diff --git a/SU2_CFD/src/numerics/NEMO/NEMO_diffusion.cpp b/SU2_CFD/src/numerics/NEMO/NEMO_diffusion.cpp
index 2553fa2915d..5a23ff0b6ce 100644
--- a/SU2_CFD/src/numerics/NEMO/NEMO_diffusion.cpp
+++ b/SU2_CFD/src/numerics/NEMO/NEMO_diffusion.cpp
@@ -169,7 +169,7 @@ CNumerics::ResidualType<> CAvgGrad_NEMO::ComputeResidual(const CConfig *config)
 
   su2double dist_ij_2[MAXNDIM] = {0.0};
   GeometryToolbox::Distance(nDim, Coord_j, Coord_i, dist_ij_2);
-  dist_ij = GeometryToolbox::SquaredNorm(nDim, dist_ij_2);
+  dist_ij = GeometryToolbox::Norm(nDim, dist_ij_2);
 
     for (auto iVar = 0ul; iVar < nVar; iVar++) {
       for (auto jVar = 0ul; jVar < nVar; jVar++) {
diff --git a/TestCases/parallel_regression.py b/TestCases/parallel_regression.py
index ff9b8c66706..69a3f2fb5ae 100755
--- a/TestCases/parallel_regression.py
+++ b/TestCases/parallel_regression.py
@@ -192,8 +192,8 @@ def main():
     visc_cone.cfg_dir = "nonequilibrium/visc_wedge"
     visc_cone.cfg_file = "axi_visccone.cfg"
     visc_cone.test_iter = 10
-    visc_cone.test_vals = [-5.215234, -5.739371, -20.559852, -20.509281, -20.408911, 1.262701, -3.205457, -0.015696, 0.093205, 32637.000000]
-    visc_cone.test_vals_aarch64 = [-5.222270, -5.746525, -20.560286, -20.510152, -20.409101, 1.255758, -3.208382, -0.016014, 0.093462, 32619.000000]
+    visc_cone.test_vals = [-5.215238, -5.739375, -20.581529, -20.580261, -20.508727, 1.262697, -3.205457, -0.015697, 0.093206, 32637.000000]
+    visc_cone.test_vals_aarch64 = [-5.215222, -5.739359, -20.581535, -20.580239, -20.508657, 1.262709, -3.205449, -0.015696, 0.093205, 32637.000000]
     test_list.append(visc_cone)
 
     # Viscous single wedge with Mutation++
@@ -209,7 +209,7 @@ def main():
     super_cat.cfg_dir = "nonequilibrium/visc_wedge"
     super_cat.cfg_file = "super_cat.cfg"
     super_cat.test_iter = 10
-    super_cat.test_vals = [-5.232595, -5.757889, -20.641415, -20.640623, -20.541670, 1.246866, -3.205258, -0.028372, 0.250647, 32440.000000]
+    super_cat.test_vals = [-5.232595, -5.757889, -20.726945, -20.749655, -20.567151, 1.246866, -3.205258, -0.028372, 0.250647, 32440.000000]
     test_list.append(super_cat)
 
     # Viscous single wedge - partially catalytic walls
@@ -225,7 +225,7 @@ def main():
     ion_gy.cfg_dir = "nonequilibrium/visc_cylinder"
     ion_gy.cfg_file = "cyl_ion_gy.cfg"
     ion_gy.test_iter = 10
-    ion_gy.test_vals = [-11.629873, -4.165562, -4.702662, -4.950351, -5.146155, -4.993878, -6.893332, 5.990109, 5.990004, -0.014849, 0.000000, 90090.000000]
+    ion_gy.test_vals = [-11.629871, -4.156741, -4.692113, -4.949447, -5.146196, -4.991664, -6.893329, 5.990109, 5.990004, -0.014849, 0.000000, 90090.000000]
     test_list.append(ion_gy)
 
     ##########################
diff --git a/TestCases/serial_regression.py b/TestCases/serial_regression.py
index aa22e30d8ae..82df6232a97 100755
--- a/TestCases/serial_regression.py
+++ b/TestCases/serial_regression.py
@@ -75,8 +75,8 @@ def main():
     visc_cone.cfg_dir = "nonequilibrium/visc_wedge"
     visc_cone.cfg_file = "axi_visccone.cfg"
     visc_cone.test_iter = 10
-    visc_cone.test_vals = [-5.215230, -5.739367, -20.560781, -20.516922, -20.406516, 1.262782, -3.205476, -0.015696, 0.093206, 32641]
-    visc_cone.test_vals_aarch64 = [-5.215250, -5.739384, -20.560917, -20.517096, -20.406630, 1.262772, -3.205492, -0.015695, 0.093205, 32641.000000]
+    visc_cone.test_vals = [-5.215233, -5.739371, -20.551239, -20.590906, -20.437427, 1.262779, -3.205476, -0.015697, 0.093207, 32641]
+    visc_cone.test_vals_aarch64 = [-5.215220, -5.739356, -20.551244, -20.590883, -20.437378, 1.262790, -3.205470, -0.015696, 0.093206, 32641.000000]
     test_list.append(visc_cone)
 
     #########################