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

Merge pull request #2966 from blattms/use-custom-handle

Browse Source 
Use a custom handle for copying doubles for above.
This commit is contained in:
Atgeirr Flø Rasmussen
2020-12-07 12:05:26 +01:00
committed by GitHub
parent c164762763 932e3c8fb9
commit 5477588a18
1 changed files with 15 additions and 2 deletions
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