実装

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

/* 体積積分の変換子 */
void WH_Fem__Isoparam3D__Tetra__J_V
(double whPt_whV_x_whV_xi[3][3],
 double *OUT__J_V)
/*
  入力引数:
    whPt_whV_x_whV_xi はテンソルで、Jacobian
  出力引数:
    OUT__J_V はスカラーで、体積積分の変換子
*/
{
  double detJ;

  WH_Tensor3D__whDet_whT_OUT_s
    (whPt_whV_x_whV_xi,
     &detJ);
  *OUT__J_V = 1.0 / 6 * detJ;
}