
    R1i#                    V   d Z ddlmZmZmZmZ ddlmZmZm	Z	m
Z
mZmZ ddlmZmZ ddlmZmZ ddlmZ ddlmZmZ dd	lmZ d
Zeez  ZdZdZdZdez  Zd Z d Z!d Z"d Z#d Z$d Z%d Z&d Z'd Z(d Z)d Z*d Z+d Z,d Z-d Z.d Z/d Z0d  Z1d! Z2d" Z3d# Z4d$ Z5d% Z6d& Z7d' Z8d( Z9d) Z:d* Z;d+ Z<d, Z=d- Z>d. Z?d/ Z@d0 ZAd1 ZBd2 ZCd3 ZDd4 ZEd5 ZFd6 ZGd7 ZHd8 ZId9 ZJd: ZKd; ZLd< ZMd= ZNd> ZOd? ZPd@ ZQdA ZRdB ZSdC ZTdD ZUdE ZVdF ZWdG ZXdH ZYdI ZZdJ Z[dK Z\dL Z]dM Z^dN Z_dO Z`dP ZadQ ZbdR ZcdS ZddT ZedU ZfdV ZgdW ZhdX ZidY ZjdZ Zkd[ Zld\ Zmd] Znd^ Zod_ Zpd` Zqda Zrdb Zsdc Ztdd Zude Zvdf Zwdg Zxdh Zydi Zzdj Z{dk Z|dl Z}dm Z~dn Zdo Zdp Zdq Zdr Zds Zdt Zdu Zdv Zdw Zdx Zdy Zdz Zd{ Zd| Zd} Zd~ Zd Zd Zd Zd Zd Zd Zd Zd Zd Zd Zd Zd Zd Zd Zd Zd Zd Zd Zd Zd Zd Zd Zd Zd Zd Zd Zd Zd Zd Zd Zd Zd Zd Zd Zd Zd Zd Zd Zd Zd Zd Zd Zd Zd Zd Zd Zd Zd Zd Zd Zd Zd Zd Zd Zd Zd Zd Zd Zd Zd Zd Zd Zd Zd Zd Zd Zd Zd ZdÄ ZdĄ Zdń ZdƄ ZdǄ ZdȄ ZdɄ Zdʄ Zd˄ Zd̄ Zd̈́ Zd΄ Zdτ ZdЄ Zdф Zd҄ Zdӄ ZdԄ ZdՄ Zy)zBAuto-generated accuracy tests vs NOVAS (see build_novas_tests.py).    )absarrayeinsummax)earthlibframelibnutationlibprecessionlibstarlibtimelib)Toposload)AU_KMAU_M)	hipparcos)BytesIO	length_of   )low_precision_ERAga]]աAg)QΠE>g߼xV42?g!tk>g      ?c               #   8   K   t        j                          y w)N)r   	timescale     e/home/cursorai/projects/iching/venv/lib/python3.12/site-packages/skyfield/tests/test_against_novas.pytsr      s     
..
s   c                     t        | d      st        |d      rt        t        | |z
              |k  sJ y t        | |z
        |k  sJ y )Nshape)hasattrr   r   )valueexpected_valueepsilons      r   comparer"      sJ    ug'.'"B3u~-./7:::5>)*g555r   c               #   &   K   t        d       y w)N	de405.bspr   r   r   r   de405r&      s     
{
s   c               #   0   K   t        d      } | d    y w)Nr$   i  r%   )ephs    r   earthr)       s     
{
C
c(Ns   c                 t    | j                   | j                  z
  }| j                   | _        || _        d| _        y )N        )tdbttwholett_fractiontdb_fraction)tdeltas     r   reduce_precisionr3   $   s.     EEADDLEeeAGAMANr   c                  X    t        t        j                  d      dz  ddt        z         y )NDDDsBA     v@gV^rJjn@ư>r"   r   earth_rotation_angle	arcsecondr   r   r   test_earth_rotation_angle_date0r;   -   s(    H))*;<uDFWy "r   c                  X    t        t        j                  d      dz  ddt        z         y )N   OBAr6   gKam@r7   r8   r   r   r   test_earth_rotation_angle_date1r>   1   '    H)))4u<>Oy "r   c                  X    t        t        j                  d      dz  ddt        z         y )N   ,BAr6   g_^q@r7   r8   r   r   r   test_earth_rotation_angle_date2rB   5   s&    H)))4u<oy "r   c                  X    t        t        j                  d      dz  ddt        z         y )N   @2BAr6   g8t@r7   r8   r   r   r   test_earth_rotation_angle_date3rE   9   r?   r   c                     t        t        j                  | j                  d            t	        d      dt
        z         y )Nr5   )gR=xq7@g&ݒr7@gNc?ga{@g΅!@h㈵>r"   r	   
earth_tilttdb_jdr   r:   r   s    r   test_earth_tilt_date0rL   =   s<    K""299->#?@uvx  CL  yLMr   c                     t        t        j                  | j                  d            t	        d      dt
        z         y )Nr=   )g3
p7@g2q7@gL0Q?gA'@gaXJ@rG   rH   rK   s    r   test_earth_tilt_date1rN   A   s;    K""299Y#78uvx  CL  yLMr   c                     t        t        j                  | j                  d            t	        d      dt
        z         y )NrA   )gtp7@gǖp7@gCgS.+gMrG   rH   rK   s    r   test_earth_tilt_date2rP   E   s?    K""299Y#78wx  {B  EN  {NOr   c                     t        t        j                  | j                  d            t	        d      dt
        z         y )NrD   )gT!p7@g;o7@g_M;MD?g/@g2M)
rG   rH   rK   s    r   test_earth_tilt_date3rR   I   s;    K""299Y#78stv}  AJ  wJKr   c                  d    t        t        j                  d      t        d      dt        z         y )Nr5   gr缉ؗҜ<r"   r	   -equation_of_the_equinoxes_complimentary_termsr   r:   r   r   r   8test_equation_of_the_equinoxes_complimentary_terms_date0rW   M   s*    KEEFWX)*,>,JLr   c                  d    t        t        j                  d      t        d      dt        z         y )Nr=   g:jݭGErT   rU   r   r   r   8test_equation_of_the_equinoxes_complimentary_terms_date1rY   Q   )    KEEiP()+=	+IKr   c                  d    t        t        j                  d      t        d      dt        z         y )NrA   gWJSE>rT   rU   r   r   r   8test_equation_of_the_equinoxes_complimentary_terms_date2r\   U   s)    KEEiP'(*<y*HJr   c                  d    t        t        j                  d      t        d      dt        z         y )NrD   g\w>GrT   rU   r   r   r   8test_equation_of_the_equinoxes_complimentary_terms_date3r^   Y   rZ   r   c                  X    t        t        j                  j                  d      dd       y )N皙?g333333?g?)g鞙?g%<S33?g\?V瞯<)r"   r   ICRS_to_J2000dotr   r   r   test_forward_frame_tiere   ]   s!    H""&&79tv{|r   c                  l    t        t        j                  j                  j	                  d      dd       y )Nr`   )gI?g!*33?g<?rb   )r"   r   rc   Trd   r   r   r   test_reverse_frame_tierh   `   s'    H""$$((9;vx}~r   c                  d    t        t        j                  d      t        d      dt        z         y )Ngڀ|ӿ)gg͡i:gyz2gg[@&.!>r"   r	   fundamental_argumentsr   r:   r   r   r    test_fundamental_arguments_date0rm   c   7    K--.ABwx  {F  IR  {RSr   c                  d    t        t        j                  d      t        d      dt        z         y )Ng2)gV3NgTwgRoąg=@^Zgz@rj   rk   r   r   r    test_fundamental_arguments_date1rp   g   s7    K--.BCwx  {F  IR  {RSr   c                  d    t        t        j                  d      t        d      dt        z         y )Nr+   )g}-@ghR@g<"2?giN:@gVwu@rj   rk   r   r   r    test_fundamental_arguments_date2rr   k   s2    K--c2qrt  CL  uLMr   c                  d    t        t        j                  d      t        d      dt        z         y )Ng! WU0?)ggsn?g-@g):?g(?g=$ rj   rk   r   r   r    test_fundamental_arguments_date3rt   o   rn   r   c                  Z    t        t        j                  d      t        ddg      d       y )Nr5   g6m2xAg	&!AMbP?r"   r	   iau2000ar   r   r   r   test_iau2000a_date0ry   s   *    K  !23%'89:ECr   c                  Z    t        t        j                  d      t        ddg      d       y )Nr=   g%BAg故IArv   rw   r   r   r   test_iau2000a_date1r|   w   s)    K  +$&7895Br   c                  Z    t        t        j                  d      t        ddg      d       y )NrA   g6gY渂rv   rw   r   r   r   test_iau2000a_date2r~   {   s)    K  +&(;<=uFr   c                  Z    t        t        j                  d      t        ddg      d       y )NrD   gjAg<Prv   rw   r   r   r   test_iau2000a_date3r      s)    K  +$&9:;UDr   c                  Z    t        t        j                  d      t        ddg      d       y )Nr5   gGtpxAg2n Arv   r"   r	   iau2000br   r   r   r   test_iau2000b_date0r      rz   r   c                  Z    t        t        j                  d      t        ddg      d       y )Nr=   gָAgyIArv   r   r   r   r   test_iau2000b_date1r      )    K  +%'9:;UDr   c                  Z    t        t        j                  d      t        ddg      d       y )NrA   gz=g.־rv   r   r   r   r   test_iau2000b_date2r      s)    K  +&(:;<eEr   c                  Z    t        t        j                  d      t        ddg      d       y )NrD   gAgPp	Prv   r   r   r   r   test_iau2000b_date3r      r   r   c                  D    t        t        j                  d      dd       y )Nr5   gW孺@r+   r"   r	   mean_obliquityr   r   r   test_mean_obliquity_date0r      s    K&&'89c#r   c                  D    t        t        j                  d      dd       y )Nr=   g|8@r+   r   r   r   r   test_mean_obliquity_date1r          K&&y1s$r   c                  D    t        t        j                  d      dd       y )NrA   g#~֙@r+   r   r   r   r   test_mean_obliquity_date2r      s    K&&y1sr   c                  D    t        t        j                  d      dd       y )NrD   gK#w@r+   r   r   r   r   test_mean_obliquity_date3r      r   r   c                     t        j                  | j                  d            }t        d|g d      }t	        d|d       y )Nr5   ij...,j...->i...r`   )g_,?g.3?g\?+=r	   compute_nutationrJ   r   r"   r   matrixresults      r   test_nutation_date0r      s;    ))"))4E*FGF&@FHEr   c                     t        j                  | j                  d            }t        d|g d      }t	        d|d       y )Nr=   r   r`   )gX9?g}_G3?g{
?r   r   r   s      r   test_nutation_date1r      s:    ))"))I*>?F&@FGEr   c                     t        j                  | j                  d            }t        d|g d      }t	        d|d       y )NrA   r   r`   )g"k2?g(13?gَS?r   r   r   s      r   test_nutation_date2r      :    ))"))I*>?F&@FHEr   c                     t        j                  | j                  d            }t        d|g d      }t	        d|d       y )NrD   r   r`   )gHY-?g.زQ3?gt>X?r   r   r   s      r   test_nutation_date3r      r   r   c                  f    t        j                  d      } t        d| g d      }t        d|d       y )Nr5   r   r`   )gcx?g2[?g@ie?rb   r
   compute_precessionr   r"   r   r   s     r   test_precession_date0r      s2    --.?@F&@FGEr   c                  f    t        j                  d      } t        d| g d      }t        d|d       y )Nr=   r   r`   )g"?g46~)?g"d֔?rb   r   r   s     r   test_precession_date1r      s1    --i8F&@FHEr   c                  f    t        j                  d      } t        d| g d      }t        d|d       y )NrA   r   r`   )ra   g233333?g?rb   r   r   s     r   test_precession_date2r      s1    --i8F&@F9Er   c                  f    t        j                  d      } t        d| g d      }t        d|d       y )NrD   r   r`   )gJ]"?g??gmnhR?rb   r   r   s     r   test_precession_date3r      s1    --i8F&@FGEr   c                      t        j                  d      j                  d      } t        t	        j
                  |       dd       y )Nr+   delta_tr5   g020@vIh%<=r   r   tt_jdr"   r   sidereal_timejds    r   test_sidereal_time_on_date0r      s4    		$	*	*+<	=BH""2&(:EBr   c                      t        j                  d      j                  ddt        z  z         } t	               5  t        t        j                  |       dd       d d d        y # 1 sw Y   y xY w)N皙X@r   r5   g|F20@r   r   r   r   
one_secondr   r"   r   r   r   s    r   0test_sidereal_time_with_nonzero_delta_t_on_date0r      sP    		%	+	+,=z@Q,Q	RB		&&r*,>F 
		   !A$$A-c                      t        j                  d      j                  d      } t        t	        j
                  |       dd       y )Nr+   r   r=   gWݦ/@r   r   r   s    r   test_sidereal_time_on_date1r      3    		$	*	*9	5BH""2&(:EBr   c                      t        j                  d      j                  ddt        z  z         } t	               5  t        t        j                  |       dd       d d d        y # 1 sw Y   y xY w)Nr   r   r=   g%Wݦ/@r   r   r   s    r   0test_sidereal_time_with_nonzero_delta_t_on_date1r      sO    		%	+	+Iz8I,I	JB		&&r*,>F 
		r   c                      t        j                  d      j                  d      } t        t	        j
                  |       dd       y )Nr+   r   rA   g(2@r   r   r   s    r   test_sidereal_time_on_date2r      s3    		$	*	*9	5BH""2&(95Ar   c                      t        j                  d      j                  ddt        z  z         } t	               5  t        t        j                  |       dd       d d d        y # 1 sw Y   y xY w)Nr   r   rA   g{&(2@r   r   r   s    r   0test_sidereal_time_with_nonzero_delta_t_on_date2r      sO    		%	+	+Iz8I,I	JB		&&r*,=uE 
		r   c                      t        j                  d      j                  d      } t        t	        j
                  |       dd       y )Nr+   r   rD   g.p>6@r   r   r   s    r   test_sidereal_time_on_date3r      r   r   c                      t        j                  d      j                  ddt        z  z         } t	               5  t        t        j                  |       dd       d d d        y # 1 sw Y   y xY w)Nr   r   rD   g?(p>6@r   r   r   s    r   0test_sidereal_time_with_nonzero_delta_t_on_date3r      sO    		%	+	+Iz8I,I	JB		&&r*,<eD 
		r   c                      t        j                  dddddd      } t        | _        | j	                          t        | j                  dd	t        z         t        | j                  d
dt        z         y )N\"w>@"L+PV@\(F@     '=
ףp=@ffffff1ra_hoursdec_degreesra_mas_per_yeardec_mas_per_yearparallax_masradial_km_per_s)gy?6Ag)N
AgO[zAg     @@)gT {gsM>?g3R%rv   )	r   Star	OLD_AU_KMau_km_compute_vectorsr"   _position_aumeter_velocity_au_per_d)stars    r   test_star_vectorr      sj    <<,(-%)5BD DJDG%K D##P5Lr   c                  Z    t        j                  ddd      } t        | ddt        z         y )N
     r+   &.>r   
