実装

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

/* 稜線上の法線方向成分指定の線積分の変換子 */
void WH_Fem__Isoparam2D__Volume__whV_J_l_Iedge
(WH_Fem__Shape2DType shapeType,
 int Iedge, 
 double whPt_whV_x_whV_xi[2][2],
 double a,
 double OUT__whV_J_l[2])
{
  switch (shapeType) {
  case WH_FEM__SHAPE2D__TRI:
    WH_Fem__Isoparam2D__Tri__whV_J_l_Iedge
      (Iedge, whPt_whV_x_whV_xi, a,
       OUT__whV_J_l);
    break;
  case WH_FEM__SHAPE2D__QUAD:
    WH_Fem__Isoparam2D__Quad__whV_J_l_Iedge
      (Iedge, whPt_whV_x_whV_xi, a,
       OUT__whV_J_l);
    break;
  default:
    assert(0);
    break;
  }
}