[libc-commits] [libc] [libc][math][c23] Implement canonicalize functions (PR #85940)

[via libc-commits]

================
@@ -73,6 +76,54 @@ LIBC_INLINE T fdim(T x, T y) {
   return (x > y ? x - y : 0);
 }
 
+template <typename T, cpp::enable_if_t<cpp::is_floating_point_v<T>, int> = 0>
+LIBC_INLINE int canonicalize(T &cx, const T &x) {
+  FPBits<T> sx(x);
+  if (LIBC_UNLIKELY(sx.is_signaling_nan())) {
+    cx = FPBits<T>::quiet_nan(sx.sign(), sx.get_explicit_mantissa()).get_val();
----------------
lntue wrote:

Does it work correctly for 80-bit extended precision?

https://github.com/llvm/llvm-project/pull/85940