public void CalculateEvents(GCLocation loc, GregorianDateTime vcStart, GregorianDateTime vcEnd)
{
//GCSunData sun = new GCSunData();
DstTypeChange ndst = 0;
int nData;
TResultEvents inEvents = this;
this.Clear();
this.EarthLocation = loc;
this.StartDateTime = new GregorianDateTime(vcStart);
this.EndDateTime = new GregorianDateTime(vcEnd);
GregorianDateTime vc = new GregorianDateTime();
GregorianDateTime vcAdd = new GregorianDateTime(), vcNext = new GregorianDateTime();
GCEarthData earth = loc.GetEarthData();
vc.Set(vcStart);
vc.TimezoneHours = loc.OffsetUtcHours;
vcAdd.Set(vc);
vcAdd.InitWeekDay();
GCHourTime sunriseData, sunsetData;
double sunRise, sunSet;
double r1, r2;
while (vcAdd.IsBeforeThis(vcEnd))
{
sunriseData = GCSunData.CalcSunrise(vcAdd, earth);
sunsetData = GCSunData.CalcSunset(vcAdd, earth);
sunRise = sunriseData.TotalDays;
sunSet = sunsetData.TotalDays;
ndst = loc.TimeZone.DetermineDaylightChange(vcAdd);
if (GCDisplaySettings.Current.getValue(GCDS.COREEVENTS_SUN) != 0)
{
ndst = loc.TimeZone.DetermineDaylightChange(vcAdd);
vcAdd.shour = sunriseData.TotalDays - 96.0 / 1440.0;
inEvents.AddEvent(vcAdd, CoreEventType.CCTYPE_S_ARUN, 0, ndst);
vcAdd.shour = sunRise;
inEvents.AddEvent(vcAdd, CoreEventType.CCTYPE_S_RISE, 0, ndst);
vcAdd.shour = (sunRise + sunSet) / 2;
inEvents.AddEvent(vcAdd, CoreEventType.CCTYPE_S_NOON, 0, ndst);
vcAdd.shour = sunSet;
inEvents.AddEvent(vcAdd, CoreEventType.CCTYPE_S_SET, 0, ndst);
}
/*if (GCDisplaySettings.Current.getValue(GCDS.COREEVENTS_ASCENDENT) != 0)
* {
* todayLongitude = sun.longitude_deg;
* vcAdd.shour = sunRise;
* todaySunriseHour = sunRise;
* if (previousLongitude < -10)
* {
* GregorianDateTime prevSunrise = new GregorianDateTime();
* prevSunrise.Set(vcAdd);
* prevSunrise.PreviousDay();
* sun.SunCalc(prevSunrise, earth);
* previousSunriseHour = sun.rise.GetDayTime() - 1;
* previousLongitude = sun.longitude_deg;
* fromTimeLimit = 0;
* }
*
* double a, b;
* double jd = vcAdd.GetJulianComplete();
* double ayan = GCAyanamsha.GetAyanamsa(jd);
* r1 = GCMath.putIn360(previousLongitude - ayan) / 30;
* r2 = GCMath.putIn360(todayLongitude - ayan) / 30;
*
* while (r2 > r1 + 13)
* {
* r2 -= 12.0;
* }
* while (r2 < r1 + 11)
* {
* r2 += 12.0;
* }
*
* a = (r2 - r1) / (todaySunriseHour - previousSunriseHour);
* b = r2 - a * todaySunriseHour;
*
* for (double tr = Math.Floor(r1) + 1.0; tr < r2; tr += 1.0)
* {
* double tm = (tr - b) / a;
* if (tm > fromTimeLimit)
* {
* vcNext.Set(vcAdd);
* vcNext.shour = tm;
* vcNext.NormalizeValues();
* inEvents.AddEvent(vcNext, CoreEventType.CCTYPE_ASCENDENT, (int)tr, ndst);
* }
* }
*
* previousLongitude = todayLongitude;
* previousSunriseHour = todaySunriseHour - 1;
* fromTimeLimit = previousSunriseHour;
* }*/
if (GCDisplaySettings.Current.getValue(GCDS.COREEVENTS_RAHUKALAM) != 0)
{
GCSunData.CalculateKala(sunRise, sunSet, vcAdd.dayOfWeek, out r1, out r2, KalaType.KT_RAHU_KALAM);
vcAdd.shour = r1;
inEvents.AddEvent(vcAdd, CoreEventType.CCTYPE_KALA_START, KalaType.KT_RAHU_KALAM, ndst);
vcAdd.shour = r2;
inEvents.AddEvent(vcAdd, CoreEventType.CCTYPE_KALA_END, KalaType.KT_RAHU_KALAM, ndst);
}
if (GCDisplaySettings.Current.getValue(GCDS.COREEVENTS_YAMAGHANTI) != 0)
{
GCSunData.CalculateKala(sunRise, sunSet, vcAdd.dayOfWeek, out r1, out r2, KalaType.KT_YAMA_GHANTI);
vcAdd.shour = r1;
inEvents.AddEvent(vcAdd, CoreEventType.CCTYPE_KALA_START, KalaType.KT_YAMA_GHANTI, ndst);
vcAdd.shour = r2;
inEvents.AddEvent(vcAdd, CoreEventType.CCTYPE_KALA_END, KalaType.KT_YAMA_GHANTI, ndst);
}
if (GCDisplaySettings.Current.getValue(GCDS.COREEVENTS_GULIKALAM) != 0)
{
GCSunData.CalculateKala(sunRise, sunSet, vcAdd.dayOfWeek, out r1, out r2, KalaType.KT_GULI_KALAM);
vcAdd.shour = r1;
inEvents.AddEvent(vcAdd, CoreEventType.CCTYPE_KALA_START, KalaType.KT_GULI_KALAM, ndst);
vcAdd.shour = r2;
inEvents.AddEvent(vcAdd, CoreEventType.CCTYPE_KALA_END, KalaType.KT_GULI_KALAM, ndst);
}
if (GCDisplaySettings.Current.getValue(GCDS.COREEVENTS_ABHIJIT_MUHURTA) != 0)
{
GCSunData.CalculateKala(sunRise, sunSet, vcAdd.dayOfWeek, out r1, out r2, KalaType.KT_ABHIJIT);
if (r1 > 0 && r2 > 0)
{
vcAdd.shour = r1;
inEvents.AddEvent(vcAdd, CoreEventType.CCTYPE_KALA_START, KalaType.KT_ABHIJIT, ndst);
vcAdd.shour = r2;
inEvents.AddEvent(vcAdd, CoreEventType.CCTYPE_KALA_END, KalaType.KT_ABHIJIT, ndst);
}
}
vcAdd.NextDay();
}
if (GCDisplaySettings.Current.getValue(GCDS.COREEVENTS_ASCENDENT) != 0)
{
GCAscendant asc = new GCAscendant();
asc.Earth = EarthLocation.GetEarthData();
asc.CurrentDateTime = new GregorianDateTime(vc);
while (asc.GetNextAscendantBefore(vcEnd))
{
ndst = loc.TimeZone.DetermineDaylightChange(vcNext);
inEvents.AddEvent(asc.CurrentDateTime, CoreEventType.CCTYPE_ASCENDENT, asc.CurrentSign, ndst);
asc.CurrentDateTime.AddHours(0.5);
}
}
if (GCDisplaySettings.Current.getValue(GCDS.COREEVENTS_TITHI) != 0)
{
vcAdd.Set(vc);
vcAdd.shour = 0.0;
while (vcAdd.IsBeforeThis(vcEnd))
{
nData = GCTithi.GetNextTithiStart(earth, vcAdd, out vcNext);
if (vcNext.GetDayInteger() < vcEnd.GetDayInteger())
{
vcNext.InitWeekDay();
ndst = loc.TimeZone.DetermineDaylightChange(vcNext);
inEvents.AddEvent(vcNext, CoreEventType.CCTYPE_TITHI, nData, ndst);
}
else
{
break;
}
vcAdd.Set(vcNext);
vcAdd.shour += 0.2;
if (vcAdd.shour >= 1.0)
{
vcAdd.shour -= 1.0;
vcAdd.NextDay();
}
}
}
if (GCDisplaySettings.Current.getValue(GCDS.COREEVENTS_NAKSATRA) != 0)
{
vcAdd.Set(vc);
vcAdd.shour = 0.0;
while (vcAdd.IsBeforeThis(vcEnd))
{
nData = GCNaksatra.GetNextNaksatra(earth, vcAdd, out vcNext);
if (vcNext.GetDayInteger() < vcEnd.GetDayInteger())
{
vcNext.InitWeekDay();
ndst = loc.TimeZone.DetermineDaylightChange(vcNext);
inEvents.AddEvent(vcNext, CoreEventType.CCTYPE_NAKS, nData, ndst);
}
else
{
break;
}
vcAdd.Set(vcNext);
vcAdd.shour += 0.2;
if (vcAdd.shour >= 1.0)
{
vcAdd.shour -= 1.0;
vcAdd.NextDay();
}
}
}
if (GCDisplaySettings.Current.getValue(GCDS.COREEVENTS_YOGA) != 0)
{
vcAdd.Set(vc);
vcAdd.shour = 0.0;
while (vcAdd.IsBeforeThis(vcEnd))
{
nData = GCYoga.GetNextYogaStart(earth, vcAdd, out vcNext);
if (vcNext.GetDayInteger() < vcEnd.GetDayInteger())
{
vcNext.InitWeekDay();
ndst = loc.TimeZone.DetermineDaylightChange(vcNext);
inEvents.AddEvent(vcNext, CoreEventType.CCTYPE_YOGA, nData, ndst);
}
else
{
break;
}
vcAdd.Set(vcNext);
vcAdd.shour += 0.2;
if (vcAdd.shour >= 1.0)
{
vcAdd.shour -= 1.0;
vcAdd.NextDay();
}
}
}
if (GCDisplaySettings.Current.getValue(GCDS.COREEVENTS_SANKRANTI) != 0)
{
vcNext = new GregorianDateTime();
vcAdd.Set(vc);
vcAdd.shour = 0.0;
while (vcAdd.IsBeforeThis(vcEnd))
{
vcNext.Set(GCSankranti.GetNextSankranti(vcAdd, earth, out nData));
if (vcNext.GetDayInteger() < vcEnd.GetDayInteger())
{
vcNext.InitWeekDay();
ndst = loc.TimeZone.DetermineDaylightChange(vcNext);
inEvents.AddEvent(vcNext, CoreEventType.CCTYPE_SANK, nData, ndst);
}
else
{
break;
}
vcAdd.Set(vcNext);
vcAdd.NextDay();
}
}
if (GCDisplaySettings.Current.getValue(GCDS.COREEVENTS_MOONRASI) != 0)
{
vcAdd.Set(vc);
vcAdd.shour = 0.0;
while (vcAdd.IsBeforeThis(vcEnd))
{
nData = GCMoonData.GetNextMoonRasi(earth, vcAdd, out vcNext);
if (vcNext.GetDayInteger() < vcEnd.GetDayInteger())
{
vcNext.InitWeekDay();
ndst = loc.TimeZone.DetermineDaylightChange(vcNext);
inEvents.AddEvent(vcNext, CoreEventType.CCTYPE_M_RASI, nData, ndst);
}
else
{
break;
}
vcAdd.Set(vcNext);
vcAdd.shour += 0.5;
vcAdd.NormalizeValues();
}
}
if (GCDisplaySettings.Current.getValue(GCDS.COREEVENTS_CONJUNCTION) != 0)
{
double dlong;
vcAdd.Set(vc);
vcAdd.shour = 0.0;
while (vcAdd.IsBeforeThis(vcEnd))
{
dlong = GCConjunction.GetNextConjunction(vcAdd, out vcNext, true, earth);
if (vcNext.GetDayInteger() < vcEnd.GetDayInteger())
{
vcNext.InitWeekDay();
ndst = loc.TimeZone.DetermineDaylightChange(vcNext);
inEvents.AddEvent(vcNext, CoreEventType.CCTYPE_CONJ, GCRasi.GetRasi(dlong, GCAyanamsha.GetAyanamsa(vcNext.GetJulianComplete())), ndst);
}
else
{
break;
}
vcAdd.Set(vcNext);
vcAdd.NextDay();
}
}
if (GCDisplaySettings.Current.getValue(GCDS.COREEVENTS_MOON) != 0)
{
vcAdd.Set(vc);
vcAdd.shour = 0.0;
while (vcAdd.IsBeforeThis(vcEnd))
{
vcNext.Set(GCMoonData.GetNextRise(earth, vcAdd, true));
inEvents.AddEvent(vcNext, CoreEventType.CCTYPE_M_RISE, 0, ndst);
vcNext.Set(GCMoonData.GetNextRise(earth, vcNext, false));
inEvents.AddEvent(vcNext, CoreEventType.CCTYPE_M_SET, 0, ndst);
vcNext.shour += 0.05;
vcNext.NormalizeValues();
vcAdd.Set(vcNext);
}
}
if (GCDisplaySettings.Current.getValue(GCDS.COREEVENTS_ASCENDENT) != 0)
{/*
* vcAdd = vc;
* vcAdd.shour = 0.0;
* while(vcAdd.IsBeforeThis(vcEnd))
* {
* nData = earth.GetNextAscendentStart(vcAdd, vcNext);
* if (vcNext.GetDayInteger() < vcEnd.GetDayInteger())
* {
* vcNext.InitWeekDay();
* ndst = TTimeZone.determineDaylightChange(vcNext, loc.m_nDST);
* inEvents.AddEvent(vcNext, CCTYPE_ASCENDENT, nData, ndst);
* }
* else
* {
* break;
* }
* vcAdd = vcNext;
* vcAdd.shour += 1/24.0;
* if (vcAdd.shour >= 1.0)
* {
* vcAdd.shour -= 1.0;
* vcAdd.NextDay();
* }
* }
*
*/
}
if (GCDisplaySettings.Current.getValue(GCDS.COREEVENTS_SORT) != 0)
{
inEvents.Sort(SORTING_BY_DATE);
}
else
{
inEvents.Sort(SORTING_BY_TYPE);
}
}
public void CalculateEvents(GCLocation loc, int nYear)
{
//GCSunData sun = new GCSunData();
//DstTypeChange ndst = DstTypeChange.DstOff;
int nData;
TResultCoreEvents inEvents = this;
this.Clear();
this.EarthLocation = loc;
this.Year = nYear;
GregorianDateTime vcStart = new GregorianDateTime(Year - 1, 12, 29);
GregorianDateTime vcEnd = new GregorianDateTime(Year + 1, 1, 2);
GregorianDateTime vc = new GregorianDateTime();
GregorianDateTime vcAdd = new GregorianDateTime(), vcNext = new GregorianDateTime();
GCEarthData earth = loc.GetEarthData();
vc.Set(vcStart);
vc.TimezoneHours = loc.OffsetUtcHours;
vcAdd.Set(vc);
vcAdd.InitWeekDay();
/*while (vcAdd.IsBeforeThis(vcEnd))
* {
* ndst = loc.TimeZone.DetermineDaylightChange(vcAdd);
* vcAdd.NextDay();
* }*/
if (Full || GCDisplaySettings.Current.getValue(GCDS.COREEVENTS_TITHI) != 0)
{
vcAdd.Set(vc);
vcAdd.shour = 0.0;
while (vcAdd.IsBeforeThis(vcEnd))
{
nData = GCTithi.GetNextTithiStart(earth, vcAdd, out vcNext);
if (vcNext.GetDayInteger() < vcEnd.GetDayInteger())
{
//vcNext.InitWeekDay();
//ndst = loc.TimeZone.DetermineDaylightChange(vcNext);
inEvents.AddEvent(vcNext, CoreEventType.CCTYPE_TITHI, nData);
}
else
{
break;
}
vcAdd.Set(vcNext);
vcAdd.shour += 0.2;
if (vcAdd.shour >= 1.0)
{
vcAdd.shour -= 1.0;
vcAdd.NextDay();
}
}
}
if (Full || GCDisplaySettings.Current.getValue(GCDS.COREEVENTS_NAKSATRA) != 0)
{
vcAdd.Set(vc);
vcAdd.shour = 0.0;
while (vcAdd.IsBeforeThis(vcEnd))
{
nData = GCNaksatra.GetNextNaksatra(earth, vcAdd, out vcNext);
if (vcNext.GetDayInteger() < vcEnd.GetDayInteger())
{
//vcNext.InitWeekDay();
//ndst = loc.TimeZone.DetermineDaylightChange(vcNext);
inEvents.AddEvent(vcNext, CoreEventType.CCTYPE_NAKS, nData);
}
else
{
break;
}
vcAdd.Set(vcNext);
vcAdd.shour += 0.2;
if (vcAdd.shour >= 1.0)
{
vcAdd.shour -= 1.0;
vcAdd.NextDay();
}
}
}
if (Full || GCDisplaySettings.Current.getValue(GCDS.COREEVENTS_YOGA) != 0)
{
vcAdd.Set(vc);
vcAdd.shour = 0.0;
while (vcAdd.IsBeforeThis(vcEnd))
{
nData = GCYoga.GetNextYogaStart(earth, vcAdd, out vcNext);
if (vcNext.GetDayInteger() < vcEnd.GetDayInteger())
{
//vcNext.InitWeekDay();
//ndst = loc.TimeZone.DetermineDaylightChange(vcNext);
inEvents.AddEvent(vcNext, CoreEventType.CCTYPE_YOGA, nData);
}
else
{
break;
}
vcAdd.Set(vcNext);
vcAdd.shour += 0.2;
if (vcAdd.shour >= 1.0)
{
vcAdd.shour -= 1.0;
vcAdd.NextDay();
}
}
}
if (Full || GCDisplaySettings.Current.getValue(GCDS.COREEVENTS_SANKRANTI) != 0)
{
vcNext = new GregorianDateTime();
vcAdd.Set(vc);
vcAdd.shour = 0.0;
while (vcAdd.IsBeforeThis(vcEnd))
{
vcNext.Set(GCSankranti.GetNextSankranti(vcAdd, earth, out nData));
if (vcNext.GetDayInteger() < vcEnd.GetDayInteger())
{
//vcNext.InitWeekDay();
//ndst = loc.TimeZone.DetermineDaylightChange(vcNext);
inEvents.AddEvent(vcNext, CoreEventType.CCTYPE_SANK, nData);
}
else
{
break;
}
vcAdd.Set(vcNext);
vcAdd.NextDay();
}
}
if (Full || GCDisplaySettings.Current.getValue(GCDS.COREEVENTS_MOONRASI) != 0)
{
vcAdd.Set(vc);
vcAdd.shour = 0.0;
while (vcAdd.IsBeforeThis(vcEnd))
{
nData = GCMoonData.GetNextMoonRasi(earth, vcAdd, out vcNext);
if (vcNext.GetDayInteger() < vcEnd.GetDayInteger())
{
//vcNext.InitWeekDay();
//ndst = loc.TimeZone.DetermineDaylightChange(vcNext);
inEvents.AddEvent(vcNext, CoreEventType.CCTYPE_M_RASI, nData);
}
else
{
break;
}
vcAdd.Set(vcNext);
vcAdd.shour += 0.5;
vcAdd.NormalizeValues();
}
}
if (Full || GCDisplaySettings.Current.getValue(GCDS.COREEVENTS_CONJUNCTION) != 0)
{
double dlong;
vcAdd.Set(vc);
vcAdd.shour = 0.0;
while (vcAdd.IsBeforeThis(vcEnd))
{
dlong = GCConjunction.GetNextConjunction(vcAdd, out vcNext, true, earth);
if (vcNext.GetDayInteger() < vcEnd.GetDayInteger())
{
vcNext.InitWeekDay();
//ndst = loc.TimeZone.DetermineDaylightChange(vcNext);
inEvents.AddEvent(vcNext, CoreEventType.CCTYPE_CONJ, GCRasi.GetRasi(dlong, GCAyanamsha.GetAyanamsa(vcNext.GetJulianComplete())));
}
else
{
break;
}
vcAdd.Set(vcNext);
vcAdd.NextDay();
}
}
inEvents.Sort();
}
/*********************************************************************/
/* */
/* */
/* */
/* */
/* */
/*********************************************************************/
public int MasaCalc(GregorianDateTime date, GCEarthData earth)
{
// SUNDATA sun;
// MOONDATA moon;
// VCTIME Conj[6], ConjA;
// double Long[6], LongA;
// int Sank[6], SankA;
const int PREV_MONTHS = 6;
double[] L = new double[8];
GregorianDateTime[] C = new GregorianDateTime[8];
int[] R = new int[8];
int n, rasi, masa = MasaId.VISNU_MASA;
int ksaya_from = -1;
int ksaya_to = -1;
date.shour = sunRise.TotalDays;
// STEP 1: calculate position of the sun and moon
// it is done by previous call of DayCalc
// and results are in argument DAYDATA day
// *DayCalc(date, earth, day, moon, sun);*
L[1] = /*Long[0] =*/ GCConjunction.GetNextConjunction(date, out C[1], false, earth);
L[0] = /*LongA =*/ GCConjunction.GetNextConjunction(C[1], out C[0], true, earth);
// on Pratipat (nTithi == 15) we need to look for previous conjunction
// but this conjunction can occur on this date before sunrise
// so we need to include this very date into looking for conjunction
// on other days we cannot include it
// and exclude it from looking for next because otherwise that will cause
// incorrect detection of Purusottama adhika masa
L[2] = GCConjunction.GetPrevConjunction(date, out C[2], false, earth);
for (n = 3; n < PREV_MONTHS; n++)
{
L[n] = GCConjunction.GetPrevConjunction(C[n - 1], out C[n], true, earth);
}
for (n = 0; n < PREV_MONTHS; n++)
{
R[n] = GCRasi.GetRasi(L[n], GCAyanamsha.GetAyanamsa(C[n].GetJulian()));
}
/* TRACE("TEST Date: %d %d %d ", date.day, date.month, date.year);
* TRACE("FOUND CONJ Date: %d %d %d rasi: %d ", C[1].day, C[1].month, C[1].year, R[1]);
* TRACE("FOUND CONJ Date: %d %d %d rasi: %d ", C[2].day, C[2].month, C[2].year, R[2]);
* TRACE("FOUND CONJ Date: %d %d %d rasi: %d ", C[3].day, C[3].month, C[3].year, R[3]);
* TRACE("FOUND CONJ Date: %d %d %d rasi: %d ", C[4].day, C[4].month, C[4].year, R[4]);
* TRACE("---
*/
// test for Adhika-Ksaya sequence
// this is like 1-2-2-4-5...
// second (2) is replaced by rasi(3)
/* if ( ((Sank[1] + 2) % 12 == SankA) && ((Sank[1] == Sank[0]) || (Sank[0] == SankA)))
* {
* Sank[0] = (Sank[1] + 1) % 12;
* }
*
* if ( ((Sank[2] + 2) % 12 == Sank[0]) && ((Sank[2] == Sank[1]) || (Sank[1] == Sank[0])))
* {
* Sank[1] = (Sank[2] + 1) % 12;
* }*/
// look for ksaya month
ksaya_from = -1;
for (n = PREV_MONTHS - 2; n >= 0; n--)
{
if ((R[n + 1] + 2) % 12 == R[n])
{
ksaya_from = n;
break;
}
}
if (ksaya_from >= 0)
{
for (n = ksaya_from; n > 0; n--)
{
if (R[n] == R[n - 1])
{
ksaya_to = n;
break;
}
}
if (ksaya_to >= 0)
{
// adhika masa found
// now correct succession of rasis
}
else
{
// adhika masa not found
// there will be some break in masa queue
ksaya_to = 0;
}
int current_rasi = R[ksaya_from + 1] + 1;
for (n = ksaya_from; n >= ksaya_to; n--)
{
R[n] = current_rasi;
current_rasi = (current_rasi + 1) % 12;
}
}
// STEP 3: test for adhika masa
// test for adhika masa
if (R[1] == R[2])
{
// it is adhika masa
masa = 12;
rasi = R[1];
}
else
{
// STEP 2. select nearest Conjunction
if (this.sunRise.Paksa == 0)
{
masa = R[1];
}
else if (this.sunRise.Paksa == 1)
{
masa = R[2];
}
rasi = masa;
}
// calculation of Gaurabda year
GaurabdaYear = date.year - 1486;
if ((rasi > 7) && (rasi < 11)) // Visnu
{
if (date.month < 6)
{
GaurabdaYear--;
}
}
Masa = masa;
return(masa);
}
