/*-
 * Copyright (C) 2003, 2005 David Schultz <das@FreeBSD.org>
 * All rights reserved.
 *
 * 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 AUTHOR 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 AUTHOR 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.
 */

/*
 * Test for scanf() floating point formats.
 */

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

#include <atf-c.h>

#define eq(type, a, b)  _eq(type##_EPSILON, (a), (b))
static int
_eq(long double epsilon, long double a, long double b)
{
        long double delta;

        delta = fabsl(a - b);
        return (delta <= epsilon);
}

ATF_TC_WITHOUT_HEAD(normalized_numbers);
ATF_TC_BODY(normalized_numbers, tc)
{
        char buf[128];
        long double ld = 0.0;
        double d = 0.0;
        float f = 0.0;

        buf[0] = '\0';
        ATF_REQUIRE(setlocale(LC_NUMERIC, ""));

        ATF_REQUIRE_EQ(1, sscanf("3.141592", "%e", &f));
        ATF_REQUIRE(eq(FLT, f, 3.141592));

        ATF_REQUIRE_EQ(1, sscanf("3.141592653589793", "%lf", &d));
        ATF_REQUIRE(eq(DBL, d, 3.141592653589793));

        ATF_REQUIRE_EQ(1, sscanf("1.234568e+06", "%E", &f));
        ATF_REQUIRE(eq(FLT, f, 1.234568e+06));

        ATF_REQUIRE_EQ(1, sscanf("-1.234568e6", "%lF", &d));
        ATF_REQUIRE(eq(DBL, d, -1.234568e6));

        ATF_REQUIRE_EQ(1, sscanf("+1.234568e-52", "%LG", &ld));
        ATF_REQUIRE(eq(LDBL, ld, 1.234568e-52L));

        ATF_REQUIRE_EQ(1, sscanf("0.1", "%la", &d));
        ATF_REQUIRE(eq(DBL, d, 0.1));

        ATF_REQUIRE_EQ(1, sscanf("00.2", "%lA", &d));
        ATF_REQUIRE(eq(DBL, d, 0.2));

        ATF_REQUIRE_EQ(2, sscanf("123456", "%5le%s", &d, buf));
        ATF_REQUIRE(eq(DBL, d, 12345.));
        ATF_REQUIRE(strcmp(buf, "6") == 0);

        ATF_REQUIRE_EQ(1, sscanf("1.0Q", "%*5le%s", buf));
        ATF_REQUIRE(strcmp(buf, "Q") == 0);

        ATF_REQUIRE_EQ(2, sscanf("-1.23e", "%e%s", &f, buf));
        ATF_REQUIRE(eq(FLT, f, -1.23));
        ATF_REQUIRE(strcmp(buf, "e") == 0);

        ATF_REQUIRE_EQ(2, sscanf("1.25e+", "%le%s", &d, buf));
        ATF_REQUIRE(eq(DBL, d, 1.25));
        ATF_REQUIRE(strcmp(buf, "e+") == 0);

        ATF_REQUIRE_EQ(2, sscanf("1.23E4E5", "%le%s", &d, buf));
        ATF_REQUIRE(eq(DBL, d, 1.23e4));
        ATF_REQUIRE(strcmp(buf, "E5") == 0);

        ATF_REQUIRE_EQ(1, sscanf("12e6", "%le", &d));
        ATF_REQUIRE(eq(DBL, d, 12e6));

        ATF_REQUIRE_EQ(2, sscanf("1.a", "%le%s", &d, buf));
        ATF_REQUIRE(eq(DBL, d, 1.0));
        ATF_REQUIRE(strcmp(buf, "a") == 0);

        ATF_REQUIRE_EQ(2, sscanf(".0p4", "%le%s", &d, buf));
        ATF_REQUIRE(eq(DBL, d, 0.0));
        ATF_REQUIRE(strcmp(buf, "p4") == 0);

        d = 0.25;
        ATF_REQUIRE_EQ(0, sscanf(".", "%le", &d));
        ATF_REQUIRE(d == 0.25);

        ATF_REQUIRE_EQ(1, sscanf("0x08", "%le", &d));
        ATF_REQUIRE(d == 0x8p0);

        ATF_REQUIRE_EQ(2, sscanf("0x90a.bcdefP+09a", "%le%s", &d, buf));
        ATF_REQUIRE(d == 0x90a.bcdefp+09);
        ATF_REQUIRE(strcmp(buf, "a") == 0);

#if (LDBL_MANT_DIG > DBL_MANT_DIG) && !defined(__i386__)
        ATF_REQUIRE_EQ(1, sscanf("3.14159265358979323846", "%Lg", &ld));
        ATF_REQUIRE(eq(LDBL, ld, 3.14159265358979323846L));

        ATF_REQUIRE_EQ(2, sscanf("  0X.0123456789abcdefffp-3g", "%Le%s", &ld, buf));
        ATF_REQUIRE(ld == 0x0.0123456789abcdefffp-3L);
        ATF_REQUIRE(strcmp(buf, "g") == 0);
#endif

        ATF_REQUIRE_EQ(2, sscanf("0xg", "%le%s", &d, buf));
        ATF_REQUIRE(d == 0.0);
        ATF_REQUIRE(strcmp(buf, "xg") == 0);

        ATF_REQUIRE(setlocale(LC_NUMERIC, "ru_RU.ISO8859-5")); /* decimalpoint==, */

        ATF_REQUIRE_EQ(2, sscanf("1.23", "%le%s", &d, buf));
        ATF_REQUIRE(d == 1.0);
        ATF_REQUIRE(strcmp(buf, ".23") == 0);

        ATF_REQUIRE_EQ(1, sscanf("1,23", "%le", &d));
        ATF_REQUIRE(d == 1.23);

        ATF_REQUIRE(setlocale(LC_NUMERIC, ""));
}

