0001 function Jdet=GetJacobianDeterminant(triaelem,xyz_list,gauss_l1l2l3),
0002
0003
0004
0005
0006
0007
0008
0009 x1=xyz_list(1,1);
0010 y1=xyz_list(1,2);
0011 z1=xyz_list(1,3);
0012 x2=xyz_list(2,1);
0013 y2=xyz_list(2,2);
0014 z2=xyz_list(2,3);
0015 x3=xyz_list(3,1);
0016 y3=xyz_list(3,2);
0017 z3=xyz_list(3,3);
0018
0019 Jdet=sqrt(3)/6.0*sqrt(((y2-y1)*(z3-z1)-(z2-z1)*(y3-y1))^2+((z2-z1)*(x3-x1)-(x2-x1)*(z3-z1))^2+((x2-x1)*(y3-y1)-(y2-y1)*(x3-x1))^2);
0020
0021 if(Jdet<0),
0022 error('negative jacobian determinant!');
0023 end