def test_axis_1_derivative_1(self):
B = self.F.simple_operators[(1,)]
BG = BOG.for_vector(self.F.grid.mesh[1],
self.F.grid.shape[0], 1, 1, 0).immigrate()
npt.assert_array_equal(B.D.data, BG.D.data)
npt.assert_equal(B, BG)
def test_axis_1_derivative_2(self):
B = self.F.simple_operators[(1,1)]
BG = BOG.for_vector(self.F.grid.mesh[1],
self.F.grid.shape[0], 2, 1, 0).immigrate()
fp(B.D.data)
print
fp(BG.D.data)
print
fp(B.D.data - BG.D.data)
print B.D.offsets
fp(B.D)
print BG.D.offsets
fp(BG.D)
npt.assert_array_equal(B.D.data, BG.D.data)
npt.assert_equal(B, BG)
