実装

要素内部発熱項ベクトルを厳密積分により求める 関数は、以下のようになる。

/* 要素内部発熱項ベクトル：三次元：シンプレックス一次 */
void WH_ThermalFem3D__Linear__Simplex__Q_body_arrayIn
(double whV_x_arrayIn[4][3],
 double q_body,
 double OUT__Q_body_arrayIn[4])
/*
  入力引数:
    whV_x_arrayIn はベクトルの節点配列で、位置（節点座標）
    q_body はスカラーで、内部発熱率
  出力引数:
    OUT__Q_body_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__Q_body_arrayIn[In] = 1.0 / 4 * V * q_body;
  }
}