ruby-calendrical 2.0.1

data/ext/calendrical/calendrical.h

@@ -0,0 +1,99 @@
+/*
+  calendrical.h: Written by Tadayoshi Funaba 1997-2000,2014
+
+  This code is in the public domain, but any use of it
+  should publically acknowledge its source.
+
+  $Id: calendrical.h,v 1.6 2014-04-08 19:27:31+09 tadf Exp $
+*/
+
+#if !defined(__GNUC__) \
+ || !( __GNUC__ > 2 || (__GNUC__ == 2 && __GNUC_MINOR__ >= 6))
+#undef __attribute__
+#define __attribute__(x)
+#endif
+
+extern int gregorian_leap_year(int year) __attribute__ ((const));
+extern int gregorian_last_day_of_month(int year, int month) __attribute__ ((const));
+extern int gregorian_to_rd(int year, int month, int day) __attribute__ ((const));
+extern void rd_to_gregorian(int date, int *ryear, int *rmonth, int *rday);
+extern int kday_on_or_before(int date, int k) __attribute__ ((const));
+extern int iso_to_rd(int year, int week, int day) __attribute__ ((const));
+extern void rd_to_iso(int date, int *ryear, int *rweek, int *rday);
+extern int julian_leap_year(int year) __attribute__ ((const));
+extern int julian_last_day_of_month(int year, int month) __attribute__ ((const));
+extern int julian_to_rd(int year, int month, int day) __attribute__ ((const));
+extern void rd_to_julian(int date, int *ryear, int *rmonth, int *rday);
+extern int islamic_leap_year(int year) __attribute__ ((const));
+extern int islamic_last_day_of_month(int year, int month) __attribute__ ((const));
+extern int islamic_to_rd(int year, int month, int day) __attribute__ ((const));
+extern void rd_to_islamic(int date, int *ryear, int *rmonth, int *rday);
+extern int hebrew_leap_year(int year) __attribute__ ((const));
+extern int hebrew_last_month_of_year(int year) __attribute__ ((const));
+extern int hebrew_last_day_of_month(int year, int month) __attribute__ ((const));
+extern int hebrew_calendar_elapsed_days(int year) __attribute__ ((const));
+extern int hebrew_days_in_year(int year) __attribute__ ((const));
+extern int long_heshvan(int year) __attribute__ ((const));
+extern int short_kislev(int year) __attribute__ ((const));
+extern int hebrew_to_rd(int year, int month, int day) __attribute__ ((const));
+extern void rd_to_hebrew(int date, int *ryear, int *rmonth, int *rday);
+#if !defined(PLUTO) && !defined(MODULE)
+extern int independence_day(int year) __attribute__ ((const));
+#endif
+#ifndef PLUTO
+extern int nth_kday(int year, int month, int n, int k) __attribute__ ((const));
+#endif
+#if !defined(PLUTO) && !defined(MODULE)
+extern int labor_day(int year) __attribute__ ((const));
+extern int memorial_day(int year) __attribute__ ((const));
+extern int daylight_savings_start(int year) __attribute__ ((const));
+extern int daylight_savings_end(int year) __attribute__ ((const));
+extern int christmas(int year) __attribute__ ((const));
+extern int advent(int year) __attribute__ ((const));
+extern int epiphany(int year) __attribute__ ((const));
+extern int eastern_orthodox_christmas(int year) __attribute__ ((const));
+#endif /* not PLUTO and not MODULE */
+extern int nicaean_rule_easter(int year) __attribute__ ((const));
+extern int easter(int year) __attribute__ ((const));
+#if !defined(PLUTO) && !defined(MODULE)
+extern int pentecost(int year) __attribute__ ((const));
+extern void islamic_date(int year, int month, int day, int date[3]);
+extern void mulad_al_nabi(int year, int date[3]);
+extern int yom_kippur(int year) __attribute__ ((const));
+extern int passover(int year) __attribute__ ((const));
+extern int purim(int year) __attribute__ ((const));
+extern int ta_anit_esther(int year) __attribute__ ((const));
+extern int tisha_b_av(int year) __attribute__ ((const));
+extern int hebrew_birthday(int birth_year, int birth_month, int birth_day, int year) __attribute__ ((const));
+extern int yahrzeit(int death_year, int death_month, int death_day, int year) __attribute__ ((const));
+#endif /* not PLUTO and not MODULE */
+extern int mayan_long_count_to_rd(int baktun, int katun, int tun, int uinal, int kin) __attribute__ ((const));
+extern void rd_to_mayan_long_count(int date, int *rbaktun, int *rkatun, int *rtun, int *ruinal, int *rkin);
+extern void rd_to_mayan_haab(int date, int *rmonth, int *rday);
+extern int mayan_haab_difference(int month1, int day1, int month2, int day2) __attribute__ ((const));
+extern int mayan_haab_on_or_before(int haab_month, int haab_day, int date) __attribute__ ((const));
+extern void rd_to_mayan_tzolkin(int date, int *rname, int *rnumber);
+extern int mayan_tzolkin_difference(int name1, int number1, int name2, int number2) __attribute__ ((const));
+extern int mayan_tzolkin_on_or_before(int name, int number, int date) __attribute__ ((const));
+extern int mayan_haab_tzolkin_on_or_before(int month, int day, int name, int number, int date) __attribute__ ((const));
+extern int french_leap_year(int year) __attribute__ ((const));
+extern int french_last_day_of_month(int year, int month) __attribute__ ((const));
+extern int french_to_rd(int year, int month, int day) __attribute__ ((const));
+extern void rd_to_french(int date, int *ryear, int *rmonth, int *rday);
+extern double solar_longitude(double days) __attribute__ ((const));
+extern double zodiac(double days) __attribute__ ((const));
+extern void rd_to_old_hindu_solar(int date, int *ryear, int *rmonth, int *rday);
+extern int old_hindu_solar_to_rd(int year, int month, int day) __attribute__ ((const));
+extern double lunar_longitude(double days) __attribute__ ((const));
+extern double lunar_phase(double days) __attribute__ ((const));
+extern double new_moon(double days) __attribute__ ((const));
+extern void rd_to_old_hindu_lunar(int date, int *ryear, int *rmonth, int *rleapmonth, int *rday);
+extern int old_hindu_lunar_precedes(int year1, int month1, int leap1, int day1, int year2, int month2, int leap2, int day2) __attribute__ ((const));
+extern int old_hindu_lunar_to_rd(int year, int month, int leapmonth, int day) __attribute__ ((const));
+
+
+/*
+Local Variables:
+c-basic-offset: 2
+End:
+*/