refractionr"   r:   rs    r   test_refraction0r     '    BD)AAsD9$%r   c                  Z    t        j                  ddd      } t        | ddt        z         y )Nr   r        @r+   r   r   r   s    r   test_refraction1r     '    BG,AAsD9$%r   c                  Z    t        j                  ddd      } t        | ddt        z         y )Nr      r   r+   r   r   r   s    r   test_refraction2r     r   r   c                  Z    t        j                  ddd      } t        | ddt        z         y )Nr   r   r   r+   r   r   r   s    r   test_refraction3r     r   r   c                  Z    t        j                  ddd      } t        | ddt        z         y )Nr   r   g|֌?r   r   r   s    r   test_refraction4r     (    BD)AA!4)#34r   c                  Z    t        j                  ddd      } t        | ddt        z         y )Nr   r   r   gN?r   r   r   s    r   test_refraction5r     (    BG,AA!4)#34r   c                  Z    t        j                  ddd      } t        | ddt        z         y )Nr   r   r   gIY6?r   r   r   s    r   test_refraction6r    r   r   c                  Z    t        j                  ddd      } t        | ddt        z         y )Nr   r   r   gr{K?r   r   r   s    r   test_refraction7r     r   r   c                  Z    t        j                  ddd      } t        | ddt        z         y )N   r   r   g>Iŵ ?r   r   r   s    r   test_refraction8r  $  s(    BD)AA"D9$45r   c                  Z    t        j                  ddd      } t        | ddt        z         y )Nr  r   r   gX?r   r   r   s    r   test_refraction9r	  (  s(    BG,AA"D9$45r   c                  Z    t        j                  ddd      } t        | ddt        z         y )Nr  r   r   g&Zpr?r   r   r   s    r   test_refraction10r  ,  s(    BD)AA#TI%56r   c                  Z    t        j                  ddd      } t        | ddt        z         y )Nr  r   r   gŋ<?r   r   r   s    r   test_refraction11r  0  s(    BG,AA#TI%56r   c                  Z    t        j                  ddd      } t        | ddt        z         y )N|V@r   r   r+   r   r   r   s    r   test_refraction12r  4  '    E2t,AAsD9$%r   c                  Z    t        j                  ddd      } t        | ddt        z         y )Nr  r   r   r+   r   r   r   s    r   test_refraction13r  8  '    E2w/AAsD9$%r   c                  Z    t        j                  ddd      } t        | ddt        z         y )Nr  r   r   r+   r   r   r   s    r   test_refraction14r  <  r  r   c                  Z    t        j                  ddd      } t        | ddt        z         y )Nr  r   r   r+   r   r   r   s    r   test_refraction15r  @  r  r   c                  Z    t        j                  ddd      } t        | ddt        z         y )Ni      $@     @g     Vr   r   refractr"   r:   alts    r   test_refract0r   D  s'    


3f
-CCy()r   c                  Z    t        j                  ddd      } t        | ddt        z         y )Nr  r  g       r   r  r  s    r   test_refract1r#  H  '    


2tV
,CCti'(r   c                  Z    t        j                  ddd      } t        | ddt        z         y )Nr   r  r  g f
ֿr   r  r  s    r   test_refract2r&  L  s(    


2tV
,CC%ti'78r   c                  Z    t        j                  ddd      } t        | ddt        z         y )Nr   r  r  g z-?r   r  r  s    r   test_refract3r(  P  s(    


1dF
+CC#TI%56r   c                  Z    t        j                  ddd      } t        | ddt        z         y )Nr   r  r  g@?r   r  r  s    r   test_refract4r*  T  (    


1dF
+CC"D9$45r   c                  Z    t        j                  ddd      } t        | ddt        z         y )N   r  r  g,	p	@r   r  r  s    r   test_refract5r.  X  r+  r   c                  Z    t        j                  ddd      } t        | ddt        z         y )N	   r  r  g@ 042"@r   r  r  s    r   test_refract6r1  \  r+  r   c                  Z    t        j                  ddd      } t        | ddt        z         y )NZ   r  r  g     V@r   r  r  s    r   test_refract7r4  `  r$  r   c                 X   t        j                  d      j                  d      }| t        ddd      z   }|j	                  |      j                  dd	
      }|j                  |      \  }}}t        |j                  ddt        z         t        |j                  ddt        z         y )Nr+   r   r5   	38.9215 N	77.0669 W      W@elevation_mgVD<Q@g~Ͽ<@alt_degrees
az_degreesepochGz(@r   p=
cL@r   r   r   r   at
from_altazradecr"   hoursr:   degreesr)   r   usnoaradecdistances          r   test_from_altaz_0rN  d  s    		$	*	*+<	=B5[d4 4D+<I[\Ab)BXBHHeTI-.CKKy 01r   c                 X   t        j                  d      j                  d      }| t        ddd      z   }|j	                  |      j                  dd	
      }|j                  |      \  }}}t        |j                  ddt        z         t        |j                  ddt        z         y )Nr+   r   r5   r6  r7  r8  r9  g91g
.e@r;  r>  r@  r   )\PrB  rH  s          r   test_from_altaz_1rQ  m  s    		$	*	*+<	=B5[d4 4D+>K]^Ab)BXBHHeTI-.CKK	!12r   c                 X   t        j                  d      j                  d      }| t        ddd      z   }|j	                  |      j                  dd	
      }|j                  |      \  }}}t        |j                  ddt        z         t        |j                  ddt        z         y )Nr+   r   r=   r6  r7  r8  r9  gh#]wP@gQp RA@r;  r>  r@  r   rA  rB  rH  s          r   test_from_altaz_2rS  v  s    		$	*	*9	5B5[d4 4D+;HZ[Ab)BXBHHeTI-.CKKy 01r   c                 X   t        j                  d      j                  d      }| t        ddd      z   }|j	                  |      j                  dd	
      }|j                  |      \  }}}t        |j                  ddt        z         t        |j                  ddt        z         y )Nr+   r   r=   r6  r7  r8  r9  gۡn2g}Me@r;  r>  r@  r   rP  rB  rH  s          r   test_from_altaz_3rU    s    		$	*	*9	5B5[d4 4D+=J\]Ab)BXBHHeTI-.CKK	!12r   c                 X   t        j                  d      j                  d      }| t        ddd      z   }|j	                  |      j                  dd	
      }|j                  |      \  }}}t        |j                  ddt        z         t        |j                  ddt        z         y )Nr+   r   rA   r6  r7  r8  r9  g]Q@gf[%Ht@r;  r>  r@  r   rA  rB  rH  s          r   test_from_altaz_4rW        		$	*	*9	5B5[d4 4D+<I[\Ab)BXBHHeTI-.CKKy 01r   c                 X   t        j                  d      j                  d      }| t        ddd      z   }|j	                  |      j                  dd	
      }|j                  |      \  }}}t        |j                  ddt        z         t        |j                  ddt        z         y )Nr+   r   rA   r6  r7  r8  r9  g,1gܿcg@r;  r>  r@  r   rP  rB  rH  s          r   test_from_altaz_5rZ    s    		$	*	*9	5B5[d4 4D+>K]^Ab)BXBHHeTI-.CKK	!12r   c                 X   t        j                  d      j                  d      }| t        ddd      z   }|j	                  |      j                  dd	
      }|j                  |      \  }}}t        |j                  ddt        z         t        |j                  ddt        z         y )Nr+   r   rD   r6  r7  r8  r9  g*®D@g\s@r;  r>  r@  r   rA  rB  rH  s          r   test_from_altaz_6r\    rX  r   c                 X   t        j                  d      j                  d      }| t        ddd      z   }|j	                  |      j                  dd	
      }|j                  |      \  }}}t        |j                  ddt        z         t        |j                  ddt        z         y )Nr+   r   rD   r6  r7  r8  r9  glwW4ZH=gޡi@r;  r>  r@  r   rP  rB  rH  s          r   test_from_altaz_7r^    s    		$	*	*9	5B5[d4 4D+>K\]Ab)BXBHHeTI-.CKK	!12r   c                  H    t        j                  d      } t        | dd       y )Nr5   g>rT   r   tdb_minus_ttr"   r   s    r   test_tdb_minus_tt_on_date0rc    s     !!"34FF+U3r   c                  H    t        j                  d      } t        | dd       y )Nr=   g1-R?rT   r`  rb  s    r   test_tdb_minus_tt_on_date1re        !!),FF)51r   c                  H    t        j                  d      } t        | dd       y )NrA   gDoM4,rT   r`  rb  s    r   test_tdb_minus_tt_on_date2rh    s    !!),FF*E2r   c                  H    t        j                  d      } t        | dd       y )NrD   g}wAUTrT   r`  rb  s    r   test_tdb_minus_tt_on_date3rj    rf  r   c                     |j                  d      }| d   j                  |      }t        |j                  j                  ddt
        z         t        |j                  j                  ddt
        z  t        z         y )NrA   r)   )g5ǿg3K!P?gx:?r   )gD4:wgLDgg^lTrG   )	r   rC  r"   positionaur   velocityau_per_dr   )r&   r   r1   es       r   test_position_and_velocityrq    so    
Ag!AAJJMMY[]`e[efAJJ!hjnqvjv  zD  kD  Er   c                 ,   |j                  d      }t        |       | d   j                  |      }| d   }t        ||j                  |      z
  j                  j
                        }t        |t        z  ddt        z         |j                  |      }|j                         \  }}}t        |j                  ddt        z         t        |j                  dd	t        z         |j                         }	|	j                         \  }}}t        |j                  d
dt        z         t        |j                  dd	t        z         |	j                  d      \  }}}t        |j                  ddt        z         t        |j                  dd	t        z         y )Nr5   r)   mercuryR>?y&1?@-C6*?U6@-C6?Ie@[U6@dater>  6ʏ@	}j6@r   r3   rC  r   rl  rm  r"   OLD_AUr   observerE  rF  ra_arcsecondrG  r:   apparent
r&   r   r1   rp  prM  astrometricrK  rL  r  s
             r   test_mercury_geocentric_date0r    sM   
"#AQg!AiA!add1g+//223HHv155=A))A,K#))+BXBHH',)>?CKK+Vi-?@##%H (BXBHH',)>?CKK+Vi-?@ V4BXBHH',)>?CKK+Vi-?@r   c                 ,   |j                  d      }t        |       | d   j                  |      }| d   }t        ||j                  |      z
  j                  j
                        }t        |t        z  ddt        z         |j                  |      }|j                         \  }}}t        |j                  ddt        z         t        |j                  dd	t        z         |j                         }	|	j                         \  }}}t        |j                  d
dt        z         t        |j                  dd	t        z         |	j                  d      \  }}}t        |j                  ddt        z         t        |j                  dd	t        z         y )Nr=   r)   rs  #)$?ru  z]@rw  Ci&@ry  ,@B&@r|  r>  熺Ȑ@$4bj&@r  r  s
             r   test_mercury_geocentric_date1r    sL   
AQg!AiA!add1g+//223HHv155=A))A,K#))+BXBHH(&<*?@CKK)6I+=>##%H (BXBHH(&<*?@CKK+Vi-?@ V4BXBHH(&<*?@CKK+Vi-?@r   c                 ,   |j                  d      }t        |       | d   j                  |      }| d   }t        ||j                  |      z
  j                  j
                        }t        |t        z  ddt        z         |j                  |      }|j                         \  }}}t        |j                  ddt        z         t        |j                  dd	t        z         |j                         }	|	j                         \  }}}t        |j                  d
dt        z         t        |j                  dd	t        z         |	j                  d      \  }}}t        |j                  ddt        z         t        |j                  dd	t        z         y )NrA   r)   rs  ?ru  <#2@rw  熢Sjk8ry  7_6ju#2@#k8r|  r>  #nRb#2@z9k8r  r  s
             r   test_mercury_geocentric_date2r    sL   
AQg!AiA!add1g+//223HHv155=A))A,K#))+BXBHH',)>?CKK+Vi-?@##%H (BXBHH',)>?CKK,fy.@A V4BXBHH(&<*?@CKK,fy.@Ar   c                 ,   |j                  d      }t        |       | d   j                  |      }| d   }t        ||j                  |      z
  j                  j
                        }t        |t        z  ddt        z         |j                  |      }|j                         \  }}}t        |j                  ddt        z         t        |j                  dd	t        z         |j                         }	|	j                         \  }}}t        |j                  d
dt        z         t        |j                  dd	t        z         |	j                  d      \  }}}t        |j                  ddt        z         t        |j                  dd	t        z         y )NrD   r)   rs  L+?ru  .􂿄"@rw  k0@ry  dȌU"@ս0@r|  r>  zcMy"@O0@r  r  s
             r   test_mercury_geocentric_date3r  
  sL   
AQg!AiA!add1g+//223HHv155=A))A,K#))+BXBHH',)>?CKK*FY,>?##%H (BXBHH',)>?CKK+Vi-?@ V4BXBHH',)>?CKK+Vi-?@r   c                 0   |j                  g d      }t        |       | d   j                  |      }| d   }t        ||j                  |      z
  j                  j
                        }t        |t        z  ddt        z         |j                  |      }|j                         \  }}}t        |j                  ddt        z         t        |j                  dd	t        z         |j                         }	|	j                         \  }}}t        |j                  d
dt        z         t        |j                  dd	t        z         |	j                  d      \  }}}t        |j                  ddt        z         t        |j                  dd	t        z         y )Nr5   r=   rA   rD   r)   rs  )rt  r  r  r  ru  )rv  r  r  r  rw  )rx  r  r  r  ry  )rz  r  r  r  )r{  r  r  r  r|  r>  )r}  r  r  r  )r~  r  r  r  r  r  s
             r   test_mercury_geocentric_date4r  "  s\   
EFAQg!AiA!add1g+//223HHvoqvy~q~))A,K#))+BXBHHceknzez{CKKfhnqzhz{##%H (BXBHHceknzez{CKKjlru~l~ V4BXBHHdflo{f{|CKKjlru~l~r   c                 ,   |j                  d      }t        |       | d   j                  |      }| d   }t        ||j                  |      z
  j                  j
                        }t        |t        z  ddt        z         |j                  |      }|j                         \  }}}t        |j                  ddt        z         t        |j                  dd	t        z         |j                         }	|	j                         \  }}}t        |j                  d
dt        z         t        |j                  dd	t        z         |	j                  d      \  }}}t        |j                  ddt        z         t        |j                  dd	t        z         y )Nr5   r)   venus~d
?ru  ,'@rw  54@ry  S'e@SRq54@r|  r>  P*@箧9*4@r  r  s
             r   test_venus_geocentric_date0r  :  sL   
"#AQg!AgA!add1g+//223HHv155=A))A,K#))+BXBHH',)>?CKK+Vi-?@##%H (BXBHH',)>?CKK+Vi-?@ V4BXBHH&(=>CKK+Vi-?@r   c                 ,   |j                  d      }t        |       | d   j                  |      }| d   }t        ||j                  |      z
  j                  j
                        }t        |t        z  ddt        z         |j                  |      }|j                         \  }}}t        |j                  ddt        z         t        |j                  dd	t        z         |j                         }	|	j                         \  }}}t        |j                  d
dt        z         t        |j                  dd	t        z         |	j                  d      \  }}}t        |j                  ddt        z         t        |j                  dd	t        z         y )Nr=   r)   r  E|r#?ru  pǖ?rw  +P@ry  Ŏ?JN@r|  r>  et?!@r  r  s
             r   test_venus_geocentric_date1r  R  sK   
AQg!AgA!add1g+//223HHv155=A))A,K#))+BXBHH',)>?CKK)6I+=>##%H (BXBHH(&<*?@CKK*FY,>? V4BXBHH(&<*?@CKK*FY,>?r   c                 ,   |j                  d      }t        |       | d   j                  |      }| d   }t        ||j                  |      z
  j                  j
                        }t        |t        z  ddt        z         |j                  |      }|j                         \  }}}t        |j                  ddt        z         t        |j                  dd	t        z         |j                         }	|	j                         \  }}}t        |j                  d
dt        z         t        |j                  dd	t        z         |	j                  d      \  }}}t        |j                  ddt        z         t        |j                  dd	t        z         y )NrA   r)   r  犹u3?ru  $t/@rw  1ls2ry  o/@hls2r|  r>  9N/@0?r2r  r  s
             r   test_venus_geocentric_date2r  j  sK   
AQg!AgA!add1g+//223HHv155=A))A,K#))+BXBHH(&<*?@CKK,fy.@A##%H (BXBHH(&<*?@CKK,fy.@A V4BXBHH(&<*?@CKK+Vi-?@r   c                 ,   |j                  d      }t        |       | d   j                  |      }| d   }t        ||j                  |      z
  j                  j
                        }t        |t        z  ddt        z         |j                  |      }|j                         \  }}}t        |j                  ddt        z         t        |j                  dd	t        z         |j                         }	|	j                         \  }}}t        |j                  d
dt        z         t        |j                  dd	t        z         |	j                  d      \  }}}t        |j                  ddt        z         t        |j                  dd	t        z         y )NrD   r)   r  ȏ-
?ru  PP̳@rw  ɴc 3@ry  RE;@i3@r|  r>  \@:3@r  r  s
             r   test_venus_geocentric_date3r    sK   
AQg!AgA!add1g+//223HHv155=A))A,K#))+BXBHH',)>?CKK+Vi-?@##%H (BXBHH',)>?CKK*FY,>? V4BXBHH',)>?CKK*FY,>?r   c                 0   |j                  g d      }t        |       | d   j                  |      }| d   }t        ||j                  |      z
  j                  j
                        }t        |t        z  ddt        z         |j                  |      }|j                         \  }}}t        |j                  ddt        z         t        |j                  dd	t        z         |j                         }	|	j                         \  }}}t        |j                  d
dt        z         t        |j                  dd	t        z         |	j                  d      \  }}}t        |j                  ddt        z         t        |j                  dd	t        z         y )Nr  r)   r  )r  r  r  r  ru  )r  r  r  r  rw  )r  r  r  r  ry  )r  r  r  r  )r  r  r  r  r|  r>  )r  r  r  r  )r  r  r  r  r  r  s
             r   test_venus_geocentric_date4r    s[   
EFAQg!AgA!add1g+//223HHvoqvy~q~))A,K#))+BXBHHceknzez{CKKhjps|j|}##%H (BXBHHdflo{f{|CKKhjps|j|} V4BXBHHceknzez{CKKgior{i{|r   c                 ,   |j                  d      }t        |       | d   j                  |      }| d   }t        ||j                  |      z
  j                  j
                        }t        |t        z  ddt        z         |j                  |      }|j                         \  }}}t        |j                  ddt        z         t        |j                  dd	t        z         |j                         }	|	j                         \  }}}t        |j                  d
dt        z         t        |j                  dd	t        z         |	j                  d      \  }}}t        |j                  ddt        z         t        |j                  dd	t        z         y )Nr5   r)   mars<=C?ru  ֗0@rw  ͋h 8ry  ><0@f 8r|  r>  d/@I8r  r  s
             r   test_mars_geocentric_date0r    sL   
"#AQg!AfA!add1g+//223HHv155=A))A,K#))+BXBHH&(=>CKK,fy.@A##%H (BXBHH',)>?CKK,fy.@A V4BXBHH',)>?CKK,fy.@Ar   c                 ,   |j                  d      }t        |       | d   j                  |      }| d   }t        ||j                  |      z
  j                  j
                        }t        |t        z  ddt        z         |j                  |      }|j                         \  }}}t        |j                  ddt        z         t        |j                  dd	t        z         |j                         }	|	j                         \  }}}t        |j                  d
dt        z         t        |j                  dd	t        z         |	j                  d      \  }}}t        |j                  ddt        z         t        |j                  dd	t        z         y )Nr=   r)   r  ޞ_N?ru  .Gg7@rw  ?Klry  n~7@|"}{/r|  r>  5n7@D+r  r  s
             r   test_mars_geocentric_date1r    sK   
AQg!AfA!add1g+//223HHv0%%-@))A,K#))+BXBHH(&<*?@CKK,fy.@A##%H (BXBHH(&<*?@CKK*FY,>? V4BXBHH(&<*?@CKK+Vi-?@r   c                 ,   |j                  d      }t        |       | d   j                  |      }| d   }t        ||j                  |      z
  j                  j
                        }t        |t        z  ddt        z         |j                  |      }|j                         \  }}}t        |j                  ddt        z         t        |j                  dd	t        z         |j                         }	|	j                         \  }}}t        |j                  d
dt        z         t        |j                  dd	t        z         |	j                  d      \  }}}t        |j                  ddt        z         t        |j                  dd	t        z         y )NrA   r)   r  I?ru  稠(6@rw  \*ry  PI6@@]*r|  r>  6@瓧Rs]*r  r  s
             r   test_mars_geocentric_date2r    sK   
AQg!AfA!add1g+//223HHv155=A))A,K#))+BXBHH(&<*?@CKK+Vi-?@##%H (BXBHH',)>?CKK,fy.@A V4BXBHH(&<*?@CKK,fy.@Ar   c                 ,   |j                  d      }t        |       | d   j                  |      }| d   }t        ||j                  |      z
  j                  j
                        }t        |t        z  ddt        z         |j                  |      }|j                         \  }}}t        |j                  ddt        z         t        |j                  dd	t        z         |j                         }	|	j                         \  }}}t        |j                  d
dt        z         t        |j                  dd	t        z         |	j                  d      \  }}}t        |j                  ddt        z         t        |j                  dd	t        z         y )NrD   r)   r  ;ܹ\C?ru  p+@rw  yb>(ry  ;+@ٿ+\>(r|  r>  ,+@=5@^(r  r  s
             r   test_mars_geocentric_date3r    sK   
AQg!AfA!add1g+//223HHv155=A))A,K#))+BXBHH(&<*?@CKK,fy.@A##%H (BXBHH(&<*?@CKK,fy.@A V4BXBHH&(=>CKK,fy.@Ar   c                 0   |j                  g d      }t        |       | d   j                  |      }| d   }t        ||j                  |      z
  j                  j
                        }t        |t        z  ddt        z         |j                  |      }|j                         \  }}}t        |j                  ddt        z         t        |j                  dd	t        z         |j                         }	|	j                         \  }}}t        |j                  d
dt        z         t        |j                  dd	t        z         |	j                  d      \  }}}t        |j                  ddt        z         t        |j                  dd	t        z         y )Nr  r)   r  )r  r  r  r  ru  )r  r  r  r  rw  )r  r  r  r  ry  )r  r  r  r  )r  r  r  r  r|  r>  )r  r  r  r  )r  r  r  r  r  r  s
             r   test_mars_geocentric_date4r    st   
EFAQg!AfA!add1g+//223HHvnpux}p}~))A,K#))+BXBHHdflo{f{|CKKlnt  xA  oA  B##%H (BXBHHdflo{f{|CKKkmsvm  A V4BXBHHceknzez{CKKlnt  xA  oA  Br   c                 ,   |j                  d      }t        |       | d   j                  |      }| d   }t        ||j                  |      z
  j                  j
                        }t        |t        z  ddt        z         |j                  |      }|j                         \  }}}t        |j                  ddt        z         t        |j                  dd	t        z         |j                         }	|	j                         \  }}}t        |j                  d
dt        z         t        |j                  dd	t        z         |	j                  d      \  }}}t        |j                  ddt        z         t        |j                  dd	t        z         y )Nr5   r)   jupiter barycenter2]y]@ru  Z{K5(@rw  ܃!?ry  #75(@a?r|  r>  U0=C'(@Z!8J?r  r  s
             r   (test_jupiter_barycenter_geocentric_date0r  *  sN   
"#AQg!A"#A!add1g+//223HHv155=A))A,K#))+BXBHH(&<*?@CKK+Vi-?@##%H (BXBHH(&<*?@CKK+Vi-?@ V4BXBHH',)>?CKK+Vi-?@r   c                 ,   |j                  d      }t        |       | d   j                  |      }| d   }t        ||j                  |      z
  j                  j
                        }t        |t        z  ddt        z         |j                  |      }|j                         \  }}}t        |j                  ddt        z         t        |j                  dd	t        z         |j                         }	|	j                         \  }}}t        |j                  d
dt        z         t        |j                  dd	t        z         |	j                  d      \  }}}t        |j                  ddt        z         t        |j                  dd	t        z         y )Nr=   r)   r  z
W4@ru  ^E@rw  qV/+7@ry  |5@nyR`+7@r|  r>  .$@Ґl.7@r  r  s
             r   (test_jupiter_barycenter_geocentric_date1r  B  sM   
AQg!A"#A!add1g+//223HHv0%%-@))A,K#))+BXBHH',)>?CKK+Vi-?@##%H (BXBHH',)>?CKK+Vi-?@ V4BXBHH',)>?CKK+Vi-?@r   c                 ,   |j                  d      }t        |       | d   j                  |      }| d   }t        ||j                  |      z
  j                  j
                        }t        |t        z  ddt        z         |j                  |      }|j                         \  }}}t        |j                  ddt        z         t        |j                  dd	t        z         |j                         }	|	j                         \  }}}t        |j                  d
