矩阵乘法的 template class 该怎么声明？

koolcoy

我的矩阵是这样定义的:

1. 
   
2. template <class Numeric, int rows, int cols>
   
3. class Matrix {
   
4.   ...
   
5. };
   

复制代码

如果要重载一个 * 运算符， 这个重载的运算符的声明该怎么写？ 这样写对吗？

1. 
   
2. template <class Numeric, int rows, int cols>
   
3. class Matrix {
   
4. public:
   
5.   template <int cols2>
   
6.   Matrix<Numeric, rows, cols2>& operator*(Matrix<Numeric, cols, cols2>&);
   
7. };
   

复制代码

[ 本帖最后由 koolcoy 于 2007-10-18 12:45 编辑 ]

koolcoy

我是这么写的， 但是有几个问题：
matrix.h:

1. 
   
2. #ifndef MATRIX_H__
   
3. #define MATRIX_H__
   
4. 
   
5. #include <iostream>
   
6. 
   
7. template <class Numeric, int rows, int cols>
   
8. class Matrix {
   
9. public:
   
10.         Numeric elements[rows][cols];
    
11.         Matrix(Numeric arr[rows][cols]) {
    
12.                 for (int i = 0; i < rows; ++i) {
    
13.                         for (int j = 0; j < cols; ++j) {
    
14.                                 elements[i][j] = arr[i][j];
    
15.                         }
    
16.                 }
    
17.         }
    
18. 
    
19.         template <int cols2>
    
20.         const Matrix<Numeric, rows, cols2>&
    
21.         operator*(Matrix<Numeric, cols, cols2>& m) {
    
22.                 Numeric arr[rows][cols2];
    
23.                 for (int i = 0; i < rows; ++i) {
    
24.                         for (int j = 0; j < cols2; ++j) {
    
25.                                 Numeric sum = 0;
    
26.                                 for (int k = 0; k < cols; ++k) {
    
27.                                         sum += elements[i][k] * m.elements[k][j];
    
28.                                 }
    
29.                                 arr[i][j] = sum;
    
30.                         }
    
31.                 }
    
32.                 return Matrix<Numeric, rows, cols2>(arr);
    
33.         }
    
34. 
    
35.         void dump() {
    
36.                 for (int i = 0; i < rows; ++i) {
    
37.                         for (int j = 0; j < cols; ++j) {
    
38.                                 std::cout << elements[i][j] << ' ';
    
39.                         }
    
40.                         std::cout << std::endl;
    
41.                 }
    
42.         }
    
43. };
    
44. 
    
45. #endif
    

复制代码

test.cpp

1. 
   
2. #include "matrix.h"
   
3. #include <iostream>
   
4. 
   
5. using namespace std;
   
6. 
   
7. int main() {
   
8.         int arr[2][3] = {{1,2,3},{4,5,6}};
   
9.         int arr2[3][2] = {{1,2},{3,4},{5,6}};
   
10.         Matrix<int, 2, 3> m1(arr);
    
11.         Matrix<int, 3, 2> m2(arr2);
    
12.         m1.dump();
    
13.         m2.dump();
    
14.         Matrix<int, 2, 2> m3 = m1 * m2;
    
15.         m3.dump();
    
16. }
    

复制代码

1. 如果我把elements这个成员声明成private的话， 代码编译不过去。
2. return Matrix<Numeric, rows, cols2>(arr); 这一行会报warning:
  matrix.h:31: warning: returning reference to temporary
为什么？ 怎么修改？

vfdff

#ifndef MATRIX_H__
#define MATRIX_H__

#include <iostream>

template <class Numeric, int rows, int cols>
class Matrix {
public:
        Numeric elements[rows][cols];
        Matrix(Numeric arr[rows][cols]) {
                for (int i = 0; i < rows; ++i) {
                        for (int j = 0; j < cols; ++j) {
                                elements[j] = arr[j];
                        }
                }
        }

        template <int cols2>
        const Matrix<Numeric, rows, cols2>&
        operator*(Matrix<Numeric, cols, cols2>& m) {
                Numeric arr[rows][cols2];
                for (int i = 0; i < rows; ++i) {
                        for (int j = 0; j < cols2; ++j) {
                                Numeric sum = 0;
                                for (int k = 0; k < cols; ++k) {
                                        sum += elements[k] * m.elements[k][j];
                                }
                                arr[j] = sum;
                        }
                }
                return Matrix<Numeric, rows, cols2>(arr);
        }

        void dump() {
                for (int i = 0; i < rows; ++i) {
                        for (int j = 0; j < cols; ++j) {
                                std::cout << elements[j] << ' ';
                        }
                        std::cout << std::endl;
                }
        }
};

#endif
上面这个算是声明还是定义呀 ？？