xref: /dragonfly/contrib/gmp/mpf/int_p.c (revision d2d4b659)
14b6a78b7SSimon Schubert /* mpf_integer_p -- test whether an mpf is an integer */
24b6a78b7SSimon Schubert 
34b6a78b7SSimon Schubert /*
44b6a78b7SSimon Schubert Copyright 2001, 2002 Free Software Foundation, Inc.
54b6a78b7SSimon Schubert 
64b6a78b7SSimon Schubert This file is part of the GNU MP Library.
74b6a78b7SSimon Schubert 
84b6a78b7SSimon Schubert The GNU MP Library is free software; you can redistribute it and/or modify
94b6a78b7SSimon Schubert it under the terms of the GNU Lesser General Public License as published by
104b6a78b7SSimon Schubert the Free Software Foundation; either version 3 of the License, or (at your
114b6a78b7SSimon Schubert option) any later version.
124b6a78b7SSimon Schubert 
134b6a78b7SSimon Schubert The GNU MP Library is distributed in the hope that it will be useful, but
144b6a78b7SSimon Schubert WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY
154b6a78b7SSimon Schubert or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU Lesser General Public
164b6a78b7SSimon Schubert License for more details.
174b6a78b7SSimon Schubert 
184b6a78b7SSimon Schubert You should have received a copy of the GNU Lesser General Public License
194b6a78b7SSimon Schubert along with the GNU MP Library.  If not, see http://www.gnu.org/licenses/.  */
204b6a78b7SSimon Schubert 
214b6a78b7SSimon Schubert #include "gmp.h"
224b6a78b7SSimon Schubert #include "gmp-impl.h"
234b6a78b7SSimon Schubert 
244b6a78b7SSimon Schubert 
254b6a78b7SSimon Schubert int
mpf_integer_p(mpf_srcptr f)26*d2d4b659SJohn Marino mpf_integer_p (mpf_srcptr f) __GMP_NOTHROW
274b6a78b7SSimon Schubert {
284b6a78b7SSimon Schubert   mp_srcptr ptr;
294b6a78b7SSimon Schubert   mp_exp_t exp;
304b6a78b7SSimon Schubert   mp_size_t size, frac, i;
314b6a78b7SSimon Schubert 
324b6a78b7SSimon Schubert   size = SIZ (f);
334b6a78b7SSimon Schubert   if (size == 0)
344b6a78b7SSimon Schubert     return 1;  /* zero is an integer */
354b6a78b7SSimon Schubert 
364b6a78b7SSimon Schubert   exp = EXP (f);
374b6a78b7SSimon Schubert   if (exp <= 0)
384b6a78b7SSimon Schubert     return 0;  /* has only fraction limbs */
394b6a78b7SSimon Schubert 
404b6a78b7SSimon Schubert   /* any fraction limbs must be zero */
414b6a78b7SSimon Schubert   frac = ABS (size) - exp;
424b6a78b7SSimon Schubert   ptr = PTR (f);
434b6a78b7SSimon Schubert   for (i = 0; i < frac; i++)
444b6a78b7SSimon Schubert     if (ptr[i] != 0)
454b6a78b7SSimon Schubert       return 0;
464b6a78b7SSimon Schubert 
474b6a78b7SSimon Schubert   return 1;
484b6a78b7SSimon Schubert }
49