Fix a bug in AMReX-Codes#4131 (AMReX-Codes#4138)

Missing sqrt.
baperry2 · Sep 6, 2024 · 87de52e · 87de52e
1 parent b454719
commit 87de52e
Showing 1 changed file with 3 additions and 3 deletions.
diff --git a/Src/EB/AMReX_EB2_3D_C.cpp b/Src/EB/AMReX_EB2_3D_C.cpp
@@ -101,9 +101,9 @@ void set_eb_data (const int i, const int j, const int k,
  bnorm(i,j,k,0) = nx;
  bnorm(i,j,k,1) = ny;
  bnorm(i,j,k,2) = nz;
- barea(i,j,k) = (nx*dapx + ny*dapy + nz*dapz) / (Math::powi<2>(nx*dx[1]*dx[2]) +
- Math::powi<2>(ny*dx[0]*dx[2]) +
- Math::powi<2>(nz*dx[0]*dx[1]));
+ barea(i,j,k) = (nx*dapx + ny*dapy + nz*dapz) / std::sqrt(Math::powi<2>(nx*dx[1]*dx[2]) +
+  Math::powi<2>(ny*dx[0]*dx[2]) +
+  Math::powi<2>(nz*dx[0]*dx[1]));
 
  Real aax = 0.5_rt*(axm+axp);
  Real aay = 0.5_rt*(aym+ayp);