dt        z         t        |j                  dd	t        z         |	j                  d      \  }}}t        |j                  ddt        z         t        |j                  dd	t        z         y )NrA   r)   r  z<|@ru  _v?rw  ﺂ/1!@ry  z)tqv?笧O1!@r|  r>  *nu?9A0!@r  r  s
             r   (test_jupiter_barycenter_geocentric_date2r  Z  sM   
AQg!A"#A!add1g+//223HHv0%%-@))A,K#))+BXBHH(&<*?@CKK*FY,>?##%H (BXBHH(&<*?@CKK)6I+=> V4BXBHH(&<*?@CKK*FY,>?r   c                 ,   |j                  d      }t        |       | d   j                  |      }| d   }t        ||j                  |      z
  j                  j
                        }t        |t        z  ddt        z         |j                  |      }|j                         \  }}}t        |j                  ddt        z         t        |j                  dd	t        z         |j                         }	|	j                         \  }}}t        |j                  d
dt        z         t        |j                  dd	t        z         |	j                  d      \  }}}t        |j                  ddt        z         t        |j                  dd	t        z         y )NrD   r)   r  篽C@ru  ݰtؖJ@rw  jT	f5@ry  ؂J@)?5b5@r|  r>  pW@Zǳ5@r  r  s
             r   (test_jupiter_barycenter_geocentric_date3r  r  sM   
AQg!A"#A!add1g+//223HHv0%%-@))A,K#))+BXBHH',)>?CKK+Vi-?@##%H (BXBHH',)>?CKK*FY,>? V4BXBHH',)>?CKK*FY,>?r   c                 0   |j                  g d      }t        |       | d   j                  |      }| d   }t        ||j                  |      z
  j                  j
                        }t        |t        z  ddt        z         |j                  |      }|j                         \  }}}t        |j                  ddt        z         t        |j                  dd	t        z         |j                         }	|	j                         \  }}}t        |j                  d
dt        z         t        |j                  dd	t        z         |	j                  d      \  }}}t        |j                  ddt        z         t        |j                  dd	t        z         y )Nr  r)   r  )r  r  r  r  ru  )r  r   r	  r  rw  )r  r  r
  r  ry  )r  r  r  r  )r  r  r  r  r|  r>  )r  r  r  r  )r  r  r  r  r  r  s
             r   (test_jupiter_barycenter_geocentric_date4r    s]   
EFAQg!A"#A!add1g+//223HHvlnsv{n{|))A,K#))+BXBHHdflo{f{|CKKhjps|j|}##%H (BXBHHdflo{f{|CKKfhnqzhz{ V4BXBHHceknzez{CKKgior{i{|r   c                 ,   |j                  d      }t        |       | d   j                  |      }| d   }t        ||j                  |      z
  j                  j
                        }t        |t        z  ddt        z         |j                  |      }|j                         \  }}}t        |j                  ddt        z         t        |j                  dd	t        z         |j                         }	|	j                         \  }}}t        |j                  d
dt        z         t        |j                  dd	t        z         |	j                  d      \  }}}t        |j                  ddt        z         t        |j                  dd	t        z         y )Nr5   r)   saturn barycenter;m"@ru  r2Qó@rw  3u(@ry  B	@:j(@r|  r>  n}@x{@^ʸ'@r  r  s
             r   'test_saturn_barycenter_geocentric_date0r$    sN   
"#AQg!A!"A!add1g+//223HHv0%%-@))A,K#))+BXBHH(&<*?@CKK+Vi-?@##%H (BXBHH',)>?CKK+Vi-?@ V4BXBHH(&<*?@CKK)6I+=>r   c                 ,   |j                  d      }t        |       | d   j                  |      }| d   }t        ||j                  |      z
  j                  j
                        }t        |t        z  ddt        z         |j                  |      }|j                         \  }}}t        |j                  ddt        z         t        |j                  dd	t        z         |j                         }	|	j                         \  }}}t        |j                  d
dt        z         t        |j                  dd	t        z         |	j                  d      \  }}}t        |j                  ddt        z         t        |j                  dd	t        z         y )Nr=   r)   r  bI"@ru  SOBr3@rw  4ry  Z!3@h4r|  r>  ;&3@$GJ4r  r  s
             r   'test_saturn_barycenter_geocentric_date1r-    sM   
AQg!A!"A!add1g+//223HHv0%%-@))A,K#))+BXBHH(&<*?@CKK,fy.@A##%H (BXBHH',)>?CKK,fy.@A V4BXBHH(&<*?@CKK,fy.@Ar   c                 ,   |j                  d      }t        |       | d   j                  |      }| d   }t        ||j                  |      z
  j                  j
                        }t        |t        z  ddt        z         |j                  |      }|j                         \  }}}t        |j                  ddt        z         t        |j                  dd	t        z         |j                         }	|	j                         \  }}}t        |j                  d
dt        z         t        |j                  dd	t        z         |	j                  d      \  }}}t        |j                  ddt        z         t        |j                  dd	t        z         y )NrA   r)   r  EC5N!@ru  n=kڬ@rw  z2;)@ry  MeC?@o7;)@r|  r>  ,vŬ@:)@r  r  s
             r   'test_saturn_barycenter_geocentric_date2r6    sM   
AQg!A!"A!add1g+//223HHv0%%-@))A,K#))+BXBHH',)>?CKK+Vi-?@##%H (BXBHH',)>?CKK+Vi-?@ V4BXBHH',)>?CKK+Vi-?@r   c                 ,   |j                  d      }t        |       | d   j                  |      }| d   }t        ||j                  |      z
  j                  j
                        }t        |t        z  ddt        z         |j                  |      }|j                         \  }}}t        |j                  ddt        z         t        |j                  dd	t        z         |j                         }	|	j                         \  }}}t        |j                  d
dt        z         t        |j                  dd	t        z         |	j                  d      \  }}}t        |j                  ddt        z         t        |j                  dd	t        z         y )NrD   r)   r  j6$@ru  ](A+@rw  Bry  i0DA+@칱.ҡr|  r>  Jk}G+@#r  r  s
             r   'test_saturn_barycenter_geocentric_date3r?    sM   
AQg!A!"A!add1g+//223HHv155=A))A,K#))+BXBHH(&<*?@CKK+Vi-?@##%H (BXBHH',)>?CKK+Vi-?@ V4BXBHH(&<*?@CKK+Vi-?@r   c                 0   |j                  g d      }t        |       | d   j                  |      }| d   }t        ||j                  |      z
  j                  j
                        }t        |t        z  ddt        z         |j                  |      }|j                         \  }}}t        |j                  ddt        z         t        |j                  dd	t        z         |j                         }	|	j                         \  }}}t        |j                  d
dt        z         t        |j                  dd	t        z         |	j                  d      \  }}}t        |j                  ddt        z         t        |j                  dd	t        z         y )Nr  r)   r  )r  r&  r/  r8  ru  )r  r'  r0  r9  rw  )r  r(  r1  r:  ry  )r   r)  r2  r;  )r!  r*  r3  r<  r|  r>  )r"  r+  r4  r=  )r#  r,  r5  r>  r  r  s
             r   'test_saturn_barycenter_geocentric_date4rA    s]   
EFAQg!A!"A!add1g+//223HHvlnsv{n{|))A,K#))+BXBHHegmp|g|}CKKjlru~l~##%H (BXBHHbdjmydyzCKKjlru~l~ V4BXBHHegmp|g|}CKKhjps|j|}r   c                 ,   |j                  d      }t        |       | d   j                  |      }| d   }t        ||j                  |      z
  j                  j
                        }t        |t        z  ddt        z         |j                  |      }|j                         \  }}}t        |j                  ddt        z         t        |j                  dd	t        z         |j                         }	|	j                         \  }}}t        |j                  d
dt        z         t        |j                  dd	t        z         |	j                  d      \  }}}t        |j                  ddt        z         t        |j                  dd	t        z         y )Nr5   r)   uranus barycenter?&2@ru  w),(@rw  DІ?ry  [n,(@g?r|  r>  vGB(@&Wx(?r  r  s
             r   'test_uranus_barycenter_geocentric_date0rK    sN   
"#AQg!A!"A!add1g+//223HHv0%%-@))A,K#))+BXBHH(&<*?@CKK,fy.@A##%H (BXBHH(&<*?@CKK,fy.@A V4BXBHH(&<*?@CKK,fy.@Ar   c                 ,   |j                  d      }t        |       | d   j                  |      }| d   }t        ||j                  |      z
  j                  j
                        }t        |t        z  ddt        z         |j                  |      }|j                         \  }}}t        |j                  ddt        z         t        |j                  dd	t        z         |j                         }	|	j                         \  }}}t        |j                  d
