On Thu, Jun 19, 2008 at 12:48 PM, Benjamin Kirk <[EMAIL PROTECTED]> wrote: >> 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 ? > > The enum LEGENDRE right now is for the radial direction in the infinite > elements, as you point out. There is not a deeper reason other than no one > has implemented them for finite elements yet. > > Adding Legendre polynomials would be pretty straightforward, though, > especially for the tensor-product elements. See for example the QUAD and > HEX sections of the case statements in src/fe/fe_lagrange_* to get a feel > for what would be involved. > > -Ben > >
Guys, I am currently using QUAD and HEX only for higher dimensions and so, this is definitely attractive for my case. I will look at the templated FE definitions for Lagrange element and create the Legendre basis on the same lines. Sometime in the future, I can create a patch with the source so that it can be useful for others who might need it. Thanks for the help. Vijay ------------------------------------------------------------------------- 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
