powered by
compute gradient of spherical harmonics functions
sph_grad_xyz(xyz, Lmax)
xyz coordinates of locations on sphere
largest degree of spherical harmonics. Current only Lmax=2 supported