Coverage Report

Created: 2024-11-21 06:47

/src/libgmp/mpz/mul.c
Line
Count
Source (jump to first uncovered line)
1
/* mpz_mul -- Multiply two integers.
2
3
Copyright 1991, 1993, 1994, 1996, 2000, 2001, 2005, 2009, 2011, 2012,
4
2015 Free Software 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 <stdio.h> /* for NULL */
33
#include "gmp-impl.h"
34
35
36
void
37
mpz_mul (mpz_ptr w, mpz_srcptr u, mpz_srcptr v)
38
11.2k
{
39
11.2k
  mp_size_t usize;
40
11.2k
  mp_size_t vsize;
41
11.2k
  mp_size_t wsize;
42
11.2k
  mp_size_t sign_product;
43
11.2k
  mp_ptr up, vp;
44
11.2k
  mp_ptr wp;
45
11.2k
  mp_ptr free_me;
46
11.2k
  size_t free_me_size;
47
11.2k
  mp_limb_t cy_limb;
48
11.2k
  TMP_DECL;
49
50
11.2k
  usize = SIZ (u);
51
11.2k
  vsize = SIZ (v);
52
11.2k
  sign_product = usize ^ vsize;
53
11.2k
  usize = ABS (usize);
54
11.2k
  vsize = ABS (vsize);
55
56
11.2k
  if (usize < vsize)
57
3.50k
    {
58
3.50k
      MPZ_SRCPTR_SWAP (u, v);
59
3.50k
      MP_SIZE_T_SWAP (usize, vsize);
60
3.50k
    }
61
62
11.2k
  if (vsize == 0)
63
122
    {
64
122
      SIZ (w) = 0;
65
122
      return;
66
122
    }
67
68
11.1k
#if HAVE_NATIVE_mpn_mul_2
69
11.1k
  if (vsize <= 2)
70
3.17k
    {
71
3.17k
      wp = MPZ_REALLOC (w, usize+vsize);
72
3.17k
      if (vsize == 1)
73
2.89k
  cy_limb = mpn_mul_1 (wp, PTR (u), usize, PTR (v)[0]);
74
281
      else
75
281
  {
76
281
    cy_limb = mpn_mul_2 (wp, PTR (u), usize, PTR (v));
77
281
    usize++;
78
281
  }
79
3.17k
      wp[usize] = cy_limb;
80
3.17k
      usize += (cy_limb != 0);
81
3.17k
      SIZ (w) = (sign_product >= 0 ? usize : -usize);
82
3.17k
      return;
83
3.17k
    }
84
#else
85
  if (vsize == 1)
86
    {
87
      wp = MPZ_REALLOC (w, usize+1);
88
      cy_limb = mpn_mul_1 (wp, PTR (u), usize, PTR (v)[0]);
89
      wp[usize] = cy_limb;
90
      usize += (cy_limb != 0);
91
      SIZ (w) = (sign_product >= 0 ? usize : -usize);
92
      return;
93
    }
94
#endif
95
96
7.96k
  TMP_MARK;
97
7.96k
  free_me = NULL;
98
7.96k
  up = PTR (u);
99
7.96k
  vp = PTR (v);
100
7.96k
  wp = PTR (w);
101
102
  /* Ensure W has space enough to store the result.  */
103
7.96k
  wsize = usize + vsize;
104
7.96k
  if (ALLOC (w) < wsize)
105
6.68k
    {
106
6.68k
      if (ALLOC (w) != 0)
107
6.48k
  {
108
6.48k
    if (wp == up || wp == vp)
109
0
      {
110
0
        free_me = wp;
111
0
        free_me_size = ALLOC (w);
112
0
      }
113
6.48k
    else
114
6.48k
      (*__gmp_free_func) (wp, (size_t) ALLOC (w) * GMP_LIMB_BYTES);
115
6.48k
  }
116
117
6.68k
      ALLOC (w) = wsize;
118
6.68k
      wp = __GMP_ALLOCATE_FUNC_LIMBS (wsize);
119
6.68k
      PTR (w) = wp;
120
6.68k
    }
121
1.28k
  else
122
1.28k
    {
123
      /* Make U and V not overlap with W.  */
124
1.28k
      if (wp == up)
125
0
  {
126
    /* W and U are identical.  Allocate temporary space for U.  */
127
0
    up = TMP_ALLOC_LIMBS (usize);
128
    /* Is V identical too?  Keep it identical with U.  */
129
0
    if (wp == vp)
130
0
      vp = up;
131
    /* Copy to the temporary space.  */
132
0
    MPN_COPY (up, wp, usize);
133
0
  }
134
1.28k
      else if (wp == vp)
135
0
  {
136
    /* W and V are identical.  Allocate temporary space for V.  */
137
0
    vp = TMP_ALLOC_LIMBS (vsize);
138
    /* Copy to the temporary space.  */
139
0
    MPN_COPY (vp, wp, vsize);
140
0
  }
141
1.28k
    }
142
143
7.96k
  if (up == vp)
144
2.88k
    {
145
2.88k
      mpn_sqr (wp, up, usize);
146
2.88k
      cy_limb = wp[wsize - 1];
147
2.88k
    }
148
5.07k
  else
149
5.07k
    {
150
5.07k
      cy_limb = mpn_mul (wp, up, usize, vp, vsize);
151
5.07k
    }
152
153
7.96k
  wsize -= cy_limb == 0;
154
155
7.96k
  SIZ (w) = sign_product < 0 ? -wsize : wsize;
156
7.96k
  if (free_me != NULL)
157
0
    (*__gmp_free_func) (free_me, free_me_size * GMP_LIMB_BYTES);
158
7.96k
  TMP_FREE;
159
7.96k
}