def test_avg_two_maps_unbiased(self):
maps2 = Maps(self.array2, Ni=self.Ni, Nj=self.Nj, Nk=self.Nk)
sigma = 2.
avg, _ = maps2.iterative_smooth_avg_var(sigma=sigma, bias=False)
maps2.smooth(sigma=sigma, inplace=True)
self.assertTrue(np.array_equal(maps2.avg().to_array(), avg.to_array()))
def test_var_two_maps_unbiased(self):
maps2 = Maps(self.array2, Ni=self.Ni, Nj=self.Nj, Nk=self.Nk)
sigma = 2.
_, var = maps2.iterative_smooth_avg_var(sigma=sigma, bias=False)
maps2.smooth(sigma=sigma, inplace=True)
self.assertTrue(
np.allclose(maps2.var(bias=False).to_array(), var.to_array()))