実装

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

/* 体積積分の変換子 */
void WH_Fem__Isoparam3D__Volume__J_V
(WH_Fem__Shape3DType shapeType,
 double whPt_whV_x_whV_xi[3][3],
 double *OUT__J_V)
{
  switch (shapeType) {
  case WH_FEM__SHAPE3D__TETRA:
    WH_Fem__Isoparam3D__Tetra__J_V
      (whPt_whV_x_whV_xi, 
       OUT__J_V);
    break;
  case WH_FEM__SHAPE3D__HEXA:
    WH_Fem__Isoparam3D__Hexa__J_V
      (whPt_whV_x_whV_xi, 
       OUT__J_V);
    break;
  default:
    assert(0);
    break;
  }
}