実装

要素整合熱容量行列を厳密積分により求める 関数は、以下のようになる。

/* 要素整合熱容量行列：三次元：シンプレックス一次 */
void WH_ThermalFem3D__Linear__Simplex__C_arrayInJn
(double whV_x_arrayIn[4][3],
 double rho,
 double c,
 double OUT__C_arrayInJn[4][4])
/*
  入力引数:
    whV_x_arrayIn はベクトルの節点配列で、位置（節点座標）
    rho はスカラーで、質量密度
    c はスカラーで、比熱
  出力引数:
    OUT__C_arrayInJn はスカラーの配列(節点、節点)で、
      要素整合熱容量行列
*/
{
  double V;
  int In, Jn;

  WH_Fem__Simplex3D__V
    (whV_x_arrayIn[0], whV_x_arrayIn[1], 
     whV_x_arrayIn[2], whV_x_arrayIn[3], 
     &V);

  for (In = 0; In < 4; In++) {
    for (Jn = 0; Jn < 4; Jn++) {
      OUT__C_arrayInJn[In][Jn] = 1.0 / 20 * V * rho * c;
    }
  }
  for (In = 0; In < 4; In++) {
    OUT__C_arrayInJn[In][In] = 2.0 / 20 * V * rho * c;
  }
}