OilfieldToolsNet/opm-common
4
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

fxing the warning related to compIdx in PTFlash

Browse Source
This commit is contained in:
Kai Bao 2022-06-27 11:19:26 +02:00
parent 41a489c9ad
commit 1728d53371
1 changed files with 13 additions and 13 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

12
opm/material/constraintsolvers/PTFlash.hpp
View File

@ -1048,7 +1048,7 @@ protected:
std::vector<double> deri(num_deri, 0.); std::vector<double> deri(num_deri, 0.);
// derivatives from P // derivatives from P
for (unsigned idx = 0; idx < num_deri; ++idx) { for (unsigned idx = 0; idx < num_deri; ++idx) {
deri[idx] = - sec_jac[compIdx][0] * p_l.derivative(idx); deri[idx] = -sec_jac[compIdx][0] * p_l.derivative(idx);
} }
for (unsigned cIdx = 0; cIdx < numComponents; ++cIdx) { for (unsigned cIdx = 0; cIdx < numComponents; ++cIdx) {
@ -1063,7 +1063,7 @@ protected:
} }
// handling y // handling y
for (unsigned idx = 0; idx < num_deri; ++idx) { for (unsigned idx = 0; idx < num_deri; ++idx) {
deri[idx] = - sec_jac[compIdx + numComponents][0]* p_v.derivative(idx); deri[idx] = -sec_jac[compIdx + numComponents][0] * p_v.derivative(idx);
} }
for (unsigned cIdx = 0; cIdx < numComponents; ++cIdx) { for (unsigned cIdx = 0; cIdx < numComponents; ++cIdx) {
const double pz = -sec_jac[compIdx + numComponents][cIdx + 1]; const double pz = -sec_jac[compIdx + numComponents][cIdx + 1];
@ -1079,10 +1079,10 @@ protected:
// handling derivatives of L // handling derivatives of L
std::vector<double> deriL(num_deri, 0.); std::vector<double> deriL(num_deri, 0.);
for (unsigned idx = 0; idx < num_deri; ++idx) { for (unsigned idx = 0; idx < num_deri; ++idx) {
deriL[idx] = - sec_jac[2*numComponents][0] * p_v.derivative(idx); deriL[idx] = -sec_jac[2 * numComponents][0] * p_v.derivative(idx);
} }
for (unsigned cIdx = 0; cIdx < numComponents; ++cIdx) { for (unsigned cIdx = 0; cIdx < numComponents; ++cIdx) {
const double pz = -sec_jac[2*numComponents][cIdx + 1]; const double pz = -sec_jac[2 * numComponents][cIdx + 1];
const auto& zi = z[cIdx]; const auto& zi = z[cIdx];
for (unsigned idx = 0; idx < num_deri; ++idx) { for (unsigned idx = 0; idx < num_deri; ++idx) {
deriL[idx] += pz * zi.derivative(idx); deriL[idx] += pz * zi.derivative(idx);
@ -1092,6 +1092,7 @@ protected:
for (unsigned idx = 0; idx < num_deri; ++idx) { for (unsigned idx = 0; idx < num_deri; ++idx) {
L_eval.setDerivative(idx, deriL[idx]); L_eval.setDerivative(idx, deriL[idx]);
} }
}
// set up the mole fractions // set up the mole fractions
for (unsigned compIdx = 0; compIdx < numComponents; ++compIdx) { for (unsigned compIdx = 0; compIdx < numComponents; ++compIdx) {
@ -1099,8 +1100,7 @@ protected:
fluid_state.setMoleFraction(FluidSystem::gasPhaseIdx, compIdx, y[compIdx]); fluid_state.setMoleFraction(FluidSystem::gasPhaseIdx, compIdx, y[compIdx]);
} }
fluid_state.setLvalue(L_eval); fluid_state.setLvalue(L_eval);
} } //end updateDerivatives
}//end updateDerivatives
// TODO: or use typename FlashFluidState::Scalar // TODO: or use typename FlashFluidState::Scalar
template <class FlashFluidState, class ComponentVector> template <class FlashFluidState, class ComponentVector>