OilfieldToolsNet/opm-simulators
4
0
Fork 0
mirror of https://github.com/OPM/opm-simulators.git synced 2025-02-25 18:55:30 -06:00
Code Issues Packages Projects Releases Wiki Activity

Transmissibility::update: use ThreadSafeMapBuilder for diffusivity_

Browse Source 
this in preparation for multithreading the loop
This commit is contained in:
Arne Morten Kvarving
2025-01-16 13:31:10 +01:00
parent 23e1654ba7
commit a502a86dc1
1 changed files with 6 additions and 3 deletions
Download Patch File Download Diff File Expand all files Collapse all files

9
opm/simulators/flow/Transmissibility_impl.hpp
View File

@@ -292,6 +292,8 @@ update(bool global, const TransUpdateQuantities update_quantities,
MapBuilderInsertionMode::Insert_Or_Assign);
ThreadSafeMapBuilder thermalHalfTrans(thermalHalfTrans_, 1,
MapBuilderInsertionMode::Insert_Or_Assign);
ThreadSafeMapBuilder diffusivity(diffusivity_, 1,
MapBuilderInsertionMode::Insert_Or_Assign);
// compute the transmissibilities for all intersections
for (const auto& elem : elements(gridView_)) {
@@ -411,7 +413,7 @@ update(bool global, const TransUpdateQuantities update_quantities,
}
if (updateDiffusivity && !onlyTrans) {
diffusivity_.insert_or_assign(details::isId(inside.elemIdx, outside.elemIdx), 0.0);
diffusivity.insert_or_assign(details::isId(inside.elemIdx, outside.elemIdx), 0.0);
}
if (updateDispersivity && !onlyTrans) {
dispersivity_.insert_or_assign(details::isId(inside.elemIdx, outside.elemIdx), 0.0);
@@ -479,8 +481,8 @@ update(bool global, const TransUpdateQuantities update_quantities,
// update the "diffusive half transmissibility" for the intersection
if (updateDiffusivity && !onlyTrans) {
diffusivity_.insert_or_assign(details::isId(inside.elemIdx, outside.elemIdx),
computeHalfMean(halfDiff, porosity_));
diffusivity.insert_or_assign(details::isId(inside.elemIdx, outside.elemIdx),
computeHalfMean(halfDiff, porosity_));
}
// update the "dispersivity half transmissibility" for the intersection
@@ -496,6 +498,7 @@ update(bool global, const TransUpdateQuantities update_quantities,
transMap.finalize();
thermalHalfTransBoundary.finalize();
thermalHalfTrans.finalize();
diffusivity.finalize();
// Potentially overwrite and/or modify transmissibilities based on input from deck
this->updateFromEclState_(global);