Coverage Report

Created: 2023-02-22 06:39

/src/gmp-6.2.1/mpz/aors.h
Line
Count
Source (jump to first uncovered line)
1
/* mpz_add, mpz_sub -- add or subtract integers.
2
3
Copyright 1991, 1993, 1994, 1996, 2000, 2001, 2011, 2012 Free Software
4
Foundation, Inc.
5
6
This file is part of the GNU MP Library.
7
8
The GNU MP Library is free software; you can redistribute it and/or modify
9
it under the terms of either:
10
11
  * the GNU Lesser General Public License as published by the Free
12
    Software Foundation; either version 3 of the License, or (at your
13
    option) any later version.
14
15
or
16
17
  * the GNU General Public License as published by the Free Software
18
    Foundation; either version 2 of the License, or (at your option) any
19
    later version.
20
21
or both in parallel, as here.
22
23
The GNU MP Library is distributed in the hope that it will be useful, but
24
WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY
25
or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License
26
for more details.
27
28
You should have received copies of the GNU General Public License and the
29
GNU Lesser General Public License along with the GNU MP Library.  If not,
30
see https://www.gnu.org/licenses/.  */
31
32
#include "gmp-impl.h"
33
34
35
#ifdef OPERATION_add
36
#define FUNCTION     mpz_add
37
#define VARIATION
38
#endif
39
#ifdef OPERATION_sub
40
#define FUNCTION     mpz_sub
41
1.62k
#define VARIATION    -
42
#endif
43
44
#ifndef FUNCTION
45
Error, need OPERATION_add or OPERATION_sub
46
#endif
47
48
49
void
50
FUNCTION (mpz_ptr w, mpz_srcptr u, mpz_srcptr v)
51
2.32k
{
52
2.32k
  mp_srcptr up, vp;
53
2.32k
  mp_ptr wp;
54
2.32k
  mp_size_t usize, vsize, wsize;
55
2.32k
  mp_size_t abs_usize;
56
2.32k
  mp_size_t abs_vsize;
57
58
2.32k
  usize = SIZ(u);
59
2.32k
  vsize = VARIATION SIZ(v);
60
2.32k
  abs_usize = ABS (usize);
61
2.32k
  abs_vsize = ABS (vsize);
62
63
2.32k
  if (abs_usize < abs_vsize)
64
1.64k
    {
65
      /* Swap U and V. */
66
1.64k
      MPZ_SRCPTR_SWAP (u, v);
67
1.64k
      MP_SIZE_T_SWAP (usize, vsize);
68
1.64k
      MP_SIZE_T_SWAP (abs_usize, abs_vsize);
69
1.64k
    }
70
71
  /* True: ABS_USIZE >= ABS_VSIZE.  */
72
73
  /* If not space for w (and possible carry), increase space.  */
74
2.32k
  wsize = abs_usize + 1;
75
2.32k
  wp = MPZ_REALLOC (w, wsize);
76
77
  /* These must be after realloc (u or v may be the same as w).  */
78
2.32k
  up = PTR(u);
79
2.32k
  vp = PTR(v);
80
81
2.32k
  if ((usize ^ vsize) < 0)
82
1.43k
    {
83
      /* U and V have different sign.  Need to compare them to determine
84
   which operand to subtract from which.  */
85
86
      /* This test is right since ABS_USIZE >= ABS_VSIZE.  */
87
1.43k
      if (abs_usize != abs_vsize)
88
969
  {
89
969
    mpn_sub (wp, up, abs_usize, vp, abs_vsize);
90
969
    wsize = abs_usize;
91
969
    MPN_NORMALIZE (wp, wsize);
92
969
    if (usize < 0)
93
821
      wsize = -wsize;
94
969
  }
95
465
      else if (mpn_cmp (up, vp, abs_usize) < 0)
96
438
  {
97
438
    mpn_sub_n (wp, vp, up, abs_usize);
98
438
    wsize = abs_usize;
99
438
    MPN_NORMALIZE (wp, wsize);
100
438
    if (usize >= 0)
101
50
      wsize = -wsize;
102
438
  }
103
27
      else
104
27
  {
105
27
    mpn_sub_n (wp, up, vp, abs_usize);
106
27
    wsize = abs_usize;
107
27
    MPN_NORMALIZE (wp, wsize);
108
27
    if (usize < 0)
109
0
      wsize = -wsize;
110
27
  }
111
1.43k
    }
112
891
  else
113
891
    {
114
      /* U and V have same sign.  Add them.  */
115
891
      mp_limb_t cy_limb = mpn_add (wp, up, abs_usize, vp, abs_vsize);
116
891
      wp[abs_usize] = cy_limb;
117
891
      wsize = abs_usize + cy_limb;
118
891
      if (usize < 0)
119
0
  wsize = -wsize;
120
891
    }
121
122
2.32k
  SIZ(w) = wsize;
123
2.32k
}
__gmpz_add
Line
Count
Source
51
697
{
52
697
  mp_srcptr up, vp;
53
697
  mp_ptr wp;
54
697
  mp_size_t usize, vsize, wsize;
55
697
  mp_size_t abs_usize;
56
697
  mp_size_t abs_vsize;
57
58
697
  usize = SIZ(u);
59
697
  vsize = VARIATION SIZ(v);
60
697
  abs_usize = ABS (usize);
61
697
  abs_vsize = ABS (vsize);
62
63
697
  if (abs_usize < abs_vsize)
64
193
    {
65
      /* Swap U and V. */
66
193
      MPZ_SRCPTR_SWAP (u, v);
67
193
      MP_SIZE_T_SWAP (usize, vsize);
68
193
      MP_SIZE_T_SWAP (abs_usize, abs_vsize);
69
193
    }
70
71
  /* True: ABS_USIZE >= ABS_VSIZE.  */
72
73
  /* If not space for w (and possible carry), increase space.  */
74
697
  wsize = abs_usize + 1;
75
697
  wp = MPZ_REALLOC (w, wsize);
76
77
  /* These must be after realloc (u or v may be the same as w).  */
78
697
  up = PTR(u);
79
697
  vp = PTR(v);
80
81
697
  if ((usize ^ vsize) < 0)
82
488
    {
83
      /* U and V have different sign.  Need to compare them to determine
84
   which operand to subtract from which.  */
85
86
      /* This test is right since ABS_USIZE >= ABS_VSIZE.  */
87
488
      if (abs_usize != abs_vsize)
88
100
  {
89
100
    mpn_sub (wp, up, abs_usize, vp, abs_vsize);
90
100
    wsize = abs_usize;
91
100
    MPN_NORMALIZE (wp, wsize);
92
100
    if (usize < 0)
93
0
      wsize = -wsize;
94
100
  }
95
388
      else if (mpn_cmp (up, vp, abs_usize) < 0)
96
388
  {
97
388
    mpn_sub_n (wp, vp, up, abs_usize);
98
388
    wsize = abs_usize;
99
388
    MPN_NORMALIZE (wp, wsize);
100
388
    if (usize >= 0)
101
0
      wsize = -wsize;
102
388
  }
103
0
      else
104
0
  {
105
0
    mpn_sub_n (wp, up, vp, abs_usize);
106
0
    wsize = abs_usize;
107
0
    MPN_NORMALIZE (wp, wsize);
108
0
    if (usize < 0)
109
0
      wsize = -wsize;
110
0
  }
111
488
    }
112
209
  else
113
209
    {
114
      /* U and V have same sign.  Add them.  */
115
209
      mp_limb_t cy_limb = mpn_add (wp, up, abs_usize, vp, abs_vsize);
116
209
      wp[abs_usize] = cy_limb;
117
209
      wsize = abs_usize + cy_limb;
118
209
      if (usize < 0)
119
0
  wsize = -wsize;
120
209
    }
121
122
697
  SIZ(w) = wsize;
123
697
}
__gmpz_sub
Line
Count
Source
51
1.62k
{
52
1.62k
  mp_srcptr up, vp;
53
1.62k
  mp_ptr wp;
54
1.62k
  mp_size_t usize, vsize, wsize;
55
1.62k
  mp_size_t abs_usize;
56
1.62k
  mp_size_t abs_vsize;
57
58
1.62k
  usize = SIZ(u);
59
1.62k
  vsize = VARIATION SIZ(v);
60
1.62k
  abs_usize = ABS (usize);
61
1.62k
  abs_vsize = ABS (vsize);
62
63
1.62k
  if (abs_usize < abs_vsize)
64
1.45k
    {
65
      /* Swap U and V. */
66
1.45k
      MPZ_SRCPTR_SWAP (u, v);
67
1.45k
      MP_SIZE_T_SWAP (usize, vsize);
68
1.45k
      MP_SIZE_T_SWAP (abs_usize, abs_vsize);
69
1.45k
    }
70
71
  /* True: ABS_USIZE >= ABS_VSIZE.  */
72
73
  /* If not space for w (and possible carry), increase space.  */
74
1.62k
  wsize = abs_usize + 1;
75
1.62k
  wp = MPZ_REALLOC (w, wsize);
76
77
  /* These must be after realloc (u or v may be the same as w).  */
78
1.62k
  up = PTR(u);
79
1.62k
  vp = PTR(v);
80
81
1.62k
  if ((usize ^ vsize) < 0)
82
946
    {
83
      /* U and V have different sign.  Need to compare them to determine
84
   which operand to subtract from which.  */
85
86
      /* This test is right since ABS_USIZE >= ABS_VSIZE.  */
87
946
      if (abs_usize != abs_vsize)
88
869
  {
89
869
    mpn_sub (wp, up, abs_usize, vp, abs_vsize);
90
869
    wsize = abs_usize;
91
869
    MPN_NORMALIZE (wp, wsize);
92
869
    if (usize < 0)
93
821
      wsize = -wsize;
94
869
  }
95
77
      else if (mpn_cmp (up, vp, abs_usize) < 0)
96
50
  {
97
50
    mpn_sub_n (wp, vp, up, abs_usize);
98
50
    wsize = abs_usize;
99
50
    MPN_NORMALIZE (wp, wsize);
100
50
    if (usize >= 0)
101
50
      wsize = -wsize;
102
50
  }
103
27
      else
104
27
  {
105
27
    mpn_sub_n (wp, up, vp, abs_usize);
106
27
    wsize = abs_usize;
107
27
    MPN_NORMALIZE (wp, wsize);
108
27
    if (usize < 0)
109
0
      wsize = -wsize;
110
27
  }
111
946
    }
112
682
  else
113
682
    {
114
      /* U and V have same sign.  Add them.  */
115
682
      mp_limb_t cy_limb = mpn_add (wp, up, abs_usize, vp, abs_vsize);
116
682
      wp[abs_usize] = cy_limb;
117
682
      wsize = abs_usize + cy_limb;
118
682
      if (usize < 0)
119
0
  wsize = -wsize;
120
682
    }
121
122
1.62k
  SIZ(w) = wsize;
123
1.62k
}