def output_catalog_tests(dates): polaris1991 = novas.make_cat_entry('11767', 'HIP', 0, 0.0, 89.26413805, 44.22, -11.74, 7.56, 0.0) polaris1991.ra = 37.94614689 # HIP uses degrees, not hours polaris = novas.transform_hip(polaris1991) for i, jd in enumerate(dates): ra, dec = call(novas.astro_star, jd, polaris) output( locals(), r""" def test_hipparcos_conversion{i}(earth): line = b'H| 11767| |02 31 47.08|+89 15 50.9| 1.97|1|H|037.94614689|+89.26413805| | 7.56| 44.22| -11.74| 0.39| 0.45| 0.48| 0.47| 0.55|-0.16| 0.05| 0.27|-0.01| 0.08| 0.05| 0.04|-0.12|-0.09|-0.36| 1| 1.22| 11767| 2.756|0.003| 2.067|0.003| | 0.636|0.003|T|0.70|0.00|L| | 2.1077|0.0021|0.014|102| | 2.09| 2.13| 3.97|P|1|A|02319+8915|I| 1| 1| | | | | | | | | |S| |P| 8890|B+88 8 | | |0.68|F7:Ib-IIv SB|G\n' df = hipparcos.load_dataframe(BytesIO(line)) star = starlib.Star.from_dataframe(df.iloc[0]) ra, dec, distance = earth.at(load.timescale().tt_jd({jd})).observe(star).radec() compare(ra.hours, {ra!r}, 0.00001 * ra_arcsecond) compare(dec.degrees, {dec!r}, 0.00001 * arcsecond) """)
def output_catalog_tests(dates): polaris1991 = novas.make_cat_entry( '11767', 'HIP', 0, 0.0, 89.26413805, 44.22, -11.74, 7.56, 0.0) polaris1991.ra = 37.94614689 # HIP uses degrees, not hours polaris = novas.transform_hip(polaris1991) for i, jd in enumerate(dates): ra, dec = call(novas.astro_star, jd, polaris) output(locals(), r""" def test_hipparcos_conversion{i}(): line = 'H| 11767| |02 31 47.08|+89 15 50.9| 1.97|1|H|037.94614689|+89.26413805| | 7.56| 44.22| -11.74| 0.39| 0.45| 0.48| 0.47| 0.55|-0.16| 0.05| 0.27|-0.01| 0.08| 0.05| 0.04|-0.12|-0.09|-0.36| 1| 1.22| 11767| 2.756|0.003| 2.067|0.003| | 0.636|0.003|T|0.70|0.00|L| | 2.1077|0.0021|0.014|102| | 2.09| 2.13| 3.97|P|1|A|02319+8915|I| 1| 1| | | | | | | | | |S| |P| 8890|B+88 8 | | |0.68|F7:Ib-IIv SB|G\n' star = hipparcos.parse(line) compare(star.ra.hours, {polaris.ra!r}, 0.001 * ra_arcsecond) compare(star.dec.degrees, {polaris.dec!r}, 0.001 * arcsecond) ra, dec, distance = de405.earth(tt={jd}).observe(star).radec() compare(ra.hours, {ra!r}, 0.001 * ra_arcsecond) compare(dec.degrees, {dec!r}, 0.001 * arcsecond) """)