Coverage Report

Created: 2023-02-22 06:39

/src/gmp-6.2.1/mpz/aors_ui.h
Line
Count
Source (jump to first uncovered line)
1
/* mpz_add_ui, mpz_sub_ui -- Add or subtract an mpz_t and an unsigned
2
   one-word integer.
3
4
Copyright 1991, 1993, 1994, 1996, 1999-2002, 2004, 2012, 2013, 2015
5
Free Software Foundation, Inc.
6
7
This file is part of the GNU MP Library.
8
9
The GNU MP Library is free software; you can redistribute it and/or modify
10
it under the terms of either:
11
12
  * the GNU Lesser General Public License as published by the Free
13
    Software Foundation; either version 3 of the License, or (at your
14
    option) any later version.
15
16
or
17
18
  * the GNU General Public License as published by the Free Software
19
    Foundation; either version 2 of the License, or (at your option) any
20
    later version.
21
22
or both in parallel, as here.
23
24
The GNU MP Library is distributed in the hope that it will be useful, but
25
WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY
26
or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License
27
for more details.
28
29
You should have received copies of the GNU General Public License and the
30
GNU Lesser General Public License along with the GNU MP Library.  If not,
31
see https://www.gnu.org/licenses/.  */
32
33
#include "gmp-impl.h"
34
35
36
#ifdef OPERATION_add_ui
37
#define FUNCTION          mpz_add_ui
38
#define FUNCTION2         mpz_add
39
#define VARIATION_CMP     >=
40
#define VARIATION_NEG
41
0
#define VARIATION_UNNEG   -
42
#endif
43
44
#ifdef OPERATION_sub_ui
45
#define FUNCTION          mpz_sub_ui
46
#define FUNCTION2         mpz_sub
47
#define VARIATION_CMP     <
48
6
#define VARIATION_NEG     -
49
#define VARIATION_UNNEG
50
#endif
51
52
#ifndef FUNCTION
53
Error, need OPERATION_add_ui or OPERATION_sub_ui
54
#endif
55
56
57
void
58
FUNCTION (mpz_ptr w, mpz_srcptr u, unsigned long int vval)
59
122
{
60
122
  mp_srcptr up;
61
122
  mp_ptr wp;
62
122
  mp_size_t usize, wsize;
63
122
  mp_size_t abs_usize;
64
65
#if BITS_PER_ULONG > GMP_NUMB_BITS  /* avoid warnings about shift amount */
66
  if (vval > GMP_NUMB_MAX)
67
    {
68
      mpz_t v;
69
      mp_limb_t vl[2];
70
      PTR(v) = vl;
71
      vl[0] = vval & GMP_NUMB_MASK;
72
      vl[1] = vval >> GMP_NUMB_BITS;
73
      SIZ(v) = 2;
74
      FUNCTION2 (w, u, v);
75
      return;
76
    }
77
#endif
78
79
122
  usize = SIZ (u);
80
122
  if (usize == 0)
81
11
    {
82
11
      MPZ_NEWALLOC (w, 1)[0] = vval;
83
11
      SIZ (w) = VARIATION_NEG (vval != 0);
84
11
      return;
85
11
    }
86
87
111
  abs_usize = ABS (usize);
88
89
  /* If not space for W (and possible carry), increase space.  */
90
111
  wp = MPZ_REALLOC (w, abs_usize + 1);
91
92
  /* These must be after realloc (U may be the same as W).  */
93
111
  up = PTR (u);
94
95
111
  if (usize VARIATION_CMP 0)
96
45
    {
97
45
      mp_limb_t cy;
98
45
      cy = mpn_add_1 (wp, up, abs_usize, (mp_limb_t) vval);
99
45
      wp[abs_usize] = cy;
100
45
      wsize = VARIATION_NEG (abs_usize + cy);
101
45
    }
102
66
  else
103
66
    {
104
      /* The signs are different.  Need exact comparison to determine
105
   which operand to subtract from which.  */
106
66
      if (abs_usize == 1 && up[0] < vval)
107
3
  {
108
3
    wp[0] = vval - up[0];
109
3
    wsize = VARIATION_NEG 1;
110
3
  }
111
63
      else
112
63
  {
113
63
    mpn_sub_1 (wp, up, abs_usize, (mp_limb_t) vval);
114
    /* Size can decrease with at most one limb.  */
115
63
    wsize = VARIATION_UNNEG (abs_usize - (wp[abs_usize - 1] == 0));
116
63
  }
117
66
    }
118
119
111
  SIZ (w) = wsize;
120
111
}
__gmpz_add_ui
Line
Count
Source
59
53
{
60
53
  mp_srcptr up;
61
53
  mp_ptr wp;
62
53
  mp_size_t usize, wsize;
63
53
  mp_size_t abs_usize;
64
65
#if BITS_PER_ULONG > GMP_NUMB_BITS  /* avoid warnings about shift amount */
66
  if (vval > GMP_NUMB_MAX)
67
    {
68
      mpz_t v;
69
      mp_limb_t vl[2];
70
      PTR(v) = vl;
71
      vl[0] = vval & GMP_NUMB_MASK;
72
      vl[1] = vval >> GMP_NUMB_BITS;
73
      SIZ(v) = 2;
74
      FUNCTION2 (w, u, v);
75
      return;
76
    }
77
#endif
78
79
53
  usize = SIZ (u);
80
53
  if (usize == 0)
81
8
    {
82
8
      MPZ_NEWALLOC (w, 1)[0] = vval;
83
8
      SIZ (w) = VARIATION_NEG (vval != 0);
84
8
      return;
85
8
    }
86
87
45
  abs_usize = ABS (usize);
88
89
  /* If not space for W (and possible carry), increase space.  */
90
45
  wp = MPZ_REALLOC (w, abs_usize + 1);
91
92
  /* These must be after realloc (U may be the same as W).  */
93
45
  up = PTR (u);
94
95
45
  if (usize VARIATION_CMP 0)
96
45
    {
97
45
      mp_limb_t cy;
98
45
      cy = mpn_add_1 (wp, up, abs_usize, (mp_limb_t) vval);
99
45
      wp[abs_usize] = cy;
100
45
      wsize = VARIATION_NEG (abs_usize + cy);
101
45
    }
102
0
  else
103
0
    {
104
      /* The signs are different.  Need exact comparison to determine
105
   which operand to subtract from which.  */
106
0
      if (abs_usize == 1 && up[0] < vval)
107
0
  {
108
0
    wp[0] = vval - up[0];
109
0
    wsize = VARIATION_NEG 1;
110
0
  }
111
0
      else
112
0
  {
113
0
    mpn_sub_1 (wp, up, abs_usize, (mp_limb_t) vval);
114
    /* Size can decrease with at most one limb.  */
115
0
    wsize = VARIATION_UNNEG (abs_usize - (wp[abs_usize - 1] == 0));
116
0
  }
117
0
    }
118
119
45
  SIZ (w) = wsize;
120
45
}
__gmpz_sub_ui
Line
Count
Source
59
69
{
60
69
  mp_srcptr up;
61
69
  mp_ptr wp;
62
69
  mp_size_t usize, wsize;
63
69
  mp_size_t abs_usize;
64
65
#if BITS_PER_ULONG > GMP_NUMB_BITS  /* avoid warnings about shift amount */
66
  if (vval > GMP_NUMB_MAX)
67
    {
68
      mpz_t v;
69
      mp_limb_t vl[2];
70
      PTR(v) = vl;
71
      vl[0] = vval & GMP_NUMB_MASK;
72
      vl[1] = vval >> GMP_NUMB_BITS;
73
      SIZ(v) = 2;
74
      FUNCTION2 (w, u, v);
75
      return;
76
    }
77
#endif
78
79
69
  usize = SIZ (u);
80
69
  if (usize == 0)
81
3
    {
82
3
      MPZ_NEWALLOC (w, 1)[0] = vval;
83
3
      SIZ (w) = VARIATION_NEG (vval != 0);
84
3
      return;
85
3
    }
86
87
66
  abs_usize = ABS (usize);
88
89
  /* If not space for W (and possible carry), increase space.  */
90
66
  wp = MPZ_REALLOC (w, abs_usize + 1);
91
92
  /* These must be after realloc (U may be the same as W).  */
93
66
  up = PTR (u);
94
95
66
  if (usize VARIATION_CMP 0)
96
0
    {
97
0
      mp_limb_t cy;
98
0
      cy = mpn_add_1 (wp, up, abs_usize, (mp_limb_t) vval);
99
0
      wp[abs_usize] = cy;
100
0
      wsize = VARIATION_NEG (abs_usize + cy);
101
0
    }
102
66
  else
103
66
    {
104
      /* The signs are different.  Need exact comparison to determine
105
   which operand to subtract from which.  */
106
66
      if (abs_usize == 1 && up[0] < vval)
107
3
  {
108
3
    wp[0] = vval - up[0];
109
3
    wsize = VARIATION_NEG 1;
110
3
  }
111
63
      else
112
63
  {
113
63
    mpn_sub_1 (wp, up, abs_usize, (mp_limb_t) vval);
114
    /* Size can decrease with at most one limb.  */
115
63
    wsize = VARIATION_UNNEG (abs_usize - (wp[abs_usize - 1] == 0));
116
63
  }
117
66
    }
118
119
66
  SIZ (w) = wsize;
120
66
}