
(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 (+ (* x (* x -0.25)) (- (log 2.0) (log x))))
double code(double x) {
return (x * (x * -0.25)) + (log(2.0) - log(x));
}
real(8) function code(x)
real(8), intent (in) :: x
code = (x * (x * (-0.25d0))) + (log(2.0d0) - log(x))
end function
public static double code(double x) {
return (x * (x * -0.25)) + (Math.log(2.0) - Math.log(x));
}
def code(x): return (x * (x * -0.25)) + (math.log(2.0) - math.log(x))
function code(x) return Float64(Float64(x * Float64(x * -0.25)) + Float64(log(2.0) - log(x))) end
function tmp = code(x) tmp = (x * (x * -0.25)) + (log(2.0) - log(x)); end
code[x_] := N[(N[(x * N[(x * -0.25), $MachinePrecision]), $MachinePrecision] + N[(N[Log[2.0], $MachinePrecision] - N[Log[x], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
x \cdot \left(x \cdot -0.25\right) + \left(\log 2 - \log x\right)
\end{array}
Initial program 99.6%
Taylor expanded in x around 0 99.7%
+-commutative99.7%
neg-mul-199.7%
associate-+l+99.7%
unpow299.7%
*-commutative99.7%
associate-*l*99.7%
+-commutative99.7%
unsub-neg99.7%
Simplified99.7%
Final simplification99.7%
(FPCore (x) :precision binary64 (log (/ (+ 1.0 (sqrt (- 1.0 (* x x)))) x)))
double code(double x) {
return log(((1.0 + sqrt((1.0 - (x * x)))) / x));
}
real(8) function code(x)
real(8), intent (in) :: x
code = log(((1.0d0 + sqrt((1.0d0 - (x * x)))) / x))
end function
public static double code(double x) {
return Math.log(((1.0 + Math.sqrt((1.0 - (x * x)))) / x));
}
def code(x): return math.log(((1.0 + math.sqrt((1.0 - (x * x)))) / x))
function code(x) return log(Float64(Float64(1.0 + sqrt(Float64(1.0 - Float64(x * x)))) / x)) end
function tmp = code(x) tmp = log(((1.0 + sqrt((1.0 - (x * x)))) / x)); end
code[x_] := N[Log[N[(N[(1.0 + N[Sqrt[N[(1.0 - N[(x * x), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] / x), $MachinePrecision]], $MachinePrecision]
\begin{array}{l}
\\
\log \left(\frac{1 + \sqrt{1 - x \cdot x}}{x}\right)
\end{array}
Initial program 99.6%
expm1-log1p-u99.6%
expm1-udef99.6%
*-un-lft-identity99.6%
div-inv99.6%
distribute-rgt-out99.6%
Applied egg-rr99.6%
expm1-def99.6%
expm1-log1p99.6%
associate-*l/99.6%
+-commutative99.6%
*-lft-identity99.6%
+-commutative99.6%
Simplified99.6%
Final simplification99.6%
(FPCore (x) :precision binary64 (log (+ (* x -0.5) (* 2.0 (/ 1.0 x)))))
double code(double x) {
return log(((x * -0.5) + (2.0 * (1.0 / x))));
}
real(8) function code(x)
real(8), intent (in) :: x
code = log(((x * (-0.5d0)) + (2.0d0 * (1.0d0 / x))))
end function
public static double code(double x) {
return Math.log(((x * -0.5) + (2.0 * (1.0 / x))));
}
def code(x): return math.log(((x * -0.5) + (2.0 * (1.0 / x))))
function code(x) return log(Float64(Float64(x * -0.5) + Float64(2.0 * Float64(1.0 / x)))) end
function tmp = code(x) tmp = log(((x * -0.5) + (2.0 * (1.0 / x)))); end
code[x_] := N[Log[N[(N[(x * -0.5), $MachinePrecision] + N[(2.0 * N[(1.0 / x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]
\begin{array}{l}
\\
\log \left(x \cdot -0.5 + 2 \cdot \frac{1}{x}\right)
\end{array}
Initial program 99.6%
Taylor expanded in x around 0 99.6%
Final simplification99.6%
(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 99.6%
Taylor expanded in x around 0 99.2%
neg-mul-199.2%
unsub-neg99.2%
Simplified99.2%
diff-log99.1%
clear-num99.1%
log-rec99.5%
div-inv99.5%
metadata-eval99.5%
Applied egg-rr99.5%
Final simplification99.5%
(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 99.6%
Taylor expanded in x around 0 99.1%
Final simplification99.1%
(FPCore (x) :precision binary64 (+ (+ 1.0 (* -0.25 (* x x))) -1.0))
double code(double x) {
return (1.0 + (-0.25 * (x * x))) + -1.0;
}
real(8) function code(x)
real(8), intent (in) :: x
code = (1.0d0 + ((-0.25d0) * (x * x))) + (-1.0d0)
end function
public static double code(double x) {
return (1.0 + (-0.25 * (x * x))) + -1.0;
}
def code(x): return (1.0 + (-0.25 * (x * x))) + -1.0
function code(x) return Float64(Float64(1.0 + Float64(-0.25 * Float64(x * x))) + -1.0) end
function tmp = code(x) tmp = (1.0 + (-0.25 * (x * x))) + -1.0; end
code[x_] := N[(N[(1.0 + N[(-0.25 * N[(x * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + -1.0), $MachinePrecision]
\begin{array}{l}
\\
\left(1 + -0.25 \cdot \left(x \cdot x\right)\right) + -1
\end{array}
Initial program 99.6%
Taylor expanded in x around 0 99.7%
+-commutative99.7%
neg-mul-199.7%
associate-+l+99.7%
unpow299.7%
*-commutative99.7%
associate-*l*99.7%
+-commutative99.7%
unsub-neg99.7%
Simplified99.7%
Taylor expanded in x around inf 99.7%
log-rec99.7%
mul-1-neg99.7%
remove-double-neg99.7%
log-div99.6%
Simplified99.6%
Taylor expanded in x around inf 2.6%
unpow22.6%
Simplified2.6%
expm1-log1p-u2.6%
expm1-udef3.1%
log1p-udef3.1%
add-exp-log3.1%
Applied egg-rr3.1%
Final simplification3.1%
(FPCore (x) :precision binary64 (* -0.25 (* x x)))
double code(double x) {
return -0.25 * (x * x);
}
real(8) function code(x)
real(8), intent (in) :: x
code = (-0.25d0) * (x * x)
end function
public static double code(double x) {
return -0.25 * (x * x);
}
def code(x): return -0.25 * (x * x)
function code(x) return Float64(-0.25 * Float64(x * x)) end
function tmp = code(x) tmp = -0.25 * (x * x); end
code[x_] := N[(-0.25 * N[(x * x), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
-0.25 \cdot \left(x \cdot x\right)
\end{array}
Initial program 99.6%
Taylor expanded in x around 0 99.7%
+-commutative99.7%
neg-mul-199.7%
associate-+l+99.7%
unpow299.7%
*-commutative99.7%
associate-*l*99.7%
+-commutative99.7%
unsub-neg99.7%
Simplified99.7%
Taylor expanded in x around inf 99.7%
log-rec99.7%
mul-1-neg99.7%
remove-double-neg99.7%
log-div99.6%
Simplified99.6%
Taylor expanded in x around inf 2.6%
unpow22.6%
Simplified2.6%
Final simplification2.6%
herbie shell --seed 2023223
(FPCore (x)
:name "Hyperbolic arc-(co)secant"
:precision binary64
(log (+ (/ 1.0 x) (/ (sqrt (- 1.0 (* x x))) x))))