tests/lib/libc/gen/t_fpsetround.c

   1 /* $NetBSD: t_fpsetround.c,v 1.6 2011/10/01 17:46:10 christos Exp $ */
   2
   3 /*-
   4  * Copyright (c) 2011 The NetBSD Foundation, Inc.
   5  * All rights reserved.
   6  *
   7  * This code is derived from software contributed to The NetBSD Foundation
   8  * by Christos Zoulas.
   9  *
  10  * Redistribution and use in source and binary forms, with or without
  11  * modification, are permitted provided that the following conditions
  12  * are met:
  13  * 1. Redistributions of source code must retain the above copyright
  14  *    notice, this list of conditions and the following disclaimer.
  15  * 2. Redistributions in binary form must reproduce the above copyright
  16  *    notice, this list of conditions and the following disclaimer in the
  17  *    documentation and/or other materials provided with the distribution.
  18  * 3. All advertising materials mentioning features or use of this software
  19  *    must display the following acknowledgement:
  20  *        This product includes software developed by the NetBSD
  21  *        Foundation, Inc. and its contributors.
  22  * 4. Neither the name of The NetBSD Foundation nor the names of its
  23  *    contributors may be used to endorse or promote products derived
  24  *    from this software without specific prior written permission.
  25  *
  26  * THIS SOFTWARE IS PROVIDED BY THE NETBSD FOUNDATION, INC. AND CONTRIBUTORS
  27  * ``AS IS'' AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED
  28  * TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR
  29  * PURPOSE ARE DISCLAIMED.  IN NO EVENT SHALL THE FOUNDATION OR CONTRIBUTORS
  30  * BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR
  31  * CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF
  32  * SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS
  33  * INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN
  34  * CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE)
  35  * ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
  36  * POSSIBILITY OF SUCH DAMAGE.
  37  */
  38 #include <sys/cdefs.h>
  39 __RCSID("$NetBSD: t_fpsetround.c,v 1.6 2011/10/01 17:46:10 christos Exp $");
  40
  41 #include <float.h>
  42 #include <math.h>
  43 #include <stdlib.h>
  44 #include <string.h>
  45 #include <stdio.h>
  46
  47 #include <atf-c.h>
  48
  49 ATF_TC(fpsetround_basic);
  50 ATF_TC_HEAD(fpsetround_basic, tc)
  51 {
  52
  53         atf_tc_set_md_var(tc, "descr",
  54             "Minimal testing of fpgetround(3) and fpsetround(3)");
  55 }
  56
  57 #ifdef _FLOAT_IEEE754
  58 #include <ieeefp.h>
  59
  60 static const struct {
  61         const char *n;
  62         int rm;
  63         int rf;
  64 } rnd[] = {
  65         { "RN", FP_RN, 1 },
  66         { "RP", FP_RP, 2 },
  67         { "RM", FP_RM, 3 },
  68         { "RZ", FP_RZ, 0 },
  69
  70 };
  71
  72 static const struct {
  73         const char *n;
  74         int v[4];
  75 } tst[] = {     /*  RN  RP  RM  RZ */
  76         {  "1.1", {  1,  1,  2,  1 } },
  77         {  "1.5", {  1,  2,  2,  1 } },
  78         {  "1.9", {  1,  2,  2,  1 } },
  79         {  "2.1", {  2,  2,  3,  2 } },
  80         {  "2.5", {  2,  2,  3,  2 } },
  81         {  "2.9", {  2,  3,  3,  2 } },
  82         { "-1.1", { -1, -1, -1, -2 } },
  83         { "-1.5", { -1, -2, -1, -2 } },
  84         { "-1.9", { -1, -2, -1, -2 } },
  85         { "-2.1", { -2, -2, -2, -3 } },
  86         { "-2.5", { -2, -2, -2, -3 } },
  87         { "-2.9", { -2, -3, -2, -3 } },
  88 };
  89
  90 static const char *
  91 getname(int r)
  92 {
  93         for (size_t i = 0; i < __arraycount(rnd); i++)
  94                 if (rnd[i].rm == r)
  95                         return rnd[i].n;
  96         return "*unknown*";
  97 }
  98
  99 static void
 100 test(int r)
 101 {
 102         int did = 0;
 103         for (size_t i = 0; i < __arraycount(tst); i++) {
 104                 double d = strtod(tst[i].n, NULL);
 105                 int g = (int)rint(d);
 106                 int e = tst[i].v[r];
 107                 ATF_CHECK_EQ(g, e);
 108                 if (g != e) {
 109                         if (!did) {
 110                                 fprintf(stderr, "Mode Value Result Expected\n");
 111                                 did = 1;
 112                         }
 113                         fprintf(stderr, "%4.4s %-5.5s %6d %8d\n", rnd[r].n,
 114                             tst[i].n, (int)rint(d), tst[i].v[r]);
 115                 }
 116         }
 117 }
 118 #endif
 119
 120
 121 ATF_TC_BODY(fpsetround_basic, tc)
 122 {
 123
 124 #ifndef _FLOAT_IEEE754
 125         atf_tc_skip("Test not applicable on this architecture.");
 126 #else
 127         int r;
 128
 129         ATF_CHECK_EQ(r = fpgetround(), FP_RN);
 130         if (FP_RN != r)
 131                 fprintf(stderr, "default expected=%s got=%s\n", getname(FP_RN),
 132                     getname(r));
 133         ATF_CHECK_EQ(FLT_ROUNDS, 1);
 134
 135         for (size_t i = 0; i < __arraycount(rnd); i++) {
 136                 const size_t j = (i + 1) & 3;
 137                 const int o = rnd[i].rm;
 138                 const int n = rnd[j].rm;
 139
 140                 ATF_CHECK_EQ(r = fpsetround(n), o);
 141                 if (o != r)
 142                         fprintf(stderr, "set %s expected=%s got=%s\n",
 143                             getname(n), getname(o), getname(r));
 144                 ATF_CHECK_EQ(r = fpgetround(), n);
 145                 if (n != r)
 146                         fprintf(stderr, "get expected=%s got=%s\n", getname(n),
 147                             getname(r));
 148                 ATF_CHECK_EQ(r = FLT_ROUNDS, rnd[j].rf);
 149                 if (r != rnd[j].rf)
 150                         fprintf(stderr, "rounds expected=%x got=%x\n",
 151                             rnd[j].rf, r);
 152                 test(r);
 153         }
 154 #endif /* _FLOAT_IEEE754 */
 155 }
 156
 157 ATF_TP_ADD_TCS(tp)
 158 {
 159
 160         ATF_TP_ADD_TC(tp, fpsetround_basic);
 161
 162         return atf_no_error();
 163 }