hjl / source-git / glibc

Forked from source-git/glibc 3 years ago
Clone

Blame time/mktime.c

Packit 6c4009
/* Convert a 'struct tm' to a time_t value.
Packit 6c4009
   Copyright (C) 1993-2018 Free Software Foundation, Inc.
Packit 6c4009
   This file is part of the GNU C Library.
Packit 6c4009
   Contributed by Paul Eggert <eggert@twinsun.com>.
Packit 6c4009
Packit 6c4009
   The GNU C Library is free software; you can redistribute it and/or
Packit 6c4009
   modify it under the terms of the GNU Lesser General Public
Packit 6c4009
   License as published by the Free Software Foundation; either
Packit 6c4009
   version 2.1 of the License, or (at your option) any later version.
Packit 6c4009
Packit 6c4009
   The GNU C Library is distributed in the hope that it will be useful,
Packit 6c4009
   but WITHOUT ANY WARRANTY; without even the implied warranty of
Packit 6c4009
   MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
Packit 6c4009
   Lesser General Public License for more details.
Packit 6c4009
Packit 6c4009
   You should have received a copy of the GNU Lesser General Public
Packit 6c4009
   License along with the GNU C Library; if not, see
Packit 6c4009
   <http://www.gnu.org/licenses/>.  */
Packit 6c4009
Packit 6c4009
/* Define this to have a standalone program to test this implementation of
Packit 6c4009
   mktime.  */
Packit 6c4009
/* #define DEBUG_MKTIME 1 */
Packit 6c4009
Packit 6c4009
#ifndef _LIBC
Packit 6c4009
# include <config.h>
Packit 6c4009
#endif
Packit 6c4009
Packit 6c4009
/* Assume that leap seconds are possible, unless told otherwise.
Packit 6c4009
   If the host has a 'zic' command with a '-L leapsecondfilename' option,
Packit 6c4009
   then it supports leap seconds; otherwise it probably doesn't.  */
Packit 6c4009
#ifndef LEAP_SECONDS_POSSIBLE
Packit 6c4009
# define LEAP_SECONDS_POSSIBLE 1
Packit 6c4009
#endif
Packit 6c4009
Packit 6c4009
#include <time.h>
Packit 6c4009
Packit 6c4009
#include <limits.h>
Packit 6c4009
Packit 6c4009
#include <string.h>		/* For the real memcpy prototype.  */
Packit 6c4009
Packit 6c4009
#if defined DEBUG_MKTIME && DEBUG_MKTIME
Packit 6c4009
# include <stdio.h>
Packit 6c4009
# include <stdlib.h>
Packit 6c4009
/* Make it work even if the system's libc has its own mktime routine.  */
Packit 6c4009
# undef mktime
Packit 6c4009
# define mktime my_mktime
Packit 6c4009
#endif /* DEBUG_MKTIME */
Packit 6c4009
Packit 6c4009
/* Some of the code in this file assumes that signed integer overflow
Packit 6c4009
   silently wraps around.  This assumption can't easily be programmed
Packit 6c4009
   around, nor can it be checked for portably at compile-time or
Packit 6c4009
   easily eliminated at run-time.
Packit 6c4009
Packit 6c4009
   Define WRAPV to 1 if the assumption is valid and if
Packit 6c4009
     #pragma GCC optimize ("wrapv")
Packit 6c4009
   does not trigger GCC bug 51793
Packit 6c4009
   <http://gcc.gnu.org/bugzilla/show_bug.cgi?id=51793>.
Packit 6c4009
   Otherwise, define it to 0; this forces the use of slower code that,
Packit 6c4009
   while not guaranteed by the C Standard, works on all production
Packit 6c4009
   platforms that we know about.  */
Packit 6c4009
#ifndef WRAPV
Packit 6c4009
# if (((__GNUC__ == 4 && 4 <= __GNUC_MINOR__) || 4 < __GNUC__) \
Packit 6c4009
      && defined __GLIBC__)
Packit 6c4009
#  pragma GCC optimize ("wrapv")
Packit 6c4009
#  define WRAPV 1
Packit 6c4009
# else
Packit 6c4009
#  define WRAPV 0
Packit 6c4009
# endif
Packit 6c4009
#endif
Packit 6c4009
Packit 6c4009
/* Verify a requirement at compile-time (unlike assert, which is runtime).  */
Packit 6c4009
#define verify(name, assertion) struct name { char a[(assertion) ? 1 : -1]; }
Packit 6c4009
Packit 6c4009
/* A signed type that is at least one bit wider than int.  */
Packit 6c4009
#if INT_MAX <= LONG_MAX / 2
Packit 6c4009
typedef long int long_int;
Packit 6c4009
#else
Packit 6c4009
typedef long long int long_int;
Packit 6c4009
#endif
Packit 6c4009
verify (long_int_is_wide_enough, INT_MAX == INT_MAX * (long_int) 2 / 2);
Packit 6c4009
Packit 6c4009
/* Shift A right by B bits portably, by dividing A by 2**B and
Packit 6c4009
   truncating towards minus infinity.  A and B should be free of side
Packit 6c4009
   effects, and B should be in the range 0 <= B <= INT_BITS - 2, where
Packit 6c4009
   INT_BITS is the number of useful bits in an int.  GNU code can
Packit 6c4009
   assume that INT_BITS is at least 32.
Packit 6c4009
Packit 6c4009
   ISO C99 says that A >> B is implementation-defined if A < 0.  Some
Packit 6c4009
   implementations (e.g., UNICOS 9.0 on a Cray Y-MP EL) don't shift
Packit 6c4009
   right in the usual way when A < 0, so SHR falls back on division if
Packit 6c4009
   ordinary A >> B doesn't seem to be the usual signed shift.  */
