1 1.1 mrg /* mpf_mul -- Multiply two floats. 2 1.1 mrg 3 1.1.1.3 mrg Copyright 1993, 1994, 1996, 2001, 2005, 2019 Free Software Foundation, Inc. 4 1.1 mrg 5 1.1 mrg This file is part of the GNU MP Library. 6 1.1 mrg 7 1.1 mrg The GNU MP Library is free software; you can redistribute it and/or modify 8 1.1.1.2 mrg it under the terms of either: 9 1.1.1.2 mrg 10 1.1.1.2 mrg * the GNU Lesser General Public License as published by the Free 11 1.1.1.2 mrg Software Foundation; either version 3 of the License, or (at your 12 1.1.1.2 mrg option) any later version. 13 1.1.1.2 mrg 14 1.1.1.2 mrg or 15 1.1.1.2 mrg 16 1.1.1.2 mrg * the GNU General Public License as published by the Free Software 17 1.1.1.2 mrg Foundation; either version 2 of the License, or (at your option) any 18 1.1.1.2 mrg later version. 19 1.1.1.2 mrg 20 1.1.1.2 mrg or both in parallel, as here. 21 1.1 mrg 22 1.1 mrg The GNU MP Library is distributed in the hope that it will be useful, but 23 1.1 mrg WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY 24 1.1.1.2 mrg or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License 25 1.1.1.2 mrg for more details. 26 1.1 mrg 27 1.1.1.2 mrg You should have received copies of the GNU General Public License and the 28 1.1.1.2 mrg GNU Lesser General Public License along with the GNU MP Library. If not, 29 1.1.1.2 mrg see https://www.gnu.org/licenses/. */ 30 1.1 mrg 31 1.1 mrg #include "gmp-impl.h" 32 1.1 mrg 33 1.1 mrg void 34 1.1 mrg mpf_mul (mpf_ptr r, mpf_srcptr u, mpf_srcptr v) 35 1.1 mrg { 36 1.1 mrg mp_size_t sign_product; 37 1.1 mrg mp_size_t prec = r->_mp_prec; 38 1.1.1.3 mrg mp_size_t rsize; 39 1.1.1.3 mrg mp_limb_t cy_limb; 40 1.1.1.3 mrg mp_ptr rp, tp; 41 1.1.1.3 mrg mp_size_t adj; 42 1.1.1.3 mrg TMP_DECL; 43 1.1 mrg 44 1.1.1.3 mrg if (u == v) 45 1.1 mrg { 46 1.1.1.3 mrg mp_srcptr up; 47 1.1.1.3 mrg mp_size_t usize; 48 1.1.1.3 mrg 49 1.1.1.3 mrg usize = u->_mp_size; 50 1.1.1.3 mrg sign_product = 0; 51 1.1.1.3 mrg 52 1.1.1.3 mrg usize = ABS (usize); 53 1.1.1.3 mrg 54 1.1.1.3 mrg up = u->_mp_d; 55 1.1.1.3 mrg if (usize > prec) 56 1.1.1.3 mrg { 57 1.1.1.3 mrg up += usize - prec; 58 1.1.1.3 mrg usize = prec; 59 1.1.1.3 mrg } 60 1.1.1.3 mrg 61 1.1.1.3 mrg if (usize == 0) 62 1.1.1.3 mrg { 63 1.1.1.3 mrg r->_mp_size = 0; 64 1.1.1.3 mrg r->_mp_exp = 0; /* ??? */ 65 1.1.1.3 mrg return; 66 1.1.1.3 mrg } 67 1.1.1.3 mrg else 68 1.1.1.3 mrg { 69 1.1.1.3 mrg TMP_MARK; 70 1.1.1.3 mrg rsize = 2 * usize; 71 1.1.1.3 mrg tp = TMP_ALLOC_LIMBS (rsize); 72 1.1.1.3 mrg 73 1.1.1.3 mrg mpn_sqr (tp, up, usize); 74 1.1.1.3 mrg cy_limb = tp[rsize - 1]; 75 1.1.1.3 mrg } 76 1.1 mrg } 77 1.1.1.3 mrg else 78 1.1 mrg { 79 1.1.1.3 mrg mp_srcptr up, vp; 80 1.1.1.3 mrg mp_size_t usize, vsize; 81 1.1.1.3 mrg 82 1.1.1.3 mrg usize = u->_mp_size; 83 1.1.1.3 mrg vsize = v->_mp_size; 84 1.1.1.3 mrg sign_product = usize ^ vsize; 85 1.1.1.3 mrg 86 1.1.1.3 mrg usize = ABS (usize); 87 1.1.1.3 mrg vsize = ABS (vsize); 88 1.1.1.3 mrg 89 1.1.1.3 mrg up = u->_mp_d; 90 1.1.1.3 mrg vp = v->_mp_d; 91 1.1.1.3 mrg if (usize > prec) 92 1.1.1.3 mrg { 93 1.1.1.3 mrg up += usize - prec; 94 1.1.1.3 mrg usize = prec; 95 1.1.1.3 mrg } 96 1.1.1.3 mrg if (vsize > prec) 97 1.1.1.3 mrg { 98 1.1.1.3 mrg vp += vsize - prec; 99 1.1.1.3 mrg vsize = prec; 100 1.1.1.3 mrg } 101 1.1.1.3 mrg 102 1.1.1.3 mrg if (usize == 0 || vsize == 0) 103 1.1.1.3 mrg { 104 1.1.1.3 mrg r->_mp_size = 0; 105 1.1.1.3 mrg r->_mp_exp = 0; 106 1.1.1.3 mrg return; 107 1.1.1.3 mrg } 108 1.1.1.3 mrg else 109 1.1.1.3 mrg { 110 1.1.1.3 mrg TMP_MARK; 111 1.1.1.3 mrg rsize = usize + vsize; 112 1.1.1.3 mrg tp = TMP_ALLOC_LIMBS (rsize); 113 1.1.1.3 mrg cy_limb = (usize >= vsize 114 1.1.1.3 mrg ? mpn_mul (tp, up, usize, vp, vsize) 115 1.1.1.3 mrg : mpn_mul (tp, vp, vsize, up, usize)); 116 1.1.1.3 mrg 117 1.1.1.3 mrg } 118 1.1 mrg } 119 1.1 mrg 120 1.1.1.3 mrg adj = cy_limb == 0; 121 1.1.1.3 mrg rsize -= adj; 122 1.1.1.3 mrg prec++; 123 1.1.1.3 mrg if (rsize > prec) 124 1.1 mrg { 125 1.1.1.3 mrg tp += rsize - prec; 126 1.1.1.3 mrg rsize = prec; 127 1.1 mrg } 128 1.1.1.3 mrg rp = r->_mp_d; 129 1.1.1.3 mrg MPN_COPY (rp, tp, rsize); 130 1.1.1.3 mrg r->_mp_exp = u->_mp_exp + v->_mp_exp - adj; 131 1.1.1.3 mrg r->_mp_size = sign_product >= 0 ? rsize : -rsize; 132 1.1.1.2 mrg 133 1.1.1.3 mrg TMP_FREE; 134 1.1 mrg } 135
Indexes created Sun May 03 00:22:47 UTC 2026