実装

内部発熱項ベクトルの 節点 $In$ における成分を計算する 関数は、以下のようになる。

/* 内部発熱項ベクトル */
void WH_ThermalFemLinear__Q_body_In
(double N_In, 
 double q_body,
 double *OUT__Q_body_In)
/*
  入力引数:
    N_In はスカラーで、節点 In の補間関数
    q_body はスカラーで、内部発熱率
  出力引数:
    OUT__Q_body_In はスカラーで、
      内部発熱項ベクトルの節点 In 成分
*/
{
  *OUT__Q_body_In = q_body * N_In;
}