> > Only with some difficulty - we don't have those elements right now, so > you'd have to create a new FEType, code the shape functions, and add a > new macroelement quadrature rule for the integration. > > None of that is easy, but it's definitely doable. Coincidentally, > that all describes my own first contributions to libMesh, in the > CLOUGH finite elements. You'd have it a bit easier than I did: the > Clough-Tocher shape functions are hairier, you've now got a > macroelement quadrature rule (QClough) to use as an example, and now > adaptivity works with arbitrary C0 elements without any special effort > on your part. > > That may be more work than you wanted, but if you do give it a try, > let us know and we'd be happy to add it in to the main library. > --- > Roy > > Sounds a difficult work. I have a question on the quadrature rule of integration, doesn't the libmesh use the unified gauss-legendre quadrature for all elements? I knew some FEA packages use the gauss-legendre for any element, or it is a different concept between the macroelement and the single element?
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
