Template:cpp/numeric/special functions/dsc sph neumann