Packit 6c4009
#define SHR(a, b)                                               \
Packit 6c4009
  ((-1 >> 1 == -1                                               \
Packit 6c4009
    && (long_int) -1 >> 1 == -1                                 \
Packit 6c4009
    && ((time_t) -1 >> 1 == -1 || ! TYPE_SIGNED (time_t)))      \
Packit 6c4009
   ? (a) >> (b)                                                 \
Packit 6c4009
   : (a) / (1 << (b)) - ((a) % (1 << (b)) < 0))
Packit 6c4009
Packit 6c4009
/* The extra casts in the following macros work around compiler bugs,
Packit 6c4009
   e.g., in Cray C 5.0.3.0.  */
Packit 6c4009
Packit 6c4009
/* True if the arithmetic type T is an integer type.  bool counts as
Packit 6c4009
   an integer.  */
Packit 6c4009
#define TYPE_IS_INTEGER(t) ((t) 1.5 == 1)
Packit 6c4009
Packit 6c4009
/* True if negative values of the signed integer type T use two's
Packit 6c4009
   complement, or if T is an unsigned integer type.  */
Packit 6c4009
#define TYPE_TWOS_COMPLEMENT(t) ((t) ~ (t) 0 == (t) -1)
Packit 6c4009
Packit 6c4009
/* True if the arithmetic type T is signed.  */
Packit 6c4009
#define TYPE_SIGNED(t) (! ((t) 0 < (t) -1))
Packit 6c4009
Packit 6c4009
/* The maximum and minimum values for the integer type T.  These
Packit 6c4009
   macros have undefined behavior if T is signed and has padding bits.
Packit 6c4009
   If this is a problem for you, please let us know how to fix it for
Packit 6c4009
   your host.  */
Packit 6c4009
#define TYPE_MINIMUM(t) \
Packit 6c4009
  ((t) (! TYPE_SIGNED (t) \
Packit 6c4009
	? (t) 0 \
Packit 6c4009
	: ~ TYPE_MAXIMUM (t)))
Packit 6c4009
#define TYPE_MAXIMUM(t) \
Packit 6c4009
  ((t) (! TYPE_SIGNED (t) \
Packit 6c4009
	? (t) -1 \
Packit 6c4009
	: ((((t) 1 << (sizeof (t) * CHAR_BIT - 2)) - 1) * 2 + 1)))
Packit 6c4009
Packit 6c4009
#ifndef TIME_T_MIN
Packit 6c4009
# define TIME_T_MIN TYPE_MINIMUM (time_t)
Packit 6c4009
#endif
Packit 6c4009
#ifndef TIME_T_MAX
Packit 6c4009
# define TIME_T_MAX TYPE_MAXIMUM (time_t)
Packit 6c4009
#endif
Packit 6c4009
#define TIME_T_MIDPOINT (SHR (TIME_T_MIN + TIME_T_MAX, 1) + 1)
Packit 6c4009
Packit 6c4009
verify (time_t_is_integer, TYPE_IS_INTEGER (time_t));
Packit 6c4009
verify (twos_complement_arithmetic,
Packit 6c4009
	(TYPE_TWOS_COMPLEMENT (int)
Packit 6c4009
	 && TYPE_TWOS_COMPLEMENT (long_int)
Packit 6c4009
	 && TYPE_TWOS_COMPLEMENT (time_t)));
Packit 6c4009
Packit 6c4009
#define EPOCH_YEAR 1970
Packit 6c4009
#define TM_YEAR_BASE 1900
Packit 6c4009
verify (base_year_is_a_multiple_of_100, TM_YEAR_BASE % 100 == 0);
Packit 6c4009
Packit 6c4009
/* Return 1 if YEAR + TM_YEAR_BASE is a leap year.  */
Packit 6c4009
static int
Packit 6c4009
leapyear (long_int year)
Packit 6c4009
{
Packit 6c4009
  /* Don't add YEAR to TM_YEAR_BASE, as that might overflow.
Packit 6c4009
     Also, work even if YEAR is negative.  */
Packit 6c4009
  return
Packit 6c4009
    ((year & 3) == 0
Packit 6c4009
     && (year % 100 != 0
Packit 6c4009
	 || ((year / 100) & 3) == (- (TM_YEAR_BASE / 100) & 3)));
Packit 6c4009
}
Packit 6c4009
Packit 6c4009
/* How many days come before each month (0-12).  */
Packit 6c4009
#ifndef _LIBC
Packit 6c4009
static
Packit 6c4009
#endif
Packit 6c4009
const unsigned short int __mon_yday[2][13] =
Packit 6c4009
  {
Packit 6c4009
    /* Normal years.  */
Packit 6c4009
    { 0, 31, 59, 90, 120, 151, 181, 212, 243, 273, 304, 334, 365 },
Packit 6c4009
    /* Leap years.  */
Packit 6c4009
    { 0, 31, 60, 91, 121, 152, 182, 213, 244, 274, 305, 335, 366 }
Packit 6c4009
  };
Packit 6c4009
Packit 6c4009
Packit 6c4009
#ifndef _LIBC
Packit 6c4009
/* Portable standalone applications should supply a <time.h> that
Packit 6c4009
   declares a POSIX-compliant localtime_r, for the benefit of older
Packit 6c4009
   implementations that lack localtime_r or have a nonstandard one.
Packit 6c4009
   See the gnulib time_r module for one way to implement this.  */
