//---------------------------------------------------------------------------
//计算4节点四面体单元的体积(节点顺序满足右手定则）
double cal_volume( Node *node1, Node *node2, Node *node3, Node *node4 )
{
double volume;
double aai = cal_det( node2, node3, node4 );
double aaj = cal_det( node3, node4, node1 );
double aak = cal_det( node4, node1, node2 );
double aal = cal_det( node1, node2, node3 );
volume = 1.0/6.0*(aai-aaj+aak-aal) ;
return volume ;
}
double cal_gconstf(double **var, const int D)
{
double gconst, tmp;
tmp = cal_det(var, D);
if (tmp == 0) {
fprintf(stderr, "WARNING : det is 0!\n");
return 0;
}
gconst = D * log(M_2PI);
gconst += tmp;
return (gconst);
}
