summaryrefslogtreecommitdiff
path: root/lib/imaxdiv.c
blob: 2502cea845492eb5d2dbfb5d96026004224b725c (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
/* imaxdiv() function: division of 'intmax_t'.
   Copyright (C) 2006, 2009 Free Software Foundation, Inc.

   This program 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.

   This program 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 this program.  If not, see <http://www.gnu.org/licenses/>.  */

#include <config.h>

/* Specification.  */
#include <inttypes.h>

#include <stdlib.h>

imaxdiv_t
imaxdiv (intmax_t numer, intmax_t denom)
{
  imaxdiv_t result;

  result.quot = numer / denom;
  result.rem = numer % denom;

  /* Verify the requirements of ISO C 99 section 6.5.5 paragraph 6:
     "When integers are divided, the result of the / operator is the
      algebraic quotient with any fractional part discarded.  (This is
      often called "truncation toward zero".)  If the quotient a/b is
      representable, the expression (a/b)*b + a%b shall equal a."  */
  if (!(denom == 0
        || (INTMAX_MIN + INTMAX_MAX < 0
            && denom == -1
            && numer < - INTMAX_MAX)))
    {
      if (!(result.quot * denom + result.rem == numer))
        /* The compiler's implementation of / and % is broken.  */
        abort ();
      if (!(numer >= 0
            ? result.rem >= 0
              && (denom >= 0
                  ? result.rem < denom
                  : /* Don't write  result.rem < - denom,
                       as it gives integer overflow if denom == INTMAX_MIN.  */
                    - result.rem > denom)
            : result.rem <= 0
              && (denom >= 0
                  ? result.rem > - denom
                  : result.rem > denom)))
        /* The compiler's implementation of / and % may be ok according to
           C89, but not to C99.  Please report this to <bug-gnulib@ngu.org>.
           This might be a big portability problem.  */
        abort ();
    }

  return result;
}