R1i$dZddlmZddlmZddlmZddlm Z ddl m Z ddl m Z dd lmZdd lmZGd d eZGd deZGddeZdZdZy)z'Vector functions and their composition.) target_names)max)C_AUDAY)reifyDeprecationError) length_of)build_position)TimeceZdZdZdZedZedZedZdZ dZ dZ d Z d Z d Zd Zd ZdZ ddZdZy)VectorFunctionz0Given a time, computes a corresponding position.Nc,t|jSN)type__name__selfs V/home/cursorai/projects/iching/venv/lib/python3.12/site-packages/skyfield/vectorlib.py vector_namezVectorFunction.vector_namesDz"""c,t|jSr) _jpl_namecenterrs r center_namezVectorFunction.center_name%%rc,t|jSr)rtargetrs r target_namezVectorFunction.target_namerrc^djt|jt|S)Nz <{0} {1}>)formatrrstrrs r__repr__zVectorFunction.__repr__s#!!$t*"5"5s4yAArcV|j|ur |jS|jSr)rr arrow_strrs r__str__zVectorFunction.__str__ s' ;;$ ## #~~rcddj|j|j|jS)Nz{0} {1} -> {2})r!rrrrs rr%zVectorFunction.arrow_str%s.&&   d..0@0@  rc|j|jk7r)|j|jk(r||}}n tdt|ddxs|f}t|ddxs|f}t |j|j||zS)Nztyou can only add two vectors if the target where one of the vectors ends is the center where the other vector startsvector_functions)rr ValueErrorgetattr VectorSum)rotherself_vfs other_vfss r__add__zVectorFunction.__add__*s ;;%,, &||t{{*#Te C 4!3T:EtgE#5t<H ellHy4HIIrct|Sr)ReversedVectorrs r__neg__zVectorFunction.__neg__:s d##rc*|j|jk7r tdt|ddxs|f}t|ddxs|f}tt |Dcgc]}| c}}t |j |j ||zScc}w)NzGyou can only subtract two vectors if they both start at the same centerr))rr*r+tuplereversedr,r)rr-r.r/vfs r__sub__zVectorFunction.__sub__=s ;;%,, &9  4!3T:EtgE#5t<H ()#<)BRC)#<=> t{{I4HII$=s B ct|tstdj||j |\}}}}|j }t |||||j}|j|_ ||_ ||_ |S)aIAt time ``t``, compute the target's position relative to the center. If ``t`` is an array of times, then the returned position object will specify as many positions as there were times. The kind of position returned depends on the value of the ``center`` attribute: * Solar System Barycenter: :class:`~skyfield.positionlib.Barycentric` * Center of the Earth: :class:`~skyfield.positionlib.Geocentric` * Difference: :class:`~skyfield.positionlib.Geometric` * Anything else: :class:`~skyfield.positionlib.ICRF` z_please provide the at() method with a Time instance as its argument, instead of the value {0!r}) isinstancer r*r!_atrr r ephemeris _ephemeris_observer_gcrs_aumessage)rtpv gcrs_positionr?rpositions ratzVectorFunction.atJs!T",,2F1I7 7(,xx{$1mW!!Q64;;?"nn%2""rc|jdk7r$tdj|jt ||S)Nrztyou can only observe() a body whose vector center is the Solar System Barycenter, but this vector has the center {0})rr*r!r_correct_for_light_travel_time)robservers r_observe_from_bcrsz!VectorFunction._observe_from_bcrsdsB ;;! C%fT%5%568 8.h==rctd)Nathe geometry_of() method has, alas, been deprecated This old method has been replaced by an improved interface. If you just need your software working again, install Skyfield 0.9.1 for a quick fix: pip install skyfield==0.9.1 Or, to update your old code, replace each operation that looks like: position = boston.geometry_of(satellite).at(t) with the vector math that was previously hiding inside the old method: position = (satellite - boston).at(t)r)rr-s r geometry_ofzVectorFunction.geometry_ofls -. .rctd)Nathe topos() method has, alas, been deprecated This old method has been replaced by an improved interface. If you just need your software working again, install Skyfield 0.9.1 for a quick fix: pip install skyfield==0.9.1 Or, to update your old code, replace each operation that looks like: boston = earth.topos(...) with the vector math that was previously hiding inside the old method: from skyfield.api import Topos boston = earth + Topos(...)r)rlatitude longitudelatitude_degreeslongitude_degrees elevation_mxys rtoposzVectorFunction.topos}s#$ $rctd)NaMthe satellite() method has, alas, been deprecated This old method has been replaced by an improved interface. If you just need your software working again, install Skyfield 0.9.1 for a quick fix: pip install skyfield==0.9.1 Or, to update your old code, replace each operation that looks like: sat = earth.satellite(tle_text) with the vector math (and the little bit of text manipulation) that was previously hiding inside the old method: from skyfield.api import EarthSatellite line1, line2 = tle_text.splitlines()[-2:] sat = earth + EarthSatellite(line1, line2)r)rtexts r satellitezVectorFunction.satellites23 3r)NNNNrXrX)r __module__ __qualname____doc__r<rrrrr#r&r%r0r3r8rErIrKrTrWrrrr s:I # # & & & &B  J $ J4>."EI@C$&3rrcNeZdZdZedZedZedZdZdZ y)r2cV|j|_|j|_||_yr)rrvector_function)rr_s r__init__zReversedVector.__init__s$%,, %,, .rc4d|jjzS)Nz Reversed )r_rrs rrzReversedVector.vector_namesT11====rc.|jjSr)r_rrs rrzReversedVector.center_name##///rc.|jjSr)r_rrs rrzReversedVector.target_namercrc|jSr)r_rs rr3zReversedVector.__neg__s###rcR|jj|\}}}}| | d|fSr)r_r;)rr@rArB_r?s rr;zReversedVector._ats4//33A61arA2tW$$rN) rrYrZr`rrrrr3r;r\rrr2r2sM/  > > 0 0 0 0$%rr2c$eZdZdZdZdZdZy)r,cb||_||_||_d|D}t|d|_y)Nc3LK|]}|jr|jywr)r<).0segments r z%VectorSum.__init__..s'-8HW!++((8Hs"$)rrr)nextr<)rrrr) ephemeridess rr`zVectorSum.__init__s5  0-8H- k40rc|j}|Dcgc]}d|jz}}djt|dj |Scc}w)N zSum of {0} vectors: {1} )r)r%r!lenjoin)rr)rlliness rr&zVectorSum.__str__s`00:JK:Jww((**:JK)00  ! IIe   LsAc$dj|S)Nz )r!rs rr#zVectorSum.__repr__s##D))rcd\}}d}|j}|D]4}|j|\}}} } |jdk(r| }||z }||z }6|djdk(rjdk(r}||| fS)N)rXrXir)r)r;r) rr@rArBrCvfsr7p2v2rgr?s rr;z VectorSum._ats1 ##B!# BAwyyC!" GA GA  q6==A "))s"2M!]G++rN)rrYrZr`r&r#r;r\rrr,r,s1 * ,rr,c|j}|j}|j}|j}|jj }|j j}|j|\}} } } t||z } d} tdD]d}| tz }|| z }tt|dkrnH|j|||z }|j|\}} } } t||z } |} ftd||z | |z ||fS)aReturn a light-time corrected astrometric position and velocity. Given an `observer` that is a `Barycentric` position somewhere in the solar system, compute where in the sky they will see the body `target`, by computing the light-time between them and figuring out where `target` was back when the light was leaving it that is now reaching the eyes or instruments of the `observer`. rX g-q=z$light-travel time failed to converge)r@tswhole tdb_fractionrDauvelocityau_per_dr;r rangerrabstdb_jdr*)rHrr@r}r~r cposition cvelocity tposition tvelocityrCr?distance light_time0i light_timedeltat2s rrGrGs  A B GGE>>L!!$$I!!**I39::a=0Iy-Y./HK 2Y' [( s5z?U "  YYulZ7 87=zz"~4 9mWY23 ?@@ y )i"7J FFrct|tst|jSt j |}| t |Sdj||S)Nz{0} {1})r:intrr_jpl_code_name_dictgetr"r!)rnames rrr sN fc "F|$$$  " "6 *D |6{   FD ))rN)r[jplephem.namesrrnumpyr constantsr descriptorlibrerrorsr functionsr positionlibr timelibr objectrr2r,rGrr\rrrsT-> $ 'V3VV3p%^%2",",H&GP*r