diff --git a/Team32/magstadyn_a.pro b/Team32/magstadyn_a.pro
index d968275e4a51fc06c7b70689587e8b39f62ef49a..50313358e26c35310eaea54ff8befb48740d93ee 100644
--- a/Team32/magstadyn_a.pro
+++ b/Team32/magstadyn_a.pro
@@ -1067,6 +1067,11 @@ Resolution {
                   Evaluate[$relaxcount=$NbrTestedFac];
                 EndIf
 
+                // To compute Residual based on PostOperation h_only
+                If (Flag_UpdateSeparated && Flag_ItLoopNRes==NLITLOOPN_POSTOPH)
+                  Generate[UpH]; Solve[UpH];
+                EndIf 
+              
                 Evaluate[$syscount = $syscount + 1 ];
                 Evaluate[$relaxcounttot=$relaxcounttot+$relaxcount];
                 Evaluate[$relax= $RelaxFac];
diff --git a/Team32/magstadyn_avs_3d.pro b/Team32/magstadyn_avs_3d.pro
index 1f0f2e1f0b2eb03eb17ccfc0bcfcec1acdc19d54..63d6b4b2a6760ec384c73a6efb15ac88696f4de2 100644
--- a/Team32/magstadyn_avs_3d.pro
+++ b/Team32/magstadyn_avs_3d.pro
@@ -1049,6 +1049,11 @@ Resolution {
                   Evaluate[$relaxcount=$NbrTestedFac];
                 EndIf
 
+                // To compute Residual based on PostOperation h_only
+                If (Flag_UpdateSeparated && Flag_ItLoopNRes==NLITLOOPN_POSTOPH)
+                  Generate[UpH]; Solve[UpH];
+                EndIf 
+                
                 Evaluate[$syscount = $syscount + 1 ];
                 Evaluate[$relaxcounttot=$relaxcounttot+$relaxcount];
                 Evaluate[$relax= $RelaxFac];