Avoid deleting derivative i.e. getting "wrong" matrix in case of zero flux

This commit is contained in:
hnil 2024-02-29 11:02:26 +01:00
parent 9b4c51de7e
commit ec11bbaac3

View File

@ -285,10 +285,6 @@ public:
distZ*g,
thpres,
problem.moduleParams());
if (pressureDifferences[phaseIdx] == 0) {
volumeFlux[phaseIdx] = 0.0;
continue;
}
const bool upwindIsInterior = (static_cast<unsigned>(upIdx[phaseIdx]) == interiorDofIdx);
const IntensiveQuantities& up = upwindIsInterior ? intQuantsIn : intQuantsEx;