実装

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

/* 形状関数(面積座標で評価) */
void WH_Fem__Simplex2D3N__N_arrayIn
(double L0, double L1, double L2,
 double OUT__N_arrayIn[3])
/*
  入力引数:
    L0, L1, L2 はスカラーで、面積座標
  出力引数:
    OUT__N_arrayIn はスカラーの節点配列で、形状関数
*/
{
  OUT__N_arrayIn[0] =  L0;
  OUT__N_arrayIn[1] =  L1;
  OUT__N_arrayIn[2] =  L2;
}