| Alternative 1 | |
|---|---|
| Accuracy | 98.2% |
| Cost | 26052 |
(FPCore (x) :precision binary32 (copysign (log (+ (fabs x) (sqrt (+ (* x x) 1.0)))) x))
(FPCore (x)
:precision binary32
(if (<= x -0.10000000149011612)
(copysign (log (/ 1.0 (- (hypot 1.0 x) x))) x)
(if (<= x 0.009999999776482582)
(copysign (+ x (* -0.16666666666666666 (pow x 3.0))) x)
(copysign (log (+ 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 tmp;
if (x <= -0.10000000149011612f) {
tmp = copysignf(logf((1.0f / (hypotf(1.0f, x) - x))), x);
} else if (x <= 0.009999999776482582f) {
tmp = copysignf((x + (-0.16666666666666666f * powf(x, 3.0f))), x);
} else {
tmp = copysignf(logf((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) tmp = Float32(0.0) if (x <= Float32(-0.10000000149011612)) tmp = copysign(log(Float32(Float32(1.0) / Float32(hypot(Float32(1.0), x) - x))), x); elseif (x <= Float32(0.009999999776482582)) tmp = copysign(Float32(x + Float32(Float32(-0.16666666666666666) * (x ^ Float32(3.0)))), x); else tmp = copysign(log(Float32(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) tmp = single(0.0); if (x <= single(-0.10000000149011612)) tmp = sign(x) * abs(log((single(1.0) / (hypot(single(1.0), x) - x)))); elseif (x <= single(0.009999999776482582)) tmp = sign(x) * abs((x + (single(-0.16666666666666666) * (x ^ single(3.0))))); else tmp = sign(x) * abs(log((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}
\mathbf{if}\;x \leq -0.10000000149011612:\\
\;\;\;\;\mathsf{copysign}\left(\log \left(\frac{1}{\mathsf{hypot}\left(1, x\right) - x}\right), x\right)\\
\mathbf{elif}\;x \leq 0.009999999776482582:\\
\;\;\;\;\mathsf{copysign}\left(x + -0.16666666666666666 \cdot {x}^{3}, x\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{copysign}\left(\log \left(x + \mathsf{hypot}\left(1, x\right)\right), x\right)\\
\end{array}
| Original | 36.6% |
|---|---|
| Target | 99.6% |
| Herbie | 99.4% |
if x < -0.100000001Initial program 52.0%
Applied egg-rr99.2%
[Start]52.0 | \[ \mathsf{copysign}\left(\log \left(\left|x\right| + \sqrt{x \cdot x + 1}\right), x\right)
\] |
|---|---|
flip-+ [=>]8.7 | \[ \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)
\] |
div-inv [=>]8.7 | \[ \mathsf{copysign}\left(\log \color{blue}{\left(\left(\left|x\right| \cdot \left|x\right| - \sqrt{x \cdot x + 1} \cdot \sqrt{x \cdot x + 1}\right) \cdot \frac{1}{\left|x\right| - \sqrt{x \cdot x + 1}}\right)}, x\right)
\] |
add-sqr-sqrt [<=]8.5 | \[ \mathsf{copysign}\left(\log \left(\left(\left|x\right| \cdot \left|x\right| - \color{blue}{\left(x \cdot x + 1\right)}\right) \cdot \frac{1}{\left|x\right| - \sqrt{x \cdot x + 1}}\right), x\right)
\] |
associate--r+ [=>]8.5 | \[ \mathsf{copysign}\left(\log \left(\color{blue}{\left(\left(\left|x\right| \cdot \left|x\right| - x \cdot x\right) - 1\right)} \cdot \frac{1}{\left|x\right| - \sqrt{x \cdot x + 1}}\right), x\right)
\] |
sqr-abs [<=]8.5 | \[ \mathsf{copysign}\left(\log \left(\left(\left(\left|x\right| \cdot \left|x\right| - \color{blue}{\left|x\right| \cdot \left|x\right|}\right) - 1\right) \cdot \frac{1}{\left|x\right| - \sqrt{x \cdot x + 1}}\right), x\right)
\] |
distribute-lft-out-- [=>]8.5 | \[ \mathsf{copysign}\left(\log \left(\left(\color{blue}{\left|x\right| \cdot \left(\left|x\right| - \left|x\right|\right)} - 1\right) \cdot \frac{1}{\left|x\right| - \sqrt{x \cdot x + 1}}\right), x\right)
\] |
add-sqr-sqrt [=>]-0.0 | \[ \mathsf{copysign}\left(\log \left(\left(\left|\color{blue}{\sqrt{x} \cdot \sqrt{x}}\right| \cdot \left(\left|x\right| - \left|x\right|\right) - 1\right) \cdot \frac{1}{\left|x\right| - \sqrt{x \cdot x + 1}}\right), x\right)
\] |
fabs-sqr [=>]-0.0 | \[ \mathsf{copysign}\left(\log \left(\left(\color{blue}{\left(\sqrt{x} \cdot \sqrt{x}\right)} \cdot \left(\left|x\right| - \left|x\right|\right) - 1\right) \cdot \frac{1}{\left|x\right| - \sqrt{x \cdot x + 1}}\right), x\right)
\] |
add-sqr-sqrt [<=]8.5 | \[ \mathsf{copysign}\left(\log \left(\left(\color{blue}{x} \cdot \left(\left|x\right| - \left|x\right|\right) - 1\right) \cdot \frac{1}{\left|x\right| - \sqrt{x \cdot x + 1}}\right), x\right)
\] |
add-sqr-sqrt [=>]-0.0 | \[ \mathsf{copysign}\left(\log \left(\left(x \cdot \left(\left|\color{blue}{\sqrt{x} \cdot \sqrt{x}}\right| - \left|x\right|\right) - 1\right) \cdot \frac{1}{\left|x\right| - \sqrt{x \cdot x + 1}}\right), x\right)
\] |
fabs-sqr [=>]-0.0 | \[ \mathsf{copysign}\left(\log \left(\left(x \cdot \left(\color{blue}{\sqrt{x} \cdot \sqrt{x}} - \left|x\right|\right) - 1\right) \cdot \frac{1}{\left|x\right| - \sqrt{x \cdot x + 1}}\right), x\right)
\] |
add-sqr-sqrt [<=]0.7 | \[ \mathsf{copysign}\left(\log \left(\left(x \cdot \left(\color{blue}{x} - \left|x\right|\right) - 1\right) \cdot \frac{1}{\left|x\right| - \sqrt{x \cdot x + 1}}\right), x\right)
\] |
add-sqr-sqrt [=>]-0.0 | \[ \mathsf{copysign}\left(\log \left(\left(x \cdot \left(x - \left|\color{blue}{\sqrt{x} \cdot \sqrt{x}}\right|\right) - 1\right) \cdot \frac{1}{\left|x\right| - \sqrt{x \cdot x + 1}}\right), x\right)
\] |
fabs-sqr [=>]-0.0 | \[ \mathsf{copysign}\left(\log \left(\left(x \cdot \left(x - \color{blue}{\sqrt{x} \cdot \sqrt{x}}\right) - 1\right) \cdot \frac{1}{\left|x\right| - \sqrt{x \cdot x + 1}}\right), x\right)
\] |
add-sqr-sqrt [<=]8.5 | \[ \mathsf{copysign}\left(\log \left(\left(x \cdot \left(x - \color{blue}{x}\right) - 1\right) \cdot \frac{1}{\left|x\right| - \sqrt{x \cdot x + 1}}\right), x\right)
\] |
Simplified99.2%
[Start]99.2 | \[ \mathsf{copysign}\left(\log \left(\left(x \cdot \left(x - x\right) - 1\right) \cdot \frac{1}{x - \mathsf{hypot}\left(1, x\right)}\right), x\right)
\] |
|---|---|
+-inverses [=>]99.2 | \[ \mathsf{copysign}\left(\log \left(\left(x \cdot \color{blue}{0} - 1\right) \cdot \frac{1}{x - \mathsf{hypot}\left(1, x\right)}\right), x\right)
\] |
mul0-rgt [=>]99.2 | \[ \mathsf{copysign}\left(\log \left(\left(\color{blue}{0} - 1\right) \cdot \frac{1}{x - \mathsf{hypot}\left(1, x\right)}\right), x\right)
\] |
metadata-eval [=>]99.2 | \[ \mathsf{copysign}\left(\log \left(\color{blue}{-1} \cdot \frac{1}{x - \mathsf{hypot}\left(1, x\right)}\right), x\right)
\] |
associate-*r/ [=>]99.2 | \[ \mathsf{copysign}\left(\log \color{blue}{\left(\frac{-1 \cdot 1}{x - \mathsf{hypot}\left(1, x\right)}\right)}, x\right)
\] |
metadata-eval [=>]99.2 | \[ \mathsf{copysign}\left(\log \left(\frac{\color{blue}{-1}}{x - \mathsf{hypot}\left(1, x\right)}\right), x\right)
\] |
metadata-eval [<=]99.2 | \[ \mathsf{copysign}\left(\log \left(\frac{\color{blue}{\frac{1}{-1}}}{x - \mathsf{hypot}\left(1, x\right)}\right), x\right)
\] |
associate-/r* [<=]99.2 | \[ \mathsf{copysign}\left(\log \color{blue}{\left(\frac{1}{-1 \cdot \left(x - \mathsf{hypot}\left(1, x\right)\right)}\right)}, x\right)
\] |
neg-mul-1 [<=]99.2 | \[ \mathsf{copysign}\left(\log \left(\frac{1}{\color{blue}{-\left(x - \mathsf{hypot}\left(1, x\right)\right)}}\right), x\right)
\] |
neg-sub0 [=>]99.2 | \[ \mathsf{copysign}\left(\log \left(\frac{1}{\color{blue}{0 - \left(x - \mathsf{hypot}\left(1, x\right)\right)}}\right), x\right)
\] |
associate--r- [=>]99.2 | \[ \mathsf{copysign}\left(\log \left(\frac{1}{\color{blue}{\left(0 - x\right) + \mathsf{hypot}\left(1, x\right)}}\right), x\right)
\] |
neg-sub0 [<=]99.2 | \[ \mathsf{copysign}\left(\log \left(\frac{1}{\color{blue}{\left(-x\right)} + \mathsf{hypot}\left(1, x\right)}\right), x\right)
\] |
mul-1-neg [<=]99.2 | \[ \mathsf{copysign}\left(\log \left(\frac{1}{\color{blue}{-1 \cdot x} + \mathsf{hypot}\left(1, x\right)}\right), x\right)
\] |
+-commutative [<=]99.2 | \[ \mathsf{copysign}\left(\log \left(\frac{1}{\color{blue}{\mathsf{hypot}\left(1, x\right) + -1 \cdot x}}\right), x\right)
\] |
mul-1-neg [=>]99.2 | \[ \mathsf{copysign}\left(\log \left(\frac{1}{\mathsf{hypot}\left(1, x\right) + \color{blue}{\left(-x\right)}}\right), x\right)
\] |
sub-neg [<=]99.2 | \[ \mathsf{copysign}\left(\log \left(\frac{1}{\color{blue}{\mathsf{hypot}\left(1, x\right) - x}}\right), x\right)
\] |
if -0.100000001 < x < 0.00999999978Initial program 19.5%
Applied egg-rr97.3%
[Start]19.5 | \[ \mathsf{copysign}\left(\log \left(\left|x\right| + \sqrt{x \cdot x + 1}\right), x\right)
\] |
|---|---|
log1p-expm1-u [=>]19.5 | \[ \mathsf{copysign}\left(\color{blue}{\mathsf{log1p}\left(\mathsf{expm1}\left(\log \left(\left|x\right| + \sqrt{x \cdot x + 1}\right)\right)\right)}, x\right)
\] |
expm1-udef [=>]19.5 | \[ \mathsf{copysign}\left(\mathsf{log1p}\left(\color{blue}{e^{\log \left(\left|x\right| + \sqrt{x \cdot x + 1}\right)} - 1}\right), x\right)
\] |
add-exp-log [<=]19.5 | \[ \mathsf{copysign}\left(\mathsf{log1p}\left(\color{blue}{\left(\left|x\right| + \sqrt{x \cdot x + 1}\right)} - 1\right), x\right)
\] |
associate--l+ [=>]97.3 | \[ \mathsf{copysign}\left(\mathsf{log1p}\left(\color{blue}{\left|x\right| + \left(\sqrt{x \cdot x + 1} - 1\right)}\right), x\right)
\] |
add-sqr-sqrt [=>]46.8 | \[ \mathsf{copysign}\left(\mathsf{log1p}\left(\left|\color{blue}{\sqrt{x} \cdot \sqrt{x}}\right| + \left(\sqrt{x \cdot x + 1} - 1\right)\right), x\right)
\] |
fabs-sqr [=>]46.8 | \[ \mathsf{copysign}\left(\mathsf{log1p}\left(\color{blue}{\sqrt{x} \cdot \sqrt{x}} + \left(\sqrt{x \cdot x + 1} - 1\right)\right), x\right)
\] |
add-sqr-sqrt [<=]97.3 | \[ \mathsf{copysign}\left(\mathsf{log1p}\left(\color{blue}{x} + \left(\sqrt{x \cdot x + 1} - 1\right)\right), x\right)
\] |
+-commutative [=>]97.3 | \[ \mathsf{copysign}\left(\mathsf{log1p}\left(x + \left(\sqrt{\color{blue}{1 + x \cdot x}} - 1\right)\right), x\right)
\] |
sqr-abs [<=]97.3 | \[ \mathsf{copysign}\left(\mathsf{log1p}\left(x + \left(\sqrt{1 + \color{blue}{\left|x\right| \cdot \left|x\right|}} - 1\right)\right), x\right)
\] |
hypot-1-def [=>]97.3 | \[ \mathsf{copysign}\left(\mathsf{log1p}\left(x + \left(\color{blue}{\mathsf{hypot}\left(1, \left|x\right|\right)} - 1\right)\right), x\right)
\] |
add-sqr-sqrt [=>]47.5 | \[ \mathsf{copysign}\left(\mathsf{log1p}\left(x + \left(\mathsf{hypot}\left(1, \left|\color{blue}{\sqrt{x} \cdot \sqrt{x}}\right|\right) - 1\right)\right), x\right)
\] |
fabs-sqr [=>]47.5 | \[ \mathsf{copysign}\left(\mathsf{log1p}\left(x + \left(\mathsf{hypot}\left(1, \color{blue}{\sqrt{x} \cdot \sqrt{x}}\right) - 1\right)\right), x\right)
\] |
add-sqr-sqrt [<=]97.3 | \[ \mathsf{copysign}\left(\mathsf{log1p}\left(x + \left(\mathsf{hypot}\left(1, \color{blue}{x}\right) - 1\right)\right), x\right)
\] |
Taylor expanded in x around 0 99.9%
if 0.00999999978 < x Initial program 53.2%
Applied egg-rr98.7%
[Start]53.2 | \[ \mathsf{copysign}\left(\log \left(\left|x\right| + \sqrt{x \cdot x + 1}\right), x\right)
\] |
|---|---|
*-un-lft-identity [=>]53.2 | \[ \mathsf{copysign}\left(\log \color{blue}{\left(1 \cdot \left(\left|x\right| + \sqrt{x \cdot x + 1}\right)\right)}, x\right)
\] |
log-prod [=>]53.2 | \[ \mathsf{copysign}\left(\color{blue}{\log 1 + \log \left(\left|x\right| + \sqrt{x \cdot x + 1}\right)}, x\right)
\] |
+-commutative [=>]53.2 | \[ \mathsf{copysign}\left(\color{blue}{\log \left(\left|x\right| + \sqrt{x \cdot x + 1}\right) + \log 1}, x\right)
\] |
add-sqr-sqrt [=>]53.2 | \[ \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 [=>]53.2 | \[ \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 [<=]53.2 | \[ \mathsf{copysign}\left(\log \left(\color{blue}{x} + \sqrt{x \cdot x + 1}\right) + \log 1, x\right)
\] |
+-commutative [=>]53.2 | \[ \mathsf{copysign}\left(\log \left(x + \sqrt{\color{blue}{1 + x \cdot x}}\right) + \log 1, x\right)
\] |
sqr-abs [<=]53.2 | \[ \mathsf{copysign}\left(\log \left(x + \sqrt{1 + \color{blue}{\left|x\right| \cdot \left|x\right|}}\right) + \log 1, x\right)
\] |
hypot-1-def [=>]98.7 | \[ \mathsf{copysign}\left(\log \left(x + \color{blue}{\mathsf{hypot}\left(1, \left|x\right|\right)}\right) + \log 1, x\right)
\] |
add-sqr-sqrt [=>]98.7 | \[ \mathsf{copysign}\left(\log \left(x + \mathsf{hypot}\left(1, \left|\color{blue}{\sqrt{x} \cdot \sqrt{x}}\right|\right)\right) + \log 1, x\right)
\] |
fabs-sqr [=>]98.7 | \[ \mathsf{copysign}\left(\log \left(x + \mathsf{hypot}\left(1, \color{blue}{\sqrt{x} \cdot \sqrt{x}}\right)\right) + \log 1, x\right)
\] |
add-sqr-sqrt [<=]98.7 | \[ \mathsf{copysign}\left(\log \left(x + \mathsf{hypot}\left(1, \color{blue}{x}\right)\right) + \log 1, x\right)
\] |
metadata-eval [=>]98.7 | \[ \mathsf{copysign}\left(\log \left(x + \mathsf{hypot}\left(1, x\right)\right) + \color{blue}{0}, x\right)
\] |
Simplified98.7%
[Start]98.7 | \[ \mathsf{copysign}\left(\log \left(x + \mathsf{hypot}\left(1, x\right)\right) + 0, x\right)
\] |
|---|---|
+-rgt-identity [=>]98.7 | \[ \mathsf{copysign}\left(\color{blue}{\log \left(x + \mathsf{hypot}\left(1, x\right)\right)}, x\right)
\] |
Final simplification99.4%
| Alternative 1 | |
|---|---|
| Accuracy | 98.2% |
| Cost | 26052 |
| Alternative 2 | |
|---|---|
| Accuracy | 98.9% |
| Cost | 9896 |
| Alternative 3 | |
|---|---|
| Accuracy | 98.2% |
| Cost | 9892 |
| Alternative 4 | |
|---|---|
| Accuracy | 98.5% |
| Cost | 6984 |
| Alternative 5 | |
|---|---|
| Accuracy | 98.3% |
| Cost | 6792 |
| Alternative 6 | |
|---|---|
| Accuracy | 98.0% |
| Cost | 6760 |
| Alternative 7 | |
|---|---|
| Accuracy | 83.8% |
| Cost | 6664 |
| Alternative 8 | |
|---|---|
| Accuracy | 97.1% |
| Cost | 6664 |
| Alternative 9 | |
|---|---|
| Accuracy | 68.9% |
| Cost | 6564 |
| Alternative 10 | |
|---|---|
| Accuracy | 62.7% |
| Cost | 6532 |
| Alternative 11 | |
|---|---|
| Accuracy | 62.7% |
| Cost | 6532 |
| Alternative 12 | |
|---|---|
| Accuracy | 54.7% |
| Cost | 3264 |
herbie shell --seed 2023138
(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))