ChOl <dZddlmZddlmZddlmZddlZ edgdZ edgd Z d Z e je jzZe j e _e j$e j&e j(e j*e j,e j.e j0e j2e j"e j4e j6g Zd Zd Zd ZdZ dZ!dZ"dZ#dZ$dZ%dZ&dZ'dZ(dZ)dZ*dZ+dZ,dZ-dfdZ.dZ/dZ0dZ1d Z2d!Z3d"\Z4Z5Z6d#Z7e jpe4e5e6e7Z9d$Z:e3e9e:Z;e;red*Z?d+Z@d,ZAd-ZBd.ZCd/ZDd0ZEd1ZFd2ZGdgd3ZHd4ZId5ZJd6ZKd7ZLd8ZMd9ZNdhd:ZOd;ZPd<ZQd=ZRd>ZSd?ZTd@ZUdAZVdBZWdCZXdDZYdEZZdFZ[dGZ\dHZ]dIZ^dJZ_e`dKk(rddlaZae dLdMdNZbe dOdPdNZce dQdRdSZde dTdUdNZee,e dVdWdXZfe,e dYd'dZZge,e d[d\d]Zhe,e dVd\d^Zie,e d_d`daZje,e d_d`dbZke,e dcdbdaZle,e dddWdeZme<e@eleee<e@emeeyy)iz9 Use Swiss ephemeris to calculate tithi, nakshatra, etc. )division)ceil) namedtupleNDate)yearmonthdayPlace)latitude longitudetimezonegpa v@cLtjtjddS)NgWf=Arswe set_sid_mode SIDM_USERd:\jyotish\support\panchanga.pyr's((8I1MrcLtjtjddS)Ng S=Arrrrrrr(sS--cmm=NPQRrcHtjtjSN)rr SIDM_LAHIRIrrrrr*sC,,S__=rcHtjtjSr)rrSIDM_FAGAN_BRADLEYrrrrr+sc..s/E/EFrcPtjdtjdtjdtjdtj dtj dtjdtjdtjd tjd i }||S) NSuryaChandraMangalaBudhaGuruShukraShaniRahuKetu) rSUNMOONMARSMERCURYJUPITERVENUSSATURN MEAN_NODEKETUPLUTO)planetnamess rget_planet_namer3.sj GGWchh 388Y KK#++vsyy( JJ&#))U[ ]% vrc||dz z|dz zS)N<ir)degsminssecss rrr5sD47NT$Y$>rclt|}||z dz}t|}t||z dzd}|||gS)Nr5)intround)degdr7mss r to_dms_precrA8sA #h! 'R$ $i! TAXOQ! Qrc<t|\}}}||t|gSr)rAr;)r=r>r?r@s rto_dmsrC?s"  '!Q QArc|}tdt|D]}||||dz ks||xxdz cc< |t|k(sJ|S)zdAdd 360 to those elements in the input list so that all elements are sorted in ascending order.h)rangelensorted)anglesresultis r unwrap_anglesrMCsX & CK a ay6!A#;q S 0  ! 6&> !" ! -rc|dk\r|dz S|SNrFrangles rrrNs#@5@rc |dzSNrFrrQs rrrQs rc|dzdzSrOrrahus rrrUs TCZ3&rctjtj|dtjd|tjtj zddz }|S)NCitra)flagrrP)rrr fixstar_ut FLG_SWIEPH FLG_SIDEREAL)pointfvals rfunctionraWsNS]]E3/ >>'5#BRBR1R STU VZ] ^D Krc~|}|}d} ||zdz }||}||}||zdk\r|}n|}||z |krn/||zdz S)Ng& .>rr) funcstartstopleftrightepsilonmiddlemidvalrtvals rbisection_searchrmesn $ % ' Ula F6lF KE ~e d   %  $,! rct|t|k(sJd}tt|D]L}d}d}tt|D]!}||k7s ||||z z}|||||z z}#||||z|z z }N|S)zOGiven two lists x and y, find the value of x = xa when y = ya, i.e., f(xa) = yarrE)rHrG)xyyatotalrLnumerdenomjs rinverse_lagrangervws Q3q6   % Q=a E E 3q6] a "qt) !A$1+  UQqT\E !!E  ,rcntj|j|j|jdS)NrY)rjuldayrrr )dates rrrsszz$))TZZ3OrcJtj|tjSr)rrevjulGREG_CALjds rrrsSZZCLL9rc tj|||||||\}}} } } } tj||| | | dtj\} }|S)zConverts local time to JD(UT1)r)r utc_time_zone utc_to_jdr|)rrr hourminutessecondsr rpr?r>hmntr@jd_etjd_ut1s rlocal_time_to_jdut1rsS((uc4'S[\!Q1c1--1aCCLLA-% -rcfd}d}t||z }|||zz }t||z }d|zd|zgS)zFGives nakshatra (1..27) and paada (1..4) in which given longitude lies竪*@竪 @rEr;)r one_starone_padaquotientreminderpadas rnakshatra_padarsN(( X% &((X--( X !$ h,D !!rcttj||tjtjz\}}t t |dS)z.to_dmssXd)%+,'$,",w6"<='#;a}Agc]; ;rz(Error calculating sunrise. Status code: N)r CALC_RISEBIT_DISC_CENTER rise_transr'r]r;print) r~placelatlontzgeoposrsmirKres_coderise rise_localrCs rsunrisers,#sB a=& ,, ,$ >>2b5j gg   nn &AY( ] AYq\d"T'>j< &*s:">"!DE FF 6xj ABD\r)i g33333&@)g?ܵL6@g cQ@@zSunrise Julian Day (Local): zSunrise Time (HH:MM:SS): rEc |\}}}tj||dz z tj||tttj z}|dd}||dz zt ||z dz|zgS)zASunset when centre of disc is at horizon for given date and placerrrErr)rrr'r; _rise_flagsCALC_SETrCr~rrrrrKsettings rsunsetrst,#sB >>"r"u*cggsCKRUR^R^D^@_ `& 1IaL' BsF FGbLB#6#;< ==rc |\}}}tj||dz z tj||tttj z}|dd}t ||z dz|zS)zCMoonrise when centre of disc is at horizon for given date and placerrrEr)rrr(r;rrrC)r~rrrrrKrs rmoonrisersf,#sB >>"r"u*chhS[SVS`S`E`Aa b& 1$ b 2% &&rc |\}}}tj||dz z tj||tttj z}|dd}t ||z dz|zS)zBMoonset when centre of disc is at horizon for given date and placerrrEr)rrr(r;rrrCrs rmoonsetrsf,#sB >>"r"u*chhS[SVS_S_E_A` a& 1IaL' 2#b( ))rc|j}t||d|dz z }t|}t|dz }|dz|z }gd}|Dcgc]}t ||zt |z dz!} }|Dcgc]}t ||zt |z dz!} }t | | D cgc] \} } | | z  } } } | }|}t|||}||z|z dz|z}t|t|g}t|dz}t|dz }||z dzdkD}|rT|dz}|dz|z }t|||}||z|z dz|jz}|dk(rdn|}|t|t|gz }|Scc}wcc}wcc} } w)zITithi at sunrise for given date and place. Also returns tithi's end time.rr ????rFrE) r r lunar_phaserlunar_longitudesolar_longitudeziprvr;rC)r~rrr moon_phasetoday degrees_leftoffsetstlunar_long_diffsolar_long_diffmoonsunrelative_motionrpro approx_endendsanswermoon_phase_tmrwtomorrow isSkipped leap_tithis rtithirs ~~" U A b ($4 * zB %j(, #'Y`bY`TUtax0?43HHCOY`/bY`bY`TUtax0?43HHCOY`/b474Y[4Y[dCdSj4Y/[! !1l3*  R 2 % *$ Jt %& q)/ /B& '(%2%))J?Z/L!!Q 5J : r !R '%.. 8DrkzJ s:t --F -3cb[s $E05$E5*E:c|\}}}t||d|dz z }gd}|Dcgc]}t||z}}t|ddzdz } t|} |} t | | | dzdz } ||z | zdz|z} t | t | g}t|ddzdz }|| z dzdkD}|rH| dz}t |||dzdz } ||z | zdz|z} | dk(rdn|}|t |t | gz }|Scc}w) zYCurrent nakshatra as of julian day (jd) 1 = Asvini, 2 = Bharani, ..., 27 = Revati rr)rYrrrrrFrrE)rrrrMrvr;rC)r~rrrrrrr longitudesnakrprorrrnak_tmrwrleap_naks r nakshatrar sL ,#sB U A c )$ ''4;c>|\}}}t||d|dz z }t|}t|}||zdz}t|dzdz } | dz|z } gd} | D cgc]} t|| zt|z dz!} } | D cgc]} t|| zt|z dz!}} t | |Dcgc] \}}||z }}}|}| }t ||| }||z|z dz|z}t | t|g}t|dz}t|dz}||zdz}t|dzdz }|| z dzdkD}|rJ| dz}|dz|z } t ||| }||z|z dz|z}| dk(rdn|}|t |t|gz }|Scc} wcc} wcc}}w) zQYoga at given jd and place. 1 = Vishkambha, 2 = Priti, ..., 27 = Vaidhrti rrrFrrrrrE)rrrrrrvr;rC)r~rrrrr lunar_long solar_longrryogrrrrrrr total_motionrprorrrlunar_long_tmrwsolar_long_tmrw total_tmrwrrleap_yogs ryogarCs  ,#sB U A c )$t$*t$*  "c )% URZ# #!E), #'Y`bY`TUtax0?43HHCOY`/bY`bY`TUtax0?43HHCOY`/b14_o1VX1V+44#:1V,X! !1l3*  b B & +$ HfTl #&$D1H-/#D1H-//1S8* *r/C' ((#~#a')QwHx(50L!!Q 5J :  "b (2 -D2Iq8H s8}fTl++F -7cbXs$F$F:Fct||d}t|}t|}||z dz}t|dz }|dz|z }t |gS)z9Returns the karana and their ending times. (from 1 to 60)rrFr:)rrrrr;)r~rrrrrrrs rkaranarusb U A $t$*t$*Z'3.* zA~ %Z', e*rc6tt|dzdzS)zFWeekday for given Julian day. 0 = Sunday, 1 = Monday,..., 6 = SaturdayrE)r;rr}s rvaarars T"q&\A  rct||d}t||d}t||d}t||d}t|}t|}||k(}|dz} | dkDr| dz} t | |gS)zbReturns lunar month and if it is adhika or not. 1 = Chaitra, 2 = Vaisakha, ..., 12 = PhalgunarrrEr)rrnew_moonraasir;) r~rticritical last_new_moon next_new_moonthis_solar_monthnext_solar_month is_leap_monthmaasas rmasars R" R  "(8R,-8R,-=)=)#'77- Q % RZ%"* e*m $$rc |dz S)Ngc!Arr}s rrrsb8mrcjt|}t|d|z dzztz }|dz }|dz}||fS)Nrrik )aharganar; sidereal_year)r~ maasa_numaharkalisakavikramas r elapsed_yearr sD "$ da)mr))] : ;$ $ 3J' trc|dk(r||z }|dk(r|d|z z}tdDcgc] }d|dz z }}|Dcgc]}t|z}}t|}t||d}|zScc}wcc}w)zReturns JDN, where opt = -1: JDN < jd such that lunar_phase(JDN) = 360 degrees opt = +1: JDN >= jd such that lunar_phase(JDN) = 360 degrees rrErrrF)rGrrMrv) r~tithi_optreoffsetrorLrpy0s rrrs  BYfBYrF{+$)"I/I&VAX I!/'()q!{519q!)A!1c""  0)s A.A3cJt|}t|}t|dz S)z-Zodiac of given jd. 1 = Mesha, ... 12 = Meenag>@)rr)r~r@solar_nirayanas rrrs&b!"2&. ns" ##rcBt|}t|}||z dz}|SrT)rr)r~rrrs rrrs+r"*r"*Z'3.* rcxt||d}|dk\r|dz dz}|dzt|dzdz dz zdz}|S) Nrir5rliPF)r r;)r~rr samvats r samvatsararsT b) $Q '$ T\dRi2-d 2ITCZ#-67 72 =& -rc|dz dzS)z0 = Vasanta,...,5 = ShishirarErcr)masa_nums rritur s Q,1 rcht||d}t||d}||z dz}|t|gS)Nrr)rrrC)r~rsrisessetdiffs r day_durationr%sA "e Q % E 1 $ ," $ t rc (|\}}}|j}tj||dz z tj||t t tj zdd}tj||dz z tj||t t tjzdd}||z }g}tddD].} |jt|| |zdz z|z dz|z0tj|dz|dz z tj||t t tj zdd}||z } tddD].} |jt|| | zdz z|z dz|z0|S)NrrrEr ) r rrr'r;rrrrGappendrC) r~rrrrr"r#day_dur end_timesrL night_durs rgauri_chogadiyar-sy,#sB ~~" ..beSWWc3s;QTQ^Q^C^?_ `ab cde f% RU CGGS#c+PSP\P\B\>] ^_` abc d$ E\') A;a VUa'kQ%66;rABFGH  .."q&BrE)377CSWZWdWdIdEe fgh ijk l%t|) A;a VTQ]a$77"<BRGHI  rc|\}}}|j}tj||dz z tj||ttj zdd}tj||dz z tj||ttj zdd}||z }t|} gdgdgdd} ||| || zz} | d |zz} | |z dz|z} | |z dz|z} t| t| gS) NrrrEr)g??rr??r)rr0rr1rr/rY)rr1rr/rYrr0)rWgulika yamagandar/) r rrr'rrrrrC) r~roptionrrrr"r#r*weekdayr start_timeend_times rtrikalamr8s,#sB ~~" ..beSWWc3{S]]?Z [\] ^_` a% RU CGGS#kCLL>X YZ[ \]^ _$ E\' "I'DCF I'w!999* %'/ )(R2%**mr !B &(  fX. //rct||dS)NrWr8r~rs rrrs xE6:rct||dS)Nr3r:r;s rrrs HR $Drct||dS)Nr2r:r;s rrrs "eX!>rcX|\}}}|j}tj||dz z tj||ttj zdd}tj|dz|dz z tj||ttj zdd}||z }tj||dz z tj||ttj zdd}||z }t|} ddgddgd d gd dgd dgd dgd dgg} ||g} ||g} | dk(r || d<|| d<ddg} ddg}tddD]]}| | |}|dk7s| || || | |zdz z| |<| ||dzdz z||<| ||z dz|z| |<|||z dz|z||<_| |gS)NrrrErg$@rYg@g!@g333333@g333333@gffffff@g@g?rcrg?) r rrr'rrrrrG)r~rrrrr#r"r,r*r5rdurbase start_timesr+rLrs r durmuhurtamrBs,#sB ~~" RU CGGS#kCLL>X YZ[ \]^ _$ .."q&BrE)377C[SVS`S`E` abc def g%t|) ..beSWWc3{S]]?Z [\] ^_` a% E\' "I'C[3Z3Z3Z3Z3Z3Z  ' '# $ \YCF$QA+!f) A;a W a F }AwQ''*:1*=!=!BBk!n ^gmb&88il$A+r1B6k!nlR'2-2il  y !!rc|\}}}|j}tj||dz z tj||ttj zdd}tj||dz z tj||ttj zdd}||z }|d|zz}|d|zz} ||z dz|z| |z dz|zgS)zhAbhijit muhurta is the 8th muhurta (middle one) of the 15 muhurtas during the day_duration (~12 hours)rrrErg?g?)r rrr'rrr) r~rrrrr"r#r*r6r7s rabhijit_muhurtarD)s,#sB ~~" ..beSWWc3{S]]?Z [\] ^_` a% RU CGGS#kCLL>X YZ[ \]^ _$ E\'v''* Vg% %( ?b 2 %2 ';b'@ AArc >||jdz z }g}tD]}|tjk7r t ||}n#t t |tj }t|dz }t|dz}|j|||t|g|S)zComputes instantaneous planetary positions (i.e., which celestial object lies in which constellation) Also gives the nakshatra-pada division rr) r planet_listrr/rketuRAHUr;rCr)r)r~rjd_ut positionsr1 nirayana_long constellation coordinatess rplanetary_positionsrN:s u~~# #%)f (7m-eSXX>?m *+M+,K fm[.:WXY rc|\}}}||dz z }ttj|||ddd}t|dz }t |dz}t ||t |gS)z8Lagna (=ascendant) calculation at any given time & placerPrErr)rr houses_exr;rCrr) r~rrrrjd_utcnirayana_lagnarLrMs r ascendantrTQsx,#sB c?&==c48;A>.nr)*-~*++ n^&D EErc@d}d|z}||z }|dz}t|dzS)zgCalculates the navamsa-sign in which given longitude falls 0 = Aries, 1 = Taurus, ..., 11 = Pisces rrrEr)r rone_sign signs_elapsed fraction_lefts rnavamsa_from_longrYcs6( (](h&-!#- ]R   rc||jdz z }g}tD]a}|tjk7r t ||}n#t t |tj }|j|t|gc|S)z!Calculates navamsa of all planetsr) r rFrr/rrGr.r)rY)r~rrRrJr1rKs rnavamsar[mss $ $&)f (8m-fcmmDEm f/ >?@  rcttjjjtt t tttt tttt tdttt tdtt dk(sJttttdtt tdgk(sJy)NrEr)rsys _getframef_codeco_namerdate2 bangalorerrrrdate4shillongrhelsinkirrr all_testsrg}s &&' "#y!"y!!$%ui #$ u  x #$ x RD () (rcttjjjt t ddd}t t ddd}t t ddd}t t ddd}t t ddd }td tttttttttttttttttt|ttt|ttt|ttt|ttt|ttttttd y) Nrcrrz tithi startz tithi end)rr^r_r`ragregorian_to_jdrrdate1rcrbdate3rfshyambirthdatehonavar)feb3apr24apr19apr20apr21s r tithi_testsrys$ &&' dAq) *$ $tQ+ ,% $tQ+ ,% $tQ+ ,% $tQ+ ,% eY eY eY eXeY dIeXeXeXng&' rcNttjjjtt t ttt tttt tttt ttyr) rr^r_r`rarrprcrbrqrdrerrrnakshatra_testsr{s[ &&' %#$ %#$ %#$ %"#rc8ttjjjt t ddd}tttttttttt|ty)Nrir]) rr^r_r`rarorrrqrcrbrf)may22s r yoga_testsrs\ &&' $tQ+ ,%UIUIUHrcttjjjt t ddd}t t ddd}t t ddd}t t ddd}t t dd d }t t dd d }td tt|ttt|ttt|ttt|ttt|ttt|ty) Nrircrir(rr'r]rmrnshyam) rr^r_r`rarorrrcrf)r~aug17aug18sep19may20may21s r masa_testsrs &&'tD!R()" $tQ+ ,% $tQ+ ,% $tQ+ ,% $tQ+ ,% $tQ+ ,%.RUIUIUIUHUHrc6ttjjjt j dddd}t|tdgdddggk(sJt j ddd d }t|td gd d dggk(sJy)Nr'rg"""""7@rc)r%rr]r.?P*@r()rmrmrj) rr^r_r`rarrxrTrcr}s rascendant_testsrs &&' zz$2{+" 2y !aq!f%= => = zz$256" 2y !aAw%? ?@ ?rc  ttjjjt j dddd}t|t}ddgddgd dgd d gdd gd d gd d gd dgddgdd gdd gg }||k(sJy)Nrr'rrr rEr]rrcrjrr:rr() rr^r_r`rarrxr[rc)r~nvexpecteds r navamsa_testsrs &&' zz$256"r9""g1v1v1v1v2w!fr2hAB!R:( xr__main__gl)@gtfS@rgl9@gnV@g(\N@g(\8@g@g镲 ,@g̜R@irririr:r'rnrkrr(riir])rrrrY)rrV)n__doc__ __future__rmathr collectionsrstructswissephrrr rrBIT_NO_REFRACTIONrr0r/r'r(r)r*r+r,r-r.URANUSNEPTUNErF revati_359_50galc_cent_mid_mularrr3from_dmsrArCrMnorm180rrGrarmrvrojd_to_gregorianrrrrrrrrr hour_valrxr~r rise_timerrrrrrrrrrrr rrrrr r%r-r8 rahu_kalamyamaganda_kalam gulika_kalamrBrDrNrTrYr[rgryr{rrrr__name__r^rcrerfrsrprbrqrdapr_8apr_10rrdeepabirthdaterrrrsS , f./w=> !!C$9$99  99ww#((CKKyy#**cmmxxS[[3 N R=F ? A $' $"P9 "==(TeS SZZeS(+ B    (1 78 %il^ 45>'*&R F/d  % $ $ &0*; D> *"X B".F$!    *  A z FFD)) 664 (( 5&$ '( '5$ '' $tQ+ ,% $tQ+ ,% $tQ* +% $tQ+ ,% $tQ* +% 4a, -&"4b!#45."4a#34. ng&'ng&''r