Re: [Libmesh-users] Legendre polynomial basis

[Vijay S. M. <vi...@gm...>]

On Thu, Jun 19, 2008 at 12:48 PM, Benjamin Kirk <ben...@na...> 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