fiat team mailing list archive

["Shi Jin" <jinzishuai@xxxxxxxxx>]

Hi there,

I have been using FIAT as a component from PETSc and I love it. Thank you.
Recently I have the need to integrate a discontinuous function in a
tetrahedral finite element (integrate only inside a sphere instead of the
whole domain) and the Gaussian quadrature rules does not work very well,
even if I put lots of quadrature points. My naive idea is that if I use the
trapezoidal quadrature rules, the accuracy of integration is guaranteed when
the number of quadrature points is large enough. Therefore, I am wondering
if it is possible to create trapezoidal quadrature rules using FIAT.

Note that I am not very familiar with the mathematics of quadratures. I
originally thought the one created by PETSc using FIAT is default to
Gaussian quadratures. But then I read the FIAT manual and it says  "FIAT
implements arbitrary-order collapsed quadrature, as discussed in Karniadakis
and Sherwin [], for the simplex of dimension one, two, or three." I am not
sure what "collapsed quadrature" is. Does it include Gaussian and
Trapezoidal quadratures? How do I choose from them? I would greatly
appreciate it if some guidance can be provided. Also, could more information
be provided about the Karniadakis and Sherwin paper, such as title, journal
and year? Maybe better quadrature rules exist for the integration of
discontinuous functions and I would love to hear some advice.

Thank you very much.

Shi Jin