実装

直交デカルト座標系への座標変換テンソルを計算する関数は 以下のようになる。

/* 座標系について、その座標変換テンソルを得る */
void WH_Tensor3D__transformationTensor
(double whCS_CS[3][3],
 double OUT__whT_P[3][3])
/*
  入力引数:
    whCS_CS は座標系
  出力引数:
    OUT__whT_P は座標変換テンソル
*/
{
  OUT__whT_P[0][0] = whCS_CS[0][0];
  OUT__whT_P[0][1] = whCS_CS[0][1];
  OUT__whT_P[0][2] = whCS_CS[0][2];

  OUT__whT_P[1][0] = whCS_CS[1][0];
  OUT__whT_P[1][1] = whCS_CS[1][1];
  OUT__whT_P[1][2] = whCS_CS[1][2];

  OUT__whT_P[2][0] = whCS_CS[2][0];
  OUT__whT_P[2][1] = whCS_CS[2][1];
  OUT__whT_P[2][2] = whCS_CS[2][2];
}