dt        z         t        |j                  dd	t        z         |	j                  d      \  }}}t        |j                  ddt        z         t        |j                  dd	t        z         y )Nr=   r)   rC  eV2@ru  D0&2@rw  ]o7ry   Z:2@]pWo7r|  r>  &˨2@ywFr7r  r  s
             r   'test_uranus_barycenter_geocentric_date1rT  2  sM   
AQg!A!"A!add1g+//223HHv155=A))A,K#))+BXBHH(&<*?@CKK,fy.@A##%H (BXBHH(&<*?@CKK,fy.@A V4BXBHH',)>?CKK,fy.@Ar   c                 ,   |j                  d      }t        |       | d   j                  |      }| d   }t        ||j                  |      z
  j                  j
                        }t        |t        z  ddt        z         |j                  |      }|j                         \  }}}t        |j                  ddt        z         t        |j                  dd	t        z         |j                         }	|	j                         \  }}}t        |j                  d
dt        z         t        |j                  dd	t        z         |	j                  d      \  }}}t        |j                  ddt        z         t        |j                  dd	t        z         y )NrA   r)   rC  o'4@ru  ˼c*5@rw  '(1ry  OO*5@A|:01r|  r>  euנ<*5@(Q31r  r  s
             r   'test_uranus_barycenter_geocentric_date2r]  J  sM   
AQg!A!"A!add1g+//223HHv155=A))A,K#))+BXBHH(&<*?@CKK,fy.@A##%H (BXBHH(&<*?@CKK,fy.@A V4BXBHH(&<*?@CKK,fy.@Ar   c                 ,   |j                  d      }t        |       | d   j                  |      }| d   }t        ||j                  |      z
  j                  j
                        }t        |t        z  ddt        z         |j                  |      }|j                         \  }}}t        |j                  ddt        z         t        |j                  dd	t        z         |j                         }	|	j                         \  }}}t        |j                  d
dt        z         t        |j                  dd	t        z         |	j                  d      \  }}}t        |j                  ddt        z         t        |j                  dd	t        z         y )NrD   r)   rC  t<3@ru  nwN?rw  篆!@ry  ,uS?t@r|  r>  #?)"&0p@r  r  s
             r   'test_uranus_barycenter_geocentric_date3rf  b  sM   
AQg!A!"A!add1g+//223HHv155=A))A,K#))+BXBHH(&<*?@CKK+Vi-?@##%H (BXBHH(&<*?@CKK*FY,>? V4BXBHH(&<*?@CKK*FY,>?r   c                 0   |j                  g d      }t        |       | d   j                  |      }| d   }t        ||j                  |      z
  j                  j
                        }t        |t        z  ddt        z         |j                  |      }|j                         \  }}}t        |j                  ddt        z         t        |j                  dd	t        z         |j                         }	|	j                         \  }}}t        |j                  d
dt        z         t        |j                  dd	t        z         |	j                  d      \  }}}t        |j                  ddt        z         t        |j                  dd	t        z         y )Nr  r)   rC  )rD  rM  rV  r_  ru  )rE  rN  rW  r`  rw  )rF  rO  rX  ra  ry  )rG  rP  rY  rb  )rH  rQ  rZ  rc  r|  r>  )rI  rR  r[  rd  )rJ  rS  r\  re  r  r  s
             r   'test_uranus_barycenter_geocentric_date4rh  z  sn   
EFAQg!A!"A!add1g+//223HHvnpux}p}~))A,K#))+BXBHHfhnq}h}~CKKlnt  xA  oA  B##%H (BXBHHfhnq}h}~CKKkmsvm  A V4BXBHHegmp|g|}CKKkmsvm  Ar   c                 ,   |j                  d      }t        |       | d   j                  |      }| d   }t        ||j                  |      z
  j                  j
                        }t        |t        z  ddt        z         |j                  |      }|j                         \  }}}t        |j                  ddt        z         t        |j                  dd	t        z         |j                         }	|	j                         \  }}}t        |j                  d
dt        z         t        |j                  dd	t        z         |	j                  d      \  }}}t        |j                  ddt        z         t        |j                  dd	t        z         y )Nr5   r)   neptune barycenterR=@ru  Pw@F/@rw   :1ry  "XF/@5721r|  r>  ퟒ7/@
zs1r  r  s
             r   (test_neptune_barycenter_geocentric_date0rr    sN   
"#AQg!A"#A!add1g+//223HHv0%%-@))A,K#))+BXBHH(&<*?@CKK+Vi-?@##%H (BXBHH',)>?CKK+Vi-?@ V4BXBHH(&<*?@CKK,fy.@Ar   c                 ,   |j                  d      }t        |       | d   j                  |      }| d   }t        ||j                  |      z
  j                  j
                        }t        |t        z  ddt        z         |j                  |      }|j                         \  }}}t        |j                  ddt        z         t        |j                  dd	t        z         |j                         }	|	j                         \  }}}t        |j                  d
dt        z         t        |j                  dd	t        z         |	j                  d      \  }}}t        |j                  ddt        z         t        |j                  dd	t        z         y )Nr=   r)   rj  
+
p}=@ru  s5F	3@rw  "5ry  OX	3@紜5r|  r>   93@,d:5r  r  s
             r   (test_neptune_barycenter_geocentric_date1r{    sM   
AQg!A"#A!add1g+//223HHv155=A))A,K#))+BXBHH',)>?CKK,fy.@A##%H (BXBHH(&<*?@CKK+Vi-?@ V4BXBHH',)>?CKK,fy.@Ar   c                 ,   |j                  d      }t        |       | d   j                  |      }| d   }t        ||j                  |      z
  j                  j
                        }t        |t        z  ddt        z         |j                  |      }|j                         \  }}}t        |j                  ddt        z         t        |j                  dd	t        z         |j                         }	|	j                         \  }}}t        |j                  d
dt        z         t        |j                  dd	t        z         |	j                  d      \  }}}t        |j                  ddt        z         t        |j                  dd	t        z         y )NrA   r)   rj  U<E?@ru  .,3\4@rw  `a63ry  Z'0\4@	C63r|  r>  焸\\4@Ԫ}63r  r  s
             r   (test_neptune_barycenter_geocentric_date2r    sM   
AQg!A"#A!add1g+//223HHv155=A))A,K#))+BXBHH(&<*?@CKK+Vi-?@##%H (BXBHH(&<*?@CKK,fy.@A V4BXBHH',)>?CKK+Vi-?@r   c                 ,   |j                  d      }t        |       | d   j                  |      }| d   }t        ||j                  |      z
  j                  j
                        }t        |t        z  ddt        z         |j                  |      }|j                         \  }}}t        |j                  ddt        z         t        |j                  dd	t        z         |j                         }	|	j                         \  }}}t        |j                  d
dt        z         t        |j                  dd	t        z         |	j                  d      \  }}}t        |j                  ddt        z         t        |j                  dd	t        z         y )NrD   r)   rj  !	(<@ru  1@6@rw  8nb'ry  .;3@6@~E['r|  r>  3C6@z_*&r  r  s
             r   (test_neptune_barycenter_geocentric_date3r    sM   
AQg!A"#A!add1g+//223HHv155=A))A,K#))+BXBHH(&<*?@CKK,fy.@A##%H (BXBHH(&<*?@CKK+Vi-?@ V4BXBHH&(=>CKK,fy.@Ar   c                 0   |j                  g d      }t        |       | d   j                  |      }| d   }t        ||j                  |      z
  j                  j
                        }t        |t        z  ddt        z         |j                  |      }|j                         \  }}}t        |j                  ddt        z         t        |j                  dd	t        z         |j                         }	|	j                         \  }}}t        |j                  d
dt        z         t        |j                  dd	t        z         |	j                  d      \  }}}t        |j                  ddt        z         t        |j                  dd	t        z         y )Nr  r)   rj  )rk  rt  r}  r  ru  )rl  ru  r~  r  rw  )rm  rv  r  r  ry  )rn  rw  r  r  )ro  rx  r  r  r|  r>  )rp  ry  r  r  )rq  rz  r  r  r  r  s
             r   (test_neptune_barycenter_geocentric_date4r    sk   
EFAQg!A"#A!add1g+//223HHvnpux}p}~))A,K#))+BXBHHegmp|g|}CKKkmsvm  A##%H (BXBHHegmp|g|}CKKjlru~l~ V4BXBHHbdjmydyzCKKlnt  xA  oA  Br   c                 ,   |j                  d      }t        |       | d   j                  |      }| d   }t        ||j                  |      z
  j                  j
                        }t        |t        z  ddt        z         |j                  |      }|j                         \  }}}t        |j                  ddt        z         t        |j                  dd	t        z         |j                         }	|	j                         \  }}}t        |j                  d
dt        z         t        |j                  dd	t        z         |	j                  d      \  }}}t        |j                  ddt        z         t        |j                  dd	t        z         y )Nr5   r)   pluto barycenter璚u(@@ru  uL(@rw  ?ܞ0@ry  ~-(@,}H|0@r|  r>  燷ҳ|'@j0@r  r  s
             r   &test_pluto_barycenter_geocentric_date0r  
  sN   
"#AQg!A !A!add1g+//223HHv155=A))A,K#))+BXBHH(&<*?@CKK+Vi-?@##%H (BXBHH',)>?CKK+Vi-?@ V4BXBHH(&<*?@CKK+Vi-?@r   c                 ,   |j                  d      }t        |       | d   j                  |      }| d   }t        ||j                  |      z
  j                  j
                        }t        |t        z  ddt        z         |j                  |      }|j                         \  }}}t        |j                  ddt        z         t        |j                  dd	t        z         |j                         }	|	j                         \  }}}t        |j                  d
dt        z         t        |j                  dd	t        z         |	j                  d      \  }}}t        |j                  ddt        z         t        |j                  dd	t        z         y )Nr=   r)   r  j<@ru  G+n.@rw  Q(NZry  52n.@eA_r|  r>  箯B7j.@ur  r  s
             r   &test_pluto_barycenter_geocentric_date1r  "  sM   
AQg!A !A!add1g+//223HHv155=A))A,K#))+BXBHH(&<*?@CKK,fy.@A##%H (BXBHH(&<*?@CKK,fy.@A V4BXBHH(&<*?@CKK,fy.@Ar   c                 ,   |j                  d      }t        |       | d   j                  |      }| d   }t        ||j                  |      z
  j                  j
                        }t        |t        z  ddt        z         |j                  |      }|j                         \  }}}t        |j                  ddt        z         t        |j                  dd	t        z         |j                         }	|	j                         \  }}}t        |j                  d
dt        z         t        |j                  dd	t        z         |	j                  d      \  }}}t        |j                  ddt        z         t        |j                  dd	t        z         y )NrA   r)   r  OSQ}?@ru   
0@rw  璨H&ry  Wi0@m>	&r|  r>  B^0@v6&r  r  s
             r   &test_pluto_barycenter_geocentric_date2r  :  sM   
AQg!A !A!add1g+//223HHv155=A))A,K#))+BXBHH(&<*?@CKK+Vi-?@##%H (BXBHH(&<*?@CKK,fy.@A V4BXBHH(&<*?@CKK+Vi-?@r   c                 ,   |j                  d      }t        |       | d   j                  |      }| d   }t        ||j                  |      z
  j                  j
                        }t        |t        z  ddt        z         |j                  |      }|j                         \  }}}t        |j                  ddt        z         t        |j                  dd	t        z         |j                         }	|	j                         \  }}}t        |j                  d
dt        z         t        |j                  dd	t        z         |	j                  d      \  }}}t        |j                  ddt        z         t        |j                  dd	t        z         y )NrD   r)   r  SS'?@ru  y}2@rw  }1c\3ry  2)}2@i#>3r|  r>  aW2@"鍊3r  r  s
             r   &test_pluto_barycenter_geocentric_date3r  R  sM   
AQg!A !A!add1g+//223HHv0%%-@))A,K#))+BXBHH(&<*?@CKK+Vi-?@##%H (BXBHH(&<*?@CKK,fy.@A V4BXBHH(&<*?@CKK,fy.@Ar   c                 0   |j                  g d      }t        |       | d   j                  |      }| d   }t        ||j                  |      z
  j                  j
                        }t        |t        z  ddt        z         |j                  |      }|j                         \  }}}t        |j                  ddt        z         t        |j                  dd	t        z         |j                         }	|	j                         \  }}}t        |j                  d
dt        z         t        |j                  dd	t        z         |	j                  d      \  }}}t        |j                  ddt        z         t        |j                  dd	t        z         y )Nr  r)   r  )r  r  r  r  ru  )r  r  r  r  rw  )r  r  r  r  ry  )r  r  r  r  )r  r  r  r  r|  r>  )r  r  r  r  )r  r  r  r  r  r  s
             r   &test_pluto_barycenter_geocentric_date4r  j  sk   
EFAQg!A !A!add1g+//223HHvnpux}p}~))A,K#))+BXBHHfhnq}h}~CKKjlru~l~##%H (BXBHHegmp|g|}CKKlnt  xA  oA  B V4BXBHHfhnq}h}~CKKkmsvm  Ar   c                 ,   |j                  d      }t        |       | d   j                  |      }| d   }t        ||j                  |      z
  j                  j
                        }t        |t        z  ddt        z         |j                  |      }|j                         \  }}}t        |j                  ddt        z         t        |j                  dd	t        z         |j                         }	|	j                         \  }}}t        |j                  d
dt        z         t        |j                  dd	t        z         |	j                  d      \  }}}t        |j                  ddt        z         t        |j                  dd	t        z         y )Nr5   r)   suncYA?ru  2f @rw  H4@ry  K7A3 @Wc4@r|  r>  I+  @玅4@r  r  s
             r   test_sun_geocentric_date0r    sL   
"#AQg!AeA!add1g+//223HHv155=A))A,K#))+BXBHH&(=>CKK+Vi-?@##%H (BXBHH&(=>CKK+Vi-?@ V4BXBHH',)>?CKK*FY,>?r   c                 ,   |j                  d      }t        |       | d   j                  |      }| d   }t        ||j                  |      z
  j                  j
                        }t        |t        z  ddt        z         |j                  |      }|j                         \  }}}t        |j                  ddt        z         t        |j                  dd	t        z         |j                         }	|	j                         \  }}}t        |j                  d
dt        z         t        |j                  dd	t        z         |	j                  d      \  }}}t        |j                  ddt        z         t        |j                  dd	t        z         y )Nr=   r)   r  yRһ0?ru  y5@rw  g3@ry  i4@<3@r|  r>  j~"@AG3@r  r  s
             r   test_sun_geocentric_date1r    sJ   
AQg!AeA!add1g+//223HHv155=A))A,K#))+BXBHH',)>?CKK+Vi-?@##%H (BXBHH',)>?CKK&9*<= V4BXBHH(&<*?@CKK+Vi-?@r   c                 ,   |j                  d      }t        |       | d   j                  |      }| d   }t        ||j                  |      z
  j                  j
                        }t        |t        z  ddt        z         |j                  |      }|j                         \  }}}t        |j                  ddt        z         t        |j                  dd	t        z         |j                         }	|	j                         \  }}}t        |j                  d
dt        z         t        |j                  dd	t        z         |	j                  d      \  }}}t        |j                  ddt        z         t        |j                  dd	t        z         y )NrA   r)   r  v.kw?ru  v2@rw  |7ry  0X2@Jz_7r|  r>  tx2@Z,Q7r  r  s
             r   test_sun_geocentric_date2r    sK   
AQg!AeA!add1g+//223HHv155=A))A,K#))+BXBHH(&<*?@CKK,fy.@A##%H (BXBHH(&<*?@CKK+Vi-?@ V4BXBHH',)>?CKK,fy.@Ar   c                 ,   |j                  d      }t        |       | d   j                  |      }| d   }t        ||j                  |      z
  j                  j
                        }t        |t        z  ddt        z         |j                  |      }|j                         \  }}}t        |j                  ddt        z         t        |j                  dd	t        z         |j                         }	|	j                         \  }}}t        |j                  d
