OpenMP for D_psi, deriv_Sb_D_Psi and gammaXY

The preconditioned and blocked versions have been omitted. This has been tested in high-statistics runs (useevenodd=no) and with the CGMMS solver. (giving binary identical results for all permutations of (no)SSE, (no)OpenMP)