実装

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

/* 面積積分の変換子 */
void WH_Element2D__Volume__J_S
(double whPt_whV_x_whV_xi[2][2],
 double a,
 double *OUT__J_S)
/*
  入力引数:
    whPt_whV_x_whV_xi はテンソルで、Jacobian
    a は厚み
  出力引数:
    OUT__J_S は面積積分の変換子
*/
{
  double detJ;

  WH_Tensor2D__whDet_whT_OUT_s
    (whPt_whV_x_whV_xi,
     &detJ);
  *OUT__J_S = 0.5 * detJ * a;
}