実装

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

/* 形状関数の自然座標についての微分 */
void WH_Fem__Shape1D__Line3N__whP_N_xi_arrayIn
(double xi,
 double OUT__whP_N_xi_arrayIn[3])
/*
  入力引数:
    xi はスカラーで、自然座標
  出力引数:
    OUT__whP_N_xi_arrayIn はスカラーの節点配列で、
      形状関数の自然座標に関する微分
*/
{
  OUT__whP_N_xi_arrayIn[0] = -0.5 + xi;
  OUT__whP_N_xi_arrayIn[1] =  0.5 + xi;
  OUT__whP_N_xi_arrayIn[2] = -2.0 * xi;
}