//  [python]
def test_orbit_method_inputvo_quantity():
    from astropy import units

    from galpy.orbit import Orbit
    from galpy.potential import MWPotential2014

    o = Orbit([1.1, 0.1, 1.1, 0.2, 0.3, 0.3])
    vo = 222.0
    assert (
        numpy.fabs(
            o.E(pot=MWPotential2014, vo=vo * units.km / units.s)
            - o.E(pot=MWPotential2014, vo=vo)
        )
        < 10.0**-8.0
    ), "Orbit method E does not return the correct value when input vo is Quantity"
    assert (
        numpy.fabs(
            o.ER(pot=MWPotential2014, vo=vo * units.km / units.s)
            - o.ER(pot=MWPotential2014, vo=vo)
        )
        < 10.0**-8.0
    ), "Orbit method ER does not return the correct value as Quantity"
    assert (
        numpy.fabs(
            o.Ez(pot=MWPotential2014, vo=vo * units.km / units.s)
            - o.Ez(pot=MWPotential2014, vo=vo)
        )
        < 10.0**-8.0
    ), "Orbit method Ez does not return the correct value when input vo is Quantity"
    assert (
        numpy.fabs(
            o.Jacobi(pot=MWPotential2014, vo=vo * units.km / units.s)
            - o.Jacobi(pot=MWPotential2014, vo=vo)
        )
        < 10.0**-8.0
    ), "Orbit method Jacobi does not return the correct value when input vo is Quantity"
    assert numpy.all(
        numpy.fabs(
            o.L(pot=MWPotential2014, vo=vo * units.km / units.s)
            - o.L(pot=MWPotential2014, vo=vo)
        )
        < 10.0**-8.0
    ), "Orbit method L does not return the correct value when input vo is Quantity"
    assert (
        numpy.fabs(
            o.rap(pot=MWPotential2014, analytic=True, vo=vo * units.km / units.s)
            - o.rap(pot=MWPotential2014, analytic=True, vo=vo)
        )
        < 10.0**-8.0
    ), "Orbit method rap does not return the correct value when input vo is Quantity"
    assert (
        numpy.fabs(
            o.rperi(pot=MWPotential2014, analytic=True, vo=vo * units.km / units.s)
            - o.rperi(pot=MWPotential2014, analytic=True, vo=vo)
        )
        < 10.0**-8.0
    ), "Orbit method rperi does not return the correct value when input vo is Quantity"
    assert (
        numpy.fabs(
            o.zmax(pot=MWPotential2014, analytic=True, vo=vo * units.km / units.s)
            - o.zmax(pot=MWPotential2014, analytic=True, vo=vo)
        )
        < 10.0**-8.0
    ), "Orbit method zmax does not return the correct value when input vo is Quantity"
    assert (
        numpy.fabs(
            o.jr(
                pot=MWPotential2014,
                type="staeckel",
                delta=0.5,
                vo=vo * units.km / units.s,
            )
            - o.jr(pot=MWPotential2014, type="staeckel", delta=0.5, vo=vo)
        )
        < 10.0**-8.0
    ), "Orbit method jr does not return the correct value when input vo is Quantity"
    assert (
        numpy.fabs(
            o.jp(
                pot=MWPotential2014,
                type="staeckel",
                delta=0.5,
                vo=vo * units.km / units.s,
            )
            - o.jp(pot=MWPotential2014, type="staeckel", delta=0.5, vo=vo)
        )
        < 10.0**-8.0
    ), "Orbit method jp does not return the correct value when input vo is Quantity"
    assert (
        numpy.fabs(
            o.jz(
                pot=MWPotential2014,
                type="staeckel",
                delta=0.5,
                vo=vo * units.km / units.s,
            )
            - o.jz(pot=MWPotential2014, type="staeckel", delta=0.5, vo=vo)
        )
        < 10.0**-8.0
    ), "Orbit method jz does not return the correct value when input vo is Quantity"
    assert (
        numpy.fabs(
            o.wr(
                pot=MWPotential2014,
                type="staeckel",
                delta=0.5,
                vo=vo * units.km / units.s,
            )
            - o.wr(pot=MWPotential2014, type="staeckel", delta=0.5, vo=vo)
        )
        < 10.0**-8.0
    ), "Orbit method wr does not return the correct value when input vo is Quantity"
    assert (
        numpy.fabs(
            o.wp(
                pot=MWPotential2014,
                type="staeckel",
                delta=0.5,
                vo=vo * units.km / units.s,
            )
            - o.wp(pot=MWPotential2014, type="staeckel", delta=0.5, vo=vo)
        )
        < 10.0**-8.0
    ), "Orbit method wp does not return the correct value when input vo is Quantity"
    assert (
        numpy.fabs(
            o.wz(
                pot=MWPotential2014,
                type="staeckel",
                delta=0.5,
                vo=vo * units.km / units.s,
            )
            - o.wz(pot=MWPotential2014, type="staeckel", delta=0.5, vo=vo)
        )
        < 10.0**-8.0
    ), "Orbit method wz does not return the correct value when input vo is Quantity"
    assert (
        numpy.fabs(
            o.Tr(
                pot=MWPotential2014,
                type="staeckel",
                delta=0.5,
                vo=vo * units.km / units.s,
            )
            - o.Tr(pot=MWPotential2014, type="staeckel", delta=0.5, vo=vo)
        )
        < 10.0**-8.0
    ), "Orbit method Tr does not return the correct value when input vo is Quantity"
    assert (
        numpy.fabs(
            o.Tp(
                pot=MWPotential2014,
                type="staeckel",
                delta=0.5,
                vo=vo * units.km / units.s,
            )
            - o.Tp(pot=MWPotential2014, type="staeckel", delta=0.5, vo=vo)
        )
        < 10.0**-8.0
    ), "Orbit method Tp does not return the correct value when input vo is Quantity"
    assert (
        numpy.fabs(
            o.Tz(
                pot=MWPotential2014,
                type="staeckel",
                delta=0.5,
                vo=vo * units.km / units.s,
            )
            - o.Tz(pot=MWPotential2014, type="staeckel", delta=0.5, vo=vo)
        )
        < 10.0**-8.0
    ), "Orbit method Tz does not return the correct value when input vo is Quantity"
    assert (
        numpy.fabs(
            o.Or(
                pot=MWPotential2014,
                type="staeckel",
                delta=0.5,
                vo=vo * units.km / units.s,
            )
            - o.Or(pot=MWPotential2014, type="staeckel", delta=0.5, vo=vo)
        )
        < 10.0**-8.0
    ), "Orbit method Or does not return the correct value when input vo is Quantity"
    assert (
        numpy.fabs(
            o.Op(
                pot=MWPotential2014,
                type="staeckel",
                delta=0.5,
                vo=vo * units.km / units.s,
            )
            - o.Op(pot=MWPotential2014, type="staeckel", delta=0.5, vo=vo)
        )
        < 10.0**-8.0
    ), "Opbit method Or does not return the correct value when input vo is Quantity"
    assert (
        numpy.fabs(
            o.Oz(
                pot=MWPotential2014,
                type="staeckel",
                delta=0.5,
                vo=vo * units.km / units.s,
            )
            - o.Oz(pot=MWPotential2014, type="staeckel", delta=0.5, vo=vo)
        )
        < 10.0**-8.0
    ), "Ozbit method Or does not return the correct value when input vo is Quantity"
    assert (
        numpy.fabs(o.time(vo=vo * units.km / units.s) - o.time(vo=vo)) < 10.0**-8.0
    ), "Orbit method time does not return the correct value when input vo is Quantity"
    assert (
        numpy.fabs(o.R(vo=vo * units.km / units.s) - o.R(vo=vo)) < 10.0**-8.0
    ), "Orbit method R does not return the correct value when input vo is Quantity"
    assert (
        numpy.fabs(o.vR(vo=vo * units.km / units.s) - o.vR(vo=vo)) < 10.0**-8.0
    ), "Orbit method vR does not return the correct value when input vo is Quantity"
    assert (
        numpy.fabs(o.vT(vo=vo * units.km / units.s) - o.vT(vo=vo)) < 10.0**-8.0
    ), "Orbit method vT does not return the correct value when input vo is Quantity"
    assert (
        numpy.fabs(o.z(vo=vo * units.km / units.s) - o.z(vo=vo)) < 10.0**-8.0
    ), "Orbit method z does not return the correct value when input vo is Quantity"
    assert (
        numpy.fabs(o.vz(vo=vo * units.km / units.s) - o.vz(vo=vo)) < 10.0**-8.0
    ), "Orbit method vz does not return the correct value when input vo is Quantity"
    assert (
        numpy.fabs(o.phi(vo=vo * units.km / units.s) - o.phi(vo=vo)) < 10.0**-8.0
    ), "Orbit method phi does not return the correct value when input vo is Quantity"
    assert (
        numpy.fabs(o.vphi(vo=vo * units.km / units.s) - o.vphi(vo=vo)) < 10.0**-8.0
    ), "Orbit method vphi does not return the correct value when input vo is Quantity"
    assert (
        numpy.fabs(o.x(vo=vo * units.km / units.s) - o.x(vo=vo)) < 10.0**-8.0
    ), "Orbit method x does not return the correct value when input vo is Quantity"
    assert (
        numpy.fabs(o.y(vo=vo * units.km / units.s) - o.y(vo=vo)) < 10.0**-8.0
    ), "Orbit method y does not return the correct value when input vo is Quantity"
    assert (
        numpy.fabs(o.vx(vo=vo * units.km / units.s) - o.vx(vo=vo)) < 10.0**-8.0
    ), "Orbit method vx does not return the correct value when input vo is Quantity"
    assert (
        numpy.fabs(o.vy(vo=vo * units.km / units.s) - o.vy(vo=vo)) < 10.0**-8.0
    ), "Orbit method vy does not return the correct value when input vo is Quantity"
    assert (
        numpy.fabs(o.theta(vo=vo * units.km / units.s) - o.theta(vo=vo)) < 10.0**-8.0
    ), "Orbit method theta does not return the correct value when input vo is Quantity"
    assert (
        numpy.fabs(o.vtheta(vo=vo * units.km / units.s) - o.vtheta(vo=vo)) < 10.0**-8.0
    ), "Orbit method vtheta does not return the correct value when input vo is Quantity"
    assert (
        numpy.fabs(o.vr(vo=vo * units.km / units.s) - o.vr(vo=vo)) < 10.0**-8.0
    ), "Orbit method vr does not return the correct value when input vo is Quantity"
    assert (
        numpy.fabs(o.ra(vo=vo * units.km / units.s) - o.ra(vo=vo)) < 10.0**-8.0
    ), "Orbit method ra does not return the correct value when input vo is Quantity"
    assert (
        numpy.fabs(o.dec(vo=vo * units.km / units.s) - o.dec(vo=vo)) < 10.0**-8.0
    ), "Orbit method dec does not return the correct value when input vo is Quantity"
    assert (
        numpy.fabs(o.ll(vo=vo * units.km / units.s) - o.ll(vo=vo)) < 10.0**-8.0
    ), "Orbit method ll does not return the correct value when input vo is Quantity"
    assert (
        numpy.fabs(o.bb(vo=vo * units.km / units.s) - o.bb(vo=vo)) < 10.0**-8.0
    ), "Orbit method bb does not return the correct value when input vo is Quantity"
    assert (
        numpy.fabs(o.dist(vo=vo * units.km / units.s) - o.dist(vo=vo)) < 10.0**-8.0
    ), "Orbit method dist does not return the correct value when input vo is Quantity"
    assert (
        numpy.fabs(o.pmra(vo=vo * units.km / units.s) - o.pmra(vo=vo)) < 10.0**-8.0
    ), "Orbit method pmra does not return the correct value when input vo is Quantity"
    assert (
        numpy.fabs(o.pmdec(vo=vo * units.km / units.s) - o.pmdec(vo=vo)) < 10.0**-8.0
    ), "Orbit method pmdec does not return the correct value when input vo is Quantity"
    assert (
        numpy.fabs(o.pmll(vo=vo * units.km / units.s) - o.pmll(vo=vo)) < 10.0**-8.0
    ), "Orbit method pmll does not return the correct value when input vo is Quantity"
    assert (
        numpy.fabs(o.pmbb(vo=vo * units.km / units.s) - o.pmbb(vo=vo)) < 10.0**-8.0
    ), "Orbit method pmbb does not return the correct value when input vo is Quantity"
    assert (
        numpy.fabs(o.vlos(vo=vo * units.km / units.s) - o.vlos(vo=vo)) < 10.0**-8.0
    ), "Orbit method vlos does not return the correct value when input vo is Quantity"
    assert (
        numpy.fabs(o.vra(vo=vo * units.km / units.s) - o.vra(vo=vo)) < 10.0**-8.0
    ), "Orbit method vra does not return the correct value when input vo is Quantity"
    assert (
        numpy.fabs(o.vdec(vo=vo * units.km / units.s) - o.vdec(vo=vo)) < 10.0**-8.0
    ), "Orbit method vdec does not return the correct value when input vo is Quantity"
    assert (
        numpy.fabs(o.vll(vo=vo * units.km / units.s) - o.vll(vo=vo)) < 10.0**-8.0
    ), "Orbit method vll does not return the correct value when input vo is Quantity"
    assert (
        numpy.fabs(o.vbb(vo=vo * units.km / units.s) - o.vbb(vo=vo)) < 10.0**-8.0
    ), "Orbit method vbb does not return the correct value when input vo is Quantity"
    assert (
        numpy.fabs(o.helioX(vo=vo * units.km / units.s) - o.helioX(vo=vo)) < 10.0**-8.0
    ), "Orbit method helioX does not return the correct value when input vo is Quantity"
    assert (
        numpy.fabs(o.helioY(vo=vo * units.km / units.s) - o.helioY(vo=vo)) < 10.0**-8.0
    ), "Orbit method helioY does not return the correct value when input vo is Quantity"
    assert (
        numpy.fabs(o.helioZ(vo=vo * units.km / units.s) - o.helioZ(vo=vo)) < 10.0**-8.0
    ), "Orbit method helioZ does not return the correct value when input vo is Quantity"
    assert (
        numpy.fabs(o.U(vo=vo * units.km / units.s) - o.U(vo=vo)) < 10.0**-8.0
    ), "Orbit method U does not return the correct value when input vo is Quantity"
    assert (
        numpy.fabs(o.V(vo=vo * units.km / units.s) - o.V(vo=vo)) < 10.0**-8.0
    ), "Orbit method V does not return the correct value when input vo is Quantity"
    assert (
        numpy.fabs(o.W(vo=vo * units.km / units.s) - o.W(vo=vo)) < 10.0**-8.0
    ), "Orbit method W does not return the correct value when input vo is Quantity"
    return None