@@ -259,13 +259,13 @@ let
259259 @test isapprox (astrom. bpn[1 ][3 ], - 0.1312227200895260194e-2 , atol= 1e-12 )
260260 @test isapprox (astrom. bpn[2 ][3 ], 0.2928082217872315680e-4 , atol= 1e-12 )
261261 @test isapprox (astrom. bpn[3 ][3 ], 0.9999991386008323373 , atol= 1e-12 )
262- @test isapprox (astrom. along, - 0.5278008060301974337 , atol= 1e-12 )
263- @test isapprox (astrom. xpl, 0.1133427418174939329e -5 , atol= 1e-17 )
264- @test isapprox (astrom. ypl, 0.1453347595745898629e -5 , atol= 1e-17 )
262+ @test isapprox (astrom. along, - 0.5278008060295995734 , atol= 1e-12 )
263+ @test isapprox (astrom. xpl, 0.1133427418130752958e -5 , atol= 1e-17 )
264+ @test isapprox (astrom. ypl, 0.1453347595780646207e -5 , atol= 1e-17 )
265265 @test isapprox (astrom. sphi, - 0.9440115679003211329 , atol= 1e-12 )
266266 @test isapprox (astrom. cphi, 0.3299123514971474711 , atol= 1e-12 )
267267 @test astrom. diurab == 0.0
268- @test isapprox (astrom. eral, 2.617608903969802566 , atol= 1e-12 )
268+ @test isapprox (astrom. eral, 2.617608903970400427 , atol= 1e-12 )
269269 @test isapprox (astrom. refa, 0.2014187790000000000e-3 , atol= 1e-15 )
270270 @test isapprox (astrom. refb, - 0.2361408310000000000e-6 , atol= 1e-18 )
271271end
@@ -308,13 +308,13 @@ let
308308 @test isapprox (astrom. bpn[1 ][3 ], - 0.001312227201745553566 , atol= 1e-12 )
309309 @test isapprox (astrom. bpn[2 ][3 ], 0.2928082218847679162e-4 , atol= 1e-12 )
310310 @test isapprox (astrom. bpn[3 ][3 ], 0.9999991386008312212 , atol= 1e-12 )
311- @test isapprox (astrom. along, - 0.5278008060301974337 , atol= 1e-12 )
312- @test isapprox (astrom. xpl, 0.1133427418174939329e -5 , atol= 1e-17 )
313- @test isapprox (astrom. ypl, 0.1453347595745898629e -5 , atol= 1e-17 )
311+ @test isapprox (astrom. along, - 0.5278008060295995733 , atol= 1e-12 )
312+ @test isapprox (astrom. xpl, 0.1133427418130752958e -5 , atol= 1e-17 )
313+ @test isapprox (astrom. ypl, 0.1453347595780646207e -5 , atol= 1e-17 )
314314 @test isapprox (astrom. sphi, - 0.9440115679003211329 , atol= 1e-12 )
315315 @test isapprox (astrom. cphi, 0.3299123514971474711 , atol= 1e-12 )
316316 @test astrom. diurab == 0.0
317- @test isapprox (astrom. eral, 2.617608909189066140 , atol= 1e-12 )
317+ @test isapprox (astrom. eral, 2.617608909189664000 , atol= 1e-12 )
318318 @test isapprox (astrom. refa, 0.2014187785940396921e-3 , atol= 1e-15 )
319319 @test isapprox (astrom. refb, - 0.2361408314943696227e-6 , atol= 1e-18 )
320320 @test isapprox (eo, - 0.003020548354802412839 , atol= 1e-14 )
@@ -430,13 +430,13 @@ let
430430
431431 astrom = iauApio (sp, theta, elong, phi, hm, xp, yp, refa, refb)
432432
433- @test isapprox (astrom. along, - 0.5278008060301974337 , atol= 1e-12 )
434- @test isapprox (astrom. xpl, 0.1133427418174939329e -5 , atol= 1e-17 )
435- @test isapprox (astrom. ypl, 0.1453347595745898629e -5 , atol= 1e-17 )
433+ @test isapprox (astrom. along, - 0.5278008060295995734 , atol= 1e-12 )
434+ @test isapprox (astrom. xpl, 0.1133427418130752958e -5 , atol= 1e-17 )
435+ @test isapprox (astrom. ypl, 0.1453347595780646207e -5 , atol= 1e-17 )
436436 @test isapprox (astrom. sphi, - 0.9440115679003211329 , atol= 1e-12 )
437437 @test isapprox (astrom. cphi, 0.3299123514971474711 , atol= 1e-12 )
438438 @test isapprox (astrom. diurab, 0.5135843661699913529e-6 , atol= 1e-12 )
439- @test isapprox (astrom. eral, 2.617608903969802566 , atol= 1e-12 )
439+ @test isapprox (astrom. eral, 2.617608903970400427 , atol= 1e-12 )
440440 @test isapprox (astrom. refa, 0.2014187790000000000e-3 , atol= 1e-15 )
441441 @test isapprox (astrom. refb, - 0.2361408310000000000e-6 , atol= 1e-18 )
442442end
@@ -458,13 +458,13 @@ let
458458 j, astrom = iauApio13 (utc1, utc2, dut1, elong, phi, hm, xp, yp,
459459 phpa, tc, rh, wl)
460460
461- @test isapprox (astrom. along, - 0.5278008060301974337 , atol= 1e-12 )
462- @test isapprox (astrom. xpl, 0.1133427418174939329e -5 , atol= 1e-17 )
463- @test isapprox (astrom. ypl, 0.1453347595745898629e -5 , atol= 1e-17 )
461+ @test isapprox (astrom. along, - 0.5278008060295995733 , atol= 1e-12 )
462+ @test isapprox (astrom. xpl, 0.1133427418130752958e -5 , atol= 1e-17 )
463+ @test isapprox (astrom. ypl, 0.1453347595780646207e -5 , atol= 1e-17 )
464464 @test isapprox (astrom. sphi, - 0.9440115679003211329 , atol= 1e-12 )
465465 @test isapprox (astrom. cphi, 0.3299123514971474711 , atol= 1e-12 )
466466 @test isapprox (astrom. diurab, 0.5135843661699913529e-6 , atol= 1e-12 )
467- @test isapprox (astrom. eral, 2.617608909189066140 , atol= 1e-12 )
467+ @test isapprox (astrom. eral, 2.617608909189664000 , atol= 1e-12 )
468468 @test isapprox (astrom. refa, 0.2014187785940396921e-3 , atol= 1e-15 )
469469 @test isapprox (astrom. refb, - 0.2361408314943696227e-6 , atol= 1e-18 )
470470end
@@ -573,11 +573,11 @@ let
573573 utc1, utc2, dut1, elong, phi, hm, xp, yp,
574574 phpa, tc, rh, wl)
575575
576- @test isapprox (aob, 0.09251774485385390973 , atol= 1e-12 )
577- @test isapprox (zob, 1.407661405256671703 , atol= 1e-12 )
578- @test isapprox (hob, - 0.09265154431430045141 , atol= 1e-12 )
579- @test isapprox (dob, 0.1716626560074556029 , atol= 1e-12 )
580- @test isapprox (rob, 2.710260453503366591 , atol= 1e-12 )
576+ @test isapprox (aob, 0.9251774485485515207e-1 , atol= 1e-12 )
577+ @test isapprox (zob, 1.407661405256499357 , atol= 1e-12 )
578+ @test isapprox (hob, - 0.9265154431529724692e-1 , atol= 1e-12 )
579+ @test isapprox (dob, 0.1716626560072526200 , atol= 1e-12 )
580+ @test isapprox (rob, 2.710260453504961012 , atol= 1e-12 )
581581 @test isapprox (eo, - 0.003020548354802412839 , atol= 1e-14 )
582582 @test j == 0
583583end
@@ -653,11 +653,11 @@ let
653653 j, aob, zob, hob, dob, rob = iauAtio13 (ri, di, utc1, utc2, dut1, elong, phi, hm,
654654 xp, yp, phpa, tc, rh, wl)
655655
656- @test isapprox (aob, 0.09233952224794989993 , atol= 1e-12 )
657- @test isapprox (zob, 1.407758704513722461 , atol= 1e-12 )
658- @test isapprox (hob, - 0.09247619879782006106 , atol= 1e-12 )
659- @test isapprox (dob, 0.1717653435758265198 , atol= 1e-12 )
660- @test isapprox (rob, 2.710085107986886201 , atol= 1e-12 )
656+ @test isapprox (aob, 0.9233952224895122499e-1 , atol= 1e-12 )
657+ @test isapprox (zob, 1.407758704513549991 , atol= 1e-12 )
658+ @test isapprox (hob, - 0.9247619879881698140e-1 , atol= 1e-12 )
659+ @test isapprox (dob, 0.1717653435756234676 , atol= 1e-12 )
660+ @test isapprox (rob, 2.710085107988480746 , atol= 1e-12 )
661661 @test j == 0
662662end
663663
@@ -682,11 +682,11 @@ let
682682
683683 aob, zob, hob, dob, rob = iauAtioq (ri, di, astrom)
684684
685- @test isapprox (aob, 0.09233952224794989993 , atol= 1e-12 )
686- @test isapprox (zob, 1.407758704513722461 , atol= 1e-12 )
687- @test isapprox (hob, - 0.09247619879782006106 , atol= 1e-12 )
688- @test isapprox (dob, 0.1717653435758265198 , atol= 1e-12 )
689- @test isapprox (rob, 2.710085107986886201 , atol= 1e-12 )
685+ @test isapprox (aob, 0.9233952224895122499e-1 , atol= 1e-12 )
686+ @test isapprox (zob, 1.407758704513549991 , atol= 1e-12 )
687+ @test isapprox (hob, - 0.9247619879881698140e-1 , atol= 1e-12 )
688+ @test isapprox (dob, 0.1717653435756234676 , atol= 1e-12 )
689+ @test isapprox (rob, 2.710085107988480746 , atol= 1e-12 )
690690end
691691
692692let
@@ -707,24 +707,24 @@ let
707707 ob2 = 0.1717653435758265198
708708 j, rc, dc = iauAtoc13 (' R' , ob1, ob2, utc1, utc2, dut1,
709709 elong, phi, hm, xp, yp, phpa, tc, rh, wl)
710- @test isapprox (rc, 2.709956744660731630 , atol= 1e-12 )
711- @test isapprox (dc, 0.1741696500896438967 , atol= 1e-12 )
710+ @test isapprox (rc, 2.709956744659136129 , atol= 1e-12 )
711+ @test isapprox (dc, 0.1741696500898471362 , atol= 1e-12 )
712712 @test j == 0
713713
714714 ob1 = - 0.09247619879782006106
715715 ob2 = 0.1717653435758265198
716716 j, rc, dc = iauAtoc13 (' H' , ob1, ob2, utc1, utc2, dut1,
717717 elong, phi, hm, xp, yp, phpa, tc, rh, wl)
718- @test isapprox (rc, 2.709956744660731630 , atol= 1e-12 )
719- @test isapprox (dc, 0.1741696500896438967 , atol= 1e-12 )
718+ @test isapprox (rc, 2.709956744659734086 , atol= 1e-12 )
719+ @test isapprox (dc, 0.1741696500898471362 , atol= 1e-12 )
720720 @test j == 0
721721
722722 ob1 = 0.09233952224794989993
723723 ob2 = 1.407758704513722461
724724 j, rc, dc = iauAtoc13 (' A' , ob1, ob2, utc1, utc2, dut1,
725725 elong, phi, hm, xp, yp, phpa, tc, rh, wl)
726- @test isapprox (rc, 2.709956744660731630 , atol= 1e-12 )
727- @test isapprox (dc, 0.1741696500896438970 , atol= 1e-12 )
726+ @test isapprox (rc, 2.709956744659734086 , atol= 1e-12 )
727+ @test isapprox (dc, 0.1741696500898471366 , atol= 1e-12 )
728728 @test j == 0
729729end
730730
@@ -746,24 +746,24 @@ let
746746 ob2 = 0.1717653435758265198
747747 j, ri, di = iauAtoi13 (' R' , ob1, ob2, utc1, utc2, dut1,
748748 elong, phi, hm, xp, yp, phpa, tc, rh, wl)
749- @test isapprox (ri, 2.710121574449135955 , atol= 1e-12 )
750- @test isapprox (di, 0.1729371839114567725 , atol= 1e-12 )
749+ @test isapprox (ri, 2.710121574447540810 , atol= 1e-12 )
750+ @test isapprox (di, 0.1729371839116608778 , atol= 1e-12 )
751751 @test j == 0
752752
753753 ob1 = - 0.09247619879782006106
754754 ob2 = 0.1717653435758265198
755755 j, ri, di = iauAtoi13 (' H' , ob1, ob2, utc1, utc2, dut1,
756756 elong, phi, hm, xp, yp, phpa, tc, rh, wl)
757- @test isapprox (ri, 2.710121574449135955 , atol= 1e-12 )
758- @test isapprox (di, 0.1729371839114567725 , atol= 1e-12 )
757+ @test isapprox (ri, 2.710121574448138676 , atol= 1e-12 )
758+ @test isapprox (di, 0.1729371839116608778 , atol= 1e-12 )
759759 @test j == 0
760760
761761 ob1 = 0.09233952224794989993
762762 ob2 = 1.407758704513722461
763763 j, ri, di = iauAtoi13 (' A' , ob1, ob2, utc1, utc2, dut1,
764764 elong, phi, hm, xp, yp, phpa, tc, rh, wl)
765- @test isapprox (ri, 2.710121574449135955 , atol= 1e-12 )
766- @test isapprox (di, 0.1729371839114567728 , atol= 1e-12 )
765+ @test isapprox (ri, 2.710121574448138676 , atol= 1e-12 )
766+ @test isapprox (di, 0.1729371839116608781 , atol= 1e-12 )
767767 @test j == 0
768768end
769769
@@ -786,20 +786,20 @@ let
786786 ob1 = 2.710085107986886201
787787 ob2 = 0.1717653435758265198
788788 ri, di = iauAtoiq (' R' , ob1, ob2, astrom)
789- @test isapprox (ri, 2.710121574449135955 , atol= 1e-12 )
790- @test isapprox (di, 0.1729371839114567725 , atol= 1e-12 )
789+ @test isapprox (ri, 2.710121574447540810 , atol= 1e-12 )
790+ @test isapprox (di, 0.17293718391166087785 , atol= 1e-12 )
791791
792792 ob1 = - 0.09247619879782006106
793793 ob2 = 0.1717653435758265198
794794 ri, di = iauAtoiq (' H' , ob1, ob2, astrom)
795- @test isapprox (ri, 2.710121574449135955 , atol= 1e-12 )
796- @test isapprox (di, 0.1729371839114567725 , atol= 1e-12 )
795+ @test isapprox (ri, 2.710121574448138676 , atol= 1e-12 )
796+ @test isapprox (di, 0.1729371839116608778 , atol= 1e-12 )
797797
798798 ob1 = 0.09233952224794989993
799799 ob2 = 1.407758704513722461
800800 ri, di = iauAtoiq (' A' , ob1, ob2, astrom)
801- @test isapprox (ri, 2.710121574449135955 , atol= 1e-12 )
802- @test isapprox (di, 0.1729371839114567728 , atol= 1e-12 )
801+ @test isapprox (ri, 2.710121574448138676 , atol= 1e-12 )
802+ @test isapprox (di, 0.1729371839116608781 , atol= 1e-12 )
803803end
804804
805805let
@@ -1544,7 +1544,7 @@ let
15441544
15451545 @test isapprox (rh, 1.767794226299947632 , atol= 1e-14 )
15461546 @test isapprox (dh, - 0.2917516070530391757 , atol= 1e-14 )
1547- @test isapprox (drh, - 0.19618741256057224e -6 ,atol= 1e-19 )
1547+ @test isapprox (drh, - 0.1961874125605721270e -6 ,atol= 1e-19 )
15481548 @test isapprox (ddh, - 0.58459905176693911e-5 , atol= 1e-19 )
15491549 @test isapprox (pxh, 0.37921 , atol= 1e-14 )
15501550 @test isapprox (rvh, - 7.6000000940000254 , atol= 1e-11 )
@@ -3040,9 +3040,9 @@ let
30403040
30413041 upv = iauPvu (2920.0 , pv)
30423042
3043- @test isapprox (upv[1 , 1 ], 126656.7598605317105 , atol= 1e-12 )
3044- @test isapprox (upv[1 , 2 ], 2118.531271155726332 , atol= 1e-12 )
3045- @test isapprox (upv[1 , 3 ], - 245216.5048590656190 , atol= 1e-12 )
3043+ @test isapprox (upv[1 , 1 ], 126656.7598605317105 , atol= 1e-6 )
3044+ @test isapprox (upv[1 , 2 ], 2118.531271155726332 , atol= 1e-8 )
3045+ @test isapprox (upv[1 , 3 ], - 245216.5048590656190 , atol= 1e-6 )
30463046
30473047 @test isapprox (upv[2 , 1 ], - 0.4051854035740713039e-2 , atol= 1e-12 )
30483048 @test isapprox (upv[2 , 2 ], - 0.6253919754866175788e-2 , atol= 1e-12 )
@@ -3055,9 +3055,9 @@ let
30553055
30563056 p = iauPvup (2920.0 , pv)
30573057
3058- @test isapprox (p[1 ], 126656.7598605317105 , atol= 1e-12 )
3059- @test isapprox (p[2 ], 2118.531271155726332 , atol= 1e-12 )
3060- @test isapprox (p[3 ], - 245216.5048590656190 , atol= 1e-12 )
3058+ @test isapprox (p[1 ], 126656.7598605317105 , atol= 1e-6 )
3059+ @test isapprox (p[2 ], 2118.531271155726332 , atol= 1e-8 )
3060+ @test isapprox (p[3 ], - 245216.5048590656190 , atol= 1e-6 )
30613061end
30623062
30633063let
0 commit comments