OpenMP parallellization of 2DM calculation