実装

形状関数を自然座標で評価する関数は 以下のようになる。 これは、全節点の形状関数について一度に計算される。

/* 形状関数 */
void WH_Fem__Shape1D__Line2N__N_arrayIn
(double xi,
 double OUT__N_arrayIn[2])
/*
  入力引数:
    xi はスカラーで、自然座標
  出力引数:
    OUT__N_arrayIn はスカラーの節点配列で、形状関数
*/
{
  OUT__N_arrayIn[0] = 0.5 * (1.0 - xi);
  OUT__N_arrayIn[1] = 0.5 * (1.0 + xi);
}