形状関数を体積座標で評価する関数は 以下のようになる。 これは、全節点の形状関数について一度に計算される。
/* 形状関数(体積座標で評価) */
void WH_Fem__Simplex3D4N__N_arrayIn
(double L0, double L1, double L2, double L3,
double OUT__N_arrayIn[4])
/*
入力引数:
L0, L1, L2, L3 はスカラーで、体積座標
出力引数:
OUT__N_arrayIn はスカラーの節点配列で、形状関数
*/
{
OUT__N_arrayIn[0] = L0;
OUT__N_arrayIn[1] = L1;
OUT__N_arrayIn[2] = L2;
OUT__N_arrayIn[3] = L3;
}