Packit 6c4009
# undef __localtime_r
Packit 6c4009
# define __localtime_r localtime_r
Packit 6c4009
# define __mktime_internal mktime_internal
Packit 6c4009
# include "mktime-internal.h"
Packit 6c4009
#endif
Packit 6c4009
Packit 6c4009
/* Return 1 if the values A and B differ according to the rules for
Packit 6c4009
   tm_isdst: A and B differ if one is zero and the other positive.  */
Packit 6c4009
static int
Packit 6c4009
isdst_differ (int a, int b)
Packit 6c4009
{
Packit 6c4009
  return (!a != !b) && (0 <= a) && (0 <= b);
Packit 6c4009
}
Packit 6c4009
Packit 6c4009
/* Return an integer value measuring (YEAR1-YDAY1 HOUR1:MIN1:SEC1) -
Packit 6c4009
   (YEAR0-YDAY0 HOUR0:MIN0:SEC0) in seconds, assuming that the clocks
Packit 6c4009
   were not adjusted between the time stamps.
Packit 6c4009
Packit 6c4009
   The YEAR values uses the same numbering as TP->tm_year.  Values
Packit 6c4009
   need not be in the usual range.  However, YEAR1 must not be less
Packit 6c4009
   than 2 * INT_MIN or greater than 2 * INT_MAX.
Packit 6c4009
Packit 6c4009
   The result may overflow.  It is the caller's responsibility to
Packit 6c4009
   detect overflow.  */
Packit 6c4009
Packit 6c4009
static time_t
Packit 6c4009
ydhms_diff (long_int year1, long_int yday1, int hour1, int min1, int sec1,
Packit 6c4009
	    int year0, int yday0, int hour0, int min0, int sec0)
Packit 6c4009
{
Packit 6c4009
  verify (C99_integer_division, -1 / 2 == 0);
Packit 6c4009
Packit 6c4009
  /* Compute intervening leap days correctly even if year is negative.
Packit 6c4009
     Take care to avoid integer overflow here.  */
Packit 6c4009
  int a4 = SHR (year1, 2) + SHR (TM_YEAR_BASE, 2) - ! (year1 & 3);
Packit 6c4009
  int b4 = SHR (year0, 2) + SHR (TM_YEAR_BASE, 2) - ! (year0 & 3);
Packit 6c4009
  int a100 = a4 / 25 - (a4 % 25 < 0);
Packit 6c4009
  int b100 = b4 / 25 - (b4 % 25 < 0);
Packit 6c4009
  int a400 = SHR (a100, 2);
Packit 6c4009
  int b400 = SHR (b100, 2);
Packit 6c4009
  int intervening_leap_days = (a4 - b4) - (a100 - b100) + (a400 - b400);
Packit 6c4009
Packit 6c4009
  /* Compute the desired time in time_t precision.  Overflow might
Packit 6c4009
     occur here.  */
Packit 6c4009
  time_t tyear1 = year1;
Packit 6c4009
  time_t years = tyear1 - year0;
Packit 6c4009
  time_t days = 365 * years + yday1 - yday0 + intervening_leap_days;
Packit 6c4009
  time_t hours = 24 * days + hour1 - hour0;
Packit 6c4009
  time_t minutes = 60 * hours + min1 - min0;
Packit 6c4009
  time_t seconds = 60 * minutes + sec1 - sec0;
Packit 6c4009
  return seconds;
Packit 6c4009
}
Packit 6c4009
Packit 6c4009
/* Return the average of A and B, even if A + B would overflow.  */
Packit 6c4009
static time_t
Packit 6c4009
time_t_avg (time_t a, time_t b)
Packit 6c4009
{
Packit 6c4009
  return SHR (a, 1) + SHR (b, 1) + (a & b & 1);
Packit 6c4009
}
Packit 6c4009
Packit 6c4009
/* Return 1 if A + B does not overflow.  If time_t is unsigned and if
Packit 6c4009
   B's top bit is set, assume that the sum represents A - -B, and
Packit 6c4009
   return 1 if the subtraction does not wrap around.  */
Packit 6c4009
static int
Packit 6c4009
time_t_add_ok (time_t a, time_t b)
Packit 6c4009
{
Packit 6c4009
  if (! TYPE_SIGNED (time_t))
Packit 6c4009
    {
Packit 6c4009
      time_t sum = a + b;
Packit 6c4009
      return (sum < a) == (TIME_T_MIDPOINT <= b);
Packit 6c4009
    }
Packit 6c4009
  else if (WRAPV)
Packit 6c4009
    {
Packit 6c4009
      time_t sum = a + b;
Packit 6c4009
      return (sum < a) == (b < 0);
Packit 6c4009
    }
Packit 6c4009
  else
Packit 6c4009
    {
Packit 6c4009
      time_t avg = time_t_avg (a, b);
Packit 6c4009
      return TIME_T_MIN / 2 <= avg && avg <= TIME_T_MAX / 2;
Packit 6c4009
    }
Packit 6c4009
}
Packit 6c4009
Packit 6c4009
/* Return 1 if A + B does not overflow.  */
Packit 6c4009
static int
Packit 6c4009
time_t_int_add_ok (time_t a, int b)
Packit 6c4009
{
Packit 6c4009
  verify (int_no_wider_than_time_t, INT_MAX <= TIME_T_MAX);
Packit 6c4009
  if (WRAPV)
Packit 6c4009
    {
Packit 6c4009
      time_t sum = a + b;
Packit 6c4009
      return (sum < a) == (b < 0);
Packit 6c4009
    }
Packit 6c4009
  else
Packit 6c4009
    {
Packit 6c4009
      int a_odd = a & 1;
Packit 6c4009
      time_t avg = SHR (a, 1) + (SHR (b, 1) + (a_odd & b));
Packit 6c4009
      return TIME_T_MIN / 2 <= avg && avg <= TIME_T_MAX / 2;
Packit 6c4009
    }
Packit 6c4009
}
Packit 6c4009
Packit 6c4009
/* Return a time_t value corresponding to (YEAR-YDAY HOUR:MIN:SEC),
Packit 6c4009
   assuming that *T corresponds to *TP and that no clock adjustments
Packit 6c4009
   occurred between *TP and the desired time.
Packit 6c4009
   If TP is null, return a value not equal to *T; this avoids false matches.
Packit 6c4009
   If overflow occurs, yield the minimal or maximal value, except do not
Packit 6c4009
   yield a value equal to *T.  */
