From 7cdfd4006ab122396283e4819991c2a989a7c635 Mon Sep 17 00:00:00 2001 From: JamesEMcclure Date: Sat, 27 Jun 2020 07:30:46 -0400 Subject: [PATCH] Eikonal solver in distance --- analysis/distance.cpp | 12 ++++++------ 1 file changed, 6 insertions(+), 6 deletions(-) diff --git a/analysis/distance.cpp b/analysis/distance.cpp index 78c344d2..3b7641d6 100644 --- a/analysis/distance.cpp +++ b/analysis/distance.cpp @@ -272,24 +272,24 @@ double Eikonal(DoubleArray &Distance, char *ID, Domain &Dm, int timesteps){ // Compute upwind derivatives for Godunov Hamiltonian if (sign < 0.0){ if (Dxp + Dxm > 0.f) Dx = Dxp*Dxp; - elseDx = Dxm*Dxm; + else Dx = Dxm*Dxm; if (Dyp + Dym > 0.f) Dy = Dyp*Dyp; - elseDy = Dym*Dym; + else Dy = Dym*Dym; if (Dzp + Dzm > 0.f) Dz = Dzp*Dzp; - elseDz = Dzm*Dzm; + else Dz = Dzm*Dzm; } else{ if (Dxp + Dxm < 0.f) Dx = Dxp*Dxp; - elseDx = Dxm*Dxm; + else Dx = Dxm*Dxm; if (Dyp + Dym < 0.f) Dy = Dyp*Dyp; - elseDy = Dym*Dym; + else Dy = Dym*Dym; if (Dzp + Dzm < 0.f) Dz = Dzp*Dzp; - elseDz = Dzm*Dzm; + else Dz = Dzm*Dzm; } //Dx = max(Dxp*Dxp,Dxm*Dxm);