public void TestBasics() { // Check that our JD computation is the same as the book's (p. 88) CalendarAstronomer astro = new CalendarAstronomer(); IBM.ICU.Util.GregorianCalendar cal3 = new IBM.ICU.Util.GregorianCalendar( IBM.ICU.Util.TimeZone.GetTimeZone("GMT"), Locale.US); IBM.ICU.Text.DateFormat d3 = IBM.ICU.Text.DateFormat.GetDateTimeInstance(cal3, IBM.ICU.Text.DateFormat.MEDIUM, IBM.ICU.Text.DateFormat.MEDIUM, Locale.US); cal3.Clear(); cal3.Set(IBM.ICU.Util.Calendar.YEAR, 1980); cal3.Set(IBM.ICU.Util.Calendar.MONTH, IBM.ICU.Util.Calendar.JULY); cal3.Set(IBM.ICU.Util.Calendar.DATE, 27); astro.SetDate(cal3.GetTime()); double jd = astro.GetJulianDay() - 2447891.5d; double exp = -3444; if (jd == exp) { Logln(d3.Format(cal3.GetTime()) + " => " + jd); } else { Errln("FAIL: " + d3.Format(cal3.GetTime()) + " => " + jd + ", expected " + exp); } // cal3.clear(); // cal3.set(cal3.YEAR, 1990); // cal3.set(cal3.MONTH, Calendar.JANUARY); // cal3.set(cal3.DATE, 1); // cal3.add(cal3.DATE, -1); // astro.setDate(cal3.getTime()); // astro.foo(); }
/// <summary> /// Compute the month and date on which this holiday falls in the year /// containing the date "date". First figure out which date Easter lands on /// in this year, and then add the offset for this holiday to get the right /// date. /// <p> /// The algorithm here is taken from the <a /// href="http://www.faqs.org/faqs/calendars/faq/">Calendar FAQ</a>. /// </summary> /// private DateTime ComputeInYear(DateTime date, GregorianCalendar cal) { if (cal == null) { cal = calendar; } lock (cal) { cal.SetTime(date); int year = cal.Get(IBM.ICU.Util.Calendar.YEAR); int g = year % 19; // "Golden Number" of year - 1 int i = 0; // # of days from 3/21 to the Paschal full moon int j = 0; // Weekday (0-based) of Paschal full moon if (ILOG.J2CsMapping.Util.DateUtil.After(cal.GetTime(), cal.GetGregorianChange())) { // We're past the Gregorian switchover, so use the Gregorian // rules. int c = year / 100; int h = (c - c / 4 - (8 * c + 13) / 25 + 19 * g + 15) % 30; i = h - (h / 28) * (1 - (h / 28) * (29 / (h + 1)) * ((21 - g) / 11)); j = (year + year / 4 + i + 2 - c + c / 4) % 7; } else { // Use the old Julian rules. i = (19 * g + 15) % 30; j = (year + year / 4 + i) % 7; } int l = i - j; int m = 3 + (l + 40) / 44; // 1-based month in which Easter falls int d = l + 28 - 31 * (m / 4); // Date of Easter within that month cal.Clear(); cal.Set(IBM.ICU.Util.Calendar.ERA, IBM.ICU.Util.GregorianCalendar.AD); cal.Set(IBM.ICU.Util.Calendar.YEAR, year); cal.Set(IBM.ICU.Util.Calendar.MONTH, m - 1); // 0-based cal.Set(IBM.ICU.Util.Calendar.DATE, d); cal.GetTime(); // JDK 1.1.2 bug workaround cal.Add(IBM.ICU.Util.Calendar.DATE, daysAfterEaster); return(cal.GetTime()); } }
public void TestLunarPosition() { IBM.ICU.Util.GregorianCalendar gc = new IBM.ICU.Util.GregorianCalendar(new SimpleTimeZone(0, "UTC")); CalendarAstronomer astro = new CalendarAstronomer(); double[][] tests = { new double[] { 1979, 2, 26, 16, 0x0, -0.3778379118188744d, -0.1399698825594198d } }; Logln(""); for (int i = 0; i < tests.Length; i++) { gc.Clear(); gc.Set((int)tests[i][0], (int)tests[i][1] - 1, (int)tests[i][2], (int)tests[i][3], (int)tests[i][4]); astro.SetDate(gc.GetTime()); IBM.ICU.Impl.CalendarAstronomer.Equatorial result = astro.GetMoonPosition(); if (result.ascension != tests[i][5]) { if ((float)result.ascension == (float)tests[i][5]) { Logln("result.ascension(" + result.ascension + ") != tests[i][5](" + tests[i][5] + ") in double for test " + i); } else { Errln("FAIL: result.ascension(" + result.ascension + ") != tests[i][5](" + tests[i][5] + ") for test " + i); } } if (result.declination != tests[i][6]) { if ((float)result.declination == (float)tests[i][6]) { Logln("result.declination(" + result.declination + ") != tests[i][6](" + tests[i][6] + ") in double for test " + i); } else { Errln("FAIL: result.declination(" + result.declination + ") != tests[i][6](" + tests[i][6] + ") for test " + i); } } } }
public void TestIsOn() { // jb 1901 SimpleHoliday sh = new SimpleHoliday(IBM.ICU.Util.Calendar.AUGUST, 15, "Doug's Day", 1958, 2058); IBM.ICU.Util.Calendar cal_0 = new IBM.ICU.Util.GregorianCalendar(); cal_0.Clear(); cal_0.Set(IBM.ICU.Util.Calendar.YEAR, 2000); cal_0.Set(IBM.ICU.Util.Calendar.MONTH, IBM.ICU.Util.Calendar.AUGUST); cal_0.Set(IBM.ICU.Util.Calendar.DAY_OF_MONTH, 15); DateTime d0 = cal_0.GetTime(); cal_0.Add(IBM.ICU.Util.Calendar.SECOND, 1); DateTime d1 = cal_0.GetTime(); cal_0.Add(IBM.ICU.Util.Calendar.SECOND, -2); DateTime d2 = cal_0.GetTime(); cal_0.Add(IBM.ICU.Util.Calendar.DAY_OF_MONTH, 1); DateTime d3 = cal_0.GetTime(); cal_0.Add(IBM.ICU.Util.Calendar.SECOND, 1); DateTime d4 = cal_0.GetTime(); cal_0.Add(IBM.ICU.Util.Calendar.SECOND, -2); cal_0.Set(IBM.ICU.Util.Calendar.YEAR, 1957); DateTime d5 = cal_0.GetTime(); cal_0.Set(IBM.ICU.Util.Calendar.YEAR, 1958); DateTime d6 = cal_0.GetTime(); cal_0.Set(IBM.ICU.Util.Calendar.YEAR, 2058); DateTime d7 = cal_0.GetTime(); cal_0.Set(IBM.ICU.Util.Calendar.YEAR, 2059); DateTime d8 = cal_0.GetTime(); DateTime[] dates = { d0, d1, d2, d3, d4, d5, d6, d7, d8 }; bool[] isOns = { true, true, false, true, false, false, true, true, false }; for (int i = 0; i < dates.Length; ++i) { DateTime d = dates[i]; Logln("\ndate: " + d); bool isOn = sh.IsOn(d); Logln("isOnDate: " + isOn); if (isOn != isOns[i]) { Errln("date: " + d + " should be on Doug's Day!"); } DateTime h = sh.FirstAfter(d); Logln("firstAfter: " + h); } }
public void TestSolarLongitude() { IBM.ICU.Util.GregorianCalendar gc = new IBM.ICU.Util.GregorianCalendar(new SimpleTimeZone(0, "UTC")); CalendarAstronomer astro = new CalendarAstronomer(); double[][] tests = { new double[] { 1980, 7, 27, 0x0, 0x0, 2.166442986535465d, 2.2070499713207730d, 0.3355704075759270d }, new double[] { 1988, 7, 27, 0x0, 0x0, 2.167484927693959d, 2.2081183335606176d, 0.3353093444275315d } }; Logln(""); for (int i = 0; i < tests.Length; i++) { gc.Clear(); gc.Set((int)tests[i][0], (int)tests[i][1] - 1, (int)tests[i][2], (int)tests[i][3], (int)tests[i][4]); astro.SetDate(gc.GetTime()); double longitude = astro.GetSunLongitude(); if (longitude != tests[i][5]) { if ((float)longitude == (float)tests[i][5]) { Logln("longitude(" + longitude + ") != tests[i][5](" + tests[i][5] + ") in double for test " + i); } else { Errln("FAIL: longitude(" + longitude + ") != tests[i][5](" + tests[i][5] + ") for test " + i); } } IBM.ICU.Impl.CalendarAstronomer.Equatorial result = astro.GetSunPosition(); if (result.ascension != tests[i][6]) { if ((float)result.ascension == (float)tests[i][6]) { Logln("result.ascension(" + result.ascension + ") != tests[i][6](" + tests[i][6] + ") in double for test " + i); } else { Errln("FAIL: result.ascension(" + result.ascension + ") != tests[i][6](" + tests[i][6] + ") for test " + i); } } if (result.declination != tests[i][7]) { if ((float)result.declination == (float)tests[i][7]) { Logln("result.declination(" + result.declination + ") != tests[i][7](" + tests[i][7] + ") in double for test " + i); } else { Errln("FAIL: result.declination(" + result.declination + ") != tests[i][7](" + tests[i][7] + ") for test " + i); } } } }
public void TestSunriseTimes() { // logln("Sunrise/Sunset times for San Jose, California, USA"); // CalendarAstronomer astro = new CalendarAstronomer(-121.55, 37.20); // TimeZone tz = TimeZone.getTimeZone("America/Los_Angeles"); // We'll use a table generated by the UNSO website as our reference // From: http://aa.usno.navy.mil/ // -Location: W079 25, N43 40 // -Rise and Set for the Sun for 2001 // -Zone: 4h West of Greenwich int[] USNO = { 6, 59, 19, 45, 6, 57, 19, 46, 6, 56, 19, 47, 6, 54, 19, 48, 6, 52, 19, 49, 6, 50, 19, 51, 6, 48, 19, 52, 6, 47, 19, 53, 6, 45, 19, 54, 6, 43, 19, 55, 6, 42, 19, 57, 6, 40, 19, 58, 6, 38, 19, 59, 6, 36, 20, 0, 6, 35, 20, 1, 6, 33, 20, 3, 6, 31, 20, 4, 6, 30, 20, 5, 6, 28, 20, 6, 6, 27, 20, 7, 6, 25, 20,8, 6, 23, 20, 10, 6, 22, 20, 11, 6, 20, 20, 12, 6, 19, 20, 13, 6, 17, 20, 14, 6, 16, 20, 16, 6, 14, 20, 17, 6, 13, 20, 18, 6, 11, 20, 19, }; Logln("Sunrise/Sunset times for Toronto, Canada"); CalendarAstronomer astro = new CalendarAstronomer(-(79 + 25 / 60), 43 + 40 / 60); // As of ICU4J 2.8 the ICU4J time zones implement pass-through // to the underlying JDK. Because of variation in the // underlying JDKs, we have to use a fixed-offset // SimpleTimeZone to get consistent behavior between JDKs. // The offset we want is [-18000000, 3600000] (raw, dst). // [aliu 10/15/03] // TimeZone tz = TimeZone.getTimeZone("America/Montreal"); IBM.ICU.Util.TimeZone tz = new SimpleTimeZone(-18000000 + 3600000, "Montreal(FIXED)"); IBM.ICU.Util.GregorianCalendar cal = new IBM.ICU.Util.GregorianCalendar(tz, Locale.US); IBM.ICU.Util.GregorianCalendar cal2 = new IBM.ICU.Util.GregorianCalendar(tz, Locale.US); cal.Clear(); cal.Set(IBM.ICU.Util.Calendar.YEAR, 2001); cal.Set(IBM.ICU.Util.Calendar.MONTH, IBM.ICU.Util.Calendar.APRIL); cal.Set(IBM.ICU.Util.Calendar.DAY_OF_MONTH, 1); cal.Set(IBM.ICU.Util.Calendar.HOUR_OF_DAY, 12); // must be near local noon for // getSunRiseSet to work IBM.ICU.Text.DateFormat df = IBM.ICU.Text.DateFormat.GetTimeInstance(cal, IBM.ICU.Text.DateFormat.MEDIUM, Locale.US); IBM.ICU.Text.DateFormat df2 = IBM.ICU.Text.DateFormat.GetDateTimeInstance(cal, IBM.ICU.Text.DateFormat.MEDIUM, IBM.ICU.Text.DateFormat.MEDIUM, Locale.US); IBM.ICU.Text.DateFormat day = IBM.ICU.Text.DateFormat.GetDateInstance(cal, IBM.ICU.Text.DateFormat.MEDIUM, Locale.US); for (int i = 0; i < 30; i++) { astro.SetDate(cal.GetTime()); DateTime sunrise = DateUtil.DateFromJavaMillis(astro.GetSunRiseSet(true)); DateTime sunset = DateUtil.DateFromJavaMillis(astro.GetSunRiseSet(false)); cal2.SetTime(cal.GetTime()); cal2.Set(IBM.ICU.Util.Calendar.SECOND, 0); cal2.Set(IBM.ICU.Util.Calendar.MILLISECOND, 0); cal2.Set(IBM.ICU.Util.Calendar.HOUR_OF_DAY, USNO[4 * i + 0]); cal2.Set(IBM.ICU.Util.Calendar.MINUTE, USNO[4 * i + 1]); DateTime exprise = cal2.GetTime(); cal2.Set(IBM.ICU.Util.Calendar.HOUR_OF_DAY, USNO[4 * i + 2]); cal2.Set(IBM.ICU.Util.Calendar.MINUTE, USNO[4 * i + 3]); DateTime expset = cal2.GetTime(); // Compute delta of what we got to the USNO data, in seconds int deltarise = Math.Abs((int)((sunrise.Ticks / 10000) - (exprise.Ticks / 10000)) / 1000); int deltaset = Math.Abs((int)((sunset.Ticks / 10000) - (expset.Ticks / 10000)) / 1000); // Allow a deviation of 0..MAX_DEV seconds // It would be nice to get down to 60 seconds, but at this // point that appears to be impossible without a redo of the // algorithm using something more advanced than Duffett-Smith. int MAX_DEV = 180; if (deltarise > MAX_DEV || deltaset > MAX_DEV) { if (deltarise > MAX_DEV) { Errln("FAIL: " + day.Format(cal.GetTime()) + ", Sunrise: " + df2.Format(sunrise) + " (USNO " + df.Format(exprise) + " d=" + deltarise + "s)"); } else { Logln(day.Format(cal.GetTime()) + ", Sunrise: " + df.Format(sunrise) + " (USNO " + df.Format(exprise) + ")"); } if (deltaset > MAX_DEV) { Errln("FAIL: " + day.Format(cal.GetTime()) + ", Sunset: " + df2.Format(sunset) + " (USNO " + df.Format(expset) + " d=" + deltaset + "s)"); } else { Logln(day.Format(cal.GetTime()) + ", Sunset: " + df.Format(sunset) + " (USNO " + df.Format(expset) + ")"); } } else { Logln(day.Format(cal.GetTime()) + ", Sunrise: " + df.Format(sunrise) + " (USNO " + df.Format(exprise) + ")" + ", Sunset: " + df.Format(sunset) + " (USNO " + df.Format(expset) + ")"); } cal.Add(IBM.ICU.Util.Calendar.DATE, 1); } // CalendarAstronomer a = new CalendarAstronomer(-(71+5/60), 42+37/60); // cal.clear(); // cal.set(cal.YEAR, 1986); // cal.set(cal.MONTH, cal.MARCH); // cal.set(cal.DATE, 10); // cal.set(cal.YEAR, 1988); // cal.set(cal.MONTH, cal.JULY); // cal.set(cal.DATE, 27); // a.setDate(cal.getTime()); // long r = a.getSunRiseSet2(true); }