Merge branch 'merge-pull-122' into lammps-icms
From Iain Bethune at EPCC in Edinborough:
This branch introduces the Dullweber-Leimkuhler McLachlan (DLM)
integrator (see http://dx.doi.org/10.1063/1.474310) into fix_nve_sphere
and fix_nh_sphere. The existing behaviour is kept as the default, to use
the new integrator add the keyword 'dlm' to the fix command. Some
convenience functions shared between the two implementations are added
to math_extra.
In our tests, DLM provides better energy conservation and small drift /
longer timesteps than the currently implemented rotational integrator in
LAMMPS.