ATF_TC_WITHOUT_HEAD(infinities_and_nans);
ATF_TC_BODY(infinities_and_nans, tc)
{
        char buf[128];
        long double ld = 0.0;
        double d = 0.0;
        float f = 0.0;

        ATF_REQUIRE(setlocale(LC_NUMERIC, "C"));

        ATF_REQUIRE_EQ(1, sscanf("-Inf", "%le", &d));
        ATF_REQUIRE(d < 0.0 && isinf(d));

        ATF_REQUIRE_EQ(2, sscanf("iNfInItY and beyond", "%le%s", &d, buf));
        ATF_REQUIRE(d > 0.0 && isinf(d));
        ATF_REQUIRE(strcmp(buf, " and beyond"));

        ATF_REQUIRE_EQ(1, sscanf("NaN", "%le", &d));
        ATF_REQUIRE(isnan(d));

        ATF_REQUIRE_EQ(2, sscanf("NAN(123Y", "%le%s", &d, buf));
        ATF_REQUIRE(isnan(d));
        ATF_REQUIRE(strcmp(buf, "(123Y") == 0);

        ATF_REQUIRE_EQ(2, sscanf("nan(f00f)plugh", "%le%s", &d, buf));
        ATF_REQUIRE(isnan(d));
        ATF_REQUIRE(strcmp(buf, "plugh") == 0);

        ATF_REQUIRE_EQ(1, sscanf("-nan", "%le", &d));
        ATF_REQUIRE(isnan(d));

        /* Only quiet NaNs should be returned. */
        ATF_REQUIRE_EQ(1, sscanf("NaN", "%e", &f));
        ATF_REQUIRE_EQ(1, sscanf("nan", "%le", &d));
        ATF_REQUIRE_EQ(1, sscanf("nan", "%Le", &ld));
        ATF_REQUIRE_EQ(0, feclearexcept(FE_ALL_EXCEPT));
        ATF_REQUIRE(f != f);
        ATF_REQUIRE(d != d);
        ATF_REQUIRE(ld != ld);
        ATF_REQUIRE(fetestexcept(FE_INVALID) == 0);
        ATF_REQUIRE_EQ(1, sscanf("nan(1234)", "%e", &f));
        ATF_REQUIRE_EQ(1, sscanf("nan(1234)", "%le", &d));
        ATF_REQUIRE_EQ(1, sscanf("nan(1234)", "%Le", &ld));
        ATF_REQUIRE_EQ(0, feclearexcept(FE_ALL_EXCEPT));
        ATF_REQUIRE(f != f);
        ATF_REQUIRE(d != d);
        ATF_REQUIRE(ld != ld);
        /* POSIX says we should only generate quiet NaNs. */
        ATF_REQUIRE(fetestexcept(FE_INVALID) == 0);
}