dt        z         t        |j                  dd	t        z         |	j                  d      \  }}}t        |j                  ddt        z         t        |j                  dd	t        z         y )NrD   r)   r  Ǌ.),?ru  CL$@rw  x{%@ry  S$@V%@r|  r>   $@
\Y`%@r  r  s
             r   test_sun_geocentric_date3r    sK   
AQg!AeA!add1g+//223HHv155=A))A,K#))+BXBHH(&<*?@CKK+Vi-?@##%H (BXBHH(&<*?@CKK+Vi-?@ V4BXBHH(&<*?@CKK+Vi-?@r   c                 0   |j                  g d      }t        |       | d   j                  |      }| d   }t        ||j                  |      z
  j                  j
                        }t        |t        z  ddt        z         |j                  |      }|j                         \  }}}t        |j                  ddt        z         t        |j                  dd	t        z         |j                         }	|	j                         \  }}}t        |j                  d
dt        z         t        |j                  dd	t        z         |	j                  d      \  }}}t        |j                  ddt        z         t        |j                  dd	t        z         y )Nr  r)   r  )r  r  r  r  ru  )r  r  r  r  rw  )r  r  r  r  ry  )r  r  r  r  )r  r  r  r  r|  r>  )r  r  r  r  )r  r  r  r  r  r  s
             r   test_sun_geocentric_date4r    s[   
EFAQg!AeA!add1g+//223HHvoqvy~q~))A,K#))+BXBHHceknzez{CKKjlru~l~##%H (BXBHHceknzez{CKKfhnqzhz{ V4BXBHHdflo{f{|CKKikqt}k}~r   c                 ,   |j                  d      }t        |       | d   j                  |      }| d   }t        ||j                  |      z
  j                  j
                        }t        |t        z  ddt        z         |j                  |      }|j                         \  }}}t        |j                  ddt        z         t        |j                  dd	t        z         |j                         }	|	j                         \  }}}t        |j                  d
dt        z         t        |j                  dd	t        z         |	j                  d      \  }}}t        |j                  ddt        z         t        |j                  dd	t        z         y )Nr5   r)   moonbMSe?ru  絏(@rw  EEټ/ry  B(@kDK/r|  r>  y:w|(@/Nr  r  s
             r   test_moon_geocentric_date0r    sL   
"#AQg!AfA!add1g+//223HHv4eemD))A,K#))+BXBHH(&<*?@CKK+Vi-?@##%H (BXBHH(&<*?@CKK+Vi-?@ V4BXBHH(&<*?@CKK+Vi-?@r   c                 ,   |j                  d      }t        |       | d   j                  |      }| d   }t        ||j                  |      z
  j                  j
                        }t        |t        z  ddt        z         |j                  |      }|j                         \  }}}t        |j                  ddt        z         t        |j                  dd	t        z         |j                         }	|	j                         \  }}}t        |j                  d
dt        z         t        |j                  dd	t        z         |	j                  d      \  }}}t        |j                  ddt        z         t        |j                  dd	t        z         y )Nr=   r)   r  h9Td?ru  Z	l+7@rw  +ev?ry  V!7@h:Ȝ?r|  r>  ִE7@<?r  r  s
             r   test_moon_geocentric_date1r    sK   
AQg!AfA!add1g+//223HHv4eemD))A,K#))+BXBHH(&<*?@CKK+Vi-?@##%H (BXBHH(&<*?@CKK*FY,>? V4BXBHH',)>?CKK+Vi-?@r   c                 ,   |j                  d      }t        |       | d   j                  |      }| d   }t        ||j                  |      z
  j                  j
                        }t        |t        z  ddt        z         |j                  |      }|j                         \  }}}t        |j                  ddt        z         t        |j                  dd	t        z         |j                         }	|	j                         \  }}}t        |j                  d
dt        z         t        |j                  dd	t        z         |	j                  d      \  }}}t        |j                  ddt        z         t        |j                  dd	t        z         y )NrA   r)   r  V^	f?ru  -t-@rw  4V %ry  xnܨ-@㬋%r|  r>  >ӽ-@vgG%r  r  s
             r   test_moon_geocentric_date2r  *  sK   
AQg!AfA!add1g+//223HHv3UU]C))A,K#))+BXBHH(&<*?@CKK,fy.@A##%H (BXBHH(&<*?@CKK+Vi-?@ V4BXBHH(&<*?@CKK,fy.@Ar   c                 ,   |j                  d      }t        |       | d   j                  |      }| d   }t        ||j                  |      z
  j                  j
                        }t        |t        z  ddt        z         |j                  |      }|j                         \  }}}t        |j                  ddt        z         t        |j                  dd	t        z         |j                         }	|	j                         \  }}}t        |j                  d
dt        z         t        |j                  dd	t        z         |	j                  d      \  }}}t        |j                  ddt        z         t        |j                  dd	t        z         y )NrD   r)   r  7+)Dd?ru  kd0@rw  珡4ry  d0@sĴz4r|  r>  碛zag0@4r  r  s
             r   test_moon_geocentric_date3r  B  sK   
AQg!AfA!add1g+//223HHv4eemD))A,K#))+BXBHH',)>?CKK+Vi-?@##%H (BXBHH',)>?CKK,fy.@A V4BXBHH',)>?CKK+Vi-?@r   c                 0   |j                  g d      }t        |       | d   j                  |      }| d   }t        ||j                  |      z
  j                  j
                        }t        |t        z  ddt        z         |j                  |      }|j                         \  }}}t        |j                  ddt        z         t        |j                  dd	t        z         |j                         }	|	j                         \  }}}t        |j                  d
dt        z         t        |j                  dd	t        z         |	j                  d      \  }}}t        |j                  ddt        z         t        |j                  dd	t        z         y )Nr  r)   r  )r  r  r  r  ru  )r  r  r  r  rw  )r  r  r  r  ry  )r  r  r  r  )r  r  r  r  r|  r>  )r  r  r  r   )r  r  r  r  r  r  s
             r   test_moon_geocentric_date4r  Z  sj   
EFAQg!AfA!add1g+//223HHvz  }B  EJ  }J  K))A,K#))+BXBHHegmp|g|}CKKjlru~l~##%H (BXBHHegmp|g|}CKKikqt}k}~ V4BXBHHdflo{f{|CKKjlru~l~r   c                    | j                  t        j                         j                  d            }t	        j
                  dddddd      }|j                  |      }|j                         \  }}}t        |j                  d	d
t        z         t        |j                  dd
t        z         |j                         }|j                         \  }}}t        |j                  dd
t        z         t        |j                  dd
t        z         |j                  d      \  }}}t        |j                  dd
t        z         t        |j                  dd
t        z         y )Nr5   r   r   r   r   r   r   r   a:@rG   PV@ .@璂PV@r|  r>  狂O @,GV@rC  r   r   r   r   r   r  rE  r"   rF  r  rG  r:   r  r)   rp  r   r  rK  rL  rM  r  s           r   test_polaris_geocentric_date0r  r  s    !''(9:;A<<,(-%)5BD ))D/K#))+BXBHH('L*@ACKK)7Y+>?##%H (BXBHH&,(>?CKK*Gi,?@ V4BXBHH('L*@ACKK*Gi,?@r   c                    | j                  t        j                         j                  d            }t	        j
                  dddddd      }|j                  |      }|j                         \  }}}t        |j                  d	d
t        z         t        |j                  dd
t        z         |j                         }|j                         \  }}}t        |j                  dd
t        z         t        |j                  dd
t        z         |j                  d      \  }}}t        |j                  dd
t        z         t        |j                  dd
t        z         y )Nr=   r   r   r   r   r   r   r   XD<@rG   \
PV@T(@3PV@r|  r>  #tҩ@J"NV@r  r  s           r   test_polaris_geocentric_date1r       !''	23A<<,(-%)5BD ))D/K#))+BXBHH('L*@ACKK*Gi,?@##%H (BXBHH'<)?@CKK*Gi,?@ V4BXBHH('L*@ACKK*Gi,?@r   c                    | j                  t        j                         j                  d            }t	        j
                  dddddd      }|j                  |      }|j                         \  }}}t        |j                  d	d
t        z         t        |j                  dd
t        z         |j                         }|j                         \  }}}t        |j                  dd
t        z         t        |j                  dd
t        z         |j                  d      \  }}}t        |j                  dd
t        z         t        |j                  dd
t        z         y )NrA   r   r   r   r   r   r   r   %	>@rG   S.PV@fŨh[@9:QV@r|  r>  }W4^@rQV@r  r  s           r   test_polaris_geocentric_date2r    r  r   c                    | j                  t        j                         j                  d            }t	        j
                  dddddd      }|j                  |      }|j                         \  }}}t        |j                  d	d
t        z         t        |j                  dd
t        z         |j                         }|j                         \  }}}t        |j                  dd
t        z         t        |j                  dd
t        z         |j                  d      \  }}}t        |j                  dd
t        z         t        |j                  dd
t        z         y )NrD   r   r   r   r   r   r   r   %M?@rG   rPV@xa7U@IPV@r|  r>  sᵨs@TV@r  r  s           r   test_polaris_geocentric_date3r'    s   !''	23A<<,(-%)5BD ))D/K#))+BXBHH'<)?@CKK*Gi,?@##%H (BXBHH'<)?@CKK*Gi,?@ V4BXBHH('L*@ACKK)7Y+>?r   c                    | j                  t        j                         j                  g d            }t	        j
                  dddddd      }|j                  |      }|j                         \  }}}t        |j                  d	d
t        z         t        |j                  dd
t        z         |j                         }|j                         \  }}}t        |j                  dd
t        z         t        |j                  dd
t        z         |j                  d      \  }}}t        |j                  dd
t        z         t        |j                  dd
t        z         y )Nr  r   r   r   r   r   r   r   )r  r  r  r!  rG   )r  r  r  r"  )r  r  r  r#  )r	  r  r  r$  r|  r>  )r
  r  r  r%  )r  r  r  r&  r  r  s           r   test_polaris_geocentric_date4r)    s,   !''(\]^A<<,(-%)5BD ))D/K#))+BXBHHegnq}g}~CKKdfmpyfyz##%H (BXBHHacjmycyzCKKegnqzgz{ V4BXBHHfhor~h~CKKdfmpyfyzr   c                    t        j                  d      j                  d      }| d   }|t        ddd      z   }|j	                  |      j                  | d	         j                         }|j                         \  }}}t        |j                  d
dt        z         t        |j                  ddt        z         |j                  d      \  }}}t        |j                  ddt        z         t        |j                  ddt        z         |j                         \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         |j                  d      \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         |j                  dd      \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         y )Nr+   r   r5   r)   r6  r7  r8  r9  rs  ȡ@Mb@?^$/U6@r|  r>  k7e@տPj6@W)G@Wwbp@standard+G@r  r  s4$+G@r   r   r   r   rC  r  r  rE  r"   rF  r  rG  r:   altaz
r&   r1   r)   rI  r  rK  rL  rM  r  azs
             r   test_mercury_topocentric_date0r9    s   s#))*;<A'NE5ktDDDwwqz!!%	"23<<>H (BXBHH(&<*?@CKK*FY,>? V4BXBHH',)>?CKK+Vi-?@ (CXCKK)6I+=>BJJ*FY,>? z2CXCKK*FY,>?BJJ*FY,>? tV4CXCKK*FY,>?BJJ*FY,>?r   c                    t        j                  d      j                  d      }| d   }|t        ddd      z   }|j	                  |      j                  | d	         j                         }|j                         \  }}}t        |j                  d
dt        z         t        |j                  ddt        z         |j                  d      \  }}}t        |j                  ddt        z         t        |j                  ddt        z         |j                         \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         |j                  d      \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         |j                  dd      \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         y )Nr+   r   r=   r)   r6  r7  r8  r9  rs  f,z@r,  K}&@r|  r>  st4+@ϊi&@|Y5W1fr@r2  r  r  r5  r7  s
             r   test_mercury_topocentric_date1rA    s   s#)))4A'NE5ktDDDwwqz!!%	"23<<>H (BXBHH',)>?CKK*FY,>? V4BXBHH',)>?CKK+Vi-?@ (CXCKK,fy.@ABJJ)6I+=> z2CXCKK,fy.@ABJJ)6I+=> tV4CXCKK,fy.@ABJJ)6I+=>r   c                    t        j                  d      j                  d      }| d   }|t        ddd      z   }|j	                  |      j                  | d	         j                         }|j                         \  }}}t        |j                  d
dt        z         t        |j                  ddt        z         |j                  d      \  }}}t        |j                  ddt        z         t        |j                  ddt        z         |j                         \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         |j                  d      \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         |j                  dd      \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         y )Nr+   r   rA   r)   r6  r7  r8  r9  rs  1{#2@r,  @|׼k8r|  r>  h#2@Z,Rk8 B?VU~^@r2   V'5?r  r   U?r5  r7  s
             r   test_mercury_topocentric_date2rK  
  s   s#)))4A'NE5ktDDDwwqz!!%	"23<<>H (BXBHH(&<*?@CKK,fy.@A V4BXBHH(&<*?@CKK,fy.@A (CXCKK-v	/ABBJJ*FY,>? z2CXCKK,fy.@ABJJ*FY,>? tV4CXCKK+Vi-?@BJJ*FY,>?r   c                    t        j                  d      j                  d      }| d   }|t        ddd      z   }|j	                  |      j                  | d	         j                         }|j                         \  }}}t        |j                  d
dt        z         t        |j                  ddt        z         |j                  d      \  }}}t        |j                  ddt        z         t        |j                  ddt        z         |j                         \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         |j                  d      \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         |j                  dd      \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         y )Nr+   r   rD   r)   r6  r7  r8  r9  rs  2G"@r,  nì0@r|  r>  	ZHk"@vW0@D[4;"vEr@r2  r  r  r5  r7  s
             r   test_mercury_topocentric_date3rS  $  s   s#)))4A'NE5ktDDDwwqz!!%	"23<<>H (BXBHH&(=>CKK*FY,>? V4BXBHH',)>?CKK+Vi-?@ (CXCKK+Vi-?@BJJ)6I+=> z2CXCKK+Vi-?@BJJ)6I+=> tV4CXCKK+Vi-?@BJJ)6I+=>r   c                 
   t        j                  d      j                  g d      }| d   }|t        ddd      z   }|j	                  |      j                  | d	         j                         }|j                         \  }}}t        |j                  d
dt        z         t        |j                  ddt        z         |j                  d      \  }}}t        |j                  ddt        z         t        |j                  ddt        z         |j                         \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         |j                  d      \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         |j                  dd      \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         y )Nr+   r   r  r)   r6  r7  r8  r9  rs  )r+  r;  rC  rM  r,  )r-  r<  rD  rN  r|  r>  )r.  r=  rE  rO  )r/  r>  rF  rP  )r0  r?  rG  rQ  )r1  r@  rH  rR  r2  )r3  r?  rI  rQ  r  r  )r4  r?  rJ  rQ  r5  r7  s
             r   test_mercury_topocentric_date4rU  >  s   s#))*^_A'NE5ktDDDwwqz!!%	"23<<>H (BXBHHceknzez{CKKgior{i{| V4BXBHHceknzez{CKKjlru~l~ (CXCKKjlru~l~BJJfhnqzhz{ z2CXCKKjlru~l~BJJfhnqzhz{ tV4CXCKKikqt}k}~BJJfhnqzhz{r   c                    t        j                  d      j                  d      }| d   }|t        ddd      z   }|j	                  |      j                  | d	         j                         }|j                         \  }}}t        |j                  d
dt        z         t        |j                  ddt        z         |j                  d      \  }}}t        |j                  ddt        z         t        |j                  ddt        z         |j                         \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         |j                  d      \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         |j                  dd      \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         y )Nr+   r   r5   r)   r6  r7  r8  r9  r  B@r,  *PNPk54@r|  r>  U@V(H*4@8"N&@Yq@r2  瀈	v&@r  r   <Y1w&@r5  r7  s
             r   test_venus_topocentric_date0r_  X  s   s#))*;<A'NE5ktDDDwwqz!!%.1::<H (BXBHH(&<*?@CKK*FY,>? V4BXBHH',)>?CKK+Vi-?@ (CXCKK+Vi-?@BJJ)6I+=> z2CXCKK*FY,>?BJJ)6I+=> tV4CXCKK+Vi-?@BJJ)6I+=>r   c                    t        j                  d      j                  d      }| d   }|t        ddd      z   }|j	                  |      j                  | d	         j                         }|j                         \  }}}t        |j                  d
dt        z         t        |j                  ddt        z         |j                  d      \  }}}t        |j                  ddt        z         t        |j                  ddt        z         |j                         \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         |j                  d      \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         |j                  dd      \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         y )Nr+   r   r=   r)   r6  r7  r8  r9  r  m/]?r,  fJ&M@r|  r>  ft?њpp@vMDAs?C1as@r2  r  r  r5  r7  s
             r   test_venus_topocentric_date1rg  r  s   s#)))4A'NE5ktDDDwwqz!!%.1::<H (BXBHH(&<*?@CKK*FY,>? V4BXBHH(&<*?@CKK*FY,>? (CXCKK,fy.@ABJJ*FY,>? z2CXCKK,fy.@ABJJ*FY,>? tV4CXCKK,fy.@ABJJ*FY,>?r   c                    t        j                  d      j                  d      }| d   }|t        ddd      z   }|j	                  |      j                  | d	         j                         }|j                         \  }}}t        |j                  d
dt        z         t        |j                  ddt        z         |j                  d      \  }}}t        |j                  ddt        z         t        |j                  ddt        z         |j                         \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         |j                  d      \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         |j                  dd      \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         y )Nr+   r   rA   r)   r6  r7  r8  r9  r  :,4y/@r,  Ʀks2r|  r>  FX/@c+AMs2Ĩ:7@ja@r2  GKD7@r  r  BrdD7@r5  r7  s
             r   test_venus_topocentric_date2rq    s   s#)))4A'NE5ktDDDwwqz!!%.1::<H (BXBHH',)>?CKK+Vi-?@ V4BXBHH',)>?CKK,fy.@A (CXCKK+Vi-?@BJJ)6I+=> z2CXCKK+Vi-?@BJJ)6I+=> tV4CXCKK+Vi-?@BJJ)6I+=>r   c                    t        j                  d      j                  d      }| d   }|t        ddd      z   }|j	                  |      j                  | d	         j                         }|j                         \  }}}t        |j                  d
dt        z         t        |j                  ddt        z         |j                  d      \  }}}t        |j                  ddt        z         t        |j                  ddt        z         |j                         \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         |j                  d      \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         |j                  dd      \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         y )Nr+   r   rD   r)   r6  r7  r8  r9  r  #ll@r,  U>K3@r|  r>  bs@@a3@爗Ĩ(\8'D8?zt@r2  r  r  r5  r7  s
             r   test_venus_topocentric_date3ry    s   s#)))4A'NE5ktDDDwwqz!!%.1::<H (BXBHH',)>?CKK*FY,>? V4BXBHH',)>?CKK+Vi-?@ (CXCKK,fy.@ABJJ(&9*<= z2CXCKK,fy.@ABJJ(&9*<= tV4CXCKK,fy.@ABJJ(&9*<=r   c                 
   t        j                  d      j                  g d      }| d   }|t        ddd      z   }|j	                  |      j                  | d	         j                         }|j                         \  }}}t        |j                  d
dt        z         t        |j                  ddt        z         |j                  d      \  }}}t        |j                  ddt        z         t        |j                  ddt        z         |j                         \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         |j                  d      \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         |j                  dd      \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         y )Nr+   r   r  r)   r6  r7  r8  r9  r  )rW  ra  ri  rs  r,  )rX  rb  rj  rt  r|  r>  )rY  rc  rk  ru  )rZ  rd  rl  rv  )r[  re  rm  rw  )r\  rf  rn  rx  r2  )r]  re  ro  rw  r  r  )r^  re  rp  rw  r5  r7  s
             r   test_venus_topocentric_date4r{    s   s#))*^_A'NE5ktDDDwwqz!!%.1::<H (BXBHHdflo{f{|CKKfhnqzhz{ V4BXBHHceknzez{CKKikqt}k}~ (CXCKKkmsvm  ABJJdfloxfxy z2CXCKKjlru~l~BJJdfloxfxy tV4CXCKKkmsvm  ABJJdfloxfxyr   c                    t        j                  d      j                  d      }| d   }|t        ddd      z   }|j	                  |      j                  | d	         j                         }|j                         \  }}}t        |j                  d
dt        z         t        |j                  ddt        z         |j                  d      \  }}}t        |j                  ddt        z         t        |j                  ddt        z         |j                         \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         |j                  d      \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         |j                  dd      \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         y )Nr+   r   r5   r)   r6  r7  r8  r9  r  's0@r,  J!8r|  r>  $/@R8瀅RZR]@r2  r  r  r5  r7  s
             r   test_mars_topocentric_date0r    s   s#))*;<A'NE5ktDDDwwqz!!%-099;H (BXBHH(&<*?@CKK,fy.@A V4BXBHH(&<*?@CKK,fy.@A (CXCKK+Vi-?@BJJ*FY,>? z2CXCKK+Vi-?@BJJ*FY,>? tV4CXCKK+Vi-?@BJJ*FY,>?r   c                    t        j                  d      j                  d      }| d   }|t        ddd      z   }|j	                  |      j                  | d	         j                         }|j                         \  }}}t        |j                  d
dt        z         t        |j                  ddt        z         |j                  d      \  }}}t        |j                  ddt        z         t        |j                  ddt        z         |j                         \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         |j                  d      \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         |j                  dd      \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         y )Nr+   r   r=   r)   r6  r7  r8  r9  r  juyv|7@r,  c!O)r|  r>  ^UCm7@	<瘒~K*/ u@r2  r  r  r5  r7  s
             r   test_mars_topocentric_date1r    s   s#)))4A'NE5ktDDDwwqz!!%-099;H (BXBHH',)>?CKK+Vi-?@ V4BXBHH',)>?CKK+Vi-?@ (CXCKK*FY,>?BJJ)6I+=> z2CXCKK*FY,>?BJJ)6I+=> tV4CXCKK*FY,>?BJJ)6I+=>r   c                    t        j                  d      j                  d      }| d   }|t        ddd      z   }|j	                  |      j                  | d	         j                         }|j                         \  }}}t        |j                  d
dt        z         t        |j                  ddt        z         |j                  d      \  }}}t        |j                  ddt        z         t        |j                  ddt        z         |j                         \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         |j                  d      \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         |j                  dd      \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         y )Nr+   r   rA   r)   r6  r7  r8  r9  r  玻6@r,  盒]*r|  r>  C6@#]*窯FsBkrS@r2  r  r  r5  r7  s
             r   test_mars_topocentric_date2r    s   s#)))4A'NE5ktDDDwwqz!!%-099;H (BXBHH(&<*?@CKK,fy.@A V4BXBHH',)>?CKK,fy.@A (CXCKK+Vi-?@BJJ)6I+=> z2CXCKK+Vi-?@BJJ)6I+=> tV4CXCKK+Vi-?@BJJ)6I+=>r   c                    t        j                  d      j                  d      }| d   }|t        ddd      z   }|j	                  |      j                  | d	         j                         }|j                         \  }}}t        |j                  d
dt        z         t        |j                  ddt        z         |j                  d      \  }}}t        |j                  ddt        z         t        |j                  ddt        z         |j                         \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         |j                  d      \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         |j                  dd      \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         y )Nr+   r   rD   r)   r6  r7  r8  r9  r  a+@r,  q).>(r|  r>  *pI+@Ɲ_()oD6@ kl@r2  簽:"6@r  r  \A"6@r5  r7  s
             r   test_mars_topocentric_date3r  (  s   s#)))4A'NE5ktDDDwwqz!!%-099;H (BXBHH&(=>CKK,fy.@A V4BXBHH(&<*?@CKK,fy.@A (CXCKK+Vi-?@BJJ)6I+=> z2CXCKK+Vi-?@BJJ)6I+=> tV4CXCKK+Vi-?@BJJ)6I+=>r   c                 
   t        j                  d      j                  g d      }| d   }|t        ddd      z   }|j	                  |      j                  | d	         j                         }|j                         \  }}}t        |j                  d
dt        z         t        |j                  ddt        z         |j                  d      \  }}}t        |j                  ddt        z         t        |j                  ddt        z         |j                         \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         |j                  d      \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         |j                  dd      \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         y )Nr+   r   r  r)   r6  r7  r8  r9  r  )r}  r  r  r  r,  )r~  r  r  r  r|  r>  )r  r  r  r  )r  r  r  r  )r  r  r  r  )r  r  r  r  r2  )r  r  r  r  r  r  )r  r  r  r  r5  r7  s
             r   test_mars_topocentric_date4r  B  s   s#))*^_A'NE5ktDDDwwqz!!%-099;H (BXBHHceknzez{CKKlnt  xA  oA  B V4BXBHHdflo{f{|CKKlnt  xA  oA  B (CXCKKhjps|j|}BJJegmpygyz z2CXCKKhjps|j|}BJJegmpygyz tV4CXCKKhjps|j|}BJJegmpygyzr   c                    t        j                  d      j                  d      }| d   }|t        ddd      z   }|j	                  |      j                  | d	         j                         }|j                         \  }}}t        |j                  d
dt        z         t        |j                  ddt        z         |j                  d      \  }}}t        |j                  ddt        z         t        |j                  ddt        z         |j                         \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         |j                  d      \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         |j                  dd      \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         y )Nr+   r   r5   r)   r6  r7  r8  r9  r  y85(@r,  y$?r|  r>  4'(@ i'H?.H@Dc@r2  犐4;յH@r  r  xٵH@r5  r7  s
             r   )test_jupiter_barycenter_topocentric_date0r  \  s   s#))*;<A'NE5ktDDDwwqz!!%(<"=>GGIH (BXBHH(&<*?@CKK+Vi-?@ V4BXBHH(&<*?@CKK*FY,>? (CXCKK*FY,>?BJJ*FY,>? z2CXCKK*FY,>?BJJ*FY,>? tV4CXCKK+Vi-?@BJJ*FY,>?r   c                    t        j                  d      j                  d      }| d   }|t        ddd      z   }|j	                  |      j                  | d	         j                         }|j                         \  }}}t        |j                  d
dt        z         t        |j                  ddt        z         |j                  d      \  }}}t        |j                  ddt        z         t        |j                  ddt        z         |j                         \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         |j                  d      \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         |j                  dd      \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         y )Nr+   r   r=   r)   r6  r7  r8  r9  r  	9b1@r,  +7@r|  r>  /~@	g.7@B C@}5p@r2  'LTC@r  r  Ko[C@r5  r7  s
             r   )test_jupiter_barycenter_topocentric_date1r  v  s   s#)))4A'NE5ktDDDwwqz!!%(<"=>GGIH (BXBHH',)>?CKK*FY,>? V4BXBHH',)>?CKK*FY,>? (CXCKK*FY,>?BJJ*FY,>? z2CXCKK*FY,>?BJJ*FY,>? tV4CXCKK*FY,>?BJJ*FY,>?r   c                    t        j                  d      j                  d      }| d   }|t        ddd      z   }|j	                  |      j                  | d	         j                         }|j                         \  }}}t        |j                  d
dt        z         t        |j                  ddt        z         |j                  d      \  }}}t        |j                  ddt        z         t        |j                  ddt        z         |j                         \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         |j                  d      \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         |j                  dd      \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         y )Nr+   r   rA   r)   r6  r7  r8  r9  r  JTQlv?r,  1!@r|  r>  sbJ}u?uf\0!@Վ=Ez=:uv@r2  r  r  r5  r7  s
             r   )test_jupiter_barycenter_topocentric_date2r    s   s#)))4A'NE5ktDDDwwqz!!%(<"=>GGIH (BXBHH(&<*?@CKK*FY,>? V4BXBHH(&<*?@CKK*FY,>? (CXCKK,fy.@ABJJ)6I+=> z2CXCKK,fy.@ABJJ)6I+=> tV4CXCKK,fy.@ABJJ)6I+=>r   c                    t        j                  d      j                  d      }| d   }|t        ddd      z   }|j	                  |      j                  | d	         j                         }|j                         \  }}}t        |j                  d
dt        z         t        |j                  ddt        z         |j                  d      \  }}}t        |j                  ddt        z         t        |j                  ddt        z         |j                         \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         |j                  d      \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         |j                  dd      \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         y )Nr+   r   rD   r)   r6  r7  r8  r9  r  竨{J@r,  9aG5@r|  r>  i-W@{Z5@XI=/֥HO@r2  r  r  r5  r7  s
             r   )test_jupiter_barycenter_topocentric_date3r    s   s#)))4A'NE5ktDDDwwqz!!%(<"=>GGIH (BXBHH&(=>CKK+Vi-?@ V4BXBHH',)>?CKK+Vi-?@ (CXCKK,fy.@ABJJ)6I+=> z2CXCKK,fy.@ABJJ)6I+=> tV4CXCKK,fy.@ABJJ)6I+=>r   c                 
   t        j                  d      j                  g d      }| d   }|t        ddd      z   }|j	                  |      j                  | d	         j                         }|j                         \  }}}t        |j                  d
dt        z         t        |j                  ddt        z         |j                  d      \  }}}t        |j                  ddt        z         t        |j                  ddt        z         |j                         \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         |j                  d      \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         |j                  dd      \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         y )Nr+   r   r  r)   r6  r7  r8  r9  r  )r  r  r  r  r,  )r  r  r  r  r|  r>  )r  r  r  r  )r  r  r  r  )r  r  r  r  )r  r  r  r  r2  )r  r  r  r  r  r  )r  r  r  r  r5  r7  s
             r   )test_jupiter_barycenter_topocentric_date4r    s   s#))*^_A'NE5ktDDDwwqz!!%(<"=>GGIH (BXBHHceknzez{CKKgior{i{| V4BXBHHdflo{f{|CKKfhnqzhz{ (CXCKKikqt}k}~BJJfhnqzhz{ z2CXCKKikqt}k}~BJJfhnqzhz{ tV4CXCKKjlru~l~BJJfhnqzhz{r   c                    t        j                  d      j                  d      }| d   }|t        ddd      z   }|j	                  |      j                  | d	         j                         }|j                         \  }}}t        |j                  d
dt        z         t        |j                  ddt        z         |j                  d      \  }}}t        |j                  ddt        z         t        |j                  ddt        z         |j                         \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         |j                  d      \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         |j                  dd      \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         y )Nr+   r   r5   r)   r6  r7  r8  r9  r  Aژ@r,  )w=S(@r|  r>   q{@%'@-٥4a
Q s@r2  r  r  r5  r7  s
             r   (test_saturn_barycenter_topocentric_date0r    s   s#))*;<A'NE5ktDDDwwqz!!%(;"<=FFHH (BXBHH(&<*?@CKK+Vi-?@ V4BXBHH(&<*?@CKK+Vi-?@ (CXCKK,fy.@ABJJ*FY,>? z2CXCKK,fy.@ABJJ*FY,>? tV4CXCKK,fy.@ABJJ*FY,>?r   c                    t        j                  d      j                  d      }| d   }|t        ddd      z   }|j	                  |      j                  | d	         j                         }|j                         \  }}}t        |j                  d
dt        z         t        |j                  ddt        z         |j                  d      \  }}}t        |j                  ddt        z         t        |j                  ddt        z         |j                         \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         |j                  d      \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         |j                  dd      \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         y )Nr+   r   r=   r)   r6  r7  r8  r9  r  .Y3@r,  `_4r|  r>  碠'3@$O4dxwHE8S@r2  r  r  r5  r7  s
             r   (test_saturn_barycenter_topocentric_date1r    s   s#)))4A'NE5ktDDDwwqz!!%(;"<=FFHH (BXBHH(&<*?@CKK,fy.@A V4BXBHH(&<*?@CKK,fy.@A (CXCKK+Vi-?@BJJ(&9*<= z2CXCKK+Vi-?@BJJ(&9*<= tV4CXCKK+Vi-?@BJJ(&9*<=r   c                    t        j                  d      j                  d      }| d   }|t        ddd      z   }|j	                  |      j                  | d	         j                         }|j                         \  }}}t        |j                  d
dt        z         t        |j                  ddt        z         |j                  d      \  }}}t        |j                  ddt        z         t        |j                  ddt        z         |j                         \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         |j                  d      \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         |j                  dd      \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         y )Nr+   r   rA   r)   r6  r7  r8  r9  r  l:@r,  r;)@r|  r>  M2@}>:)@Z>@BWSu@r2  r  r  r5  r7  s
             r   (test_saturn_barycenter_topocentric_date2r  	  s   s#)))4A'NE5ktDDDwwqz!!%(;"<=FFHH (BXBHH(&<*?@CKK+Vi-?@ V4BXBHH',)>?CKK+Vi-?@ (CXCKK,fy.@ABJJ*FY,>? z2CXCKK,fy.@ABJJ*FY,>? tV4CXCKK,fy.@ABJJ*FY,>?r   c                    t        j                  d      j                  d      }| d   }|t        ddd      z   }|j	                  |      j                  | d	         j                         }|j                         \  }}}t        |j                  d
dt        z         t        |j                  ddt        z         |j                  d      \  }}}t        |j                  ddt        z         t        |j                  ddt        z         |j                         \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         |j                  d      \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         |j                  dd      \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         y )Nr+   r   rD   r)   r6  r7  r8  r9  r  a:~\A+@r,  b{r|  r>  #ă|G+@DH|q|}6@Ak3m@r2  PJ:M7@r  r  [Xg7@r5  r7  s
             r   (test_saturn_barycenter_topocentric_date3r  ,	  s   s#)))4A'NE5ktDDDwwqz!!%(;"<=FFHH (BXBHH(&<*?@CKK+Vi-?@ V4BXBHH(&<*?@CKK+Vi-?@ (CXCKK*FY,>?BJJ*FY,>? z2CXCKK+Vi-?@BJJ*FY,>? tV4CXCKK+Vi-?@BJJ*FY,>?r   c                 
   t        j                  d      j                  g d      }| d   }|t        ddd      z   }|j	                  |      j                  | d	         j                         }|j                         \  }}}t        |j                  d
