Definition of custom types required by the code.
real dhRoma(real x, real h)
Discrete delta function from Roma et al. (1999).
double real
Is a float or a double depending on the machine precision.
real delta(real x, real y, real h)
Two-dimensional discrete delta function.