Changes in Step2 of integrator_omelyan 

error bounds in iteration loop and some tidying up of varaibles and definitions.
This is the normal Newton iteration method. The rotation matrix inverse is added
as well since the torque in the lab frame has to be rotated opposite to the rotation
of the molecule to get its representation right in the body fixed frame.