From 4cd8fcb495dc01a08dcf7943d05a7cb2ee7e9573 Mon Sep 17 00:00:00 2001 From: Chun Shen Date: Wed, 23 Sep 2026 07:51:22 -0400 Subject: [PATCH] fixed the bug based on Claude's suggestion --- src/framework/SurfaceFinder.cc | 6 +++--- 1 file changed, 3 insertions(+), 3 deletions(-) diff --git a/src/framework/SurfaceFinder.cc b/src/framework/SurfaceFinder.cc index 98a5a108..76e76650 100644 --- a/src/framework/SurfaceFinder.cc +++ b/src/framework/SurfaceFinder.cc @@ -505,9 +505,9 @@ SurfaceCellInfo SurfaceFinder::PrepareASurfaceCell( temp_cell.mu_Q = fluid_cell.mu_C; temp_cell.mu_S = fluid_cell.mu_S; - double u0 = - sqrt(1. + fluid_cell.vx * fluid_cell.vx + fluid_cell.vy * fluid_cell.vy + - fluid_cell.vz * fluid_cell.vz); + double v2 = (fluid_cell.vx * fluid_cell.vx + fluid_cell.vy * fluid_cell.vy + + fluid_cell.vz * fluid_cell.vz); + double u0 = 1. / sqrt(1. - v2); double uz = u0 * fluid_cell.vz; temp_cell.umu[0] = u0 * cosh(eta) - uz * sinh(eta); temp_cell.umu[1] = u0 * fluid_cell.vx;