METHOD=CHARGE_CONSTRAIN with PURIFICATION=NONE must be set, otherwise Error message
IF for switching between charge_constrained ADMM and normal basis_projection in the correct file, so that there is no double coding.
Tr(DS) implemented, code compiles, but crashes when calculating the Trace.