実装

もし、 ベクトル未知数の自由度が各軸方向の3成分により表現される場合、 補間関数ベクトルを計算する 関数は以下のようになる。 これは、全節点の全自由度について一度に計算される。

/* 補間関数ベクトル */
void WH_Fem__Isoparam3D__whV_N_arrayInId
(int nNodes, 
 double N_arrayIn[/* nNodes */],
 double OUT__whV_N_arrayInId[/* nNodes */][3][3])
/*
  入力引数:
    nNodes は整数で、要素節点数
    N_arrayIn はスカラーの節点配列で、形状関数
  出力引数:
    OUT__whV_N_arrayInId はベクトルの節点自由度配列で、
      補間関数ベクトル
*/
{
  int In;

  for (In = 0; In < nNodes; In++) {
    OUT__whV_N_arrayInId[In][0][0] = N_arrayIn[In];
    OUT__whV_N_arrayInId[In][0][1] = 0.0;
    OUT__whV_N_arrayInId[In][0][2] = 0.0;

    OUT__whV_N_arrayInId[In][1][0] = 0.0;
    OUT__whV_N_arrayInId[In][1][1] = N_arrayIn[In];
    OUT__whV_N_arrayInId[In][1][2] = 0.0;

    OUT__whV_N_arrayInId[In][2][0] = 0.0;
    OUT__whV_N_arrayInId[In][2][1] = 0.0;
    OUT__whV_N_arrayInId[In][2][2] = N_arrayIn[In];
  }
}