実装

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

/* ２階テンソルのテンソル積 */
void WH_Tensor2D__whT_otimes_whT_OUT_whT4
(double whT_X[2][2],
 double whT_Y[2][2],
 double OUT__whT4_A[2][2][2][2])
/*
  入力引数:
    whT_X, whT_Y はテンソル
  出力引数:
    OUT__whT4_Aは４階テンソルで、whT_X と whT_Y のテンソル積
*/
{
  int i, j, k, l;

  for (i = 0; i < 2; i++) {
    for (j = 0; j < 2; j++) {
      for (k = 0; k < 2; k++) {
        for (l = 0; l < 2; l++) {
          OUT__whT4_A[i][j][k][l] = whT_X[i][j] * whT_Y[k][l];
        }
      }
    }
  }
}