Skip to content
Snippets Groups Projects

[velocityoutput] fix missing cast to NumEqVector

Merged Bernd Flemisch requested to merge fix/missing-cast-to-numeqvector into master
1 file
+ 2
1
Compare changes
  • Side-by-side
  • Inline
@@ -305,7 +305,8 @@ public:
{
// check if we have Neumann no flow, we can just use 0
const auto neumannFlux = problem_.neumann(element, fvGeometry, elemVolVars, scvf);
if (Dune::FloatCmp::eq<std::decay_t<decltype(neumannFlux)>, Dune::FloatCmp::CmpStyle::absolute>(neumannFlux, 0.0, 1e-30))
using NumEqVector = std::decay_t<decltype(neumannFlux)>;
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())
Loading