実装

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

/* 稜線上の線積分の変換子 */
void WH_Fem__Isoparam2D__Volume__J_l_Iedge
(WH_Fem__Shape2DType shapeType,
 int Iedge, 
 double whPt_whV_x_whV_xi[2][2],
 double a,
 double *OUT__J_l)
{
  switch (shapeType) {
  case WH_FEM__SHAPE2D__TRI:
    WH_Fem__Isoparam2D__Tri__J_l_Iedge
      (Iedge, whPt_whV_x_whV_xi, a,
       OUT__J_l);
    break;
  case WH_FEM__SHAPE2D__QUAD:
    WH_Fem__Isoparam2D__Quad__J_l_Iedge
      (Iedge, whPt_whV_x_whV_xi, a,
       OUT__J_l);
    break;
  default:
    assert(0);
    break;
  }
}