def test_field2d_d_x2(): fld = fls.Field2D(2, 1, 2, 1, 10, 1, int(5)) assert np.allclose(fld.d_x2().toarray(), [[-2, 1, 0, 0], [1, -2, 1, 0], [0, 1, -2, 1], [0, 0, 1, -2]])