模板矩阵 - 矩阵乘法c ++

时间:2016-04-01 23:21:51

标签: c++

我正在尝试使用模板进行Matrix-Matrix乘法,并且我不断收到以下错误。 (我试图乘以非方矩阵)

错误1错误C2593:'运营商*'含糊不清

任何人都可以就如何解决这个问题给我一些建议吗?

//Matrix.h
#pragma once
#include <iostream>
#include <vector>
using namespace std;

template<class T, int m, int n>
class Matrix;

template<class T, int m, int n, int l>
Matrix<T, m, n> operator*(const Matrix<T, m, n>&, const Matrix<T, n, l>&);

template<class T, int m, int n>
class Matrix
{
vector<vector<T>> elements;
int nrow;
int ncol;

public:
    Matrix();
    ~Matrix();
    void print();
    template<int l>
    friend Matrix<T, m, l> operator*<>(const Matrix<T, m, n>&, const Matrix<T, n, l>&);

};

template<class T, int m, int n>
Matrix<T, m, n>::Matrix() : nrow(m), ncol(n)
{
    for (int i = 0; i < nrow; i++){
        vector<T> row(ncol, i);
        elements.push_back(row);
    }
}

template<class T, int m, int n>
Matrix<T, m, n>::~Matrix(){}

template<class T, int m, int n>
void Matrix<T, m, n>::print()
{
    for (int i = 0; i < nrow; i++){
        for (int j = 0; j < ncol; j++){
            cout << elements[i][j] << " ";
        }
    cout << endl;
    }
}

template<class T, int m, int n, int l> 
Matrix<T, m, l> operator*(const Matrix<T, m, n>& m1, const Matrix<T, n, l>& m2){
    int nrow = m1.nrow;
    int ncol = m2.ncol;
    Matrix<T, m, l> m3;
    for (int i = 0; i < nrow; ++i){
        for (int j = 0; j < ncol; ++j){
            m3.elements[i][j] = 0;
            for (int k = 0; k < m1.ncol; k++){
                T temp = m1.elements[i][k] * m2.elements[k][j];
                m3.elements[i][j] = temp + m3.elements[i][j];
            }
        }
    }
return m3;
}

//main.cpp
#include "Matrix.h"
using namespace std;

int main()
{
Matrix<int, 3, 2> a;
Matrix<int, 2, 1> b;
Matrix<int, 3, 1> c;
c = a*b;

c.print();
}

问题出现在矩阵乘法中,可能是由于模板中的编码错误。

2 个答案:

答案 0 :(得分:1)

错误在于:

./matrix.cpp:48:28: error: function template partial specialization is not allowed
    friend Matrix<T, m, l> operator*<>(const Matrix<T, m, n>&, const Matrix<T, n, l>&);
                           ^        ~~
1 error generated.

将其更改为:

friend Matrix<T, m, l> operator*(const Matrix<T, m, n>&, const Matrix<T, n, l>&);

答案 1 :(得分:1)

我必须改变Richard改变的内容,但我还必须更改operator*friend的声明,如下所示:

template<class T, int m, int n, int l>
Matrix<T, m, l> operator*(const Matrix<T, m, n>&, const Matrix<T, n, l>&);
          // ^ here

template<class _T, int _m, int _n, int l>
friend Matrix<_T, _m, l> operator*(const Matrix<_T, _m, _n>&, const Matrix<_T, _n, l>&);

如果我没有更改这两个中的第一个,我得到了“operator*含糊不清”的错误,因为前向声明与进一步向下的实例化不匹配。

目前正在输出:

0
1
2

这似乎不太正确,但我还没有进一步调试。