/* glpmpl05.c */

/***********************************************************************

*  This code is part of GLPK (GNU Linear Programming Kit).

*

*  Authors: Andrew Makhorin <mao@gnu.org>

*           Heinrich Schuchardt <xypron.glpk@gmx.de>

*

*  Copyright (C) 2000, 2001, 2002, 2003, 2004, 2005, 2006, 2007, 2008,

*  2009, 2010 Andrew Makhorin, Department for Applied Informatics,

*  Moscow Aviation Institute, Moscow, Russia. All rights reserved.

*  E-mail: <mao@gnu.org>.

*

*  GLPK is free software: you can redistribute it and/or modify it

*  under the terms of the GNU General Public License as published by

*  the Free Software Foundation, either version 3 of the License, or

*  (at your option) any later version.

*

*  GLPK is distributed in the hope that it will be useful, but WITHOUT

*  ANY WARRANTY; without even the implied warranty of MERCHANTABILITY

*  or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public

*  License for more details.

*

*  You should have received a copy of the GNU General Public License

*  along with GLPK. If not, see <http://www.gnu.org/licenses/>.

***********************************************************************/

#define _GLPSTD_STDIO

#define _GLPSTD_TIME

#include "glpmpl.h"

double fn_gmtime(MPL *mpl)

{     /* obtain the current calendar time (UTC) */

      time_t timer;

      struct tm *tm;

      int j;

      time(&timer);

      if (timer == (time_t)(-1))

err:     error(mpl, "gmtime(); unable to obtain current calendar time");

      tm = gmtime(&timer);

      if (tm == NULL) goto err;

      j = jday(tm->tm_mday, tm->tm_mon + 1, 1900 + tm->tm_year);

      if (j < 0) goto err;

      return (((double)(j - jday(1, 1, 1970)) * 24.0 +

         (double)tm->tm_hour) * 60.0 + (double)tm->tm_min) * 60.0 +

         (double)tm->tm_sec;

}

static char *week[] = { "Monday", "Tuesday", "Wednesday", "Thursday",

      "Friday", "Saturday", "Sunday" };

static char *moon[] = { "January", "February", "March", "April", "May",

      "June", "July", "August", "September", "October", "November",

      "December" };

static void error1(MPL *mpl, const char *str, const char *s,

      const char *fmt, const char *f, const char *msg)

{     xprintf("Input string passed to str2time:\n");

      xprintf("%s\n", str);

      xprintf("%*s\n", (s - str) + 1, "^");

      xprintf("Format string passed to str2time:\n");

      xprintf("%s\n", fmt);

      xprintf("%*s\n", (f - fmt) + 1, "^");

      error(mpl, "%s", msg);

      /* no return */

}

double fn_str2time(MPL *mpl, const char *str, const char *fmt)

