Final bug fix for the input of (Dxx, Dyy, Dzz, Dxy, Dxz, Dyz) ellipsoidal diffusion parameters.
The numpy eigenvalue decomposition has been replaced by SVD as the eigenvalue decomposition could
generate a proper rotation matrix when 2 eigenvalues are the same. The test suite now passes fully.