diff --git a/dumux/porousmediumflow/velocityoutput.hh b/dumux/porousmediumflow/velocityoutput.hh
index b6e4205e16d0e29b432b128e4e572f4fabb305f1..b29aa21021a5f66aed10d21d12df3670d819b75c 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())