Coverage Report

Created: 2024-11-25 06:29

/src/gmp/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
5.81k
{
39
5.81k
  mp_size_t usize;
40
5.81k
  mp_size_t vsize;
41
5.81k
  mp_size_t wsize;
42
5.81k
  mp_size_t sign_product;
43
5.81k
  mp_ptr up, vp;
44
5.81k
  mp_ptr wp;
45
5.81k
  mp_ptr free_me;
46
5.81k
  size_t free_me_size;
47
5.81k
  mp_limb_t cy_limb;
48
5.81k
  TMP_DECL;
49
50
5.81k
  usize = SIZ (u);
51
5.81k
  vsize = SIZ (v);
52
5.81k
  sign_product = usize ^ vsize;
53
5.81k
  usize = ABS (usize);
54
5.81k
  vsize = ABS (vsize);
55
56
5.81k
  if (usize < vsize)
57
0
    {
58
0
      MPZ_SRCPTR_SWAP (u, v);
59
0
      MP_SIZE_T_SWAP (usize, vsize);
60
0
    }
61
62
5.81k
  if (vsize == 0)
63
12
    {
64
12
      SIZ (w) = 0;
65
12
      return;
66
12
    }
67
68
#if HAVE_NATIVE_mpn_mul_2
69
  if (vsize <= 2)
70
    {
71
      wp = MPZ_REALLOC (w, usize+vsize);
72
      if (vsize == 1)
73
  cy_limb = mpn_mul_1 (wp, PTR (u), usize, PTR (v)[0]);
74
      else
75
  {
76
    cy_limb = mpn_mul_2 (wp, PTR (u), usize, PTR (v));
77
    usize++;
78
  }
79
      wp[usize] = cy_limb;
80
      usize += (cy_limb != 0);
81
      SIZ (w) = (sign_product >= 0 ? usize : -usize);
82
      return;
83
    }
84
#else
85
5.80k
  if (vsize == 1)
86
211
    {
87
211
      wp = MPZ_REALLOC (w, usize+1);
88
211
      cy_limb = mpn_mul_1 (wp, PTR (u), usize, PTR (v)[0]);
89
211
      wp[usize] = cy_limb;
90
211
      usize += (cy_limb != 0);
91
211
      SIZ (w) = (sign_product >= 0 ? usize : -usize);
92
211
      return;
93
211
    }
94
5.59k
#endif
95
96
5.59k
  TMP_MARK;
97
5.59k
  free_me = NULL;
98
5.59k
  up = PTR (u);
99
5.59k
  vp = PTR (v);
100
5.59k
  wp = PTR (w);
101
102
  /* Ensure W has space enough to store the result.  */
103
5.59k
  wsize = usize + vsize;
104
5.59k
  if (ALLOC (w) < wsize)
105
5.59k
    {
106
5.59k
      if (ALLOC (w) != 0)
107
476
  {
108
476
    if (wp == up || wp == vp)
109
476
      {
110
476
        free_me = wp;
111
476
        free_me_size = ALLOC (w);
112
476
      }
113
0
    else
114
0
      (*__gmp_free_func) (wp, (size_t) ALLOC (w) * GMP_LIMB_BYTES);
115
476
  }
116
117
5.59k
      ALLOC (w) = wsize;
118
5.59k
      wp = __GMP_ALLOCATE_FUNC_LIMBS (wsize);
119
5.59k
      PTR (w) = wp;
120
5.59k
    }
121
1
  else
122
1
    {
123
      /* Make U and V not overlap with W.  */
124
1
      if (wp == up)
125
1
  {
126
    /* W and U are identical.  Allocate temporary space for U.  */
127
1
    up = TMP_ALLOC_LIMBS (usize);
128
    /* Is V identical too?  Keep it identical with U.  */
129
1
    if (wp == vp)
130
0
      vp = up;
131
    /* Copy to the temporary space.  */
132
1
    MPN_COPY (up, wp, usize);
133
1
  }
134
0
      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
    }
142
143
5.59k
  if (up == vp)
144
1.04k
    {
145
1.04k
      mpn_sqr (wp, up, usize);
146
1.04k
      cy_limb = wp[wsize - 1];
147
1.04k
    }
148
4.54k
  else
149
4.54k
    {
150
4.54k
      cy_limb = mpn_mul (wp, up, usize, vp, vsize);
151
4.54k
    }
152
153
5.59k
  wsize -= cy_limb == 0;
154
155
5.59k
  SIZ (w) = sign_product < 0 ? -wsize : wsize;
156
5.59k
  if (free_me != NULL)
157
476
    (*__gmp_free_func) (free_me, free_me_size * GMP_LIMB_BYTES);
158
5.59k
  TMP_FREE;
159
5.59k
}