From f407a1c8768470626939aabc8ae8dc613a03c864 Mon Sep 17 00:00:00 2001
From: Timo Koch <timo.koch@iws.uni-stuttgart.de>
Date: Tue, 27 Nov 2018 18:03:28 +0100
Subject: [PATCH] [bugfix] Fix typo and missing include

---
 dumux/porousmediumflow/velocityoutput.hh | 3 ++-
 1 file changed, 2 insertions(+), 1 deletion(-)

diff --git a/dumux/porousmediumflow/velocityoutput.hh b/dumux/porousmediumflow/velocityoutput.hh
index b6e4205e16..b29aa21021 100644
--- a/dumux/porousmediumflow/velocityoutput.hh
+++ b/dumux/porousmediumflow/velocityoutput.hh
@@ -24,6 +24,7 @@
 #ifndef DUMUX_POROUSMEDIUMFLOW_VELOCITYOUTPUT_HH
 #define DUMUX_POROUSMEDIUMFLOW_VELOCITYOUTPUT_HH
 
+#include <dune/common/float_cmp.hh>
 #include <dune/geometry/referenceelements.hh>
 
 #include <dumux/common/parameters.hh>
@@ -306,7 +307,7 @@ public:
                         // check if we have Neumann no flow, we can just use 0
                         const auto neumannFlux = problem_.neumann(element, fvGeometry, elemVolVars, scvf);
                         using NumEqVector = std::decay_t<decltype(neumannFlux)>;
-                        if (Dune::FloatCmp::eq<NumEqVector>, Dune::FloatCmp::CmpStyle::absolute>(neumannFlux, NumEqVector(0.0), 1e-30))
+                        if (Dune::FloatCmp::eq<NumEqVector, Dune::FloatCmp::CmpStyle::absolute>(neumannFlux, NumEqVector(0.0), 1e-30))
                             scvfFluxes[scvfIndexInInside[localScvfIdx]] = 0;
                         // cubes
                         else if (dim == 1 || geomType.isCube())
-- 
GitLab