def METHOD_NAME(axis, order, use_dask):
    mc_hdu = moment_cube()
    sc = SpectralCube.read(mc_hdu, use_dask=use_dask)
    cwise = sc.moment(axis=axis, order=order, how='cube')
    swise = sc.moment(axis=axis, order=order, how='slice')
    rwise = sc.moment(axis=axis, order=order, how='ray')
    assert_allclose(cwise, swise, rtol=rtol, atol=atol)
    assert_allclose(cwise, rwise, rtol=rtol, atol=atol)