tests/test-printf-frexp.c

   1 /* Test of splitting a double into fraction and mantissa.
   2    Copyright (C) 2007 Free Software Foundation, Inc.
   3
   4    This program is free software; you can redistribute it and/or modify
   5    it under the terms of the GNU General Public License as published by
   6    the Free Software Foundation; either version 2, or (at your option)
   7    any later version.
   8
   9    This program is distributed in the hope that it will be useful,
  10    but WITHOUT ANY WARRANTY; without even the implied warranty of
  11    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
  12    GNU General Public License for more details.
  13
  14    You should have received a copy of the GNU General Public License
  15    along with this program; if not, write to the Free Software Foundation,
  16    Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301, USA.  */
  17
  18 /* Written by Bruno Haible <bruno@clisp.org>, 2007.  */
  19
  20 #include <config.h>
  21
  22 #include "printf-frexp.h"
  23
  24 #include <float.h>
  25 #include <stdlib.h>
  26
  27 #define ASSERT(expr) if (!(expr)) abort ();
  28
  29 static double
  30 my_ldexp (double x, int d)
  31 {
  32   for (; d > 0; d--)
  33     x *= 2.0;
  34   for (; d < 0; d++)
  35     x *= 0.5;
  36   return x;
  37 }
  38
  39 int
  40 main ()
  41 {
  42   int i;
  43   double x;
  44
  45   for (i = 1, x = 1.0; i <= DBL_MAX_EXP; i++, x *= 2.0)
  46     {
  47       int exp = -9999;
  48       double mantissa = printf_frexp (x, &exp);
  49       ASSERT (exp == i - 1);
  50       ASSERT (mantissa == 1.0);
  51     }
  52   for (i = 1, x = 1.0; i >= DBL_MIN_EXP; i--, x *= 0.5)
  53     {
  54       int exp = -9999;
  55       double mantissa = printf_frexp (x, &exp);
  56       ASSERT (exp == i - 1);
  57       ASSERT (mantissa == 1.0);
  58     }
  59   for (; i >= DBL_MIN_EXP - 100 && x > 0.0; i--, x *= 0.5)
  60     {
  61       int exp = -9999;
  62       double mantissa = printf_frexp (x, &exp);
  63       ASSERT (exp == DBL_MIN_EXP - 1);
  64       ASSERT (mantissa == my_ldexp (1.0, i - DBL_MIN_EXP));
  65     }
  66
  67   for (i = 1, x = 1.01; i <= DBL_MAX_EXP; i++, x *= 2.0)
  68     {
  69       int exp = -9999;
  70       double mantissa = printf_frexp (x, &exp);
  71       ASSERT (exp == i - 1);
  72       ASSERT (mantissa == 1.01);
  73     }
  74   for (i = 1, x = 1.01; i >= DBL_MIN_EXP; i--, x *= 0.5)
  75     {
  76       int exp = -9999;
  77       double mantissa = printf_frexp (x, &exp);
  78       ASSERT (exp == i - 1);
  79       ASSERT (mantissa == 1.01);
  80     }
  81   for (; i >= DBL_MIN_EXP - 100 && x > 0.0; i--, x *= 0.5)
  82     {
  83       int exp = -9999;
  84       double mantissa = printf_frexp (x, &exp);
  85       ASSERT (exp == DBL_MIN_EXP - 1);
  86       ASSERT (mantissa >= my_ldexp (1.0, i - DBL_MIN_EXP));
  87       ASSERT (mantissa <= my_ldexp (2.0, i - DBL_MIN_EXP));
  88       ASSERT (mantissa == my_ldexp (x, - exp));
  89     }
  90
  91   for (i = 1, x = 1.73205; i <= DBL_MAX_EXP; i++, x *= 2.0)
  92     {
  93       int exp = -9999;
  94       double mantissa = printf_frexp (x, &exp);
  95       ASSERT (exp == i - 1);
  96       ASSERT (mantissa == 1.73205);
  97     }
  98   for (i = 1, x = 1.73205; i >= DBL_MIN_EXP; i--, x *= 0.5)
  99     {
 100       int exp = -9999;
 101       double mantissa = printf_frexp (x, &exp);
 102       ASSERT (exp == i - 1);
 103       ASSERT (mantissa == 1.73205);
 104     }
 105   for (; i >= DBL_MIN_EXP - 100 && x > 0.0; i--, x *= 0.5)
 106     {
 107       int exp = -9999;
 108       double mantissa = printf_frexp (x, &exp);
 109       ASSERT (exp == DBL_MIN_EXP - 1);
 110       ASSERT (mantissa >= my_ldexp (1.0, i - DBL_MIN_EXP));
 111       ASSERT (mantissa <= my_ldexp (2.0, i - DBL_MIN_EXP));
 112       ASSERT (mantissa == my_ldexp (x, - exp));
 113     }
 114
 115   return 0;
 116 }