| Alternative 1 | |
|---|---|
| Accuracy | 99.5% |
| Cost | 39144 |
(FPCore (x) :precision binary32 (copysign (log (+ (fabs x) (sqrt (+ (* x x) 1.0)))) x))
(FPCore (x)
:precision binary32
(let* ((t_0 (copysign (log (+ (fabs x) (sqrt (+ (* x x) 1.0)))) x)))
(if (<= t_0 -0.5)
(copysign (- (log (- (hypot 1.0 x) x))) x)
(if (<= t_0 0.11999999731779099)
(copysign
(+
(* -0.16666666666666666 (pow x 3.0))
(+ (* 0.075 (pow x 5.0)) (+ x (* -0.044642857142857144 (pow x 7.0)))))
x)
(copysign (log (+ (fabs x) (hypot 1.0 x))) x)))))float code(float x) {
return copysignf(logf((fabsf(x) + sqrtf(((x * x) + 1.0f)))), x);
}
float code(float x) {
float t_0 = copysignf(logf((fabsf(x) + sqrtf(((x * x) + 1.0f)))), x);
float tmp;
if (t_0 <= -0.5f) {
tmp = copysignf(-logf((hypotf(1.0f, x) - x)), x);
} else if (t_0 <= 0.11999999731779099f) {
tmp = copysignf(((-0.16666666666666666f * powf(x, 3.0f)) + ((0.075f * powf(x, 5.0f)) + (x + (-0.044642857142857144f * powf(x, 7.0f))))), x);
} else {
tmp = copysignf(logf((fabsf(x) + hypotf(1.0f, x))), x);
}
return tmp;
}
function code(x) return copysign(log(Float32(abs(x) + sqrt(Float32(Float32(x * x) + Float32(1.0))))), x) end
function code(x) t_0 = copysign(log(Float32(abs(x) + sqrt(Float32(Float32(x * x) + Float32(1.0))))), x) tmp = Float32(0.0) if (t_0 <= Float32(-0.5)) tmp = copysign(Float32(-log(Float32(hypot(Float32(1.0), x) - x))), x); elseif (t_0 <= Float32(0.11999999731779099)) tmp = copysign(Float32(Float32(Float32(-0.16666666666666666) * (x ^ Float32(3.0))) + Float32(Float32(Float32(0.075) * (x ^ Float32(5.0))) + Float32(x + Float32(Float32(-0.044642857142857144) * (x ^ Float32(7.0)))))), x); else tmp = copysign(log(Float32(abs(x) + hypot(Float32(1.0), x))), x); end return tmp end
function tmp = code(x) tmp = sign(x) * abs(log((abs(x) + sqrt(((x * x) + single(1.0)))))); end
function tmp_2 = code(x) t_0 = sign(x) * abs(log((abs(x) + sqrt(((x * x) + single(1.0)))))); tmp = single(0.0); if (t_0 <= single(-0.5)) tmp = sign(x) * abs(-log((hypot(single(1.0), x) - x))); elseif (t_0 <= single(0.11999999731779099)) tmp = sign(x) * abs(((single(-0.16666666666666666) * (x ^ single(3.0))) + ((single(0.075) * (x ^ single(5.0))) + (x + (single(-0.044642857142857144) * (x ^ single(7.0))))))); else tmp = sign(x) * abs(log((abs(x) + hypot(single(1.0), x)))); end tmp_2 = tmp; end
\mathsf{copysign}\left(\log \left(\left|x\right| + \sqrt{x \cdot x + 1}\right), x\right)
\begin{array}{l}
t_0 := \mathsf{copysign}\left(\log \left(\left|x\right| + \sqrt{x \cdot x + 1}\right), x\right)\\
\mathbf{if}\;t_0 \leq -0.5:\\
\;\;\;\;\mathsf{copysign}\left(-\log \left(\mathsf{hypot}\left(1, x\right) - x\right), x\right)\\
\mathbf{elif}\;t_0 \leq 0.11999999731779099:\\
\;\;\;\;\mathsf{copysign}\left(-0.16666666666666666 \cdot {x}^{3} + \left(0.075 \cdot {x}^{5} + \left(x + -0.044642857142857144 \cdot {x}^{7}\right)\right), x\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{copysign}\left(\log \left(\left|x\right| + \mathsf{hypot}\left(1, x\right)\right), x\right)\\
\end{array}
| Original | 37.8% |
|---|---|
| Target | 99.6% |
| Herbie | 99.5% |
if (copysign.f32 (log.f32 (+.f32 (fabs.f32 x) (sqrt.f32 (+.f32 (*.f32 x x) 1)))) x) < -0.5Initial program 50.9%
Applied egg-rr17.3%
[Start]50.9 | \[ \mathsf{copysign}\left(\log \left(\left|x\right| + \sqrt{x \cdot x + 1}\right), x\right)
\] |
|---|---|
flip-+ [=>]9.4 | \[ \mathsf{copysign}\left(\log \color{blue}{\left(\frac{\left|x\right| \cdot \left|x\right| - \sqrt{x \cdot x + 1} \cdot \sqrt{x \cdot x + 1}}{\left|x\right| - \sqrt{x \cdot x + 1}}\right)}, x\right)
\] |
frac-2neg [=>]9.4 | \[ \mathsf{copysign}\left(\log \color{blue}{\left(\frac{-\left(\left|x\right| \cdot \left|x\right| - \sqrt{x \cdot x + 1} \cdot \sqrt{x \cdot x + 1}\right)}{-\left(\left|x\right| - \sqrt{x \cdot x + 1}\right)}\right)}, x\right)
\] |
log-div [=>]9.5 | \[ \mathsf{copysign}\left(\color{blue}{\log \left(-\left(\left|x\right| \cdot \left|x\right| - \sqrt{x \cdot x + 1} \cdot \sqrt{x \cdot x + 1}\right)\right) - \log \left(-\left(\left|x\right| - \sqrt{x \cdot x + 1}\right)\right)}, x\right)
\] |
sqr-abs [=>]9.5 | \[ \mathsf{copysign}\left(\log \left(-\left(\color{blue}{x \cdot x} - \sqrt{x \cdot x + 1} \cdot \sqrt{x \cdot x + 1}\right)\right) - \log \left(-\left(\left|x\right| - \sqrt{x \cdot x + 1}\right)\right), x\right)
\] |
add-sqr-sqrt [<=]9.4 | \[ \mathsf{copysign}\left(\log \left(-\left(x \cdot x - \color{blue}{\left(x \cdot x + 1\right)}\right)\right) - \log \left(-\left(\left|x\right| - \sqrt{x \cdot x + 1}\right)\right), x\right)
\] |
fma-def [=>]9.4 | \[ \mathsf{copysign}\left(\log \left(-\left(x \cdot x - \color{blue}{\mathsf{fma}\left(x, x, 1\right)}\right)\right) - \log \left(-\left(\left|x\right| - \sqrt{x \cdot x + 1}\right)\right), x\right)
\] |
add-sqr-sqrt [=>]-0.0 | \[ \mathsf{copysign}\left(\log \left(-\left(x \cdot x - \mathsf{fma}\left(x, x, 1\right)\right)\right) - \log \left(-\left(\left|\color{blue}{\sqrt{x} \cdot \sqrt{x}}\right| - \sqrt{x \cdot x + 1}\right)\right), x\right)
\] |
fabs-sqr [=>]-0.0 | \[ \mathsf{copysign}\left(\log \left(-\left(x \cdot x - \mathsf{fma}\left(x, x, 1\right)\right)\right) - \log \left(-\left(\color{blue}{\sqrt{x} \cdot \sqrt{x}} - \sqrt{x \cdot x + 1}\right)\right), x\right)
\] |
add-sqr-sqrt [<=]17.3 | \[ \mathsf{copysign}\left(\log \left(-\left(x \cdot x - \mathsf{fma}\left(x, x, 1\right)\right)\right) - \log \left(-\left(\color{blue}{x} - \sqrt{x \cdot x + 1}\right)\right), x\right)
\] |
+-commutative [=>]17.3 | \[ \mathsf{copysign}\left(\log \left(-\left(x \cdot x - \mathsf{fma}\left(x, x, 1\right)\right)\right) - \log \left(-\left(x - \sqrt{\color{blue}{1 + x \cdot x}}\right)\right), x\right)
\] |
Simplified98.4%
[Start]17.3 | \[ \mathsf{copysign}\left(\log \left(-\left(x \cdot x - \mathsf{fma}\left(x, x, 1\right)\right)\right) - \log \left(-\left(x - \mathsf{hypot}\left(1, x\right)\right)\right), x\right)
\] |
|---|---|
neg-mul-1 [=>]17.3 | \[ \mathsf{copysign}\left(\log \color{blue}{\left(-1 \cdot \left(x \cdot x - \mathsf{fma}\left(x, x, 1\right)\right)\right)} - \log \left(-\left(x - \mathsf{hypot}\left(1, x\right)\right)\right), x\right)
\] |
fma-udef [=>]17.3 | \[ \mathsf{copysign}\left(\log \left(-1 \cdot \left(x \cdot x - \color{blue}{\left(x \cdot x + 1\right)}\right)\right) - \log \left(-\left(x - \mathsf{hypot}\left(1, x\right)\right)\right), x\right)
\] |
associate--r+ [=>]47.5 | \[ \mathsf{copysign}\left(\log \left(-1 \cdot \color{blue}{\left(\left(x \cdot x - x \cdot x\right) - 1\right)}\right) - \log \left(-\left(x - \mathsf{hypot}\left(1, x\right)\right)\right), x\right)
\] |
+-inverses [=>]98.4 | \[ \mathsf{copysign}\left(\log \left(-1 \cdot \left(\color{blue}{0} - 1\right)\right) - \log \left(-\left(x - \mathsf{hypot}\left(1, x\right)\right)\right), x\right)
\] |
metadata-eval [=>]98.4 | \[ \mathsf{copysign}\left(\log \left(-1 \cdot \color{blue}{-1}\right) - \log \left(-\left(x - \mathsf{hypot}\left(1, x\right)\right)\right), x\right)
\] |
metadata-eval [=>]98.4 | \[ \mathsf{copysign}\left(\log \color{blue}{1} - \log \left(-\left(x - \mathsf{hypot}\left(1, x\right)\right)\right), x\right)
\] |
metadata-eval [=>]98.4 | \[ \mathsf{copysign}\left(\color{blue}{0} - \log \left(-\left(x - \mathsf{hypot}\left(1, x\right)\right)\right), x\right)
\] |
neg-sub0 [<=]98.4 | \[ \mathsf{copysign}\left(\color{blue}{-\log \left(-\left(x - \mathsf{hypot}\left(1, x\right)\right)\right)}, x\right)
\] |
neg-sub0 [=>]98.4 | \[ \mathsf{copysign}\left(-\log \color{blue}{\left(0 - \left(x - \mathsf{hypot}\left(1, x\right)\right)\right)}, x\right)
\] |
associate--r- [=>]98.4 | \[ \mathsf{copysign}\left(-\log \color{blue}{\left(\left(0 - x\right) + \mathsf{hypot}\left(1, x\right)\right)}, x\right)
\] |
neg-sub0 [<=]98.4 | \[ \mathsf{copysign}\left(-\log \left(\color{blue}{\left(-x\right)} + \mathsf{hypot}\left(1, x\right)\right), x\right)
\] |
+-commutative [<=]98.4 | \[ \mathsf{copysign}\left(-\log \color{blue}{\left(\mathsf{hypot}\left(1, x\right) + \left(-x\right)\right)}, x\right)
\] |
sub-neg [<=]98.4 | \[ \mathsf{copysign}\left(-\log \color{blue}{\left(\mathsf{hypot}\left(1, x\right) - x\right)}, x\right)
\] |
if -0.5 < (copysign.f32 (log.f32 (+.f32 (fabs.f32 x) (sqrt.f32 (+.f32 (*.f32 x x) 1)))) x) < 0.119999997Initial program 21.2%
Applied egg-rr21.3%
[Start]21.2 | \[ \mathsf{copysign}\left(\log \left(\left|x\right| + \sqrt{x \cdot x + 1}\right), x\right)
\] |
|---|---|
*-un-lft-identity [=>]21.2 | \[ \mathsf{copysign}\left(\log \color{blue}{\left(1 \cdot \left(\left|x\right| + \sqrt{x \cdot x + 1}\right)\right)}, x\right)
\] |
*-commutative [=>]21.2 | \[ \mathsf{copysign}\left(\log \color{blue}{\left(\left(\left|x\right| + \sqrt{x \cdot x + 1}\right) \cdot 1\right)}, x\right)
\] |
log-prod [=>]21.2 | \[ \mathsf{copysign}\left(\color{blue}{\log \left(\left|x\right| + \sqrt{x \cdot x + 1}\right) + \log 1}, x\right)
\] |
add-sqr-sqrt [=>]12.9 | \[ \mathsf{copysign}\left(\log \left(\left|\color{blue}{\sqrt{x} \cdot \sqrt{x}}\right| + \sqrt{x \cdot x + 1}\right) + \log 1, x\right)
\] |
fabs-sqr [=>]12.9 | \[ \mathsf{copysign}\left(\log \left(\color{blue}{\sqrt{x} \cdot \sqrt{x}} + \sqrt{x \cdot x + 1}\right) + \log 1, x\right)
\] |
add-sqr-sqrt [<=]21.4 | \[ \mathsf{copysign}\left(\log \left(\color{blue}{x} + \sqrt{x \cdot x + 1}\right) + \log 1, x\right)
\] |
+-commutative [=>]21.4 | \[ \mathsf{copysign}\left(\log \left(x + \sqrt{\color{blue}{1 + x \cdot x}}\right) + \log 1, x\right)
\] |
hypot-1-def [=>]21.3 | \[ \mathsf{copysign}\left(\log \left(x + \color{blue}{\mathsf{hypot}\left(1, x\right)}\right) + \log 1, x\right)
\] |
metadata-eval [=>]21.3 | \[ \mathsf{copysign}\left(\log \left(x + \mathsf{hypot}\left(1, x\right)\right) + \color{blue}{0}, x\right)
\] |
Simplified21.3%
[Start]21.3 | \[ \mathsf{copysign}\left(\log \left(x + \mathsf{hypot}\left(1, x\right)\right) + 0, x\right)
\] |
|---|---|
+-rgt-identity [=>]21.3 | \[ \mathsf{copysign}\left(\color{blue}{\log \left(x + \mathsf{hypot}\left(1, x\right)\right)}, x\right)
\] |
Taylor expanded in x around 0 100.0%
if 0.119999997 < (copysign.f32 (log.f32 (+.f32 (fabs.f32 x) (sqrt.f32 (+.f32 (*.f32 x x) 1)))) x) Initial program 49.2%
Simplified98.7%
[Start]49.2 | \[ \mathsf{copysign}\left(\log \left(\left|x\right| + \sqrt{x \cdot x + 1}\right), x\right)
\] |
|---|---|
+-commutative [=>]49.2 | \[ \mathsf{copysign}\left(\log \left(\left|x\right| + \sqrt{\color{blue}{1 + x \cdot x}}\right), x\right)
\] |
hypot-1-def [=>]98.7 | \[ \mathsf{copysign}\left(\log \left(\left|x\right| + \color{blue}{\mathsf{hypot}\left(1, x\right)}\right), x\right)
\] |
Final simplification99.2%
| Alternative 1 | |
|---|---|
| Accuracy | 99.5% |
| Cost | 39144 |
| Alternative 2 | |
|---|---|
| Accuracy | 99.5% |
| Cost | 10120 |
| Alternative 3 | |
|---|---|
| Accuracy | 98.8% |
| Cost | 9896 |
| Alternative 4 | |
|---|---|
| Accuracy | 99.5% |
| Cost | 9896 |
| Alternative 5 | |
|---|---|
| Accuracy | 98.0% |
| Cost | 6760 |
| Alternative 6 | |
|---|---|
| Accuracy | 98.1% |
| Cost | 6760 |
| Alternative 7 | |
|---|---|
| Accuracy | 97.3% |
| Cost | 6696 |
| Alternative 8 | |
|---|---|
| Accuracy | 83.8% |
| Cost | 6664 |
| Alternative 9 | |
|---|---|
| Accuracy | 97.1% |
| Cost | 6664 |
| Alternative 10 | |
|---|---|
| Accuracy | 68.8% |
| Cost | 6564 |
| Alternative 11 | |
|---|---|
| Accuracy | 62.6% |
| Cost | 6532 |
| Alternative 12 | |
|---|---|
| Accuracy | 54.5% |
| Cost | 3264 |
herbie shell --seed 2023164
(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))