def METHOD_NAME(self, space='system'):
    a = np.random.rand(2,3,4,5)
    a = a.astype(np.float64)
    b = a.transpose(0,3,2,1).copy()
    c = bf.zeros(a.shape, dtype=a.dtype, space='system')
    c[...] = a
    c = c.copy(space=space)
    d = c.transpose(0,3,2,1).copy(space='system')
    # Use ctypes to directly access the memory
    b_data = ctypes.cast(b.ctypes.data, ctypes.POINTER(ctypes.c_double))
    b_data = np.array([b_data[i] for i in range(b.size)])
    d_data = ctypes.cast(d.ctypes.data, ctypes.POINTER(ctypes.c_double))
    d_data = np.array([d_data[i] for i in range(d.size)])
    np.testing.assert_equal(d_data, b_data)