dt        z         t        |j                  ddt        z         |j                  d      \  }}}t        |j                  ddt        z         t        |j                  ddt        z         |j                         \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         |j                  d      \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         |j                  dd      \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         y )Nr+   r   r  r)   r6  r7  r8  r9  r  )r  r  r  r  r,  )r  r  r  r  r|  r>  )r  r  r  r  )r  r  r  r  )r  r  r  r  )r  r  r  r  r2  )r  r  r  r  r  r  )r  r  r  r  r5  r7  s
             r   (test_saturn_barycenter_topocentric_date4r  F	  s   s#))*^_A'NE5ktDDDwwqz!!%(;"<=FFHH (BXBHHfhnq}h}~CKKjlru~l~ V4BXBHHegmp|g|}CKKjlru~l~ (CXCKKjlru~l~BJJfhnqzhz{ z2CXCKKkmsvm  ABJJfhnqzhz{ tV4CXCKKkmsvm  ABJJfhnqzhz{r   c                    t        j                  d      j                  d      }| d   }|t        ddd      z   }|j	                  |      j                  | d	         j                         }|j                         \  }}}t        |j                  d
dt        z         t        |j                  ddt        z         |j                  d      \  }}}t        |j                  ddt        z         t        |j                  ddt        z         |j                         \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         |j                  d      \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         |j                  dd      \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         y )Nr+   r   r5   r)   r6  r7  r8  r9  rC  @q,(@r,  c?r|  r>  C(@UFV'?缌U0H@m/Zc@r2  0p6H@r  r  K]H@r5  r7  s
             r   (test_uranus_barycenter_topocentric_date0r  `	  s   s#))*;<A'NE5ktDDDwwqz!!%(;"<=FFHH (BXBHH(&<*?@CKK,fy.@A V4BXBHH(&<*?@CKK,fy.@A (CXCKK*FY,>?BJJ*FY,>? z2CXCKK+Vi-?@BJJ*FY,>? tV4CXCKK*FY,>?BJJ*FY,>?r   c                    t        j                  d      j                  d      }| d   }|t        ddd      z   }|j	                  |      j                  | d	         j                         }|j                         \  }}}t        |j                  d
dt        z         t        |j                  ddt        z         |j                  d      \  }}}t        |j                  ddt        z         t        |j                  ddt        z         |j                         \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         |j                  d      \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         |j                  dd      \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         y )Nr+   r   r=   r)   r6  r7  r8  r9  rC  :2@r,  nao7r|  r>  -&̨2@
ZQr7RBޭV@r2  r  r  r5  r7  s
             r   (test_uranus_barycenter_topocentric_date1r  z	  s   s#)))4A'NE5ktDDDwwqz!!%(;"<=FFHH (BXBHH(&<*?@CKK+Vi-?@ V4BXBHH',)>?CKK,fy.@A (CXCKK*FY,>?BJJ)6I+=> z2CXCKK*FY,>?BJJ)6I+=> tV4CXCKK*FY,>?BJJ)6I+=>r   c                    t        j                  d      j                  d      }| d   }|t        ddd      z   }|j	                  |      j                  | d	         j                         }|j                         \  }}}t        |j                  d
dt        z         t        |j                  ddt        z         |j                  d      \  }}}t        |j                  ddt        z         t        |j                  ddt        z         |j                         \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         |j                  d      \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         |j                  dd      \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         y )Nr+   r   rA   r)   r6  r7  r8  r9  rC  珡ZO*5@r,  t21r|  r>  H<*5@dj61Pj,=G_6V@r2  r  r  r5  r7  s
             r   (test_uranus_barycenter_topocentric_date2r  	  s   s#)))4A'NE5ktDDDwwqz!!%(;"<=FFHH (BXBHH',)>?CKK,fy.@A V4BXBHHov'<=CKK,fy.@A (CXCKK,fy.@ABJJ)6I+=> z2CXCKK,fy.@ABJJ)6I+=> tV4CXCKK,fy.@ABJJ)6I+=>r   c                    t        j                  d      j                  d      }| d   }|t        ddd      z   }|j	                  |      j                  | d	         j                         }|j                         \  }}}t        |j                  d
dt        z         t        |j                  ddt        z         |j                  d      \  }}}t        |j                  ddt        z         t        |j                  ddt        z         |j                         \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         |j                  d      \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         |j                  dd      \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         y )Nr+   r   rD   r)   r6  r7  r8  r9  rC  )S?r,  	[L'@r|  r>  )?5ep@XU-c}YR@r2  r  r  r5  r7  s
             r   (test_uranus_barycenter_topocentric_date3r  	  s   s#)))4A'NE5ktDDDwwqz!!%(;"<=FFHH (BXBHH)6L+@ACKK*FY,>? V4BXBHH(&<*?@CKK+Vi-?@ (CXCKK*FY,>?BJJ)6I+=> z2CXCKK*FY,>?BJJ)6I+=> tV4CXCKK*FY,>?BJJ)6I+=>r   c                 
   t        j                  d      j                  g d      }| d   }|t        ddd      z   }|j	                  |      j                  | d	         j                         }|j                         \  }}}t        |j                  d
dt        z         t        |j                  ddt        z         |j                  d      \  }}}t        |j                  ddt        z         t        |j                  ddt        z         |j                         \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         |j                  d      \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         |j                  dd      \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         y )Nr+   r   r  r)   r6  r7  r8  r9  rC  )r  r  r  r  r,  )r  r  r  r  r|  r>  )r  r  r  r  )r  r  r   r  )r  r  r  r	  )r  r  r  r
  r2  )r  r  r  r	  r  r  )r  r  r  r	  r5  r7  s
             r   (test_uranus_barycenter_topocentric_date4r  	  s   s#))*^_A'NE5ktDDDwwqz!!%(;"<=FFHH (BXBHHfhnq}h}~CKKjlru~l~ V4BXBHHbdjmydyzCKKlnt  xA  oA  B (CXCKKgior{i{|BJJegmpygyz z2CXCKKhjps|j|}BJJegmpygyz tV4CXCKKgior{i{|BJJegmpygyzr   c                    t        j                  d      j                  d      }| d   }|t        ddd      z   }|j	                  |      j                  | d	         j                         }|j                         \  }}}t        |j                  d
dt        z         t        |j                  ddt        z         |j                  d      \  }}}t        |j                  ddt        z         t        |j                  ddt        z         |j                         \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         |j                  d      \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         |j                  dd      \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         y )Nr+   r   r5   r)   r6  r7  r8  r9  rj  &XF/@r,  ʉ41r|  r>  pY7/@k\u1pmP.>z@AR]@r2  ilD @r  r  dS!@r5  r7  s
             r   )test_neptune_barycenter_topocentric_date0r  	  s   s#))*;<A'NE5ktDDDwwqz!!%(<"=>GGIH (BXBHH(&<*?@CKK,fy.@A V4BXBHH(&<*?@CKK,fy.@A (CXCKK)6I+=>BJJ*FY,>? z2CXCKK*FY,>?BJJ*FY,>? tV4CXCKK*FY,>?BJJ*FY,>?r   c                    t        j                  d      j                  d      }| d   }|t        ddd      z   }|j	                  |      j                  | d	         j                         }|j                         \  }}}t        |j                  d
dt        z         t        |j                  ddt        z         |j                  d      \  }}}t        |j                  ddt        z         t        |j                  ddt        z         |j                         \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         |j                  d      \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         |j                  dd      \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         y )Nr+   r   r=   r)   r6  r7  r8  r9  rj  \jY	3@r,  ;5r|  r>  ^.YI3@&C5(r7DJqRk,U@r2  r  r  r5  r7  s
             r   )test_neptune_barycenter_topocentric_date1r  	  s   s#)))4A'NE5ktDDDwwqz!!%(<"=>GGIH (BXBHH(&<*?@CKK,fy.@A V4BXBHH(&<*?@CKK,fy.@A (CXCKK+Vi-?@BJJ)6I+=> z2CXCKK+Vi-?@BJJ)6I+=> tV4CXCKK+Vi-?@BJJ)6I+=>r   c                    t        j                  d      j                  d      }| d   }|t        ddd      z   }|j	                  |      j                  | d	         j                         }|j                         \  }}}t        |j                  d
dt        z         t        |j                  ddt        z         |j                  d      \  }}}t        |j                  ddt        z         t        |j                  ddt        z         |j                         \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         |j                  d      \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         |j                  dd      \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         y )Nr+   r   rA   r)   r6  r7  r8  r9  rj  ۬f\4@r,  0ϲ63r|  r>  $䒸\4@&u̗63&5.cX@r2  r  r  r5  r7  s
             r   )test_neptune_barycenter_topocentric_date2r'  
  s   s#)))4A'NE5ktDDDwwqz!!%(<"=>GGIH (BXBHH(&<*?@CKK,fy.@A V4BXBHH',)>?CKK+Vi-?@ (CXCKK,fy.@ABJJ)6I+=> z2CXCKK,fy.@ABJJ)6I+=> tV4CXCKK,fy.@ABJJ)6I+=>r   c                    t        j                  d      j                  d      }| d   }|t        ddd      z   }|j	                  |      j                  | d	         j                         }|j                         \  }}}t        |j                  d
dt        z         t        |j                  ddt        z         |j                  d      \  }}}t        |j                  ddt        z         t        |j                  ddt        z         |j                         \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         |j                  d      \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         |j                  dd      \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         y )Nr+   r   rD   r)   r6  r7  r8  r9  rj  4@6@r,  V`'r|  r>  TC6@`Dz&`FU@-NʳZ@r2  ` nm_@r  r  @7@d@r5  r7  s
             r   )test_neptune_barycenter_topocentric_date3r1  0
  s   s#)))4A'NE5ktDDDwwqz!!%(<"=>GGIH (BXBHH(&<*?@CKK,fy.@A V4BXBHH',)>?CKK,fy.@A (CXCKK)6I+=>BJJ)6I+=> z2CXCKK+Vi-?@BJJ)6I+=> tV4CXCKK+Vi-?@BJJ)6I+=>r   c                 
   t        j                  d      j                  g d      }| d   }|t        ddd      z   }|j	                  |      j                  | d	         j                         }|j                         \  }}}t        |j                  d
