GetJacobianDeterminant

PURPOSE ^

GETJACOBIANDETERMINANT - computes the determinant of the jacobian for a triaelem

SYNOPSIS ^

function Jdet=GetJacobianDeterminant(triaelem,xyz_list,gauss_l1l2l3),

DESCRIPTION ^

GETJACOBIANDETERMINANT - computes the determinant of the jacobian for a triaelem

   Usage:
      Jdet=GetJacobianDeterminant(triaelem,xyz_list,gauss_coord)

   See also GETJACOBIAN, GETJACOBIANINVERT

CROSS-REFERENCE INFORMATION ^

This function calls: This function is called by:

SOURCE CODE ^

0001 function Jdet=GetJacobianDeterminant(triaelem,xyz_list,gauss_l1l2l3),
0002 %GETJACOBIANDETERMINANT - computes the determinant of the jacobian for a triaelem
0003 %
0004 %   Usage:
0005 %      Jdet=GetJacobianDeterminant(triaelem,xyz_list,gauss_coord)
0006 %
0007 %   See also GETJACOBIAN, GETJACOBIANINVERT
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

Generated on Sun 29-Mar-2009 20:22:55 by m2html © 2003