From ffefd5faa8aa73efc0e4c776c097f8d7dbaf3252 Mon Sep 17 00:00:00 2001
From: David Wells <drwells@email.unc.edu>
Date: Wed, 8 Jun 2022 12:10:34 -0400
Subject: [PATCH] Fix some CGAL tests.

1. Avoid an implicit muParser dependency
2. Fix compilation
---
 tests/cgal/cgal_surface_triangulation_03.cc | 10 ++++++++--
 tests/cgal/implicit_triangulation_01.cc     | 17 ++++++++++++++---
 tests/cgal/implicit_triangulation_02.cc     | 17 ++++++++++++++---
 3 files changed, 36 insertions(+), 8 deletions(-)

diff --git a/tests/cgal/cgal_surface_triangulation_03.cc b/tests/cgal/cgal_surface_triangulation_03.cc
index d5c2f4076e..98e9e348a3 100644
--- a/tests/cgal/cgal_surface_triangulation_03.cc
+++ b/tests/cgal/cgal_surface_triangulation_03.cc
@@ -47,7 +47,10 @@ main()
     std::ifstream input_sm(SOURCE_DIR "/input_grids/tripod.off");
     input_sm >> sm;
     cgal_surface_mesh_to_dealii_triangulation(sm, tria);
-    go.write_msh(tria, "temp.msh");
+    {
+      std::ofstream out("temp.msh");
+      go.write_msh(tria, out);
+    }
     remove("temp.msh");
     // If we got here, everything is fine.
     deallog << "OK" << std::endl;
@@ -61,7 +64,10 @@ main()
     std::ifstream input_p(SOURCE_DIR "/input_grids/tripod.off");
     input_p >> P;
     cgal_surface_mesh_to_dealii_triangulation(P, tria);
-    go.write_msh(tria, "temp.msh");
+    {
+      std::ofstream out("temp.msh");
+      go.write_msh(tria, out);
+    }
     remove("temp.msh");
     // If we got here, everything is fine.
     deallog << "OK" << std::endl;
diff --git a/tests/cgal/implicit_triangulation_01.cc b/tests/cgal/implicit_triangulation_01.cc
index 43c8201681..054fc51ea2 100644
--- a/tests/cgal/implicit_triangulation_01.cc
+++ b/tests/cgal/implicit_triangulation_01.cc
@@ -17,7 +17,7 @@
 
 #include <deal.II/base/config.h>
 
-#include <deal.II/base/function_parser.h>
+#include <deal.II/base/function.h>
 
 #include <deal.II/grid/grid_generator.h>
 #include <deal.II/grid/grid_out.h>
@@ -27,13 +27,24 @@
 
 using namespace CGAL::parameters;
 
+class ImplicitFunction : public Function<3>
+{
+public:
+  virtual double
+  value(const Point<3> &p, const unsigned int component = 0) const override
+  {
+    return std::pow(1 - std::sqrt(p[0] * p[0] + p[1] * p[1]), 2) + p[2] * p[2] -
+           .25;
+  }
+};
+
 int
 main()
 {
   initlog();
   // Build a deal.II triangulation
-  Triangulation<3>  tria;
-  FunctionParser<3> implicit_function("(1-sqrt(x^2+y^2))^2+z^2-.25");
+  Triangulation<3>                tria;
+  ImplicitFunction                implicit_function;
   CGALWrappers::AdditionalData<3> data;
   data.cell_size = 0.4;
   GridGenerator::implicit_function(
diff --git a/tests/cgal/implicit_triangulation_02.cc b/tests/cgal/implicit_triangulation_02.cc
index a0a76b0008..5e8456b7ff 100644
--- a/tests/cgal/implicit_triangulation_02.cc
+++ b/tests/cgal/implicit_triangulation_02.cc
@@ -17,7 +17,7 @@
 
 #include <deal.II/base/config.h>
 
-#include <deal.II/base/function_parser.h>
+#include <deal.II/base/function.h>
 
 #include <deal.II/grid/grid_generator.h>
 #include <deal.II/grid/grid_out.h>
@@ -27,13 +27,24 @@
 
 using namespace CGAL::parameters;
 
+class ImplicitFunction : public Function<3>
+{
+public:
+  virtual double
+  value(const Point<3> &p, const unsigned int component = 0) const override
+  {
+    return std::pow(1 - std::sqrt(p[0] * p[0] + p[1] * p[1]), 2) + p[2] * p[2] -
+           .25;
+  }
+};
+
 int
 main()
 {
   initlog();
   // Build a deal.II triangulation
-  Triangulation<2, 3> tria;
-  FunctionParser<3>   implicit_function("(1-sqrt(x^2+y^2))^2+z^2-.25");
+  Triangulation<2, 3>             tria;
+  ImplicitFunction                implicit_function;
   CGALWrappers::AdditionalData<2> data;
   data.angular_bound  = 30.;
   data.radius_bound   = .1;
-- 
2.39.5