実装

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

/* ベクトルのテンソル積 */
void WH_Tensor2D__whV_otimes_whV_OUT_whT
(double whV_a[2],
 double whV_b[2],
 double OUT__whT_X[2][2])
/*
  入力引数:
    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[1][0] = whV_a[1] * whV_b[0];
  OUT__whT_X[1][1] = whV_a[1] * whV_b[1];
}