为下/上三角矩阵动态分配内存并使其成为 C/C++ 中的一维数组

如何解决为下/上三角矩阵动态分配内存并使其成为 C/C++ 中的一维数组

我正在做一个项目，需要我创建 15 个维度为 100000*100000 的双数组。它们都是下三角矩阵。因此，仅仅为 1 个这样的数组耗尽内存并给出错误 terminate called after throwing an instance of 'std::bad_alloc' what(): std::bad_alloc

也就不足为奇了

我应该如何删除这些矩阵中的零并像这样编码它们

#define WIDTH 100000
#define HEIGHT 100000
int jimmy [HEIGHT * WIDTH];
int n,m;

int main ()
{
  for (n=0; n<HEIGHT; n++)
    for (m=0; m<WIDTH; m++)
   {
     jimmy[n*WIDTH+m]=(n+1)*(m+1);
  }
}

这是我正在处理的代码片段。我使用了以下未显示的矩阵。它们都是下三角形。我不需要所有这些。

double T  = 0.3;   // final time.
double iN = 100; // samples(I want it 100000) 
double h  = T/N;   // step-size

double a_i11[iN][iN],a_i12[iN][iN],a_i13[iN][iN];
double a_i21[iN][iN],a_i22[iN][iN],a_i23[iN][iN];
double a_i31[iN][iN],a_i32[iN][iN],a_i33[iN][iN];
double a_f[iN][iN];

double b_i12[iN][iN],b_i13[iN][iN];
double b_i22[iN][iN],b_i23[iN][iN];
double b_p[iN][iN];

This answer 展示了存储下三角矩阵的最佳方法。谁能告诉我结合这两种方法的最有效方法？我的最后一个目标是将此 C++ 代码转换为 CUDA 代码。