28 const FVElementGeometry& fvGeometry)
31 using VelocityVector =
typename FVElementGeometry::GridGeometry::GlobalCoordinate;
32 VelocityVector result(0.0);
36 return std::find_if(vector.begin(), vector.end(), [](
const auto& x) { return std::abs(x) > 1e-8; } ) - vector.begin();
39 for (
const auto& scvf : scvfs(fvGeometry))
42 result[dirIdx] += 0.5*getFaceVelocity(fvGeometry, scvf)[dirIdx];