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 Marinompf_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