tools/regression/lib/msun/test-ilogb.c

   1 /*-
   2  * Copyright (c) 2004 Stefan Farfeleder
   3  * All rights reserved.
   4  *
   5  * Redistribution and use in source and binary forms, with or without
   6  * modification, are permitted provided that the following conditions
   7  * are met:
   8  * 1. Redistributions of source code must retain the above copyright
   9  *    notice, this list of conditions and the following disclaimer.
  10  * 2. Redistributions in binary form must reproduce the above copyright
  11  *    notice, this list of conditions and the following disclaimer in the
  12  *    documentation and/or other materials provided with the distribution.
  13  *
  14  * THIS SOFTWARE IS PROVIDED BY THE AUTHOR AND CONTRIBUTORS ``AS IS'' AND
  15  * ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
  16  * IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
  17  * ARE DISCLAIMED.  IN NO EVENT SHALL THE AUTHOR OR CONTRIBUTORS BE LIABLE
  18  * FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL
  19  * DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS
  20  * OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION)
  21  * HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT
  22  * LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY
  23  * OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF
  24  * SUCH DAMAGE.
  25  *
  26  * $FreeBSD$
  27  */
  28
  29 #include <assert.h>
  30 #include <float.h>
  31 #include <limits.h>
  32 #include <math.h>
  33 #include <stdio.h>
  34 #include <stdlib.h>
  35
  36 int
  37 main(void)
  38 {
  39         char buf[128], *end;
  40         double d;
  41         float f;
  42         long double ld;
  43         int e, i;
  44
  45         printf("1..3\n");
  46         assert(ilogb(0) == FP_ILOGB0);
  47         assert(ilogb(NAN) == FP_ILOGBNAN);
  48         assert(ilogb(INFINITY) == INT_MAX);
  49         for (e = DBL_MIN_EXP - DBL_MANT_DIG; e < DBL_MAX_EXP; e++) {
  50                 snprintf(buf, sizeof(buf), "0x1.p%d", e);
  51                 d = strtod(buf, &end);
  52                 assert(*end == '\0');
  53                 i = ilogb(d);
  54                 assert(i == e);
  55         }
  56         printf("ok 1 - ilogb\n");
  57
  58         assert(ilogbf(0) == FP_ILOGB0);
  59         assert(ilogbf(NAN) == FP_ILOGBNAN);
  60         assert(ilogbf(INFINITY) == INT_MAX);
  61         for (e = FLT_MIN_EXP - FLT_MANT_DIG; e < FLT_MAX_EXP; e++) {
  62                 snprintf(buf, sizeof(buf), "0x1.p%d", e);
  63                 f = strtof(buf, &end);
  64                 assert(*end == '\0');
  65                 i = ilogbf(f);
  66                 assert(i == e);
  67         }
  68         printf("ok 2 - ilogbf\n");
  69
  70         assert(ilogbl(0) == FP_ILOGB0);
  71         assert(ilogbl(NAN) == FP_ILOGBNAN);
  72         assert(ilogbl(INFINITY) == INT_MAX);
  73         for (e = LDBL_MIN_EXP - LDBL_MANT_DIG; e < LDBL_MAX_EXP; e++) {
  74                 snprintf(buf, sizeof(buf), "0x1.p%d", e);
  75                 ld = strtold(buf, &end);
  76                 assert(*end == '\0');
  77                 i = ilogbl(ld);
  78                 assert(i == e);
  79         }
  80         printf("ok 3 - ilogbl\n");
  81
  82         return (0);
  83 }