mirrors/cosmopolitan
1
0
Fork 0
mirror of https://github.com/jart/cosmopolitan.git synced 2025-05-23 05:42:29 +00:00
Code Issues Projects Releases Packages Wiki Activity Actions 1

Use ARM's faster math functions on non-tiny builds

Browse source
This commit is contained in:
Justine Tunney 2022-07-11 18:34:10 -07:00
parent 3c10fb5580
commit 4814b6bdf8
58 changed files with 3760 additions and 361 deletions
Download patch file Download diff file Expand all files Collapse all files

10
test/libc/tinymath/ceil_test.c
View file

@ -18,6 +18,7 @@
╚─────────────────────────────────────────────────────────────────────────────*/
#include "libc/math.h"
#include "libc/runtime/gc.internal.h"
#include "libc/testlib/ezbench.h"
#include "libc/testlib/testlib.h"
#include "libc/x/x.h"
@ -53,3 +54,12 @@ TEST(ceill, test) {
EXPECT_STREQ("INFINITY", gc(xdtoal(ceill(INFINITY))));
EXPECT_STREQ("-INFINITY", gc(xdtoal(ceill(-INFINITY))));
}
BENCH(ceill, bench) {
double _ceil(double) asm("ceil");
float _ceilf(float) asm("ceilf");
long double _ceill(long double) asm("ceill");
EZBENCH2("ceil", donothing, _ceil(.7)); /* ~3ns */
EZBENCH2("ceilf", donothing, _ceilf(.7)); /* ~3ns */
EZBENCH2("ceill", donothing, _ceill(.7)); /* ~9ns */
}

10
test/libc/tinymath/exp2_test.c
View file

@ -18,6 +18,7 @@
╚─────────────────────────────────────────────────────────────────────────────*/
#include "libc/math.h"
#include "libc/runtime/gc.internal.h"
#include "libc/testlib/ezbench.h"
#include "libc/testlib/testlib.h"
#include "libc/x/x.h"
@ -54,3 +55,12 @@ TEST(exp2f, test) {
EXPECT_STREQ("0", gc(xdtoaf(exp2f(-132098844872390))));
EXPECT_STREQ("INFINITY", gc(xdtoaf(exp2f(132098844872390))));
}
BENCH(exp2l, bench) {
double _exp2(double) asm("exp2");
float _exp2f(float) asm("exp2f");
long double _exp2l(long double) asm("exp2l");
EZBENCH2("exp2", donothing, _exp2(.7)); /* ~6ns */
EZBENCH2("exp2f", donothing, _exp2f(.7)); /* ~5ns */
EZBENCH2("exp2l", donothing, _exp2l(.7)); /* ~30ns */
}

10
test/libc/tinymath/exp_test.c
View file

@ -19,6 +19,7 @@
#include "libc/math.h"
#include "libc/runtime/gc.internal.h"
#include "libc/stdio/stdio.h"
#include "libc/testlib/ezbench.h"
#include "libc/testlib/testlib.h"
#include "libc/x/x.h"
@ -61,3 +62,12 @@ TEST(exp, fun) {
ASSERT_STREQ("6.389056", gc(xasprintf("%f", expm1(2.0))));
ASSERT_STREQ("6.389056", gc(xasprintf("%f", exp(2.0) - 1.0)));
}
BENCH(expl, bench) {
double _exp(double) asm("exp");
float _expf(float) asm("expf");
long double _expl(long double) asm("expl");
EZBENCH2("exp", donothing, _exp(.7)); /* ~6ns */
EZBENCH2("expf", donothing, _expf(.7)); /* ~5ns */
EZBENCH2("expl", donothing, _expl(.7)); /* ~20ns */
}

10
test/libc/tinymath/floor_test.c
View file

@ -18,6 +18,7 @@
╚─────────────────────────────────────────────────────────────────────────────*/
#include "libc/math.h"
#include "libc/runtime/gc.internal.h"
#include "libc/testlib/ezbench.h"
#include "libc/testlib/testlib.h"
#include "libc/x/x.h"
@ -53,3 +54,12 @@ TEST(floorl, test) {
EXPECT_STREQ("INFINITY", gc(xdtoal(floorl(INFINITY))));
EXPECT_STREQ("-INFINITY", gc(xdtoal(floorl(-INFINITY))));
}
BENCH(floorl, bench) {
double _floor(double) asm("floor");
float _floorf(float) asm("floorf");
long double _floorl(long double) asm("floorl");
EZBENCH2("floor", donothing, _floor(.7)); /* ~3ns */
EZBENCH2("floorf", donothing, _floorf(.7)); /* ~3ns */
EZBENCH2("floorl", donothing, _floorl(.7)); /* ~9ns */
}

66
test/libc/tinymath/log10_test.c Normal file
View file

@ -0,0 +1,66 @@
/*-*- mode:c;indent-tabs-mode:nil;c-basic-offset:2;tab-width:8;coding:utf-8 -*-│
│vi: set net ft=c ts=2 sts=2 sw=2 fenc=utf-8 :vi│
╞══════════════════════════════════════════════════════════════════════════════╡
│ Copyright 2021 Justine Alexandra Roberts Tunney │
│ │
│ Permission to use, copy, modify, and/or distribute this software for │
│ any purpose with or without fee is hereby granted, provided that the │
│ above copyright notice and this permission notice appear in all copies. │
│ │
│ THE SOFTWARE IS PROVIDED "AS IS" AND THE AUTHOR DISCLAIMS ALL │
│ WARRANTIES WITH REGARD TO THIS SOFTWARE INCLUDING ALL IMPLIED │
│ WARRANTIES OF MERCHANTABILITY AND FITNESS. IN NO EVENT SHALL THE │
│ AUTHOR BE LIABLE FOR ANY SPECIAL, DIRECT, INDIRECT, OR CONSEQUENTIAL │
│ DAMAGES OR ANY DAMAGES WHATSOEVER RESULTING FROM LOSS OF USE, DATA OR │
│ PROFITS, WHETHER IN AN ACTION OF CONTRACT, NEGLIGENCE OR OTHER │
│ TORTIOUS ACTION, ARISING OUT OF OR IN CONNECTION WITH THE USE OR │
│ PERFORMANCE OF THIS SOFTWARE. │
╚─────────────────────────────────────────────────────────────────────────────*/
#include "libc/math.h"
#include "libc/runtime/gc.internal.h"
#include "libc/testlib/ezbench.h"
#include "libc/testlib/testlib.h"
#include "libc/x/x.h"
TEST(log10l, test) {
EXPECT_STREQ("1", gc(xdtoal(log10l(10))));
EXPECT_STREQ("NAN", gc(xdtoal(log10l(NAN))));
EXPECT_STREQ("0", gc(xdtoal(log10l(1))));
EXPECT_STREQ("INFINITY", gc(xdtoal(log10l(INFINITY))));
EXPECT_STREQ("-INFINITY", gc(xdtoal(log10l(0))));
EXPECT_STREQ("-INFINITY", gc(xdtoal(log10l(-0.))));
EXPECT_STREQ("-NAN", gc(xdtoal(log10l(-1))));
EXPECT_STREQ("-NAN", gc(xdtoal(log10l(-2))));
}
TEST(log10, test) {
EXPECT_STREQ(".301029995663981", gc(xdtoa(log10(2))));
EXPECT_STREQ("2", gc(xdtoa(log10(100))));
EXPECT_STREQ("NAN", gc(xdtoa(log10(NAN))));
EXPECT_STREQ("0", gc(xdtoa(log10(1))));
EXPECT_STREQ("INFINITY", gc(xdtoa(log10(INFINITY))));
EXPECT_STREQ("-INFINITY", gc(xdtoa(log10(0))));
EXPECT_STREQ("-INFINITY", gc(xdtoa(log10(-0.))));
EXPECT_STREQ("-NAN", gc(xdtoa(log10(-1))));
EXPECT_STREQ("-NAN", gc(xdtoa(log10(-2))));
}
TEST(log10f, test) {
EXPECT_STREQ("1", gc(xdtoaf(log10f(10))));
EXPECT_STREQ("NAN", gc(xdtoaf(log10f(NAN))));
EXPECT_STREQ("0", gc(xdtoaf(log10f(1))));
EXPECT_STREQ("INFINITY", gc(xdtoaf(log10f(INFINITY))));
EXPECT_STREQ("-INFINITY", gc(xdtoaf(log10f(0))));
EXPECT_STREQ("-INFINITY", gc(xdtoaf(log10f(-0.))));
EXPECT_STREQ("-NAN", gc(xdtoaf(log10f(-1))));
EXPECT_STREQ("-NAN", gc(xdtoaf(log10f(-2))));
}
BENCH(log10, bench) {
double _log10(double) asm("log10");
float _log10f(float) asm("log10f");
long double _log10l(long double) asm("log10l");
EZBENCH2("log10", donothing, _log10(.7)); /* ~20ns */
EZBENCH2("log10f", donothing, _log10f(.7)); /* ~21ns */
EZBENCH2("log10l", donothing, _log10l(.7)); /* ~21ns */
}

