> > Gauss-Legendre rules of high enough order give exact results for > polynomial non-macroelements, but on a macroelement the reduced > smoothness between subelements hurts you. The easiest (albeit not the > most efficient) solution for piecewise-polynomial macroelements is to > just combine a (translated and reweighted) Gauss rule on each > subelement. > > libMesh has a bunch of non-Gauss-Legendre rules for general elements > available too; some for verification and underintegration, some just > because there are multiple ways to do quadrature in 2D/3D with > different tradeoffs. Check the source code first if you're curious; > John has most of it pretty well documented. > --- > Roy > Interesting, I appreciate your clear explanation, yes I just know something about basic non-macroelement from the textbook. it is obvious that libmesh is not only a software but also a wonderful modern FEA tutorial. there are already lots of libmesh source codes on my list. Thanks!
Liang ------------------------------------------------------------------------------ Download Intel® Parallel Studio Eval Try the new software tools for yourself. Speed compiling, find bugs proactively, and fine-tune applications for parallel performance. See why Intel Parallel Studio got high marks during beta. http://p.sf.net/sfu/intel-sw-dev _______________________________________________ Libmesh-users mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-users
