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

Added a parallel version for computing the global reductions.

Browse Source
This commit is contained in:
Markus Blatt 2015-01-29 11:32:39 +01:00
parent 4527ce8ffd
commit ff9b8d790d
1 changed files with 40 additions and 8 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

48
opm/autodiff/FullyImplicitBlackoilSolver_impl.hpp
View File

@ -29,6 +29,7 @@
#include <opm/core/grid.h> #include <opm/core/grid.h>
#include <opm/core/linalg/LinearSolverInterface.hpp> #include <opm/core/linalg/LinearSolverInterface.hpp>
#include <opm/core/linalg/ParallelIstlInformation.hpp>
#include <opm/core/props/rock/RockCompressibility.hpp> #include <opm/core/props/rock/RockCompressibility.hpp>
#include <opm/core/simulator/BlackoilState.hpp> #include <opm/core/simulator/BlackoilState.hpp>
#include <opm/core/utility/ErrorMacros.hpp> #include <opm/core/utility/ErrorMacros.hpp>
@ -1898,16 +1899,47 @@ namespace {
int nc) const int nc) const
{ {
// Do the global reductions // Do the global reductions
for ( int idx=0; idx<MaxNumPhases; ++idx ) #if HAVE_MPI
if(linsolver_.parallelInformation().type()==typeid(ParallelISTLInformation))
{ {
if (active_[idx]) { const ParallelISTLInformation& info =
B_avg[idx] = B.col(idx).sum()/nc; boost::any_cast<const ParallelISTLInformation&>(linsolver_.parallelInformation());
maxCoeff[idx]=tempV.col(idx).maxCoeff(); for ( int idx=0; idx<MaxNumPhases; ++idx )
R_sum[idx] = R.col(idx).sum();
}
else
{ {
R_sum[idx] = B_avg[idx] = maxCoeff[idx] =0.; if (active_[idx]) {
auto values = std::tuple<double,double,double>(0.,0.,0.);
auto containers = std::make_tuple(B.col(idx),
tempV.col(idx),
R.col(idx));
auto operators = std::make_tuple(Opm::Reduction::makeGlobalSumFunctor<double>(),
Opm::Reduction::makeGlobalMaxFunctor<double>(),
Opm::Reduction::makeGlobalSumFunctor<double>());
nc=info.communicator().sum(nc);
info.computeReduction(containers,operators,values);
B_avg[idx] = std::get<0>(values)/nc;
maxCoeff[idx] = std::get<1>(values);
R_sum[idx] = std::get<2>(values);
}
else
{
R_sum[idx] = B_avg[idx] = maxCoeff[idx] = 0.;
}
}
}
else
#endif
{
for ( int idx=0; idx<MaxNumPhases; ++idx )
{
if (active_[idx]) {
B_avg[idx] = B.col(idx).sum()/nc;
maxCoeff[idx]=tempV.col(idx).maxCoeff();
R_sum[idx] = R.col(idx).sum();
}
else
{
R_sum[idx] = B_avg[idx] = maxCoeff[idx] =0.;
}
} }
} }
// Compute total pore volume // Compute total pore volume