dt        z         t        |j                  ddt        z         |j                  d      \  }}}t        |j                  ddt        z         t        |j                  ddt        z         |j                         \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         |j                  d      \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         |j                  dd      \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         y )Nr+   r   r  r)   r6  r7  r8  r9  rj  )r  r  r!  r)  r,  )r  r  r"  r*  r|  r>  )r  r  r#  r+  )r  r  r$  r,  )r  r  r%  r-  )r  r  r&  r.  r2  )r  r  r%  r/  r  r  )r  r  r%  r0  r5  r7  s
             r   )test_neptune_barycenter_topocentric_date4r3  J
  s   s#))*^_A'NE5ktDDDwwqz!!%(<"=>GGIH (BXBHHfhnq}h}~CKKmou  yB  pB  C V4BXBHHdflo{f{|CKKlnt  xA  oA  B (CXCKKfhnqzhz{BJJegmpygyz z2CXCKKikqt}k}~BJJegmpygyz tV4CXCKKikqt}k}~BJJegmpygyzr   c                    t        j                  d      j                  d      }| d   }|t        ddd      z   }|j	                  |      j                  | d	         j                         }|j                         \  }}}t        |j                  d
dt        z         t        |j                  ddt        z         |j                  d      \  }}}t        |j                  ddt        z         t        |j                  ddt        z         |j                         \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         |j                  d      \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         |j                  dd      \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         y )Nr+   r   r5   r)   r6  r7  r8  r9  r  BW(@r,  Iz0@r|  r>  r}'@Au70@@.P@qʁfb@r2  ߜ/P@r  r   /P@r5  r7  s
             r   'test_pluto_barycenter_topocentric_date0r=  d
  s   s#))*;<A'NE5ktDDDwwqz!!%(:";<EEGH (BXBHH(&<*?@CKK+Vi-?@ V4BXBHH(&<*?@CKK+Vi-?@ (CXCKK*FY,>?BJJ)6I+=> z2CXCKK*FY,>?BJJ)6I+=> tV4CXCKK*FY,>?BJJ)6I+=>r   c                    t        j                  d      j                  d      }| d   }|t        ddd      z   }|j	                  |      j                  | d	         j                         }|j                         \  }}}t        |j                  d
dt        z         t        |j                  ddt        z         |j                  d      \  }}}t        |j                  ddt        z         t        |j                  ddt        z         |j                         \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         |j                  d      \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         |j                  dd      \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         y )Nr+   r   r=   r)   r6  r7  r8  r9  r  \n.@r,  u7_r|  r>  8?
j.@7A/c;0@8_^YZ@r2  DI0@r  r  \z{J0@r5  r7  s
             r   'test_pluto_barycenter_topocentric_date1rG  ~
  s   s#)))4A'NE5ktDDDwwqz!!%(:";<EEGH (BXBHH(&<*?@CKK+Vi-?@ V4BXBHH(&<*?@CKK,fy.@A (CXCKK+Vi-?@BJJ)6I+=> z2CXCKK*FY,>?BJJ)6I+=> tV4CXCKK+Vi-?@BJJ)6I+=>r   c                    t        j                  d      j                  d      }| d   }|t        ddd      z   }|j	                  |      j                  | d	         j                         }|j                         \  }}}t        |j                  d
dt        z         t        |j                  ddt        z         |j                  d      \  }}}t        |j                  ddt        z         t        |j                  ddt        z         |j                         \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         |j                  d      \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         |j                  dd      \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         y )Nr+   r   rA   r)   r6  r7  r8  r9  r  ̩J0@r,  ^&r|  r>  縵|0@#?&,kht6@&_@r2  @d6@r  r  ~6@r5  r7  s
             r   'test_pluto_barycenter_topocentric_date2rQ  
  s   s#)))4A'NE5ktDDDwwqz!!%(:";<EEGH (BXBHH(&<*?@CKK,fy.@A V4BXBHH',)>?CKK+Vi-?@ (CXCKK+Vi-?@BJJ*FY,>? z2CXCKK+Vi-?@BJJ*FY,>? tV4CXCKK*FY,>?BJJ*FY,>?r   c                    t        j                  d      j                  d      }| d   }|t        ddd      z   }|j	                  |      j                  | d	         j                         }|j                         \  }}}t        |j                  d
