'Test predictions using erfa/plan94.
Accuracies are maximum deviations listed in erfa/plan94.c, for Jupiter and
Mercury, and that quoted in Meeus "Astronomical Algorithms" (1998) for the Moon.'

@pytest.mark.parametrize((u'body', u'sep_tol', u'dist_tol'), ((u'mercury', (7.0 * u.arcsec), (1000 * u.km)), (u'jupiter', (78.0 * u.arcsec), (76000 * u.km)), (u'moon', (20.0 * u.arcsec), (80 * u.km)), (u'sun', (5.0 * u.arcsec), (11.0 * u.km))))
def test_erfa_planet(self, body, sep_tol, dist_tol):

    astropy = get_body(body, self.t, ephemeris=u'builtin')
    horizons = self.horizons[body]
    astropy = _apparent_position_in_true_coordinates(astropy)
    assert (astropy.separation(horizons) < sep_tol)
    assert_quantity_allclose(astropy.distance, horizons.distance, atol=dist_tol)
