Replacing 1e-6 tolerance for f.s check with std::numeric_limits<double>::epsilon()

Author: aabrown100-git

tests/unitTests/material_model_tests/test_material_holzapfel_ogden.cpp

@@ -29,6 +29,7 @@
  */
 
 #include "test_material_holzapfel_ogden.h"
+#include <limits>
 
 
 
@@ -91,7 +92,7 @@ class HolzapfelOgdenTest : public :: MaterialTestFixture {
 
         // Check f.s = 0
         double dot_fs = params.f.dot(params.s);
-        if (fabs(dot_fs) > 1e-6) {
+        if (fabs(dot_fs) > std::numeric_limits<double>::epsilon()) {
             std::cout << "f.s = " << dot_fs << std::endl;
             std::cout << "f = [" << params.f[0] << ", " << params.f[1] << ", " << params.f[2] << "]" << std::endl;
             std::cout << "s = [" << params.s[0] << ", " << params.s[1] << ", " << params.s[2] << "]" << std::endl;

tests/unitTests/material_model_tests/test_material_holzapfel_ogden_MA.cpp

@@ -29,6 +29,7 @@
  */
 
 #include "test_material_holzapfel_ogden_MA.h"
+#include <limits>
 
 // ----------------------------------------------------------------------------
 // ------------- Holzapfel-Ogden (Modified Anisotropy) Material  --------------
@@ -87,7 +88,7 @@ class HolzapfelOgdenMATest : public MaterialTestFixture {
 
         // Check f.s = 0
         double dot_fs = params.f.dot(params.s);
-        if (fabs(dot_fs) > 1e-6) {
+        if (fabs(dot_fs) > std::numeric_limits<double>::epsilon()) {
             std::cout << "f.s = " << dot_fs << std::endl;
             std::cout << "f = [" << params.f[0] << ", " << params.f[1] << ", " << params.f[2] << "]" << std::endl;
             std::cout << "s = [" << params.s[0] << ", " << params.s[1] << ", " << params.s[2] << "]" << std::endl;