--- trunk/OOPSE-2.0/src/math/SquareMatrix3.hpp	2005/09/27 20:02:57	2332
+++ trunk/OOPSE-2.0/src/math/SquareMatrix3.hpp	2005/09/28 16:55:57	2336
@@ -179,7 +179,7 @@ namespace oopse {
       Real ad1, ad2, ad3;    
       t = this->data_[0][0] + this->data_[1][1] + this->data_[2][2] + 1.0;
 
-      if( t > 0.0 ){
+      if( t > NumericConstant::epsilon ){
 
 	s = 0.5 / sqrt( t );
 	q[0] = 0.25 / s;