Ich würde versuchen, anstelle eines eigenen COMPLEX-Typs den vordefinierten Typ std::complex<T> zu verwenden, in deinem Fall als std::complex<double>. Dann brauchst du die Rechenregeln für komplexe Zahlen nicht noch mal nachzuprogrammieren, das bringt nämlich keinen wissenschaftlichen Gewinn und ist nur eine weitere Fehlerquelle.
Danke - ja das werde ich auf alle fälle - das sollte nur ein versuch sein. Aber prinzipiell würde mich interessieren ob das so in Ordnung ginge oder ob man das doch anders machen sollte....danke