vasilito
facf0c92e0
Red Bear OS is a full fork. All sources must be available from git clone with zero network access. Removed gitignore rules that excluded fetched source trees under recipes/*/source/, local/recipes/kde/*/source/, local/recipes/qt/*/source/, and vendor source trees. Build artifacts (target/, build/, source.tar, *.o, *.so) remain excluded. 127291 files added — kernel, relibc, base, bootloader, pkgar, all KDE/Qt frameworks, mesa, wayland, DRM drivers, and every other recipe source.
243 lines
7.2 KiB
C
243 lines
7.2 KiB
C
/* Test file for mpfr_frexp.
|
|
|
|
Copyright 2011-2025 Free Software Foundation, Inc.
|
|
Contributed by the Pascaline and Caramba projects, INRIA.
|
|
|
|
This file is part of the GNU MPFR Library.
|
|
|
|
The GNU MPFR Library is free software; you can redistribute it and/or modify
|
|
it under the terms of the GNU Lesser General Public License as published by
|
|
the Free Software Foundation; either version 3 of the License, or (at your
|
|
option) any later version.
|
|
|
|
The GNU 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 Lesser General Public
|
|
License for more details.
|
|
|
|
You should have received a copy of the GNU Lesser General Public License
|
|
along with the GNU MPFR Library; see the file COPYING.LESSER.
|
|
If not, see <https://www.gnu.org/licenses/>. */
|
|
|
|
#include "mpfr-test.h"
|
|
|
|
static void
|
|
check_special (void)
|
|
{
|
|
mpfr_t x, y;
|
|
int inex;
|
|
mpfr_exp_t exp;
|
|
|
|
mpfr_init2 (x, 53);
|
|
mpfr_init2 (y, 53);
|
|
|
|
mpfr_set_nan (x);
|
|
inex = mpfr_frexp (&exp, y, x, MPFR_RNDN);
|
|
if (mpfr_nan_p (y) == 0 || inex != 0)
|
|
{
|
|
printf ("Error for mpfr_frexp(NaN)\n");
|
|
exit (1);
|
|
}
|
|
|
|
mpfr_set_inf (x, 1);
|
|
inex = mpfr_frexp (&exp, y, x, MPFR_RNDN);
|
|
if (mpfr_inf_p (y) == 0 || mpfr_sgn (y) <= 0 || inex != 0)
|
|
{
|
|
printf ("Error for mpfr_frexp(+Inf)\n");
|
|
exit (1);
|
|
}
|
|
|
|
mpfr_set_inf (x, -1);
|
|
inex = mpfr_frexp (&exp, y, x, MPFR_RNDN);
|
|
if (mpfr_inf_p (y) == 0 || mpfr_sgn (y) >= 0 || inex != 0)
|
|
{
|
|
printf ("Error for mpfr_frexp(-Inf)\n");
|
|
exit (1);
|
|
}
|
|
|
|
mpfr_set_zero (x, 1);
|
|
inex = mpfr_frexp (&exp, y, x, MPFR_RNDN);
|
|
if (mpfr_zero_p (y) == 0 || mpfr_signbit (y) != 0 || inex != 0 || exp != 0)
|
|
{
|
|
printf ("Error for mpfr_frexp(+0)\n");
|
|
exit (1);
|
|
}
|
|
|
|
mpfr_set_zero (x, -1);
|
|
inex = mpfr_frexp (&exp, y, x, MPFR_RNDN);
|
|
if (mpfr_zero_p (y) == 0 || mpfr_signbit (y) == 0 || inex != 0 || exp != 0)
|
|
{
|
|
printf ("Error for mpfr_frexp(-0)\n");
|
|
exit (1);
|
|
}
|
|
|
|
mpfr_set_ui (x, 17, MPFR_RNDN);
|
|
inex = mpfr_frexp (&exp, y, x, MPFR_RNDN);
|
|
/* 17 = 17/32*2^5 */
|
|
if (mpfr_cmp_ui_2exp (y, 17, -5) != 0 || inex != 0 || exp != 5)
|
|
{
|
|
printf ("Error for mpfr_frexp(17)\n");
|
|
exit (1);
|
|
}
|
|
|
|
mpfr_set_si (x, -17, MPFR_RNDN);
|
|
inex = mpfr_frexp (&exp, y, x, MPFR_RNDN);
|
|
if (mpfr_cmp_si_2exp (y, -17, -5) != 0 || inex != 0 || exp != 5)
|
|
{
|
|
printf ("Error for mpfr_frexp(-17)\n");
|
|
exit (1);
|
|
}
|
|
|
|
/* now reduce the precision of y */
|
|
mpfr_set_prec (y, 4);
|
|
mpfr_set_ui (x, 17, MPFR_RNDN);
|
|
inex = mpfr_frexp (&exp, y, x, MPFR_RNDN);
|
|
/* 17 -> 16/32*2^5 */
|
|
if (mpfr_cmp_ui_2exp (y, 16, -5) != 0 || inex >= 0 || exp != 5)
|
|
{
|
|
printf ("Error for mpfr_frexp(17) with prec=4, RNDN\n");
|
|
exit (1);
|
|
}
|
|
|
|
mpfr_set_ui (x, 17, MPFR_RNDN);
|
|
inex = mpfr_frexp (&exp, y, x, MPFR_RNDZ);
|
|
if (mpfr_cmp_ui_2exp (y, 16, -5) != 0 || inex >= 0 || exp != 5)
|
|
{
|
|
printf ("Error for mpfr_frexp(17) with prec=4, RNDZ\n");
|
|
exit (1);
|
|
}
|
|
|
|
mpfr_set_ui (x, 17, MPFR_RNDN);
|
|
inex = mpfr_frexp (&exp, y, x, MPFR_RNDD);
|
|
if (mpfr_cmp_ui_2exp (y, 16, -5) != 0 || inex >= 0 || exp != 5)
|
|
{
|
|
printf ("Error for mpfr_frexp(17) with prec=4, RNDD\n");
|
|
exit (1);
|
|
}
|
|
|
|
mpfr_set_ui (x, 17, MPFR_RNDN);
|
|
inex = mpfr_frexp (&exp, y, x, MPFR_RNDU);
|
|
if (mpfr_cmp_ui_2exp (y, 18, -5) != 0 || inex <= 0 || exp != 5)
|
|
{
|
|
printf ("Error for mpfr_frexp(17) with prec=4, RNDU\n");
|
|
exit (1);
|
|
}
|
|
|
|
mpfr_clear (y);
|
|
mpfr_clear (x);
|
|
}
|
|
|
|
static void check1 (void)
|
|
{
|
|
mpfr_exp_t emin, emax, e;
|
|
mpfr_t x, y1, y2;
|
|
int r, neg, red;
|
|
|
|
emin = mpfr_get_emin ();
|
|
emax = mpfr_get_emax ();
|
|
set_emin (MPFR_EMIN_MIN);
|
|
set_emax (MPFR_EMAX_MAX);
|
|
|
|
mpfr_init2 (x, 7);
|
|
mpfr_inits2 (4, y1, y2, (mpfr_ptr) 0);
|
|
|
|
mpfr_set_ui_2exp (x, 1, -2, MPFR_RNDN);
|
|
while (mpfr_regular_p (x))
|
|
{
|
|
/* Test the exponents up to 3 and with the maximum exponent
|
|
(to check potential intermediate overflow). */
|
|
if (MPFR_GET_EXP (x) == 4)
|
|
mpfr_set_exp (x, MPFR_EMAX_MAX);
|
|
e = MPFR_GET_EXP (x);
|
|
for (neg = 0; neg < 2; neg++)
|
|
{
|
|
RND_LOOP (r)
|
|
{
|
|
int inex1, inex2;
|
|
mpfr_exp_t e1, e2;
|
|
mpfr_flags_t flags1, flags2;
|
|
|
|
for (red = 0; red < 2; red++)
|
|
{
|
|
if (red)
|
|
{
|
|
/* e1: exponent of the rounded value of x. */
|
|
MPFR_ASSERTN (e1 == e || e1 == e + 1);
|
|
set_emin (e);
|
|
set_emax (e);
|
|
mpfr_clear_flags ();
|
|
inex1 = e1 < 0 ?
|
|
mpfr_mul_2ui (y1, x, -e1, (mpfr_rnd_t) r) :
|
|
mpfr_div_2ui (y1, x, e1, (mpfr_rnd_t) r);
|
|
flags1 = __gmpfr_flags;
|
|
}
|
|
else
|
|
{
|
|
inex1 = mpfr_set (y1, x, (mpfr_rnd_t) r);
|
|
e1 = MPFR_IS_INF (y1) ? e + 1 : MPFR_GET_EXP (y1);
|
|
flags1 = inex1 != 0 ? MPFR_FLAGS_INEXACT : 0;
|
|
}
|
|
mpfr_clear_flags ();
|
|
inex2 = mpfr_frexp (&e2, y2, x, (mpfr_rnd_t) r);
|
|
flags2 = __gmpfr_flags;
|
|
set_emin (MPFR_EMIN_MIN);
|
|
set_emax (MPFR_EMAX_MAX);
|
|
if ((!red || e == 0) &&
|
|
(! mpfr_regular_p (y2) || MPFR_GET_EXP (y2) != 0))
|
|
{
|
|
printf ("Error in check1 for %s, red = %d, x = ",
|
|
mpfr_print_rnd_mode ((mpfr_rnd_t) r), red);
|
|
mpfr_dump (x);
|
|
printf ("Expected 1/2 <= |y| < 1, got y = ");
|
|
mpfr_dump (y2);
|
|
exit (1);
|
|
}
|
|
if (!red)
|
|
{
|
|
if (e2 > 0)
|
|
mpfr_mul_2ui (y2, y2, e2, MPFR_RNDN);
|
|
else if (e2 < 0)
|
|
mpfr_div_2ui (y2, y2, -e2, MPFR_RNDN);
|
|
}
|
|
if (! (SAME_SIGN (inex1, inex2) &&
|
|
mpfr_equal_p (y1, y2) &&
|
|
flags1 == flags2))
|
|
{
|
|
printf ("Error in check1 for %s, red = %d, x = ",
|
|
mpfr_print_rnd_mode ((mpfr_rnd_t) r), red);
|
|
mpfr_dump (x);
|
|
printf ("Expected y1 = ");
|
|
mpfr_dump (y1);
|
|
printf ("Got y2 = ");
|
|
mpfr_dump (y2);
|
|
printf ("Expected inex ~= %d, got %d\n", inex1, inex2);
|
|
printf ("Expected flags:");
|
|
flags_out (flags1);
|
|
printf ("Got flags: ");
|
|
flags_out (flags2);
|
|
exit (1);
|
|
}
|
|
}
|
|
}
|
|
mpfr_neg (x, x, MPFR_RNDN);
|
|
}
|
|
mpfr_nextabove (x);
|
|
}
|
|
|
|
mpfr_clears (x, y1, y2, (mpfr_ptr) 0);
|
|
set_emin (emin);
|
|
set_emax (emax);
|
|
}
|
|
|
|
int
|
|
main (int argc, char *argv[])
|
|
{
|
|
tests_start_mpfr ();
|
|
|
|
check_special ();
|
|
check1 ();
|
|
|
|
tests_end_mpfr ();
|
|
return 0;
|
|
}
|