From 0ee91f785084fe6be1b1164cb2d92a857a26404b Mon Sep 17 00:00:00 2001 From: Magnus Ulimoen Date: Tue, 3 Sep 2019 19:44:54 +0200 Subject: [PATCH] get h from SBPoperator --- src/maxwell.rs | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/src/maxwell.rs b/src/maxwell.rs index 89ed152..3cb4313 100644 --- a/src/maxwell.rs +++ b/src/maxwell.rs @@ -104,7 +104,7 @@ impl System { let ny = y.0.shape()[0]; let nx = y.0.shape()[1]; - let h = 49.0 / 144.0 / (nx - 1) as f32; // TODO: Get from schema + let h = SBP::h()[0] / (nx - 1) as f32; let hinv = 1.0 / h; // East boundary @@ -130,7 +130,7 @@ impl System { k[i].2[(j, 0)] += tau * hinv * (-0.5 * (v.1 - g.1) - 0.5 * (v.2 - g.2)); } - let h = 49.0 / 144.0 / (ny - 1) as f32; // TODO: same as above + let h = SBP::h()[0] / (ny - 1) as f32; let hinv = 1.0 / h; // North boundary