1272343Sngie/* $NetBSD: t_fpclassify.c,v 1.3 2011/10/01 21:47:08 christos Exp $ */ 2272343Sngie 3272343Sngie/*- 4272343Sngie * Copyright (c) 2011 The NetBSD Foundation, Inc. 5272343Sngie * All rights reserved. 6272343Sngie * 7272343Sngie * Redistribution and use in source and binary forms, with or without 8272343Sngie * modification, are permitted provided that the following conditions 9272343Sngie * are met: 10272343Sngie * 1. Redistributions of source code must retain the above copyright 11272343Sngie * notice, this list of conditions and the following disclaimer. 12272343Sngie * 2. Redistributions in binary form must reproduce the above copyright 13272343Sngie * notice, this list of conditions and the following disclaimer in the 14272343Sngie * documentation and/or other materials provided with the distribution. 15272343Sngie * 16272343Sngie * THIS SOFTWARE IS PROVIDED BY THE NETBSD FOUNDATION, INC. AND CONTRIBUTORS 17272343Sngie * ``AS IS'' AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED 18272343Sngie * TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR 19272343Sngie * PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE FOUNDATION OR CONTRIBUTORS 20272343Sngie * BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR 21272343Sngie * CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF 22272343Sngie * SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS 23272343Sngie * INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN 24272343Sngie * CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) 25272343Sngie * ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE 26272343Sngie * POSSIBILITY OF SUCH DAMAGE. 27272343Sngie */ 28272343Sngie 29272343Sngie#include <atf-c.h> 30272343Sngie 31272343Sngie#include <float.h> 32272343Sngie#include <math.h> 33272343Sngie#include <stdio.h> 34272343Sngie#include <string.h> 35272343Sngie 36272343Sngie#ifndef _FLOAT_IEEE754 37272343Sngie 38272343SngieATF_TC(no_test); 39272343SngieATF_TC_HEAD(no_test, tc) 40272343Sngie{ 41272343Sngie atf_tc_set_md_var(tc, "descr", "Dummy test"); 42272343Sngie} 43272343Sngie 44272343SngieATF_TC_BODY(no_test,tc) 45272343Sngie{ 46272343Sngie atf_tc_skip("Test not available on this architecture"); 47272343Sngie} 48272343Sngie 49272343Sngie#else /* defined(_FLOAT_IEEE754) */ 50272343Sngie 51272343SngieATF_TC(fpclassify_float); 52272343SngieATF_TC_HEAD(fpclassify_float, tc) 53272343Sngie{ 54272343Sngie 55272343Sngie atf_tc_set_md_var(tc, "descr", "Test float operations"); 56272343Sngie} 57272343Sngie 58272343SngieATF_TC_BODY(fpclassify_float, tc) 59272343Sngie{ 60272343Sngie float d0, d1, d2, f, ip; 61272343Sngie int e, i; 62272343Sngie 63272343Sngie d0 = FLT_MIN; 64272343Sngie ATF_REQUIRE_EQ(fpclassify(d0), FP_NORMAL); 65272343Sngie f = frexpf(d0, &e); 66272343Sngie ATF_REQUIRE_EQ(e, FLT_MIN_EXP); 67272343Sngie ATF_REQUIRE_EQ(f, 0.5); 68272343Sngie d1 = d0; 69272343Sngie 70272343Sngie /* shift a "1" bit through the mantissa (skip the implicit bit) */ 71272343Sngie for (i = 1; i < FLT_MANT_DIG; i++) { 72272343Sngie d1 /= 2; 73272343Sngie ATF_REQUIRE_EQ(fpclassify(d1), FP_SUBNORMAL); 74272343Sngie ATF_REQUIRE(d1 > 0 && d1 < d0); 75272343Sngie 76272343Sngie d2 = ldexpf(d0, -i); 77272343Sngie ATF_REQUIRE_EQ(d2, d1); 78272343Sngie 79272343Sngie d2 = modff(d1, &ip); 80272343Sngie ATF_REQUIRE_EQ(d2, d1); 81272343Sngie ATF_REQUIRE_EQ(ip, 0); 82272343Sngie 83272343Sngie f = frexpf(d1, &e); 84272343Sngie ATF_REQUIRE_EQ(e, FLT_MIN_EXP - i); 85272343Sngie ATF_REQUIRE_EQ(f, 0.5); 86272343Sngie } 87272343Sngie 88272343Sngie d1 /= 2; 89272343Sngie ATF_REQUIRE_EQ(fpclassify(d1), FP_ZERO); 90272343Sngie f = frexpf(d1, &e); 91272343Sngie ATF_REQUIRE_EQ(e, 0); 92272343Sngie ATF_REQUIRE_EQ(f, 0); 93272343Sngie} 94272343Sngie 95272343SngieATF_TC(fpclassify_double); 96272343SngieATF_TC_HEAD(fpclassify_double, tc) 97272343Sngie{ 98272343Sngie 99272343Sngie atf_tc_set_md_var(tc, "descr", "Test double operations"); 100272343Sngie} 101272343Sngie 102272343SngieATF_TC_BODY(fpclassify_double, tc) 103272343Sngie{ 104272343Sngie double d0, d1, d2, f, ip; 105272343Sngie int e, i; 106272343Sngie 107272343Sngie d0 = DBL_MIN; 108272343Sngie ATF_REQUIRE_EQ(fpclassify(d0), FP_NORMAL); 109272343Sngie f = frexp(d0, &e); 110272343Sngie ATF_REQUIRE_EQ(e, DBL_MIN_EXP); 111272343Sngie ATF_REQUIRE_EQ(f, 0.5); 112272343Sngie d1 = d0; 113272343Sngie 114272343Sngie /* shift a "1" bit through the mantissa (skip the implicit bit) */ 115272343Sngie for (i = 1; i < DBL_MANT_DIG; i++) { 116272343Sngie d1 /= 2; 117272343Sngie ATF_REQUIRE_EQ(fpclassify(d1), FP_SUBNORMAL); 118272343Sngie ATF_REQUIRE(d1 > 0 && d1 < d0); 119272343Sngie 120272343Sngie d2 = ldexp(d0, -i); 121272343Sngie ATF_REQUIRE_EQ(d2, d1); 122272343Sngie 123272343Sngie d2 = modf(d1, &ip); 124272343Sngie ATF_REQUIRE_EQ(d2, d1); 125272343Sngie ATF_REQUIRE_EQ(ip, 0); 126272343Sngie 127272343Sngie f = frexp(d1, &e); 128272343Sngie ATF_REQUIRE_EQ(e, DBL_MIN_EXP - i); 129272343Sngie ATF_REQUIRE_EQ(f, 0.5); 130272343Sngie } 131272343Sngie 132272343Sngie d1 /= 2; 133272343Sngie ATF_REQUIRE_EQ(fpclassify(d1), FP_ZERO); 134272343Sngie f = frexp(d1, &e); 135272343Sngie ATF_REQUIRE_EQ(e, 0); 136272343Sngie ATF_REQUIRE_EQ(f, 0); 137272343Sngie} 138272343Sngie 139272343Sngie/* 140272343Sngie * XXX NetBSD doesn't have long-double flavors of frexp, ldexp, and modf, 141272343Sngie * XXX so this test is disabled. 142272343Sngie */ 143272343Sngie 144272343Sngie#ifdef TEST_LONG_DOUBLE 145272343Sngie 146272343SngieATF_TC(fpclassify_long_double); 147272343SngieATF_TC_HEAD(fpclassify_long_double, tc) 148272343Sngie{ 149272343Sngie 150272343Sngie atf_tc_set_md_var(tc, "descr", "Test long double operations"); 151272343Sngie} 152272343Sngie 153272343SngieATF_TC_BODY(fpclassify_long_double, tc) 154272343Sngie{ 155272343Sngie long double d0, d1, d2, f, ip; 156272343Sngie int e, i; 157272343Sngie 158272343Sngie d0 = LDBL_MIN; 159272343Sngie ATF_REQUIRE_EQ(fpclassify(d0), FP_NORMAL); 160272343Sngie f = frexpl(d0, &e); 161272343Sngie ATF_REQUIRE_EQ(e, LDBL_MIN_EXP); 162272343Sngie ATF_REQUIRE_EQ(f, 0.5); 163272343Sngie d1 = d0; 164272343Sngie 165272343Sngie /* shift a "1" bit through the mantissa (skip the implicit bit) */ 166272343Sngie for (i = 1; i < LDBL_MANT_DIG; i++) { 167272343Sngie d1 /= 2; 168272343Sngie ATF_REQUIRE_EQ(fpclassify(d1), FP_SUBNORMAL); 169272343Sngie ATF_REQUIRE(d1 > 0 && d1 < d0); 170272343Sngie 171272343Sngie d2 = ldexpl(d0, -i); 172272343Sngie ATF_REQUIRE_EQ(d2, d1); 173272343Sngie 174272343Sngie d2 = modfl(d1, &ip); 175272343Sngie ATF_REQUIRE_EQ(d2, d1); 176272343Sngie ATF_REQUIRE_EQ(ip, 0); 177272343Sngie 178272343Sngie f = frexpl(d1, &e); 179272343Sngie ATF_REQUIRE_EQ(e, LDBL_MIN_EXP - i); 180272343Sngie ATF_REQUIRE_EQ(f, 0.5); 181272343Sngie } 182272343Sngie 183272343Sngie d1 /= 2; 184272343Sngie ATF_REQUIRE_EQ(fpclassify(d1), FP_ZERO); 185272343Sngie f = frexpl(d1, &e); 186272343Sngie ATF_REQUIRE_EQ(e, 0); 187272343Sngie ATF_REQUIRE_EQ(f, 0); 188272343Sngie} 189272343Sngie#endif /* TEST_LONG_DOUBLE */ 190272343Sngie#endif /* _FLOAT_IEEE754 */ 191272343Sngie 192272343SngieATF_TP_ADD_TCS(tp) 193272343Sngie{ 194272343Sngie 195272343Sngie#ifndef _FLOAT_IEEE754 196272343Sngie ATF_TP_ADD_TC(tp, no_test); 197272343Sngie#else 198272343Sngie ATF_TP_ADD_TC(tp, fpclassify_float); 199272343Sngie ATF_TP_ADD_TC(tp, fpclassify_double); 200272343Sngie#ifdef TEST_LONG_DOUBLE 201272343Sngie ATF_TP_ADD_TC(tp, fpclassify_long_double); 202272343Sngie#endif /* TEST_LONG_DOUBLE */ 203272343Sngie#endif /* _FLOAT_IEEE754 */ 204272343Sngie 205272343Sngie return atf_no_error(); 206272343Sngie} 207