Say I have a function $f : \mathbf{R}^3 \to \mathbf{R}$ that I wish to integrate over a tetrahedron $T \subset \mathbf{R}^3$. If $f$ was arbitrary, Gauss quadrature would be a good solution, but I happen to know that $f$ is harmonic. How much can Gauss quadrature be accelerated using this information?

For example, if $T$ was instead a sphere, evaluating $f$ once at the center of the sphere gives the exact answer by the mean value property.

A search turned up the following paper, which is interesting but generalizes the sphere case in a different direction (to polyharmonic instead of away from spheres):

Bojanov and Dimitrov, Gaussian extended cubature formulae for polyharmonic functions


1 Answer 1


I found something that might be interesting. http://www.math.kth.se/~gbjorn/exact.pdf

I hope this helps, Tom

  • $\begingroup$ That's an interesting paper, but it looks like it and its references only treat integrals of differential operators of harmonic functions. Do you know if they can be used for straight integrals? $\endgroup$ Mar 4, 2013 at 17:38
  • $\begingroup$ I'm wondering if introducing a quadrature formula with the so-called "Poisson kernel" (en.wikipedia.org/wiki/Poisson_kernel) could help... Otherwise I know that some xfem techiques use harmonic functions to enrich the FE space, and therefore should use specific quadrature methods for integrating the variational forms (?). $\endgroup$
    – Tom
    Mar 4, 2013 at 20:46

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge you have read our privacy policy.

Not the answer you're looking for? Browse other questions tagged or ask your own question.