実装

要素集中化熱容量ベクトルを厳密積分により求める 関数は、以下のようになる。

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

  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++) {
    OUT__C_arrayIn[In] = 1.0 / 4 * V * rho * c;
  }
}