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

avoiding creating dx_V unnecessarily

Browse Source 
by using dx.array() in function solveWellEq().
This commit is contained in:
Kai Bao 2015-06-16 14:32:47 +02:00
parent 919286d700
commit f2089f5a1b
1 changed files with 1 additions and 4 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

5
opm/autodiff/BlackoilModelBase_impl.hpp
View File

@ -1266,10 +1266,7 @@ namespace detail {
const std::vector<M>& Jn = total_residual.derivative();
const Eigen::SparseLU< M > solver(Jn[0]);
const Eigen::VectorXd& dx = solver.solve(total_residual.value().matrix());
const int numeq = well_state.numWells()*(well_state.numPhases()+1);
V dx_V = V(numeq);
std::copy_n(dx.data(), numeq, dx_V.data());
updateWellState(dx_V, well_state);
updateWellState(dx.array(), well_state);
updateWellControls(well_state);
} while (it < 15);