Hi everybody,

Does deal.II provide a way to use Gauss-Lobatto-Chebychev quadratures when 
using FE_DGQArbitraryNodes<2>?
Specifically, I want to integrate functions f(x) using Lagrange polynomials 
on the GLC support points.
The quadrature there is usually defined with a weighting function, as \int 
f(x) = \int w(x) \Phi(x) ~= sum a_i^GLC \Phi(x_i).

In order to use this in 2D, I need weights that are scaled with the 
weighting function, as a_i=a_i^GLC / w(x_i), otherwise they will not sum up 
to 1.

Is this only possible via defining a new quadrature?
Even when defining the weights a_i in a new quadrature, their sum would 
only appoximate 1, preventing the quadrature from being extended to 2D.

What is the best way to approach this issue in deal.II, and in general?

Best,
Andreas

-- 
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see 
https://groups.google.com/d/forum/dealii?hl=en
--- 
You received this message because you are subscribed to the Google Groups 
"deal.II User Group" group.
To unsubscribe from this group and stop receiving emails from it, send an email 
to dealii+unsubscr...@googlegroups.com.
For more options, visit https://groups.google.com/d/optout.

Reply via email to