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