def set_jupiter_transits(self):
if 'guru-saGkrAntiH' not in self.rules_collection.name_to_rule:
return
jd_end = self.panchaanga.jd_start + self.panchaanga.duration + 13
check_window = 400 # Max t between two Jupiter transits is ~396 (checked across 180y)
# Let's check for transitions in a relatively large window
# to finalise what is the FINAL transition post retrograde movements
transits = Graha.singleton(Graha.JUPITER).get_transits(self.panchaanga.jd_start - 13, jd_end + check_window, anga_type=AngaType.RASHI,
ayanaamsha_id=self.ayanaamsha_id)
if len(transits) > 0:
for i, transit in enumerate(transits):
(jd_transit, rashi1, rashi2) = (transit.jd, transit.value_1, transit.value_2)
if self.panchaanga.jd_start - 13 < jd_transit < jd_end:
fday = int(floor(jd_transit) - floor(self.daily_panchaangas[0].julian_day_start))
fest = TransitionFestivalInstance(name='guru-saGkrAntiH', status_1_hk=names.NAMES['RASHI_NAMES']['sa']['hk'][rashi1], status_2_hk=names.NAMES['RASHI_NAMES']['sa']['hk'][rashi2])
self.panchaanga.add_festival_instance(festival_instance=fest, date=self.daily_panchaangas[fday].date)
if (rashi1 % 12 + 1) == rashi2 and ((transits[i + 1].value_1 % 12) + 1) == transits[i + 1].value_2:
# Considering only non-retrograde transits for pushkara computations
# logging.debug('Non-retrograde transit; we have a pushkaram!')
(madhyanha_start, madhyaahna_end) = interval.get_interval(self.daily_panchaangas[fday].jd_sunrise,
self.daily_panchaangas[fday].jd_sunset, 2, 5).to_tuple()
if jd_transit < madhyaahna_end:
fday_pushkara = fday
else:
fday_pushkara = fday + 1
self.panchaanga.add_festival(
fest_id='%s-Adi-puSkara-ArambhaH' % names.NAMES['PUSHKARA_NAMES']['sa']['hk'][rashi2], date=self.daily_panchaangas[fday_pushkara].date)
self.panchaanga.add_festival(
fest_id='%s-Adi-puSkara-samApanam' % names.NAMES['PUSHKARA_NAMES']['sa']['hk'][rashi2], date=self.daily_panchaangas[fday_pushkara].date + 11)
self.panchaanga.add_festival(
fest_id='%s-antya-puSkara-samApanam' % names.NAMES['PUSHKARA_NAMES']['sa']['hk'][rashi1], date=self.daily_panchaangas[fday_pushkara].date - 1)
self.panchaanga.add_festival(
fest_id='%s-antya-puSkara-ArambhaH' % names.NAMES['PUSHKARA_NAMES']['sa']['hk'][rashi1], date=self.daily_panchaangas[fday_pushkara].date - 12)
def compute_tb_muhuurtas(self):
""" Computes muhuurta-s according to taittiriiya brAhmaNa.
"""
if getattr(self, "jd_sunrise", None) is None:
self.compute_sun_moon_transitions()
tb_muhuurtas = []
for muhuurta_id in range(0, 15):
(jd_start,
jd_end) = interval.get_interval(start_jd=self.jd_sunrise,
end_jd=self.jd_sunset,
part_index=muhuurta_id,
num_parts=15).to_tuple()
from jyotisha.panchaanga.temporal.interval import TbSayanaMuhuurta
tb_muhuurtas.append(
TbSayanaMuhuurta(jd_start=jd_start,
jd_end=jd_end,
muhuurta_id=muhuurta_id))
self.day_length_based_periods.tb_muhuurtas = tb_muhuurtas
def assign_shraaddha_tithi(self, debug_shraaddha_tithi=False):
tithi_days = [{z: [] for z in range(0, 32)} for _x in range(13)]
lunar_tithi_days = {}
daily_panchaangas = self.daily_panchaangas
def _assign(self, fday, tithi):
if daily_panchaangas[
fday].shraaddha_tithi == [] or daily_panchaangas[
fday].shraaddha_tithi == [tithi]:
daily_panchaangas[fday].shraaddha_tithi = [tithi]
else:
daily_panchaangas[fday].shraaddha_tithi.append(tithi)
if daily_panchaangas[fday -
1].shraaddha_tithi.count(tithi) == 1:
daily_panchaangas[fday - 1].shraaddha_tithi.remove(tithi)
for d in range(self.panchaanga.duration_prior_padding,
self.panchaanga.duration + 1):
[y, m, dt,
t] = time.jd_to_utc_gregorian(self.panchaanga.jd_start + d -
1).to_date_fractional_hour_tuple()
(d0_angas, d1_angas) = get_2_day_interval_boundary_angas(
kaala="aparaahna",
anga_type=AngaType.TITHI,
p0=self.daily_panchaangas[d],
p1=self.daily_panchaangas[d + 1])
angam_start = d0_angas.start
next_anga = angam_start + 1
nnext_anga = next_anga + 1
# Calc vyaaptis
t_start_d, t_end_d = interval.get_interval(
daily_panchaangas[d].jd_sunrise,
daily_panchaangas[d].jd_sunset, 3, 5).to_tuple()
span_1 = t_end_d - t_start_d
span_2 = 0
for tithi_span in daily_panchaangas[
d].sunrise_day_angas.tithis_with_ends:
tithi_end = tithi_span.jd_end
if tithi_end is None:
pass
elif t_start_d < tithi_end < t_end_d:
span_1 = tithi_end - t_start_d
span_2 = t_end_d - tithi_end
t_start_d1, t_end_d1 = interval.get_interval(
daily_panchaangas[d + 1].jd_sunrise,
daily_panchaangas[d + 1].jd_sunset, 3, 5).to_tuple()
vyapti_3 = t_end_d1 - t_start_d1
for tithi_span in daily_panchaangas[
d + 1].sunrise_day_angas.tithis_with_ends:
tithi_end = tithi_span.jd_end
if tithi_end is None:
pass
elif t_start_d1 < tithi_end < t_end_d1:
vyapti_3 = tithi_end - t_start_d1
# Combinations
# <a> 1 1 1 1 - d + 1: 1
# <b> 1 1 2 2 - d: 1
# <f> 1 1 2 3 - d: 1, d+1: 2
# <e> 1 1 1 2 - d, or vyApti (just in case next day aparaahna is slightly longer): 1
# <d> 1 1 3 3 - d: 1, 2
# <h> 1 2 3 3 - d: 2
# <c> 1 2 2 2 - d + 1: 2
# <g> 1 2 2 3 - vyApti: 2
fday = -1
reason = '?'
# if sunrise_day_angas[1] == angam_start:
# logging.debug('Pre-emptively assign %2d to %3d, can be removed tomorrow if need be.' % (angam_start, d))
# _assign(self, d, angam_start)
if d1_angas.end == angam_start: # <a>
# Full aparaahnas on both days, so second day
fday = d + 1
s_tithi = angam_start
reason = '%2d incident on consecutive days; paraviddhA' % s_tithi.index
elif (d0_angas.end == angam_start) and (d1_angas.start
== next_anga): # <b>/<f>
fday = d
s_tithi = d0_angas.start
reason = '%2d not incident on %3d' % (s_tithi.index, d + 1)
if d1_angas.end == nnext_anga: # <f>
if debug_shraaddha_tithi:
logging.debug('%03d [%4d-%02d-%02d]: %s' % (
d, y, m, dt,
'Need to assign %2d to %3d as it is present only at start of aparAhna tomorrow!)'
% (next_anga.index, d + 1)))
_assign(self, d + 1, next_anga)
elif d1_angas.start == angam_start: # <e>
if span_1 > vyapti_3:
# Most likely
fday = d
s_tithi = d1_angas.start
reason = '%2d has more vyApti on day %3d (%f ghatikAs; full?) compared to day %3d (%f ghatikAs)' % (
s_tithi.index, d, span_1 * 60, d + 1, vyapti_3 * 60)
else:
fday = d + 1
s_tithi = d1_angas.start
reason = '%2d has more vyApti on day %3d (%f ghatikAs) compared to day %3d (%f ghatikAs) --- unusual!' % (
s_tithi.index, d + 1, vyapti_3 * 60, d, span_1 * 60)
elif d1_angas.start == nnext_anga: # <d>/<h>
if d0_angas.end == next_anga: # <h>
fday = d
s_tithi = d0_angas.end
reason = '%2d has some vyApti on day %3d; not incident on day %3d at all' % (
s_tithi.index, d, d + 1)
else: # <d>
s_tithi = angam_start
fday = d
reason = '%2d is incident fully at aparAhna today (%3d), and not incident tomorrow (%3d)!' % (
s_tithi.index, d, d + 1)
# Need to check vyApti of next_anga in sAyaMkAla: if it's nearly entire sAyaMkAla ie 5-59-30 or more!
if debug_shraaddha_tithi:
logging.debug('%03d [%4d-%02d-%02d]: %s' % (
d, y, m, dt,
'%2d not incident at aparAhna on either day (%3d/%3d); picking second day %3d!'
% (next_anga, d, d + 1, d + 1)))
_assign(self, d + 1, next_anga)
# logging.debug(reason)
elif d0_angas.end == d1_angas.start == d1_angas.end == next_anga: # <c>
s_tithi = next_anga
fday = d + 1
reason = '%2d has more vyApti on %3d (full) compared to %3d (part)' % (
s_tithi.index, d + 1, d)
elif d0_angas.end == d1_angas.start == next_anga: # <g>
s_tithi = d1_angas.start
if span_2 > vyapti_3:
# Most likely
fday = d
reason = '%2d has more vyApti on day %3d (%f ghatikAs) compared to day %3d (%f ghatikAs)' % (
s_tithi.index, d, span_2 * 60, d + 1, vyapti_3 * 60)
else:
fday = d + 1
reason = '%2d has more vyApti on day %3d (%f ghatikAs) compared to day %3d (%f ghatikAs)' % (
s_tithi.index, d + 1, vyapti_3 * 60, d, span_2 * 60
) # Examine for greater vyApti
else:
logging.error('Should not reach here ever! %s %s',
d0_angas.to_tuple(), d1_angas.to_tuple())
reason = '?'
if debug_shraaddha_tithi:
logging.debug(
'%03d [%4d-%02d-%02d]: Assigning tithi %2d to %3d (%s).' %
(d, y, m, dt, s_tithi.index, fday, reason))
_assign(self, fday, s_tithi)
if debug_shraaddha_tithi:
logging.debug(self.panchaanga.shraaddha_tithi)
for z in set([x.lunar_month_sunrise for x in self.daily_panchaangas]):
lunar_tithi_days[z] = {}
for d in range(self.panchaanga.duration_prior_padding,
self.panchaanga.duration + 1):
for t in daily_panchaangas[d].shraaddha_tithi:
lunar_tithi_days[daily_panchaangas[d].lunar_month_sunrise][
t.index] = d
# Following this primary assignment, we must now "clean" for Sankranti, and repetitions
# If there are two tithis, take second. However, if the second has sankrAnti dushtam, take
# first. If both have sankrAnti dushtam, take second.
for d in range(self.panchaanga.duration_prior_padding,
self.panchaanga.duration + 1):
if daily_panchaangas[d].shraaddha_tithi != []:
if daily_panchaangas[
d].solar_sidereal_date_sunset.month_transition is not None:
if debug_shraaddha_tithi:
logging.debug(
(d, daily_panchaangas[d].
solar_sidereal_date_sunset.month_transition))
aparaahna_start, aparaahna_end = interval.get_interval(
daily_panchaangas[d].jd_sunrise,
daily_panchaangas[d].jd_sunset, 3, 5).to_tuple()
m1 = daily_panchaangas[
d -
1].solar_sidereal_date_sunset.month # Previous month
m2 = daily_panchaangas[
d].solar_sidereal_date_sunset.month # Current month
if aparaahna_start < daily_panchaangas[
d].solar_sidereal_date_sunset.month_transition < aparaahna_end:
if debug_shraaddha_tithi:
logging.debug(
'Sankranti in aparaahna! Assigning to both months!'
)
assert daily_panchaangas[
d].solar_sidereal_date_sunset.day == 1
for t in daily_panchaangas[d].shraaddha_tithi:
# Assigning to both months --- should get eliminated because of a second occurrence
tithi_days[m1][t.index].extend([d, '*'])
tithi_days[m2][t.index].extend([d, '*'])
if daily_panchaangas[
d].solar_sidereal_date_sunset.month_transition < aparaahna_start:
if debug_shraaddha_tithi:
logging.debug('Sankranti before aparaahna!')
assert daily_panchaangas[
d].solar_sidereal_date_sunset.day == 1
for t in daily_panchaangas[d].shraaddha_tithi:
tithi_days[m2][t.index].extend([d, '*'])
if aparaahna_end < daily_panchaangas[
d].solar_sidereal_date_sunset.month_transition:
if debug_shraaddha_tithi:
logging.debug('Sankranti after aparaahna!')
# Depending on whether sankranti is before or after sunset, m2 may or may not be equal to m1
# In any case, we wish to assign this tithi to the previous month, where it really occurs.
for t in daily_panchaangas[d].shraaddha_tithi:
tithi_days[m1][t.index].extend([d, '*'])
else:
for t in daily_panchaangas[d].shraaddha_tithi:
tithi_days[daily_panchaangas[d].
solar_sidereal_date_sunset.month][
t.index].append(d)
# We have now assigned all tithis. Now, remove duplicates based on the above-mentioned rules.
# TODO: This is not the best way to clean. Need to examine one month at a time.
for m in range(1, 13):
for t in range(1, 31):
if len(tithi_days[m][t]) == 1:
continue
elif len(tithi_days[m][t]) == 2:
if tithi_days[m][t][1] == '*':
# Only one tithi available!
if debug_shraaddha_tithi:
logging.debug(
'Only one %2d tithi in month %2d, on day %3d, despite sankrAnti dushtam!'
% (t, m, tithi_days[m][t][0]))
del tithi_days[m][t][1]
tithi_days[m][t][0] = '%d::%d' % (tithi_days[m][t][0],
m)
if debug_shraaddha_tithi:
logging.debug('Note %s' % str(tithi_days[m][t]))
else:
daily_panchaangas[
tithi_days[m][t][0]].shraaddha_tithi = [
Anga.get_cached(
index=0, anga_type_id=AngaType.TITHI.name)
] # Shunya
if debug_shraaddha_tithi:
logging.debug('Removed %d' % tithi_days[m][t][0])
del tithi_days[m][t][0]
if debug_shraaddha_tithi:
logging.debug(
'Two %2d tithis in month %2d: retaining second on %2d!'
% (t, m, tithi_days[m][t][0]))
elif len(tithi_days[m][t]) == 3:
if debug_shraaddha_tithi:
logging.debug('Two %2d tithis in month %2d: %s' %
(t, m, str(tithi_days[m][t])))
if tithi_days[m][t][1] == '*':
daily_panchaangas[
tithi_days[m][t][0]].shraaddha_tithi = [
Anga.get_cached(
index=0, anga_type_id=AngaType.TITHI.name)
] # Shunya
if debug_shraaddha_tithi:
logging.debug('Removed %d' % tithi_days[m][t][0])
del tithi_days[m][t][:2]
elif tithi_days[m][t][2] == '*':
daily_panchaangas[
tithi_days[m][t][1]].shraaddha_tithi = [
Anga.get_cached(
index=0, anga_type_id=AngaType.TITHI.name)
] # Shunya
if debug_shraaddha_tithi:
logging.debug('Removed %d' % tithi_days[m][t][1])
del tithi_days[m][t][1:]
if debug_shraaddha_tithi:
logging.debug(' Retaining non-dushta: %s' %
(str(tithi_days[m][t])))
elif len(tithi_days[m][t]) == 4:
if debug_shraaddha_tithi:
logging.debug(
'Two dushta %2d tithis in month %2d: %s' %
(t, m, str(tithi_days[m][t])))
daily_panchaangas[tithi_days[m][t][0]].shraaddha_tithi = [
Anga.get_cached(index=0,
anga_type_id=AngaType.TITHI.name)
] # Shunya
if debug_shraaddha_tithi:
logging.debug('Removed %d' % tithi_days[m][t][0])
tithi_days[m][t][3] = str(m)
del tithi_days[m][t][:2]
if debug_shraaddha_tithi:
logging.debug(' Retaining: %s' %
(str(tithi_days[m][t])))
tithi_days[m][t][0] = '%d::%d' % (tithi_days[m][t][0], m)
if debug_shraaddha_tithi:
logging.debug('Note %s' % str(tithi_days[m][t]))
elif len(tithi_days[m][t]) == 0:
logging.warning(
'Rare issue. No tithi %d in this solar month (%d). Therefore use lunar tithi.'
% (t, m))
# सौरमासे तिथ्यलाभे चान्द्रमानेन कारयेत्
# tithi_days[m][t] = lunar_tithi_days[m][t]
else:
logging.error(
'Something weird. len(tithi_days[m][t]) is not in 1:4!! : %s (m=%d, t=%d)',
str(tithi_days[m][t]), m, t)
if debug_shraaddha_tithi:
logging.debug(tithi_days)
def test_get_interval():
assert interval.get_interval(start_jd=1,
end_jd=11,
part_index=[0],
num_parts=1) == Interval(jd_start=1,
jd_end=11)
