/*      $NetBSD: t_strtod.c,v 1.39 2025/04/07 02:28:00 riastradh Exp $ */

/*-
 * Copyright (c) 2011 The NetBSD Foundation, Inc.
 * All rights reserved.
 *
 * This code is derived from software contributed to The NetBSD Foundation
 * by Jukka Ruohonen.
 *
 * Redistribution and use in source and binary forms, with or without
 * modification, are permitted provided that the following conditions
 * are met:
 * 1. Redistributions of source code must retain the above copyright
 *    notice, this list of conditions and the following disclaimer.
 * 2. Redistributions in binary form must reproduce the above copyright
 *    notice, this list of conditions and the following disclaimer in the
 *    documentation and/or other materials provided with the distribution.
 *
 * THIS SOFTWARE IS PROVIDED BY THE NETBSD FOUNDATION, INC. AND CONTRIBUTORS
 * ``AS IS'' AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED
 * TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR
 * PURPOSE ARE DISCLAIMED.  IN NO EVENT SHALL THE FOUNDATION OR CONTRIBUTORS
 * BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR
 * CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF
 * SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS
 * INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN
 * CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE)
 * ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
 * POSSIBILITY OF SUCH DAMAGE.
 */

/* Public domain, Otto Moerbeek <otto@drijf.net>, 2006. */

#include <sys/cdefs.h>
__RCSID("$NetBSD: t_strtod.c,v 1.39 2025/04/07 02:28:00 riastradh Exp $");

#include <errno.h>
#include <fenv.h>
#include <float.h>
#include <math.h>
#include <stdio.h>
#include <stdlib.h>
#include <string.h>

#include <atf-c.h>

static const char * const inf_strings[] =
    { "Inf", "INF", "-Inf", "-INF", "Infinity", "+Infinity",
      "INFINITY", "-INFINITY", "InFiNiTy", "+InFiNiTy" };
const char * const nan_string = "NaN(x)y";

ATF_TC(strtod_basic);
ATF_TC_HEAD(strtod_basic, tc)
{
        atf_tc_set_md_var(tc, "descr", "A basic test of strtod(3)");
}

ATF_TC_BODY(strtod_basic, tc)
{
        static const size_t n = 1024 * 1000;

        for (size_t i = 1; i < n; i = i + 1024) {
                char buf[512];
                (void)snprintf(buf, sizeof(buf), "%zu.%zu", i, i + 1);

                errno = 0;
                double d = strtod(buf, NULL);

                ATF_CHECK_MSG(d > 0, "i=%zu buf=\"%s\" d=%g errno=%d",
                    i, buf, d, errno);
                ATF_CHECK_EQ_MSG(errno, 0, "i=%zu buf=\"%s\" d=%g errno=%d",
                    i, buf, d, errno);
        }
}

ATF_TC(strtold_basic);
ATF_TC_HEAD(strtold_basic, tc)
{
        atf_tc_set_md_var(tc, "descr", "Some examples of strtold(3)");
}

ATF_TC_BODY(strtold_basic, tc)
{
        static const struct {
                const char *str;
                long double val;
        } testcases[] = {
                { "0x0p0", 0.0L },
                { "0x1.234p0", 0x1234 / 4096.0L },
                { "0x1.234p0", 0x1234 / 4096.0L },
#if FLT_RADIX == 2 && LDBL_MAX_EXP == 16384 && LDBL_MANT_DIG == 64
                { "2.16", 0x8.a3d70a3d70a3d71p-2L },
#endif
        };

        for (size_t i = 0, n = __arraycount(testcases); i < n; i++) {
                char *end;
                errno = 0;
                long double val = strtold(testcases[i].str, &end);

                ATF_CHECK_MSG(
                    errno == 0 && *end == '\0' && val == testcases[i].val,
                    "'%s' want %La have %La errno %d end '%s'",
                    testcases[i].str, testcases[i].val, val, errno, end);
        }
}

ATF_TC(strtod_hex);
ATF_TC_HEAD(strtod_hex, tc)
{
        atf_tc_set_md_var(tc, "descr", "A strtod(3) with hexadecimals");
}

ATF_TC_BODY(strtod_hex, tc)
{
        const char *str;
        char *end;
        volatile double d;

        str = "-0x0";
        d = strtod(str, &end);  /* -0.0 */

        ATF_CHECK_EQ_MSG(end, str + 4, "str=%p end=%p", str, end);
        ATF_CHECK_MSG(signbit(d) != 0, "d=%g=%a signbit=%d", d, d, signbit(d));
        ATF_CHECK_EQ_MSG(fabs(d), 0, "d=%g=%a, fabs(d)=%g=%a",
            d, d, fabs(d), fabs(d));

        str = "-0x";
        d = strtod(str, &end);  /* -0.0 */

        ATF_CHECK_EQ_MSG(end, str + 2, "str=%p end=%p", str, end);
        ATF_CHECK_MSG(signbit(d) != 0, "d=%g=%a signbit=%d", d, d, signbit(d));
        ATF_CHECK_EQ_MSG(fabs(d), 0, "d=%g=%a fabs(d)=%g=%a",
            d, d, fabs(d), fabs(d));
}

