--- trunk/OOPSE-2.0/src/math/RectMatrix.hpp 2004/10/13 22:24:59 1564 +++ trunk/OOPSE-2.0/src/math/RectMatrix.hpp 2004/10/13 23:53:40 1567 @@ -34,9 +34,11 @@ #ifndef MATH_RECTMATRIX_HPP #define MATH_RECTMATRIX_HPP +#include #include "Vector.hpp" namespace oopse { + const double epsilon = 0.000001; template inline bool equal(T e1, T e2) { @@ -45,12 +47,12 @@ namespace oopse { template<> inline bool equal(float e1, float e2) { - return e1 == e2; + return fabs(e1 - e2) < epsilon; } template<> inline bool equal(double e1, double e2) { - return e1 == e2; + return fabs(e1 - e2) < epsilon; } /**