From a488c2126bb2dba4667466ab414dd9f6d6bc86cc Mon Sep 17 00:00:00 2001
From: DennisGlaeser <dennis.glaeser@iws.uni-stuttgart.de>
Date: Thu, 25 Feb 2016 23:43:11 +0100
Subject: [PATCH] [LocalResidual] Adjust evalStorage method to new structure
---
dumux/implicit/localresidual.hh | 14 +++++++-------
1 file changed, 7 insertions(+), 7 deletions(-)
diff --git a/dumux/implicit/localresidual.hh b/dumux/implicit/localresidual.hh
index 2dbfb1b1c0..26882fb074 100644
--- a/dumux/implicit/localresidual.hh
+++ b/dumux/implicit/localresidual.hh
@@ -276,13 +276,13 @@ protected:
// calculate the amount of conservation each quantity inside
// all sub control volumes
- for (int scvIdx = 0; scvIdx < fvGeometry_().numScv; scvIdx++) {
- Valgrind::SetUndefined(storageTerm_[scvIdx]);
- asImp_().computeStorage(storageTerm_[scvIdx], scvIdx, /*isOldSol=*/false);
- storageTerm_[scvIdx] *=
- fvGeometry_().subContVol[scvIdx].volume
- * curVolVars_(scvIdx).extrusionFactor();
- Valgrind::CheckDefined(storageTerm_[scvIdx]);
+ for (auto&& scv : fvGeometry_().scvs())
+ {
+ int scvIdx = scv.indexInElement();
+
+ storageTerm_[scvIdx] = asImp_().computeStorage(scv, /*isOldSol=*/false);
+ storageTerm_[scvIdx] *= scv.volume();
+ storageTerm_[scvIdx] *= problem_().model().curVolVars(scv).extrusionFactor();
}
}
--
GitLab