On Thu, Jun 19, 2008 at 1:42 PM, Vijay S. Mahadevan <[EMAIL PROTECTED]> wrote: > Hi, > > I was wondering if there is something special that needs to be > specified in order to use LEGENDRE FE basis functions for say > diffusion problems. I am particularly interested in basis functions > that are orthogonal to each other by nature since it simplifies the > mass matrix structure. Hence, when I try to create elements with > legendre basis, I end up with a runtime error because all the > functions in FEInterface class never have a switch statement for > LEGENDRE option.
They've never been implemented, though they are quite useful in certain special situations. I think they're only orthogonal on affine elements (?). For steady diffusion problems there are also "integrated legendre" polynomials so that the diffusion matrix instead of the mass matrix is simplified ... > So my question is, can I use LEGENDRE basis at all for non-infinite > elements ? Would it require a lot of work to add this support ? Or is > there a deeper reason why this was intentionally left out of the > implementation and was written only for infinite elements alone ? It depends how you define "a lot". I would definitely call it non-trivial though. We have some interesting templated FE definitions that take a while to get used to, but other than that it shouldn't be too bad. -- John ------------------------------------------------------------------------- Check out the new SourceForge.net Marketplace. It's the best place to buy or sell services for just about anything Open Source. http://sourceforge.net/services/buy/index.php _______________________________________________ Libmesh-users mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-users