Packit 6c4009
static time_t
Packit 6c4009
guess_time_tm (long_int year, long_int yday, int hour, int min, int sec,
Packit 6c4009
	       const time_t *t, const struct tm *tp)
Packit 6c4009
{
Packit 6c4009
  if (tp)
Packit 6c4009
    {
Packit 6c4009
      time_t d = ydhms_diff (year, yday, hour, min, sec,
Packit 6c4009
			     tp->tm_year, tp->tm_yday,
Packit 6c4009
			     tp->tm_hour, tp->tm_min, tp->tm_sec);
Packit 6c4009
      if (time_t_add_ok (*t, d))
Packit 6c4009
	return *t + d;
Packit 6c4009
    }
Packit 6c4009
Packit 6c4009
  /* Overflow occurred one way or another.  Return the nearest result
Packit 6c4009
     that is actually in range, except don't report a zero difference
Packit 6c4009
     if the actual difference is nonzero, as that would cause a false
Packit 6c4009
     match; and don't oscillate between two values, as that would
Packit 6c4009
     confuse the spring-forward gap detector.  */
Packit 6c4009
  return (*t < TIME_T_MIDPOINT
Packit 6c4009
	  ? (*t <= TIME_T_MIN + 1 ? *t + 1 : TIME_T_MIN)
Packit 6c4009
	  : (TIME_T_MAX - 1 <= *t ? *t - 1 : TIME_T_MAX));
Packit 6c4009
}
Packit 6c4009
Packit 6c4009
/* Use CONVERT to convert *T to a broken down time in *TP.
Packit 6c4009
   If *T is out of range for conversion, adjust it so that
Packit 6c4009
   it is the nearest in-range value and then convert that.  */
Packit 6c4009
static struct tm *
Packit 6c4009
ranged_convert (struct tm *(*convert) (const time_t *, struct tm *),
Packit 6c4009
		time_t *t, struct tm *tp)
Packit 6c4009
{
Packit 6c4009
  struct tm *r = convert (t, tp);
Packit 6c4009
Packit 6c4009
  if (!r && *t)
Packit 6c4009
    {
Packit 6c4009
      time_t bad = *t;
Packit 6c4009
      time_t ok = 0;
Packit 6c4009
Packit 6c4009
      /* BAD is a known unconvertible time_t, and OK is a known good one.
Packit 6c4009
	 Use binary search to narrow the range between BAD and OK until
Packit 6c4009
	 they differ by 1.  */
Packit 6c4009
      while (bad != ok + (bad < 0 ? -1 : 1))
Packit 6c4009
	{
Packit 6c4009
	  time_t mid = *t = time_t_avg (ok, bad);
Packit 6c4009
	  r = convert (t, tp);
Packit 6c4009
	  if (r)
Packit 6c4009
	    ok = mid;
Packit 6c4009
	  else
Packit 6c4009
	    bad = mid;
Packit 6c4009
	}
Packit 6c4009
Packit 6c4009
      if (!r && ok)
Packit 6c4009
	{
Packit 6c4009
	  /* The last conversion attempt failed;
Packit 6c4009
	     revert to the most recent successful attempt.  */
Packit 6c4009
	  *t = ok;
Packit 6c4009
	  r = convert (t, tp);
Packit 6c4009
	}
Packit 6c4009
    }
Packit 6c4009
Packit 6c4009
  return r;
Packit 6c4009
}
Packit 6c4009
Packit 6c4009
Packit 6c4009
/* Convert *TP to a time_t value, inverting
Packit 6c4009
   the monotonic and mostly-unit-linear conversion function CONVERT.
Packit 6c4009
   Use *OFFSET to keep track of a guess at the offset of the result,
Packit 6c4009
   compared to what the result would be for UTC without leap seconds.
Packit 6c4009
   If *OFFSET's guess is correct, only one CONVERT call is needed.
Packit 6c4009
   This function is external because it is used also by timegm.c.  */
Packit 6c4009
time_t
Packit 6c4009
__mktime_internal (struct tm *tp,
Packit 6c4009
		   struct tm *(*convert) (const time_t *, struct tm *),
Packit 6c4009
		   time_t *offset)
