実装

面積を求める 関数は以下のようになる。

/* 面積 */
void WH_Fem__Simplex2D__S
(double whV_x0[2], 
 double whV_x1[2], 
 double whV_x2[2],
 double *OUT__S)
/*
  入力引数:
    whV_x0, whV_x1, whV_x2 はベクトルで、各頂点の位置
  出力引数:
    OUT__S はスカラーで、面積
*/
{
  double x0 = whV_x0[0];
  double x1 = whV_x1[0];
  double x2 = whV_x2[0];
  double y0 = whV_x0[1];
  double y1 = whV_x1[1];
  double y2 = whV_x2[1];

  *OUT__S 
    = 0.5 * ((x1 - x0) * (y2 - y0) - (x2 - x0) * (y1 - y0));
}