実装

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

/* ベクトル三重積 */
void WH_Tensor3D__whV_times_whV_times_whV_OUT_whV
(double whV_a[3],
 double whV_b[3],
 double whV_c[3],
 double OUT__whV_d[3])
/*
  入力引数:
    whV_a, whV_b, whV_c はベクトル
  出力引数:
    OUT__whV_dはベクトルで、whV_a, whV_b, whV_c のベクトル三重積
*/
{
  double whV_ab[3];

  WH_Tensor3D__whV_times_whV_OUT_whV 
    (whV_a, whV_b,
     whV_ab);
  WH_Tensor3D__whV_times_whV_OUT_whV 
    (whV_ab, whV_c,
     OUT__whV_d);
}