(M3R + i*M3I) = (M1R + i*M1I)*D*(M2RT - i*M2IT)M3R = M1R*D*M1RT + M1I*D*M1IT; M3I = -M1R*D*M1IT + M1I*D*M1RT