実装

行列は二次元配列変数として実装される。

/*
  行列変数 whM_A は、以下のように定義される。

#define NCOLUMNS 1000
#define NROWS 1000

double whM_A[NROWS][NCOLUMNS];

  配列 whM_A はサイズがNROWS(行)、NCOLUMNS(列)の行列である。
*/

/* 行列用のマクロ */

#define WH_MATRIX__AT(A, nRows, nColumns, i, j) \
  ((double*)(A))[(nColumns) * (i) + (j)]