{     /* convert character string to the calendar time */

      int j, year, month, day, hh, mm, ss, zone;

      const char *s, *f;

      year = month = day = hh = mm = ss = -1, zone = INT_MAX;

      s = str;

      for (f = fmt; *f != '\0'; f++)

      {  if (*f == '%')

         {  f++;

            if (*f == 'b' || *f == 'h')

            {  /* the abbreviated month name */

               int k;

               char *name;

               if (month >= 0)

                  error1(mpl, str, s, fmt, f, "month multiply specified"

                     );

               while (*s == ' ') s++;

               for (month = 1; month <= 12; month++)

               {  name = moon[month-1];

                  for (k = 0; k <= 2; k++)

                  {  if (toupper((unsigned char)s[k]) !=

                         toupper((unsigned char)name[k])) goto next;

                  }

                  s += 3;

                  for (k = 3; name[k] != '\0'; k++)

                  {  if (toupper((unsigned char)*s) !=

                         toupper((unsigned char)name[k])) break;

                     s++;

                  }

                  break;

next:             ;

               }

               if (month > 12)

                  error1(mpl, str, s, fmt, f, "abbreviated month name m"

                     "issing or invalid");

            }

            else if (*f == 'd')

            {  /* the day of the month as a decimal number (01..31) */

               if (day >= 0)

                  error1(mpl, str, s, fmt, f, "day multiply specified");

               while (*s == ' ') s++;

               if (!('0' <= *s && *s <= '9'))

                  error1(mpl, str, s, fmt, f, "day missing or invalid");

               day = (*s++) - '0';

               if ('0' <= *s && *s <= '9')

                  day = 10 * day + ((*s++) - '0');

               if (!(1 <= day && day <= 31))

                  error1(mpl, str, s, fmt, f, "day out of range");

            }

            else if (*f == 'H')

            {  /* the hour as a decimal number, using a 24-hour clock

                  (00..23) */

               if (hh >= 0)

                  error1(mpl, str, s, fmt, f, "hour multiply specified")

                     ;

               while (*s == ' ') s++;

               if (!('0' <= *s && *s <= '9'))

                  error1(mpl, str, s, fmt, f, "hour missing or invalid")

                     ;

               hh = (*s++) - '0';

               if ('0' <= *s && *s <= '9')

                  hh = 10 * hh + ((*s++) - '0');

               if (!(0 <= hh && hh <= 23))

                  error1(mpl, str, s, fmt, f, "hour out of range");

            }

            else if (*f == 'm')

            {  /* the month as a decimal number (01..12) */

               if (month >= 0)

                  error1(mpl, str, s, fmt, f, "month multiply specified"

                     );

               while (*s == ' ') s++;

               if (!('0' <= *s && *s <= '9'))

                  error1(mpl, str, s, fmt, f, "month missing or invalid"

                     );

               month = (*s++) - '0';

               if ('0' <= *s && *s <= '9')

                  month = 10 * month + ((*s++) - '0');

               if (!(1 <= month && month <= 12))

                  error1(mpl, str, s, fmt, f, "month out of range");

            }

            else if (*f == 'M')

            {  /* the minute as a decimal number (00..59) */

               if (mm >= 0)

                  error1(mpl, str, s, fmt, f, "minute multiply specifie"

                     "d");

               while (*s == ' ') s++;

               if (!('0' <= *s && *s <= '9'))

                  error1(mpl, str, s, fmt, f, "minute missing or invali"

                     "d");

               mm = (*s++) - '0';

               if ('0' <= *s && *s <= '9')

                  mm = 10 * mm + ((*s++) - '0');

               if (!(0 <= mm && mm <= 59))

                  error1(mpl, str, s, fmt, f, "minute out of range");

            }

            else if (*f == 'S')

            {  /* the second as a decimal number (00..60) */

               if (ss >= 0)

                  error1(mpl, str, s, fmt, f, "second multiply specifie"

                     "d");

               while (*s == ' ') s++;

               if (!('0' <= *s && *s <= '9'))

                  error1(mpl, str, s, fmt, f, "second missing or invali"

                     "d");

               ss = (*s++) - '0';

               if ('0' <= *s && *s <= '9')

                  ss = 10 * ss + ((*s++) - '0');

               if (!(0 <= ss && ss <= 60))

                  error1(mpl, str, s, fmt, f, "second out of range");

            }

            else if (*f == 'y')

            {  /* the year without a century as a decimal number

                  (00..99); the values 00 to 68 mean the years 2000 to

                  2068 while the values 69 to 99 mean the years 1969 to

                  1999 */

               if (year >= 0)

                  error1(mpl, str, s, fmt, f, "year multiply specified")

                     ;

               while (*s == ' ') s++;

               if (!('0' <= *s && *s <= '9'))

                  error1(mpl, str, s, fmt, f, "year missing or invalid")

                     ;

               year = (*s++) - '0';

               if ('0' <= *s && *s <= '9')

                  year = 10 * year + ((*s++) - '0');

               year += (year >= 69 ? 1900 : 2000);

            }

            else if (*f == 'Y')

            {  /* the year as a decimal number, using the Gregorian

                  calendar */

               if (year >= 0)

                  error1(mpl, str, s, fmt, f, "year multiply specified")

                     ;

               while (*s == ' ') s++;

               if (!('0' <= *s && *s <= '9'))

                  error1(mpl, str, s, fmt, f, "year missing or invalid")

                     ;

               year = 0;

               for (j = 1; j <= 4; j++)

               {  if (!('0' <= *s && *s <= '9')) break;

                  year = 10 * year + ((*s++) - '0');

               }

               if (!(1 <= year && year <= 4000))

                  error1(mpl, str, s, fmt, f, "year out of range");

            }

            else if (*f == 'z')

            {  /* time zone offset in the form zhhmm */

               int z, hh, mm;

               if (zone != INT_MAX)

                  error1(mpl, str, s, fmt, f, "time zone offset multipl"

                     "y specified");

               while (*s == ' ') s++;

               if (*s == 'Z')

               {  z = hh = mm = 0, s++;

                  goto skip;

               }

               if (*s == '+')

                  z = +1, s++;

               else if (*s == '-')

                  z = -1, s++;

               else

                  error1(mpl, str, s, fmt, f, "time zone offset sign mi"

                     "ssing");

               hh = 0;

               for (j = 1; j <= 2; j++)

               {  if (!('0' <= *s && *s <= '9'))

err1:                error1(mpl, str, s, fmt, f, "time zone offset valu"

                        "e incomplete or invalid");

                  hh = 10 * hh + ((*s++) - '0');

               }

               if (hh > 23)

err2:             error1(mpl, str, s, fmt, f, "time zone offset value o"

                     "ut of range");

               if (*s == ':')

               {  s++;

                  if (!('0' <= *s && *s <= '9')) goto err1;

               }

               mm = 0;

               if (!('0' <= *s && *s <= '9')) goto skip;

               for (j = 1; j <= 2; j++)

               {  if (!('0' <= *s && *s <= '9')) goto err1;

                  mm = 10 * mm + ((*s++) - '0');

               }

               if (mm > 59) goto err2;

skip:          zone = z * (60 * hh + mm);

            }

            else if (*f == '%')

            {  /* literal % character */

               goto test;

            }

            else

               error1(mpl, str, s, fmt, f, "invalid conversion specifie"

                  "r");

         }

         else if (*f == ' ')

            ;

         else

test:    {  /* check a matching character in the input string */

            if (*s != *f)

               error1(mpl, str, s, fmt, f, "character mismatch");

            s++;

         }

      }

      if (year < 0) year = 1970;

      if (month < 0) month = 1;

      if (day < 0) day = 1;

      if (hh < 0) hh = 0;

      if (mm < 0) mm = 0;

      if (ss < 0) ss = 0;

      if (zone == INT_MAX) zone = 0;

      j = jday(day, month, year);

      xassert(j >= 0);

      return (((double)(j - jday(1, 1, 1970)) * 24.0 + (double)hh) *

         60.0 + (double)mm) * 60.0 + (double)ss - 60.0 * (double)zone;

}

