C++ Operatorüberladung für eine Matrixklasse mit reellen und komplexen Matrizen

Question

Ich versuche, eine Matrixklasse zu implementieren, die sowohl echte als auch komplexe Matrizen verarbeiten kann. Ich habe Probleme, wenn ich versuche, den Multiplikationsoperator zu überladen. Genauer gesagt, wenn ich versuche, eine doppelte Matrix mit einer komplexen Matrix zu multiplizieren (in dieser Reihenfolge). Das Ergebnis sollte komplex sein, aber der Operator * ist ein Mitglied einer Doppelmatrix. In diesem Fall sehe ich nicht, wie man eine komplexe Matrix zurückgibt (ich habe versucht, einen Freundoperator zu verwenden, das scheint auch nicht zu funktionieren). Hier ist der relevante Code-Schnipsel:

template <class V> Matrix<T> operator *(const Matrix<V> &b) 
{ 
    long i, j, k; T temp; Matrix<T> c(mRows, b.Cols()); 

    for (i = 0; i < mRows; i++) 
     for (j = 0; j < c.Cols(); j++) 
     { 
      for (temp = 0, k = 0; k < mCols; k++) 
       temp += this->Element(i, k)*b.Element(k, j); 

      c(i, j) = temp; 
     } 

     return c; 
} 

Also, wenn A real ist, und B und C sind komplex, C = B A funktioniert gut, aber C = A B nicht. Insbesondere weist der Compilerfehler die Zeile "temp + =" darauf hin, dass kein globaler Operator gefunden wurde, der den Typ komplex (oder keine akzeptable Konvertierung) akzeptiert. Ich möchte vermeiden, die Vorlage für reale und komplexe Typen vollständig zu spezialisieren, gibt es einen Weg um dies zu umgehen?

Vielen Dank im Voraus.

Answer 1

ich denke Problem darin: this-> Element (...)

ich weiß nicht genau, was für Fest tun. aber ich befehle, beide Matrix zu ändern. (A < -> B)

Answer 2

Wie in den Kommentaren erwähnt, in diesem Fall ist es besser, eine Dritt Überlastung der operator* zu verwenden (die mit der Matrix-Klasse be- friend ED).

Aber das ist nicht der Grund, warum es nicht funktioniert. Das Problem hier ist die Erklärung von

T temp; Matrix<T> c; 

in Ihrem in-class operator*. Sie sollten jede T mit std::common_type_t<T,U> ersetzen:

using C = std::common_type_t<T,U>; 
C temp; Matrix<C> c; 

Ansonsten ist das Element-Typ der zurückgegebenen Matrix wird die gleiche wie die der anrufenden Matrix immer sein. Das heißt, wenn Sie C = A*B aufrufen, die als C = A.operator*(B) aufgelöst wird, ist der Elementtyp C derselbe wie A.

Dann, wenn A ein Matrix<double> ist und Sie multiplizieren es mit einem Matrix<complex<double> >, würden Sie einen Matrix<double> als Rückkehr bekommen - was offensichtlich falsch ist.