1 /* mpf_random2 -- Generate a positive random mpf_t of specified size, with
2 long runs of consecutive ones and zeros in the binary representation.
3 Intended for testing of other MP routines.
5 Copyright 1995, 1996, 2001, 2002, 2003 Free Software Foundation, Inc.
7 This file is part of the GNU MP Library.
9 The GNU MP Library is free software; you can redistribute it and/or modify
10 it under the terms of the GNU Lesser General Public License as published by
11 the Free Software Foundation; either version 3 of the License, or (at your
12 option) any later version.
14 The GNU MP Library is distributed in the hope that it will be useful, but
15 WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY
16 or FITNESS FOR A PARTICULAR PURPOSE. See the GNU Lesser General Public
17 License for more details.
19 You should have received a copy of the GNU Lesser General Public License
20 along with the GNU MP Library. If not, see http://www.gnu.org/licenses/. */
27 mpf_random2 (mpf_ptr x
, mp_size_t xs
, mp_exp_t exp
)
46 /* General random mantissa. */
47 mpn_random2 (PTR(x
), xn
);
49 /* Generate random exponent. */
50 _gmp_rand (&elimb
, RANDS
, GMP_NUMB_BITS
);
52 exp
= elimb
% (2 * exp
+ 1) - exp
;
55 SIZ(x
) = xs
< 0 ? -xn
: xn
;