use libm functions

Author: malezjaa

library/std/src/num/f16.rs

@@ -832,13 +832,7 @@ impl f16 {
     #[unstable(feature = "f16", issue = "116909")]
     #[must_use = "method returns a new number and does not mutate the original value"]
     pub fn asinh(self) -> f16 {
-        let ax = self.abs();
-        if ax >= (1u16 << f16::MANTISSA_DIGITS / 2) as f16 {
-            return (ax.ln() + consts::LN_2).copysign(self);
-        }
-
-        let ix = 1.0 / ax;
-        (ax + (ax / (Self::hypot(1.0, ix) + ix))).ln_1p().copysign(self)
+        cmath::asinhf(self as f32) as f16
     }
 
     /// Inverse hyperbolic cosine function.
@@ -869,13 +863,7 @@ impl f16 {
     #[unstable(feature = "f16", issue = "116909")]
     #[must_use = "method returns a new number and does not mutate the original value"]
     pub fn acosh(self) -> f16 {
-        if self < 1.0 {
-            Self::NAN
-        } else if self >= (1u16 << f16::MANTISSA_DIGITS / 2) as f16 {
-            self.ln() + consts::LN_2
-        } else {
-            (self + ((self - 1.0).sqrt() * (self + 1.0).sqrt())).ln()
-        }
+        cmath::acoshf(self as f32) as f16
     }
 
     /// Inverse hyperbolic tangent function.

library/std/src/num/f32.rs

@@ -1091,13 +1091,7 @@ impl f32 {
     #[stable(feature = "rust1", since = "1.0.0")]
     #[inline]
     pub fn asinh(self) -> f32 {
-        let ax = self.abs();
-        if ax >= (1u32 << f32::MANTISSA_DIGITS / 2) as f32 {
-            return (ax.ln() + consts::LN_2).copysign(self);
-        }
-
-        let ix = 1.0 / ax;
-        (ax + (ax / (Self::hypot(1.0, ix) + ix))).ln_1p().copysign(self)
+        cmath::asinhf(self)
     }
 
     /// Inverse hyperbolic cosine function.
@@ -1123,13 +1117,7 @@ impl f32 {
     #[stable(feature = "rust1", since = "1.0.0")]
     #[inline]
     pub fn acosh(self) -> f32 {
-        if self < 1.0 {
-            Self::NAN
-        } else if self >= (1u32 << f32::MANTISSA_DIGITS / 2) as f32 {
-            self.ln() + consts::LN_2
-        } else {
-            (self + ((self - 1.0).sqrt() * (self + 1.0).sqrt())).ln()
-        }
+        cmath::acoshf(self)
     }
 
     /// Inverse hyperbolic tangent function.

library/std/src/num/f64.rs

@@ -1091,13 +1091,7 @@ impl f64 {
     #[stable(feature = "rust1", since = "1.0.0")]
     #[inline]
     pub fn asinh(self) -> f64 {
-        let ax = self.abs();
-        if ax >= (1u64 << f64::MANTISSA_DIGITS / 2) as f64 {
-            return (ax.ln() + consts::LN_2).copysign(self);
-        }
-
-        let ix = 1.0 / ax;
-        (ax + (ax / (Self::hypot(1.0, ix) + ix))).ln_1p().copysign(self)
+        cmath::asinh(self)
     }
 
     /// Inverse hyperbolic cosine function.
@@ -1123,13 +1117,7 @@ impl f64 {
     #[stable(feature = "rust1", since = "1.0.0")]
     #[inline]
     pub fn acosh(self) -> f64 {
-        if self < 1.0 {
-            Self::NAN
-        } else if self >= (1u64 << f64::MANTISSA_DIGITS / 2) as f64 {
-            self.ln() + consts::LN_2
-        } else {
-            (self + ((self - 1.0).sqrt() * (self + 1.0).sqrt())).ln()
-        }
+        cmath::acosh(self)
     }
 
     /// Inverse hyperbolic tangent function.

library/std/src/sys/cmath.rs

@@ -4,7 +4,9 @@
 // or by `compiler-builtins` on unsupported platforms.
 unsafe extern "C" {
     pub safe fn acos(n: f64) -> f64;
+    pub safe fn acosh(n: f64) -> f64;
     pub safe fn asin(n: f64) -> f64;
+    pub safe fn asinh(n: f64) -> f64;
     pub safe fn atan(n: f64) -> f64;
     pub safe fn atan2(a: f64, b: f64) -> f64;
     pub safe fn cosh(n: f64) -> f64;
@@ -57,6 +59,16 @@ cfg_select! {
             f64::acos(n as f64) as f32
         }
 
+        #[inline]
+        pub fn acoshf(n: f32) -> f32 {
+            f64::acosh(n as f64) as f32
+        }
+
+        #[inline]
+        pub fn asinhf(n: f32) -> f32 {
+            f64::asinh(n as f64) as f32
+        }
+
         #[inline]
         pub fn asinf(n: f32) -> f32 {
             f64::asin(n as f64) as f32
@@ -95,7 +107,9 @@ cfg_select! {
     _ => {
         unsafe extern "C" {
             pub safe fn acosf(n: f32) -> f32;
+            pub safe fn acoshf(n: f32) -> f32;
             pub safe fn asinf(n: f32) -> f32;
+            pub safe fn asinhf(n: f32) -> f32;
             pub safe fn atan2f(a: f32, b: f32) -> f32;
             pub safe fn atanf(n: f32) -> f32;
             pub safe fn coshf(n: f32) -> f32;