How to choose Gauss-Lobatto Quadrature?

1

Hi all,
I want to use Gauss-Lobatto when assembling a 1D problem.
I can do it by hacking FIAT/quadrature.py to change make_quadrature
as
elif ref_el.get_shape() == reference_element.LINE:
#return GaussJacobiQuadratureLineRule(ref_el, m)
return GaussLobattoLegendreQuadratureLineRule(ref_el, m)
but that is unsatisfactory. Any better suggestions?
Thanks,

2

@bleyerj has created a monkey patch for using Gauss-Lobatto here: https://gitlab.enpc.fr/navier-fenics/fenics-optim/-/blob/bcd1882f34ffbb9e90fb1df50c184d28daab9162/fenics_optim/custom_quadrature_schemes.py

2 Likes
3

Thanks!
Will follow that lead.

4

Yep - that worked. Thanks @dokken + @bleyerj