ATF_TC(strtod_inf);
ATF_TC_HEAD(strtod_inf, tc)
{
        atf_tc_set_md_var(tc, "descr", "A strtod(3) with INF (PR lib/33262)");
}

ATF_TC_BODY(strtod_inf, tc)
{

        if (!isinf(INFINITY))
                atf_tc_skip("no infinities on this architecture");

        for (size_t i = 0; i < __arraycount(inf_strings); i++) {
                volatile double d = strtod(inf_strings[i], NULL);
                ATF_CHECK_MSG(isinf(d), "inf_strings[%zu]=\"%s\" d=%g=%a",
                    i, inf_strings[i], d, d);
        }
}

ATF_TC(strtof_inf);
ATF_TC_HEAD(strtof_inf, tc)
{
        atf_tc_set_md_var(tc, "descr", "A strtof(3) with INF (PR lib/33262)");
}

ATF_TC_BODY(strtof_inf, tc)
{

        if (!isinf(INFINITY))
                atf_tc_skip("no infinities on this architecture");

        for (size_t i = 0; i < __arraycount(inf_strings); i++) {
                volatile float f = strtof(inf_strings[i], NULL);
                ATF_CHECK_MSG(isinf(f), "inf_strings[%zu]=\"%s\" f=%g=%a",
                    i, inf_strings[i], f, f);
                ATF_CHECK_MSG(isinff(f), "inf_strings[%zu]=\"%s\" f=%g=%a",
                    i, inf_strings[i], f, f);
        }
}

ATF_TC(strtold_inf);
ATF_TC_HEAD(strtold_inf, tc)
{
        atf_tc_set_md_var(tc, "descr", "A strtold(3) with INF (PR lib/33262)");
}

ATF_TC_BODY(strtold_inf, tc)
{

        if (!isinf(INFINITY))
                atf_tc_skip("no infinities on this architecture");

        for (size_t i = 0; i < __arraycount(inf_strings); i++) {
                volatile long double ld = strtold(inf_strings[i], NULL);
                ATF_CHECK_MSG(isinf(ld), "inf_strings[%zu]=\"%s\" ld=%Lg=%La",
                    i, inf_strings[i], ld, ld);
        }
}

ATF_TC(strtod_nan);
ATF_TC_HEAD(strtod_nan, tc)
{
        atf_tc_set_md_var(tc, "descr", "A strtod(3) with NaN");
}

ATF_TC_BODY(strtod_nan, tc)
{
        char *end;

#ifndef NAN
        atf_tc_skip("no NaNs on this architecture");
#endif

        volatile double d = strtod(nan_string, &end);
        ATF_CHECK_MSG(isnan(d), "nan_string=\"%s\" d=%g=%a", nan_string, d, d);
        ATF_CHECK_MSG(strcmp(end, "y") == 0, "nan_string=\"%s\"@%p end=%p",
            nan_string, nan_string, end);
}

ATF_TC(strtof_nan);
ATF_TC_HEAD(strtof_nan, tc)
{
        atf_tc_set_md_var(tc, "descr", "A strtof(3) with NaN");
}

ATF_TC_BODY(strtof_nan, tc)
{
        char *end;

#ifndef NAN
        atf_tc_skip("no NaNs on this architecture");
#endif

        volatile float f = strtof(nan_string, &end);
        ATF_CHECK_MSG(isnan(f), "nan_string=\"%s\" f=%g=%a", nan_string, f, f);
        ATF_CHECK_MSG(isnanf(f), "nan_string=\"%s\" f=%g=%a", nan_string,
            f, f);
        ATF_CHECK_MSG(strcmp(end, "y") == 0, "nan_string=\"%s\"@%p end=%p",
            nan_string, nan_string, end);
}

ATF_TC(strtold_nan);
ATF_TC_HEAD(strtold_nan, tc)
{
        atf_tc_set_md_var(tc, "descr", "A strtold(3) with NaN (PR lib/45020)");
}

ATF_TC_BODY(strtold_nan, tc)
{
        char *end;

#ifndef NAN
        atf_tc_skip("no NaNs on this architecture");
#endif

        volatile long double ld = strtold(nan_string, &end);
        ATF_CHECK_MSG(isnan(ld), "nan_string=\"%s\" ld=%Lg=%La",
            nan_string, ld, ld);
        ATF_CHECK_MSG(isnanf(ld), "nan_string=\"%s\" ld=%Lg=%La",
            nan_string, ld, ld);
        ATF_CHECK_MSG(strcmp(end, "y") == 0, "nan_string=\"%s\"@%p end=%p",
            nan_string, nan_string, end);
}

ATF_TC(strtod_round);
ATF_TC_HEAD(strtod_round, tc)
{
        atf_tc_set_md_var(tc, "descr", "Test rounding in strtod(3)");
}

