fixed bug in scaling of rhs for cpr

This commit is contained in:
hnil 2019-03-18 15:34:49 +01:00 committed by Atgeirr Flø Rasmussen
parent 31f4c01938
commit b27e56e386

View File

@ -769,11 +769,11 @@ protected:
}
block[pressureEqnIndex] = neweq;
}
BlockVector newrhs(0.0);
Scalar newrhs(0.0);
for (std::size_t ii = 0; ii < brhs.size(); ii++) {
newrhs += bweights[ii]*brhs[ii];
}
brhs = newrhs;
brhs[pressureEqnIndex] = newrhs;
}
}