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

tried to get openmp work

Browse Source
This commit is contained in:
hnil 2022-06-21 09:16:26 +02:00 committed by Atgeirr Flø Rasmussen
parent df32014109
commit 08c0848e0d
3 changed files with 9 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

4
opm/models/blackoil/blackoilintensivequantitiessimple.hh
View File

@ -133,8 +133,8 @@ public:
void update(const Problem& problem,const PrimaryVariables& primaryVars,unsigned globalSpaceIdx, unsigned timeIdx)
{
//ParentType::update(elemCtx, dofIdx, timeIdx);//only used for extrusion factor
//const auto& materialParams = problem.materialLawParams(globalSpaceIdx);
const auto& materialParams = problem.materialLawParams(0);
const auto& materialParams = problem.materialLawParams(globalSpaceIdx);
//const auto& materialParams = problem.materialLawParams(0);//NB improve speed
Scalar RvMax;
if (FluidSystem::enableVaporizedOil()) {
RvMax = problem.maxOilVaporizationFactor(timeIdx, globalSpaceIdx);

4
opm/models/discretization/common/fvbasediscretization.hh
View File

@ -801,11 +801,15 @@ public:
{
//invalidateIntensiveQuantitiesCache(timeIdx);
size_t numGridDof = primaryVars.size();
#ifdef _OPENMP
#pragma omp parallel
#endif
for (unsigned dofIdx = 0; dofIdx < numGridDof; ++dofIdx) {
const auto& primaryVar = primaryVars[dofIdx];
auto& intquant = intensiveQuantityCache_[timeIdx][dofIdx];
intquant.update(problem, primaryVar, dofIdx, timeIdx);
}
std::fill(intensiveQuantityCacheUpToDate_[timeIdx].begin(),
intensiveQuantityCacheUpToDate_[timeIdx].end(),
/*value=*/true);

4
opm/models/discretization/common/linearizertpfa.hh
View File

@ -453,6 +453,9 @@ private:
const bool well_local = false;
resetSystem_();
unsigned numCells = model_().numTotalDof();
#ifdef _OPENMP
#pragma omp parallel
#endif
for(unsigned globI = 0; globI < numCells; globI++){
const auto& neighbours = neighbours_[globI];// this is a set but should maybe be changed
// accumulation term
@ -507,7 +510,6 @@ private:
setResAndJacobi(res, bMat, adres);
residual_[globI] += res;
jacobian_->addToBlock(globI, globI, bMat);
//residual_[globJ] += res;
bMat *= -1.0;
jacobian_->addToBlock(globJ, globI, bMat);
}