
(FPCore (x) :precision binary64 (log (+ (/ 1.0 x) (/ (sqrt (- 1.0 (* x x))) x))))
double code(double x) {
return log(((1.0 / x) + (sqrt((1.0 - (x * x))) / x)));
}
real(8) function code(x)
real(8), intent (in) :: x
code = log(((1.0d0 / x) + (sqrt((1.0d0 - (x * x))) / x)))
end function
public static double code(double x) {
return Math.log(((1.0 / x) + (Math.sqrt((1.0 - (x * x))) / x)));
}
def code(x): return math.log(((1.0 / x) + (math.sqrt((1.0 - (x * x))) / x)))
function code(x) return log(Float64(Float64(1.0 / x) + Float64(sqrt(Float64(1.0 - Float64(x * x))) / x))) end
function tmp = code(x) tmp = log(((1.0 / x) + (sqrt((1.0 - (x * x))) / x))); end
code[x_] := N[Log[N[(N[(1.0 / x), $MachinePrecision] + N[(N[Sqrt[N[(1.0 - N[(x * x), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] / x), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]
\begin{array}{l}
\\
\log \left(\frac{1}{x} + \frac{\sqrt{1 - x \cdot x}}{x}\right)
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 7 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (x) :precision binary64 (log (+ (/ 1.0 x) (/ (sqrt (- 1.0 (* x x))) x))))
double code(double x) {
return log(((1.0 / x) + (sqrt((1.0 - (x * x))) / x)));
}
real(8) function code(x)
real(8), intent (in) :: x
code = log(((1.0d0 / x) + (sqrt((1.0d0 - (x * x))) / x)))
end function
public static double code(double x) {
return Math.log(((1.0 / x) + (Math.sqrt((1.0 - (x * x))) / x)));
}
def code(x): return math.log(((1.0 / x) + (math.sqrt((1.0 - (x * x))) / x)))
function code(x) return log(Float64(Float64(1.0 / x) + Float64(sqrt(Float64(1.0 - Float64(x * x))) / x))) end
function tmp = code(x) tmp = log(((1.0 / x) + (sqrt((1.0 - (x * x))) / x))); end
code[x_] := N[Log[N[(N[(1.0 / x), $MachinePrecision] + N[(N[Sqrt[N[(1.0 - N[(x * x), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] / x), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]
\begin{array}{l}
\\
\log \left(\frac{1}{x} + \frac{\sqrt{1 - x \cdot x}}{x}\right)
\end{array}
(FPCore (x) :precision binary64 (- (log (/ x (+ 1.0 (sqrt (- 1.0 (pow x 2.0))))))))
double code(double x) {
return -log((x / (1.0 + sqrt((1.0 - pow(x, 2.0))))));
}
real(8) function code(x)
real(8), intent (in) :: x
code = -log((x / (1.0d0 + sqrt((1.0d0 - (x ** 2.0d0))))))
end function
public static double code(double x) {
return -Math.log((x / (1.0 + Math.sqrt((1.0 - Math.pow(x, 2.0))))));
}
def code(x): return -math.log((x / (1.0 + math.sqrt((1.0 - math.pow(x, 2.0))))))
function code(x) return Float64(-log(Float64(x / Float64(1.0 + sqrt(Float64(1.0 - (x ^ 2.0))))))) end
function tmp = code(x) tmp = -log((x / (1.0 + sqrt((1.0 - (x ^ 2.0)))))); end
code[x_] := (-N[Log[N[(x / N[(1.0 + N[Sqrt[N[(1.0 - N[Power[x, 2.0], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision])
\begin{array}{l}
\\
-\log \left(\frac{x}{1 + \sqrt{1 - {x}^{2}}}\right)
\end{array}
Initial program 100.0%
add-sqr-sqrt100.0%
log-prod100.0%
*-un-lft-identity100.0%
div-inv100.0%
distribute-rgt-out100.0%
pow2100.0%
*-un-lft-identity100.0%
div-inv100.0%
distribute-rgt-out100.0%
Applied egg-rr100.0%
count-2100.0%
associate-*l/100.0%
*-lft-identity100.0%
Simplified100.0%
*-commutative100.0%
add-log-exp100.0%
exp-to-pow100.0%
pow2100.0%
add-sqr-sqrt100.0%
clear-num100.0%
log-div100.0%
metadata-eval100.0%
Applied egg-rr100.0%
neg-sub0100.0%
Simplified100.0%
(FPCore (x) :precision binary64 (log (+ (/ 1.0 x) (/ (sqrt (- 1.0 (* x x))) x))))
double code(double x) {
return log(((1.0 / x) + (sqrt((1.0 - (x * x))) / x)));
}
real(8) function code(x)
real(8), intent (in) :: x
code = log(((1.0d0 / x) + (sqrt((1.0d0 - (x * x))) / x)))
end function
public static double code(double x) {
return Math.log(((1.0 / x) + (Math.sqrt((1.0 - (x * x))) / x)));
}
def code(x): return math.log(((1.0 / x) + (math.sqrt((1.0 - (x * x))) / x)))
function code(x) return log(Float64(Float64(1.0 / x) + Float64(sqrt(Float64(1.0 - Float64(x * x))) / x))) end
function tmp = code(x) tmp = log(((1.0 / x) + (sqrt((1.0 - (x * x))) / x))); end
code[x_] := N[Log[N[(N[(1.0 / x), $MachinePrecision] + N[(N[Sqrt[N[(1.0 - N[(x * x), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] / x), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]
\begin{array}{l}
\\
\log \left(\frac{1}{x} + \frac{\sqrt{1 - x \cdot x}}{x}\right)
\end{array}
Initial program 100.0%
(FPCore (x) :precision binary64 (log (/ (+ 2.0 (* (* x x) -0.5)) x)))
double code(double x) {
return log(((2.0 + ((x * x) * -0.5)) / x));
}
real(8) function code(x)
real(8), intent (in) :: x
code = log(((2.0d0 + ((x * x) * (-0.5d0))) / x))
end function
public static double code(double x) {
return Math.log(((2.0 + ((x * x) * -0.5)) / x));
}
def code(x): return math.log(((2.0 + ((x * x) * -0.5)) / x))
function code(x) return log(Float64(Float64(2.0 + Float64(Float64(x * x) * -0.5)) / x)) end
function tmp = code(x) tmp = log(((2.0 + ((x * x) * -0.5)) / x)); end
code[x_] := N[Log[N[(N[(2.0 + N[(N[(x * x), $MachinePrecision] * -0.5), $MachinePrecision]), $MachinePrecision] / x), $MachinePrecision]], $MachinePrecision]
\begin{array}{l}
\\
\log \left(\frac{2 + \left(x \cdot x\right) \cdot -0.5}{x}\right)
\end{array}
Initial program 100.0%
Taylor expanded in x around 0 99.5%
*-commutative99.5%
Simplified99.5%
unpow299.5%
Applied egg-rr99.5%
(FPCore (x) :precision binary64 (- (log (* x 0.5))))
double code(double x) {
return -log((x * 0.5));
}
real(8) function code(x)
real(8), intent (in) :: x
code = -log((x * 0.5d0))
end function
public static double code(double x) {
return -Math.log((x * 0.5));
}
def code(x): return -math.log((x * 0.5))
function code(x) return Float64(-log(Float64(x * 0.5))) end
function tmp = code(x) tmp = -log((x * 0.5)); end
code[x_] := (-N[Log[N[(x * 0.5), $MachinePrecision]], $MachinePrecision])
\begin{array}{l}
\\
-\log \left(x \cdot 0.5\right)
\end{array}
Initial program 100.0%
Taylor expanded in x around 0 98.8%
clear-num98.8%
log-div98.9%
metadata-eval98.9%
div-inv98.9%
metadata-eval98.9%
Applied egg-rr98.9%
neg-sub098.9%
Simplified98.9%
(FPCore (x) :precision binary64 (log (/ 2.0 x)))
double code(double x) {
return log((2.0 / x));
}
real(8) function code(x)
real(8), intent (in) :: x
code = log((2.0d0 / x))
end function
public static double code(double x) {
return Math.log((2.0 / x));
}
def code(x): return math.log((2.0 / x))
function code(x) return log(Float64(2.0 / x)) end
function tmp = code(x) tmp = log((2.0 / x)); end
code[x_] := N[Log[N[(2.0 / x), $MachinePrecision]], $MachinePrecision]
\begin{array}{l}
\\
\log \left(\frac{2}{x}\right)
\end{array}
Initial program 100.0%
Taylor expanded in x around 0 98.8%
(FPCore (x) :precision binary64 (- (log1p -1.0)))
double code(double x) {
return -log1p(-1.0);
}
public static double code(double x) {
return -Math.log1p(-1.0);
}
def code(x): return -math.log1p(-1.0)
function code(x) return Float64(-log1p(-1.0)) end
code[x_] := (-N[Log[1 + -1.0], $MachinePrecision])
\begin{array}{l}
\\
-\mathsf{log1p}\left(-1\right)
\end{array}
Initial program 100.0%
Taylor expanded in x around 0 98.8%
clear-num98.8%
log-div98.9%
metadata-eval98.9%
div-inv98.9%
metadata-eval98.9%
Applied egg-rr98.9%
neg-sub098.9%
Simplified98.9%
add-sqr-sqrt0.0%
sqrt-unprod1.5%
sqr-neg1.5%
sqrt-unprod1.5%
add-sqr-sqrt1.5%
neg-log1.5%
*-commutative1.5%
associate-/r*1.5%
metadata-eval1.5%
log1p-expm1-u1.5%
expm1-define1.5%
add-exp-log1.5%
sub-neg1.5%
Applied egg-rr6.3%
Taylor expanded in x around 0 3.1%
(FPCore (x) :precision binary64 0.0)
double code(double x) {
return 0.0;
}
real(8) function code(x)
real(8), intent (in) :: x
code = 0.0d0
end function
public static double code(double x) {
return 0.0;
}
def code(x): return 0.0
function code(x) return 0.0 end
function tmp = code(x) tmp = 0.0; end
code[x_] := 0.0
\begin{array}{l}
\\
0
\end{array}
Initial program 100.0%
Taylor expanded in x around 0 98.8%
clear-num98.8%
log-div98.9%
metadata-eval98.9%
div-inv98.9%
metadata-eval98.9%
Applied egg-rr98.9%
neg-sub098.9%
Simplified98.9%
Applied egg-rr3.1%
*-inverses3.1%
Simplified3.1%
metadata-eval3.1%
Applied egg-rr3.1%
Final simplification3.1%
herbie shell --seed 2024123
(FPCore (x)
:name "Hyperbolic arc-(co)secant"
:precision binary64
(log (+ (/ 1.0 x) (/ (sqrt (- 1.0 (* x x))) x))))