Fixes wrong indexing for residual_.material_balance_eq.

Curiously, this indexes differently from  other stuff. It
uses the phase index and not the index within the phases
that are present.
This commit is contained in:
Markus Blatt 2015-01-28 19:11:13 +01:00
parent cb9048d928
commit 3fcb99c27a

View File

@ -1942,7 +1942,7 @@ namespace {
const int pos = pu.phase_pos[idx];
const ADB& tempB = rq_[pos].b;
B.col(pos) = 1./tempB.value();
R.col(pos) = residual_.material_balance_eq[pos].value();
R.col(pos) = residual_.material_balance_eq[idx].value();
tempV.col(pos) = R.col(pos).abs()/pv;
}
}