実装

熱流束ベクトルを求める 関数は、２次元と３次元問題についてそれぞれ以下のようになる。

/* 熱流束ベクトル：２次元 */
void WH_Thermal2D__whV_q
(double whT_lambda[2][2],
 double whP_T_whV_X[2],
 double OUT__whV_q[2])
/*
  入力引数:
    whT_lambda はテンソルで、熱伝導率テンソル
    whP_T_whV_X はベクトルで、温度勾配
  出力引数:
    OUT__whV_q はベクトルで、熱流束
*/
{
  WH_Tensor2D__whT_cdot_whV_OUT_whV
    (whT_lambda, whP_T_whV_X,
     OUT__whV_q);
}

/* 熱流束ベクトル：３次元 */
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);
}