I think it would be useful to have high order Legendre basis function for high order discontinuous Galerkin discretizations. But I would need them also for tets and prisms... I know that in this case it is necessary to introduce a change of coordinates in order to keep the polynomial defined in the [-1, 1] interval. I don't think it's a big complication, is it?
Lorenzo Il giorno 19/giu/08, alle ore 21:09, Vijay S. Mahadevan ha scritto: > 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 ------------------------------------------------------------------------- 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
