Blob Blame History Raw
/* random.c -- Handle seed for random numbers.

Copyright (C) 2013 INRIA

This file is part of GNU MPC.

GNU MPC 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.

GNU MPC 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 this program. If not, see http://www.gnu.org/licenses/ .
*/

#include "mpc-tests.h"

/* set MPFR flags to random values */
void
set_mpfr_flags (int counter)
{
  if (counter & 1)
    mpfr_set_underflow ();
  else
    mpfr_clear_underflow ();
  if (counter & 2)
    mpfr_set_overflow ();
  else
    mpfr_clear_overflow ();
  /* the divide-by-0 flag was added in MPFR 3.1.0 */
#ifdef mpfr_set_divby0
  if (counter & 4)
    mpfr_set_divby0 ();
  else
    mpfr_clear_divby0 ();
#endif
  if (counter & 8)
    mpfr_set_nanflag ();
  else
    mpfr_clear_nanflag ();
  if (counter & 16)
    mpfr_set_inexflag ();
  else
    mpfr_clear_inexflag ();
  if (counter & 32)
    mpfr_set_erangeflag ();
  else
    mpfr_clear_erangeflag ();
}

/* Check MPFR flags: we allow that some flags are set internally by MPC,
   for example if MPC does internal computations (using MPFR) which yield
   an overflow, even if the final MPC result fits in the exponent range.
   However we don't allow MPC to *clear* the MPFR flags */
void
check_mpfr_flags (int counter)
{
  int old, neu;

  old = (counter & 1) != 0;
  neu = mpfr_underflow_p () != 0;
  if (old && (neu == 0))
    {
      printf ("Error, underflow flag has been modified from %d to %d\n",
              old, neu);
      exit (1);
    }
  old = (counter & 2) != 0;
  neu = mpfr_overflow_p () != 0;
  if (old && (neu == 0))
    {
      printf ("Error, overflow flag has been modified from %d to %d\n",
              old, neu);
      exit (1);
    }
#ifdef mpfr_divby0_p
  old = (counter & 4) != 0;
  neu = mpfr_divby0_p () != 0;
  if (old && (neu == 0))
    {
      printf ("Error, divby0 flag has been modified from %d to %d\n",
              old, neu);
      exit (1);
    }
#endif
  old = (counter & 8) != 0;
  neu = mpfr_nanflag_p () != 0;
  if (old && (neu == 0))
    {
      printf ("Error, nanflag flag has been modified from %d to %d\n",
              old, neu);
      exit (1);
    }
  old = (counter & 16) != 0;
  neu = mpfr_inexflag_p () != 0;
  if (old && (neu == 0))
    {
      printf ("Error, inexflag flag has been modified from %d to %d\n",
              old, neu);
      exit (1);
    }
  old = (counter & 32) != 0;
  neu = mpfr_erangeflag_p () != 0;
  if (old && (neu == 0))
    {
      printf ("Error, erangeflag flag has been modified from %d to %d\n",
              old, neu);
      exit (1);
    }
}