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()))