Packit 6c4009
{
Packit 6c4009
  time_t t, gt, t0, t1, t2;
Packit 6c4009
  struct tm tm;
Packit 6c4009
Packit 6c4009
  /* The maximum number of probes (calls to CONVERT) should be enough
Packit 6c4009
     to handle any combinations of time zone rule changes, solar time,
Packit 6c4009
     leap seconds, and oscillations around a spring-forward gap.
Packit 6c4009
     POSIX.1 prohibits leap seconds, but some hosts have them anyway.  */
Packit 6c4009
  int remaining_probes = 6;
Packit 6c4009
Packit 6c4009
  /* Time requested.  Copy it in case CONVERT modifies *TP; this can
Packit 6c4009
     occur if TP is localtime's returned value and CONVERT is localtime.  */
Packit 6c4009
  int sec = tp->tm_sec;
Packit 6c4009
  int min = tp->tm_min;
Packit 6c4009
  int hour = tp->tm_hour;
Packit 6c4009
  int mday = tp->tm_mday;
Packit 6c4009
  int mon = tp->tm_mon;
Packit 6c4009
  int year_requested = tp->tm_year;
Packit 6c4009
  int isdst = tp->tm_isdst;
Packit 6c4009
Packit 6c4009
  /* 1 if the previous probe was DST.  */
Packit 6c4009
  int dst2;
Packit 6c4009
Packit 6c4009
  /* Ensure that mon is in range, and set year accordingly.  */
Packit 6c4009
  int mon_remainder = mon % 12;
Packit 6c4009
  int negative_mon_remainder = mon_remainder < 0;
Packit 6c4009
  int mon_years = mon / 12 - negative_mon_remainder;
Packit 6c4009
  long_int lyear_requested = year_requested;
Packit 6c4009
  long_int year = lyear_requested + mon_years;
Packit 6c4009
Packit 6c4009
  /* The other values need not be in range:
Packit 6c4009
     the remaining code handles minor overflows correctly,
Packit 6c4009
     assuming int and time_t arithmetic wraps around.
Packit 6c4009
     Major overflows are caught at the end.  */
Packit 6c4009
Packit 6c4009
  /* Calculate day of year from year, month, and day of month.
Packit 6c4009
     The result need not be in range.  */
Packit 6c4009
  int mon_yday = ((__mon_yday[leapyear (year)]
Packit 6c4009
		   [mon_remainder + 12 * negative_mon_remainder])
Packit 6c4009
		  - 1);
Packit 6c4009
  long_int lmday = mday;
Packit 6c4009
  long_int yday = mon_yday + lmday;
Packit 6c4009
Packit 6c4009
  time_t guessed_offset = *offset;
Packit 6c4009
Packit 6c4009
  int sec_requested = sec;
Packit 6c4009
Packit 6c4009
  if (LEAP_SECONDS_POSSIBLE)
Packit 6c4009
    {
Packit 6c4009
      /* Handle out-of-range seconds specially,
Packit 6c4009
	 since ydhms_tm_diff assumes every minute has 60 seconds.  */
Packit 6c4009
      if (sec < 0)
Packit 6c4009
	sec = 0;
Packit 6c4009
      if (59 < sec)
Packit 6c4009
	sec = 59;
Packit 6c4009
    }
Packit 6c4009
Packit 6c4009
  /* Invert CONVERT by probing.  First assume the same offset as last
Packit 6c4009
     time.  */
Packit 6c4009
Packit 6c4009
  t0 = ydhms_diff (year, yday, hour, min, sec,
Packit 6c4009
		   EPOCH_YEAR - TM_YEAR_BASE, 0, 0, 0, - guessed_offset);
Packit 6c4009
Packit 6c4009
  if (TIME_T_MAX / INT_MAX / 366 / 24 / 60 / 60 < 3)
Packit 6c4009
    {
Packit 6c4009
      /* time_t isn't large enough to rule out overflows, so check
Packit 6c4009
	 for major overflows.  A gross check suffices, since if t0
Packit 6c4009
	 has overflowed, it is off by a multiple of TIME_T_MAX -
Packit 6c4009
	 TIME_T_MIN + 1.  So ignore any component of the difference
Packit 6c4009
	 that is bounded by a small value.  */
Packit 6c4009
Packit 6c4009
      /* Approximate log base 2 of the number of time units per
Packit 6c4009
	 biennium.  A biennium is 2 years; use this unit instead of
Packit 6c4009
	 years to avoid integer overflow.  For example, 2 average
Packit 6c4009
	 Gregorian years are 2 * 365.2425 * 24 * 60 * 60 seconds,
Packit 6c4009
	 which is 63113904 seconds, and rint (log2 (63113904)) is
Packit 6c4009
	 26.  */
Packit 6c4009
      int ALOG2_SECONDS_PER_BIENNIUM = 26;
Packit 6c4009
      int ALOG2_MINUTES_PER_BIENNIUM = 20;
Packit 6c4009
      int ALOG2_HOURS_PER_BIENNIUM = 14;
Packit 6c4009
      int ALOG2_DAYS_PER_BIENNIUM = 10;
Packit 6c4009
      int LOG2_YEARS_PER_BIENNIUM = 1;
Packit 6c4009
Packit 6c4009
      int approx_requested_biennia =
Packit 6c4009
	(SHR (year_requested, LOG2_YEARS_PER_BIENNIUM)
Packit 6c4009
	 - SHR (EPOCH_YEAR - TM_YEAR_BASE, LOG2_YEARS_PER_BIENNIUM)
Packit 6c4009
	 + SHR (mday, ALOG2_DAYS_PER_BIENNIUM)
Packit 6c4009
	 + SHR (hour, ALOG2_HOURS_PER_BIENNIUM)
Packit 6c4009
	 + SHR (min, ALOG2_MINUTES_PER_BIENNIUM)
Packit 6c4009
	 + (LEAP_SECONDS_POSSIBLE
Packit 6c4009
	    ? 0
Packit 6c4009
	    : SHR (sec, ALOG2_SECONDS_PER_BIENNIUM)));
Packit 6c4009
Packit 6c4009
      int approx_biennia = SHR (t0, ALOG2_SECONDS_PER_BIENNIUM);
Packit 6c4009
      int diff = approx_biennia - approx_requested_biennia;
Packit 6c4009
      int approx_abs_diff = diff < 0 ? -1 - diff : diff;
Packit 6c4009
Packit 6c4009
      /* IRIX 4.0.5 cc miscalculates TIME_T_MIN / 3: it erroneously
Packit 6c4009
	 gives a positive value of 715827882.  Setting a variable
Packit 6c4009
	 first then doing math on it seems to work.
Packit 6c4009
	 (ghazi@caip.rutgers.edu) */
Packit 6c4009
      time_t time_t_max = TIME_T_MAX;
Packit 6c4009
      time_t time_t_min = TIME_T_MIN;
Packit 6c4009
      time_t overflow_threshold =
Packit 6c4009
	(time_t_max / 3 - time_t_min / 3) >> ALOG2_SECONDS_PER_BIENNIUM;
Packit 6c4009
Packit 6c4009
      if (overflow_threshold < approx_abs_diff)
Packit 6c4009
	{
Packit 6c4009
	  /* Overflow occurred.  Try repairing it; this might work if
Packit 6c4009
	     the time zone offset is enough to undo the overflow.  */
Packit 6c4009
	  time_t repaired_t0 = -1 - t0;
Packit 6c4009
	  approx_biennia = SHR (repaired_t0, ALOG2_SECONDS_PER_BIENNIUM);
Packit 6c4009
	  diff = approx_biennia - approx_requested_biennia;
Packit 6c4009
	  approx_abs_diff = diff < 0 ? -1 - diff : diff;
Packit 6c4009
	  if (overflow_threshold < approx_abs_diff)
Packit 6c4009
	    return -1;
Packit 6c4009
	  guessed_offset += repaired_t0 - t0;
Packit 6c4009
	  t0 = repaired_t0;
Packit 6c4009
	}
Packit 6c4009
    }
Packit 6c4009
Packit 6c4009
  /* Repeatedly use the error to improve the guess.  */
Packit 6c4009
Packit 6c4009
  for (t = t1 = t2 = t0, dst2 = 0;
Packit 6c4009
       (gt = guess_time_tm (year, yday, hour, min, sec, &t,
Packit 6c4009
			    ranged_convert (convert, &t, &tm)),
Packit 6c4009
	t != gt);
Packit 6c4009
       t1 = t2, t2 = t, t = gt, dst2 = tm.tm_isdst != 0)
Packit 6c4009
    if (t == t1 && t != t2
Packit 6c4009
	&& (tm.tm_isdst < 0
Packit 6c4009
	    || (isdst < 0
Packit 6c4009
		? dst2 <= (tm.tm_isdst != 0)
Packit 6c4009
		: (isdst != 0) != (tm.tm_isdst != 0))))
Packit 6c4009
      /* We can't possibly find a match, as we are oscillating
Packit 6c4009
	 between two values.  The requested time probably falls
Packit 6c4009
	 within a spring-forward gap of size GT - T.  Follow the common
Packit 6c4009
	 practice in this case, which is to return a time that is GT - T
Packit 6c4009
	 away from the requested time, preferring a time whose
Packit 6c4009
	 tm_isdst differs from the requested value.  (If no tm_isdst
Packit 6c4009
	 was requested and only one of the two values has a nonzero
Packit 6c4009
	 tm_isdst, prefer that value.)  In practice, this is more
Packit 6c4009
	 useful than returning -1.  */
Packit 6c4009
      goto offset_found;
Packit 6c4009
    else if (--remaining_probes == 0)
Packit 6c4009
      return -1;
Packit 6c4009
Packit 6c4009
  /* We have a match.  Check whether tm.tm_isdst has the requested
Packit 6c4009
     value, if any.  */
Packit 6c4009
  if (isdst_differ (isdst, tm.tm_isdst))
Packit 6c4009
    {
Packit 6c4009
      /* tm.tm_isdst has the wrong value.  Look for a neighboring
Packit 6c4009
	 time with the right value, and use its UTC offset.
Packit 6c4009
Packit 6c4009
	 Heuristic: probe the adjacent timestamps in both directions,
Packit 6c4009
	 looking for the desired isdst.  This should work for all real
Packit 6c4009
	 time zone histories in the tz database.  */
Packit 6c4009
Packit 6c4009
      /* Distance between probes when looking for a DST boundary.  In
Packit 6c4009
	 tzdata2003a, the shortest period of DST is 601200 seconds
Packit 6c4009
	 (e.g., America/Recife starting 2000-10-08 01:00), and the
Packit 6c4009
	 shortest period of non-DST surrounded by DST is 694800
Packit 6c4009
	 seconds (Africa/Tunis starting 1943-04-17 01:00).  Use the
Packit 6c4009
	 minimum of these two values, so we don't miss these short
Packit 6c4009
	 periods when probing.  */
Packit 6c4009
      int stride = 601200;
Packit 6c4009
Packit 6c4009
      /* The longest period of DST in tzdata2003a is 536454000 seconds
Packit 6c4009
	 (e.g., America/Jujuy starting 1946-10-01 01:00).  The longest
Packit 6c4009
	 period of non-DST is much longer, but it makes no real sense
Packit 6c4009
	 to search for more than a year of non-DST, so use the DST
Packit 6c4009
	 max.  */
Packit 6c4009
      int duration_max = 536454000;
Packit 6c4009
Packit 6c4009
      /* Search in both directions, so the maximum distance is half
Packit 6c4009
	 the duration; add the stride to avoid off-by-1 problems.  */
Packit 6c4009
      int delta_bound = duration_max / 2 + stride;
Packit 6c4009
Packit 6c4009
      int delta, direction;
Packit 6c4009
Packit 6c4009
      for (delta = stride; delta < delta_bound; delta += stride)
Packit 6c4009
	for (direction = -1; direction <= 1; direction += 2)
Packit 6c4009
	  if (time_t_int_add_ok (t, delta * direction))
Packit 6c4009
	    {
Packit 6c4009
	      time_t ot = t + delta * direction;
Packit 6c4009
	      struct tm otm;
Packit 6c4009
	      ranged_convert (convert, &ot, &otm;;
Packit 6c4009
	      if (! isdst_differ (isdst, otm.tm_isdst))
Packit 6c4009
		{
Packit 6c4009
		  /* We found the desired tm_isdst.
Packit 6c4009
		     Extrapolate back to the desired time.  */
Packit 6c4009
		  t = guess_time_tm (year, yday, hour, min, sec, &ot, &otm;;
Packit 6c4009
		  ranged_convert (convert, &t, &tm;;
Packit 6c4009
		  goto offset_found;
Packit 6c4009
		}
Packit 6c4009
	    }
Packit 6c4009
    }
Packit 6c4009
Packit 6c4009
 offset_found:
Packit 6c4009
  *offset = guessed_offset + t - t0;
Packit 6c4009
Packit 6c4009
  if (LEAP_SECONDS_POSSIBLE && sec_requested != tm.tm_sec)
Packit 6c4009
    {
Packit 6c4009
      /* Adjust time to reflect the tm_sec requested, not the normalized value.
Packit 6c4009
	 Also, repair any damage from a false match due to a leap second.  */
Packit 6c4009
      int sec_adjustment = (sec == 0 && tm.tm_sec == 60) - sec;
Packit 6c4009
      if (! time_t_int_add_ok (t, sec_requested))
Packit 6c4009
	return -1;
Packit 6c4009
      t1 = t + sec_requested;
Packit 6c4009
      if (! time_t_int_add_ok (t1, sec_adjustment))
Packit 6c4009
	return -1;
Packit 6c4009
      t2 = t1 + sec_adjustment;
Packit 6c4009
      if (! convert (&t2, &tm))
Packit 6c4009
	return -1;
Packit 6c4009
      t = t2;
Packit 6c4009
    }
Packit 6c4009
Packit 6c4009
  *tp = tm;
Packit 6c4009
  return t;
Packit 6c4009
}
Packit 6c4009
Packit 6c4009
Packit 6c4009
/* FIXME: This should use a signed type wide enough to hold any UTC
Packit 6c4009
   offset in seconds.  'int' should be good enough for GNU code.  We
Packit 6c4009
   can't fix this unilaterally though, as other modules invoke
Packit 6c4009
   __mktime_internal.  */
Packit 6c4009
static time_t localtime_offset;
Packit 6c4009
Packit 6c4009
/* Convert *TP to a time_t value.  */
Packit 6c4009
time_t
Packit 6c4009
mktime (struct tm *tp)
Packit 6c4009
{
Packit 6c4009
#ifdef _LIBC
Packit 6c4009
  /* POSIX.1 8.1.1 requires that whenever mktime() is called, the
Packit 6c4009
     time zone names contained in the external variable 'tzname' shall
Packit 6c4009
     be set as if the tzset() function had been called.  */
Packit 6c4009
  __tzset ();
Packit 6c4009
#endif
Packit 6c4009
Packit 6c4009
  return __mktime_internal (tp, __localtime_r, &localtime_offset);
Packit 6c4009
}
Packit 6c4009
Packit 6c4009
#ifdef weak_alias
Packit 6c4009
weak_alias (mktime, timelocal)
Packit 6c4009
#endif
Packit 6c4009
Packit 6c4009
#ifdef _LIBC
Packit 6c4009
libc_hidden_def (mktime)
Packit 6c4009
libc_hidden_weak (timelocal)
Packit 6c4009
#endif
Packit 6c4009

Packit 6c4009
#if defined DEBUG_MKTIME && DEBUG_MKTIME
Packit 6c4009
Packit 6c4009
static int
Packit 6c4009
not_equal_tm (const struct tm *a, const struct tm *b)
Packit 6c4009
{
Packit 6c4009
  return ((a->tm_sec ^ b->tm_sec)
Packit 6c4009
	  | (a->tm_min ^ b->tm_min)
Packit 6c4009
	  | (a->tm_hour ^ b->tm_hour)
Packit 6c4009
	  | (a->tm_mday ^ b->tm_mday)
Packit 6c4009
	  | (a->tm_mon ^ b->tm_mon)
Packit 6c4009
	  | (a->tm_year ^ b->tm_year)
Packit 6c4009
	  | (a->tm_yday ^ b->tm_yday)
Packit 6c4009
	  | isdst_differ (a->tm_isdst, b->tm_isdst));
Packit 6c4009
}
Packit 6c4009
Packit 6c4009
static void
Packit 6c4009
print_tm (const struct tm *tp)
Packit 6c4009
{
Packit 6c4009
  if (tp)
Packit 6c4009
    printf ("%04d-%02d-%02d %02d:%02d:%02d yday %03d wday %d isdst %d",
Packit 6c4009
	    tp->tm_year + TM_YEAR_BASE, tp->tm_mon + 1, tp->tm_mday,
Packit 6c4009
	    tp->tm_hour, tp->tm_min, tp->tm_sec,
Packit 6c4009
	    tp->tm_yday, tp->tm_wday, tp->tm_isdst);
Packit 6c4009
  else
Packit 6c4009
    printf ("0");
Packit 6c4009
}
Packit 6c4009
Packit 6c4009
static int
Packit 6c4009
check_result (time_t tk, struct tm tmk, time_t tl, const struct tm *lt)
Packit 6c4009
{
Packit 6c4009
  if (tk != tl || !lt || not_equal_tm (&tmk, lt))
Packit 6c4009
    {
Packit 6c4009
      printf ("mktime (");
Packit 6c4009
      print_tm (lt);
Packit 6c4009
      printf (")\nyields (");
Packit 6c4009
      print_tm (&tmk);
Packit 6c4009
      printf (") == %ld, should be %ld\n", (long int) tk, (long int) tl);
Packit 6c4009
      return 1;
Packit 6c4009
    }
Packit 6c4009
Packit 6c4009
  return 0;
Packit 6c4009
}
Packit 6c4009
Packit 6c4009
int
Packit 6c4009
main (int argc, char **argv)
Packit 6c4009
{
Packit 6c4009
  int status = 0;
Packit 6c4009
  struct tm tm, tmk, tml;
Packit 6c4009
  struct tm *lt;
Packit 6c4009
  time_t tk, tl, tl1;
Packit 6c4009
  char trailer;
Packit 6c4009
Packit 6c4009
  if ((argc == 3 || argc == 4)
Packit 6c4009
      && (sscanf (argv[1], "%d-%d-%d%c",
Packit 6c4009
		  &tm.tm_year, &tm.tm_mon, &tm.tm_mday, &trailer)
Packit 6c4009
	  == 3)
Packit 6c4009
      && (sscanf (argv[2], "%d:%d:%d%c",
Packit 6c4009
		  &tm.tm_hour, &tm.tm_min, &tm.tm_sec, &trailer)
Packit 6c4009
	  == 3))
Packit 6c4009
    {
Packit 6c4009
      tm.tm_year -= TM_YEAR_BASE;
Packit 6c4009
      tm.tm_mon--;
Packit 6c4009
      tm.tm_isdst = argc == 3 ? -1 : atoi (argv[3]);
Packit 6c4009
      tmk = tm;
Packit 6c4009
      tl = mktime (&tmk);
Packit 6c4009
      lt = localtime (&tl;;
Packit 6c4009
      if (lt)
Packit 6c4009
	{
Packit 6c4009
	  tml = *lt;
Packit 6c4009
	  lt = &tm;;
Packit 6c4009
	}
Packit 6c4009
      printf ("mktime returns %ld == ", (long int) tl);
Packit 6c4009
      print_tm (&tmk);
Packit 6c4009
      printf ("\n");
Packit 6c4009
      status = check_result (tl, tmk, tl, lt);
Packit 6c4009
    }
Packit 6c4009
  else if (argc == 4 || (argc == 5 && strcmp (argv[4], "-") == 0))
Packit 6c4009
    {
Packit 6c4009
      time_t from = atol (argv[1]);
Packit 6c4009
      time_t by = atol (argv[2]);
Packit 6c4009
      time_t to = atol (argv[3]);
Packit 6c4009
Packit 6c4009
      if (argc == 4)
Packit 6c4009
	for (tl = from; by < 0 ? to <= tl : tl <= to; tl = tl1)
Packit 6c4009
	  {
Packit 6c4009
	    lt = localtime (&tl;;
Packit 6c4009
	    if (lt)
Packit 6c4009
	      {
Packit 6c4009
		tmk = tml = *lt;
Packit 6c4009
		tk = mktime (&tmk);
Packit 6c4009
		status |= check_result (tk, tmk, tl, &tml);
Packit 6c4009
	      }
Packit 6c4009
	    else
Packit 6c4009
	      {
Packit 6c4009
		printf ("localtime (%ld) yields 0\n", (long int) tl);
Packit 6c4009
		status = 1;
Packit 6c4009
	      }
Packit 6c4009
	    tl1 = tl + by;
Packit 6c4009
	    if ((tl1 < tl) != (by < 0))
Packit 6c4009
	      break;
Packit 6c4009
	  }
Packit 6c4009
      else
Packit 6c4009
	for (tl = from; by < 0 ? to <= tl : tl <= to; tl = tl1)
Packit 6c4009
	  {
Packit 6c4009
	    /* Null benchmark.  */
Packit 6c4009
	    lt = localtime (&tl;;
Packit 6c4009
	    if (lt)
Packit 6c4009
	      {
Packit 6c4009
		tmk = tml = *lt;
Packit 6c4009
		tk = tl;
Packit 6c4009
		status |= check_result (tk, tmk, tl, &tml);
Packit 6c4009
	      }
Packit 6c4009
	    else
Packit 6c4009
	      {
Packit 6c4009
		printf ("localtime (%ld) yields 0\n", (long int) tl);
Packit 6c4009
		status = 1;
Packit 6c4009
	      }
Packit 6c4009
	    tl1 = tl + by;
Packit 6c4009
	    if ((tl1 < tl) != (by < 0))
Packit 6c4009
	      break;
Packit 6c4009
	  }
Packit 6c4009
    }
Packit 6c4009
  else
Packit 6c4009
    printf ("Usage:\
Packit 6c4009
\t%s YYYY-MM-DD HH:MM:SS [ISDST] # Test given time.\n\
Packit 6c4009
\t%s FROM BY TO # Test values FROM, FROM+BY, ..., TO.\n\
Packit 6c4009
\t%s FROM BY TO - # Do not test those values (for benchmark).\n",
Packit 6c4009
	    argv[0], argv[0], argv[0]);
Packit 6c4009
Packit 6c4009
  return status;
Packit 6c4009
}
Packit 6c4009
Packit 6c4009
#endif /* DEBUG_MKTIME */
Packit 6c4009

Packit 6c4009
/*
Packit 6c4009
Local Variables:
Packit 6c4009
compile-command: "gcc -DDEBUG_MKTIME -I. -Wall -W -O2 -g mktime.c -o mktime"
Packit 6c4009
End:
Packit 6c4009
*/