diff --git a/src/math.f90 b/src/math.f90 index c6a358407..863661c2c 100644 --- a/src/math.f90 +++ b/src/math.f90 @@ -1380,10 +1380,10 @@ subroutine math_selfTest() if (any(dNeq0(math_eye(3),math_inv33(math_I3)))) & error stop 'math_inv33(math_I3)' - if (any(dNeq(t33,math_symmetric33(t33)+math_skew33(t33)))) & + if (any(dNeq(t33,math_symmetric33(t33)+math_skew33(t33),1.0e-10_pReal))) & error stop 'math_symmetric/skew' - if (any(dNeq(t33,math_spherical33(t33)+math_deviatoric33(t33)))) & + if (any(dNeq(t33,math_spherical33(t33)+math_deviatoric33(t33),1.0e-10_pReal))) & error stop 'math_spherical/deviatoric' do while(abs(math_det33(t33))<1.0e-9_pREAL)