実装：マトリクス標記

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

/* 熱流束の列ベクトル：三次元 */
void WH_Thermal3D__whCv_q
(double whM_lambda[3][3],
 double whCv_whP_T_whV_x[3],
 double OUT__whCv_q[3])
/*
  入力引数:
    whM_lambda は行列で、熱伝導率の行列
    whCv_whP_T_whV_x は列ベクトルで、温度勾配の列ベクトル
  出力引数:
    OUT__whCv_q は列ベクトルで、熱流束の列ベクトル
*/
{
  WH_Matrix__whM_mul_whCv_OUT_whCv
    (3, 3, whM_lambda, 
     3, whCv_whP_T_whV_x,
     3, OUT__whCv_q);
  WH_Matrix__whCv_mul_s_OUT_whCv
    (3, OUT__whCv_q,
     -1.0,
     3, OUT__whCv_q);
}