Use GridKit::Testing::isEqual in Enzyme examples.

Author: nkoukpaizan

examples/Enzyme/Library/Scalar/CMakeLists.txt

@@ -1,4 +1,4 @@
 add_executable(EnzymeLibScalarCheck EnzymeScalar.cpp ScalarModel.cpp)
-target_link_libraries(EnzymeLibScalarCheck ClangEnzymeFlags)
+target_link_libraries(EnzymeLibScalarCheck ClangEnzymeFlags GRIDKIT::Utilities)
 
 add_test(NAME "EnzymeLibScalarCheck" COMMAND ${CMAKE_CURRENT_BINARY_DIR}/EnzymeLibScalarCheck)

examples/Enzyme/Library/Scalar/EnzymeScalar.cpp

@@ -1,7 +1,7 @@
 #include <iostream>
-#include <limits>
 
 #include "ScalarModel.hpp"
+#include <Utilities/Testing.hpp>
 
 /**
  * @brief Example that computes the derivative of a library function
@@ -23,7 +23,7 @@ int main()
   double dsq = scalar_model.getDerivativeValue();
 
   std::cout << "x = " << var << ", x^2 = " << sq << ", d(x^2)/dx = " << dsq << "\n";
-  if (std::abs(dsq - 2.0 * var) > std::numeric_limits<double>::epsilon())
+  if (!GridKit::Testing::isEqual(dsq, 2.0 * var))
   {
     fail++;
     std::cout << "Result incorrect\n";

examples/Enzyme/Library/Vector/CMakeLists.txt

@@ -1,4 +1,4 @@
 add_executable(EnzymeLibVectorCheck EnzymeVector.cpp VectorModel.cpp)
-target_link_libraries(EnzymeLibVectorCheck ClangEnzymeFlags GRIDKIT::DenseMatrix)
+target_link_libraries(EnzymeLibVectorCheck ClangEnzymeFlags GRIDKIT::DenseMatrix GRIDKIT::Utilities)
 
 add_test(NAME "EnzymeLibVectorCheck" COMMAND ${CMAKE_CURRENT_BINARY_DIR}/EnzymeLibVectorCheck)

examples/Enzyme/Library/Vector/EnzymeVector.cpp

@@ -1,7 +1,7 @@
 #include <iostream>
-#include <limits>
 
 #include "VectorModel.hpp"
+#include <Utilities/Testing.hpp>
 
 /**
  * @brief Example that computes the Jacobian of a vector-valued residual
@@ -63,7 +63,7 @@ int main()
   {
     for (size_t idx = 0; idx < var.size(); ++idx)
     {
-      if (std::abs(jac.getValue(idx, idy) - jac_ref.getValue(idx, idy)) > std::numeric_limits<double>::epsilon())
+      if (!GridKit::Testing::isEqual(jac.getValue(idx, idy), jac_ref.getValue(idx, idy)))
       {
         fail++;
         if (verbose)

examples/Enzyme/PowerElectronics/CMakeLists.txt

@@ -1,5 +1,5 @@
 add_executable(EnzymePowerElectronicsCheck main.cpp)
 target_compile_options(EnzymePowerElectronicsCheck PUBLIC -fno-vectorize -ffast-math -fno-unroll-loops)
-target_link_libraries(EnzymePowerElectronicsCheck ClangEnzymeFlags GRIDKIT::DenseMatrix GRIDKIT::power_elec_disgen)
+target_link_libraries(EnzymePowerElectronicsCheck ClangEnzymeFlags GRIDKIT::DenseMatrix GRIDKIT::power_elec_disgen GRIDKIT::Utilities)
 
 add_test(NAME "EnzymePowerElectronicsCheck" COMMAND ${CMAKE_CURRENT_BINARY_DIR}/EnzymePowerElectronicsCheck)

examples/Enzyme/PowerElectronics/main.cpp

@@ -1,9 +1,9 @@
 #include <iostream>
-#include <limits>
 
 #include <LinearAlgebra/DenseMatrix/DenseMatrix.hpp>
 #include <Model/PowerElectronics/DistributedGenerator/DistributedGenerator.hpp>
 #include <Model/PowerElectronics/SystemModelPowerElectronics.hpp>
+#include <Utilities/Testing.hpp>
 
 /**
  * @brief Standalone example that computes the Jacobian associated with the
@@ -177,7 +177,7 @@ int main()
     {
       double jac_value     = jac_autodiff.getValue(idx, idy);
       double jac_ref_value = jac_ref_dense.getValue(idx, idy);
-      if (std::abs(jac_value - jac_ref_value) > 10 * std::numeric_limits<double>::epsilon())
+      if (!GridKit::Testing::isEqual(jac_value, jac_ref_value))
       {
         fail++;
         if (verbose)

examples/Enzyme/Standalone/CMakeLists.txt

@@ -1,12 +1,12 @@
 add_executable(EnzymeStandaloneScalarCheck EnzymeScalar.cpp)
-target_link_libraries(EnzymeStandaloneScalarCheck ClangEnzymeFlags)
+target_link_libraries(EnzymeStandaloneScalarCheck ClangEnzymeFlags GRIDKIT::Utilities)
 
 add_executable(EnzymeStandaloneVectorCheck EnzymeVector.cpp)
-target_link_libraries(EnzymeStandaloneVectorCheck ClangEnzymeFlags GRIDKIT::DenseMatrix)
+target_link_libraries(EnzymeStandaloneVectorCheck ClangEnzymeFlags GRIDKIT::DenseMatrix GRIDKIT::Utilities)
 
 add_executable(EnzymeStandaloneSparseCheck EnzymeSparse.cpp)
 target_compile_options(EnzymeStandaloneSparseCheck PUBLIC -mllvm -enzyme-auto-sparsity=1)
-target_link_libraries(EnzymeStandaloneSparseCheck ClangEnzymeFlags GRIDKIT::SparseMatrix)
+target_link_libraries(EnzymeStandaloneSparseCheck ClangEnzymeFlags GRIDKIT::SparseMatrix GRIDKIT::Utilities)
 
 add_test(NAME "EnzymeStandaloneScalarCheck" COMMAND ${CMAKE_CURRENT_BINARY_DIR}/EnzymeStandaloneScalarCheck)
 add_test(NAME "EnzymeStandaloneVectorCheck" COMMAND ${CMAKE_CURRENT_BINARY_DIR}/EnzymeStandaloneVectorCheck)

examples/Enzyme/Standalone/EnzymeScalar.cpp

@@ -1,5 +1,5 @@
 #include <iostream>
-#include <limits>
+#include <Utilities/Testing.hpp>
 
 /**
  * @brief Standalone example that computes the derivative of a scalar function
@@ -27,7 +27,7 @@ int main()
   double sq   = square(var);
   double dsq  = dsquare(var);
   std::cout << "x = " << var << ", x^2 = " << sq << ", d(x^2)/dx = " << dsq << "\n";
-  if (std::abs(dsq - 2.0 * var) > std::numeric_limits<double>::epsilon())
+  if (!GridKit::Testing::isEqual(dsq, 2.0 * var))
   {
     fail++;
     std::cout << "Result incorrect\n";

examples/Enzyme/Standalone/EnzymeSparse.cpp

@@ -1,10 +1,10 @@
 #include <assert.h>
 #include <cmath>
-#include <limits>
 #include <stdio.h>
 #include <vector>
 
 #include <LinearAlgebra/SparseMatrix/COO_Matrix.hpp>
+#include <Utilities/Testing.hpp>
 
 /**
  * @brief Standalone example that computes the sparse Jacobian of a vector-valued function
@@ -165,7 +165,7 @@ void check(SparseMatrix matrix_1, SparseMatrix matrix_2, int& fail)
       fail++;
     if (ccord_1[ind] != ccord_2[ind])
       fail++;
-    if (std::abs(vals_1[ind] - vals_2[ind]) > std::numeric_limits<double>::epsilon())
+    if (!GridKit::Testing::isEqual(vals_1[ind], vals_2[ind]))
       fail++;
   }
 }

examples/Enzyme/Standalone/EnzymeVector.cpp

@@ -1,8 +1,8 @@
 #include <iostream>
-#include <limits>
 #include <vector>
 
 #include <LinearAlgebra/DenseMatrix/DenseMatrix.hpp>
+#include <Utilities/Testing.hpp>
 
 /**
  * @brief Standalone example that computes the Jacobian of a vector-valued function
@@ -110,7 +110,7 @@ int main()
   {
     for (size_t idx = 0; idx < x.size(); ++idx)
     {
-      if (std::abs(dsq.getValue(idx, idy) - dsq_ref.getValue(idx, idy)) > std::numeric_limits<double>::epsilon())
+      if (!GridKit::Testing::isEqual(dsq.getValue(idx, idy), dsq_ref.getValue(idx, idy)))
       {
         fail++;
         if (verbose)