mul_i.h
上传用户:qaz666999
上传日期:2022-08-06
资源大小:2570k
文件大小:3k
源码类别:

数学计算

开发平台:

Unix_Linux

  1. /* mpz_mul_ui/si (product, multiplier, small_multiplicand) -- Set PRODUCT to
  2.    MULTIPLICATOR times SMALL_MULTIPLICAND.
  3. Copyright 1991, 1993, 1994, 1996, 2000, 2001, 2002, 2005, 2008 Free Software
  4. Foundation, Inc.
  5. This file is part of the GNU MP Library.
  6. The GNU MP Library is free software; you can redistribute it and/or modify
  7. it under the terms of the GNU Lesser General Public License as published by
  8. the Free Software Foundation; either version 3 of the License, or (at your
  9. option) any later version.
  10. The GNU MP Library is distributed in the hope that it will be useful, but
  11. WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY
  12. or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU Lesser General Public
  13. License for more details.
  14. You should have received a copy of the GNU Lesser General Public License
  15. along with the GNU MP Library.  If not, see http://www.gnu.org/licenses/.  */
  16. #include "gmp.h"
  17. #include "gmp-impl.h"
  18. #ifdef OPERATION_mul_si
  19. #define FUNCTION               mpz_mul_si
  20. #define MULTIPLICAND_UNSIGNED
  21. #define MULTIPLICAND_ABS(x)    ((unsigned long) ABS(x))
  22. #endif
  23. #ifdef OPERATION_mul_ui
  24. #define FUNCTION               mpz_mul_ui
  25. #define MULTIPLICAND_UNSIGNED  unsigned
  26. #define MULTIPLICAND_ABS(x)    x
  27. #endif
  28. #ifndef FUNCTION
  29. Error, error, unrecognised OPERATION
  30. #endif
  31. void
  32. FUNCTION (mpz_ptr prod, mpz_srcptr mult,
  33.           MULTIPLICAND_UNSIGNED long int small_mult)
  34. {
  35.   mp_size_t size = SIZ(mult);
  36.   mp_size_t sign_product = size;
  37.   mp_limb_t sml;
  38.   mp_limb_t cy;
  39.   mp_ptr pp;
  40.   if (size == 0 || small_mult == 0)
  41.     {
  42.       SIZ(prod) = 0;
  43.       return;
  44.     }
  45.   size = ABS (size);
  46.   sml = MULTIPLICAND_ABS (small_mult);
  47.   if (sml <= GMP_NUMB_MAX)
  48.     {
  49.       MPZ_REALLOC (prod, size + 1);
  50.       pp = PTR(prod);
  51.       cy = mpn_mul_1 (pp, PTR(mult), size, sml & GMP_NUMB_MASK);
  52.       pp[size] = cy;
  53.       size += cy != 0;
  54.     }
  55. #if GMP_NAIL_BITS != 0
  56.   else
  57.     {
  58.       /* Operand too large for the current nails size.  Use temporary for
  59.  intermediate products, to allow prod and mult being identical.  */
  60.       mp_ptr tp;
  61.       TMP_DECL;
  62.       TMP_MARK;
  63.       tp = TMP_ALLOC_LIMBS (size + 2);
  64.       cy = mpn_mul_1 (tp, PTR(mult), size, sml & GMP_NUMB_MASK);
  65.       tp[size] = cy;
  66.       cy = mpn_addmul_1 (tp + 1, PTR(mult), size, sml >> GMP_NUMB_BITS);
  67.       tp[size + 1] = cy;
  68.       size += 2;
  69.       MPN_NORMALIZE_NOT_ZERO (tp, size); /* too general, need to trim one or two limb */
  70.       MPZ_REALLOC (prod, size);
  71.       pp = PTR(prod);
  72.       MPN_COPY (pp, tp, size);
  73.       TMP_FREE;
  74.     }
  75. #endif
  76.   SIZ(prod) = ((sign_product < 0) ^ (small_mult < 0)) ? -size : size;
  77. }