| 78 | 
  | 
    SMat3x3::jacobi(a, w1L, v); | 
| 79 | 
  | 
 | 
| 80 | 
  | 
    CPPUNIT_ASSERT_DOUBLES_EQUAL(w1L[0], w1R[0], 0.0001); | 
| 81 | 
< | 
    CPPUNIT_ASSERT_DOUBLES_EQUAL(w1L[1], w1R[1], oopse::epsilon); | 
| 82 | 
< | 
    CPPUNIT_ASSERT_DOUBLES_EQUAL(w1L[2], w1R[2], oopse::epsilon); | 
| 81 | 
> | 
    CPPUNIT_ASSERT_DOUBLES_EQUAL(w1L[1], w1R[1], OpenMD::NumericConstant::epsilon); | 
| 82 | 
> | 
    CPPUNIT_ASSERT_DOUBLES_EQUAL(w1L[2], w1R[2], OpenMD::NumericConstant::epsilon); | 
| 83 | 
  | 
 | 
| 84 | 
  | 
} | 
| 85 | 
  | 
 | 
| 86 | 
  | 
void SquareMatrixTestCase::testTrace() { | 
| 87 | 
< | 
    CPPUNIT_ASSERT_DOUBLES_EQUAL(identMat.trace(), 3.0, oopse::epsilon); | 
| 88 | 
< | 
    CPPUNIT_ASSERT_DOUBLES_EQUAL(symMat.trace(), 4.0, oopse::epsilon); | 
| 87 | 
> | 
    CPPUNIT_ASSERT_DOUBLES_EQUAL(identMat.trace(), 3.0, OpenMD::NumericConstant::epsilon); | 
| 88 | 
> | 
    CPPUNIT_ASSERT_DOUBLES_EQUAL(symMat.trace(), 4.0, OpenMD::NumericConstant::epsilon); | 
| 89 | 
  | 
} | 
| 90 | 
  | 
void SquareMatrixTestCase::testIsSymmertric(){ | 
| 91 | 
  | 
    CPPUNIT_ASSERT(identMat.isSymmetric()); |