Hello Stefano,

FE_Q(QGaussLobatto(n)) works just fine. All quadratures in deal.II are defined on [0,1].

Best Regards,
Markus



Am 12.12.2011 11:35, schrieb Stefano Zampini:
Hi,

I'm a new user. I'm wondering how to use spectral quadrilateral elements based on tensor product lagrangian polinomials on Gauss Lobatto Legendre points. Searching the mailing list and the classes, I found two different constructors to do this, namely: fe_Q(GLLpoints) (with GLL points on [-1,1] which must be shifted to [0,1]) or fe_DGQArbtrarynodes(...)? The latter solution seems to be the recomended one. Why do I need to resort to DG methods since I need a conforming triangulation?

Thanks,
--
Stefano


_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii

Reply via email to