static void error2(MPL *mpl, const char *fmt, const char *f,

      const char *msg)

{     xprintf("Format string passed to time2str:\n");

      xprintf("%s\n", fmt);

      xprintf("%*s\n", (f - fmt) + 1, "^");

      error(mpl, "%s", msg);

      /* no return */

}

static int weekday(int j)

{     /* determine weekday number (1 = Mon, ..., 7 = Sun) */

      return (j + jday(1, 1, 1970)) % 7 + 1;

}

static int firstday(int year)

{     /* determine the first day of the first week for a specified year

         according to ISO 8601 */

      int j;

      /* if 1 January is Monday, Tuesday, Wednesday or Thursday, it is

         in week 01; if 1 January is Friday, Saturday or Sunday, it is

         in week 52 or 53 of the previous year */

      j = jday(1, 1, year) - jday(1, 1, 1970);

      switch (weekday(j))

      {  case 1: /* 1 Jan is Mon */ j += 0; break;

         case 2: /* 1 Jan is Tue */ j -= 1; break;

         case 3: /* 1 Jan is Wed */ j -= 2; break;

         case 4: /* 1 Jan is Thu */ j -= 3; break;

         case 5: /* 1 Jan is Fri */ j += 3; break;

         case 6: /* 1 Jan is Sat */ j += 2; break;

         case 7: /* 1 Jan is Sun */ j += 1; break;

         default: xassert(j != j);

      }

      /* the first day of the week must be Monday */

      xassert(weekday(j) == 1);

      return j;

}

void fn_time2str(MPL *mpl, char *str, double t, const char *fmt)

