Introduction of a New Theory for the Calculation of Magnetic Coupling Based on Spin–Flip Constricted Variational Density Functional Theory. Application to Trinuclear Copper Complexes which Model the Native Intermediate in Multicopper Oxidases