Add a method to get the second derivatives of the directors

with respect to the quaternion coordinates.  Or, in other words,
the second derivative of the map from quaternions to 
orthogonal matrices.

[[Imported from SVN: r8407]]