10
test/libc/tinymath/log2_test.c
View file

@ -18,6 +18,7 @@
╚─────────────────────────────────────────────────────────────────────────────*/
#include "libc/math.h"
#include "libc/runtime/gc.internal.h"
#include "libc/testlib/ezbench.h"
#include "libc/testlib/testlib.h"
#include "libc/x/x.h"
@ -54,3 +55,12 @@ TEST(log2f, test) {
EXPECT_STREQ("-NAN", gc(xdtoaf(log2f(-1))));
EXPECT_STREQ("-NAN", gc(xdtoaf(log2f(-2))));
}
BENCH(log2, bench) {
double _log2(double) asm("log2");
float _log2f(float) asm("log2f");
long double _log2l(long double) asm("log2l");
EZBENCH2("log2", donothing, _log2(.7)); /* ~9ns */
EZBENCH2("log2f", donothing, _log2f(.7)); /* ~6ns */
EZBENCH2("log2l", donothing, _log2l(.7)); /* ~21ns */
}

10
test/libc/tinymath/log_test.c
View file

@ -18,6 +18,7 @@
╚─────────────────────────────────────────────────────────────────────────────*/
#include "libc/math.h"
#include "libc/runtime/gc.internal.h"
#include "libc/testlib/ezbench.h"
#include "libc/testlib/testlib.h"
#include "libc/x/x.h"
@ -54,3 +55,12 @@ TEST(logf, test) {
EXPECT_STREQ("-NAN", gc(xdtoaf(logf(-1))));
EXPECT_STREQ("-NAN", gc(xdtoaf(logf(-2))));
}
BENCH(logl, bench) {
double _log(double) asm("log");
float _logf(float) asm("logf");
long double _logl(long double) asm("logl");
EZBENCH2("log", donothing, _log(.7)); /* ~9ns */
EZBENCH2("logf", donothing, _logf(.7)); /* ~20ns */
EZBENCH2("logl", donothing, _logl(.7)); /* ~20ns */
}

536
test/libc/tinymath/powl_test.c
View file