data/ext/calendrical/calendrical2.c

@@ -0,0 +1,802 @@
+/*
+  calendrical2.c (beta): Written by Tadayoshi Funaba 1997,1999-2002,2014
+
+  This code is in the public domain, but any use of it
+  should publically acknowledge its source.
+
+  $Id: calendrical2.c,v 1.10 2014-04-27 09:42:02+09 tadf Exp $
+*/
+
+#define MODULE
+
+#ifndef NULL
+#define NULL 0
+#endif
+
+#include <stdlib.h>
+#include <math.h>
+#include "calendrical.h"
+#include "calendrical2.h"
+
+static div_t
+divmod(int m, int n)
+{
+  div_t d;
+
+  d = div(m, n);
+  if ((m < 0) != (n < 0) && d.rem != 0) {
+    d.quot -= 1;
+    d.rem += n;
+  }
+  return d;
+}
+
+#undef quotient
+#define quotient(m, n) i_quotient(m, n)
+static int
+i_quotient(int m, int n)
+{
+  return divmod(m, n).quot;
+}
+#undef mod
+#define mod(m, n) i_mod(m, n)
+static int
+i_mod(int m, int n)
+{
+  return divmod(m, n).rem;
+}
+#undef oddp
+#define oddp(n) ((n) % 2)
+
+static int
+approximate(int x, int y)
+{
+  return quotient(x, (x < 0) ? (y - 1) : y);
+}
+
+int
+world_leap_year(int year)
+{
+  return
+    mod(year, 4) == 0 &&
+    !(mod(year, 400) == 100 ||
+      mod(year, 400) == 200 ||
+      mod(year, 400) == 300);
+}
+
+int
+world_last_day_of_month(int year, int month)
+{
+  if (month == 6 && world_leap_year(year))
+    return 31;
+  switch (month) {
+  case  1: return 31;
+  case  2: return 30;
+  case  3: return 30;
+  case  4: return 31;
+  case  5: return 30;
+  case  6: return 30;
+  case  7: return 31;
+  case  8: return 30;
+  case  9: return 30;
+  case 10: return 31;
+  case 11: return 30;
+  case 12: return 31;
+  default: return 0;
+  }
+}
+
+int
+world_to_rd(int year, int month, int day)
+{
+  int sumres;
+
+  {
+    int temp, m;
+    for (temp = 0, m = 1;
+	 (m < month);
+	 temp = temp + world_last_day_of_month(year, m), m++)
+      ;
+    sumres = temp;
+  }
+  return day
+    + sumres
+    + (365 * (year - 1))
+    + quotient(year - 1, 4)
+    - quotient(year - 1, 100)
+    + quotient(year - 1, 400);
+}
+
+void
+rd_to_world(int date, int *ryear, int *rmonth, int *rday)
+{
+  int approx, month, day, year;
+  int sumres1, sumres2;
+
+  approx = approximate(date, 366);
+  {
+    int temp, y;
+    for (temp = 0, y = approx;
+	 (date >= world_to_rd(1 + y, 1, 1));
+	 temp = temp + 1, y++)
+      ;
+    sumres1 = temp;
+  }
+  year = approx
+    + sumres1;
+  {
+    int temp, m;
+    for (temp = 0, m = 1;
+	 (date > world_to_rd
+	  (year, m, world_last_day_of_month(year, m)));
+	 temp = temp + 1, m++)
+      ;
+    sumres2 = temp;
+  }
+  month = 1
+    + sumres2;
+  day = date
+    - (world_to_rd(year, month, 1) - 1);
+  if (rmonth)
+    *rmonth = month;
+  if (rday)
+    *rday = day;
+  if (ryear)
+    *ryear = year;
+}
+
+int
+rd_to_world_day_of_week(int date)
+{
+  int month, day, year;
+
+  rd_to_world(date, &year, &month, &day);
+  if (month == 12 && day == 31)
+    return 7;
+  else if (month == 6 && day == 31)
+    return 8;
+  else
+    switch ((month - 1) % 3) {
+    case 0: return (day - 1) % 7;
+    case 1: return (day + 2) % 7;
+    case 2: return (day + 4) % 7;
+    }
+  /*NOTREACHED*/
+  return -1;
+}
+
+int
+coptic_leap_year(int year)
+{
+  return mod(year, 4) == 3;
+}
+
+int
+coptic_last_day_of_month(int year, int month)
+{
+  if (month == 13 && coptic_leap_year(year))
+    return 6;
+  switch (month) {
+  case  1: return 30;
+  case  2: return 30;
+  case  3: return 30;
+  case  4: return 30;
+  case  5: return 30;
+  case  6: return 30;
+  case  7: return 30;
+  case  8: return 30;
+  case  9: return 30;
+  case 10: return 30;
+  case 11: return 30;
+  case 12: return 30;
+  case 13: return 5;
+  default: return 0;
+  }
+}
+
+#define EPOCH_COPTIC 103605
+#define BEFORE_COPTIC (103605 - 1)
+
+int
+coptic_to_rd(int year, int month, int day)
+{
+  int sumres;
+
+  {
+    int temp, m;
+    for (temp = 0, m = 1;
+	 (m < month);
+	 temp = temp + coptic_last_day_of_month(year, m), m++)
+      ;
+    sumres = temp;
+  }
+  return day
+    + sumres
+    + (365 * (year - 1))
+    + quotient(year, 4)
+    + BEFORE_COPTIC;
+}
+
+void
+rd_to_coptic(int date, int *ryear, int *rmonth, int *rday)
+{
+  int approx, month, day, year;
+  int sumres1, sumres2;
+
+  approx = approximate(date - BEFORE_COPTIC, 366);
+  {
+    int temp, y;
+    for (temp = 0, y = approx;
+	 (date >= coptic_to_rd(1 + y, 1, 1));
+	 temp = temp + 1, y++)
+      ;
+    sumres1 = temp;
+  }
+  year = approx
+    + sumres1;
+  {
+    int temp, m;
+    for (temp = 0, m = 1;
+	 (date > coptic_to_rd
+	  (year, m, coptic_last_day_of_month(year, m)));
+	 temp = temp + 1, m++)
+      ;
+    sumres2 = temp;
+  }
+  month = 1
+    + sumres2;
+  day = date - (coptic_to_rd(year, month, 1) - 1);
+  if (rmonth)
+    *rmonth = month;
+  if (rday)
+    *rday = day;
+  if (ryear)
+    *ryear = year;
+}
+
+int
+ethiopian_leap_year(int year)
+{
+  return coptic_leap_year(year);
+}
+
+int
+ethiopian_last_day_of_month(int year, int month)
+{
+  return coptic_last_day_of_month(year, month);
+}
+
+#define EPOCH_ETHIOPIAN 2796
+#define BEFORE_ETHIOPIAN (2796 - 1)
+
+int
+ethiopian_to_rd(int year, int month, int day)
+{
+  return
+    coptic_to_rd(year, month, day) - (EPOCH_COPTIC - EPOCH_ETHIOPIAN);
+}
+
+void
+rd_to_ethiopian(int date, int *ryear, int *rmonth, int *rday)
+{
+  rd_to_coptic
+    (date + (EPOCH_COPTIC - EPOCH_ETHIOPIAN),
+     ryear, rmonth, rday);
+}
+
+int
+jalali_leap_year(int year)
+{
+    return mod(29 + (8 * year), 33) < 8;
+}
+
+int
+jalali_last_day_of_month(int year, int month)
+{
+  if (month == 12 && jalali_leap_year(year))
+      return 30;
+  switch (month) {
+  case  1: return 31;
+  case  2: return 31;
+  case  3: return 31;
+  case  4: return 31;
+  case  5: return 31;
+  case  6: return 31;
+  case  7: return 30;
+  case  8: return 30;
+  case  9: return 30;
+  case 10: return 30;
+  case 11: return 30;
+  case 12: return 29;
+  default: return 0;
+  }
+}
+
+#define BEFORE_JALALI (226895 - 1)
+
+int
+jalali_to_rd(int year, int month, int day)
+{
+  int sumres;
+
+  {
+    int temp, m;
+    for (temp = 0, m = 1;
+	 (m < month);
+	 temp = temp + jalali_last_day_of_month(year, m), m++)
+      ;
+    sumres = temp;
+  }
+  return day
+    + sumres
+    + (365 * (year - 1))
+    + quotient(21 + (8 * year), 33)
+    + BEFORE_JALALI;
+}
+
+void
+rd_to_jalali(int date, int *ryear, int *rmonth, int *rday)
+{
+  int approx, month, day, year;
+  int sumres1, sumres2;
+
+  approx = approximate(date - BEFORE_JALALI, 366);
+  {
+    int temp, y;
+    for (temp = 0, y = approx;
+	 (date >= jalali_to_rd(1 + y, 1, 1));
+	 temp = temp + 1, y++)
+      ;
+    sumres1 = temp;
+  }
+  year = approx
+    + sumres1;
+  {
+    int temp, m;
+    for (temp = 0, m = 1;
+	 (date > jalali_to_rd
+	  (year, m, jalali_last_day_of_month(year, m)));
+	 temp = temp + 1, m++)
+      ;
+    sumres2 = temp;
+  }
+  month = 1
+    + sumres2;
+  day = date - (jalali_to_rd(year, month, 1) - 1);
+  if (rmonth)
+    *rmonth = month;
+  if (rday)
+    *rday = day;
+  if (ryear)
+    *ryear = year;
+}
+
+int
+bahai_leap_year(int year)
+{
+    return gregorian_leap_year(year + 1844);
+}
+
+int
+bahai_last_day_of_month(int year, int month)
+{
+  if (month == 19) {
+    if (bahai_leap_year(year))
+      return 5;
+    else
+      return 4;
+  }
+  return 19;
+}
+
+#define BEFORE_BAHAI (673222 - 1)
+
+int
+bahai_to_rd(int year, int month, int day)
+{
+  int sumres;
+
+  {
+    int temp, m;
+    for (temp = 0, m = 1;
+	 (m < month);
+	 temp = temp + bahai_last_day_of_month(year, m), m++)
+      ;
+    sumres = temp;
+  }
+  return day
+    + sumres
+    + gregorian_to_rd(year + 1843, 3, 20);
+}
+
+void
+rd_to_bahai(int date, int *ryear, int *rmonth, int *rday)
+{
+  int approx, month, day, year;
+  int sumres1, sumres2;
+
+  approx = approximate(date - BEFORE_BAHAI, 366);
+  {
+    int temp, y;
+    for (temp = 0, y = approx;
+	 (date >= bahai_to_rd(1 + y, 1, 1));
+	 temp = temp + 1, y++)
+      ;
+    sumres1 = temp;
+  }
+  year = approx
+    + sumres1;
+  {
+    int temp, m;
+    for (temp = 0, m = 1;
+	 (date > bahai_to_rd
+	  (year, m, bahai_last_day_of_month(year, m)));
+	 temp = temp + 1, m++)
+      ;
+    sumres2 = temp;
+  }
+  month = 1
+    + sumres2;
+  day = date - (bahai_to_rd(year, month, 1) - 1);
+  if (rmonth)
+    *rmonth = month;
+  if (rday)
+    *rday = day;
+  if (ryear)
+    *ryear = year;
+}
+
+void
+bahai_year_to_vahid(int year, int *rkull, int *rvahid, int *ryear)
+{
+  int r;
+  *rkull = quotient(year - 1,  361) + 1;
+  r = mod(year - 1, 361);
+  *rvahid = quotient(r, 19) + 1;
+  *ryear = mod(r, 19) + 1;
+}
+
+int
+bahai_vahid_to_year(int kull, int vahid, int year)
+{
+  return (kull - 1) * 361 + (vahid - 1) * 19 + year;
+}
+
+int
+indian_national_leap_year(int year)
+{
+    return gregorian_leap_year(year + 78);
+}
+
+int
+indian_national_last_day_of_month(int year, int month)
+{
+  if (month == 1 && indian_national_leap_year(year))
+      return 31;
+  switch (month) {
+  case  1: return 30;
+  case  2: return 31;
+  case  3: return 31;
+  case  4: return 31;
+  case  5: return 31;
+  case  6: return 31;
+  case  7: return 30;
+  case  8: return 30;
+  case  9: return 30;
+  case 10: return 30;
+  case 11: return 30;
+  case 12: return 30;
+  default: return 0;
+  }
+}
+
+int
+indian_national_to_rd(int year, int month, int day)
+{
+  int sumres;
+
+  {
+    int temp, m;
+    for (temp = 0, m = 1;
+	 (m < month);
+	 temp = temp + indian_national_last_day_of_month(year, m), m++)
+      ;
+    sumres = temp;
+  }
+  return day
+    + sumres
+    + gregorian_to_rd(year + 78, 1, 1) + 79;
+}
+
+#define BEFORE_INDIAN_NATIONAL (28205 - 1)
+
+void
+rd_to_indian_national(int date, int *ryear, int *rmonth, int *rday)
+{
+  int approx, month, day, year;
+  int sumres1, sumres2;
+
+  approx = approximate(date - BEFORE_INDIAN_NATIONAL, 366);
+  {
+    int temp, y;
+    for (temp = -1, y = approx;
+	 (date >= indian_national_to_rd(y, 1, 1));
+	 temp = temp + 1, y++)
+      ;
+    sumres1 = temp;
+  }
+  year = approx
+    + sumres1;
+  {
+    int temp, m;
+    for (temp = 0, m = 1;
+	 (date > indian_national_to_rd
+	  (year, m, indian_national_last_day_of_month(year, m)));
+	 temp = temp + 1, m++)
+      ;
+    sumres2 = temp;
+  }
+  month = 1
+    + sumres2;
+  day = date - (indian_national_to_rd(year, month, 1) - 1);
+  if (rmonth)
+    *rmonth = month;
+  if (rday)
+    *rday = day;
+  if (ryear)
+    *ryear = year;
+}
+
+int
+bengali_leap_year(int year)
+{
+    return gregorian_leap_year(year + 594);
+}
+
+int
+bengali_last_day_of_month(int year, int month)
+{
+  if (month == 11 && bengali_leap_year(year))
+      return 31;
+  switch (month) {
+  case  1: return 31;
+  case  2: return 31;
+  case  3: return 31;
+  case  4: return 31;
+  case  5: return 31;
+  case  6: return 30;
+  case  7: return 30;
+  case  8: return 30;
+  case  9: return 30;
+  case 10: return 30;
+  case 11: return 30;
+  case 12: return 30;
+  default: return 0;
+  }
+}
+
+int
+bengali_to_rd(int year, int month, int day)
+{
+  int sumres;
+
+  {
+    int temp, m;
+    for (temp = 0, m = 1;
+	 (m < month);
+	 temp = temp + bengali_last_day_of_month(year, m), m++)
+      ;
+    sumres = temp;
+  }
+  return day
+    + sumres
+    + gregorian_to_rd(year + 593, 4, 13);
+}
+
+#define BEFORE_BENGALI (216693 - 1)
+
+void
+rd_to_bengali(int date, int *ryear, int *rmonth, int *rday)
+{
+  int approx, month, day, year;
+  int sumres1, sumres2;
+
+  approx = approximate(date - BEFORE_BENGALI, 366);
+  {
+    int temp, y;
+    for (temp = 0, y = approx;
+	 (date >= bengali_to_rd(1 + y, 1, 1));
+	 temp = temp + 1, y++)
+      ;
+    sumres1 = temp;
+  }
+  year = approx
+    + sumres1;
+  {
+    int temp, m;
+    for (temp = 0, m = 1;
+	 (date > bengali_to_rd
+	  (year, m, bengali_last_day_of_month(year, m)));
+	 temp = temp + 1, m++)
+      ;
+    sumres2 = temp;
+  }
+  month = 1
+    + sumres2;
+  day = date - (bengali_to_rd(year, month, 1) - 1);
+  if (rmonth)
+    *rmonth = month;
+  if (rday)
+    *rday = day;
+  if (ryear)
+    *ryear = year;
+}
+
+int
+nanakshahi_leap_year(int year)
+{
+    return gregorian_leap_year(year + 1469);
+}
+
+int
+nanakshahi_last_day_of_month(int year, int month)
+{
+  if (month == 12 && nanakshahi_leap_year(year))
+      return 31;
+  switch (month) {
+  case  1: return 31;
+  case  2: return 31;
+  case  3: return 31;
+  case  4: return 31;
+  case  5: return 31;
+  case  6: return 30;
+  case  7: return 30;
+  case  8: return 30;
+  case  9: return 30;
+  case 10: return 30;
+  case 11: return 30;
+  case 12: return 30;
+  default: return 0;
+  }
+}
+
+int
+nanakshahi_to_rd(int year, int month, int day)
+{
+  int sumres;
+
+  {
+    int temp, m;
+    for (temp = 0, m = 1;
+	 (m < month);
+	 temp = temp + nanakshahi_last_day_of_month(year, m), m++)
+      ;
+    sumres = temp;
+  }
+  return day
+    + sumres
+    + gregorian_to_rd(year + 1468, 3, 13);
+}
+
+#define BEFORE_NANAKSHAHI (536258 - 1)
+
+void
+rd_to_nanakshahi(int date, int *ryear, int *rmonth, int *rday)
+{
+  int approx, month, day, year;
+  int sumres1, sumres2;
+
+  approx = approximate(date - BEFORE_NANAKSHAHI, 366);
+  {
+    int temp, y;
+    for (temp = 0, y = approx;
+	 (date >= nanakshahi_to_rd(1 + y, 1, 1));
+	 temp = temp + 1, y++)
+      ;
+    sumres1 = temp;
+  }
+  year = approx
+    + sumres1;
+  {
+    int temp, m;
+    for (temp = 0, m = 1;
+	 (date > nanakshahi_to_rd
+	  (year, m, nanakshahi_last_day_of_month(year, m)));
+	 temp = temp + 1, m++)
+      ;
+    sumres2 = temp;
+  }
+  month = 1
+    + sumres2;
+  day = date - (nanakshahi_to_rd(year, month, 1) - 1);
+  if (rmonth)
+    *rmonth = month;
+  if (rday)
+    *rday = day;
+  if (ryear)
+    *ryear = year;
+}
+
+#ifndef PLUTO
+int
+ordinal_to_rd(int year, int day)
+{
+  return gregorian_to_rd(year, 1, 1) + day - 1;
+}
+
+void
+rd_to_ordinal(int date, int *ryear, int *rday)
+{
+  int month, day, year, date2;
+
+  rd_to_gregorian(date, &year, &month, &day);
+  date2 = gregorian_to_rd(year, 1, 1);
+  if (rday)
+    *rday = date - date2 + 1;
+  if (ryear)
+    *ryear = year;
+}
+#endif
+
+int
+jd_to_rd(int date) /* chronological */
+{
+  return date - 1721425;
+}
+
+int
+rd_to_jd(int date) /* chronological */
+{
+  return date + 1721425;
+}
+
+#ifndef PLUTO
+int
+mjd_to_rd(int date) /* chronological */
+{
+  return date + 678576;
+}
+
+int
+rd_to_mjd(int date) /* chronological */
+{
+  return date - 678576;
+}
+
+int
+mjd_to_jd(int date) /* chronological */
+{
+  return date + 2400001;
+}
+
+int
+jd_to_mjd(int date) /* chronological */
+{
+  return date - 2400001;
+}
+
+int
+ld_to_rd(int date)
+{
+  return date + 577735;
+}
+
+int
+rd_to_ld(int date)
+{
+  return date - 577735;
+}
+
+int
+rd_to_day_of_week(int date)
+{
+  return mod(date, 7);
+}
+#endif /* not PLUTO */
+
+
+/*
+Local Variables:
+c-basic-offset: 2
+End:
+*/