@zigc/lib 0.16.0-test.1 → 0.17.0-dev.27

package/libc/musl/src/math/i386/rintf.c

@@ -1,7 +0,0 @@
-#include <math.h>
-
-float rintf(float x)
-{
-	__asm__ ("frndint" : "+t"(x));
-	return x;
-}

package/libc/musl/src/math/lrint.c

@@ -1,72 +0,0 @@
-#include <limits.h>
-#include <fenv.h>
-#include <math.h>
-#include "libm.h"
-
-/*
-If the result cannot be represented (overflow, nan), then
-lrint raises the invalid exception.
-
-Otherwise if the input was not an integer then the inexact
-exception is raised.
-
-C99 is a bit vague about whether inexact exception is
-allowed to be raised when invalid is raised.
-(F.9 explicitly allows spurious inexact exceptions, F.9.6.5
-does not make it clear if that rule applies to lrint, but
-IEEE 754r 7.8 seems to forbid spurious inexact exception in
-the ineger conversion functions)
-
-So we try to make sure that no spurious inexact exception is
-raised in case of an overflow.
-
-If the bit size of long > precision of double, then there
-cannot be inexact rounding in case the result overflows,
-otherwise LONG_MAX and LONG_MIN can be represented exactly
-as a double.
-*/
-
-#if LONG_MAX < 1U<<53 && defined(FE_INEXACT)
-#include <float.h>
-#include <stdint.h>
-#if FLT_EVAL_METHOD==0 || FLT_EVAL_METHOD==1
-#define EPS DBL_EPSILON
-#elif FLT_EVAL_METHOD==2
-#define EPS LDBL_EPSILON
-#endif
-#ifdef __GNUC__
-/* avoid stack frame in lrint */
-__attribute__((noinline))
-#endif
-static long lrint_slow(double x)
-{
-	#pragma STDC FENV_ACCESS ON
-	int e;
-
-	e = fetestexcept(FE_INEXACT);
-	x = rint(x);
-	if (!e && (x > LONG_MAX || x < LONG_MIN))
-		feclearexcept(FE_INEXACT);
-	/* conversion */
-	return x;
-}
-
-long lrint(double x)
-{
-	uint32_t abstop = asuint64(x)>>32 & 0x7fffffff;
-	uint64_t sign = asuint64(x) & (1ULL << 63);
-
-	if (abstop < 0x41dfffff) {
-		/* |x| < 0x7ffffc00, no overflow */
-		double_t toint = asdouble(asuint64(1/EPS) | sign);
-		double_t y = x + toint - toint;
-		return (long)y;
-	}
-	return lrint_slow(x);
-}
-#else
-long lrint(double x)
-{
-	return rint(x);
-}
-#endif

package/libc/musl/src/math/lrintf.c

@@ -1,8 +0,0 @@
-#include <math.h>
-
-/* uses LONG_MAX > 2^24, see comments in lrint.c */
-
-long lrintf(float x)
-{
-	return rintf(x);
-}

package/libc/musl/src/math/powerpc64/lrint.c

@@ -1,16 +0,0 @@
-#include <math.h>
-
-#ifdef _ARCH_PWR5X
-
-long lrint(double x)
-{
-	long n;
-	__asm__ ("fctid %0, %1" : "=d"(n) : "d"(x));
-	return n;
-}
-
-#else
-
-#include "../lrint.c"
-
-#endif

package/libc/musl/src/math/powerpc64/lrintf.c

@@ -1,16 +0,0 @@
-#include <math.h>
-
-#ifdef _ARCH_PWR5X
-
-long lrintf(float x)
-{
-	long n;
-	__asm__ ("fctid %0, %1" : "=d"(n) : "f"(x));
-	return n;
-}
-
-#else
-
-#include "../lrintf.c"
-
-#endif

package/libc/musl/src/math/rintf.c

@@ -1,30 +0,0 @@
-#include <float.h>
-#include <math.h>
-#include <stdint.h>
-
-#if FLT_EVAL_METHOD==0
-#define EPS FLT_EPSILON
-#elif FLT_EVAL_METHOD==1
-#define EPS DBL_EPSILON
-#elif FLT_EVAL_METHOD==2
-#define EPS LDBL_EPSILON
-#endif
-static const float_t toint = 1/EPS;
-
-float rintf(float x)
-{
-	union {float f; uint32_t i;} u = {x};
-	int e = u.i>>23 & 0xff;
-	int s = u.i>>31;
-	float_t y;
-
-	if (e >= 0x7f+23)
-		return x;
-	if (s)
-		y = x - toint + toint;
-	else
-		y = x + toint - toint;
-	if (y == 0)
-		return s ? -0.0f : 0.0f;
-	return y;
-}

package/libc/musl/src/math/s390x/rintf.c

@@ -1,15 +0,0 @@
-#include <math.h>
-
-#if defined(__HTM__) || __ARCH__ >= 9
-
-float rintf(float x)
-{
-	__asm__ ("fiebr %0, 0, %1" : "=f"(x) : "f"(x));
-	return x;
-}
-
-#else
-
-#include "../rintf.c"
-
-#endif

package/libc/musl/src/math/sincosl.c

