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

avoid potentially uninitialized warnings in serial

Browse Source 
this was actually broken by an upstream change from return-by-reference
to return-by-value, but we were lucky, the return value was
kept in memory until end of block it seems.
This commit is contained in:
Arne Morten Kvarving 2020-03-19 10:03:44 +01:00
parent 563e1d9fbc
commit 42a2d22d89
1 changed files with 5 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

8
ebos/ecltransmissibility.hh
View File

@ -161,13 +161,15 @@ public:
// compute the axis specific "centroids" used for the transmissibilities. for
// consistency with the flow simulator, we use the element centers as
// computed by opm-parser's Opm::EclipseGrid class for all axes.
const double* centroid;
std::array<double, 3> centroid;
if (vanguard_.gridView().comm().rank() == 0) {
const auto& eclGrid = eclState.getInputGrid();
unsigned cartesianCellIdx = cartMapper.cartesianIndex(elemIdx);
centroid = &eclGrid.getCellCenter(cartesianCellIdx)[0];
centroid = eclGrid.getCellCenter(cartesianCellIdx);
} else
centroid = &centroids[centroidIdx * dimWorld];
std::copy(centroids.begin() + centroidIdx * dimWorld,
centroids.begin() + (centroidIdx + 1) * dimWorld,
centroid.begin());
for (unsigned axisIdx = 0; axisIdx < dimWorld; ++axisIdx)
for (unsigned dimIdx = 0; dimIdx < dimWorld; ++dimIdx)