Merge pull request #4655 from akva2/fix_stdwell_mobilility_scalar

fixed: copy the mobility values into the temporary Eval vector
This commit is contained in:
Kai Bao 2023-05-15 15:34:29 +02:00 committed by GitHub
commit 6021388045
No known key found for this signature in database
GPG Key ID: 4AEE18F83AFDEB23

View File

@ -955,6 +955,8 @@ namespace Opm
// as a result, the polymer and water share the same viscosity
if constexpr (!Base::has_polymermw) {
std::vector<EvalWell> mob_eval(this->num_components_, {this->primary_variables_.numWellEq() + Indices::numEq, 0.});
for (size_t i = 0; i < mob.size(); ++i)
mob_eval[i].setValue(mob[i]);
updateWaterMobilityWithPolymer(ebosSimulator, perf, mob_eval, deferred_logger);
for (size_t i = 0; i < mob.size(); ++i) {
mob[i] = getValue(mob_eval[i]);