/*********************************************************************/
/* Finds date of Pratipat, Krsna Paksa, Visnu Masa */
/* */
/* EARTHDATA earth - location */
/* int nYear - Gregorian year */
/* */
/*********************************************************************/
public static GregorianDateTime GetFirstDayOfYear(GCEarthData earth, int nYear)
{
int[] a = { 2, 15, 3, 1, 3, 15, 4, 1, 4, 15 };
GregorianDateTime d = new GregorianDateTime();
GCAstroData day = new GCAstroData();
int gy, j, masa;
int step;
UInt32 tmp;
if (nYear >= 1950 && nYear < 2058)
{
tmp = gGaurBeg[(nYear - 1950) * 26 + 22];
d.month = (int)(tmp & 0x3e0) >> 5;
d.day = (int)(tmp & 0x1f);
d.year = nYear;
d.NextDay();
a[0] = d.month;
a[1] = d.day;
}
for (int i = 0; i < 10; i += 2)
{
d.year = nYear;
d.month = a[i];
d.day = a[i + 1];
day.DayCalc(d, earth);
masa = day.MasaCalc(d, earth);
gy = day.GaurabdaYear;
if (masa == 11) // visnu masa
{
do
{
// shifts date
step = day.sunRise.Tithi / 2;
step = (step > 0) ? step : 1;
for (j = step; j > 0; j--)
{
d.PreviousDay();
}
// try new time
day.DayCalc(d, earth);
}while (day.sunRise.Tithi < 28);
d.NextDay();
d.TimezoneHours = earth.OffsetUtcHours;
d.shour = day.sunRise.TotalDays;
return(d);
}
}
d.year = -1;
d.month = -1;
d.day = -1;
d.TimezoneHours = earth.OffsetUtcHours;
d.shour = day.sunRise.TotalDays;
return(d);
}
public static XmlDocument GetSankrantiListXml(GCLocation loc, GregorianDateTime vcStart, GregorianDateTime vcEnd)
{
GregorianDateTime d = new GregorianDateTime();
int zodiac;
XmlDocument doc = new XmlDocument();
XmlElement e1 = doc.CreateElement("xml");
doc.AppendChild(e1);
XmlElement e2, e3;
GCEarthData earth = loc.GetEarthData();
// open file
d.Set(vcStart);
e2 = doc.CreateElement("request");
e1.AppendChild(e2);
e2.SetAttribute("name", "Sankranti");
e2.SetAttribute("version", GCStrings.getString(130));
e2.SetAttribute("longitude", loc.Longitude.ToString());
e2.SetAttribute("latitude", loc.Latitude.ToString());
e2.SetAttribute("timezone", loc.OffsetUtcHours.ToString());
e2.SetAttribute("startdate", vcStart.ToString());
e2.SetAttribute("enddate", vcEnd.ToString());
e2 = doc.CreateElement("result");
e2.SetAttribute("name", "SankrantiList");
e1.AppendChild(e2);
while (d.IsBeforeThis(vcEnd))
{
d.Set(GCSankranti.GetNextSankranti(d, earth, out zodiac));
d.InitWeekDay();
e3 = doc.CreateElement("sank");
e2.AppendChild(e3);
e3.SetAttribute("date", d.ToString());
e3.SetAttribute("dayweekid", d.dayOfWeek.ToString());
e3.SetAttribute("dayweek", GCCalendar.GetWeekdayName(d.dayOfWeek));
e3.SetAttribute("time", d.LongTimeString());
e3.SetAttribute("rasi", zodiac.ToString());
e3.SetAttribute("rasiName", GCRasi.GetName(zodiac));
e3.SetAttribute("rasiNameEn", GCRasi.GetNameEn(zodiac));
d.NextDay();
d.NextDay();
}
return(doc);
}
private int CalculateYoga(GCLocation loc, GregorianDateTime vcEnd, GregorianDateTime vcAdd, GCEarthData earth, GCConfigRatedEvents rec)
{
int nData = 0;
GregorianDateTime vcNext = new GregorianDateTime();
int ndst;
vcAdd.SubtractDays(1);
while (vcAdd.IsBeforeThis(vcEnd))
{
nData = GCYoga.GetNextYogaStart(earth, vcAdd, out vcNext);
if (vcNext.GetDayInteger() < vcEnd.GetDayInteger())
{
vcNext.InitWeekDay();
vcNext.AddHours(loc.TimeZone.GetBiasMinutesForDay(vcNext) / 60.0);
//AddEvent(vcNext, CoreEventType.CCTYPE_YOGA, nData, ndst);
AddRating(vcNext, rec.rateYoga[nData], rec.rateYoga[Prev(nData, 27)]);
}
else
{
break;
}
vcAdd.Set(vcNext);
vcAdd.shour += 0.2;
if (vcAdd.shour >= 1.0)
{
vcAdd.shour -= 1.0;
vcAdd.NextDay();
}
}
return(nData);
}
private void CalculateSunRasi(GCLocation loc, GregorianDateTime vcEnd, GregorianDateTime vcAdd, GCConfigRatedEvents rec)
{
int nData;
GregorianDateTime vcNext = new GregorianDateTime();
int ndst;
GCEarthData earth = loc.GetEarthData();
vcAdd.SubtractDays(30);
while (vcAdd.IsBeforeThis(vcEnd))
{
vcNext.Set(GCSankranti.GetNextSankranti(vcAdd, earth, out nData));
if (vcNext.GetDayInteger() < vcEnd.GetDayInteger())
{
vcNext.InitWeekDay();
vcNext.AddHours(loc.TimeZone.GetBiasMinutesForDay(vcNext) / 60.0);
//AddEvent(vcNext, CoreEventType.CCTYPE_SANK, nData, ndst);
AddRating(vcNext, rec.rateGrahaRasi[0, nData], rec.rateGrahaRasi[0, Prev(nData, 12)]);
}
else
{
break;
}
vcAdd.Set(vcNext);
vcAdd.NextDay();
}
}
public static int GetNextYogaStart(GCEarthData ed, GregorianDateTime startDate, out GregorianDateTime nextDate)
{
double phi = 40.0 / 3.0;
double l1, l2, longitudeSun;
double jday = startDate.GetJulianComplete();
double xj;
double longitudeMoon;
GregorianDateTime d = new GregorianDateTime();
d.Set(startDate);
GregorianDateTime xd = new GregorianDateTime();
double scan_step = 0.5;
int prev_tit = 0;
int new_tit = -1;
double ayanamsha = GCAyanamsha.GetAyanamsa(jday);
longitudeMoon = GCCoreAstronomy.GetMoonLongitude(d, ed);
longitudeSun = GCCoreAstronomy.GetSunLongitude(d, ed);
l1 = GCMath.putIn360(longitudeMoon + longitudeSun - 2 * ayanamsha);
prev_tit = Convert.ToInt32(Math.Floor(l1 / phi));
int counter = 0;
while (counter < 20)
{
xj = jday;
xd.Set(d);
jday += scan_step;
d.shour += scan_step;
if (d.shour > 1.0)
{
d.shour -= 1.0;
d.NextDay();
}
longitudeMoon = GCCoreAstronomy.GetMoonLongitude(d, ed);
longitudeSun = GCCoreAstronomy.GetSunLongitude(d, ed);
l2 = GCMath.putIn360(longitudeMoon + longitudeSun - 2 * ayanamsha);
new_tit = Convert.ToInt32(Math.Floor(l2 / phi));
if (prev_tit != new_tit)
{
jday = xj;
d.Set(xd);
scan_step *= 0.5;
counter++;
continue;
}
else
{
l1 = l2;
}
}
nextDate = d;
return(new_tit);
}
public static int GetNextMoonRasi(GCEarthData ed, GregorianDateTime startDate, out GregorianDateTime nextDate)
{
double phi = 30.0;
double l1, l2, longitudeMoon;
double jday = startDate.GetJulianComplete();
GregorianDateTime d = new GregorianDateTime();
d.Set(startDate);
double ayanamsa = GCAyanamsha.GetAyanamsa(jday);
double scan_step = 0.5;
int prev_naks = 0;
int new_naks = -1;
double xj;
GregorianDateTime xd = new GregorianDateTime();
longitudeMoon = GCCoreAstronomy.GetMoonLongitude(d, ed);
l1 = GCMath.putIn360(longitudeMoon - ayanamsa);
prev_naks = GCMath.IntFloor(l1 / phi);
int counter = 0;
while (counter < 20)
{
xj = jday;
xd.Set(d);
jday += scan_step;
d.shour += scan_step;
if (d.shour > 1.0)
{
d.shour -= 1.0;
d.NextDay();
}
longitudeMoon = GCCoreAstronomy.GetMoonLongitude(d, ed);
l2 = GCMath.putIn360(longitudeMoon - ayanamsa);
new_naks = GCMath.IntFloor(l2 / phi);
if (prev_naks != new_naks)
{
jday = xj;
d.Set(xd);
scan_step *= 0.5;
counter++;
continue;
}
else
{
l1 = l2;
}
}
nextDate = new GregorianDateTime();
nextDate.Set(d);
return(new_naks);
}
private int CalculateNaksatras(GCLocation loc, GregorianDateTime vcEnd, GregorianDateTime vcAdd, GCEarthData earth, GCConfigRatedEvents rec)
{
int nData = 0;
GregorianDateTime vcNext;
int ndst;
bool prevNaksatraValid = false;
GregorianDateTime prevNaksatra = new GregorianDateTime();
vcAdd.SubtractDays(1);
while (vcAdd.IsBeforeThis(vcEnd))
{
nData = GCNaksatra.GetNextNaksatra(earth, vcAdd, out vcNext);
if (vcNext.GetDayInteger() < vcEnd.GetDayInteger())
{
vcNext.InitWeekDay();
vcNext.AddHours(loc.TimeZone.GetBiasMinutesForDay(vcNext) / 60.0);
//AddEvent(vcNext, CoreEventType.CCTYPE_NAKS, nData, ndst);
AddRating(vcNext, rec.rateNaksatra[nData], rec.rateNaksatra[Prev(nData, 27)]);
if (prevNaksatraValid)
{
double padaLength = (vcNext.GetJulianComplete() - prevNaksatra.GetJulianComplete()) / 4.0;
for (int j = 0; j < 4; j++)
{
//AddEvent(prevNaksatra, CoreEventType.CCTYPE_NAKS_PADA1 + j, nData, ndst);
int prevPada = (nData * 4 + j + 107) % 108;
AddRating(vcNext, rec.rateNaksatraPada[nData, j], rec.rateNaksatraPada[prevPada / 4, prevPada % 4]);
prevNaksatra.shour += padaLength;
prevNaksatra.NormalizeValues();
}
}
prevNaksatra.Set(vcNext);
prevNaksatraValid = true;
}
else
{
break;
}
vcAdd.Set(vcNext);
vcAdd.shour += 0.2;
if (vcAdd.shour >= 1.0)
{
vcAdd.shour -= 1.0;
vcAdd.NextDay();
}
}
return(nData);
}
/*********************************************************************/
/* Finds next time when rasi is changed */
/* */
/* startDate - starting date and time, timezone member must be valid */
/* zodiac [out] - found zodiac sign into which is changed */
/* */
/*********************************************************************/
public static GregorianDateTime GetNextSankranti(GregorianDateTime startDate, GCEarthData earth, out int zodiac)
{
GregorianDateTime d = new GregorianDateTime();
double step = 1.0;
int count = 0;
double ld, prev;
int prev_rasi, new_rasi;
GregorianDateTime prevday;
d.Set(startDate);
//d.ChangeTimeZone(0.0);
//d.shour = 0.0;
zodiac = 0;
prev = GCMath.putIn360(GCCoreAstronomy.GetSunLongitude(d, earth) - GCAyanamsha.GetAyanamsa(d.GetJulian()));
prev_rasi = GCMath.IntFloor(prev / 30.0);
while (count < 20)
{
prevday = new GregorianDateTime();
prevday.Set(d);
d.shour += step;
if (d.shour > 1.0)
{
d.shour -= 1.0;
d.NextDay();
}
ld = GCMath.putIn360(GCCoreAstronomy.GetSunLongitude(d, earth) - GCAyanamsha.GetAyanamsa(d.GetJulian()));
new_rasi = GCMath.IntFloor(ld / 30.0);
if (prev_rasi != new_rasi)
{
zodiac = new_rasi;
//v uplynulom dni je sankranti
step *= 0.5;
d.Set(prevday);
count++;
continue;
}
}
return(d);
}
public int GetDaylightTimeStartDate(int nYear, ref GregorianDateTime vcStart)
{
vcStart.day = 1;
vcStart.month = Convert.ToInt32(StartDst.Month);
vcStart.year = nYear;
if (StartDst.Type == 1)
{
vcStart.day = StartDst.Day;
}
else
{
if (StartDst.Month == 5)
{
vcStart.day = GregorianDateTime.GetMonthMaxDays(nYear, StartDst.Month);
vcStart.InitWeekDay();
while (vcStart.dayOfWeek != StartDst.Day)
{
vcStart.PreviousDay();
vcStart.dayOfWeek = (vcStart.dayOfWeek + 6) % 7;
}
}
else
{
vcStart.day = 1;
vcStart.InitWeekDay();
while (vcStart.dayOfWeek != StartDst.Day)
{
vcStart.NextDay();
vcStart.dayOfWeek = (vcStart.dayOfWeek + 1) % 7;
}
vcStart.day += StartDst.Week * 7;
}
}
vcStart.shour = 1 / 24.0;
return(0);
}
public static int writeGaurabdaTithiXml(string fileName, GCLocation loc, GaurabdaDate vaStart, GaurabdaDate vaEnd)
{
int gyearA = vaStart.gyear;
int gyearB = vaEnd.gyear;
int gmasa = vaStart.masa;
int gpaksa = vaStart.tithi / 15;
int gtithi = vaStart.tithi % 15;
if (gyearB < gyearA)
{
gyearB = gyearA;
}
using (StreamWriter xml = new StreamWriter(fileName))
{
xml.Write("<xml>\n");
xml.Write("\t<request name=\"Tithi\" version=\"");
xml.Write(GCStrings.getString(130));
xml.Write("\">\n");
xml.Write("\t\t<arg name=\"longitude\" val=\"");
xml.Write(loc.Longitude);
xml.Write("\" />\n");
xml.Write("\t\t<arg name=\"latitude\" val=\"");
xml.Write(loc.Latitude);
xml.Write("\" />\n");
xml.Write("\t\t<arg name=\"timezone\" val=\"");
xml.Write(loc.OffsetUtcHours);
xml.Write("\" />\n");
if (gyearA > 1500)
{
xml.Write("\t\t<arg name=\"year-start\" val=\"");
xml.Write(gyearA);
xml.Write("\" />\n");
xml.Write("\t\t<arg name=\"year-end\" val=\"");
xml.Write(gyearB);
xml.Write("\" />\n");
}
else
{
xml.Write("\t\t<arg name=\"gaurabdayear-start\" val=\"");
xml.Write(gyearA);
xml.Write("\" />\n");
xml.Write("\t\t<arg name=\"gaurabdayear-end\" val=\"");
xml.Write(gyearB);
xml.Write("\" />\n");
}
xml.Write("\t\t<arg name=\"masa\" val=\"");
xml.Write(gmasa);
xml.Write("\" />\n");
xml.Write("\t\t<arg name=\"paksa\" val=\"");
xml.Write(gpaksa);
xml.Write("\" />\n");
xml.Write("\t\t<arg name=\"tithi\" val=\"");
xml.Write(gtithi);
xml.Write("\" />\n");
xml.Write("\t</request>\n");
xml.Write("\t<result name=\"Tithi\">\n");
GCEarthData earth = loc.GetEarthData();
GregorianDateTime vcs = new GregorianDateTime(), vce = new GregorianDateTime(), today = new GregorianDateTime();
GCHourTime sun = new GCHourTime();
int A, B;
double sunrise;
GCAstroData day = new GCAstroData();
int oTithi, oPaksa, oMasa, oYear;
if (gyearA > 1500)
{
A = gyearA - 1487;
B = gyearB - 1485;
}
else
{
A = gyearA;
B = gyearB;
}
for (; A <= B; A++)
{
vcs.Set(GCTithi.CalcTithiEnd(A, gmasa, gpaksa, gtithi, earth, out vce));
if (gyearA > 1500)
{
if ((vcs.year < gyearA) || (vcs.year > gyearB))
{
continue;
}
}
oTithi = gpaksa * 15 + gtithi;
oMasa = gmasa;
oPaksa = gpaksa;
oYear = 0;
xml.Write("\t<celebration\n");
// xml.Write("\t\t<tithi\n";
xml.Write("\t\trtithi=\"");
xml.Write(GCTithi.GetName(oTithi));
xml.Write("\"\n");
xml.Write("\t\trmasa=\"");
xml.Write(GCMasa.GetName(oMasa));
xml.Write("\"\n");
xml.Write("\t\trpaksa=\"");
xml.Write(oPaksa != 0 ? "Gaura" : "Krsna");
xml.Write("\"\n");
// test ci je ksaya
today.Set(vcs);
today.shour = 0.5;
sun = GCSunData.CalcSunrise(today, earth);
sunrise = sun.TotalDays;
if (sunrise < vcs.shour)
{
today.Set(vce);
sun = GCSunData.CalcSunrise(today, earth);
sunrise = sun.TotalDays;
if (sunrise < vce.shour)
{
// normal type
vcs.NextDay();
xml.Write("\t\ttype=\"normal\"\n");
}
else
{
// ksaya
vcs.NextDay();
day.DayCalc(vcs, earth);
oTithi = day.sunRise.Tithi;
oPaksa = day.sunRise.Paksa;
oMasa = day.MasaCalc(vcs, earth);
oYear = day.GaurabdaYear;
xml.Write("\t\ttype=\"ksaya\"\n");
}
}
else
{
// normal, alebo prvy den vriddhi
today.Set(vce);
sun = GCSunData.CalcSunrise(today, earth);
if (sun.TotalDays < vce.shour)
{
// first day of vriddhi type
xml.Write("\t\ttype=\"vriddhi\"\n");
}
else
{
// normal
xml.Write("\t\ttype=\"normal\"\n");
}
}
xml.Write("\t\tdate=\"");
xml.Write(vcs);
xml.Write("\"\n");
xml.Write("\t\totithi=\"");
xml.Write(GCTithi.GetName(oTithi));
xml.Write("\"\n");
xml.Write("\t\tomasa=\"");
xml.Write(GCMasa.GetName(oMasa));
xml.Write("\"\n");
xml.Write("\t\topaksa=\"");
xml.Write(GCPaksa.GetName(oPaksa));
xml.Write("\"\n");
xml.Write("\t/>\n");
}
xml.Write("\t</result>\n");
xml.Write("</xml>\n");
}
return(1);
}
/*********************************************************************/
/* Calculates Date of given Tithi */
/*********************************************************************/
public static GregorianDateTime CalcTithiDate(int nGYear, int nMasa, int nPaksa, int nTithi, GCEarthData earth)
{
int i = 0, gy = 0;
GregorianDateTime d = new GregorianDateTime();
GCAstroData day = new GCAstroData();
int tithi = 0;
int counter = 0;
UInt32 tmp = 0;
if (nGYear >= 464 && nGYear < 572)
{
tmp = GCAstroData.gGaurBeg[(nGYear - 464) * 26 + nMasa * 2 + nPaksa];
d.month = Convert.ToInt32((tmp & 0x3e0) >> 5);
d.day = Convert.ToInt32(tmp & 0x1f);
d.year = Convert.ToInt32(tmp & 0xfffc00) >> 10;
d.TimezoneHours = earth.OffsetUtcHours;
d.NextDay();
day.DayCalc(d, earth);
day.Masa = day.MasaCalc(d, earth);
gy = day.GaurabdaYear;
}
else
{
//d = GetFirstDayOfYear(earth, nGYear + 1486);
d.day = 15;
d.month = 2 + nMasa;
d.year = nGYear + 1486;
if (d.month > 12)
{
d.month -= 12;
d.year++;
}
d.shour = 0.5;
d.TimezoneHours = earth.OffsetUtcHours;
i = 0;
do
{
d.AddDays(13);
day.DayCalc(d, earth);
day.Masa = day.MasaCalc(d, earth);
gy = day.GaurabdaYear;
i++;
}while (((day.sunRise.Paksa != nPaksa) || (day.Masa != nMasa)) && (i <= 30));
}
if (i >= 30)
{
d.year = d.month = d.day = -1;
return(d);
}
// we found masa and paksa
// now we have to find tithi
tithi = nTithi + nPaksa * 15;
if (day.sunRise.Tithi == tithi)
{
// loc1
// find tithi juncts in this day and according to that times,
// look in previous or next day for end and start of this tithi
d.PreviousDay();
day.DayCalc(d, earth);
if ((day.sunRise.Tithi > tithi) && (day.sunRise.Paksa != nPaksa))
{
d.NextDay();
}
return(d);
}
if (day.sunRise.Tithi < tithi)
{
// do increment of date until nTithi == tithi
// but if nTithi > tithi
// then do decrement of date
counter = 0;
while (counter < 16)
{
d.NextDay();
day.DayCalc(d, earth);
if (day.sunRise.Tithi == tithi)
{
return(d);
}
if ((day.sunRise.Tithi < tithi) && (day.sunRise.Paksa != nPaksa))
{
return(d);
}
if (day.sunRise.Tithi > tithi)
{
return(d);
}
counter++;
}
// somewhere is error
d.year = d.month = d.day = 0;
return(d);
}
else
{
// do decrement of date until nTithi <= tithi
counter = 0;
while (counter < 16)
{
d.PreviousDay();
day.DayCalc(d, earth);
if (day.sunRise.Tithi == tithi)
{
return(d);
}
if ((day.sunRise.Tithi > tithi) && (day.sunRise.Paksa != nPaksa))
{
d.NextDay();
return(d);
}
if (day.sunRise.Tithi < tithi)
{
d.NextDay();
return(d);
}
counter++;
}
// somewhere is error
d.year = d.month = d.day = 0;
return(d);
}
// now we know the type of day-accurancy
// nType = 0 means, that we dont found a day
// nType = 1 means, we find day, when tithi was present at sunrise
// nType = 2 means, we found day, when tithi started after sunrise
// but ended before next sunrise
//
return(d);
}
public static GregorianDateTime CalcTithiEndEx(GregorianDateTime vcStart, int GYear, int nMasa, int nPaksa, int nTithi, GCEarthData earth, out GregorianDateTime endTithi)
{
int i, gy, nType;
GregorianDateTime d = new GregorianDateTime();
GCAstroData day = new GCAstroData();
int tithi;
int counter;
GregorianDateTime start = new GregorianDateTime(), end = new GregorianDateTime();
// SUNDATA sun;
// MOONDATA moon;
double sunrise;
start.shour = -1.0;
end.shour = -1.0;
start.day = start.month = start.year = -1;
end.day = end.month = end.year = -1;
/* d = GetFirstDayOfYear(earth, nGYear + 1486);
* d.shour = 0.5;
* d.TimeZone = earth.tzone;
*/
d.Set(vcStart);
i = 0;
do
{
d.AddDays(13);
day.DayCalc(d, earth);
day.Masa = day.MasaCalc(d, earth);
gy = day.GaurabdaYear;
i++;
}while (((day.sunRise.Paksa != nPaksa) || (day.Masa != nMasa)) && (i <= 30));
if (i >= 30)
{
d.year = d.month = d.day = -1;
endTithi = null;
return(d);
}
// we found masa and paksa
// now we have to find tithi
tithi = nTithi + nPaksa * 15;
if (day.sunRise.Tithi == tithi)
{
// loc1
// find tithi juncts in this day and according to that times,
// look in previous or next day for end and start of this tithi
nType = 1;
}
else
{
if (day.sunRise.Tithi < tithi)
{
// do increment of date until nTithi == tithi
// but if nTithi > tithi
// then do decrement of date
counter = 0;
while (counter < 30)
{
d.NextDay();
day.DayCalc(d, earth);
if (day.sunRise.Tithi == tithi)
{
goto cont_2;
}
if ((day.sunRise.Tithi < tithi) && (day.sunRise.Paksa != nPaksa))
{
d.PreviousDay();
goto cont_2;
}
if (day.sunRise.Tithi > tithi)
{
d.PreviousDay();
goto cont_2;
}
counter++;
}
// somewhere is error
d.year = d.month = d.day = 0;
nType = 0;
}
else
{
// do decrement of date until nTithi <= tithi
counter = 0;
while (counter < 30)
{
d.PreviousDay();
day.DayCalc(d, earth);
if (day.sunRise.Tithi == tithi)
{
goto cont_2;
}
if ((day.sunRise.Tithi > tithi) && (day.sunRise.Paksa != nPaksa))
{
goto cont_2;
}
if (day.sunRise.Tithi < tithi)
{
goto cont_2;
}
counter++;
}
// somewhere is error
d.year = d.month = d.day = 0;
nType = 0;
}
cont_2:
if (day.sunRise.Tithi == tithi)
{
// do the same as in loc1
nType = 1;
}
else
{
// nTithi != tithi and nTithi < tithi
// but on next day is nTithi > tithi
// that means we will find start and the end of tithi
// in this very day or on next day before sunrise
nType = 2;
}
}
// now we know the type of day-accurancy
// nType = 0 means, that we dont found a day
// nType = 1 means, we find day, when tithi was present at sunrise
// nType = 2 means, we found day, when tithi started after sunrise
// but ended before next sunrise
//
sunrise = day.sunRise.TotalDays;
endTithi = new GregorianDateTime();
if (nType == 1)
{
GregorianDateTime d1, d2;
d.shour = sunrise;
GCTithi.GetPrevTithiStart(earth, d, out d1);
//d = d1;
//d.shour += 0.02;
GCTithi.GetNextTithiStart(earth, d, out d2);
endTithi.Set(d2);
return(d1);
}
else if (nType == 2)
{
GregorianDateTime d1, d2;
d.shour = sunrise;
GCTithi.GetNextTithiStart(earth, d, out d1);
d.Set(d1);
d.shour += 0.1;
d.NormalizeValues();
GCTithi.GetNextTithiStart(earth, d, out d2);
endTithi.Set(d2);
return(d1);
}
// if nType == 0, then this algoritmus has some failure
if (nType == 0)
{
d.year = 0;
d.month = 0;
d.day = 0;
d.shour = 0.0;
endTithi.Set(d);
}
else
{
d.Set(start);
endTithi.Set(end);
}
return(d);
}
public int GetNextAscendentStart(GregorianDateTime startDate, out GregorianDateTime nextDate)
{
int countOfElements = 12;
double phi = 360.0 / countOfElements;
double l1, l2;
double jday = startDate.GetJulianDetailed();
double xj;
GregorianDateTime d = new GregorianDateTime(startDate);
// we are calculating in UTC timezone
d.shour -= startDate.TimezoneHours / 24.0;
d.TimezoneHours = 0.0;
GregorianDateTime xd = new GregorianDateTime();
double scan_step = 1.0 / 48.0;
int prev_tit = 0;
int new_tit = -1;
l1 = GetAscendantDegrees(jday);
prev_tit = Convert.ToInt32(Math.Floor(l1 / phi));
//Debugger.Log(0, "", "=== FROM " + startDate.LongTimeString() + " ===\n");
//Debugger.Log(0, "", string.Format("{0:00}: {1} {2} {3} {4} \n", 99, l1, prev_tit, jday, d.LongTimeString()));
int counter = 0;
while (counter < 20)
{
xj = jday;
xd.Set(d);
jday += scan_step;
d.shour += scan_step;
if (d.shour > 1.0)
{
d.shour -= 1.0;
d.NextDay();
}
l2 = GetAscendantDegrees(jday);
new_tit = Convert.ToInt32(Math.Floor(l2 / phi));
//Debugger.Log(0, "", string.Format("{0:00}: {1} {2} {3} {4} \n", counter, l1, prev_tit, jday, d.LongTimeString()));
if (prev_tit != new_tit)
{
jday = xj;
d.Set(xd);
scan_step *= 0.5;
counter++;
//Debugger.Log(0, "", string.Format(" Going back to {0}\n", d.LongTimeString()));
continue;
}
else
{
l1 = l2;
}
}
// date D was calculated in timezone UTC+0h, therefore
// we have to convert it into timezone of input datetime
nextDate = new GregorianDateTime(d);
nextDate.shour += startDate.TimezoneHours / 24.0;
nextDate.TimezoneHours = startDate.TimezoneHours;
nextDate.NormalizeValues();
return(new_tit);
}
/*********************************************************************/
/* */
/* */
/* */
/* */
/* */
/*********************************************************************/
public static double GetNextConjunction(GregorianDateTime test_date, out GregorianDateTime found, bool this_conj, GCEarthData earth)
{
double phi = 12.0;
double l1, l2, longitudeSun, longitudeMoon;
if (this_conj)
{
test_date.shour += 0.2;
test_date.NormalizeValues();
}
double jday = test_date.GetJulianComplete();
double xj;
GregorianDateTime d = new GregorianDateTime();
d.Set(test_date);
GregorianDateTime xd = new GregorianDateTime();
double scan_step = 1.0;
int prev_tit = 0;
int new_tit = -1;
longitudeMoon = GCCoreAstronomy.GetMoonLongitude(d, earth);
longitudeSun = GCSunData.GetSunLongitude(d);
l1 = GCMath.putIn180(longitudeMoon - longitudeSun);
prev_tit = GCMath.IntFloor(l1 / phi);
int counter = 0;
while (counter < 20)
{
xj = jday;
xd.Set(d);
jday += scan_step;
d.shour += scan_step;
if (d.shour > 1.0)
{
d.shour -= 1.0;
d.NextDay();
}
longitudeMoon = GCCoreAstronomy.GetMoonLongitude(d, earth);
longitudeSun = GCSunData.GetSunLongitude(d);
l2 = GCMath.putIn180(longitudeMoon - longitudeSun);
new_tit = GCMath.IntFloor(l2 / phi);
if (prev_tit < 0 && new_tit >= 0)
{
jday = xj;
d.Set(xd);
scan_step *= 0.5;
counter++;
continue;
}
else
{
l1 = l2;
}
prev_tit = new_tit;
}
found = new GregorianDateTime();
found.Set(d);
return(longitudeSun);
}
///////////////////////////////////////////////////////////////////////
// GET PREVIOUS CONJUNCTION OF THE SUN AND MOON
//
// THIS IS HELP FUNCTION FOR GetPrevConjunction(VCTIME test_date,
// VCTIME &found, bool this_day, EARTHDATA earth)
//
// looking for previous sun-moon conjunction
// it starts from input day from 12:00 AM (noon) UTC
// so output can be the same day
// if output is the same day, then conjunction occurs between 00:00 AM and noon of this day
//
// date - input / output
// earth - input
// return value - sun longitude in degs
//
// error is when return value is greater than 999.0 deg
public static double GetPrevConjunction(ref GregorianDateTime date, GCEarthData earth)
{
int bCont = 32;
double prevSun = 0.0, prevMoon = 0.0, prevDiff = 0.0;
double nowSun = 0.0, nowMoon = 0.0, nowDiff = 0.0;
// SUNDATA sun;
double jd, longitudeMoon;
GregorianDateTime d = new GregorianDateTime();
d.Set(date);
d.shour = 0.5;
d.TimezoneHours = 0.0;
jd = d.GetJulian();//GetJulianDay(d.year, d.month, d.day);
// set initial data for input day
// NOTE: for grenwich
//SunPosition(d, earth, sun);
longitudeMoon = GCCoreAstronomy.GetMoonLongitude(d, earth);
prevSun = GCSunData.GetSunLongitude(d);
prevMoon = longitudeMoon;
prevDiff = GCMath.putIn180(prevSun - prevMoon);
do
{
// shift to day before
d.PreviousDay();
jd -= 1.0;
// calculation
//SunPosition(d, earth, sun);
longitudeMoon = GCCoreAstronomy.GetMoonLongitude(d, earth);
nowSun = GCSunData.GetSunLongitude(d);
nowMoon = longitudeMoon;
nowDiff = GCMath.putIn180(nowSun - nowMoon);
// if difference of previous has another sign than now calculated
// then it is the case that moon was faster than sun and he
//printf(" prevsun=%f\nprevmoon=%f\nnowsun=%f\nnowmoon=%f\n", prevSun, prevMoon, nowSun, nowMoon);
if (IsConjunction(nowMoon, nowSun, prevSun, prevMoon))
{
// now it calculates actual time and zodiac of conjunction
double x;
if (prevDiff == nowDiff)
{
return(0);
}
x = Math.Abs(nowDiff) / Math.Abs(prevDiff - nowDiff);
if (x < 0.5)
{
d.shour = x + 0.5;
}
else
{
d.NextDay();
d.shour = x - 0.5;
}
date.Set(d);
double other = GCSunData.GetSunLongitude(d);
// double other2 = nowSun + (prevSun - nowSun)*x;
GCMath.putIn360(prevSun);
GCMath.putIn360(nowSun);
if (Math.Abs(prevSun - nowSun) > 10.0)
{
return(GCMath.putIn180(nowSun) + (GCMath.putIn180(prevSun) - GCMath.putIn180(nowSun)) * x);
}
else
{
return(nowSun + (prevSun - nowSun) * x);
}
// return other2;
}
prevSun = nowSun;
prevMoon = nowMoon;
prevDiff = nowDiff;
bCont--;
}while (bCont >= 0);
return(1000.0);
}
/// <summary>
/// Main function for CALENDAR calculations
/// </summary>
/// <param name="loc"></param>
/// <param name="begDate"></param>
/// <param name="iCount"></param>
/// <returns></returns>
public int CalculateCalendar(GCLocation loc, GregorianDateTime begDate, int iCount)
{
int i, weekday;
int nTotalCount = BEFORE_DAYS + iCount + BEFORE_DAYS;
GregorianDateTime date;
GCEarthData earth;
int lastMasa = 0;
int lastGYear = 0;
String tempStr;
bool bCalcMoon = (GCDisplaySettings.Current.getValue(4) > 0 || GCDisplaySettings.Current.getValue(5) > 0);
m_nCount = 0;
m_Location = loc;
m_vcStart = new GregorianDateTime(begDate);
m_vcCount = iCount;
earth = loc.GetEarthData();
// alokacia pola
m_pData = new VAISNAVADAY[nTotalCount];
// inicializacia poctovych premennych
m_nCount = nTotalCount;
m_PureCount = iCount;
date = new GregorianDateTime();
date.Set(begDate);
date.shour = 0.0;
date.TimezoneHours = loc.OffsetUtcHours;
date.SubtractDays(BEFORE_DAYS);
date.InitWeekDay();
weekday = date.dayOfWeek;
GCFestivalSpecialExecutor exec = new GCFestivalSpecialExecutor(this);
long utcDayStart = -1;
long utcDayEnd = -1;
// 1
// initialization of days
for (i = 0; i < nTotalCount; i++)
{
m_pData[i] = new VAISNAVADAY();
m_pData[i].date = new GregorianDateTime(date);
m_pData[i].date.dayOfWeek = weekday;
date.NextDay();
weekday = (weekday + 1) % 7;
m_pData[i].moonrise.SetValue(-1);
m_pData[i].moonset.SetValue(-1);
if (utcDayStart < 0)
{
utcDayStart = Convert.ToInt64(m_pData[i].date.GetJulianComplete() * 86400);
}
m_pData[i].UtcDayStart = utcDayStart;
utcDayStart += 86400;
}
for (i = 0; i < nTotalCount; i++)
{
m_pData[i].Previous = i > 0 ? m_pData[i - 1] : null;
m_pData[i].Next = i < nTotalCount - 1 ? m_pData[i + 1] : null;
}
TResultCoreEvents recentEvents = null;
// 2
// calculating DST data and core events
foreach (VAISNAVADAY t in m_pData)
{
t.BiasMinutes = loc.TimeZone.GetBiasMinutesForDay(t.date);
if (t.Previous != null)
{
t.DstDayType = (t.BiasMinutes == 0
? (t.Previous.BiasMinutes == 0 ? DstTypeChange.DstOff : DstTypeChange.DstEnd)
: (t.Previous.BiasMinutes == 0 ? DstTypeChange.DstStart : DstTypeChange.DstOn));
}
utcDayStart = t.UtcDayStart - (t.Previous != null ? t.Previous.BiasMinutes * 60 : 0);
utcDayEnd = t.UtcDayStart + 86400 - t.BiasMinutes * 60;
if (recentEvents == null || recentEvents.p_events.Count == 0 || recentEvents.Year != t.date.year)
{
recentEvents = GCCoreAstronomy.GetCoreEventsYear(loc, t.date.year);
}
recentEvents.GetCoreEvents(t.coreEvents, utcDayStart, utcDayEnd);
utcDayStart = utcDayEnd;
}
// 3
if (bCalcMoon)
{
foreach (VAISNAVADAY t in m_pData)
{
GCMoonData.CalcMoonTimes(earth, t.date, Convert.ToDouble(t.BiasMinutes / 60.0), out t.moonrise, out t.moonset);
if (!EkadasiOnly && GCDisplaySettings.Current.getValue(GCDS.CAL_MOON_RISE) != 0 && t.moonrise.hour >= 0)
{
tempStr = t.Format(GCStrings.Localized("Moonrise {moonRiseTime} ({dstSig})"));
t.AddEvent(DisplayPriorities.PRIO_MOON, GCDS.CAL_MOON_RISE, tempStr);
}
if (!EkadasiOnly && GCDisplaySettings.Current.getValue(GCDS.CAL_MOON_SET) != 0 && t.moonset.hour >= 0)
{
tempStr = t.Format(GCStrings.Localized("Moonset {moonSetTime} ({dstSig})"));
t.AddEvent(DisplayPriorities.PRIO_MOON, GCDS.CAL_MOON_SET, tempStr);
}
}
}
// 4
// init of astro data
foreach (VAISNAVADAY t in m_pData)
{
t.astrodata = new GCAstroData();
t.astrodata.DayCalc(t.date, earth);
t.coreEvents.InsertByTime(new TCoreEvent()
{
Time = t.UtcDayStart + (long)t.astrodata.sunArunodaya.TotalSeconds,
nType = CoreEventType.CCTYPE_S_ARUN
});
t.coreEvents.InsertByTime(new TCoreEvent()
{
Time = t.UtcDayStart + (long)t.astrodata.sunRise.TotalSeconds,
nType = CoreEventType.CCTYPE_S_RISE
});
t.coreEvents.InsertByTime(new TCoreEvent()
{
Time = t.UtcDayStart + (long)t.astrodata.sunNoon.TotalSeconds,
nType = CoreEventType.CCTYPE_S_NOON
});
t.coreEvents.InsertByTime(new TCoreEvent()
{
Time = t.UtcDayStart + (long)t.astrodata.sunSet.TotalSeconds,
nType = CoreEventType.CCTYPE_S_SET
});
}
// 5
// init of masa
foreach (VAISNAVADAY t in m_pData)
{
if (t.Previous == null ||
t.astrodata.sunRise.Paksa != t.Previous.astrodata.sunRise.Paksa)
{
t.astrodata.MasaCalc(t.date, earth);
lastMasa = t.astrodata.Masa;
lastGYear = t.astrodata.GaurabdaYear;
}
t.astrodata.Masa = lastMasa;
t.astrodata.GaurabdaYear = lastGYear;
if (!EkadasiOnly && GCDisplaySettings.Current.getValue(GCDS.CAL_SUN_LONG) != 0)
{
tempStr = string.Format("{0}: {1} (*)", GCStrings.Localized("Sun Longitude")
, t.astrodata.sunRise.longitude);
t.AddEvent(DisplayPriorities.PRIO_ASTRO, GCDS.CAL_SUN_LONG, tempStr);
}
if (!EkadasiOnly && GCDisplaySettings.Current.getValue(GCDS.CAL_MOON_LONG) != 0)
{
tempStr = string.Format("{0}: {1} (*)", GCStrings.Localized("Moon Longitude")
, t.astrodata.sunRise.longitudeMoon);
t.AddEvent(DisplayPriorities.PRIO_ASTRO, GCDS.CAL_MOON_LONG, tempStr);
}
if (!EkadasiOnly && GCDisplaySettings.Current.getValue(GCDS.CAL_AYANAMSHA) != 0)
{
tempStr = string.Format("{0} {1} ({2}) (*)", GCStrings.Localized("Ayanamsha")
, t.astrodata.Ayanamsa
, GCAyanamsha.GetAyanamsaName(GCAyanamsha.GetAyanamsaType()));
t.AddEvent(DisplayPriorities.PRIO_ASTRO, GCDS.CAL_AYANAMSHA, tempStr);
}
if (!EkadasiOnly && GCDisplaySettings.Current.getValue(GCDS.CAL_JULIAN) != 0)
{
tempStr = string.Format("{0} {1} (*)", GCStrings.Localized("Julian Time")
, t.astrodata.JulianDay);
t.AddEvent(DisplayPriorities.PRIO_ASTRO, GCDS.CAL_JULIAN, tempStr);
}
}
if (!EkadasiOnly && GCDisplaySettings.Current.getValue(GCDS.CAL_MASA_CHANGE) != 0)
{
String str;
foreach (VAISNAVADAY t in m_pData)
{
if (t.Previous != null && t.Previous.astrodata.Masa != t.astrodata.Masa)
{
str = t.Format(GCStrings.Localized("First day of {masaName} Masa"));
t.AddEvent(DisplayPriorities.PRIO_MASA_CHANGE, GCDS.CAL_MASA_CHANGE, str);
}
if (t.Next != null && t.Next.astrodata.Masa != t.astrodata.Masa)
{
str = t.Format(GCStrings.Localized("Last day of {masaName} Masa"));
t.AddEvent(DisplayPriorities.PRIO_MASA_CHANGE, GCDS.CAL_MASA_CHANGE, str);
}
}
}
if (!EkadasiOnly && GCDisplaySettings.Current.getValue(GCDS.CAL_DST_CHANGE) != 0)
{
foreach (VAISNAVADAY t in m_pData)
{
if (t.Previous != null && t.Previous.BiasMinutes == 0 && t.BiasMinutes != 0)
{
t.AddEvent(DisplayPriorities.PRIO_DST_CHANGE, GCDS.CAL_DST_CHANGE, GCStrings.Localized("First day of Daylight Saving Time"));
}
else if (t.Next != null && t.BiasMinutes != 0 && t.Next.BiasMinutes == 0)
{
t.AddEvent(DisplayPriorities.PRIO_DST_CHANGE, GCDS.CAL_DST_CHANGE, GCStrings.Localized("Last day of Daylight Saving Time"));
}
}
}
// 6
// init of mahadvadasis
foreach (VAISNAVADAY t in m_pData)
{
t.MahadvadasiCalc(earth);
}
//
// init for Ekadasis
foreach (VAISNAVADAY t in m_pData)
{
if (t.Previous == null)
{
continue;
}
t.EkadasiCalc(earth);
if (t.Previous.nFastID == FastType.FAST_EKADASI)
{
t.CalculateEParana(earth);
}
}
if (EkadasiOnly)
{
return(1);
}
// init ksaya data
// init of second day of vriddhi
CalculateKsayaVriddhiTithis();
//
// calculate sankrantis
CalculateSankrantis();
// 7
// init of festivals
AddMatchedFestivals(exec);
//
// apply daylight saving time
ApplyDaylightSavingHours();
//
// resolve festivals fasting
for (i = BEFORE_DAYS; i < m_PureCount + BEFORE_DAYS; i++)
{
ResolveFestivalsFasting(i);
if (GCDisplaySettings.Current.getValue(GCDS.CAL_SUN_RISE) != 0)
{
tempStr = string.Format("{0}-{1}-{2} {3} - {4} - {5} ({6})",
GCStrings.Localized("Sunrise"),
GCStrings.Localized("Noon"),
GCStrings.Localized("Sunset"),
m_pData[i].astrodata.sunRise.ToShortTimeString(),
m_pData[i].astrodata.sunNoon.ToShortTimeString(),
m_pData[i].astrodata.sunSet.ToShortTimeString(),
GCStrings.GetDSTSignature(m_pData[i].BiasMinutes));
m_pData[i].AddEvent(DisplayPriorities.PRIO_SUN, GCDS.CAL_SUN_RISE, tempStr);
}
if (GCDisplaySettings.Current.getValue(GCDS.CAL_SUN_SANDHYA) != 0)
{
tempStr = string.Format("{0}: {1} | {2} | {3} ({4})",
GCStrings.Localized("Sandhyas"),
m_pData[i].astrodata.sunRise.ShortSandhyaString(),
m_pData[i].astrodata.sunNoon.ShortSandhyaString(),
m_pData[i].astrodata.sunSet.ShortSandhyaString(),
GCStrings.GetDSTSignature(m_pData[i].BiasMinutes));
m_pData[i].AddEvent(DisplayPriorities.PRIO_SUN, GCDS.CAL_SUN_SANDHYA, tempStr);
}
if (GCDisplaySettings.Current.getValue(GCDS.CAL_BRAHMA_MUHURTA) != 0)
{
tempStr = string.Format("{0}: {1} ({2})",
GCStrings.Localized("Brahma Muhurta"),
m_pData[i].astrodata.sunRise.ShortMuhurtaString(-2),
GCStrings.GetDSTSignature(m_pData[i].BiasMinutes));
m_pData[i].AddEvent(DisplayPriorities.PRIO_SUN, GCDS.CAL_BRAHMA_MUHURTA, tempStr);
}
}
if (GCDisplaySettings.Current.getValue(GCDS.CAL_COREEVENTS) != 0)
{
int number = 1;
foreach (VAISNAVADAY today in m_pData)
{
foreach (TCoreEvent tde in today.coreEvents)
{
m_pData[i].AddEvent(DisplayPriorities.PRIO_ASTRO + number, GCDS.CAL_COREEVENTS,
tde.TypeString + " " + today.GetCoreEventTime(tde));
number++;
}
}
}
// sorting day events according priority
VAISNAVAEVENTComparer vec = new VAISNAVAEVENTComparer();
foreach (VAISNAVADAY t in m_pData)
{
t.dayEvents.Sort(vec);
}
return(1);
}
public static int writeGaurabdaNextTithiXml(string fileName, GCLocation loc, GregorianDateTime vcStart, GaurabdaDate vaStart)
{
int gmasa, gpaksa, gtithi;
GCHourTime sunRise;
using (StreamWriter xml = new StreamWriter(fileName))
{
gmasa = vaStart.masa;
gpaksa = vaStart.tithi / 15;
gtithi = vaStart.tithi % 15;
xml.Write("<xml>\n");
xml.Write("\t<request name=\"Tithi\" version=\"");
xml.Write(GCStrings.getString(130));
xml.Write("\">\n");
xml.Write("\t\t<arg name=\"longitude\" val=\"");
xml.Write(loc.Longitude);
xml.Write("\" />\n");
xml.Write("\t\t<arg name=\"latitude\" val=\"");
xml.Write(loc.Latitude);
xml.Write("\" />\n");
xml.Write("\t\t<arg name=\"timezone\" val=\"");
xml.Write(loc.OffsetUtcHours);
xml.Write("\" />\n");
xml.Write("\t\t<arg name=\"start date\" val=\"");
xml.Write(vcStart);
xml.Write("\" />\n");
xml.Write("\t\t<arg name=\"masa\" val=\"");
xml.Write(gmasa);
xml.Write("\" />\n");
xml.Write("\t\t<arg name=\"paksa\" val=\"");
xml.Write(gpaksa);
xml.Write("\" />\n");
xml.Write("\t\t<arg name=\"tithi\" val=\"");
xml.Write(gtithi);
xml.Write("\" />\n");
xml.Write("\t</request>\n");
xml.Write("\t<result name=\"Tithi\">\n");
GCEarthData earth = loc.GetEarthData();
GregorianDateTime vcs = new GregorianDateTime(), vce = new GregorianDateTime(), today = new GregorianDateTime();
//GCSunData sun = new GCSunData();
int A;
double sunrise;
GCAstroData day = new GCAstroData();
int oTithi, oPaksa, oMasa, oYear;
today.Set(vcStart);
today.PreviousDay();
vcStart.SubtractDays(15);
for (A = 0; A <= 3; A++)
{
vcs.Set(GCTithi.CalcTithiEndEx(vcStart, 0, gmasa, gpaksa, gtithi, earth, out vce));
if (!vcs.IsBeforeThis(today))
{
oTithi = gpaksa * 15 + gtithi;
oMasa = gmasa;
oPaksa = gpaksa;
oYear = 0;
xml.Write("\t<celebration\n");
// xml.Write("\t\t<tithi\n";
xml.Write("\t\trtithi=\"");
xml.Write(GCTithi.GetName(oTithi));
xml.Write("\"\n");
xml.Write("\t\trmasa=\"");
xml.Write(GCMasa.GetName(oMasa));
xml.Write("\"\n");
xml.Write("\t\trpaksa=\"");
xml.Write((oPaksa != 0 ? "Gaura" : "Krsna"));
xml.Write("\"\n");
// test ci je ksaya
today.Set(vcs);
today.shour = 0.5;
sunRise = GCSunData.CalcSunrise(today, earth);
sunrise = sunRise.TotalDays;
if (sunrise < vcs.shour)
{
today.Set(vce);
sunRise = GCSunData.CalcSunrise(today, earth);
sunrise = sunRise.TotalDays;
if (sunrise < vce.shour)
{
// normal type
vcs.NextDay();
xml.Write("\t\ttype=\"normal\"\n");
}
else
{
// ksaya
vcs.NextDay();
day.DayCalc(vcs, earth);
oTithi = day.sunRise.Tithi;
oPaksa = day.sunRise.Paksa;
oMasa = day.MasaCalc(vcs, earth);
oYear = day.GaurabdaYear;
xml.Write("\t\ttype=\"ksaya\"\n");
}
}
else
{
// normal, alebo prvy den vriddhi
today.Set(vce);
sunRise = GCSunData.CalcSunrise(today, earth);
if (sunRise.TotalDays < vce.shour)
{
// first day of vriddhi type
xml.Write("\t\ttype=\"vriddhi\"\n");
}
else
{
// normal
xml.Write("\t\ttype=\"normal\"\n");
}
}
xml.Write("\t\tdate=\"");
xml.Write(vcs);
xml.Write("\"\n");
xml.Write("\t\totithi=\"");
xml.Write(GCTithi.GetName(oTithi));
xml.Write("\"\n");
xml.Write("\t\tomasa=\"");
xml.Write(GCMasa.GetName(oMasa));
xml.Write("\"\n");
xml.Write("\t\topaksa=\"");
xml.Write(GCPaksa.GetName(oPaksa));
xml.Write("\"\n");
xml.Write("\t/>\n");
break;
}
else
{
vcStart.Set(vcs);
vcs.NextDay();
}
}
xml.Write("\t</result>\n");
xml.Write("</xml>\n");
}
return(1);
}
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);
}
}
private void CalculateEvents(GCLocation loc, GregorianDateTime vcStart, GregorianDateTime vcEnd, GCConfigRatedEvents rec)
{
GregorianDateTime vc = new GregorianDateTime(vcStart);
double biasHours = 0;
ShowAboveOnly = rec.useAcceptLimit;
ShowLongerThan = rec.useMinPeriodLength;
RatingsAboveOnly = rec.acceptLimit;
MinIntervalLength = rec.minPeriodLength;
//bool hasAscendent = rec.RequiredGrahaRasi(9);
bool hasRahuKalam = rec.rateKalas[KalaType.KT_RAHU_KALAM].Usable;
bool hasYamaghanti = rec.rateKalas[KalaType.KT_YAMA_GHANTI].Usable;
bool hasGulikalam = rec.rateKalas[KalaType.KT_GULI_KALAM].Usable;
bool hasAbhijit = rec.rateKalas[KalaType.KT_ABHIJIT].Usable;
bool hasMuhurta = rec.RequiredMuhurta();
p_events.Clear();
EarthLocation = new GCLocation(loc);
DateStart = new GregorianDateTime(vcStart);
DateEnd = new GregorianDateTime(vcEnd);
bool hasPreviousSunset = false;
GregorianDateTime previousSunset = new GregorianDateTime();
GregorianDateTime vcAdd = new GregorianDateTime(vcStart);
GregorianDateTime vcNext = new GregorianDateTime();
GCEarthData earth = loc.GetEarthData();
vcAdd.InitWeekDay();
GCHourTime sunRise, sunSet;
double muhurtaLength = 48.0 / 1440.0;
GregorianDateTime muhurtaDate = new GregorianDateTime();
while (vcAdd.IsBeforeThis(vcEnd))
{
biasHours = loc.TimeZone.GetBiasMinutesForDay(vcAdd) / 60.0;
sunRise = GCSunData.CalcSunrise(vcAdd, earth);
sunSet = GCSunData.CalcSunset(vcAdd, earth);
sunRise.AddMinutes(biasHours * 60);
sunSet.AddMinutes(biasHours * 60);
vcAdd.shour = sunRise.TotalDays;
vcAdd.InitWeekDay();
//AddEvent(vcAdd, CoreEventType.CCTYPE_S_RISE, 0, ndst);
AddRating(vcAdd, rec.rateDayHours[0], rec.rateDayHours[3]);
AddRating(vcAdd, rec.rateDay[0], rec.rateDayHours[1]);
//AddEvent(vcAdd, CoreEventType.CCTYPE_DAY_OF_WEEK, vcAdd.dayOfWeek, ndst);
AddRating(vcAdd, rec.weekday[vcAdd.dayOfWeek], rec.weekday[Prev(vcAdd.dayOfWeek, 7)]);
if (hasPreviousSunset)
{
previousSunset.shour += (vcAdd.shour + 1.0 - previousSunset.shour) / 2;
previousSunset.NormalizeValues();
//AddEvent(previousSunset, CoreEventType.CCTYPE_S_MIDNIGHT, 0, ndst);
AddRating(previousSunset, rec.rateDayHours[3], rec.rateDayHours[2]);
}
if (hasMuhurta)
{
muhurtaDate.Set(vcAdd);
muhurtaDate.shour -= 2 * muhurtaLength;
muhurtaDate.NormalizeValues();
for (int j = 0; j < 30; j++)
{
//AddEvent(muhurtaDate, CoreEventType.CCTYPE_DAY_MUHURTA, (j + 28) % 30, ndst);
int mi = (j + 28) % 30;
AddRating(muhurtaDate, rec.rateMuhurta[mi], rec.rateMuhurta[Prev(mi, 30)]);
}
}
vcAdd.shour = (sunRise.TotalDays + sunSet.TotalDays) / 2;
//AddEvent(vcAdd, CoreEventType.CCTYPE_S_NOON, 0, ndst);
AddRating(vcAdd, rec.rateDayHours[1], rec.rateDayHours[0]);
vcAdd.shour = sunSet.TotalDays;
//AddEvent(vcAdd, CoreEventType.CCTYPE_S_SET, 0, ndst);
AddRating(vcAdd, rec.rateDayHours[2], rec.rateDayHours[1]);
AddRating(vcAdd, rec.rateDay[1], rec.rateDayHours[0]);
previousSunset.Set(vcAdd);
hasPreviousSunset = true;
if (hasRahuKalam)
{
CalculateKalam(vcAdd, sunRise.TotalDays, sunSet.TotalDays, KalaType.KT_RAHU_KALAM, rec);
}
if (hasYamaghanti)
{
CalculateKalam(vcAdd, sunRise.TotalDays, sunSet.TotalDays, KalaType.KT_YAMA_GHANTI, rec);
}
if (hasGulikalam)
{
CalculateKalam(vcAdd, sunRise.TotalDays, sunSet.TotalDays, KalaType.KT_GULI_KALAM, rec);
}
if (hasAbhijit)
{
CalculateKalam(vcAdd, sunRise.TotalDays, sunSet.TotalDays, KalaType.KT_ABHIJIT, rec);
}
vcAdd.NextDay();
}
if (rec.RequiredGrahaRasi(9))
{
GCAscendant asc = new GCAscendant();
asc.Earth = EarthLocation.GetEarthData();
asc.CurrentDateTime = new GregorianDateTime(vc);
while (asc.GetNextAscendantBefore(vcEnd))
{
AddRating(asc.CurrentDateTime.TimeWithOffset(biasHours), rec.rateGrahaRasi[9, asc.CurrentSign], rec.rateGrahaRasi[9, Prev(asc.CurrentSign, 12)]);
asc.CurrentDateTime.AddHours(0.5);
}
}
if (rec.RequiredTithi())
{
vcAdd.Set(vc);
vcAdd.shour = 0.0;
CalculateTithis(loc, vcEnd, vcAdd, earth, rec);
}
if (rec.RequiredNaksatra())
{
vcAdd.Set(vc);
vcAdd.shour = 0.0;
CalculateNaksatras(loc, vcEnd, vcAdd, earth, rec);
}
if (rec.RequiredYoga())
{
vcAdd.Set(vc);
vcAdd.shour = 0.0;
CalculateYoga(loc, vcEnd, vcAdd, earth, rec);
}
for (int i = 0; i < 9; i++)
{
if (rec.RequiredGrahaRasi(i))
{
vcAdd.Set(vc);
vcAdd.shour = 0.0;
CalculatePlanetRasi(i, loc, vcEnd, vcAdd, rec);
//CalculateSunRasi(loc, vcEnd, vcAdd, rec);
}
if (rec.RequiredGrahaHouse(i))
{
vcAdd.Set(vc);
vcAdd.shour = 0.0;
CalculatePlanetHouse(i, loc, vcEnd, vcAdd, rec);
//CalculateSunRasi(loc, vcEnd, vcAdd, rec);
}
}
/*if (rec.RequiredGrahaRasi(1))
* {
* vcAdd.Set(vc);
* vcAdd.shour = 0.0;
* CalculateMoonRasi(loc, vcEnd, vcAdd, earth, rec);
*
* }*/
if (rec.RequiredMoonTimes())
{
vcAdd.Set(vc);
vcAdd.shour = 0.0;
CalculateMoonTimes(vcEnd, vcAdd, earth, rec);
}
IComparer <GCRatedMoment> C = new GCRatedMoment.ComparerClass();
p_ratings.Sort(C);
}
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();
}
/// <summary>
/// Calculation of Masa List
/// </summary>
/// <param name="loc">Location</param>
/// <param name="nYear">Starting year</param>
/// <param name="nCount">Number of years</param>
/// <returns></returns>
public int CalculateMasaList(GCLocation loc, int nYear, int nCount)
{
GCAstroData day = new GCAstroData();
GregorianDateTime d = new GregorianDateTime(), de = new GregorianDateTime(), t = new GregorianDateTime();
int lm = -1;
TResultMasaList mlist = this;
GCEarthData earth = loc.GetEarthData();
mlist.n_startYear = nYear;
mlist.n_countYears = nCount;
mlist.vc_start = new GregorianDateTime();
mlist.vc_end = new GregorianDateTime();
mlist.vc_start.Set(GCAstroData.GetFirstDayOfYear(earth, nYear));
mlist.vc_end.Set(GCAstroData.GetFirstDayOfYear(earth, nYear + nCount));
mlist.m_location = loc;
d.Set(mlist.vc_start);
de.Set(mlist.vc_end);
int i = 0;
int prev_paksa = -1;
int current = 0;
while (d.IsBeforeThis(de))
{
day.DayCalc(d, earth);
if (prev_paksa != day.sunRise.Paksa)
{
day.Masa = day.MasaCalc(d, earth);
if (lm != day.Masa)
{
if (lm >= 0)
{
t.Set(d);
t.PreviousDay();
if (mlist.arr.Count <= current)
{
mlist.arr.Add(new TResultMasa());
}
mlist.arr[current].vc_end = new GregorianDateTime(t);
current++;
}
lm = day.Masa;
if (mlist.arr.Count <= current)
{
mlist.arr.Add(new TResultMasa());
}
mlist.arr[current].masa = day.Masa;
mlist.arr[current].year = day.GaurabdaYear;
mlist.arr[current].vc_start = new GregorianDateTime(d);
}
}
prev_paksa = day.sunRise.Paksa;
d.NextDay();
i++;
}
mlist.arr[current].vc_end = new GregorianDateTime(d);
current++;
mlist.n_countMasa = current;
return(1);
}
