実装

ベクトルのテンソル積を関数として以下のように実装する。

/* ベクトルのテンソル積 */
void WH_Tensor3D__whV_otimes_whV_OUT_whT
(double whV_a[3],
 double whV_b[3],
 double OUT__whT_X[3][3])
/*
  入力引数:
    whV_a、whV_b はベクトル
  出力引数:
    OUT__whT_Xはテンソルで、whV_a と whV_b のテンソル積
*/
{
  OUT__whT_X[0][0] = whV_a[0] * whV_b[0];
  OUT__whT_X[0][1] = whV_a[0] * whV_b[1];
  OUT__whT_X[0][2] = whV_a[0] * whV_b[2];

  OUT__whT_X[1][0] = whV_a[1] * whV_b[0];
  OUT__whT_X[1][1] = whV_a[1] * whV_b[1];
  OUT__whT_X[1][2] = whV_a[1] * whV_b[2];

  OUT__whT_X[2][0] = whV_a[2] * whV_b[0];
  OUT__whT_X[2][1] = whV_a[2] * whV_b[1];
  OUT__whT_X[2][2] = whV_a[2] * whV_b[2];
}