From d89cbfdda6ddb48b054dc3a7120d17936eee5865 Mon Sep 17 00:00:00 2001 From: Pavel Tomin Date: Tue, 15 Oct 2024 14:06:25 -0500 Subject: [PATCH] fix: Bug fix for residual computation introduced in https://github.com/GEOS-DEV/GEOS/pull/3364 --- .../fluidFlow/ThermalCompositionalMultiphaseBaseKernels.hpp | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/src/coreComponents/physicsSolvers/fluidFlow/ThermalCompositionalMultiphaseBaseKernels.hpp b/src/coreComponents/physicsSolvers/fluidFlow/ThermalCompositionalMultiphaseBaseKernels.hpp index a27d879c215..d39b7f51ddc 100644 --- a/src/coreComponents/physicsSolvers/fluidFlow/ThermalCompositionalMultiphaseBaseKernels.hpp +++ b/src/coreComponents/physicsSolvers/fluidFlow/ThermalCompositionalMultiphaseBaseKernels.hpp @@ -911,7 +911,7 @@ class ResidualNormKernel : public solverBaseKernels::ResidualNormKernelBase< 3 > // step 2: volume residual real64 const valVol = LvArray::math::abs( m_localResidual[stack.localRow + m_numComponents] ) / volumeNormalizer; - if( valVol > stack.localValue[0] ) + if( valVol > stack.localValue[1] ) { stack.localValue[1] = valVol; } @@ -919,7 +919,7 @@ class ResidualNormKernel : public solverBaseKernels::ResidualNormKernelBase< 3 > // step 3: energy residual real64 const valEnergy = LvArray::math::abs( m_localResidual[stack.localRow + m_numComponents + 1] ) / energyNormalizer; - if( valEnergy > stack.localValue[1] ) + if( valEnergy > stack.localValue[2] ) { stack.localValue[2] = valEnergy; }