| 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 |
|
} |
