実装

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

/* 表面上の面積積分の変換子 */
void WH_Fem__Isoparam3D__Volume__J_S_Iface
(WH_Fem__Shape3DType shapeType,
 int Iface, 
 double whPt_whV_x_whV_xi[3][3],
 double *OUT__J_S)
{
  switch (shapeType) {
  case WH_FEM__SHAPE3D__TETRA:
    WH_Fem__Isoparam3D__Tetra__J_S_Iface
      (Iface, whPt_whV_x_whV_xi,
       OUT__J_S);
    break;
  case WH_FEM__SHAPE3D__HEXA:
    WH_Fem__Isoparam3D__Hexa__J_S_Iface
      (Iface, whPt_whV_x_whV_xi,
       OUT__J_S);
    break;
  default:
    assert(0);
    break;
  }
}