実装

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

/* ３階テンソルとベクトルのドット積 */
void WH_Tensor3D__whT3_cdot_whV_OUT_whT
(double whT3_B[3][3][3],
 double whV_f[3],
 double OUT__whT_X[3][3])
/*
  入力引数:
    whT3_B は３階テンソル
    whV_f はベクトル
  出力引数:
    OUT__whT_Xはテンソルで、whT3_B と whV_f のドット積
*/
{
  int i, j, k;

  for (i = 0; i < 3; i++) {
    for (j = 0; j < 3; j++) {
      double tmp = 0;

      for (k = 0; k < 3; k++) {
        tmp += whT3_B[i][j][k] * whV_f[k];
      }
      OUT__whT_X[i][j] = tmp;
    }
  }
}

/* ベクトルと３階テンソルのドット積 */
void WH_Tensor3D__whV_cdot_whT3_OUT_whT
(double whV_f[3],
 double whT3_B[3][3][3],
 double OUT__whT_X[3][3])
/*
  入力引数:
    whV_f はベクトル
    whT3_B は３階テンソル
  出力引数:
    OUT__whT_Xはテンソルで、whV_f と whT3_B のドット積
*/
{
  int i, j, k;

  for (j = 0; j < 3; j++) {
    for (k = 0; k < 3; k++) {
      double tmp = 0;

      for (i = 0; i < 3; i++) {
        tmp += whV_f[i] * whT3_B[i][j][k];
      }
      OUT__whT_X[j][k] = tmp;
    }
  }
}