absolute value distance to fix fpe error

zhichen3 · Nov 13, 2024 · 942be8a · 942be8a
1 parent 76380d9
commit 942be8a
Show file tree

Hide file tree

Showing 2 changed files with 2 additions and 2 deletions.
diff --git a/Source/driver/Castro_util.H b/Source/driver/Castro_util.H
@@ -174,7 +174,7 @@ Real distance(GeometryData const& geomdata, GpuArray<Real, 3>& loc)
         return std::sqrt(loc[0]*loc[0] + loc[1]*loc[1]);
     }
 
-    return loc[0];
+    return std::abs(loc[0]);
 }
 
 

diff --git a/Source/driver/Derive.cpp b/Source/driver/Derive.cpp
@@ -661,7 +661,7 @@ extern "C"
           der(i,j,k,0) = (-dat(i,j,k,1)*loc[1] + dat(i,j,k,2)*loc[0]) / (dat(i,j,k,0)*r);
 #endif
         } else {
-            Real r = distance(geomdata, loc);
+          Real r = distance(geomdata, loc);
 
           // we really mean just the velocity component that is
           // perpendicular to radial, and in general 3-d (e.g. a
-Original file line number
+Diff line change
@@ Expand Up @@
             return std::sqrt(loc[0]*loc[0] + loc[1]*loc[1]);
         }
-        return loc[0];
+        return std::abs(loc[0]);
     }
@@ Expand Down @@