実装

スカラー量を形状関数で補間する 関数は以下のようになる。

/* スカラー量の補間 */
void WH_Fem__Isoparam3D__f
(int nNodes,
 double f_arrayIn[/* nNodes */],
 double N_arrayIn[/* nNodes */],
 double *OUT__f)
/*
  入力引数:
    nNodes は整数で、要素節点数
    f_arrayIn はスカラーの節点配列で、スカラー量
    N_arrayIn はスカラーの節点配列で、形状関数
  出力引数:
    OUT__f はスカラーで、スカラー量
*/
{
  int In;

  *OUT__f = 0;
  for (In = 0; In < nNodes; In++) {
    *OUT__f += f_arrayIn[In] * N_arrayIn[In];
  }
}