mirror of
https://github.com/boykopovar/AnyPS5.git
synced 2026-10-22 00:22:57 -05:00
hypermistandClaude Opus 5.5
f8b611b7da
Soft-float arithmetic and comparisons, float/double/x87 conversions to and from 32-, 64- and 128-bit integers, __extendsfdf2, __truncdfsf2, complex multiplication and division and __powixf2, matching the libc.prx bundled with PPSA08427: soft-float results ignore MXCSR and raise no flags, division flushes subnormal quotients, NaN payloads are moved unshifted and float-to-integer conversions do not saturate. Co-Authored-By: Claude Opus 5.5 <noreply@anthropic.com> Claude-Session: https://claude.ai/code/session_01Q8tinbP1MUpUaDsbNWnGnu
229 lines
13 KiB
C++
229 lines
13 KiB
C++
#include "prx/libc/include/general/VabiMacros.hpp"
|
|
#include <bit>
|
|
#include <cmath>
|
|
#include <cstdint>
|
|
#include <cstdlib>
|
|
#include <cstring>
|
|
#include <xmmintrin.h>
|
|
|
|
using Int128 = __int128;
|
|
using Uint128 = unsigned __int128;
|
|
|
|
extern "C" {
|
|
double APS5_VABI __adddf3_nid_postfix(double, double);
|
|
float APS5_VABI __addsf3_nid_postfix(float, float);
|
|
double APS5_VABI __subdf3_nid_postfix(double, double);
|
|
float APS5_VABI __subsf3_nid_postfix(float, float);
|
|
double APS5_VABI __muldf3_nid_postfix(double, double);
|
|
float APS5_VABI __mulsf3_nid_postfix(float, float);
|
|
double APS5_VABI __divdf3_nid_postfix(double, double);
|
|
float APS5_VABI __divsf3_nid_postfix(float, float);
|
|
double APS5_VABI __negdf2_nid_postfix(double);
|
|
float APS5_VABI __negsf2_nid_postfix(float);
|
|
int APS5_VABI __eqdf2_nid_postfix(double, double);
|
|
int APS5_VABI __nedf2_nid_postfix(double, double);
|
|
int APS5_VABI __ltdf2_nid_postfix(double, double);
|
|
int APS5_VABI __ledf2_nid_postfix(double, double);
|
|
int APS5_VABI __gtdf2_nid_postfix(double, double);
|
|
int APS5_VABI __gedf2_nid_postfix(double, double);
|
|
int APS5_VABI __unorddf2_nid_postfix(double, double);
|
|
int APS5_VABI __eqsf2_nid_postfix(float, float);
|
|
int APS5_VABI __ltsf2_nid_postfix(float, float);
|
|
int APS5_VABI __gesf2_nid_postfix(float, float);
|
|
int APS5_VABI __unordsf2_nid_postfix(float, float);
|
|
double APS5_VABI __extendsfdf2_nid_postfix(float);
|
|
float APS5_VABI __truncdfsf2_nid_postfix(double);
|
|
std::int32_t APS5_VABI __fixsfsi_nid_postfix(float);
|
|
std::int64_t APS5_VABI __fixsfdi_nid_postfix(float);
|
|
Int128 APS5_VABI __fixsfti_nid_postfix(float);
|
|
std::int32_t APS5_VABI __fixdfsi_nid_postfix(double);
|
|
std::int64_t APS5_VABI __fixdfdi_nid_postfix(double);
|
|
Int128 APS5_VABI __fixdfti_nid_postfix(double);
|
|
std::uint32_t APS5_VABI __fixunssfsi_nid_postfix(float);
|
|
std::uint64_t APS5_VABI __fixunssfdi_nid_postfix(float);
|
|
Uint128 APS5_VABI __fixunssfti_nid_postfix(float);
|
|
std::uint32_t APS5_VABI __fixunsdfsi_nid_postfix(double);
|
|
std::uint64_t APS5_VABI __fixunsdfdi_nid_postfix(double);
|
|
Uint128 APS5_VABI __fixunsdfti_nid_postfix(double);
|
|
std::int64_t APS5_VABI __fixxfdi_nid_postfix(long double);
|
|
Int128 APS5_VABI __fixxfti_nid_postfix(long double);
|
|
std::uint32_t APS5_VABI __fixunsxfsi_nid_postfix(long double);
|
|
std::uint64_t APS5_VABI __fixunsxfdi_nid_postfix(long double);
|
|
Uint128 APS5_VABI __fixunsxfti_nid_postfix(long double);
|
|
double APS5_VABI __floatsidf_nid_postfix(std::int32_t);
|
|
double APS5_VABI __floatunsidf_nid_postfix(std::uint32_t);
|
|
float APS5_VABI __floatsisf_nid_postfix(std::int32_t);
|
|
float APS5_VABI __floatunsisf_nid_postfix(std::uint32_t);
|
|
double APS5_VABI __floatdidf_nid_postfix(std::int64_t);
|
|
float APS5_VABI __floatdisf_nid_postfix(std::int64_t);
|
|
long double APS5_VABI __floatdixf_nid_postfix(std::int64_t);
|
|
double APS5_VABI __floatundidf_nid_postfix(std::uint64_t);
|
|
float APS5_VABI __floatundisf_nid_postfix(std::uint64_t);
|
|
long double APS5_VABI __floatundixf_nid_postfix(std::uint64_t);
|
|
double APS5_VABI __floattidf_nid_postfix(Int128);
|
|
float APS5_VABI __floattisf_nid_postfix(Int128);
|
|
long double APS5_VABI __floattixf_nid_postfix(Int128);
|
|
double APS5_VABI __floatuntidf_nid_postfix(Uint128);
|
|
float APS5_VABI __floatuntisf_nid_postfix(Uint128);
|
|
long double APS5_VABI __floatuntixf_nid_postfix(Uint128);
|
|
__complex__ float APS5_VABI __mulsc3_nid_postfix(float, float, float, float);
|
|
__complex__ double APS5_VABI __muldc3_nid_postfix(double, double, double, double);
|
|
__complex__ long double APS5_VABI __mulxc3_nid_postfix(long double, long double, long double, long double);
|
|
__complex__ float APS5_VABI __divsc3_nid_postfix(float, float, float, float);
|
|
__complex__ double APS5_VABI __divdc3_nid_postfix(double, double, double, double);
|
|
__complex__ long double APS5_VABI __divxc3_nid_postfix(long double, long double, long double, long double);
|
|
long double APS5_VABI __powixf2_nid_postfix(long double, int);
|
|
}
|
|
|
|
static void Require(bool value) { if (!value) std::abort(); }
|
|
|
|
static std::uint64_t Bits(double value) { return std::bit_cast<std::uint64_t>(value); }
|
|
static std::uint32_t Bits(float value) { return std::bit_cast<std::uint32_t>(value); }
|
|
static double Double(std::uint64_t bits) { return std::bit_cast<double>(bits); }
|
|
static float Float(std::uint32_t bits) { return std::bit_cast<float>(bits); }
|
|
static Uint128 Make(std::uint64_t high, std::uint64_t low) { return (static_cast<Uint128>(high) << 64) | low; }
|
|
|
|
static constexpr unsigned GuestMxcsr = 0x1f80 | 0x6000 | 0x8040;
|
|
|
|
static void CheckArithmetic() {
|
|
Require(__adddf3_nid_postfix(1., 2.) == 3. && __subsf3_nid_postfix(1.f, 2.f) == -1.f);
|
|
Require(__muldf3_nid_postfix(3., -4.) == -12. && __divsf3_nid_postfix(1.f, 4.f) == .25f);
|
|
Require(Bits(__adddf3_nid_postfix(INFINITY, -INFINITY)) == 0x7ff8000000000000ull);
|
|
Require(Bits(__muldf3_nid_postfix(0., INFINITY)) == 0x7ff8000000000000ull);
|
|
Require(Bits(__divdf3_nid_postfix(0., 0.)) == 0x7ff8000000000000ull);
|
|
Require(Bits(__subsf3_nid_postfix(INFINITY, INFINITY)) == 0x7fc00000u);
|
|
Require(Bits(__mulsf3_nid_postfix(0.f, -INFINITY)) == 0x7fc00000u);
|
|
Require(Bits(__adddf3_nid_postfix(Double(0xfff0000000000001ull), 1.)) == 0xfff8000000000001ull);
|
|
Require(Bits(__adddf3_nid_postfix(Double(0x7ff8000000000002ull), Double(0x7ff0000000000003ull))) == 0x7ff8000000000002ull);
|
|
Require(Bits(__muldf3_nid_postfix(2., Double(0x7ff0000000000003ull))) == 0x7ff8000000000003ull);
|
|
Require(Bits(__subdf3_nid_postfix(1., Double(0x7ff8000000000004ull))) == 0xfff8000000000004ull);
|
|
Require(Bits(__negdf2_nid_postfix(Double(0x7ff8000000000004ull))) == 0xfff8000000000004ull);
|
|
Require(Bits(__negsf2_nid_postfix(0.f)) == 0x80000000u);
|
|
Require(Bits(__divdf3_nid_postfix(-0x1p-1000, 0x1p30)) == 0x8000000000000000ull);
|
|
Require(Bits(__divsf3_nid_postfix(0x1p-120f, 0x1p10f)) == 0u);
|
|
Require(Bits(__divdf3_nid_postfix(0x1p-1022, 0x1.0000000000001p0)) == 0u);
|
|
Require(__divdf3_nid_postfix(0x1p-1000, 0x1p-22) == 0x1p-978 && __divsf3_nid_postfix(1.f, 3.f) == 1.f / 3.f);
|
|
const unsigned saved = _mm_getcsr();
|
|
_mm_setcsr(GuestMxcsr);
|
|
const double sum = __adddf3_nid_postfix(1., 0x1.8p-53);
|
|
const float denormal = __mulsf3_nid_postfix(0x1p-126f, .5f);
|
|
const double quotient = __divdf3_nid_postfix(1., 3.);
|
|
const unsigned after = _mm_getcsr();
|
|
_mm_setcsr(saved);
|
|
Require(sum == 1. + 0x1p-52);
|
|
Require(Bits(denormal) == 0x00400000u);
|
|
Require(Bits(quotient) == 0x3fd5555555555555ull);
|
|
Require(after == GuestMxcsr);
|
|
}
|
|
|
|
static void CheckComparisons() {
|
|
Require(__ledf2_nid_postfix(1., 2.) == -1 && __ledf2_nid_postfix(2., 1.) == 1 && __ledf2_nid_postfix(2., 2.) == 0);
|
|
Require(__ltdf2_nid_postfix(-2., -1.) == -1 && __gtdf2_nid_postfix(-1., -2.) == 1);
|
|
Require(__eqdf2_nid_postfix(0., -0.) == 0 && __nedf2_nid_postfix(1., -1.) != 0);
|
|
Require(__ledf2_nid_postfix(NAN, 1.) == 1 && __gedf2_nid_postfix(1., NAN) == -1 && __gtdf2_nid_postfix(NAN, NAN) == -1);
|
|
Require(__unorddf2_nid_postfix(1., NAN) == 1 && __unorddf2_nid_postfix(1., INFINITY) == 0);
|
|
Require(__eqsf2_nid_postfix(-0.f, 0.f) == 0 && __ltsf2_nid_postfix(-INFINITY, -1e30f) == -1);
|
|
Require(__gesf2_nid_postfix(NAN, 0.f) == -1 && __unordsf2_nid_postfix(NAN, 0.f) == 1);
|
|
Require(__ledf2_nid_postfix(Double(0x0000000000000001ull), 0.) == 1);
|
|
}
|
|
|
|
static void CheckFloatConversions() {
|
|
Require(Bits(__extendsfdf2_nid_postfix(Float(0x7f800001u))) == 0x7ff0000000000001ull);
|
|
Require(Bits(__extendsfdf2_nid_postfix(Float(0xffc00001u))) == 0xfff8000000000001ull);
|
|
Require(__extendsfdf2_nid_postfix(Float(1u)) == std::ldexp(1., -149));
|
|
Require(__extendsfdf2_nid_postfix(Float(0x80400000u)) == -std::ldexp(1., -127));
|
|
Require(__extendsfdf2_nid_postfix(1.5f) == 1.5 && __extendsfdf2_nid_postfix(-INFINITY) == -INFINITY);
|
|
Require(Bits(__truncdfsf2_nid_postfix(Double(0x7ff0000020000001ull))) == 0x7fc00001u);
|
|
Require(Bits(__truncdfsf2_nid_postfix(Double(0xfff4000000123456ull))) == 0xffd23456u);
|
|
const unsigned saved = _mm_getcsr();
|
|
_mm_setcsr(GuestMxcsr);
|
|
const float third = __truncdfsf2_nid_postfix(1. / 3.);
|
|
const float tiny = __truncdfsf2_nid_postfix(0x1p-140);
|
|
_mm_setcsr(saved);
|
|
Require(Bits(third) == 0x3eaaaaabu && Bits(tiny) == 0x00000200u);
|
|
}
|
|
|
|
static void CheckFloatToInteger() {
|
|
Require(__fixdfsi_nid_postfix(-3.7) == -3 && __fixsfsi_nid_postfix(.99f) == 0);
|
|
Require(__fixdfsi_nid_postfix(0x1p31) == INT32_MIN && __fixdfsi_nid_postfix(1e10) == 1410065408);
|
|
Require(__fixsfsi_nid_postfix(Float(0x7fc00000u)) == INT32_MIN);
|
|
Require(__fixdfdi_nid_postfix(0x1p63) == INT64_MIN && __fixdfdi_nid_postfix(-0x1p62) == -0x4000000000000000);
|
|
Require(__fixsfdi_nid_postfix(-0x1p40f) == -(1ll << 40));
|
|
Require(__fixdfti_nid_postfix(0x1p100) == static_cast<Int128>(Make(1ull << 36, 0)));
|
|
Require(__fixsfti_nid_postfix(-0x1p70f) == -static_cast<Int128>(Make(1ull << 6, 0)));
|
|
Require(__fixdfdi_nid_postfix(0x1p64) == 0 && __fixdfdi_nid_postfix(0x1p120) == 1ll << 56);
|
|
Require(__fixdfti_nid_postfix(0x1p200) == static_cast<Int128>(1) << 72);
|
|
Require(__fixunsdfsi_nid_postfix(-1.) == 0 && __fixunsdfsi_nid_postfix(4294967295.) == 0xffffffffu);
|
|
Require(__fixunsdfsi_nid_postfix(0x1p32) == 1);
|
|
Require(__fixunssfsi_nid_postfix(0x1p32f) == 0 && __fixunssfsi_nid_postfix(0x1p55f) == 1u << 23 && __fixunssfdi_nid_postfix(0x1p63f) == 1ull << 63);
|
|
Require(__fixunsdfdi_nid_postfix(18446744073709549568.) == 18446744073709549568ull);
|
|
Require(__fixunssfti_nid_postfix(0x1p127f) == Make(1ull << 63, 0) && __fixunsdfti_nid_postfix(-2.) == 0);
|
|
Require(__fixxfdi_nid_postfix(-5.5L) == -5 && __fixxfdi_nid_postfix(0x1p63L) == INT64_MIN);
|
|
Require(__fixunsxfdi_nid_postfix(0x1p64L) == 1 && __fixunsxfdi_nid_postfix(18446744073709551615.L) == ~0ull);
|
|
Require(__fixxfti_nid_postfix(-0x1p70L) == -static_cast<Int128>(Make(1ull << 6, 0)));
|
|
Require(__fixunsxfti_nid_postfix(0x1.8p100L) == Make(3ull << 35, 0) && __fixunsxfti_nid_postfix(-1.L) == 0);
|
|
Require(__fixunsxfsi_nid_postfix(3e9L) == 3000000000u && __fixunsxfsi_nid_postfix(0x1p32L) == 1);
|
|
}
|
|
|
|
static void CheckIntegerToFloat() {
|
|
Require(__floatsidf_nid_postfix(INT32_MIN) == -0x1p31 && __floatunsidf_nid_postfix(0xffffffffu) == 4294967295.);
|
|
Require(__floatdidf_nid_postfix(-3) == -3. && __floatdisf_nid_postfix(1ll << 40) == 0x1p40f);
|
|
Require(__floatundidf_nid_postfix(~0ull) == 0x1p64 && __floatundisf_nid_postfix(~0ull) == 0x1p64f);
|
|
Require(__floatdixf_nid_postfix(INT64_MIN) == -0x1p63L && __floatundixf_nid_postfix(~0ull) == 18446744073709551615.L);
|
|
const unsigned saved = _mm_getcsr();
|
|
_mm_setcsr(GuestMxcsr);
|
|
const float tie = __floatsisf_nid_postfix(16777217);
|
|
const float up = __floatunsisf_nid_postfix(16777219u);
|
|
const float negative = __floatsisf_nid_postfix(-16777219);
|
|
const double wide = __floattidf_nid_postfix(static_cast<Int128>(Make(0, (1ull << 63) | 0x601)));
|
|
_mm_setcsr(saved);
|
|
Require(tie == 16777216.f && up == 16777220.f && negative == -16777220.f);
|
|
Require(wide == 0x1p63 + 0x800);
|
|
Require(__floattidf_nid_postfix(static_cast<Int128>(Make(1, 1))) == 0x1p64);
|
|
Require(__floattisf_nid_postfix(-static_cast<Int128>(Make(1ull << 40, 0))) == -0x1p104f);
|
|
Require(__floatuntisf_nid_postfix(Make(~0ull, ~0ull)) == INFINITY);
|
|
Require(__floatuntidf_nid_postfix(Make(~0ull, ~0ull)) == 0x1p128);
|
|
Require(__floatuntidf_nid_postfix(0) == 0. && !std::signbit(__floattidf_nid_postfix(0)));
|
|
Require(__floattixf_nid_postfix(-static_cast<Int128>(Make(1, 1))) == -0x1p64L);
|
|
Require(__floattixf_nid_postfix(-static_cast<Int128>(Make(1, 3))) == -(0x1p64L + 4));
|
|
Require(__floatuntixf_nid_postfix(Make(~0ull, ~0ull)) == 0x1p128L);
|
|
Require(__floatuntixf_nid_postfix(Make(0, ~0ull)) == 18446744073709551615.L);
|
|
}
|
|
|
|
static void CheckComplex() {
|
|
const __complex__ float product = __mulsc3_nid_postfix(1.f, 2.f, 3.f, 4.f);
|
|
Require(__real__ product == -5.f && __imag__ product == 10.f);
|
|
const __complex__ double infinite = __muldc3_nid_postfix(INFINITY, NAN, 1., 0.);
|
|
Require(__real__ infinite == INFINITY);
|
|
const __complex__ long double extended = __mulxc3_nid_postfix(1.L, 1.L, 1.L, -1.L);
|
|
Require(__real__ extended == 2.L && __imag__ extended == 0.L);
|
|
const __complex__ double quotient = __divdc3_nid_postfix(6., 8., 3., 4.);
|
|
Require(__real__ quotient == 2. && __imag__ quotient == 0.);
|
|
const __complex__ float single = __divsc3_nid_postfix(1.f, 0.f, 0.f, 1.f);
|
|
Require(__real__ single == 0.f && __imag__ single == -1.f);
|
|
const __complex__ double byZero = __divdc3_nid_postfix(1., 0., 0., 0.);
|
|
Require(__real__ byZero == INFINITY);
|
|
const __complex__ long double large = __divxc3_nid_postfix(1e4000L, 1e4000L, 1e4000L, 1e4000L);
|
|
Require(__real__ large == 1.L && __imag__ large == 0.L);
|
|
const __complex__ double byInfinity = __divdc3_nid_postfix(1., 1., INFINITY, INFINITY);
|
|
Require(__real__ byInfinity == 0. && __imag__ byInfinity == 0.);
|
|
}
|
|
|
|
int main() {
|
|
CheckArithmetic();
|
|
CheckComparisons();
|
|
CheckFloatConversions();
|
|
CheckFloatToInteger();
|
|
CheckIntegerToFloat();
|
|
CheckComplex();
|
|
Require(__powixf2_nid_postfix(2.L, 10) == 1024.L && __powixf2_nid_postfix(2.L, -2) == .25L);
|
|
Require(__powixf2_nid_postfix(-3.L, 3) == -27.L && __powixf2_nid_postfix(5.L, 0) == 1.L);
|
|
long double signaling;
|
|
const unsigned char signalingBytes[10] = {1, 0, 0, 0, 0, 0, 0, 0x80, 0xff, 0x7f};
|
|
std::memset(&signaling, 0, sizeof(signaling));
|
|
std::memcpy(&signaling, signalingBytes, sizeof(signalingBytes));
|
|
const long double same = __powixf2_nid_postfix(signaling, 1);
|
|
Require(std::memcmp(&same, signalingBytes, sizeof(signalingBytes)) == 0);
|
|
return 0;
|
|
}
|