実装

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

/* 形状関数 */
void WH_Fem__Shape2D__Quad8N__N_arrayIn
(double whV_xi[2],
 double OUT__N_arrayIn[8])
/*
  入力引数:
    whV_xi はベクトルで、自然座標
  出力引数:
    OUT__N_arrayIn はスカラーの節点配列で、形状関数
*/
{
  double xi = whV_xi[0];
  double eta = whV_xi[1];

  OUT__N_arrayIn[0] 
    = 0.25 * (1.0 - xi) * (1.0 - eta) * (-1.0 - xi - eta);
  OUT__N_arrayIn[1] 
    = 0.25 * (1.0 + xi) * (1.0 - eta) * (-1.0 + xi - eta);
  OUT__N_arrayIn[2] 
    = 0.25 * (1.0 + xi) * (1.0 + eta) * (-1.0 + xi + eta);
  OUT__N_arrayIn[3] 
    = 0.25 * (1.0 - xi) * (1.0 + eta) * (-1.0 - xi + eta);

  OUT__N_arrayIn[4] = 0.5 * (1.0 - xi * xi) * (1.0 - eta);
  OUT__N_arrayIn[5] = 0.5 * (1.0 + xi) * (1.0 - eta * eta);
  OUT__N_arrayIn[6] = 0.5 * (1.0 - xi * xi) * (1.0 + eta);
  OUT__N_arrayIn[7] = 0.5 * (1.0 - xi) * (1.0 - eta * eta);
}