@ -40,6 +40,10 @@ void SetUp(void) {
rando = rand() & 0xffff;
}
static double POW(double x, double y) {
return powl(x, y);
}
static char *fmtd(double x) {
sprintf(buf, "%.15g", x);
return buf;
@ -209,356 +213,356 @@ TEST(powf, test) {
}
TEST(powl, errors) {
EXPECT_STREQ("1", fmtd(pow(0., 0.)));
EXPECT_STREQ("1", fmtd(pow(0., -0.)));
EXPECT_STREQ("0", fmtd(pow(0., .5)));
EXPECT_STREQ("1", fmtd(POW(0., 0.)));
EXPECT_STREQ("1", fmtd(POW(0., -0.)));
EXPECT_STREQ("0", fmtd(POW(0., .5)));
errno = 0;
EXPECT_STREQ("inf", fmtd(pow(0., -.5)));
EXPECT_STREQ("inf", fmtd(POW(0., -.5)));
EXPECT_EQ(ERANGE, errno);
EXPECT_STREQ("0", fmtd(pow(0., 1.)));
EXPECT_STREQ("inf", fmtd(pow(0., -1.)));
EXPECT_STREQ("0", fmtd(pow(0., 1.5)));
EXPECT_STREQ("0", fmtd(POW(0., 1.)));
EXPECT_STREQ("inf", fmtd(POW(0., -1.)));
EXPECT_STREQ("0", fmtd(POW(0., 1.5)));
errno = 0;
EXPECT_STREQ("inf", fmtd(pow(0., -1.5)));
EXPECT_STREQ("inf", fmtd(POW(0., -1.5)));
EXPECT_EQ(ERANGE, errno);
EXPECT_STREQ("0", fmtd(pow(0., 2.)));
EXPECT_STREQ("0", fmtd(POW(0., 2.)));
errno = 0;
EXPECT_STREQ("inf", fmtd(pow(0., -2.)));
EXPECT_STREQ("inf", fmtd(POW(0., -2.)));
EXPECT_EQ(ERANGE, errno);
EXPECT_TRUE(isnan(pow(0., NAN)));
EXPECT_TRUE(isnan(pow(0., -NAN)));
EXPECT_STREQ("0", fmtd(pow(0., INFINITY)));
EXPECT_STREQ("inf", fmtd(pow(0., -INFINITY)));
EXPECT_STREQ("0", fmtd(pow(0., __DBL_MIN__)));
EXPECT_STREQ("0", fmtd(pow(0., __DBL_MAX__)));
EXPECT_STREQ("1", fmtd(pow(-0., 0.)));
EXPECT_STREQ("1", fmtd(pow(-0., -0.)));
EXPECT_STREQ("0", fmtd(pow(-0., .5)));
EXPECT_TRUE(isnan(POW(0., NAN)));
EXPECT_TRUE(isnan(POW(0., -NAN)));
EXPECT_STREQ("0", fmtd(POW(0., INFINITY)));
EXPECT_STREQ("inf", fmtd(POW(0., -INFINITY)));
EXPECT_STREQ("0", fmtd(POW(0., __DBL_MIN__)));
EXPECT_STREQ("0", fmtd(POW(0., __DBL_MAX__)));
EXPECT_STREQ("1", fmtd(POW(-0., 0.)));
EXPECT_STREQ("1", fmtd(POW(-0., -0.)));
EXPECT_STREQ("0", fmtd(POW(-0., .5)));
errno = 0;
EXPECT_STREQ("inf", fmtd(pow(-0., -.5)));
EXPECT_STREQ("inf", fmtd(POW(-0., -.5)));
EXPECT_EQ(ERANGE, errno);
EXPECT_STREQ("-0", fmtd(pow(-0., 1.)));
EXPECT_STREQ("-inf", fmtd(pow(-0., -1.)));
EXPECT_STREQ("0", fmtd(pow(-0., 1.5)));
EXPECT_STREQ("-0", fmtd(POW(-0., 1.)));
EXPECT_STREQ("-inf", fmtd(POW(-0., -1.)));
EXPECT_STREQ("0", fmtd(POW(-0., 1.5)));
errno = 0;
EXPECT_STREQ("inf", fmtd(pow(-0., -1.5)));
EXPECT_STREQ("inf", fmtd(POW(-0., -1.5)));
EXPECT_EQ(ERANGE, errno);
EXPECT_STREQ("0", fmtd(pow(-0., 2.)));
EXPECT_STREQ("0", fmtd(POW(-0., 2.)));
errno = 0;
EXPECT_STREQ("inf", fmtd(pow(-0., -2.)));
EXPECT_STREQ("inf", fmtd(POW(-0., -2.)));
EXPECT_EQ(ERANGE, errno);
EXPECT_TRUE(isnan(pow(-0., NAN)));
EXPECT_TRUE(isnan(pow(-0., -NAN)));
EXPECT_STREQ("0", fmtd(pow(-0., INFINITY)));
EXPECT_STREQ("inf", fmtd(pow(-0., -INFINITY)));
EXPECT_STREQ("0", fmtd(pow(-0., __DBL_MIN__)));
EXPECT_STREQ("0", fmtd(pow(-0., __DBL_MAX__)));
EXPECT_STREQ("1", fmtd(pow(.5, 0.)));
EXPECT_STREQ("1", fmtd(pow(.5, -0.)));
EXPECT_STREQ("0.707106781186548", fmtd(pow(.5, .5)));
EXPECT_STREQ("1.4142135623731", fmtd(pow(.5, -.5)));
EXPECT_STREQ("0.5", fmtd(pow(.5, 1.)));
EXPECT_STREQ("2", fmtd(pow(.5, -1.)));
EXPECT_STREQ("0.353553390593274", fmtd(pow(.5, 1.5)));
EXPECT_STREQ("2.82842712474619", fmtd(pow(.5, -1.5)));
EXPECT_STREQ("0.25", fmtd(pow(.5, 2.)));
EXPECT_STREQ("4", fmtd(pow(.5, -2.)));
EXPECT_TRUE(isnan(pow(.5, NAN)));
EXPECT_TRUE(isnan(pow(.5, -NAN)));
EXPECT_STREQ("0", fmtd(pow(.5, INFINITY)));
EXPECT_STREQ("inf", fmtd(pow(.5, -INFINITY)));
EXPECT_STREQ("1", fmtd(pow(.5, __DBL_MIN__)));
EXPECT_TRUE(isnan(POW(-0., NAN)));
EXPECT_TRUE(isnan(POW(-0., -NAN)));
EXPECT_STREQ("0", fmtd(POW(-0., INFINITY)));
EXPECT_STREQ("inf", fmtd(POW(-0., -INFINITY)));
EXPECT_STREQ("0", fmtd(POW(-0., __DBL_MIN__)));
EXPECT_STREQ("0", fmtd(POW(-0., __DBL_MAX__)));
EXPECT_STREQ("1", fmtd(POW(.5, 0.)));
EXPECT_STREQ("1", fmtd(POW(.5, -0.)));
EXPECT_STREQ("0.707106781186548", fmtd(POW(.5, .5)));
EXPECT_STREQ("1.4142135623731", fmtd(POW(.5, -.5)));
EXPECT_STREQ("0.5", fmtd(POW(.5, 1.)));
EXPECT_STREQ("2", fmtd(POW(.5, -1.)));
EXPECT_STREQ("0.353553390593274", fmtd(POW(.5, 1.5)));
EXPECT_STREQ("2.82842712474619", fmtd(POW(.5, -1.5)));
EXPECT_STREQ("0.25", fmtd(POW(.5, 2.)));
EXPECT_STREQ("4", fmtd(POW(.5, -2.)));
EXPECT_TRUE(isnan(POW(.5, NAN)));
EXPECT_TRUE(isnan(POW(.5, -NAN)));
EXPECT_STREQ("0", fmtd(POW(.5, INFINITY)));
EXPECT_STREQ("inf", fmtd(POW(.5, -INFINITY)));
EXPECT_STREQ("1", fmtd(POW(.5, __DBL_MIN__)));
errno = 0;
EXPECT_STREQ("0", fmtd(pow(.5, __DBL_MAX__)));
EXPECT_STREQ("0", fmtd(POW(.5, __DBL_MAX__)));
/* EXPECT_EQ(ERANGE, errno); */
EXPECT_STREQ("1", fmtd(pow(-.5, 0.)));
EXPECT_STREQ("1", fmtd(pow(-.5, -0.)));
EXPECT_STREQ("1", fmtd(POW(-.5, 0.)));
EXPECT_STREQ("1", fmtd(POW(-.5, -0.)));
errno = 0;
EXPECT_TRUE(isnan(pow(-.5, .5)));
EXPECT_TRUE(isnan(POW(-.5, .5)));
EXPECT_EQ(EDOM, errno);
errno = 0;
EXPECT_TRUE(isnan(pow(-.5, -.5)));
EXPECT_TRUE(isnan(POW(-.5, -.5)));
EXPECT_EQ(EDOM, errno);
EXPECT_STREQ("-0.5", fmtd(pow(-.5, 1.)));
EXPECT_STREQ("-2", fmtd(pow(-.5, -1.)));
EXPECT_STREQ("-0.5", fmtd(POW(-.5, 1.)));
EXPECT_STREQ("-2", fmtd(POW(-.5, -1.)));
errno = 0;
EXPECT_TRUE(isnan(pow(-.5, 1.5)));
EXPECT_TRUE(isnan(POW(-.5, 1.5)));
EXPECT_EQ(EDOM, errno);
errno = 0;
EXPECT_TRUE(isnan(pow(-.5, -1.5)));
EXPECT_TRUE(isnan(POW(-.5, -1.5)));
EXPECT_EQ(EDOM, errno);
EXPECT_STREQ("0.25", fmtd(pow(-.5, 2.)));
EXPECT_STREQ("4", fmtd(pow(-.5, -2.)));
EXPECT_TRUE(isnan(pow(-.5, NAN)));
EXPECT_TRUE(isnan(pow(-.5, -NAN)));
EXPECT_STREQ("0", fmtd(pow(-.5, INFINITY)));
EXPECT_STREQ("inf", fmtd(pow(-.5, -INFINITY)));
EXPECT_STREQ("0.25", fmtd(POW(-.5, 2.)));
EXPECT_STREQ("4", fmtd(POW(-.5, -2.)));
EXPECT_TRUE(isnan(POW(-.5, NAN)));
EXPECT_TRUE(isnan(POW(-.5, -NAN)));
EXPECT_STREQ("0", fmtd(POW(-.5, INFINITY)));
EXPECT_STREQ("inf", fmtd(POW(-.5, -INFINITY)));
errno = 0;
EXPECT_TRUE(isnan(pow(-.5, __DBL_MIN__)));
EXPECT_TRUE(isnan(POW(-.5, __DBL_MIN__)));
EXPECT_EQ(EDOM, errno);
errno = 0;
EXPECT_STREQ("0", fmtd(pow(-.5, __DBL_MAX__)));
EXPECT_STREQ("0", fmtd(POW(-.5, __DBL_MAX__)));
/* EXPECT_EQ(ERANGE, errno); */
EXPECT_STREQ("1", fmtd(pow(1., 0.)));
EXPECT_STREQ("1", fmtd(pow(1., -0.)));
EXPECT_STREQ("1", fmtd(pow(1., .5)));
EXPECT_STREQ("1", fmtd(pow(1., -.5)));
EXPECT_STREQ("1", fmtd(pow(1., 1.)));
EXPECT_STREQ("1", fmtd(pow(1., -1.)));
EXPECT_STREQ("1", fmtd(pow(1., 1.5)));
EXPECT_STREQ("1", fmtd(pow(1., -1.5)));
EXPECT_STREQ("1", fmtd(pow(1., 2.)));
EXPECT_STREQ("1", fmtd(pow(1., -2.)));
EXPECT_STREQ("1", fmtd(pow(1., NAN)));
EXPECT_STREQ("1", fmtd(pow(1., -NAN)));
EXPECT_STREQ("1", fmtd(pow(1., INFINITY)));
EXPECT_STREQ("1", fmtd(pow(1., -INFINITY)));
EXPECT_STREQ("1", fmtd(pow(1., __DBL_MIN__)));
EXPECT_STREQ("1", fmtd(pow(1., __DBL_MAX__)));
EXPECT_STREQ("1", fmtd(pow(-1., 0.)));
EXPECT_STREQ("1", fmtd(pow(-1., -0.)));
EXPECT_STREQ("1", fmtd(POW(1., 0.)));
EXPECT_STREQ("1", fmtd(POW(1., -0.)));
EXPECT_STREQ("1", fmtd(POW(1., .5)));
EXPECT_STREQ("1", fmtd(POW(1., -.5)));
EXPECT_STREQ("1", fmtd(POW(1., 1.)));
EXPECT_STREQ("1", fmtd(POW(1., -1.)));
EXPECT_STREQ("1", fmtd(POW(1., 1.5)));
EXPECT_STREQ("1", fmtd(POW(1., -1.5)));
EXPECT_STREQ("1", fmtd(POW(1., 2.)));
EXPECT_STREQ("1", fmtd(POW(1., -2.)));
EXPECT_STREQ("1", fmtd(POW(1., NAN)));
EXPECT_STREQ("1", fmtd(POW(1., -NAN)));
EXPECT_STREQ("1", fmtd(POW(1., INFINITY)));
EXPECT_STREQ("1", fmtd(POW(1., -INFINITY)));
EXPECT_STREQ("1", fmtd(POW(1., __DBL_MIN__)));
EXPECT_STREQ("1", fmtd(POW(1., __DBL_MAX__)));
EXPECT_STREQ("1", fmtd(POW(-1., 0.)));
EXPECT_STREQ("1", fmtd(POW(-1., -0.)));
errno = 0;
EXPECT_TRUE(isnan(pow(-1., .5)));
EXPECT_TRUE(isnan(POW(-1., .5)));
EXPECT_EQ(EDOM, errno);
errno = 0;
EXPECT_TRUE(isnan(pow(-1., -.5)));
EXPECT_TRUE(isnan(POW(-1., -.5)));
EXPECT_EQ(EDOM, errno);
EXPECT_STREQ("-1", fmtd(pow(-1., 1.)));
EXPECT_STREQ("-1", fmtd(pow(-1., -1.)));
EXPECT_STREQ("-1", fmtd(POW(-1., 1.)));
EXPECT_STREQ("-1", fmtd(POW(-1., -1.)));
errno = 0;
EXPECT_TRUE(isnan(pow(-1., 1.5)));
EXPECT_TRUE(isnan(POW(-1., 1.5)));
EXPECT_EQ(EDOM, errno);
errno = 0;
EXPECT_TRUE(isnan(pow(-1., -1.5)));
EXPECT_TRUE(isnan(POW(-1., -1.5)));
EXPECT_EQ(EDOM, errno);
EXPECT_STREQ("1", fmtd(pow(-1., 2.0)));
EXPECT_STREQ("1", fmtd(pow(-1., -2.0)));
EXPECT_TRUE(isnan(pow(-1., NAN)));
EXPECT_TRUE(isnan(pow(-1., -NAN)));
EXPECT_STREQ("1", fmtd(pow(-1., INFINITY)));
EXPECT_STREQ("1", fmtd(pow(-1., -INFINITY)));
EXPECT_STREQ("1", fmtd(POW(-1., 2.0)));
EXPECT_STREQ("1", fmtd(POW(-1., -2.0)));
EXPECT_TRUE(isnan(POW(-1., NAN)));
EXPECT_TRUE(isnan(POW(-1., -NAN)));
EXPECT_STREQ("1", fmtd(POW(-1., INFINITY)));
EXPECT_STREQ("1", fmtd(POW(-1., -INFINITY)));
errno = 0;
EXPECT_TRUE(isnan(pow(-1., __DBL_MIN__)));
EXPECT_TRUE(isnan(POW(-1., __DBL_MIN__)));
EXPECT_EQ(EDOM, errno);
EXPECT_STREQ("1", fmtd(pow(-1., __DBL_MAX__)));
EXPECT_STREQ("1", fmtd(pow(1.5, 0.)));
EXPECT_STREQ("1", fmtd(pow(1.5, -0.)));
EXPECT_STREQ("1.22474487139159", fmtd(pow(1.5, .5)));
EXPECT_STREQ("0.816496580927726", fmtd(pow(1.5, -.5)));
EXPECT_STREQ("1.5", fmtd(pow(1.5, 1.)));
EXPECT_STREQ("0.666666666666667", fmtd(pow(1.5, -1.)));
EXPECT_STREQ("1.83711730708738", fmtd(pow(1.5, 1.5)));
EXPECT_STREQ("0.544331053951817", fmtd(pow(1.5, -1.5)));
EXPECT_STREQ("2.25", fmtd(pow(1.5, 2.0)));
EXPECT_STREQ("0.444444444444444", fmtd(pow(1.5, -2.0)));
EXPECT_TRUE(isnan(pow(1.5, NAN)));
EXPECT_TRUE(isnan(pow(1.5, -NAN)));
EXPECT_STREQ("inf", fmtd(pow(1.5, INFINITY)));
EXPECT_STREQ("0", fmtd(pow(1.5, -INFINITY)));
EXPECT_STREQ("1", fmtd(pow(1.5, __DBL_MIN__)));
EXPECT_STREQ("1", fmtd(POW(-1., __DBL_MAX__)));
EXPECT_STREQ("1", fmtd(POW(1.5, 0.)));
EXPECT_STREQ("1", fmtd(POW(1.5, -0.)));
EXPECT_STREQ("1.22474487139159", fmtd(POW(1.5, .5)));
EXPECT_STREQ("0.816496580927726", fmtd(POW(1.5, -.5)));
EXPECT_STREQ("1.5", fmtd(POW(1.5, 1.)));
EXPECT_STREQ("0.666666666666667", fmtd(POW(1.5, -1.)));
EXPECT_STREQ("1.83711730708738", fmtd(POW(1.5, 1.5)));
EXPECT_STREQ("0.544331053951817", fmtd(POW(1.5, -1.5)));
EXPECT_STREQ("2.25", fmtd(POW(1.5, 2.0)));
EXPECT_STREQ("0.444444444444444", fmtd(POW(1.5, -2.0)));
EXPECT_TRUE(isnan(POW(1.5, NAN)));
EXPECT_TRUE(isnan(POW(1.5, -NAN)));
EXPECT_STREQ("inf", fmtd(POW(1.5, INFINITY)));
EXPECT_STREQ("0", fmtd(POW(1.5, -INFINITY)));
EXPECT_STREQ("1", fmtd(POW(1.5, __DBL_MIN__)));
errno = 0;
EXPECT_STREQ("inf", fmtd(pow(1.5, __DBL_MAX__)));
EXPECT_STREQ("inf", fmtd(POW(1.5, __DBL_MAX__)));
/* EXPECT_EQ(ERANGE, errno); */
EXPECT_STREQ("1", fmtd(pow(-1.5, 0.)));
EXPECT_STREQ("1", fmtd(pow(-1.5, -0.)));
EXPECT_STREQ("1", fmtd(POW(-1.5, 0.)));
EXPECT_STREQ("1", fmtd(POW(-1.5, -0.)));
errno = 0;
EXPECT_TRUE(isnan(pow(-1.5, .5)));
EXPECT_TRUE(isnan(POW(-1.5, .5)));
EXPECT_EQ(EDOM, errno);
errno = 0;
EXPECT_TRUE(isnan(pow(-1.5, -.5)));
EXPECT_TRUE(isnan(POW(-1.5, -.5)));
EXPECT_EQ(EDOM, errno);
EXPECT_STREQ("-1.5", fmtd(pow(-1.5, 1.)));
EXPECT_STREQ("-0.666666666666667", fmtd(pow(-1.5, -1.)));
EXPECT_STREQ("-1.5", fmtd(POW(-1.5, 1.)));
EXPECT_STREQ("-0.666666666666667", fmtd(POW(-1.5, -1.)));
errno = 0;
EXPECT_TRUE(isnan(pow(-1.5, 1.5)));
EXPECT_TRUE(isnan(POW(-1.5, 1.5)));
EXPECT_EQ(EDOM, errno);
errno = 0;
EXPECT_TRUE(isnan(pow(-1.5, -1.5)));
EXPECT_TRUE(isnan(POW(-1.5, -1.5)));
EXPECT_EQ(EDOM, errno);
EXPECT_STREQ("2.25", fmtd(pow(-1.5, 2.0)));
EXPECT_STREQ("0.444444444444444", fmtd(pow(-1.5, -2.0)));
EXPECT_TRUE(isnan(pow(-1.5, NAN)));
EXPECT_TRUE(isnan(pow(-1.5, -NAN)));
EXPECT_STREQ("inf", fmtd(pow(-1.5, INFINITY)));
EXPECT_STREQ("0", fmtd(pow(-1.5, -INFINITY)));
EXPECT_STREQ("2.25", fmtd(POW(-1.5, 2.0)));
EXPECT_STREQ("0.444444444444444", fmtd(POW(-1.5, -2.0)));
EXPECT_TRUE(isnan(POW(-1.5, NAN)));
EXPECT_TRUE(isnan(POW(-1.5, -NAN)));
EXPECT_STREQ("inf", fmtd(POW(-1.5, INFINITY)));
EXPECT_STREQ("0", fmtd(POW(-1.5, -INFINITY)));
errno = 0;
EXPECT_TRUE(isnan(pow(-1.5, __DBL_MIN__)));
EXPECT_TRUE(isnan(POW(-1.5, __DBL_MIN__)));
EXPECT_EQ(EDOM, errno);
errno = 0;
EXPECT_STREQ("inf", fmtd(pow(-1.5, __DBL_MAX__)));
EXPECT_STREQ("inf", fmtd(POW(-1.5, __DBL_MAX__)));
/* EXPECT_EQ(ERANGE, errno); */
EXPECT_STREQ("1", fmtd(pow(+2.0, 0.)));
EXPECT_STREQ("1", fmtd(pow(+2.0, -0.)));
EXPECT_STREQ("1.4142135623731", fmtd(pow(+2.0, .5)));
EXPECT_STREQ("0.707106781186548", fmtd(pow(+2.0, -.5)));
EXPECT_STREQ("2", fmtd(pow(+2.0, 1.)));
EXPECT_STREQ("0.5", fmtd(pow(+2.0, -1.)));
EXPECT_STREQ("2.82842712474619", fmtd(pow(+2.0, 1.5)));
EXPECT_STREQ("0.353553390593274", fmtd(pow(+2.0, -1.5)));
EXPECT_STREQ("4", fmtd(pow(+2.0, 2.0)));
EXPECT_STREQ("0.25", fmtd(pow(+2.0, -2.0)));
EXPECT_TRUE(isnan(pow(+2.0, NAN)));
EXPECT_TRUE(isnan(pow(+2.0, -NAN)));
EXPECT_STREQ("inf", fmtd(pow(+2.0, INFINITY)));
EXPECT_STREQ("0", fmtd(pow(+2.0, -INFINITY)));
EXPECT_STREQ("1", fmtd(pow(+2.0, __DBL_MIN__)));
EXPECT_STREQ("1", fmtd(POW(+2.0, 0.)));
EXPECT_STREQ("1", fmtd(POW(+2.0, -0.)));
EXPECT_STREQ("1.4142135623731", fmtd(POW(+2.0, .5)));
EXPECT_STREQ("0.707106781186548", fmtd(POW(+2.0, -.5)));
EXPECT_STREQ("2", fmtd(POW(+2.0, 1.)));
EXPECT_STREQ("0.5", fmtd(POW(+2.0, -1.)));
EXPECT_STREQ("2.82842712474619", fmtd(POW(+2.0, 1.5)));
EXPECT_STREQ("0.353553390593274", fmtd(POW(+2.0, -1.5)));
EXPECT_STREQ("4", fmtd(POW(+2.0, 2.0)));
EXPECT_STREQ("0.25", fmtd(POW(+2.0, -2.0)));
EXPECT_TRUE(isnan(POW(+2.0, NAN)));
EXPECT_TRUE(isnan(POW(+2.0, -NAN)));
EXPECT_STREQ("inf", fmtd(POW(+2.0, INFINITY)));
EXPECT_STREQ("0", fmtd(POW(+2.0, -INFINITY)));
EXPECT_STREQ("1", fmtd(POW(+2.0, __DBL_MIN__)));
errno = 0;
EXPECT_STREQ("inf", fmtd(pow(+2.0, __DBL_MAX__)));
EXPECT_STREQ("inf", fmtd(POW(+2.0, __DBL_MAX__)));
/* EXPECT_EQ(ERANGE, errno); */
EXPECT_STREQ("1", fmtd(pow(-2.0, 0.)));
EXPECT_STREQ("1", fmtd(pow(-2.0, -0.)));
EXPECT_STREQ("1", fmtd(POW(-2.0, 0.)));
EXPECT_STREQ("1", fmtd(POW(-2.0, -0.)));
errno = 0;
EXPECT_TRUE(isnan(pow(-2.0, .5)));
EXPECT_TRUE(isnan(POW(-2.0, .5)));
EXPECT_EQ(EDOM, errno);
errno = 0;
EXPECT_TRUE(isnan(pow(-2.0, -.5)));
EXPECT_TRUE(isnan(POW(-2.0, -.5)));
EXPECT_EQ(EDOM, errno);
EXPECT_STREQ("-2", fmtd(pow(-2.0, 1.)));
EXPECT_STREQ("-0.5", fmtd(pow(-2.0, -1.)));
EXPECT_STREQ("-2", fmtd(POW(-2.0, 1.)));
EXPECT_STREQ("-0.5", fmtd(POW(-2.0, -1.)));
errno = 0;
EXPECT_TRUE(isnan(pow(-2.0, 1.5)));
EXPECT_TRUE(isnan(POW(-2.0, 1.5)));
EXPECT_EQ(EDOM, errno);
errno = 0;
EXPECT_TRUE(isnan(pow(-2.0, -1.5)));
EXPECT_TRUE(isnan(POW(-2.0, -1.5)));
EXPECT_EQ(EDOM, errno);
EXPECT_STREQ("4", fmtd(pow(-2.0, 2.0)));
EXPECT_STREQ("0.25", fmtd(pow(-2.0, -2.0)));
EXPECT_TRUE(isnan(pow(-2.0, NAN)));
EXPECT_TRUE(isnan(pow(-2.0, -NAN)));
EXPECT_STREQ("inf", fmtd(pow(-2.0, INFINITY)));
EXPECT_STREQ("0", fmtd(pow(-2.0, -INFINITY)));
EXPECT_STREQ("4", fmtd(POW(-2.0, 2.0)));
EXPECT_STREQ("0.25", fmtd(POW(-2.0, -2.0)));
EXPECT_TRUE(isnan(POW(-2.0, NAN)));
EXPECT_TRUE(isnan(POW(-2.0, -NAN)));
EXPECT_STREQ("inf", fmtd(POW(-2.0, INFINITY)));
EXPECT_STREQ("0", fmtd(POW(-2.0, -INFINITY)));
errno = 0;
EXPECT_TRUE(isnan(pow(-2.0, __DBL_MIN__)));
EXPECT_TRUE(isnan(POW(-2.0, __DBL_MIN__)));
EXPECT_EQ(EDOM, errno);
errno = 0;
EXPECT_STREQ("inf", fmtd(pow(-2.0, __DBL_MAX__)));
EXPECT_STREQ("inf", fmtd(POW(-2.0, __DBL_MAX__)));
/* EXPECT_EQ(ERANGE, errno); */
EXPECT_STREQ("1", fmtd(pow(NAN, 0.)));
EXPECT_STREQ("1", fmtd(pow(NAN, -0.)));
EXPECT_TRUE(isnan(pow(NAN, .5)));
EXPECT_TRUE(isnan(pow(NAN, -.5)));
EXPECT_TRUE(isnan(pow(NAN, 1.)));
EXPECT_TRUE(isnan(pow(NAN, -1.)));
EXPECT_TRUE(isnan(pow(NAN, 1.5)));
EXPECT_TRUE(isnan(pow(NAN, -1.5)));
EXPECT_TRUE(isnan(pow(NAN, 2.0)));
EXPECT_TRUE(isnan(pow(NAN, -2.0)));
EXPECT_TRUE(isnan(pow(NAN, NAN)));
EXPECT_TRUE(isnan(pow(NAN, -NAN)));
EXPECT_TRUE(isnan(pow(NAN, INFINITY)));
EXPECT_TRUE(isnan(pow(NAN, -INFINITY)));
EXPECT_TRUE(isnan(pow(NAN, __DBL_MIN__)));
EXPECT_TRUE(isnan(pow(NAN, __DBL_MAX__)));
EXPECT_STREQ("1", fmtd(pow(-NAN, 0.)));
EXPECT_STREQ("1", fmtd(pow(-NAN, -0.)));
EXPECT_TRUE(isnan(pow(-NAN, .5)));
EXPECT_TRUE(isnan(pow(-NAN, -.5)));
EXPECT_TRUE(isnan(pow(-NAN, 1.)));
EXPECT_TRUE(isnan(pow(-NAN, -1.)));
EXPECT_TRUE(isnan(pow(-NAN, 1.5)));
EXPECT_TRUE(isnan(pow(-NAN, -1.5)));
EXPECT_TRUE(isnan(pow(-NAN, 2.0)));
EXPECT_TRUE(isnan(pow(-NAN, -2.0)));
EXPECT_TRUE(isnan(pow(-NAN, NAN)));
EXPECT_TRUE(isnan(pow(-NAN, -NAN)));
EXPECT_TRUE(isnan(pow(-NAN, INFINITY)));
EXPECT_TRUE(isnan(pow(-NAN, -INFINITY)));
EXPECT_TRUE(isnan(pow(-NAN, __DBL_MIN__)));
EXPECT_TRUE(isnan(pow(-NAN, __DBL_MAX__)));
EXPECT_STREQ("1", fmtd(pow(INFINITY, 0.)));
EXPECT_STREQ("1", fmtd(pow(INFINITY, -0.)));
EXPECT_STREQ("inf", fmtd(pow(INFINITY, .5)));
EXPECT_STREQ("0", fmtd(pow(INFINITY, -.5)));
EXPECT_STREQ("inf", fmtd(pow(INFINITY, 1.)));
EXPECT_STREQ("0", fmtd(pow(INFINITY, -1.)));
EXPECT_STREQ("inf", fmtd(pow(INFINITY, 1.5)));
EXPECT_STREQ("0", fmtd(pow(INFINITY, -1.5)));
EXPECT_STREQ("inf", fmtd(pow(INFINITY, 2.0)));
EXPECT_STREQ("0", fmtd(pow(INFINITY, -2.0)));
EXPECT_TRUE(isnan(pow(INFINITY, NAN)));
EXPECT_TRUE(isnan(pow(INFINITY, -NAN)));
EXPECT_STREQ("inf", fmtd(pow(INFINITY, INFINITY)));
EXPECT_STREQ("0", fmtd(pow(INFINITY, -INFINITY)));
EXPECT_STREQ("inf", fmtd(pow(INFINITY, __DBL_MIN__)));
EXPECT_STREQ("inf", fmtd(pow(INFINITY, __DBL_MAX__)));
EXPECT_STREQ("1", fmtd(pow(-INFINITY, 0.)));
EXPECT_STREQ("1", fmtd(pow(-INFINITY, -0.)));
EXPECT_STREQ("inf", fmtd(pow(-INFINITY, .5)));
EXPECT_STREQ("0", fmtd(pow(-INFINITY, -.5)));
EXPECT_STREQ("-inf", fmtd(pow(-INFINITY, 1.)));
EXPECT_STREQ("-0", fmtd(pow(-INFINITY, -1.)));
EXPECT_STREQ("inf", fmtd(pow(-INFINITY, 1.5)));
EXPECT_STREQ("0", fmtd(pow(-INFINITY, -1.5)));
EXPECT_STREQ("inf", fmtd(pow(-INFINITY, 2.0)));
EXPECT_STREQ("0", fmtd(pow(-INFINITY, -2.0)));
EXPECT_TRUE(isnan(pow(-INFINITY, NAN)));
EXPECT_TRUE(isnan(pow(-INFINITY, -NAN)));
EXPECT_STREQ("inf", fmtd(pow(-INFINITY, INFINITY)));
EXPECT_STREQ("0", fmtd(pow(-INFINITY, -INFINITY)));
EXPECT_STREQ("inf", fmtd(pow(-INFINITY, __DBL_MIN__)));
EXPECT_STREQ("inf", fmtd(pow(-INFINITY, __DBL_MAX__)));
EXPECT_STREQ("1", fmtd(pow(__DBL_MIN__, 0.)));
EXPECT_STREQ("1", fmtd(pow(__DBL_MIN__, -0.)));
EXPECT_STREQ("1.49166814624004e-154", fmtd(pow(__DBL_MIN__, .5)));
EXPECT_STREQ("6.7039039649713e+153", fmtd(pow(__DBL_MIN__, -.5)));
EXPECT_STREQ("2.2250738585072e-308", fmtd(pow(__DBL_MIN__, 1.)));
EXPECT_STREQ("4.49423283715579e+307", fmtd(pow(__DBL_MIN__, -1.)));
EXPECT_STREQ("1", fmtd(POW(NAN, 0.)));
EXPECT_STREQ("1", fmtd(POW(NAN, -0.)));
EXPECT_TRUE(isnan(POW(NAN, .5)));
EXPECT_TRUE(isnan(POW(NAN, -.5)));
EXPECT_TRUE(isnan(POW(NAN, 1.)));
EXPECT_TRUE(isnan(POW(NAN, -1.)));
EXPECT_TRUE(isnan(POW(NAN, 1.5)));
EXPECT_TRUE(isnan(POW(NAN, -1.5)));
EXPECT_TRUE(isnan(POW(NAN, 2.0)));
EXPECT_TRUE(isnan(POW(NAN, -2.0)));
EXPECT_TRUE(isnan(POW(NAN, NAN)));
EXPECT_TRUE(isnan(POW(NAN, -NAN)));
EXPECT_TRUE(isnan(POW(NAN, INFINITY)));
EXPECT_TRUE(isnan(POW(NAN, -INFINITY)));
EXPECT_TRUE(isnan(POW(NAN, __DBL_MIN__)));
EXPECT_TRUE(isnan(POW(NAN, __DBL_MAX__)));
EXPECT_STREQ("1", fmtd(POW(-NAN, 0.)));
EXPECT_STREQ("1", fmtd(POW(-NAN, -0.)));
EXPECT_TRUE(isnan(POW(-NAN, .5)));
EXPECT_TRUE(isnan(POW(-NAN, -.5)));
EXPECT_TRUE(isnan(POW(-NAN, 1.)));
EXPECT_TRUE(isnan(POW(-NAN, -1.)));
EXPECT_TRUE(isnan(POW(-NAN, 1.5)));
EXPECT_TRUE(isnan(POW(-NAN, -1.5)));
EXPECT_TRUE(isnan(POW(-NAN, 2.0)));
EXPECT_TRUE(isnan(POW(-NAN, -2.0)));
EXPECT_TRUE(isnan(POW(-NAN, NAN)));
EXPECT_TRUE(isnan(POW(-NAN, -NAN)));
EXPECT_TRUE(isnan(POW(-NAN, INFINITY)));
EXPECT_TRUE(isnan(POW(-NAN, -INFINITY)));
EXPECT_TRUE(isnan(POW(-NAN, __DBL_MIN__)));
EXPECT_TRUE(isnan(POW(-NAN, __DBL_MAX__)));
EXPECT_STREQ("1", fmtd(POW(INFINITY, 0.)));
EXPECT_STREQ("1", fmtd(POW(INFINITY, -0.)));
EXPECT_STREQ("inf", fmtd(POW(INFINITY, .5)));
EXPECT_STREQ("0", fmtd(POW(INFINITY, -.5)));
EXPECT_STREQ("inf", fmtd(POW(INFINITY, 1.)));
EXPECT_STREQ("0", fmtd(POW(INFINITY, -1.)));
EXPECT_STREQ("inf", fmtd(POW(INFINITY, 1.5)));
EXPECT_STREQ("0", fmtd(POW(INFINITY, -1.5)));
EXPECT_STREQ("inf", fmtd(POW(INFINITY, 2.0)));
EXPECT_STREQ("0", fmtd(POW(INFINITY, -2.0)));
EXPECT_TRUE(isnan(POW(INFINITY, NAN)));
EXPECT_TRUE(isnan(POW(INFINITY, -NAN)));
EXPECT_STREQ("inf", fmtd(POW(INFINITY, INFINITY)));
EXPECT_STREQ("0", fmtd(POW(INFINITY, -INFINITY)));
EXPECT_STREQ("inf", fmtd(POW(INFINITY, __DBL_MIN__)));
EXPECT_STREQ("inf", fmtd(POW(INFINITY, __DBL_MAX__)));
EXPECT_STREQ("1", fmtd(POW(-INFINITY, 0.)));
EXPECT_STREQ("1", fmtd(POW(-INFINITY, -0.)));
EXPECT_STREQ("inf", fmtd(POW(-INFINITY, .5)));
EXPECT_STREQ("0", fmtd(POW(-INFINITY, -.5)));
EXPECT_STREQ("-inf", fmtd(POW(-INFINITY, 1.)));
EXPECT_STREQ("-0", fmtd(POW(-INFINITY, -1.)));
EXPECT_STREQ("inf", fmtd(POW(-INFINITY, 1.5)));
EXPECT_STREQ("0", fmtd(POW(-INFINITY, -1.5)));
EXPECT_STREQ("inf", fmtd(POW(-INFINITY, 2.0)));
EXPECT_STREQ("0", fmtd(POW(-INFINITY, -2.0)));
EXPECT_TRUE(isnan(POW(-INFINITY, NAN)));
EXPECT_TRUE(isnan(POW(-INFINITY, -NAN)));
EXPECT_STREQ("inf", fmtd(POW(-INFINITY, INFINITY)));
EXPECT_STREQ("0", fmtd(POW(-INFINITY, -INFINITY)));
EXPECT_STREQ("inf", fmtd(POW(-INFINITY, __DBL_MIN__)));
EXPECT_STREQ("inf", fmtd(POW(-INFINITY, __DBL_MAX__)));
EXPECT_STREQ("1", fmtd(POW(__DBL_MIN__, 0.)));
EXPECT_STREQ("1", fmtd(POW(__DBL_MIN__, -0.)));
EXPECT_STREQ("1.49166814624004e-154", fmtd(POW(__DBL_MIN__, .5)));
EXPECT_STREQ("6.7039039649713e+153", fmtd(POW(__DBL_MIN__, -.5)));
EXPECT_STREQ("2.2250738585072e-308", fmtd(POW(__DBL_MIN__, 1.)));
EXPECT_STREQ("4.49423283715579e+307", fmtd(POW(__DBL_MIN__, -1.)));
errno = 0;
EXPECT_STREQ("0", fmtd(pow(__DBL_MIN__, 1.5)));
EXPECT_STREQ("0", fmtd(POW(__DBL_MIN__, 1.5)));
/* EXPECT_EQ(ERANGE, errno); */
errno = 0;
EXPECT_STREQ("inf", fmtd(pow(__DBL_MIN__, -1.5)));
EXPECT_STREQ("inf", fmtd(POW(__DBL_MIN__, -1.5)));
/* EXPECT_EQ(ERANGE, errno); */
EXPECT_STREQ("0", fmtd(pow(__DBL_MIN__, 2.0)));
EXPECT_STREQ("inf", fmtd(pow(__DBL_MIN__, -2.0)));
EXPECT_TRUE(isnan(pow(__DBL_MIN__, NAN)));
EXPECT_TRUE(isnan(pow(__DBL_MIN__, -NAN)));
EXPECT_STREQ("0", fmtd(pow(__DBL_MIN__, INFINITY)));
EXPECT_STREQ("inf", fmtd(pow(__DBL_MIN__, -INFINITY)));
EXPECT_STREQ("1", fmtd(pow(__DBL_MIN__, __DBL_MIN__)));
EXPECT_STREQ("0", fmtd(POW(__DBL_MIN__, 2.0)));
EXPECT_STREQ("inf", fmtd(POW(__DBL_MIN__, -2.0)));
EXPECT_TRUE(isnan(POW(__DBL_MIN__, NAN)));
EXPECT_TRUE(isnan(POW(__DBL_MIN__, -NAN)));
EXPECT_STREQ("0", fmtd(POW(__DBL_MIN__, INFINITY)));
EXPECT_STREQ("inf", fmtd(POW(__DBL_MIN__, -INFINITY)));
EXPECT_STREQ("1", fmtd(POW(__DBL_MIN__, __DBL_MIN__)));
/* errno = 0; */ /* wut */
/* EXPECT_STREQ("0", fmtd(pow(__DBL_MIN__, __DBL_MAX__))); */
/* EXPECT_STREQ("0", fmtd(POW(__DBL_MIN__, __DBL_MAX__))); */
/* EXPECT_EQ(ERANGE, errno); */
EXPECT_STREQ("1", fmtd(pow(__DBL_MAX__, 0.)));
EXPECT_STREQ("1", fmtd(pow(__DBL_MAX__, -0.)));
EXPECT_STREQ("1.34078079299426e+154", fmtd(pow(__DBL_MAX__, .5)));
EXPECT_STREQ("7.45834073120021e-155", fmtd(pow(__DBL_MAX__, -.5)));
EXPECT_STREQ("1.79769313486232e+308", fmtd(pow(__DBL_MAX__, 1.)));
EXPECT_STREQ("5.562684646268e-309", fmtd(pow(__DBL_MAX__, -1.)));
EXPECT_STREQ("1", fmtd(POW(__DBL_MAX__, 0.)));
EXPECT_STREQ("1", fmtd(POW(__DBL_MAX__, -0.)));
EXPECT_STREQ("1.34078079299426e+154", fmtd(POW(__DBL_MAX__, .5)));
EXPECT_STREQ("7.45834073120021e-155", fmtd(POW(__DBL_MAX__, -.5)));
EXPECT_STREQ("1.79769313486232e+308", fmtd(POW(__DBL_MAX__, 1.)));
EXPECT_STREQ("5.562684646268e-309", fmtd(POW(__DBL_MAX__, -1.)));
errno = 0;
EXPECT_STREQ("inf", fmtd(pow(__DBL_MAX__, 1.5)));
EXPECT_STREQ("inf", fmtd(POW(__DBL_MAX__, 1.5)));
/* EXPECT_EQ(ERANGE, errno); */
errno = 0;
EXPECT_STREQ("0", fmtd(pow(__DBL_MAX__, -1.5)));
EXPECT_STREQ("0", fmtd(POW(__DBL_MAX__, -1.5)));
/* EXPECT_EQ(ERANGE, errno); */
EXPECT_STREQ("inf", fmtd(pow(__DBL_MAX__, 2.0)));
EXPECT_STREQ("inf", fmtd(POW(__DBL_MAX__, 2.0)));
errno = 0;
EXPECT_STREQ("0", fmtd(pow(__DBL_MAX__, -2.0)));
EXPECT_STREQ("0", fmtd(POW(__DBL_MAX__, -2.0)));
/* EXPECT_EQ(ERANGE, errno); */
EXPECT_TRUE(isnan(pow(__DBL_MAX__, NAN)));
EXPECT_TRUE(isnan(pow(__DBL_MAX__, -NAN)));
EXPECT_STREQ("inf", fmtd(pow(__DBL_MAX__, INFINITY)));
EXPECT_STREQ("0", fmtd(pow(__DBL_MAX__, -INFINITY)));
EXPECT_STREQ("1", fmtd(pow(__DBL_MAX__, __DBL_MIN__)));
EXPECT_TRUE(isnan(POW(__DBL_MAX__, NAN)));
EXPECT_TRUE(isnan(POW(__DBL_MAX__, -NAN)));
EXPECT_STREQ("inf", fmtd(POW(__DBL_MAX__, INFINITY)));
EXPECT_STREQ("0", fmtd(POW(__DBL_MAX__, -INFINITY)));
EXPECT_STREQ("1", fmtd(POW(__DBL_MAX__, __DBL_MIN__)));
errno = 0;
EXPECT_STREQ("inf", fmtd(pow(__DBL_MAX__, __DBL_MAX__)));
EXPECT_STREQ("inf", fmtd(POW(__DBL_MAX__, __DBL_MAX__)));
/* EXPECT_EQ(ERANGE, errno); */
EXPECT_STREQ("1", gc(xasprintf("%.15g", pow(0., 0))));
EXPECT_STREQ("1", gc(xasprintf("%.15g", pow(-0., 0))));
EXPECT_STREQ("-0", gc(xasprintf("%.15g", pow(-0., 1))));
EXPECT_STREQ("-0", gc(xasprintf("%.15g", pow(-0., 11))));
EXPECT_STREQ("-0", gc(xasprintf("%.15g", pow(-0., 111))));
EXPECT_STREQ("0", gc(xasprintf("%.15g", pow(-0., 2))));
EXPECT_STREQ("0", gc(xasprintf("%.15g", pow(-0., 22))));
EXPECT_STREQ("0", gc(xasprintf("%.15g", pow(-0., 222))));
EXPECT_STREQ("0", gc(xasprintf("%.15g", pow(-0., 2.5))));
EXPECT_STREQ("1", gc(xasprintf("%.15g", POW(0., 0))));
EXPECT_STREQ("1", gc(xasprintf("%.15g", POW(-0., 0))));
EXPECT_STREQ("-0", gc(xasprintf("%.15g", POW(-0., 1))));
EXPECT_STREQ("-0", gc(xasprintf("%.15g", POW(-0., 11))));
EXPECT_STREQ("-0", gc(xasprintf("%.15g", POW(-0., 111))));
EXPECT_STREQ("0", gc(xasprintf("%.15g", POW(-0., 2))));
EXPECT_STREQ("0", gc(xasprintf("%.15g", POW(-0., 22))));
EXPECT_STREQ("0", gc(xasprintf("%.15g", POW(-0., 222))));
EXPECT_STREQ("0", gc(xasprintf("%.15g", POW(-0., 2.5))));
}
BENCH(powl, bench) {
double _pow(double, double) asm("pow");
float _powf(float, float) asm("powf");
long double _powl(long double, long double) asm("powl");
EZBENCH2("pow", donothing, _pow(.7, .2)); /* ~56ns */
EZBENCH2("pow", donothing, _pow(.7, .2)); /* ~18ns */
EZBENCH2("powf", donothing, _powf(.7, .2)); /* ~56ns */
EZBENCH2("powl", donothing, _powl(.7, .2)); /* ~56ns */
}

10
test/libc/tinymath/trunc_test.c
View file

@ -18,6 +18,7 @@
╚─────────────────────────────────────────────────────────────────────────────*/
#include "libc/math.h"
#include "libc/runtime/gc.internal.h"
#include "libc/testlib/ezbench.h"
#include "libc/testlib/testlib.h"
#include "libc/x/x.h"
@ -53,3 +54,12 @@ TEST(truncl, test) {
EXPECT_STREQ("INFINITY", gc(xdtoal(truncl(INFINITY))));
EXPECT_STREQ("-INFINITY", gc(xdtoal(truncl(-INFINITY))));
}
BENCH(truncl, bench) {
double _trunc(double) asm("trunc");
float _truncf(float) asm("truncf");
long double _truncl(long double) asm("truncl");
EZBENCH2("trunc", donothing, _trunc(.7)); /* ~2ns */
EZBENCH2("truncf", donothing, _truncf(.7)); /* ~2ns */
EZBENCH2("truncl", donothing, _truncl(.7)); /* ~9ns */
}