実装：テンソル標記

熱流束ベクトルを求める 関数は、以下のようになる。

/* 熱流束ベクトル：三次元 */
void WH_Thermal3D__whV_q
(double whT_lambda[3][3],
 double whP_T_whV_x[3],
 double OUT__whV_q[3])
/*
  入力引数:
    whT_lambda はテンソルで、熱伝導率テンソル
    whP_T_whV_x はベクトルで、温度勾配ベクトル
  出力引数:
    OUT__whV_q はベクトルで、熱流束ベクトル
*/
{
  WH_Tensor3D__whT_cdot_whV_OUT_whV
    (whT_lambda, whP_T_whV_x,
     OUT__whV_q);
  WH_Tensor3D__whV_mul_s_OUT_whV
    (OUT__whV_q, -1.0,
     OUT__whV_q);
}