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

Use a custom handle for copying doubles for above.

Browse Source 
The one previously used is unfortunately protected upstream (not noticed with gcc <=8).
This commit is contained in:
Markus Blatt 2020-12-07 11:44:10 +01:00
parent 0269a36865
commit 932e3c8fb9
1 changed files with 15 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

17
opm/simulators/wells/ParallelWellInfo.cpp
View File

@ -90,6 +90,20 @@ void CommunicateAbove::endReset()
#endif
}
struct CopyGatherScatter
{
static const double& gather(const double* a, std::size_t i)
{
return a[i];
}
static void scatter(double* a, const double& v, std::size_t i)
{
a[i] = v;
}
};
std::vector<double> CommunicateAbove::communicate(double first_above,
const double* current,
std::size_t size)
@ -99,13 +113,12 @@ std::vector<double> CommunicateAbove::communicate(double first_above,
#if HAVE_MPI
if (comm_.size() > 1)
{
using Handle = Dune::OwnerOverlapCopyCommunication<int,int>::CopyGatherScatter<double*>;
auto aboveData = above.data();
// Ugly const_cast needed as my compiler says, that
// passing const double*& and double* as parameter is
// incompatible with function decl template<Data> forward(const Data&, Data&))
// That would need the first argument to be double* const&
communicator_.forward<Handle>(const_cast<double*>(current), aboveData);
communicator_.forward<CopyGatherScatter>(const_cast<double*>(current), aboveData);
}
else
#endif