製作一個具有兩個固定長度和一個可以變化的3D矢量

Question

我希望構建一個3D結構，它將由一個具有i = 1..N個網站的數組組成，在每個網站中都有一個2D矢量，該向量有一個有限的dim k = 1 ... L，並且在模擬過程中發生變化，從M = 0開始，使用push_back命令更新L元素。

我已經添加了一個解釋結構的圖，該如何寫這個東西，請注意，重要的是我可以通過索引來訪問這個結構，比如struct [i] [j] [k] 。

注意：這個結構存儲了N個實現的模擬數據，每個實現中有M個時間步長可以不同（在運行期間發生這種情況），並且在每個時間步中有一個L個點的點陣是更新。

我能做到它在2D情況下，但不是3D，我至今是：

array <vector <double>, N> Array; 
array <array <vector <double>, N>, L> Brray; 
for (int i=0; i<N; i++) 
{ 
    // calculate cal1 for 0'th index 
    Array[i].push_back(i) = cal1; 
    j=1; 
    while (j<M) 
    { 
     // calculate cal1 
     Array[i].push_back(Array[i][j-1]+cal1); 
     for (int k=0; k<L; j++) 
     { 
      // calculate cal2 
      Brray[i][j][k] = cal2; // this is the line I have a problem 
     } 
     //perform cal3 that decides whether there will be another iteration 
     if (cal3==1) 
     { 
      j++; 
     } 
    } 
} 

Answer 1

你需要：

• 含有
  • 含有
    • 含有常量位置L的陣列雙值
  可變大小爲M的向量固定大小N的陣列

您應該使用array<vector<array<double, L> >, N>。但是您需要在每個迭代中分配一個新的array<double, L>。

您的代碼將成爲（僅Brray部分編輯）：

array <vector <double>, N> Array; 
array <vector <array<double, L> >, N> Brray; // CHANGED 
for (int i=0; i<N; i++) 
{ 
    // calculate cal1 for 0'th index 
    Array[i].push_back(i) = cal1; 
    M = value; // value is different for every i 
    j=1; 
    while (j<M) 
    { 
     Brray[i].push_back(array<double, L>()); //ADDED: add a new array to vector 
     // calculate cal1 
     Array[i].push_back(Array[i][j-1]+cal1); 
     for (int k=0; k<L; j++) 
     { 
      // calculate cal2 
      Brray[i][j][k] = cal2; // this is the line I have a problem - should work now 
     } 
     j++; 
    } 
} 

我無法測試上面的代碼，但在C++ 11模式（鐺）這一個正常工作：

#include <iostream> 
#include <array> 
#include <vector> 


int main() { 
    constexpr int N=5, L=8; 
    std::array<std::vector<std::array<double, L> >, N> arr; 
    for (int i=0; i<N; i++) { 
     const int M = 2*i; 
     for (int j=0; j<M; j++) { 
      arr[i].push_back(std::array<double, L>()); 
      for (int k=0; k< L; k++) { 
       arr[i][j][k] = i+j+k; 
      } 
     } 
    } 

    for (int i=0; i<N; i++) { 
     const int M = 2*i; 
     std::cout << i << std::endl; 
     for (int j=0; j<M; j++) { 
      std::cout << "\t" << j << " :"; 
      for (int k=0; k< L; k++) { 
       std::cout << " " << arr[i][j][k]; 
      } 
      std::cout << std::endl; 
     } 
    } 

    return 0; 
} 

Answer 2

vector<vector<vector<double>>>，或者，由於兩個尺寸是編譯時間常數，array<vector<array<double, L>>, N>。

編輯

要使用此，你會做沿着這些路線的東西：

array<vector<array<double, L>>, N> dat; 

此創建初始數組來保存所有的數據。在每次實現開始時，您已經有一個空向量來保存M個結果。你既可以使用

dat[i].resize(M); 

，如果你在你的循環開始知道M或

dat[i].emplace_back(); 

添加一個當您需要。然後這將給你你的陣列與L網站，你可以使用索引（dat[i][j][k]）或迭代器填充它。