実装

表面上の法線方向成分指定の面積積分の変換子 $whV{J^S}$ を評価する 関数は以下のようになる。

/* 表面上の法線方向成分指定の面積積分の変換子 */
void WH_Fem__Isoparam3D__Volume__whV_J_S_Iface
(WH_Fem__Shape3DType shapeType,
 int Iface, 
 double whPt_whV_x_whV_xi[3][3],
 double OUT__whV_J_S[3])
{
  switch (shapeType) {
  case WH_FEM__SHAPE3D__TETRA:
    WH_Fem__Isoparam3D__Tetra__whV_J_S_Iface
      (Iface, whPt_whV_x_whV_xi,
       OUT__whV_J_S);
    break;
  case WH_FEM__SHAPE3D__HEXA:
    WH_Fem__Isoparam3D__Hexa__whV_J_S_Iface
      (Iface, whPt_whV_x_whV_xi,
       OUT__whV_J_S);
    break;
  default:
    assert(0);
    break;
  }
}