ATF_TC_WITHOUT_HEAD(rounding_tests);
ATF_TC_BODY(rounding_tests, tc)
{
        long double ld = 0.0;
        double d = 0.0;

        ATF_REQUIRE(setlocale(LC_NUMERIC, "C"));

        ATF_REQUIRE_EQ(0, fesetround(FE_DOWNWARD));

        ATF_REQUIRE_EQ(1, sscanf("1.999999999999999999999999999999999", "%le", &d));
        ATF_REQUIRE(d < 2.0);
        ATF_REQUIRE_EQ(1, sscanf("0x1.ffffffffffffffp0", "%le", &d));
        ATF_REQUIRE(d < 2.0);
        ATF_REQUIRE_EQ(1, sscanf("1.999999999999999999999999999999999", "%Le", &ld));
        ATF_REQUIRE(ld < 2.0);

        ATF_REQUIRE_EQ(1, sscanf("1.0571892669084007", "%le", &d));
        ATF_REQUIRE(d == 0x1.0ea3f4af0dc59p0);
        ATF_REQUIRE_EQ(1, sscanf("-1.0571892669084007", "%le", &d));
        ATF_REQUIRE(d == -0x1.0ea3f4af0dc5ap0);
        ATF_REQUIRE_EQ(1, sscanf("1.0571892669084010", "%le", &d));
        ATF_REQUIRE(d == 0x1.0ea3f4af0dc5ap0);

        ATF_REQUIRE_EQ(1, sscanf("0x1.23p-5000", "%le", &d));
        ATF_REQUIRE(d == 0.0);

        ATF_REQUIRE_EQ(1, sscanf("0x1.2345678p-1050", "%le", &d));
        ATF_REQUIRE(d == 0x1.234567p-1050);

        ATF_REQUIRE_EQ(0, fesetround(FE_UPWARD));

        ATF_REQUIRE_EQ(1, sscanf("1.0571892669084007", "%le", &d));
        ATF_REQUIRE(d == 0x1.0ea3f4af0dc5ap0);
        ATF_REQUIRE_EQ(1, sscanf("-1.0571892669084007", "%le", &d));
        ATF_REQUIRE(d == -0x1.0ea3f4af0dc59p0);
        ATF_REQUIRE_EQ(1, sscanf("1.0571892669084010", "%le", &d));
        ATF_REQUIRE(d == 0x1.0ea3f4af0dc5bp0);

        ATF_REQUIRE_EQ(1, sscanf("0x1.23p-5000", "%le", &d));
        ATF_REQUIRE(d == 0x1p-1074);

        ATF_REQUIRE_EQ(1, sscanf("0x1.2345678p-1050", "%le", &d));
        ATF_REQUIRE(d == 0x1.234568p-1050);

        ATF_REQUIRE_EQ(0, fesetround(FE_TOWARDZERO));

        ATF_REQUIRE_EQ(1, sscanf("1.0571892669084007", "%le", &d));
        ATF_REQUIRE(d == 0x1.0ea3f4af0dc59p0);
        ATF_REQUIRE_EQ(1, sscanf("-1.0571892669084007", "%le", &d));
        ATF_REQUIRE(d == -0x1.0ea3f4af0dc59p0);
        ATF_REQUIRE_EQ(1, sscanf("1.0571892669084010", "%le", &d));
        ATF_REQUIRE(d == 0x1.0ea3f4af0dc5ap0);

        ATF_REQUIRE_EQ(1, sscanf("0x1.23p-5000", "%le", &d));
        ATF_REQUIRE(d == 0.0);

        ATF_REQUIRE_EQ(1, sscanf("0x1.2345678p-1050", "%le", &d));
        ATF_REQUIRE(d == 0x1.234567p-1050);

        ATF_REQUIRE_EQ(0, fesetround(FE_TONEAREST));

        /* 1.0571892669084007 is slightly closer to 0x1.0ea3f4af0dc59p0 */
        ATF_REQUIRE_EQ(1, sscanf("1.0571892669084007", "%le", &d));
        ATF_REQUIRE(d == 0x1.0ea3f4af0dc59p0);
        ATF_REQUIRE_EQ(1, sscanf("-1.0571892669084007", "%le", &d));
        ATF_REQUIRE(d == -0x1.0ea3f4af0dc59p0);
        ATF_REQUIRE_EQ(1, sscanf("1.0571892669084010", "%le", &d));
        ATF_REQUIRE(d == 0x1.0ea3f4af0dc5bp0);

        /* strtod() should round small numbers to 0. */
        ATF_REQUIRE_EQ(1, sscanf("0x1.23p-5000", "%le", &d));
        ATF_REQUIRE(d == 0.0);

        /* Extra digits in a denormal shouldn't break anything. */
        ATF_REQUIRE_EQ(1, sscanf("0x1.2345678p-1050", "%le", &d));
        ATF_REQUIRE(d == 0x1.234568p-1050);
}

ATF_TC_WITHOUT_HEAD(strtod);
ATF_TC_BODY(strtod, tc)
{
        char *endp;

        ATF_REQUIRE(setlocale(LC_NUMERIC, "C"));

        ATF_REQUIRE(strtod("0xy", &endp) == 0);
        ATF_REQUIRE(strcmp("xy", endp) == 0);

        /* This used to cause an infinite loop and round the wrong way. */
        ATF_REQUIRE_EQ(0, fesetround(FE_DOWNWARD));
        ATF_REQUIRE(strtof("3.5e38", &endp) == FLT_MAX);
        ATF_REQUIRE(strtod("2e308", &endp) == DBL_MAX);
        ATF_REQUIRE_EQ(0, fesetround(FE_UPWARD));
        ATF_REQUIRE(strtof("3.5e38", &endp) == INFINITY);
        ATF_REQUIRE(strtod("2e308", &endp) == INFINITY);
        ATF_REQUIRE_EQ(0, fesetround(FE_TOWARDZERO));
        ATF_REQUIRE(strtof("3.5e38", &endp) == FLT_MAX);
        ATF_REQUIRE(strtod("2e308", &endp) == DBL_MAX);
        ATF_REQUIRE_EQ(0, fesetround(FE_TONEAREST));
        ATF_REQUIRE(strtof("3.5e38", &endp) == INFINITY);
        ATF_REQUIRE(strtod("2e308", &endp) == INFINITY);
}

ATF_TP_ADD_TCS(tp)
{

        ATF_TP_ADD_TC(tp, normalized_numbers);
        ATF_TP_ADD_TC(tp, infinities_and_nans);
        ATF_TP_ADD_TC(tp, rounding_tests);
        ATF_TP_ADD_TC(tp, strtod);

        return (atf_no_error());
}