{     /* convert the calendar time to character string */

      int j, year, month, day, hh, mm, ss, len;

      double temp;

      const char *f;

      char buf[MAX_LENGTH+1];

      if (!(-62135596800.0 <= t && t <= 64092211199.0))

         error(mpl, "time2str(%.*g,...); argument out of range",

            DBL_DIG, t);

      t = floor(t + 0.5);

      temp = fabs(t) / 86400.0;

      j = (int)floor(temp);

      if (t < 0.0)

      {  if (temp == floor(temp))

            j = - j;

         else

            j = - (j + 1);

      }

      xassert(jdate(j + jday(1, 1, 1970), &day, &month, &year) == 0);

      ss = (int)(t - 86400.0 * (double)j);

      xassert(0 <= ss && ss < 86400);

      mm = ss / 60, ss %= 60;

      hh = mm / 60, mm %= 60;

      len = 0;

      for (f = fmt; *f != '\0'; f++)

      {  if (*f == '%')

         {  f++;

            if (*f == 'a')

            {  /* the abbreviated weekday name */

               memcpy(buf, week[weekday(j)-1], 3), buf[3] = '\0';

            }

            else if (*f == 'A')

            {  /* the full weekday name */

               strcpy(buf, week[weekday(j)-1]);

            }

            else if (*f == 'b' || *f == 'h')

            {  /* the abbreviated month name */

               memcpy(buf, moon[month-1], 3), buf[3] = '\0';

            }

            else if (*f == 'B')

            {  /* the full month name */

               strcpy(buf, moon[month-1]);

            }

            else if (*f == 'C')

            {  /* the century of the year */

               sprintf(buf, "%02d", year / 100);

            }

            else if (*f == 'd')

            {  /* the day of the month as a decimal number (01..31) */

               sprintf(buf, "%02d", day);

            }

            else if (*f == 'D')

            {  /* the date using the format %m/%d/%y */

               sprintf(buf, "%02d/%02d/%02d", month, day, year % 100);

            }

            else if (*f == 'e')

            {  /* the day of the month like with %d, but padded with

                  blank (1..31) */

               sprintf(buf, "%2d", day);

            }

            else if (*f == 'F')

            {  /* the date using the format %Y-%m-%d */

               sprintf(buf, "%04d-%02d-%02d", year, month, day);

            }

            else if (*f == 'g')

            {  /* the year corresponding to the ISO week number, but

                  without the century (range 00 through 99); this has

                  the same format and value as %y, except that if the

                  ISO week number (see %V) belongs to the previous or

                  next year, that year is used instead */

               int iso;

               if (j < firstday(year))

                  iso = year - 1;

               else if (j < firstday(year + 1))

                  iso = year;

               else

                  iso = year + 1;

               sprintf(buf, "%02d", iso % 100);

            }

            else if (*f == 'G')

            {  /* the year corresponding to the ISO week number; this

                  has the same format and value as %Y, excepth that if

                  the ISO week number (see %V) belongs to the previous

                  or next year, that year is used instead */

               int iso;

               if (j < firstday(year))

                  iso = year - 1;

               else if (j < firstday(year + 1))

                  iso = year;

               else

                  iso = year + 1;

               sprintf(buf, "%04d", iso);

            }

            else if (*f == 'H')

            {  /* the hour as a decimal number, using a 24-hour clock

                  (00..23) */

               sprintf(buf, "%02d", hh);

            }

            else if (*f == 'I')

            {  /* the hour as a decimal number, using a 12-hour clock

                  (01..12) */

               sprintf(buf, "%02d",

                  hh == 0 ? 12 : hh <= 12 ? hh : hh - 12);

            }

            else if (*f == 'j')

            {  /* the day of the year as a decimal number (001..366) */

               sprintf(buf, "%03d",

                  jday(day, month, year) - jday(1, 1, year) + 1);

            }

            else if (*f == 'k')

            {  /* the hour as a decimal number, using a 24-hour clock

                  like %H, but padded with blank (0..23) */

               sprintf(buf, "%2d", hh);

            }

            else if (*f == 'l')

            {  /* the hour as a decimal number, using a 12-hour clock

                  like %I, but padded with blank (1..12) */

               sprintf(buf, "%2d",

                  hh == 0 ? 12 : hh <= 12 ? hh : hh - 12);

            }

            else if (*f == 'm')

            {  /* the month as a decimal number (01..12) */

               sprintf(buf, "%02d", month);

            }

            else if (*f == 'M')

            {  /* the minute as a decimal number (00..59) */

               sprintf(buf, "%02d", mm);

            }

            else if (*f == 'p')

            {  /* either AM or PM, according to the given time value;

                  noon is treated as PM and midnight as AM */

               strcpy(buf, hh <= 11 ? "AM" : "PM");

            }

            else if (*f == 'P')

            {  /* either am or pm, according to the given time value;

                  noon is treated as pm and midnight as am */

               strcpy(buf, hh <= 11 ? "am" : "pm");

            }

            else if (*f == 'r')

            {  /* the calendar time using the format %I:%M:%S %p */

               sprintf(buf, "%02d:%02d:%02d %s",

                  hh == 0 ? 12 : hh <= 12 ? hh : hh - 12,

                  mm, ss, hh <= 11 ? "AM" : "PM");

            }

            else if (*f == 'R')

            {  /* the hour and minute using the format %H:%M */

               sprintf(buf, "%02d:%02d", hh, mm);

            }

            else if (*f == 'S')

            {  /* the second as a decimal number (00..59) */

               sprintf(buf, "%02d", ss);

            }

            else if (*f == 'T')

            {  /* the time of day using the format %H:%M:%S */

               sprintf(buf, "%02d:%02d:%02d", hh, mm, ss);

            }

            else if (*f == 'u')

            {  /* the day of the week as a decimal number (1..7),

                  Monday being 1 */

               sprintf(buf, "%d", weekday(j));

            }

            else if (*f == 'U')

            {  /* the week number of the current year as a decimal

                  number (range 00 through 53), starting with the first

                  Sunday as the first day of the first week; days

                  preceding the first Sunday in the year are considered

                  to be in week 00 */

#if 1 /* 09/I-2009 */

#undef sun

/* causes compilation error in SunOS */

#endif

               int sun;

               /* sun = the first Sunday of the year */

               sun = jday(1, 1, year) - jday(1, 1, 1970);

               sun += (7 - weekday(sun));

               sprintf(buf, "%02d", (j + 7 - sun) / 7);

            }

            else if (*f == 'V')

            {  /* the ISO week number as a decimal number (range 01

                  through 53); ISO weeks start with Monday and end with

                  Sunday; week 01 of a year is the first week which has

                  the majority of its days in that year; week 01 of

                  a year can contain days from the previous year; the

                  week before week 01 of a year is the last week (52 or

                  53) of the previous year even if it contains days

                  from the new year */

               int iso;

               if (j < firstday(year))

                  iso = j - firstday(year - 1);

               else if (j < firstday(year + 1))

                  iso = j - firstday(year);

               else

                  iso = j - firstday(year + 1);

               sprintf(buf, "%02d", iso / 7 + 1);

            }

            else if (*f == 'w')

            {  /* the day of the week as a decimal number (0..6),

                  Sunday being 0 */

               sprintf(buf, "%d", weekday(j) % 7);

            }

            else if (*f == 'W')

            {  /* the week number of the current year as a decimal

                  number (range 00 through 53), starting with the first

                  Monday as the first day of the first week; days

                  preceding the first Monday in the year are considered

                  to be in week 00 */

               int mon;

               /* mon = the first Monday of the year */

               mon = jday(1, 1, year) - jday(1, 1, 1970);

               mon += (8 - weekday(mon)) % 7;

               sprintf(buf, "%02d", (j + 7 - mon) / 7);

            }

            else if (*f == 'y')

            {  /* the year without a century as a decimal number

                  (00..99) */

               sprintf(buf, "%02d", year % 100);

            }

            else if (*f == 'Y')

            {  /* the year as a decimal number, using the Gregorian

                  calendar */

               sprintf(buf, "%04d", year);

            }

            else if (*f == '%')

            {  /* a literal % character */

               buf[0] = '%', buf[1] = '\0';

            }

            else

               error2(mpl, fmt, f, "invalid conversion specifier");

         }

         else

            buf[0] = *f, buf[1] = '\0';

         if (len + strlen(buf) > MAX_LENGTH)

            error(mpl, "time2str; output string length exceeds %d chara"

               "cters", MAX_LENGTH);

         memcpy(str+len, buf, strlen(buf));

         len += strlen(buf);

      }

      str[len] = '\0';

      return;

}

/* eof */