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

数学计算

开发平台:

Unix_Linux

  1. /* zdiv_round() -- divide integers, round to nearest */
  2. /*
  3. Copyright 1999 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. void
  17. zdiv_round (mpz_t rop, mpz_t n, mpz_t d)
  18. {
  19.   mpf_t f_n, f_d;
  20.   mpf_init (f_n);
  21.   mpf_init (f_d);
  22.   mpf_set_z (f_d, d);
  23.   mpf_set_z (f_n, n);
  24.   mpf_div (f_n, f_n, f_d);
  25.   mpf_set_d (f_d, .5);
  26.   if (mpf_sgn (f_n) < 0)
  27.     mpf_neg (f_d, f_d);
  28.   mpf_add (f_n, f_n, f_d);
  29.   mpz_set_f (rop, f_n);
  30.   mpf_clear (f_n);
  31.   mpf_clear (f_d);
  32.   return;
  33. }