c++ 線形代数 Eigenが使えない場合

Question

c++　で行列の計算を頻繁に行うプログラムを書いています．

#現状
今までは，古の文化を大事に， new deleteを用いて書いていました．

  int ***data = new int**[N];
    for(i=0;i<N;i++){
      data[i] = new int*[M];
      for(j=0;j<M;j++){
        data[i][j] = new int[M];
      }
    }

  for(i=0;i<N;i++){
	for(j=0;j<M;j++){
      delete [] data[i][j];
    }
    delete [] data[i];
  }
  delete [] data;

#改善
このやり方に限界を感じ，Eigen というライブラリがあるのを見つけそちらを使おうと思ったのですが，
Eigenでは

data(M, M)

のサイズは扱えても，初心者の私が調べた限りでは

data(N, M, M);

の3次元は扱えませんでした．
行列計算は，私が扱う中で[256][4][4]
くらいのサイズが出てきて，4*4の行列の積を256回行ったり，LU分解や転置なども出てきます．

大変恐縮ではありますが，3次元の行列計算を扱える形での最も良い方法を教えていただければと思います．

私の調べた中では，Eigenを諦めてvectorを使い自分で計算するプログラムは書くべきなのかと考えています．

Comments

[yuki23]

3次元行列の演算がしたいのですか？それとも通常の行列演算がたくさんあるだけですか？

[hromj]

後者です。
分かりづらい例で申し訳ないのですが、
4×4の正方行列の積をたくさん行います。

Answer 1

4×4の正方行列の積をたくさん行います。

それなら単に Matrix を vector などでたくさん用意すればいいかと思います

例

std::vector<Eigen::MatrixXd> data;

Comments

[hromj]

ご回答ありがとうございました。こちらの方法でやってみたいと思います。