get h from SBPoperator

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