実装

面積積分の変換子 $J^S$ を評価する 関数は以下のようになる。

/* 面積積分の変換子 */
void WH_Fem__Isoparam2D__Volume__J_S
(WH_Fem__Shape2DType shapeType,
 double whPt_whV_x_whV_xi[2][2],
 double a,
 double *OUT__J_S)
{
  switch (shapeType) {
  case WH_FEM__SHAPE2D__TRI:
    WH_Fem__Isoparam2D__Tri__J_S
      (whPt_whV_x_whV_xi, a,
       OUT__J_S);
    break;
  case WH_FEM__SHAPE2D__QUAD:
    WH_Fem__Isoparam2D__Quad__J_S
      (whPt_whV_x_whV_xi, a,
       OUT__J_S);
    break;
  default:
    assert(0);
    break;
  }
}