38 |
|
#include "Vector.hpp" |
39 |
|
|
40 |
|
namespace oopse { |
41 |
– |
const double epsilon = 0.000001; |
42 |
– |
|
43 |
– |
template<typename T> |
44 |
– |
inline bool equal(T e1, T e2) { |
45 |
– |
return e1 == e2; |
46 |
– |
} |
47 |
– |
|
48 |
– |
template<> |
49 |
– |
inline bool equal(float e1, float e2) { |
50 |
– |
return fabs(e1 - e2) < epsilon; |
51 |
– |
} |
52 |
– |
|
53 |
– |
template<> |
54 |
– |
inline bool equal(double e1, double e2) { |
55 |
– |
return fabs(e1 - e2) < epsilon; |
56 |
– |
} |
41 |
|
|
42 |
|
/** |
43 |
|
* @class RectMatrix RectMatrix.hpp "math/RectMatrix.hpp" |
412 |
|
for (unsigned int i = 0; i < Row; i++) |
413 |
|
for (unsigned int j = 0; j < Col; j++) |
414 |
|
for (unsigned int k = 0; k < SameDim; k++) |
415 |
< |
result(i, j) = m1(i, k) * m2(k, j); |
415 |
> |
result(i, j) += m1(i, k) * m2(k, j); |
416 |
|
|
417 |
|
return result; |
418 |
|
} |