Merge pull request #626 from totto82/updateIntensive

Make a method that invalidate and updates the intensive quantities
OPM · Sep 28, 2020 · f79c543 · f79c543
2 parents c89f660 + d662902
commit f79c543
Showing 1 changed file with 20 additions and 0 deletions.
diff --git a/opm/models/discretization/common/fvbasediscretization.hh b/opm/models/discretization/common/fvbasediscretization.hh
@@ -765,6 +765,26 @@ public:
         }
     }
 
+    void invalidateAndUpdateIntensiveQuantities(unsigned timeIdx) const
+    {
+        invalidateIntensiveQuantitiesCache(timeIdx);
+
+        // loop over all elements...
+        ThreadedEntityIterator<GridView, /*codim=*/0> threadedElemIt(gridView_);
+#ifdef _OPENMP
+#pragma omp parallel
+#endif
+        {
+            ElementContext elemCtx(simulator_);
+            ElementIterator elemIt = threadedElemIt.beginParallel();
+            for (; !threadedElemIt.isFinished(elemIt); elemIt = threadedElemIt.increment()) {
+                const Element& elem = *elemIt;
+                elemCtx.updatePrimaryStencil(elem);
+                elemCtx.updatePrimaryIntensiveQuantities(/*timeIdx=*/0);
+            }
+        }
+    }
+
     /*!
      * \brief Move the intensive quantities for a given time index to the back.
      *