Home | History | Annotate | Line # | Download | only in tests 
      1 /* tnorm -- test file for mpc_norm.
      2 
      3 Copyright (C) 2008, 2011, 2012, 2013 INRIA
      4 
      5 This file is part of GNU MPC.
      6 
      7 GNU MPC is free software; you can redistribute it and/or modify it under
      8 the terms of the GNU Lesser General Public License as published by the
      9 Free Software Foundation; either version 3 of the License, or (at your
     10 option) any later version.
     11 
     12 GNU MPC is distributed in the hope that it will be useful, but WITHOUT ANY
     13 WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS
     14 FOR A PARTICULAR PURPOSE. See the GNU Lesser General Public License for
     15 more details.
     16 
     17 You should have received a copy of the GNU Lesser General Public License
     18 along with this program. If not, see http://www.gnu.org/licenses/ .
     19 */
     20 
     21 #include "mpc-tests.h"
     22 
     23 static void
     24 test_underflow (void)
     25 {
     26   mpfr_exp_t emin = mpfr_get_emin ();
     27   mpc_t z;
     28   mpfr_t f;
     29   int inex;
     30 
     31   mpfr_set_emin (-1); /* smallest positive number is 0.5*2^emin = 0.25 */
     32   mpc_init2 (z, 10);
     33   mpfr_set_ui_2exp (mpc_realref (z), 1023, -11, MPFR_RNDN); /* exact */
     34   mpfr_set_ui_2exp (mpc_imagref (z), 1023, -11, MPFR_RNDN); /* exact */
     35   mpfr_init2 (f, 10);
     36 
     37   inex = mpc_norm (f, z, MPFR_RNDZ); /* should give 511/1024 */
     38   if (inex >= 0)
     39     {
     40       printf ("Error in underflow case (1)\n");
     41       printf ("expected inex < 0, got %d\n", inex);
     42       exit (1);
     43     }
     44   if (mpfr_cmp_ui_2exp (f, 511, -10) != 0)
     45     {
     46       printf ("Error in underflow case (1)\n");
     47       printf ("got      ");
     48       mpfr_dump (f);
     49       printf ("expected ");
     50       mpfr_set_ui_2exp (f, 511, -10, MPFR_RNDZ);
     51       mpfr_dump (f);
     52       exit (1);
     53     }
     54 
     55   inex = mpc_norm (f, z, MPFR_RNDN); /* should give 511/1024 */
     56   if (inex >= 0)
     57     {
     58       printf ("Error in underflow case (2)\n");
     59       printf ("expected inex < 0, got %d\n", inex);
     60       exit (1);
     61     }
     62   if (mpfr_cmp_ui_2exp (f, 511, -10) != 0)
     63     {
     64       printf ("Error in underflow case (2)\n");
     65       printf ("got      ");
     66       mpfr_dump (f);
     67       printf ("expected ");
     68       mpfr_set_ui_2exp (f, 511, -10, MPFR_RNDZ);
     69       mpfr_dump (f);
     70       exit (1);
     71     }
     72 
     73   inex = mpc_norm (f, z, MPFR_RNDU); /* should give 1023/2048 */
     74   if (inex <= 0)
     75     {
     76       printf ("Error in underflow case (3)\n");
     77       printf ("expected inex > 0, got %d\n", inex);
     78       exit (1);
     79     }
     80   if (mpfr_cmp_ui_2exp (f, 1023, -11) != 0)
     81     {
     82       printf ("Error in underflow case (3)\n");
     83       printf ("got      ");
     84       mpfr_dump (f);
     85       printf ("expected ");
     86       mpfr_set_ui_2exp (f, 1023, -11, MPFR_RNDZ);
     87       mpfr_dump (f);
     88       exit (1);
     89     }
     90 
     91   mpc_clear (z);
     92   mpfr_clear (f);
     93   mpfr_set_emin (emin);
     94 }
     95 
     96 #define MPC_FUNCTION_CALL                                               \
     97   P[0].mpfr_inex = mpc_norm (P[1].mpfr, P[2].mpc, P[3].mpfr_rnd)
     98 
     99 #include "data_check.tpl"
    100 #include "tgeneric.tpl"
    101 
    102 int
    103 main (void)
    104 {
    105   test_start ();
    106 
    107   data_check_template ("norm.dsc", "norm.dat");
    108 
    109   tgeneric_template ("norm.dsc", 2, 1024, 1, 4096);
    110 
    111   test_underflow ();
    112 
    113   test_end ();
    114 
    115   return 0;
    116 }
    117