Need to make sure m4/Makefile.am exists before running gnulib-tool.
[pspp] / src / calendar.c
1 #include <config.h>
2 #include "calendar.h"
3 #include <assert.h>
4 #include <stdbool.h>
5 #include "settings.h"
6 #include "val.h"
7
8 #include "gettext.h"
9 #define _(msgid) gettext (msgid)
10
11 /* 14 Oct 1582. */
12 #define EPOCH (-577734)
13
14 /* Calculates and returns floor(a/b) for integer b > 0. */
15 static int
16 floor_div (int a, int b) 
17 {
18   assert (b > 0);
19   return (a >= 0 ? a : a - b + 1) / b;
20 }
21
22 /* Calculates floor(a/b) and the corresponding remainder and
23    stores them into *Q and *R. */
24 static void
25 floor_divmod (int a, int b, int *q, int *r) 
26 {
27   *q = floor_div (a, b);
28   *r = a - b * *q;
29 }
30
31 /* Returns true if Y is a leap year, false otherwise. */
32 static bool
33 is_leap_year (int y) 
34 {
35   return y % 4 == 0 && (y % 100 != 0 || y % 400 == 0);
36 }
37
38 static int
39 raw_gregorian_to_offset (int y, int m, int d) 
40 {
41   return (EPOCH - 1
42           + 365 * (y - 1)
43           + floor_div (y - 1, 4)
44           - floor_div (y - 1, 100)
45           + floor_div (y - 1, 400)
46           + floor_div (367 * m - 362, 12)
47           + (m <= 2 ? 0 : (m >= 2 && is_leap_year (y) ? -1 : -2))
48           + d);
49 }
50
51 /* Returns the number of days from 14 Oct 1582 to (Y,M,D) in the
52    Gregorian calendar.  Returns SYSMIS for dates before 14 Oct
53    1582. */
54 double
55 calendar_gregorian_to_offset (int y, int m, int d,
56                               calendar_error_func *error, void *aux)
57 {
58   /* Normalize year. */
59   if (y >= 0 && y < 100) 
60     {
61       int epoch = get_epoch ();
62       int century = epoch / 100 + (y < epoch % 100);
63       y += century * 100;
64     }
65
66   /* Normalize month. */
67   if (m < 1 || m > 12) 
68     {
69       if (m == 0) 
70         {
71           y--;
72           m = 12;
73         }
74       else if (m == 13) 
75         {
76           y++;
77           m = 1;
78         }
79       else
80         {
81           error (aux, _("Month %d is not in acceptable range of 0 to 13."), m);
82           return SYSMIS;
83         }
84     }
85
86   /* Normalize day. */
87   if (d < 0 || d > 31) 
88     {
89       error (aux, _("Day %d is not in acceptable range of 0 to 31."), d);
90       return SYSMIS;
91     }
92
93   /* Validate date. */
94   if (y < 1582 || (y == 1582 && (m < 10 || (m == 10 && d < 15)))) 
95     {
96       error (aux, _("Date %04d-%d-%d is before the earliest acceptable "
97                     "date of 1582-10-15."), y, m, d);
98       return SYSMIS;
99     }
100
101   /* Calculate offset. */
102   return raw_gregorian_to_offset (y, m, d);
103 }
104
105 /* Returns the number of days in the given YEAR from January 1 up
106    to (but not including) the first day of MONTH. */
107 static int
108 cum_month_days (int year, int month) 
109 {
110   static const int cum_month_days[12] = 
111     {
112       0,
113       31, /* Jan */
114       31 + 28, /* Feb */
115       31 + 28 + 31, /* Mar */
116       31 + 28 + 31 + 30, /* Apr */
117       31 + 28 + 31 + 30 + 31, /* May */
118       31 + 28 + 31 + 30 + 31 + 30, /* Jun */
119       31 + 28 + 31 + 30 + 31 + 30 + 31, /* Jul */
120       31 + 28 + 31 + 30 + 31 + 30 + 31 + 31, /* Aug */
121       31 + 28 + 31 + 30 + 31 + 30 + 31 + 31 + 30, /* Sep */
122       31 + 28 + 31 + 30 + 31 + 30 + 31 + 31 + 30 + 31, /* Oct */
123       31 + 28 + 31 + 30 + 31 + 30 + 31 + 31 + 30 + 31 + 30, /* Nov */
124     };
125
126   assert (month >= 1 && month <= 12);
127   return cum_month_days[month - 1] + (month >= 3 && is_leap_year (year));
128 }
129
130 /* Takes a count of days from 14 Oct 1582 and returns the
131    Gregorian calendar year it is in.  Dates both before and after
132    the epoch are supported. */
133 int
134 calendar_offset_to_year (int ofs) 
135 {
136   int d0;
137   int n400, d1;
138   int n100, d2;
139   int n4, d3;
140   int n1;
141   int y;
142
143   d0 = ofs - EPOCH;
144   floor_divmod (d0, 365 * 400 + 100 - 3, &n400, &d1);
145   floor_divmod (d1, 365 * 100 + 25 - 1, &n100, &d2);
146   floor_divmod (d2, 365 * 4 + 1, &n4, &d3);
147   n1 = floor_div (d3, 365);
148   y = 400 * n400 + 100 * n100 + 4 * n4 + n1;
149   if (n100 != 4 && n1 != 4)
150     y++;
151
152   return y;
153 }
154
155 /* Takes a count of days from 14 Oct 1582 and translates it into
156    a Gregorian calendar date in (*Y,*M,*D).  Dates both before
157    and after the epoch are supported. */
158 void
159 calendar_offset_to_gregorian (int ofs, int *y, int *m, int *d)
160 {
161   int year = *y = calendar_offset_to_year (ofs);
162   int january1 = raw_gregorian_to_offset (year, 1, 1);
163   int yday = ofs - january1 + 1;
164   int march1 = january1 + cum_month_days (year, 3);
165   int correction = ofs < march1 ? 0 : (is_leap_year (year) ? 1 : 2);
166   int month = *m = (12 * (yday - 1 + correction) + 373) / 367;
167   *d = yday - cum_month_days (year, month);
168 }
169
170 /* Takes a count of days from 14 Oct 1582 and returns the 1-based
171    year-relative day number, that is, the number of days from the
172    beginning of the year. */
173 int
174 calendar_offset_to_yday (int ofs)
175 {
176   int year = calendar_offset_to_year (ofs);
177   int january1 = raw_gregorian_to_offset (year, 1, 1);
178   int yday = ofs - january1 + 1;
179   return yday;
180 }
181
182 /* Takes a count of days from 14 Oct 1582 and returns the
183    corresponding weekday 1...7, with 1=Sunday. */
184 int
185 calendar_offset_to_wday (int ofs)
186 {
187   int wday = (ofs - EPOCH + 1) % 7 + 1;
188   if (wday <= 0)
189     wday += 7;
190   return wday;
191 }
192
193 /* Takes a count of days from 14 Oct 1582 and returns the month
194    it is in. */
195 int
196 calendar_offset_to_month (int ofs) 
197 {
198   int y, m, d;
199   calendar_offset_to_gregorian (ofs, &y, &m, &d);
200   return m;
201 }
202
203 /* Takes a count of days from 14 Oct 1582 and returns the
204    corresponding day of the month. */
205 int
206 calendar_offset_to_mday (int ofs) 
207 {
208   int y, m, d;
209   calendar_offset_to_gregorian (ofs, &y, &m, &d);
210   return d;
211 }