def test_dd_perp_delta(self): den = self.zz den[self.nx/2, self.ny/2] = 1.0 den_diff = dens_diff_perp(den=den, idxs=range(self.nx/8), outer_len=1.0) # print # print den_diff ok_(not np.allclose(den_diff[:,1], 0.0))
def test_dd_perp_const(self): den = self.zz den.fill(0.1) den_diff = dens_diff_perp(den=den, idxs=range(10), outer_len=1.0) ok_(np.allclose(den_diff[:,1], 0.0))
def test_dd_perp_zeros(self): den_diff = dens_diff_perp(den=self.zz, idxs=range(10), outer_len=1.0) ok_(np.allclose(den_diff[:,1], 0.0))