GitLab

angles are in degrees (unlike *every* other angle force out there), and
stripUnits automatically reduces all angles to radians if they come in with
units.

The approach here is to *slightly* change the API, so that
AmoebaAngleForce.addAngle interprets input angles <2*pi as radians, and >2*pi as
degrees. This is heuristic, but should work in every case out in the wild so
far.

I've also updated the documentation to reflect this behavior, and fixed the
units attached to the return value of AmoebaAngleForce.getAngleParameters() to
return degrees instead of radians.

o Replaced Cartesians in the reciprocal and self PME terms.
o Added non-PME spherical harmonics field generation and force evaluation routines.
o Added spherical harmonic potential evaluation, to allow removal of the Cartesian moments.

o Permanent and induced moment energies, forces and torques are implemented, for PME direct space.
o The code is activated by a preprocessor macro at the top of AmoebaReferenceMultipoleForce.h, and can be completely disabled.
o Currently, the Cartesian momenents are stored alongside the spherical harmonic moments, but these can eventually be removed when the reciprocal and self terms have been adapted to use spherical harmonics (in progress) and the non-PME routines have been coded.