Skip to content

feat(libc): implement extended-precision math exports #2271

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Closed
obojeoboje wants to merge 2 commits into from
Closed

feat(libc): implement extended-precision math exports #2271

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
2 commits
Select commit Hold shift + click to select a range
1b40c04
feat(libc): implement extended-precision math exports
Oct 9, 2026
91f4371
fix(libc): preserve extended extrema signs and neighbor exceptions
Oct 9, 2026
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
16 changes: 16 additions & 0 deletions core/libs/prx/libc/CMakeLists.txt
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -11,6 +11,7 @@ target_sources(libc PRIVATE
${CMAKE_CURRENT_SOURCE_DIR}/src/Multibyte.cpp
${CMAKE_CURRENT_SOURCE_DIR}/src/StringExtensions.cpp
${CMAKE_CURRENT_SOURCE_DIR}/src/Math.cpp
${CMAKE_CURRENT_SOURCE_DIR}/src/ExtendedMath.cpp
${CMAKE_CURRENT_SOURCE_DIR}/src/Time.cpp
${CMAKE_CURRENT_SOURCE_DIR}/src/FileAndHeap.cpp
${CMAKE_CURRENT_SOURCE_DIR}/src/StdioExtras.cpp
Expand All @@ -37,6 +38,9 @@ target_sources(libc PRIVATE
${CMAKE_CURRENT_SOURCE_DIR}/src/exception/Exports.cpp
)

set_source_files_properties(${CMAKE_CURRENT_SOURCE_DIR}/src/ExtendedMath.cpp
TARGET_DIRECTORY libc PROPERTIES COMPILE_OPTIONS "-frounding-math")

target_include_directories(libc PRIVATE
${LIBS_INCLUDE_DIR}
)
Expand All @@ -46,6 +50,18 @@ if(WIN32)
endif()

if(BUILD_TESTING)
add_test_executable(guest_extended_math_tests ${LIBS_INCLUDE_DIR}/tests/GuestExtendedMath.cpp)
target_include_directories(guest_extended_math_tests PRIVATE ${LIBS_INCLUDE_DIR})
target_link_libraries(guest_extended_math_tests PRIVATE libc)
configure_windows_unwind(guest_extended_math_tests)
set_source_files_properties(${LIBS_INCLUDE_DIR}/tests/GuestExtendedMath.cpp
PROPERTIES COMPILE_OPTIONS "-frounding-math")
add_test(NAME guest_extended_math COMMAND guest_extended_math_tests)
if(WIN32)
set_tests_properties(guest_extended_math PROPERTIES ENVIRONMENT_MODIFICATION
"PATH=path_list_prepend:$<TARGET_FILE_DIR:libc>")
endif()

add_test_executable(guest_time_conversion_tests ${LIBS_INCLUDE_DIR}/tests/GuestTimeConversion.cpp)
target_include_directories(guest_time_conversion_tests PRIVATE ${LIBS_INCLUDE_DIR})
target_link_libraries(guest_time_conversion_tests PRIVATE libc)
Expand Down
130 changes: 130 additions & 0 deletions core/libs/prx/libc/src/ExtendedMath.cpp
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,130 @@
#include "prx/libc/include/general/VabiMacros.hpp"

#include <algorithm>
#include <cerrno>
#include <cfenv>
#include <cmath>
#include <cstdint>
#include <limits>

static_assert(sizeof(long double) == 16 && std::numeric_limits<long double>::digits == 64);
static_assert(sizeof(long long) == sizeof(std::int64_t));

namespace {

long double DomainError_nid_no_patch() {
errno = EDOM;
std::feraiseexcept(FE_INVALID);
return std::numeric_limits<long double>::quiet_NaN();
}

long double CheckRange_nid_no_patch(long double input, long double result) {
if (std::isfinite(input)) {
if (std::isinf(result)) {
errno = ERANGE;
std::feraiseexcept(FE_OVERFLOW | FE_INEXACT);
} else if (result == 0.L && input != 0.L) {
errno = ERANGE;
std::feraiseexcept(FE_UNDERFLOW | FE_INEXACT);
}
}
return result;
}

long double CheckNeighbor_nid_no_patch(long double from, long double to, long double result) {
if (from == to || std::isnan(result)) return result;
if (std::isfinite(from) && std::isinf(result)) {
errno = ERANGE;
std::feraiseexcept(FE_OVERFLOW | FE_INEXACT);
} else if (result == 0.L || std::fpclassify(result) == FP_SUBNORMAL) {
errno = ERANGE;
std::feraiseexcept(FE_UNDERFLOW | FE_INEXACT);
}
return result;
}

}

extern "C" {

long double APS5_VABI acosl_nid_postfix(long double value) { return std::acos(value); }
long double APS5_VABI asinl_nid_postfix(long double value) { return std::asin(value); }
long double APS5_VABI atanl_nid_postfix(long double value) { return std::atan(value); }
long double APS5_VABI atan2l_nid_postfix(long double y, long double x) { return std::atan2(y, x); }
long double APS5_VABI tanl_nid_postfix(long double value) { return std::tan(value); }

long double APS5_VABI acoshl_nid_postfix(long double value) {
return value < 1.L ? DomainError_nid_no_patch() : std::acosh(value);
}
long double APS5_VABI asinhl_nid_postfix(long double value) { return std::asinh(value); }
long double APS5_VABI atanhl_nid_postfix(long double value) {
if (std::fabs(value) > 1.L) {
return DomainError_nid_no_patch();
}
if (std::fabs(value) == 1.L) {
errno = ERANGE;
std::feraiseexcept(FE_DIVBYZERO);
return std::copysign(std::numeric_limits<long double>::infinity(), value);
}
return std::atanh(value);
}
long double APS5_VABI coshl_nid_postfix(long double value) { return CheckRange_nid_no_patch(value, std::cosh(value)); }
long double APS5_VABI sinhl_nid_postfix(long double value) { return CheckRange_nid_no_patch(value, std::sinh(value)); }
long double APS5_VABI tanhl_nid_postfix(long double value) { return std::tanh(value); }

long double APS5_VABI expl_nid_postfix(long double value) { return CheckRange_nid_no_patch(value, std::exp(value)); }
long double APS5_VABI exp2l_nid_postfix(long double value) { return std::exp2(value); }
long double APS5_VABI expm1l_nid_postfix(long double value) { return CheckRange_nid_no_patch(value, std::expm1(value)); }
long double APS5_VABI log10l_nid_postfix(long double value) { return std::log10(value); }
long double APS5_VABI log1pl_nid_postfix(long double value) { return std::log1p(value); }
long double APS5_VABI log2l_nid_postfix(long double value) { return std::log2(value); }
long double APS5_VABI logbl_nid_postfix(long double value) { return std::logb(value); }

long double APS5_VABI sqrtl_nid_postfix(long double value) {
return value < 0.L ? DomainError_nid_no_patch() : std::sqrt(value);
}
long double APS5_VABI cbrtl_nid_postfix(long double value) { return std::cbrt(value); }
long double APS5_VABI hypotl_nid_postfix(long double x, long double y) { return std::hypot(x, y); }
long double APS5_VABI fabsl_nid_postfix(long double value) { return std::fabs(value); }
long double APS5_VABI copysignl_nid_postfix(long double magnitude, long double sign) { return std::copysign(magnitude, sign); }
long double APS5_VABI fdiml_nid_postfix(long double x, long double y) { return std::fdim(x, y); }
long double APS5_VABI fmaxl_nid_postfix(long double x, long double y) {
if (x == 0.L && y == 0.L) return std::signbit(x) ? y : x;
return std::fmax(x, y);
}
long double APS5_VABI fminl_nid_postfix(long double x, long double y) {
if (x == 0.L && y == 0.L) return std::signbit(x) ? x : y;
return std::fmin(x, y);
}
long double APS5_VABI fmal_nid_postfix(long double x, long double y, long double z) { return std::fma(x, y, z); }

long double APS5_VABI ceill_nid_postfix(long double value) { return std::ceil(value); }
long double APS5_VABI floorl_nid_postfix(long double value) { return std::floor(value); }
long double APS5_VABI roundl_nid_postfix(long double value) { return std::round(value); }
std::int64_t APS5_VABI lroundl_nid_postfix(long double value) { return std::llround(value); }
std::int64_t APS5_VABI llroundl_nid_postfix(long double value) { return std::llround(value); }
long double APS5_VABI nearbyintl_nid_postfix(long double value) { return std::nearbyint(value); }

long double APS5_VABI fmodl_nid_postfix(long double x, long double y) { return std::fmod(x, y); }
long double APS5_VABI remainderl_nid_postfix(long double x, long double y) { return std::remainder(x, y); }
long double APS5_VABI remquol_nid_postfix(long double x, long double y, int* quotient) { return std::remquo(x, y, quotient); }
long double APS5_VABI nextafterl_nid_postfix(long double from, long double to) { return CheckNeighbor_nid_no_patch(from, to, std::nextafter(from, to)); }
long double APS5_VABI nexttowardl_nid_postfix(long double from, long double to) { return CheckNeighbor_nid_no_patch(from, to, std::nexttoward(from, to)); }

long double APS5_VABI frexpl_nid_postfix(long double value, int* exponent) { return std::frexp(value, exponent); }
long double APS5_VABI ldexpl_nid_postfix(long double value, int exponent) { return std::ldexp(value, exponent); }
long double APS5_VABI modfl_nid_postfix(long double value, long double* integral) {
if (std::isinf(value)) {
*integral = value;
return std::copysign(0.L, value);
}
return std::modf(value, integral);
}
long double APS5_VABI scalbnl_nid_postfix(long double value, int exponent) { return std::scalbn(value, exponent); }
long double APS5_VABI scalblnl_nid_postfix(long double value, std::int64_t exponent) {
const auto bounded = std::clamp(exponent, static_cast<std::int64_t>(std::numeric_limits<int>::min()),
static_cast<std::int64_t>(std::numeric_limits<int>::max()));
return std::scalbn(value, static_cast<int>(bounded));
}

}
Loading
Loading