double nablaw(double r) { final h = smoothingRadius; final nabla = (945.0 / (32.0 * math.pi * math.pow(h, 9))) * (h * h - r * r) * (7 * r * r - 3 * h * h); return nabla; }