dt        z         t        |j                  ddt        z         |j                  d      \  }}}t        |j                  ddt        z         t        |j                  ddt        z         |j                         \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         |j                  d      \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         |j                  dd      \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         y )Nr+   r   rD   r)   r6  r7  r8  r9  r  mzQ}2@r,  \vA3r|  r>  /X2@?(3?V<@pQ?c@r2  )^<@r  r  n^<@r5  r7  s
             r   'test_pluto_barycenter_topocentric_date3r[  
  s   s#)))4A'NE5ktDDDwwqz!!%(:";<EEGH (BXBHH(&<*?@CKK,fy.@A V4BXBHH(&<*?@CKK,fy.@A (CXCKK*FY,>?BJJ*FY,>? z2CXCKK+Vi-?@BJJ*FY,>? tV4CXCKK+Vi-?@BJJ*FY,>?r   c                 
   t        j                  d      j                  g d      }| d   }|t        ddd      z   }|j	                  |      j                  | d	         j                         }|j                         \  }}}t        |j                  d
dt        z         t        |j                  ddt        z         |j                  d      \  }}}t        |j                  ddt        z         t        |j                  ddt        z         |j                         \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         |j                  d      \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         |j                  dd      \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         y )Nr+   r   r  r)   r6  r7  r8  r9  r  )r5  r?  rI  rS  r,  )r6  r@  rJ  rT  r|  r>  )r7  rA  rK  rU  )r8  rB  rL  rV  )r9  rC  rM  rW  )r:  rD  rN  rX  r2  )r;  rE  rO  rY  r  r  )r<  rF  rP  rZ  r5  r7  s
             r   'test_pluto_barycenter_topocentric_date4r]  
  s   s#))*^_A'NE5ktDDDwwqz!!%(:";<EEGH (BXBHHfhnq}h}~CKKkmsvm  A V4BXBHHegmp|g|}CKKkmsvm  A (CXCKKgior{i{|BJJfhnqzhz{ z2CXCKKgior{i{|BJJfhnqzhz{ tV4CXCKKgior{i{|BJJfhnqzhz{r   c                    t        j                  d      j                  d      }| d   }|t        ddd      z   }|j	                  |      j                  | d	         j                         }|j                         \  }}}t        |j                  d
dt        z         t        |j                  ddt        z         |j                  d      \  }}}t        |j                  ddt        z         t        |j                  ddt        z         |j                         \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         |j                  d      \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         |j                  dd      \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         y )Nr+   r   r5   r)   r6  r7  r8  r9  r  t' @r,  S&4@r|  r>  2q  @B14@_!\G@V(p@r2  bs^G@r  r  Ҭ'^G@r5  r7  s
             r   test_sun_topocentric_date0rg  
  s   s#))*;<A'NE5ktDDDwwqz!!%,/88:H (BXBHH&(=>CKK+Vi-?@ V4BXBHH',)>?CKK+Vi-?@ (CXCKK*FY,>?BJJ)6I+=> z2CXCKK*FY,>?BJJ)6I+=> tV4CXCKK*FY,>?BJJ)6I+=>r   c                    t        j                  d      j                  d      }| d   }|t        ddd      z   }|j	                  |      j                  | d	         j                         }|j                         \  }}}t        |j                  d
dt        z         t        |j                  ddt        z         |j                  d      \  }}}t        |j                  ddt        z         t        |j                  ddt        z         |j                         \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         |j                  d      \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         |j                  dd      \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         y )Nr+   r   r=   r)   r6  r7  r8  r9  r  Oh4@r,  a3@r|  r>  : !@/
3@@S]FM_r@r2  @@r  r  @@r5  r7  s
             r   test_sun_topocentric_date1rq     s   s#)))4A'NE5ktDDDwwqz!!%,/88:H (BXBHH(&<*?@CKK*FY,>? V4BXBHH(&<*?@CKK*FY,>? (CXCKK+Vi-?@BJJ*FY,>? z2CXCKK+Vi-?@BJJ*FY,>? tV4CXCKK*FY,>?BJJ*FY,>?r   c                    t        j                  d      j                  d      }| d   }|t        ddd      z   }|j	                  |      j                  | d	         j                         }|j                         \  }}}t        |j                  d
dt        z         t        |j                  ddt        z         |j                  d      \  }}}t        |j                  ddt        z         t        |j                  ddt        z         |j                         \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         |j                  d      \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         |j                  dd      \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         y )Nr+   r   rA   r)   r6  r7  r8  r9  r  :e2@r,  QZ
	7r|  r>  >q2@焏x7@fy. C|\@r2  r  r  r5  r7  s
             r   test_sun_topocentric_date2ry       s#)))4A'NE5ktDDDwwqz!!%,/88:H (BXBHH(&<*?@CKK+Vi-?@ V4BXBHH(&<*?@CKK+Vi-?@ (CXCKK+Vi-?@BJJ*FY,>? z2CXCKK+Vi-?@BJJ*FY,>? tV4CXCKK+Vi-?@BJJ*FY,>?r   c                    t        j                  d      j                  d      }| d   }|t        ddd      z   }|j	                  |      j                  | d	         j                         }|j                         \  }}}t        |j                  d
dt        z         t        |j                  ddt        z         |j                  d      \  }}}t        |j                  ddt        z         t        |j                  ddt        z         |j                         \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         |j                  d      \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         |j                  dd      \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         y )Nr+   r   rD   r)   r6  r7  r8  r9  r  W$@r,  :%@r|  r>  (F$@i_%@@JBچq@r2  r  r  r5  r7  s
             r   test_sun_topocentric_date3r  4  rz  r   c                 
   t        j                  d      j                  g d      }| d   }|t        ddd      z   }|j	                  |      j                  | d	         j                         }|j                         \  }}}t        |j                  d
dt        z         t        |j                  ddt        z         |j                  d      \  }}}t        |j                  ddt        z         t        |j                  ddt        z         |j                         \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         |j                  d      \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         |j                  dd      \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         y )Nr+   r   r  r)   r6  r7  r8  r9  r  )r_  ri  rs  r|  r,  )r`  rj  rt  r}  r|  r>  )ra  rk  ru  r~  )rb  rl  rv  r  )rc  rm  rw  r  )rd  rn  rx  r  r2  )re  ro  rw  r  r  r  )rf  rp  rw  r  r5  r7  s
             r   test_sun_topocentric_date4r  N  s   s#))*^_A'NE5ktDDDwwqz!!%,/88:H (BXBHHdflo{f{|CKKhjps|j|} V4BXBHHegmp|g|}CKKhjps|j|} (CXCKKhjps|j|}BJJgior{i{| z2CXCKKhjps|j|}BJJgior{i{| tV4CXCKKgior{i{|BJJgior{i{|r   c                    t        j                  d      j                  d      }| d   }|t        ddd      z   }|j	                  |      j                  | d	         j                         }|j                         \  }}}t        |j                  d
dt        z         t        |j                  ddt        z         |j                  d      \  }}}t        |j                  ddt        z         t        |j                  ddt        z         |j                         \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         |j                  d      \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         |j                  dd      \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         y )Nr+   r   r5   r)   r6  r7  r8  r9  r  usm(@r,  m/8Br|  r>  t~(@sk5E!D@ oNb@r2  _%D@r  r  cC+D@r5  r7  s
             r   test_moon_topocentric_date0r  h  s   s#))*;<A'NE5ktDDDwwqz!!%-099;H (BXBHH(&<*?@CKK+Vi-?@ V4BXBHH(&<*?@CKK+Vi-?@ (CXCKK*FY,>?BJJ*FY,>? z2CXCKK+Vi-?@BJJ*FY,>? tV4CXCKK+Vi-?@BJJ*FY,>?r   c                    t        j                  d      j                  d      }| d   }|t        ddd      z   }|j	                  |      j                  | d	         j                         }|j                         \  }}}t        |j                  d
dt        z         t        |j                  ddt        z         |j                  d      \  }}}t        |j                  ddt        z         t        |j                  ddt        z         |j                         \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         |j                  d      \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         |j                  dd      \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         y )Nr+   r   r=   r)   r6  r7  r8  r9  r  bc٩7@r,  j,s?r|  r>  @4̧7@燤?$y_G3= "u@r2  r  r  r5  r7  s
             r   test_moon_topocentric_date1r    s   s#)))4A'NE5ktDDDwwqz!!%-099;H (BXBHH(&<*?@CKK*FY,>? V4BXBHH(&<*?@CKK+Vi-?@ (CXCKK+Vi-?@BJJ*FY,>? z2CXCKK+Vi-?@BJJ*FY,>? tV4CXCKK+Vi-?@BJJ*FY,>?r   c                    t        j                  d      j                  d      }| d   }|t        ddd      z   }|j	                  |      j                  | d	         j                         }|j                         \  }}}t        |j                  d
dt        z         t        |j                  ddt        z         |j                  d      \  }}}t        |j                  ddt        z         t        |j                  ddt        z         |j                         \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         |j                  d      \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         |j                  dd      \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         y )Nr+   r   rA   r)   r6  r7  r8  r9  r  ?-@r,  l0.'r|  r>  vx,ް-@E-'$-O0B@U큉c@r2  9?>3B@r  r  ?3B@r5  r7  s
             r   test_moon_topocentric_date2r    s   s#)))4A'NE5ktDDDwwqz!!%-099;H (BXBHH(&<*?@CKK,fy.@A V4BXBHH(&<*?@CKK+Vi-?@ (CXCKK+Vi-?@BJJ)6I+=> z2CXCKK*FY,>?BJJ)6I+=> tV4CXCKK+Vi-?@BJJ)6I+=>r   c                    t        j                  d      j                  d      }| d   }|t        ddd      z   }|j	                  |      j                  | d	         j                         }|j                         \  }}}t        |j                  d
dt        z         t        |j                  ddt        z         |j                  d      \  }}}t        |j                  ddt        z         t        |j                  ddt        z         |j                         \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         |j                  d      \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         |j                  dd      \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         y )Nr+   r   rD   r)   r6  r7  r8  r9  r  ^`g|a0@r,  ^3]5r|  r>  璾d0@)Bp5vѺ{p<@nmpg@r2  \r2x<@r  r  lFx<@r5  r7  s
             r   test_moon_topocentric_date3r    s   s#)))4A'NE5ktDDDwwqz!!%-099;H (BXBHH(&<*?@CKK+Vi-?@ V4BXBHH(&<*?@CKK+Vi-?@ (CXCKK+Vi-?@BJJ*FY,>? z2CXCKK*FY,>?BJJ*FY,>? tV4CXCKK*FY,>?BJJ*FY,>?r   c                 
   t        j                  d      j                  g d      }| d   }|t        ddd      z   }|j	                  |      j                  | d	         j                         }|j                         \  }}}t        |j                  d
dt        z         t        |j                  ddt        z         |j                  d      \  }}}t        |j                  ddt        z         t        |j                  ddt        z         |j                         \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         |j                  d      \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         |j                  dd      \  }}	}t        |j                  ddt        z         t        |	j                  ddt        z         y )Nr+   r   r  r)   r6  r7  r8  r9  r  )r  r  r  r  r,  )r  r  r  r  r|  r>  )r  r  r  r  )r  r  r  r  )r  r  r  r  )r  r  r  r  r2  )r  r  r  r  r  r  )r  r  r  r  r5  r7  s
             r   test_moon_topocentric_date4r    s   s#))*^_A'NE5ktDDDwwqz!!%-099;H (BXBHHfhnq}h}~CKKikqt}k}~ V4BXBHHfhnq}h}~CKKikqt}k}~ (CXCKKhjps|j|}BJJgior{i{| z2CXCKKgior{i{|BJJgior{i{| tV4CXCKKhjps|j|}BJJgior{i{|r   c                    d}t        j                  t        |            }t        j                  j                  |j                  d         }| j                  t        j                         j                  d            j                  |      j                         \  }}}t        |j                  ddt        z         t        |j                   ddt"        z         y )N  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
r   r5   ke:@rG   PV@r   load_dataframer   r   r   from_dataframeilocrC  r   r   r   r  rE  r"   rF  r  rG  r:   r)   linedfr   rK  rL  rM  s          r   test_hipparcos_conversion0r    s     RD		!	!'$-	0B<<&&rwwqz2D!1!7!78I!JKSSTXY__aBXBHH('L*@ACKK*Gi,?@r   c                    d}t        j                  t        |            }t        j                  j                  |j                  d         }| j                  t        j                         j                  d            j                  |      j                         \  }}}t        |j                  ddt        z         t        |j                   ddt"        z         y )Nr  r   r=   <@rG   bFPV@r  r  s          r   test_hipparcos_conversion1r    s     RD		!	!'$-	0B<<&&rwwqz2D!1!7!7	!BCKKDQWWYBXBHH'<)?@CKK*Gi,?@r   c                    d}t        j                  t        |            }t        j                  j                  |j                  d         }| j                  t        j                         j                  d            j                  |      j                         \  }}}t        |j                  ddt        z         t        |j                   ddt"        z         y )Nr  r   rA   Vn	>@rG   5o/PV@r  r  s          r   test_hipparcos_conversion2r         RD		!	!'$-	0B<<&&rwwqz2D!1!7!7	!BCKKDQWWYBXBHH('L*@ACKK*Gi,?@r   c                    d}t        j                  t        |            }t        j                  j                  |j                  d         }| j                  t        j                         j                  d            j                  |      j                         \  }}}t        |j                  ddt        z         t        |j                   ddt"        z         y )Nr  r   rD   B&O?@rG   j>PV@r  r  s          r   test_hipparcos_conversion3r    r  r   c                    d}t        j                  t        |            }t        j                  j                  |j                  d         }| j                  t        j                         j                  g d            j                  |      j                         \  }}}t        |j                  ddt        z         t        |j                   ddt"        z         y )Nr  r   r  )r  r  r  r  rG   )r  r  r  r  r  r  s          r   test_hipparcos_conversion4r  
  s     RD		!	!'$-	0B<<&&rwwqz2D!1!7!78l!mnvvw{|  C  C  EBXBHHegnq}g}~CKKegnqzgz{r   N)__doc__numpyr   r   r   r   skyfieldr   r   r	   r
   r   r   skyfield.apir   r   skyfield.constantsr   r   skyfield.datar   skyfield.functionsr   r   fixesr   r   r  r   r:   r  r   r   r"   r&   r)   r3   r;   r>   rB   rE   rL   rN   rP   rR   rW   rY   r\   r^   re   rh   rm   rp   rr   rt   ry   r|   r~   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r   r  r  r  r	  r  r  r  r  r  r  r   r#  r&  r(  r*  r.  r1  r4  rN  rQ  rS  rU  rW  rZ  r\  r^  rc  re  rh  rj  rq  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r$  r-  r6  r?  rA  rK  rT  r]  rf  rh  rr  r{  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r'  r)  r9  rA  rK  rS  rU  r_  rg  rq  ry  r{  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r'  r1  r3  r=  rG  rQ  r[  r]  rg  rq  ry  r  r  r  r  r  r  r  r  r  r  r  r  r   r   r   <module>r     s   D ) )7 7 $ * # 1 $			%
	)d
6""""MMOKLKJK}SSMSCBFDCDED#$$CG
CG
BF
CE
&&&&55556677&&&&*)97666)232323234232EA0A0B0A0@0A0@0A0@0}0B0A0B0B0B0A0A0@0@0}0?0B0A0A0~0B0B0B0@0A0B0B0A0B0B0A0B0A0B0A0@0A0B0A00A0A0B0A0@0A(A(A(@({(@4?4@4?4|4?4@4?4>4z4@4?4?4?4{4@4@4?4?4|4@4>4@4@4|4@4?4?4?4{4@4?4?4?4{4?4?4@4@4|4?4@4@4@4}4@4@4?4@4}4AAAA|r   