Major changes:

• New feature: Second derivatives of internal coordinates w/r.t. Cartesian coordinates and transformation of Hessians from Cartesian to internal coordinates

• Improved performance and stability of constrained optimizations that set --enforce