> 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 ------------------------------------------------------------------------- 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