@@ -1,60 +0,0 @@
-#define _GNU_SOURCE
-#include "libm.h"
-
-#if LDBL_MANT_DIG == 53 && LDBL_MAX_EXP == 1024
-void sincosl(long double x, long double *sin, long double *cos)
-{
-	double sind, cosd;
-	sincos(x, &sind, &cosd);
-	*sin = sind;
-	*cos = cosd;
-}
-#elif (LDBL_MANT_DIG == 64 || LDBL_MANT_DIG == 113) && LDBL_MAX_EXP == 16384
-void sincosl(long double x, long double *sin, long double *cos)
-{
-	union ldshape u = {x};
-	unsigned n;
-	long double y[2], s, c;
-
-	u.i.se &= 0x7fff;
-	if (u.i.se == 0x7fff) {
-		*sin = *cos = x - x;
-		return;
-	}
-	if (u.f < M_PI_4) {
-		if (u.i.se < 0x3fff - LDBL_MANT_DIG) {
-			/* raise underflow if subnormal */
-			if (u.i.se == 0) FORCE_EVAL(x*0x1p-120f);
-			*sin = x;
-			/* raise inexact if x!=0 */
-			*cos = 1.0 + x;
-			return;
-		}
-		*sin = __sinl(x, 0, 0);
-		*cos = __cosl(x, 0);
-		return;
-	}
-	n = __rem_pio2l(x, y);
-	s = __sinl(y[0], y[1], 1);
-	c = __cosl(y[0], y[1]);
-	switch (n & 3) {
-	case 0:
-		*sin = s;
-		*cos = c;
-		break;
-	case 1:
-		*sin = c;
-		*cos = -s;
-		break;
-	case 2:
-		*sin = -s;
-		*cos = -c;
-		break;
-	case 3:
-	default:
-		*sin = -c;
-		*cos = s;
-		break;
-	}
-}
-#endif

package/libc/musl/src/math/sinl.c

@@ -1,41 +0,0 @@
-#include "libm.h"
-
-#if LDBL_MANT_DIG == 53 && LDBL_MAX_EXP == 1024
-long double sinl(long double x)
-{
-	return sin(x);
-}
-#elif (LDBL_MANT_DIG == 64 || LDBL_MANT_DIG == 113) && LDBL_MAX_EXP == 16384
-long double sinl(long double x)
-{
-	union ldshape u = {x};
-	unsigned n;
-	long double y[2], hi, lo;
-
-	u.i.se &= 0x7fff;
-	if (u.i.se == 0x7fff)
-		return x - x;
-	if (u.f < M_PI_4) {
-		if (u.i.se < 0x3fff - LDBL_MANT_DIG/2) {
-			/* raise inexact if x!=0 and underflow if subnormal */
-			FORCE_EVAL(u.i.se == 0 ? x*0x1p-120f : x+0x1p120f);
-			return x;
-		}
-		return __sinl(x, 0.0, 0);
-	}
-	n = __rem_pio2l(x, y);
-	hi = y[0];
-	lo = y[1];
-	switch (n & 3) {
-	case 0:
-		return __sinl(hi, lo, 1);
-	case 1:
-		return __cosl(hi, lo);
-	case 2:
-		return -__sinl(hi, lo, 1);
-	case 3:
-	default:
-		return -__cosl(hi, lo);
-	}
-}
-#endif

package/libc/musl/src/math/tanl.c

@@ -1,29 +0,0 @@
-#include "libm.h"
-
-#if LDBL_MANT_DIG == 53 && LDBL_MAX_EXP == 1024
-long double tanl(long double x)
-{
-	return tan(x);
-}
-#elif (LDBL_MANT_DIG == 64 || LDBL_MANT_DIG == 113) && LDBL_MAX_EXP == 16384
-long double tanl(long double x)
-{
-	union ldshape u = {x};
-	long double y[2];
-	unsigned n;
-
-	u.i.se &= 0x7fff;
-	if (u.i.se == 0x7fff)
-		return x - x;
-	if (u.f < M_PI_4) {
-		if (u.i.se < 0x3fff - LDBL_MANT_DIG/2) {
-			/* raise inexact if x!=0 and underflow if subnormal */
-			FORCE_EVAL(u.i.se == 0 ? x*0x1p-120f : x+0x1p120f);
-			return x;
-		}
-		return __tanl(x, 0, 0);
-	}
-	n = __rem_pio2l(x, y);
-	return __tanl(y[0], y[1], n&1);
-}
-#endif

package/libc/musl/src/math/x32/lrint.s

@@ -1,5 +0,0 @@
-.global lrint
-.type lrint,@function
-lrint:
-	cvtsd2si %xmm0,%rax
-	ret

package/libc/musl/src/math/x32/lrintf.s

@@ -1,5 +0,0 @@
-.global lrintf
-.type lrintf,@function
-lrintf:
-	cvtss2si %xmm0,%rax
-	ret

package/libc/musl/src/math/x86_64/lrint.c

@@ -1,8 +0,0 @@
-#include <math.h>
-
-long lrint(double x)
-{
-	long r;
-	__asm__ ("cvtsd2si %1, %0" : "=r"(r) : "x"(x));
-	return r;
-}

package/libc/musl/src/math/x86_64/lrintf.c

@@ -1,8 +0,0 @@
-#include <math.h>
-
-long lrintf(float x)
-{
-	long r;
-	__asm__ ("cvtss2si %1, %0" : "=r"(r) : "x"(x));
-	return r;
-}

package/libc/wasi/libc-bottom-half/sources/reallocarray.c

@@ -1,14 +0,0 @@
-#include <stdlib.h>
-#include <errno.h>
-
-void *__reallocarray(void *ptr, size_t nmemb, size_t size) {
-    size_t bytes;
-    if (__builtin_umull_overflow(nmemb, size, &bytes)) {
-        errno = ENOMEM;
-        return NULL;
-    }
-    return realloc(ptr, bytes);
-}
-
-void *reallocarray(void *ptr, size_t nmemb, size_t size)
-    __attribute__((__weak__, __alias__("__reallocarray")));