import math


def sph_kernel(r,h):
    """Evaluates the sph kernel used in gadget."""
    if r > h:
        return 0
    elif r > h/2:
        return 2*(1-r/h)**3*8/math.pi/h**3
    else:
        return 8/math.pi/h**3*(1-6*(r/h)**2+6*(r/h)**3)