[BACK]Return to set_si.c CVS log [TXT][DIR] Up to [local] / OpenXM_contrib / gmp / mpfr

File: [local] / OpenXM_contrib / gmp / mpfr / Attic / set_si.c (download)

Revision 1.1.1.1 (vendor branch), Sat Sep 9 14:12:19 2000 UTC (23 years, 8 months ago) by maekawa
Branch: GMP
CVS Tags: maekawa-ipv6, VERSION_3_1_1, VERSION_3_1, RELEASE_1_2_2, RELEASE_1_2_1, RELEASE_1_1_3
Changes since 1.1: +0 -0 lines

Import gmp 3.1

/* mpfr_set_si, mpfr_set_ui -- set a MPFR number from a machine integer

Copyright (C) 1999 PolKA project, Inria Lorraine and Loria

This file is part of the MPFR Library.

The MPFR Library is free software; you can redistribute it and/or modify
it under the terms of the GNU Library General Public License as published by
the Free Software Foundation; either version 2 of the License, or (at your
option) any later version.

The MPFR Library is distributed in the hope that it will be useful, but
WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY
or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU Library General Public
License for more details.

You should have received a copy of the GNU Library General Public License
along with the MPFR Library; see the file COPYING.LIB.  If not, write to
the Free Software Foundation, Inc., 59 Temple Place - Suite 330, Boston,
MA 02111-1307, USA. */

#include <stdio.h>
#include "gmp.h"
#include "gmp-impl.h"
#include "longlong.h"
#include "mpfr.h"

void
#if __STDC__
mpfr_set_si(mpfr_ptr x, long i, unsigned char rnd_mode)
#else
mpfr_set_si(x, i, rnd_mode)
     mpfr_ptr x;
     long i;
     unsigned char rnd_mode;
#endif
{
  unsigned long xn, cnt; mp_limb_t ai;

  if (i==0) { SET_ZERO(x); return; }
  xn = (PREC(x)-1)/BITS_PER_MP_LIMB;
  ai = ABS(i); 

  count_leading_zeros(cnt, ai); 

  x -> _mp_d[xn] = ai << cnt;
  /* don't forget to put zero in lower limbs */
  MPN_ZERO(MANT(x), xn);
  x -> _mp_exp = BITS_PER_MP_LIMB - cnt;
  /* warning: don't change the precision of x! */
  if (i*SIGN(x) < 0) CHANGE_SIGN(x);

  return; 
}

void
#if __STDC__
mpfr_set_ui(mpfr_ptr x, unsigned long i, unsigned char rnd_mode)
#else
mpfr_set_ui(x, i, rnd_mode)
     mpfr_ptr x;
     long i;
     unsigned char rnd_mode;
#endif  
{
  unsigned int xn, cnt;

  if (i==0) { SET_ZERO(x); return; }
  xn = (PREC(x)-1)/BITS_PER_MP_LIMB;
  count_leading_zeros(cnt, (mp_limb_t) i); 

  x -> _mp_d[xn] = ((mp_limb_t) i) << cnt; 
  /* don't forget to put zero in lower limbs */
  MPN_ZERO(MANT(x), xn);
  x -> _mp_exp = BITS_PER_MP_LIMB - cnt;
  /* warning: don't change the precision of x! */
  if (SIGN(x) < 0) CHANGE_SIGN(x);

  return; 
}