wrote a test to make sure it gives NaN outside of interpolation range

eas342 · Feb 6, 2021 · 31289f0 · 31289f0
1 parent 430e7f4
commit 31289f0
Showing 1 changed file with 5 additions and 1 deletion.
diff --git a/tests/test_mie.py b/tests/test_mie.py
@@ -41,7 +41,11 @@ def test_negative_r_for_dist(self):
         with self.assertRaisesRegex(ValueError,'Negative or zero radius not allowed') as cm:
             calc_mie.get_mie_coeff_distribution(1.0,-0.3)
 
-
+    def test_outside_interpolation_range(self):
+        wave = np.array([1e-5,2.5])
+        qext, qsca, qback, g = calc_mie.get_mie_coeff(wave,r=1.0,material='Fe2SiO4')
+        self.assertEqual(np.isfinite(qext)[0],False)
+
     def test_large_r(self):
         res = calc_mie.get_mie_coeff_distribution(0.8,r=1000.,material='Mg2SiO4')