実装

３つのベクトルのトリアドを関数として以下のように実装する。

/* ３つのベクトルのトリアド */
void WH_Tensor3D__whV_otimes_whV_otimes_whV_OUT_whT3
(double whV_a[3],
 double whV_b[3],
 double whV_c[3],
 double OUT__whT3_A[3][3][3])
/*
  入力引数:
    whV_a, whV_b, whV_c はベクトル
  出力引数:
    OUT__whT3_Aは３階テンソルで、whV_a, whV_b, whV_c のトリアド
*/
{
  int i, j, k;

  for (i = 0; i < 3; i++) {
    for (j = 0; j < 3; j++) {
      for (k = 0; k < 3; k++) {
        OUT__whT3_A[i][j][k] = whV_a[i] * whV_b[j] * whV_c[k];
      }
    }
  }
}