An efficient routine for computing symmetric real spherical harmonics for high orders of expansion

A numerically efficient method of constructing symmetric real spherical harmonics is presented. Symmetric spherical harmonics are real spherical harmonics with built-in invariance with respect to rotations or inversions. Such symmetry-invariant spherical harmonics are linear combinations of non-symmetric ones. They are obtained as eigenvectors of an appropriate operator, depending on symmetry. This approach allows for fast and stable computation up to very high order symmetric harmonic bases, which can be used in e.g. averaging of non-crystallographic symmetry in protein crystallography or refinement of large viruses in electron microscopy.