-
Notifications
You must be signed in to change notification settings - Fork 62
/
MatrixOpe.cpp
101 lines (101 loc) · 2.59 KB
/
MatrixOpe.cpp
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
#include <MatrixOpe.h>
#include <stdlib.h>
namespace MLL
{
Matrix operator+(const Matrix &matrix1, const Matrix &matrix2)
{
if(matrix1._col != matrix2._col || matrix1._row != matrix2._row)
{
std::cout << matrix1._row << "*" << matrix1._col << std::endl;
std::cout << matrix2._row << "*" << matrix2._col << std::endl;
std::cout << "addData data1 data2 is no" << std::endl;
exit(-1);
}
size_t i = 0, j = 0;
Matrix add(matrix1._row, matrix1._col, 0);
for(i = 0; i < matrix1._row; i++)
{
for(j = 0; j < matrix1._col; j++)
{
add._data[i][j] = matrix1._data[i][j] + matrix2._data[i][j];
}
}
return add;
}
Matrix operator-(const Matrix &matrix1, const Matrix &matrix2)
{
if(matrix1._col != matrix2._col || matrix1._row != matrix2._row)
{
std::cout << matrix1._row<< "*" << matrix1._col << std::endl;
std::cout << matrix2._row<< "*" << matrix2._col << std::endl;
std::cout << "addData data1 data2 is no" << std::endl;
exit(-1);
}
size_t i = 0, j = 0;
Matrix sub(matrix1._row, matrix1._col, 0);
for(i = 0; i < matrix1._row; i++)
{
for(j = 0; j < matrix1._col; j++)
{
sub._data[i][j] = matrix1._data[i][j] - matrix2._data[i][j];
}
}
return sub;
}
Matrix operator*(const Matrix &matrix1, const Matrix &matrix2)
{
if(matrix1._col != matrix2._row)
{
std::cout << matrix1._row << "*" << matrix1._col << std::endl;
std::cout << matrix2._row << "*" << matrix2._col << std::endl;
std::cout << "multsData error" << std::endl;
exit(-1);
}
size_t i = 0, j = 0, k = 0;
Matrix mults(matrix1._row,matrix2._col,0);
for(i = 0; i < matrix1._row; i++)
{
for(j = 0; j < matrix2._col; j++)
{
mults._data[i][j] = 0;
}
}
for(i = 0; i < matrix1._row; i++)
{
for(j = 0; j < matrix2._col; j++)
{
for(k = 0; k < matrix1._col; k++)
{
mults._data[i][j] += matrix1._data[i][k] * matrix2._data[k][j];
}
}
}
return mults;
}
Matrix operator*(const double &alpha, const Matrix &matrix1)
{
size_t i = 0, j = 0;
Matrix mults(matrix1._row,matrix1._col,0);
for(i = 0; i < matrix1._row; i++)
{
for(j = 0; j < matrix1._col; j++)
{
mults._data[i][j] = alpha * matrix1._data[i][j];
}
}
return mults;
}
Matrix operator/(const double &alpha, const Matrix &matrix1)
{
size_t i = 0, j = 0;
Matrix mults(matrix1._row,matrix1._col,0);
for(i = 0 ; i < matrix1._row; i++)
{
for(j = 0; j < matrix1._col; j++)
{
mults._data[i][j] = matrix1._data[i][j] / alpha;
}
}
return mults;
}
}