2

我正在创建一个矩阵类（主要用于自我练习，我知道 Eigen）。在其中，我将矩阵乘法定义为人们所期望的 - 该方法采用对 Matrix 对象的 const 引用并返回一个 Matrix 对象。

Matrix Matrix::operator *(const Matrix& rhs) const
{
    Matrix temp(mRows,rhs.mCols);
    ... //do Matrix Multiplication
    return temp;
}

然后我派生了一个 DCM 类（方向余弦矩阵——只是一种特殊的矩阵）。 

class DCM: public Matrix
{
    .... // class definition here, but does NOT contain definition for
         // matrix multiplication    
};

我可以创建两个 DCM 对象并将它们相乘没有问题，但返回的对象是 Matrix 类型。

DCM Rvb_q(q);

DCM Rvb_p(p);

DCM Rvb_pq(Rvb_p*Rvb_q); // error because Rvb_p*Rvb_q returns a Matrix object

有没有办法让该函数返回一个 DCM 对象而不必在派生类中重新编码该函数？一种方法是将此构造函数添加到派生类：

    DCM(const Matrix &M):Matrix(M) {}

但这似乎真的很低效（为乘法创建一个新对象，然后在制作 DCM 对象时复制它）和限制（任何时候我将两个 DCM 对象相乘，我都必须为 Matrix 对象创建一个新的 DCM 对象复制到——我不能只使用返回的对象作为 DCM）。有没有办法重用基类函数，但让它返回派生类类型？

4

2 回答 2

2

您可以使用Curiously Recurrent Template Pattern来实现。CRTP 是一个习惯用法，其中一个类从模板类继承，并将其作为模板参数。例如：

template<typename T>
struct CRTP_base {};

struct foo : public CRTP_base<foo> {};

这个习语的强大之处在于基类知道从它派生什么类。您可以使用该属性来解决您的问题：

template<typename T>
struct matrix_base
{
    T& operator*=(const T& rhs)
    {
        /* Multiplication code */
        return static_cast<T&>(*this);
    }
};

struct DCM : public matrix_base<DCM> {};

int main()
{
    DCM a,b c;

    a *= b; //OK, call to matrix_base<DCM>::operator*, which expects a DCM as parameter,
            //and returns a reference to a DCM.
}
于 2013-10-05T11:08:16.567 回答
1

如果从矩阵继承的任何东西乘以相同的类型（在乘法下关闭），那么您可以将类似的东西组合在一起：

template <class TMatrix>
TMatrix Matrix::operator *(const TMatrix& rhs)
{
    TMatrix temp(mRows, rhs.mCols);
    // multiply
    return temp;
}

但我不建议这样做。在这种情况下，我认为最好明确将 DCM 乘以 DCM 会得到 DCM：

DCM DCM::operator *(const DCM& rhs)
{
    return Matrix::operator*(rhs);
}

在这种情况下，您还应该声明（可能是私有的）DCM(const Matrix&) and DCM(Matrix&&)。

于 2013-10-05T05:28:19.430 回答