ATF_TC_BODY(strtod_round, tc)
{
#ifdef __HAVE_FENV

        /*
         * Test that strtod(3) honors the current rounding mode.
         * The used value is somewhere near 1 + DBL_EPSILON + FLT_EPSILON.
         */
        const char *val =
            "1.00000011920928977282585492503130808472633361816406";

        (void)fesetround(FE_UPWARD);
        volatile double d1 = strtod(val, NULL);
        (void)fesetround(FE_DOWNWARD);
        volatile double d2 = strtod(val, NULL);
        ATF_CHECK_MSG(d1 > d2, "d1=%g=%a d2=%g=%a", d1, d1, d2, d2);
#else
        atf_tc_skip("Requires <fenv.h> support");
#endif
}

ATF_TC(strtod_underflow);
ATF_TC_HEAD(strtod_underflow, tc)
{
        atf_tc_set_md_var(tc, "descr", "Test underflow in strtod(3)");
}

ATF_TC_BODY(strtod_underflow, tc)
{

        const char *tmp =
            "0.0000000000000000000000000000000000000000000000000000"
            "000000000000000000000000000000000000000000000000000000"
            "000000000000000000000000000000000000000000000000000000"
            "000000000000000000000000000000000000000000000000000000"
            "000000000000000000000000000000000000000000000000000000"
            "000000000000000000000000000000000000000000000000000000"
            "000000000000000000000000000000000000000000000000000000"
            "000000000000000002";

        errno = 0;
        volatile double d = strtod(tmp, NULL);

        if (d != 0 || errno != ERANGE)
                atf_tc_fail("strtod(3) did not detect underflow");
}

/*
 * Bug found by Geza Herman.
 * See
 * http://www.exploringbinary.com/a-bug-in-the-bigcomp-function-of-david-gays-strtod/
 */
ATF_TC(strtod_gherman_bug);
ATF_TC_HEAD(strtod_gherman_bug, tc)
{
        atf_tc_set_md_var(tc, "descr", "Test a bug found by Geza Herman");
}

ATF_TC_BODY(strtod_gherman_bug, tc)
{
        /*
         * Input is _just barely below_ halfway from one binary64
         * (p=53) floating-point number 0x1.d34fd8378ea83p+0 to the
         * next one 0x1.d34fd8378ea84p+0:
         *
         *                                                       v
         * 1.110100110100111111011000001101111000111010101000001101111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111110..._2
         */
        const char *volatile str =
            "1.8254370818746402660437411213933955878019332885742187";
        int error;

        errno = 0;
        volatile double d = strtod(str, NULL);
        error = errno;
        ATF_CHECK_MSG(error == 0, "errno=%d (%s)", error, strerror(error));

#if DBL_MANT_DIG == 53
        ATF_CHECK_EQ_MSG(d, 0x1.d34fd8378ea83p+0, "d=%g=%a", d, d);
#elif DBL_MANT_DIG == 56
        /* a.k.a. 0xe.9a7ec1bc7541cp-3 */
        ATF_CHECK_EQ_MSG(d, 0x1.d34fd8378ea838p+0, "d=%g=%a", d, d);
#else
#  error Unknown DBL_MANT_DIG value!
#endif

#if DBL_MANT_DIG >= 56
        /*
         * Same deal, but VAX D (p=56) between 0xe.9a7ec1bc7541fp-3 and
         * 0xe.9a7ec1bc75420p-3 (a.k.a. 0x1.d34fd8378ea83ep+0 and
         * 0x1.d34fd8378ea840p+0, respectively):
         *
         *                                                          v
         * 1.11010011010011111101100000110111100011101010100000111110111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111110..._2
         *
         * (Not actually sure if this would have triggered the bug!
         * Maybe someone with a VAX can try and find out.)
         */
        str = "1.8254370818746403631882557760945928748697042465209960";

        errno = 0;
        d = strtod(str, NULL);
        error = errno;
        ATF_CHECK_MSG(error == 0, "errno=%d (%s)", error, strerror(error));

        /* a.k.a. 0x1.d34fd8378ea83ep+0 */
        ATF_CHECK_EQ_MSG(d, 0xe.9a7ec1bc7541fp-3, "d=%g=%a", d, d);
#endif
}

ATF_TP_ADD_TCS(tp)
{

        ATF_TP_ADD_TC(tp, strtod_basic);
        ATF_TP_ADD_TC(tp, strtold_basic);
        ATF_TP_ADD_TC(tp, strtod_hex);
        ATF_TP_ADD_TC(tp, strtod_inf);
        ATF_TP_ADD_TC(tp, strtof_inf);
        ATF_TP_ADD_TC(tp, strtold_inf);
        ATF_TP_ADD_TC(tp, strtod_nan);
        ATF_TP_ADD_TC(tp, strtof_nan);
        ATF_TP_ADD_TC(tp, strtold_nan);
        ATF_TP_ADD_TC(tp, strtod_round);
        ATF_TP_ADD_TC(tp, strtod_underflow);
        ATF_TP_ADD_TC(tp, strtod_gherman_bug);

        return atf_no_error();
}