
    R1i                     J    d dl mZ d dlmZmZ d dlmZ d Zd Zd Z	d Z
d Zy	)
    )assert_raises)Toposload)
Geocentricc                      t        d      } | d   }| d   }t        t        d      5  ||z    d d d        y # 1 sw Y   y xY w)N	de421.bspearthmarsz(the center where the other vector starts)r   r   
ValueError)planetsr	   r
   s      _/home/cursorai/projects/iching/venv/lib/python3.12/site-packages/skyfield/tests/test_vectors.pytest_bad_additionr      s>    ;GGE6?D	z#M	N 
O	N	Ns   5>c                      t        d      } | d   }t        ddd      }t        t        d      5  ||z
   d d d        y # 1 sw Y   y xY w)Nr   r	   	38.9215 N	77.0669 W      W@elevation_mz%if they both start at the same center)r   r   r   r   )r   r	   usnos      r   test_bad_subtractionr      sC    ;GGEkt<D	z#J	K 
L	K	Ks	   >Ac                  p    t        d      } | d   | d   z
  }t        |      dk(  sJ t        |      dk(  sJ y )Nr   zearth barycentersunzSum of 2 vectors:
 Reversed 'de421.bsp' segment 10 SUN -> 0 SOLAR SYSTEM BARYCENTER
 'de421.bsp' segment 0 SOLAR SYSTEM BARYCENTER -> 3 EARTH BARYCENTERz<VectorSum of 2 vectors:
 Reversed 'de421.bsp' segment 10 SUN -> 0 SOLAR SYSTEM BARYCENTER
 'de421.bsp' segment 0 SOLAR SYSTEM BARYCENTER -> 3 EARTH BARYCENTER>)r   strrepr)r   vs     r   test_chebyshev_subtractionr      s`    ;G"#gen4Aq6 H H H H
 7 I I I I    c                     t        j                         } | j                  ddddd      }t        ddd	      }| }|j	                  |      }|j	                  |      }|j
                  j                  |j
                  j                   k(  j                         sJ |j                  j                  |j                  j                   k(  j                         sJ | }||u sJ y )
Ni              r   r   r   r   )
r   	timescaleutcr   atpositionauallvelocityau_per_d)tstr   negp1p2s         r   test_negationr0   #   s    		B
tQB"Akt<D%C	B	BKKNN..33555KK  bkk&:&:$::??AAA $C$;;r   c                  f   t        j                         } | j                  ddddd      }t        d      }|d   }|d   }|}t        |      d	k(  sJ t	        |      d
k(  sJ t        |j                  |            dk(  sJ |t        ddd      z   }t        |      dk(  sJ t	        |      dk(  sJ t        |j                  |            dk(  sJ ||z
  }t        |      dk(  sJ t	        |      dk(  sJ t        |j                  |            dk(  sJ t        g d      }|j                  dk(  sJ y )Ni        
   ,   r   r	   r
   zSum of 2 vectors:
 'de421.bsp' segment 0 SOLAR SYSTEM BARYCENTER -> 3 EARTH BARYCENTER
 'de421.bsp' segment 3 EARTH BARYCENTER -> 399 EARTHz<VectorSum of 2 vectors:
 'de421.bsp' segment 0 SOLAR SYSTEM BARYCENTER -> 3 EARTH BARYCENTER
 'de421.bsp' segment 3 EARTH BARYCENTER -> 399 EARTH>zF<Barycentric BCRS position and velocity at date t center=0 target=399>r   r   r   r   zSum of 3 vectors:
 'de421.bsp' segment 0 SOLAR SYSTEM BARYCENTER -> 3 EARTH BARYCENTER
 'de421.bsp' segment 3 EARTH BARYCENTER -> 399 EARTH
 Geodetic 399 EARTH -> IERS2010 latitude +38.9215 N longitude -77.0669 E elevation 92.0 mz<VectorSum of 3 vectors:
 'de421.bsp' segment 0 SOLAR SYSTEM BARYCENTER -> 3 EARTH BARYCENTER
 'de421.bsp' segment 3 EARTH BARYCENTER -> 399 EARTH
 Geodetic 399 EARTH -> IERS2010 latitude +38.9215 N longitude -77.0669 E elevation 92.0 m>z<Barycentric BCRS position and velocity at date t center=0 target=IERS2010 latitude +38.9215 N longitude -77.0669 E elevation 92.0 m>a  Sum of 4 vectors:
 Reversed 'de421.bsp' segment 499 MARS -> 4 MARS BARYCENTER
 Reversed 'de421.bsp' segment 4 MARS BARYCENTER -> 0 SOLAR SYSTEM BARYCENTER
 'de421.bsp' segment 0 SOLAR SYSTEM BARYCENTER -> 3 EARTH BARYCENTER
 'de421.bsp' segment 3 EARTH BARYCENTER -> 399 EARTHa  <VectorSum of 4 vectors:
 Reversed 'de421.bsp' segment 499 MARS -> 4 MARS BARYCENTER
 Reversed 'de421.bsp' segment 4 MARS BARYCENTER -> 0 SOLAR SYSTEM BARYCENTER
 'de421.bsp' segment 0 SOLAR SYSTEM BARYCENTER -> 3 EARTH BARYCENTER
 'de421.bsp' segment 3 EARTH BARYCENTER -> 399 EARTH>z<<ICRF position and velocity at date t center=499 target=399>)r   r   r   i  )	r   r#   ttr   r   r%   r   r   center)r+   r,   r   r	   r
   r   
geocentrics          r   test_vectorsr9   1   s   		B
dAr2r"A;GGE6?DAq6 8 8 8 8
 7 9 9 9 9
 qttAw< H H H H 	k;DAAAq6 ] ] ] ] 7 ^ ^ ^ ^ qttAw< I I I I 	Aq6 8 8 8 8 7 9 9 9 9 qttAw< > > > > G$J###r   N)assayr   skyfield.apir   r   skyfield.positionlibr   r   r   r   r0   r9    r   r   <module>r>      s*      $ +I<$r   