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

数学计算

开发平台:

Unix_Linux

  1. /* mpq_set_si(dest,ulong_num,ulong_den) -- Set DEST to the rational number
  2.    ULONG_NUM/ULONG_DEN.
  3. Copyright 1991, 1994, 1995, 2001, 2003 Free Software Foundation, Inc.
  4. This file is part of the GNU MP Library.
  5. The GNU MP Library is free software; you can redistribute it and/or modify
  6. it under the terms of the GNU Lesser General Public License as published by
  7. the Free Software Foundation; either version 3 of the License, or (at your
  8. option) any later version.
  9. The GNU MP Library is distributed in the hope that it will be useful, but
  10. WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY
  11. or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU Lesser General Public
  12. License for more details.
  13. You should have received a copy of the GNU Lesser General Public License
  14. along with the GNU MP Library.  If not, see http://www.gnu.org/licenses/.  */
  15. #include "gmp.h"
  16. #include "gmp-impl.h"
  17. void
  18. mpq_set_si (MP_RAT *dest, signed long int num, unsigned long int den)
  19. {
  20.   unsigned long int abs_num;
  21.   if (GMP_NUMB_BITS < BITS_PER_ULONG)
  22.     {
  23.       if (num == 0)  /* Canonicalize 0/d to 0/1.  */
  24.         den = 1;
  25.       mpz_set_si (mpq_numref (dest), num);
  26.       mpz_set_ui (mpq_denref (dest), den);
  27.       return;
  28.     }
  29.   abs_num = ABS (num);
  30.   if (num == 0)
  31.     {
  32.       /* Canonicalize 0/d to 0/1.  */
  33.       den = 1;
  34.       dest->_mp_num._mp_size = 0;
  35.     }
  36.   else
  37.     {
  38.       dest->_mp_num._mp_d[0] = abs_num;
  39.       dest->_mp_num._mp_size = num > 0 ? 1 : -1;
  40.     }
  41.   dest->_mp_den._mp_d[0] = den;
  42.   dest->_mp_den._mp_size = (den != 0);
  43. }