On Mon, Sep 20, 2010 at 9:48 AM, yunfei zhu <[email protected]> wrote: > > I noted one of the comments in quadrature_gauss_3D.C, the tetrahedral > quadrature rules are taken from pg. 222 of "The finite element method" vol.1 > ed.5 by Zienkiewicz & Taylor, > so I checked the details in the book. I find that for each of the quadrature > point, there are 4 tetrahedral coordinates and the sum of the weights are > equal to 1.
Those are probably barycentric coordinates, one of which is not independent, i.e. (x,y,z,1-x-y-z). We have simply converted those to (x,y,z). I don't have the book in front of me at the moment, so I'm not totally sure, but you can see more about the barycentric coordinates in fe_lagrange_shape_3D.C, around line 78. -- John ------------------------------------------------------------------------------ Start uncovering the many advantages of virtual appliances and start using them to simplify application deployment and accelerate your shift to cloud computing. http://p.sf.net/sfu/novell-sfdev2dev _______________________________________________ Libmesh-users mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-users
