C++ 矩阵加法与乘法的实现

#include <iostream>using namespace std;//将矩阵a和b相加得到ctemplate <class T>void matrixAdd(T **a, T **b, T **c, int numberOfRows, int numberOfColumns){for (int i = 0; i < numberOfRows; i++)for (int j = 0; j < numberOfColumns; j++)c[i][j] = a[i][j] + b[i][j];}//两个n*n矩阵的乘法;将矩阵a和b相乘得到ctemplate <class T>void squareMatrixMultiply(T **a, T **b, T **c, int n){for (int i = 0; i < n; i++){for (int j = 0; j < n; j++){T sum = 0;for (int k = 0; k < n; k++){sum += a[i][k] * b[k][j];}c[i][j] = sum;}}}//m*n矩阵与n*p矩阵的乘法template <class T>void matrixMultiply(T **a, T **b, T **c, int n,int p){for (int i = 0; i < m; i++){for (int j = 0; j < p; j++){for (int k = 0; k < n; k++){sum += a[i][k] * b[k][j];}a[i][j] = sum;}}}int main(){int arr1[3] = { 1, 2, 3 }; //定义一维数组int arr2[3] = { 4, 5, 6 };int arr3[3] = { 7, 8, 9 };int arr4[3] = { 10, 11, 12 };int *a[3] = { arr1, arr2, arr3 }; //指针数组int *b[3] = { arr1, arr2, arr3 };int *c[4] = { arr1, arr2, arr3, arr4 };int *result1[3];int *result2[3];int *result3[4];int **f = a;matrixAdd(a, b, result1, 3, 3);squareMatrixMultiply(a, b, result2, 3);matrixMultiply(a, b, result3, 3, 4);return 0;}