fix: extract_mesh (#110)

hpc-maths · Jul 11, 2023 · 3c2c05e · 3c2c05e
1 parent 7be28c4
commit 3c2c05e
Show file tree

Hide file tree

Showing 2 changed files with 10 additions and 0 deletions.
diff --git a/include/samurai/utils.hpp b/include/samurai/utils.hpp
@@ -23,6 +23,9 @@ namespace samurai
     template <class F, class... CT>
     class field_function;
 
+    template <template <class T> class OP, class... CT>
+    class field_operator_function;
+
     namespace detail
     {
 
@@ -151,6 +154,11 @@ namespace samurai
         {
         };
 
+        template <template <class T> class OP, class... CT>
+        struct is_field_function<field_operator_function<OP, CT...>> : std::true_type
+        {
+        };
+
         template <class... T>
         auto& extract_mesh(const std::tuple<T...>& t)
         {

diff --git a/test/test_utils.cpp b/test/test_utils.cpp
@@ -4,6 +4,7 @@
 
 #include <samurai/box.hpp>
 #include <samurai/field.hpp>
+#include <samurai/stencil_field.hpp>
 #include <samurai/uniform_mesh.hpp>
 #include <samurai/utils.hpp>
 
@@ -18,5 +19,6 @@ namespace samurai
         auto u = make_field<double, 1>("u", mesh);
 
         static_assert(detail::is_field_function<decltype(5 + u)>::value);
+        static_assert(detail::is_field_function<decltype(upwind(1, u))>::value);
     }
 }