実装

整合熱容量行列の 節点 $In$ 節点 $Jn$ における成分を計算する 関数は、以下のようになる。

/* 整合熱容量行列 */
void WH_ThermalFemLinear__C_InJn
(double N_In, double N_Jn,
 double rho,
 double c,
 double *OUT__C_InJn)
/*
  入力引数:
    N_In, N_Jn はスカラーで、それぞれ節点 In, Jn の補間関数
    rho はスカラーで、質量密度
    c はスカラーで、比熱
  出力引数:
    OUT__C_InJn はスカラーで、
      整合熱容量行列の節点 In 節点 Jn 成分
*/
{
  *OUT__C_InJn = rho * c * N_In * N_Jn;
}