
(FPCore (x) :precision binary32 (asinh x))
float code(float x) {
return asinhf(x);
}
function code(x) return asinh(x) end
function tmp = code(x) tmp = asinh(x); end
\begin{array}{l}
\\
\sinh^{-1} x
\end{array}
Sampling outcomes in binary32 precision:
Herbie found 13 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (x) :precision binary32 (copysign (log (+ (fabs x) (sqrt (+ (* x x) 1.0)))) x))
float code(float x) {
return copysignf(logf((fabsf(x) + sqrtf(((x * x) + 1.0f)))), x);
}
function code(x) return copysign(log(Float32(abs(x) + sqrt(Float32(Float32(x * x) + Float32(1.0))))), x) end
function tmp = code(x) tmp = sign(x) * abs(log((abs(x) + sqrt(((x * x) + single(1.0)))))); end
\begin{array}{l}
\\
\mathsf{copysign}\left(\log \left(\left|x\right| + \sqrt{x \cdot x + 1}\right), x\right)
\end{array}
(FPCore (x)
:precision binary32
(if (<= x -10.0)
(copysign (log (/ -0.5 x)) x)
(if (<= x 0.05000000074505806)
(copysign (+ x (* (pow x 3.0) -0.16666666666666666)) x)
(copysign (log (+ x (hypot 1.0 x))) x))))
float code(float x) {
float tmp;
if (x <= -10.0f) {
tmp = copysignf(logf((-0.5f / x)), x);
} else if (x <= 0.05000000074505806f) {
tmp = copysignf((x + (powf(x, 3.0f) * -0.16666666666666666f)), x);
} else {
tmp = copysignf(logf((x + hypotf(1.0f, x))), x);
}
return tmp;
}
function code(x) tmp = Float32(0.0) if (x <= Float32(-10.0)) tmp = copysign(log(Float32(Float32(-0.5) / x)), x); elseif (x <= Float32(0.05000000074505806)) tmp = copysign(Float32(x + Float32((x ^ Float32(3.0)) * Float32(-0.16666666666666666))), x); else tmp = copysign(log(Float32(x + hypot(Float32(1.0), x))), x); end return tmp end
function tmp_2 = code(x) tmp = single(0.0); if (x <= single(-10.0)) tmp = sign(x) * abs(log((single(-0.5) / x))); elseif (x <= single(0.05000000074505806)) tmp = sign(x) * abs((x + ((x ^ single(3.0)) * single(-0.16666666666666666)))); else tmp = sign(x) * abs(log((x + hypot(single(1.0), x)))); end tmp_2 = tmp; end
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq -10:\\
\;\;\;\;\mathsf{copysign}\left(\log \left(\frac{-0.5}{x}\right), x\right)\\
\mathbf{elif}\;x \leq 0.05000000074505806:\\
\;\;\;\;\mathsf{copysign}\left(x + {x}^{3} \cdot -0.16666666666666666, x\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{copysign}\left(\log \left(x + \mathsf{hypot}\left(1, x\right)\right), x\right)\\
\end{array}
\end{array}
if x < -10Initial program 48.3%
Taylor expanded in x around -inf 99.9%
sub-neg99.9%
neg-mul-199.9%
unsub-neg99.9%
associate-*r/99.9%
metadata-eval99.9%
distribute-neg-frac99.9%
metadata-eval99.9%
Simplified99.9%
Taylor expanded in x around 0 -0.0%
neg-mul-1-0.0%
unsub-neg-0.0%
log-div99.0%
Simplified99.0%
if -10 < x < 0.0500000007Initial program 19.3%
Applied egg-rr19.6%
+-commutative19.6%
fma-udef19.6%
unpow219.6%
+-commutative19.6%
associate-+l+19.6%
log1p-def97.3%
unpow297.3%
distribute-lft-out--97.3%
Simplified97.3%
Taylor expanded in x around 0 99.3%
*-commutative99.3%
Simplified99.3%
if 0.0500000007 < x Initial program 61.9%
*-un-lft-identity61.9%
*-commutative61.9%
log-prod61.9%
+-commutative61.9%
hypot-1-def98.1%
add-sqr-sqrt98.0%
fabs-sqr98.0%
add-sqr-sqrt98.1%
metadata-eval98.1%
Applied egg-rr98.1%
+-rgt-identity98.1%
Simplified98.1%
Final simplification98.9%
(FPCore (x) :precision binary32 (if (<= (copysign (log (+ (fabs x) (sqrt (+ (* x x) 1.0)))) x) -2.0) (copysign (log (+ (fabs x) (hypot 1.0 x))) x) (copysign (log1p (+ x (+ (hypot 1.0 x) -1.0))) x)))
float code(float x) {
float tmp;
if (copysignf(logf((fabsf(x) + sqrtf(((x * x) + 1.0f)))), x) <= -2.0f) {
tmp = copysignf(logf((fabsf(x) + hypotf(1.0f, x))), x);
} else {
tmp = copysignf(log1pf((x + (hypotf(1.0f, x) + -1.0f))), x);
}
return tmp;
}
function code(x) tmp = Float32(0.0) if (copysign(log(Float32(abs(x) + sqrt(Float32(Float32(x * x) + Float32(1.0))))), x) <= Float32(-2.0)) tmp = copysign(log(Float32(abs(x) + hypot(Float32(1.0), x))), x); else tmp = copysign(log1p(Float32(x + Float32(hypot(Float32(1.0), x) + Float32(-1.0)))), x); end return tmp end
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\mathsf{copysign}\left(\log \left(\left|x\right| + \sqrt{x \cdot x + 1}\right), x\right) \leq -2:\\
\;\;\;\;\mathsf{copysign}\left(\log \left(\left|x\right| + \mathsf{hypot}\left(1, x\right)\right), x\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{copysign}\left(\mathsf{log1p}\left(x + \left(\mathsf{hypot}\left(1, x\right) + -1\right)\right), x\right)\\
\end{array}
\end{array}
if (copysign.f32 (log.f32 (+.f32 (fabs.f32 x) (sqrt.f32 (+.f32 (*.f32 x x) 1)))) x) < -2Initial program 48.3%
+-commutative48.3%
hypot-1-def100.0%
Simplified100.0%
if -2 < (copysign.f32 (log.f32 (+.f32 (fabs.f32 x) (sqrt.f32 (+.f32 (*.f32 x x) 1)))) x) Initial program 33.2%
*-un-lft-identity33.2%
add-log-exp18.4%
exp-prod18.2%
log-pow32.9%
+-commutative32.9%
hypot-1-def44.8%
add-sqr-sqrt38.7%
fabs-sqr38.7%
add-sqr-sqrt44.9%
Applied egg-rr44.9%
log1p-expm1-u44.9%
expm1-udef44.9%
rem-log-exp45.3%
add-exp-log45.3%
*-rgt-identity45.3%
Applied egg-rr45.3%
associate--l+98.0%
Simplified98.0%
Final simplification98.5%
(FPCore (x) :precision binary32 (if (<= x -10.0) (copysign (log (+ (- (fabs x) x) (/ -0.5 x))) x) (copysign (log1p (+ x (+ (hypot 1.0 x) -1.0))) x)))
float code(float x) {
float tmp;
if (x <= -10.0f) {
tmp = copysignf(logf(((fabsf(x) - x) + (-0.5f / x))), x);
} else {
tmp = copysignf(log1pf((x + (hypotf(1.0f, x) + -1.0f))), x);
}
return tmp;
}
function code(x) tmp = Float32(0.0) if (x <= Float32(-10.0)) tmp = copysign(log(Float32(Float32(abs(x) - x) + Float32(Float32(-0.5) / x))), x); else tmp = copysign(log1p(Float32(x + Float32(hypot(Float32(1.0), x) + Float32(-1.0)))), x); end return tmp end
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq -10:\\
\;\;\;\;\mathsf{copysign}\left(\log \left(\left(\left|x\right| - x\right) + \frac{-0.5}{x}\right), x\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{copysign}\left(\mathsf{log1p}\left(x + \left(\mathsf{hypot}\left(1, x\right) + -1\right)\right), x\right)\\
\end{array}
\end{array}
if x < -10Initial program 48.3%
Taylor expanded in x around -inf 99.9%
sub-neg99.9%
neg-mul-199.9%
unsub-neg99.9%
associate-*r/99.9%
metadata-eval99.9%
distribute-neg-frac99.9%
metadata-eval99.9%
Simplified99.9%
if -10 < x Initial program 33.2%
*-un-lft-identity33.2%
add-log-exp18.4%
exp-prod18.2%
log-pow32.9%
+-commutative32.9%
hypot-1-def44.8%
add-sqr-sqrt38.7%
fabs-sqr38.7%
add-sqr-sqrt44.9%
Applied egg-rr44.9%
log1p-expm1-u44.9%
expm1-udef44.9%
rem-log-exp45.3%
add-exp-log45.3%
*-rgt-identity45.3%
Applied egg-rr45.3%
associate--l+98.0%
Simplified98.0%
Final simplification98.5%
(FPCore (x) :precision binary32 (if (<= x -10.0) (copysign (log (- (fabs x) (+ x (/ 0.5 x)))) x) (copysign (log1p (+ x (+ (hypot 1.0 x) -1.0))) x)))
float code(float x) {
float tmp;
if (x <= -10.0f) {
tmp = copysignf(logf((fabsf(x) - (x + (0.5f / x)))), x);
} else {
tmp = copysignf(log1pf((x + (hypotf(1.0f, x) + -1.0f))), x);
}
return tmp;
}
function code(x) tmp = Float32(0.0) if (x <= Float32(-10.0)) tmp = copysign(log(Float32(abs(x) - Float32(x + Float32(Float32(0.5) / x)))), x); else tmp = copysign(log1p(Float32(x + Float32(hypot(Float32(1.0), x) + Float32(-1.0)))), x); end return tmp end
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq -10:\\
\;\;\;\;\mathsf{copysign}\left(\log \left(\left|x\right| - \left(x + \frac{0.5}{x}\right)\right), x\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{copysign}\left(\mathsf{log1p}\left(x + \left(\mathsf{hypot}\left(1, x\right) + -1\right)\right), x\right)\\
\end{array}
\end{array}
if x < -10Initial program 48.3%
Taylor expanded in x around -inf 99.9%
neg-mul-199.9%
associate-*r/99.9%
metadata-eval99.9%
Simplified99.9%
if -10 < x Initial program 33.2%
*-un-lft-identity33.2%
add-log-exp18.4%
exp-prod18.2%
log-pow32.9%
+-commutative32.9%
hypot-1-def44.8%
add-sqr-sqrt38.7%
fabs-sqr38.7%
add-sqr-sqrt44.9%
Applied egg-rr44.9%
log1p-expm1-u44.9%
expm1-udef44.9%
rem-log-exp45.3%
add-exp-log45.3%
*-rgt-identity45.3%
Applied egg-rr45.3%
associate--l+98.0%
Simplified98.0%
Final simplification98.5%
(FPCore (x) :precision binary32 (if (<= x -400.0) (copysign (log (/ -0.5 x)) x) (copysign (log1p (+ x (+ (hypot 1.0 x) -1.0))) x)))
float code(float x) {
float tmp;
if (x <= -400.0f) {
tmp = copysignf(logf((-0.5f / x)), x);
} else {
tmp = copysignf(log1pf((x + (hypotf(1.0f, x) + -1.0f))), x);
}
return tmp;
}
function code(x) tmp = Float32(0.0) if (x <= Float32(-400.0)) tmp = copysign(log(Float32(Float32(-0.5) / x)), x); else tmp = copysign(log1p(Float32(x + Float32(hypot(Float32(1.0), x) + Float32(-1.0)))), x); end return tmp end
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq -400:\\
\;\;\;\;\mathsf{copysign}\left(\log \left(\frac{-0.5}{x}\right), x\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{copysign}\left(\mathsf{log1p}\left(x + \left(\mathsf{hypot}\left(1, x\right) + -1\right)\right), x\right)\\
\end{array}
\end{array}
if x < -400Initial program 46.8%
Taylor expanded in x around -inf 100.0%
sub-neg100.0%
neg-mul-1100.0%
unsub-neg100.0%
associate-*r/100.0%
metadata-eval100.0%
distribute-neg-frac100.0%
metadata-eval100.0%
Simplified100.0%
Taylor expanded in x around 0 -0.0%
neg-mul-1-0.0%
unsub-neg-0.0%
log-div99.9%
Simplified99.9%
if -400 < x Initial program 33.9%
*-un-lft-identity33.9%
add-log-exp18.8%
exp-prod18.6%
log-pow33.6%
+-commutative33.6%
hypot-1-def45.4%
add-sqr-sqrt38.3%
fabs-sqr38.3%
add-sqr-sqrt45.3%
Applied egg-rr45.3%
log1p-expm1-u45.3%
expm1-udef45.3%
rem-log-exp45.6%
add-exp-log45.7%
*-rgt-identity45.7%
Applied egg-rr45.7%
associate--l+97.8%
Simplified97.8%
Final simplification98.4%
(FPCore (x)
:precision binary32
(if (<= x -10.0)
(copysign (log (/ -0.5 x)) x)
(if (<= x 1.0)
(copysign (+ x (* (pow x 3.0) -0.16666666666666666)) x)
(copysign (log1p (+ x (+ x (+ (/ 0.5 x) -1.0)))) x))))
float code(float x) {
float tmp;
if (x <= -10.0f) {
tmp = copysignf(logf((-0.5f / x)), x);
} else if (x <= 1.0f) {
tmp = copysignf((x + (powf(x, 3.0f) * -0.16666666666666666f)), x);
} else {
tmp = copysignf(log1pf((x + (x + ((0.5f / x) + -1.0f)))), x);
}
return tmp;
}
function code(x) tmp = Float32(0.0) if (x <= Float32(-10.0)) tmp = copysign(log(Float32(Float32(-0.5) / x)), x); elseif (x <= Float32(1.0)) tmp = copysign(Float32(x + Float32((x ^ Float32(3.0)) * Float32(-0.16666666666666666))), x); else tmp = copysign(log1p(Float32(x + Float32(x + Float32(Float32(Float32(0.5) / x) + Float32(-1.0))))), x); end return tmp end
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq -10:\\
\;\;\;\;\mathsf{copysign}\left(\log \left(\frac{-0.5}{x}\right), x\right)\\
\mathbf{elif}\;x \leq 1:\\
\;\;\;\;\mathsf{copysign}\left(x + {x}^{3} \cdot -0.16666666666666666, x\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{copysign}\left(\mathsf{log1p}\left(x + \left(x + \left(\frac{0.5}{x} + -1\right)\right)\right), x\right)\\
\end{array}
\end{array}
if x < -10Initial program 48.3%
Taylor expanded in x around -inf 99.9%
sub-neg99.9%
neg-mul-199.9%
unsub-neg99.9%
associate-*r/99.9%
metadata-eval99.9%
distribute-neg-frac99.9%
metadata-eval99.9%
Simplified99.9%
Taylor expanded in x around 0 -0.0%
neg-mul-1-0.0%
unsub-neg-0.0%
log-div99.0%
Simplified99.0%
if -10 < x < 1Initial program 21.6%
Applied egg-rr21.8%
+-commutative21.8%
fma-udef21.7%
unpow221.7%
+-commutative21.7%
associate-+l+21.8%
log1p-def97.2%
unpow297.2%
distribute-lft-out--97.2%
Simplified97.2%
Taylor expanded in x around 0 98.3%
*-commutative98.3%
Simplified98.3%
if 1 < x Initial program 59.6%
*-un-lft-identity59.6%
add-log-exp9.7%
exp-prod9.7%
log-pow59.6%
+-commutative59.6%
hypot-1-def98.3%
add-sqr-sqrt98.3%
fabs-sqr98.3%
add-sqr-sqrt98.3%
Applied egg-rr98.3%
log1p-expm1-u98.3%
expm1-udef98.3%
rem-log-exp98.3%
add-exp-log98.3%
*-rgt-identity98.3%
Applied egg-rr98.3%
associate--l+98.3%
Simplified98.3%
Taylor expanded in x around inf 97.1%
associate--l+97.1%
associate-*r/97.1%
metadata-eval97.1%
Simplified97.1%
Final simplification98.2%
(FPCore (x)
:precision binary32
(if (<= x -10.0)
(copysign (log (/ -0.5 x)) x)
(if (<= x 1.0)
(copysign (+ x (* (pow x 3.0) -0.16666666666666666)) x)
(copysign (log (* x 2.0)) x))))
float code(float x) {
float tmp;
if (x <= -10.0f) {
tmp = copysignf(logf((-0.5f / x)), x);
} else if (x <= 1.0f) {
tmp = copysignf((x + (powf(x, 3.0f) * -0.16666666666666666f)), x);
} else {
tmp = copysignf(logf((x * 2.0f)), x);
}
return tmp;
}
function code(x) tmp = Float32(0.0) if (x <= Float32(-10.0)) tmp = copysign(log(Float32(Float32(-0.5) / x)), x); elseif (x <= Float32(1.0)) tmp = copysign(Float32(x + Float32((x ^ Float32(3.0)) * Float32(-0.16666666666666666))), x); else tmp = copysign(log(Float32(x * Float32(2.0))), x); end return tmp end
function tmp_2 = code(x) tmp = single(0.0); if (x <= single(-10.0)) tmp = sign(x) * abs(log((single(-0.5) / x))); elseif (x <= single(1.0)) tmp = sign(x) * abs((x + ((x ^ single(3.0)) * single(-0.16666666666666666)))); else tmp = sign(x) * abs(log((x * single(2.0)))); end tmp_2 = tmp; end
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq -10:\\
\;\;\;\;\mathsf{copysign}\left(\log \left(\frac{-0.5}{x}\right), x\right)\\
\mathbf{elif}\;x \leq 1:\\
\;\;\;\;\mathsf{copysign}\left(x + {x}^{3} \cdot -0.16666666666666666, x\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{copysign}\left(\log \left(x \cdot 2\right), x\right)\\
\end{array}
\end{array}
if x < -10Initial program 48.3%
Taylor expanded in x around -inf 99.9%
sub-neg99.9%
neg-mul-199.9%
unsub-neg99.9%
associate-*r/99.9%
metadata-eval99.9%
distribute-neg-frac99.9%
metadata-eval99.9%
Simplified99.9%
Taylor expanded in x around 0 -0.0%
neg-mul-1-0.0%
unsub-neg-0.0%
log-div99.0%
Simplified99.0%
if -10 < x < 1Initial program 21.6%
Applied egg-rr21.8%
+-commutative21.8%
fma-udef21.7%
unpow221.7%
+-commutative21.7%
associate-+l+21.8%
log1p-def97.2%
unpow297.2%
distribute-lft-out--97.2%
Simplified97.2%
Taylor expanded in x around 0 98.3%
*-commutative98.3%
Simplified98.3%
if 1 < x Initial program 59.6%
Applied egg-rr36.5%
+-commutative36.5%
fma-udef36.5%
unpow236.5%
+-commutative36.5%
associate-+l+36.5%
log1p-def36.5%
unpow236.5%
distribute-lft-out--36.5%
Simplified36.5%
Taylor expanded in x around inf 94.7%
associate--l+95.0%
distribute-rgt-out--96.4%
metadata-eval96.4%
*-commutative96.4%
+-commutative96.4%
mul-1-neg96.4%
log-rec96.5%
remove-double-neg96.5%
log-prod95.5%
Simplified95.5%
Final simplification97.8%
(FPCore (x) :precision binary32 (if (<= x -10.0) (copysign (log (- x)) x) (if (<= x 1.0) (copysign x x) (copysign (log (* x 2.0)) x))))
float code(float x) {
float tmp;
if (x <= -10.0f) {
tmp = copysignf(logf(-x), x);
} else if (x <= 1.0f) {
tmp = copysignf(x, x);
} else {
tmp = copysignf(logf((x * 2.0f)), x);
}
return tmp;
}
function code(x) tmp = Float32(0.0) if (x <= Float32(-10.0)) tmp = copysign(log(Float32(-x)), x); elseif (x <= Float32(1.0)) tmp = copysign(x, x); else tmp = copysign(log(Float32(x * Float32(2.0))), x); end return tmp end
function tmp_2 = code(x) tmp = single(0.0); if (x <= single(-10.0)) tmp = sign(x) * abs(log(-x)); elseif (x <= single(1.0)) tmp = sign(x) * abs(x); else tmp = sign(x) * abs(log((x * single(2.0)))); end tmp_2 = tmp; end
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq -10:\\
\;\;\;\;\mathsf{copysign}\left(\log \left(-x\right), x\right)\\
\mathbf{elif}\;x \leq 1:\\
\;\;\;\;\mathsf{copysign}\left(x, x\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{copysign}\left(\log \left(x \cdot 2\right), x\right)\\
\end{array}
\end{array}
if x < -10Initial program 48.3%
Taylor expanded in x around -inf 99.9%
sub-neg99.9%
neg-mul-199.9%
unsub-neg99.9%
associate-*r/99.9%
metadata-eval99.9%
distribute-neg-frac99.9%
metadata-eval99.9%
Simplified99.9%
Taylor expanded in x around inf -0.0%
mul-1-neg-0.0%
log-rec-0.0%
remove-double-neg-0.0%
log-prod44.7%
neg-mul-144.7%
Simplified44.7%
if -10 < x < 1Initial program 21.6%
Applied egg-rr21.8%
+-commutative21.8%
fma-udef21.7%
unpow221.7%
+-commutative21.7%
associate-+l+21.8%
log1p-def97.2%
unpow297.2%
distribute-lft-out--97.2%
Simplified97.2%
Taylor expanded in x around 0 96.6%
if 1 < x Initial program 59.6%
Applied egg-rr36.5%
+-commutative36.5%
fma-udef36.5%
unpow236.5%
+-commutative36.5%
associate-+l+36.5%
log1p-def36.5%
unpow236.5%
distribute-lft-out--36.5%
Simplified36.5%
Taylor expanded in x around inf 94.7%
associate--l+95.0%
distribute-rgt-out--96.4%
metadata-eval96.4%
*-commutative96.4%
+-commutative96.4%
mul-1-neg96.4%
log-rec96.5%
remove-double-neg96.5%
log-prod95.5%
Simplified95.5%
Final simplification82.7%
(FPCore (x) :precision binary32 (if (<= x -10.0) (copysign (log (/ -0.5 x)) x) (if (<= x 1.0) (copysign x x) (copysign (log (* x 2.0)) x))))
float code(float x) {
float tmp;
if (x <= -10.0f) {
tmp = copysignf(logf((-0.5f / x)), x);
} else if (x <= 1.0f) {
tmp = copysignf(x, x);
} else {
tmp = copysignf(logf((x * 2.0f)), x);
}
return tmp;
}
function code(x) tmp = Float32(0.0) if (x <= Float32(-10.0)) tmp = copysign(log(Float32(Float32(-0.5) / x)), x); elseif (x <= Float32(1.0)) tmp = copysign(x, x); else tmp = copysign(log(Float32(x * Float32(2.0))), x); end return tmp end
function tmp_2 = code(x) tmp = single(0.0); if (x <= single(-10.0)) tmp = sign(x) * abs(log((single(-0.5) / x))); elseif (x <= single(1.0)) tmp = sign(x) * abs(x); else tmp = sign(x) * abs(log((x * single(2.0)))); end tmp_2 = tmp; end
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq -10:\\
\;\;\;\;\mathsf{copysign}\left(\log \left(\frac{-0.5}{x}\right), x\right)\\
\mathbf{elif}\;x \leq 1:\\
\;\;\;\;\mathsf{copysign}\left(x, x\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{copysign}\left(\log \left(x \cdot 2\right), x\right)\\
\end{array}
\end{array}
if x < -10Initial program 48.3%
Taylor expanded in x around -inf 99.9%
sub-neg99.9%
neg-mul-199.9%
unsub-neg99.9%
associate-*r/99.9%
metadata-eval99.9%
distribute-neg-frac99.9%
metadata-eval99.9%
Simplified99.9%
Taylor expanded in x around 0 -0.0%
neg-mul-1-0.0%
unsub-neg-0.0%
log-div99.0%
Simplified99.0%
if -10 < x < 1Initial program 21.6%
Applied egg-rr21.8%
+-commutative21.8%
fma-udef21.7%
unpow221.7%
+-commutative21.7%
associate-+l+21.8%
log1p-def97.2%
unpow297.2%
distribute-lft-out--97.2%
Simplified97.2%
Taylor expanded in x around 0 96.6%
if 1 < x Initial program 59.6%
Applied egg-rr36.5%
+-commutative36.5%
fma-udef36.5%
unpow236.5%
+-commutative36.5%
associate-+l+36.5%
log1p-def36.5%
unpow236.5%
distribute-lft-out--36.5%
Simplified36.5%
Taylor expanded in x around inf 94.7%
associate--l+95.0%
distribute-rgt-out--96.4%
metadata-eval96.4%
*-commutative96.4%
+-commutative96.4%
mul-1-neg96.4%
log-rec96.5%
remove-double-neg96.5%
log-prod95.5%
Simplified95.5%
Final simplification97.0%
(FPCore (x) :precision binary32 (if (<= x -10.0) (copysign (log (- x)) x) (copysign (log1p x) x)))
float code(float x) {
float tmp;
if (x <= -10.0f) {
tmp = copysignf(logf(-x), x);
} else {
tmp = copysignf(log1pf(x), x);
}
return tmp;
}
function code(x) tmp = Float32(0.0) if (x <= Float32(-10.0)) tmp = copysign(log(Float32(-x)), x); else tmp = copysign(log1p(x), x); end return tmp end
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq -10:\\
\;\;\;\;\mathsf{copysign}\left(\log \left(-x\right), x\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{copysign}\left(\mathsf{log1p}\left(x\right), x\right)\\
\end{array}
\end{array}
if x < -10Initial program 48.3%
Taylor expanded in x around -inf 99.9%
sub-neg99.9%
neg-mul-199.9%
unsub-neg99.9%
associate-*r/99.9%
metadata-eval99.9%
distribute-neg-frac99.9%
metadata-eval99.9%
Simplified99.9%
Taylor expanded in x around inf -0.0%
mul-1-neg-0.0%
log-rec-0.0%
remove-double-neg-0.0%
log-prod44.7%
neg-mul-144.7%
Simplified44.7%
if -10 < x Initial program 33.2%
Taylor expanded in x around 0 25.6%
log1p-def78.4%
Simplified78.4%
expm1-log1p-u78.4%
expm1-udef25.6%
add-sqr-sqrt20.5%
fabs-sqr20.5%
add-sqr-sqrt25.7%
Applied egg-rr25.7%
expm1-def78.4%
expm1-log1p78.4%
Simplified78.4%
Final simplification69.6%
(FPCore (x) :precision binary32 (if (<= x 3.200000047683716) (copysign x x) (copysign (log x) x)))
float code(float x) {
float tmp;
if (x <= 3.200000047683716f) {
tmp = copysignf(x, x);
} else {
tmp = copysignf(logf(x), x);
}
return tmp;
}
function code(x) tmp = Float32(0.0) if (x <= Float32(3.200000047683716)) tmp = copysign(x, x); else tmp = copysign(log(x), x); end return tmp end
function tmp_2 = code(x) tmp = single(0.0); if (x <= single(3.200000047683716)) tmp = sign(x) * abs(x); else tmp = sign(x) * abs(log(x)); end tmp_2 = tmp; end
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq 3.200000047683716:\\
\;\;\;\;\mathsf{copysign}\left(x, x\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{copysign}\left(\log x, x\right)\\
\end{array}
\end{array}
if x < 3.20000005Initial program 31.0%
Applied egg-rr16.8%
+-commutative16.8%
fma-udef16.8%
unpow216.8%
+-commutative16.8%
associate-+l+16.8%
log1p-def66.4%
unpow266.4%
distribute-lft-out--66.4%
Simplified66.4%
Taylor expanded in x around 0 67.3%
if 3.20000005 < x Initial program 58.9%
Taylor expanded in x around inf 44.0%
mul-1-neg44.0%
log-rec44.0%
remove-double-neg44.0%
Simplified44.0%
Final simplification62.1%
(FPCore (x) :precision binary32 (if (<= x 1.0) (copysign x x) (copysign (log1p x) x)))
float code(float x) {
float tmp;
if (x <= 1.0f) {
tmp = copysignf(x, x);
} else {
tmp = copysignf(log1pf(x), x);
}
return tmp;
}
function code(x) tmp = Float32(0.0) if (x <= Float32(1.0)) tmp = copysign(x, x); else tmp = copysign(log1p(x), x); end return tmp end
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq 1:\\
\;\;\;\;\mathsf{copysign}\left(x, x\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{copysign}\left(\mathsf{log1p}\left(x\right), x\right)\\
\end{array}
\end{array}
if x < 1Initial program 30.6%
Applied egg-rr16.4%
+-commutative16.4%
fma-udef16.3%
unpow216.3%
+-commutative16.3%
associate-+l+16.4%
log1p-def66.3%
unpow266.3%
distribute-lft-out--66.3%
Simplified66.3%
Taylor expanded in x around 0 67.5%
if 1 < x Initial program 59.6%
Taylor expanded in x around 0 43.9%
log1p-def43.9%
Simplified43.9%
expm1-log1p-u43.9%
expm1-udef43.9%
add-sqr-sqrt43.9%
fabs-sqr43.9%
add-sqr-sqrt43.9%
Applied egg-rr43.9%
expm1-def43.9%
expm1-log1p43.9%
Simplified43.9%
Final simplification62.1%
(FPCore (x) :precision binary32 (copysign x x))
float code(float x) {
return copysignf(x, x);
}
function code(x) return copysign(x, x) end
function tmp = code(x) tmp = sign(x) * abs(x); end
\begin{array}{l}
\\
\mathsf{copysign}\left(x, x\right)
\end{array}
Initial program 37.2%
Applied egg-rr20.9%
+-commutative20.9%
fma-udef20.9%
unpow220.9%
+-commutative20.9%
associate-+l+20.9%
log1p-def59.5%
unpow259.5%
distribute-lft-out--59.5%
Simplified59.5%
Taylor expanded in x around 0 54.9%
Final simplification54.9%
(FPCore (x) :precision binary32 (let* ((t_0 (/ 1.0 (fabs x)))) (copysign (log1p (+ (fabs x) (/ (fabs x) (+ (hypot 1.0 t_0) t_0)))) x)))
float code(float x) {
float t_0 = 1.0f / fabsf(x);
return copysignf(log1pf((fabsf(x) + (fabsf(x) / (hypotf(1.0f, t_0) + t_0)))), x);
}
function code(x) t_0 = Float32(Float32(1.0) / abs(x)) return copysign(log1p(Float32(abs(x) + Float32(abs(x) / Float32(hypot(Float32(1.0), t_0) + t_0)))), x) end
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \frac{1}{\left|x\right|}\\
\mathsf{copysign}\left(\mathsf{log1p}\left(\left|x\right| + \frac{\left|x\right|}{\mathsf{hypot}\left(1, t_0\right) + t_0}\right), x\right)
\end{array}
\end{array}
herbie shell --seed 2024013
(FPCore (x)
:name "Rust f32::asinh"
:precision binary32
:herbie-target
(copysign (log1p (+ (fabs x) (/ (fabs x) (+ (hypot 1.0 (/ 1.0 (fabs x))) (/ 1.0 (fabs x)))))) x)
(copysign (log (+ (fabs x) (sqrt (+ (* x x) 1.0)))) x))