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

convertSats(): Transpose loop nest for locality of data

Browse Source
This commit is contained in:
Bård Skaflestad 2015-09-04 15:08:05 +02:00 committed by Andreas Lauser
parent 2cbf99fcde
commit b08f89e68a
1 changed files with 13 additions and 5 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

18
opm/core/simulator/initStateEquil_impl.hpp
View File

@ -780,14 +780,22 @@ namespace Opm
inline std::vector<double> inline std::vector<double>
convertSats(const std::vector< std::vector<double> >& sat) convertSats(const std::vector< std::vector<double> >& sat)
{ {
const int np = sat.size(); const auto np = sat.size();
const int nc = sat[0].size(); const auto nc = sat[0].size();
std::vector<double> s(np * nc); std::vector<double> s(np * nc);
for (int c = 0; c < nc; ++c) {
for (int p = 0; p < np; ++p) { for (decltype(sat.size()) p = 0; p < np; ++p) {
s[np*c + p] = sat[p][c]; const auto& sat_p = sat[p];
double* sp = & s[0*nc + p];
for (decltype(sat[0].size()) c = 0;
c < nc; ++c, sp += np)
{
*sp = sat_p[c];
} }
} }
return s; return s;
} }
} // namespace Details } // namespace Details