# Using Gauss quadrature for a discontinuous integrand

Suppose I have the following integral: $$\int_{-1}^1 \int_{-1}^1 C_if(x,y)dxdy$$ where $$C_i= \begin{cases} C_1 \quad \text{in }\Omega_1\\ C_2 \quad \text{in } \Omega_2\\ \end{cases}$$ So how can I compute it using Gauss quadrature without having to splitting the unit cell into two seperate parts or re-transforming them into another 2 unit cells. I mean, since the integrand is discontinuous over the unit cell, it seems to be not possible to use Guass quadrature directly, so what other information do I need or what treatment should I take to deal with it? Is it possible to accomplish this goal?

• So, you want to have a method that does not imply splitting the integration domain but leads to high precision? Jan 15 '16 at 22:16
• Yes, @nicoguaro Jan 16 '16 at 1:24
• I think that those two conditions compete. Jan 16 '16 at 1:38
• Hi, @nicoguaro, what do you mean "compete"? Jan 16 '16 at 1:47

Your question may sound like "one can not use Gauss quadrature for discontinuous integrand" that is not necessary a case. The problem is that a straightforward application of it (without taking into account the discontinuity) can give unsatisfactory accuracy (e.g. adding more quadrature points decrease only slowly the quadrature error).

To get an answer that will fit your goal, you should specify in your question how your sets $\Omega_i$ are described. I can give a reference for one particular case. I am aware of such integration when using level set methods, where an interface can split a domain Omega into two parts, so the sets $\Omega_i$ are given implicitly by a sign of the level set function, because the interface itself is given by the zero level set.

So if you can construct a smooth function $\phi$ (especially that $\nabla \phi$ is defined almost everywhere) such that $\Omega_1 = \{ x : \phi(x) < 0 \}$ and $\Omega_2 = \{x : \phi(x)>0\}$ then search in Google or in scientific databases for terms like "quadrature methods implicitly defined volumes". I can not recommend a specific publication, because I have not practical experiences with this problem.

Similarly the so called Volume of Fluid Method has to deal with such problem without using a level set function, but a "characteristic functions" for $\Omega_i$.