行列データはすべて double へのポインタとして実現した。 ある行列を 2 次元の配列と表現したい場合と、1 次元のベクトルとして 表現したい場合とがあって煩わしいので、すべて自前で処理することにした。 具体的には matrix_util.h の中で
#ifndef scalar #define scalar double #endif typedef scalar *vector, *matrix;