Forum Message

 

 

We have moved the forum to https://forum.ngsolve.org . This is an archived version of the topics until 05/05/23. All the topics were moved to the new forum and conversations can be continued there. This forum is just kept as legacy to not invalidate old links. If you want to continue a conversation just look for the topic in the new forum.

Notice

The forum is in read only mode.

Regarding structure of L2 basis functions

More
3 years 2 months ago #3506 by alex_sch
I am currently using ngsolve for my DG solver. For certain post-processing, I am looking for the structure of basis functions of order p on elements in terms of monomials i.e.

(Basis_function) = \sum( coeff * x^m y^n) where m+n <= p.


where Basis_function represents one of the basis functions for order p, coeff represents the coefficient of monomial x^m y^n and \sum represents the summation over monomials.

Currently, I have basis functions and their derivative at the quadrature points. Since many commonly used basis functions such as Legendre polynomials are generated via recursion and thus, gives us the possibility to look into their structure. I wonder that if a similar approach can be taken with ngsolve?
More
3 years 2 months ago #3510 by mneunteufel
Hi alex_sch,

for the triangular L2 shape functions (see L2HighOrderFE_Shape<ET_TRIG>::T_CalcShape, line ~385 in fem/l2hofe_impl.hpp) NGSolve uses a Dubiner basis (class DubinerBasis in fem/recursive_pol.hpp line ~2584). The evaluation is done via a recursion involving Legendre and Jacobi polynomials.

Best,
Michael
Time to create page: 0.171 seconds