From 48db4062b58b65024b65158c50add45f4218d60d Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?Atgeirr=20Fl=C3=B8=20Rasmussen?= Date: Fri, 23 Nov 2018 21:37:52 +0100 Subject: [PATCH] Count lower/upper entries simpler and faster. --- opm/autodiff/ParallelOverlappingILU0.hpp | 10 +++++----- 1 file changed, 5 insertions(+), 5 deletions(-) diff --git a/opm/autodiff/ParallelOverlappingILU0.hpp b/opm/autodiff/ParallelOverlappingILU0.hpp index ce44781b5..9b63dbba4 100644 --- a/opm/autodiff/ParallelOverlappingILU0.hpp +++ b/opm/autodiff/ParallelOverlappingILU0.hpp @@ -446,14 +446,14 @@ namespace Opm size_type numLower = 0; size_type numUpper = 0; const auto endi = A.end(); - for (auto i=A.begin(); i!=endi; ++i) { + for (auto i = A.begin(); i != endi; ++i) { const size_type iIndex = i.index(); + size_type numLowerRow = 0; for (auto j = (*i).begin(); j.index() < iIndex; ++j) { - ++numLower; - } - for (auto j = (*i).beforeEnd(); j.index() > iIndex; --j) { - ++numUpper; + ++numLowerRow; } + numLower += numLowerRow; + numUpper += (*i).size() - numLowerRow - 1; } assert(numLower + numUpper + A.N() == A.nonzeroes());