double spherical harmonics method