
(FPCore (x) :precision binary32 (acosh x))
float code(float x) {
return acoshf(x);
}
function code(x) return acosh(x) end
function tmp = code(x) tmp = acosh(x); end
\begin{array}{l}
\\
\cosh^{-1} x
\end{array}
Sampling outcomes in binary32 precision:
Herbie found 7 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (x) :precision binary32 (log (+ x (sqrt (- (* x x) 1.0)))))
float code(float x) {
return logf((x + sqrtf(((x * x) - 1.0f))));
}
real(4) function code(x)
real(4), intent (in) :: x
code = log((x + sqrt(((x * x) - 1.0e0))))
end function
function code(x) return log(Float32(x + sqrt(Float32(Float32(x * x) - Float32(1.0))))) end
function tmp = code(x) tmp = log((x + sqrt(((x * x) - single(1.0))))); end
\begin{array}{l}
\\
\log \left(x + \sqrt{x \cdot x - 1}\right)
\end{array}
(FPCore (x) :precision binary32 (log (- (* x 2.0) (/ 0.5 x))))
float code(float x) {
return logf(((x * 2.0f) - (0.5f / x)));
}
real(4) function code(x)
real(4), intent (in) :: x
code = log(((x * 2.0e0) - (0.5e0 / x)))
end function
function code(x) return log(Float32(Float32(x * Float32(2.0)) - Float32(Float32(0.5) / x))) end
function tmp = code(x) tmp = log(((x * single(2.0)) - (single(0.5) / x))); end
\begin{array}{l}
\\
\log \left(x \cdot 2 - \frac{0.5}{x}\right)
\end{array}
Initial program 56.9%
Taylor expanded in x around inf 98.1%
*-commutative98.1%
associate-*r/98.1%
metadata-eval98.1%
Simplified98.1%
Final simplification98.1%
(FPCore (x) :precision binary32 (log (+ x x)))
float code(float x) {
return logf((x + x));
}
real(4) function code(x)
real(4), intent (in) :: x
code = log((x + x))
end function
function code(x) return log(Float32(x + x)) end
function tmp = code(x) tmp = log((x + x)); end
\begin{array}{l}
\\
\log \left(x + x\right)
\end{array}
Initial program 56.9%
Taylor expanded in x around inf 97.2%
Final simplification97.2%
(FPCore (x) :precision binary32 (log x))
float code(float x) {
return logf(x);
}
real(4) function code(x)
real(4), intent (in) :: x
code = log(x)
end function
function code(x) return log(x) end
function tmp = code(x) tmp = log(x); end
\begin{array}{l}
\\
\log x
\end{array}
Initial program 56.9%
pow1/256.9%
add-cube-cbrt56.9%
pow356.9%
pow-pow56.9%
fma-neg56.9%
metadata-eval56.9%
metadata-eval56.9%
Applied egg-rr56.9%
Taylor expanded in x around inf 44.1%
Final simplification44.1%
(FPCore (x) :precision binary32 (* x (+ x -1.0)))
float code(float x) {
return x * (x + -1.0f);
}
real(4) function code(x)
real(4), intent (in) :: x
code = x * (x + (-1.0e0))
end function
function code(x) return Float32(x * Float32(x + Float32(-1.0))) end
function tmp = code(x) tmp = x * (x + single(-1.0)); end
\begin{array}{l}
\\
x \cdot \left(x + -1\right)
\end{array}
Initial program 56.9%
pow1/256.9%
add-cube-cbrt56.9%
pow356.9%
pow-pow56.9%
fma-neg56.9%
metadata-eval56.9%
metadata-eval56.9%
Applied egg-rr56.9%
Taylor expanded in x around 0 -0.0%
Simplified1.7%
+-commutative1.7%
add-sqr-sqrt-0.0%
difference-of-sqr--1-0.0%
add-sqr-sqrt-0.0%
sqrt-unprod-0.0%
mul-1-neg-0.0%
mul-1-neg-0.0%
sqr-neg-0.0%
sqrt-unprod-0.0%
add-sqr-sqrt-0.0%
sqrt-prod-0.0%
add-sqr-sqrt-0.0%
Applied egg-rr8.7%
Taylor expanded in x around inf 8.8%
unpow28.8%
distribute-rgt-in8.8%
Simplified8.8%
Final simplification8.8%
(FPCore (x) :precision binary32 (- (* x x) x))
float code(float x) {
return (x * x) - x;
}
real(4) function code(x)
real(4), intent (in) :: x
code = (x * x) - x
end function
function code(x) return Float32(Float32(x * x) - x) end
function tmp = code(x) tmp = (x * x) - x; end
\begin{array}{l}
\\
x \cdot x - x
\end{array}
Initial program 56.9%
pow1/256.9%
add-cube-cbrt56.9%
pow356.9%
pow-pow56.9%
fma-neg56.9%
metadata-eval56.9%
metadata-eval56.9%
Applied egg-rr56.9%
Taylor expanded in x around 0 -0.0%
Simplified1.7%
+-commutative1.7%
add-sqr-sqrt-0.0%
difference-of-sqr--1-0.0%
add-sqr-sqrt-0.0%
sqrt-unprod-0.0%
mul-1-neg-0.0%
mul-1-neg-0.0%
sqr-neg-0.0%
sqrt-unprod-0.0%
add-sqr-sqrt-0.0%
sqrt-prod-0.0%
add-sqr-sqrt-0.0%
Applied egg-rr8.7%
Taylor expanded in x around inf 8.8%
unpow28.8%
+-commutative8.8%
mul-1-neg8.8%
unsub-neg8.8%
Simplified8.8%
Final simplification8.8%
(FPCore (x) :precision binary32 (* x x))
float code(float x) {
return x * x;
}
real(4) function code(x)
real(4), intent (in) :: x
code = x * x
end function
function code(x) return Float32(x * x) end
function tmp = code(x) tmp = x * x; end
\begin{array}{l}
\\
x \cdot x
\end{array}
Initial program 56.9%
pow1/256.9%
add-cube-cbrt56.9%
pow356.9%
pow-pow56.9%
fma-neg56.9%
metadata-eval56.9%
metadata-eval56.9%
Applied egg-rr56.9%
Taylor expanded in x around 0 -0.0%
Simplified1.7%
+-commutative1.7%
add-sqr-sqrt-0.0%
difference-of-sqr--1-0.0%
add-sqr-sqrt-0.0%
sqrt-unprod-0.0%
mul-1-neg-0.0%
mul-1-neg-0.0%
sqr-neg-0.0%
sqrt-unprod-0.0%
add-sqr-sqrt-0.0%
sqrt-prod-0.0%
add-sqr-sqrt-0.0%
Applied egg-rr8.7%
Taylor expanded in x around inf 8.8%
unpow28.8%
Simplified8.8%
Final simplification8.8%
(FPCore (x) :precision binary32 -1.0)
float code(float x) {
return -1.0f;
}
real(4) function code(x)
real(4), intent (in) :: x
code = -1.0e0
end function
function code(x) return Float32(-1.0) end
function tmp = code(x) tmp = single(-1.0); end
\begin{array}{l}
\\
-1
\end{array}
Initial program 56.9%
pow1/256.9%
add-cube-cbrt56.9%
pow356.9%
pow-pow56.9%
fma-neg56.9%
metadata-eval56.9%
metadata-eval56.9%
Applied egg-rr56.9%
Taylor expanded in x around 0 -0.0%
Simplified1.7%
Taylor expanded in x around 0 3.1%
Final simplification3.1%
(FPCore (x) :precision binary32 (log (+ x (* (sqrt (- x 1.0)) (sqrt (+ x 1.0))))))
float code(float x) {
return logf((x + (sqrtf((x - 1.0f)) * sqrtf((x + 1.0f)))));
}
real(4) function code(x)
real(4), intent (in) :: x
code = log((x + (sqrt((x - 1.0e0)) * sqrt((x + 1.0e0)))))
end function
function code(x) return log(Float32(x + Float32(sqrt(Float32(x - Float32(1.0))) * sqrt(Float32(x + Float32(1.0)))))) end
function tmp = code(x) tmp = log((x + (sqrt((x - single(1.0))) * sqrt((x + single(1.0)))))); end
\begin{array}{l}
\\
\log \left(x + \sqrt{x - 1} \cdot \sqrt{x + 1}\right)
\end{array}
herbie shell --seed 2023199
(FPCore (x)
:name "Rust f32::acosh"
:precision binary32
:pre (>= x 1.0)
:herbie-target
(log (+ x (* (sqrt (- x 1.0)) (sqrt (+ x 1.0)))))
(log (+ x (sqrt (- (* x x) 1.0)))))