実装

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

/* ベクトルの外積 */
void WH_Tensor3D__whV_times_whV_OUT_whV
(double whV_a[3],
 double whV_b[3],
 double OUT__whV_c[3])
/*
  入力引数:
    whV_a, whV_b はベクトル
  出力引数:
    OUT__whV_cはベクトルで、whV_a と whV_b の外積
*/
{
  OUT__whV_c[0] = whV_a[1] * whV_b[2] - whV_a[2] * whV_b[1];
  OUT__whV_c[1] = whV_a[2] * whV_b[0] - whV_a[0] * whV_b[2];
  OUT__whV_c[2] = whV_a[0] * whV_b[1] - whV_a[1] * whV_b[0];
}