def test_fresnel_integral(self):
v = np.array([-1, 0, 1, 2.4])
o = np.array([
-20 * np.log(1.12), -20 * np.log(0.5),
-20 * np.log(0.4 - (0.1184 - 0.28**2)**0.5),
-20 * np.log(0.225 / 2.4)
])
npt.assert_almost_equal(geo.fresnel_integral(v), o)
def test_get_fresnel_single(self):
five = shapely.geometry.LineString([[527990,183005,0],[528020,183005,15]])
expected=geo.fresnel_integral([5 *( 2 / (0.1903 * 10))**0.5])
self.assertAlmostEqual(geo.get_fresnel(five,self.buildings,[10]),expected[0])
def test_fresnel(self):
a = geo.fresnel(self.map_box, self.rays_box)
five = geo.fresnel_integral([5 *( 2 / (0.1903 * 10))**0.5])[0]
b = pd.Series(data=[five,0],index=[1,2])
pdt.assert_series_equal(a, b)