実装：マトリクス標記

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

/* 熱流束の列ベクトル：三次元：等方性材料 */
void WH_Thermal3D__whCv_q_IF_isotropic
(double lambda,
 double whCv_whP_T_whV_x[3],
 double OUT__whCv_q[3])
/*
  入力引数:
    lambda はスカラーで、熱伝導率
    whCv_whP_T_whV_x は列ベクトルで、温度勾配の列ベクトル
  出力引数:
    OUT__whCv_q は列ベクトルで、熱流束の列ベクトル
*/
{
  WH_Matrix__s_mul_whCv_OUT_whCv
    (-1.0 * lambda, 
     3, whCv_whP_T_whV_x,
     3, OUT__whCv_q);
}