Before revision:  This code involves a three-dimensional inner product between two functions, one of which requires an integral, on a non-rectangular sparse domain .
After revision:  This code involves a three-dimensional inner product between two functions, one of which requires an integral, on a non-rectangular domain containing a sparse grid .