From 147f3040c361e64759ec89c3100f9b7c7b42433a Mon Sep 17 00:00:00 2001 From: Magnus Ulimoen Date: Fri, 10 Apr 2020 01:07:57 +0200 Subject: [PATCH] divide instead of multiply h --- sbp/src/euler.rs | 8 ++++---- 1 file changed, 4 insertions(+), 4 deletions(-) diff --git a/sbp/src/euler.rs b/sbp/src/euler.rs index b7507bc..66ed772 100644 --- a/sbp/src/euler.rs +++ b/sbp/src/euler.rs @@ -630,7 +630,7 @@ fn SAT_characteristics( ) { // North boundary { - let hi = (k.ny() - 1) as Float * SBP::h()[0]; + let hi = (k.ny() - 1) as Float / SBP::h()[0]; let sign = -1.0; let tau = 1.0; let slice = s![y.ny() - 1, ..]; @@ -648,7 +648,7 @@ fn SAT_characteristics( } // South boundary { - let hi = (k.ny() - 1) as Float * SBP::h()[0]; + let hi = (k.ny() - 1) as Float / SBP::h()[0]; let sign = 1.0; let tau = -1.0; let slice = s![0, ..]; @@ -666,7 +666,7 @@ fn SAT_characteristics( } // West Boundary { - let hi = (k.nx() - 1) as Float * SBP::h()[0]; + let hi = (k.nx() - 1) as Float / SBP::h()[0]; let sign = 1.0; let tau = -1.0; let slice = s![.., 0]; @@ -684,7 +684,7 @@ fn SAT_characteristics( } // East Boundary { - let hi = (k.nx() - 1) as Float * SBP::h()[0]; + let hi = (k.nx() - 1) as Float / SBP::h()[0]; let sign = -1.0; let tau = 1.0; let slice = s![.., y.nx() - 1];