implemented pairing matrix with new npdm code, called by setting npdm_order=-1