
(FPCore (x) :precision binary64 (* (/ 1.0 2.0) (log (/ (+ 1.0 x) (- 1.0 x)))))
double code(double x) {
return (1.0 / 2.0) * log(((1.0 + x) / (1.0 - x)));
}
real(8) function code(x)
real(8), intent (in) :: x
code = (1.0d0 / 2.0d0) * log(((1.0d0 + x) / (1.0d0 - x)))
end function
public static double code(double x) {
return (1.0 / 2.0) * Math.log(((1.0 + x) / (1.0 - x)));
}
def code(x): return (1.0 / 2.0) * math.log(((1.0 + x) / (1.0 - x)))
function code(x) return Float64(Float64(1.0 / 2.0) * log(Float64(Float64(1.0 + x) / Float64(1.0 - x)))) end
function tmp = code(x) tmp = (1.0 / 2.0) * log(((1.0 + x) / (1.0 - x))); end
code[x_] := N[(N[(1.0 / 2.0), $MachinePrecision] * N[Log[N[(N[(1.0 + x), $MachinePrecision] / N[(1.0 - x), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{1}{2} \cdot \log \left(\frac{1 + x}{1 - x}\right)
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 5 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (x) :precision binary64 (* (/ 1.0 2.0) (log (/ (+ 1.0 x) (- 1.0 x)))))
double code(double x) {
return (1.0 / 2.0) * log(((1.0 + x) / (1.0 - x)));
}
real(8) function code(x)
real(8), intent (in) :: x
code = (1.0d0 / 2.0d0) * log(((1.0d0 + x) / (1.0d0 - x)))
end function
public static double code(double x) {
return (1.0 / 2.0) * Math.log(((1.0 + x) / (1.0 - x)));
}
def code(x): return (1.0 / 2.0) * math.log(((1.0 + x) / (1.0 - x)))
function code(x) return Float64(Float64(1.0 / 2.0) * log(Float64(Float64(1.0 + x) / Float64(1.0 - x)))) end
function tmp = code(x) tmp = (1.0 / 2.0) * log(((1.0 + x) / (1.0 - x))); end
code[x_] := N[(N[(1.0 / 2.0), $MachinePrecision] * N[Log[N[(N[(1.0 + x), $MachinePrecision] / N[(1.0 - x), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{1}{2} \cdot \log \left(\frac{1 + x}{1 - x}\right)
\end{array}
(FPCore (x) :precision binary64 (* (+ (* -2.0 x) (+ (* -0.6666666666666666 (pow x 3.0)) (* -0.4 (pow x 5.0)))) -0.5))
double code(double x) {
return ((-2.0 * x) + ((-0.6666666666666666 * pow(x, 3.0)) + (-0.4 * pow(x, 5.0)))) * -0.5;
}
real(8) function code(x)
real(8), intent (in) :: x
code = (((-2.0d0) * x) + (((-0.6666666666666666d0) * (x ** 3.0d0)) + ((-0.4d0) * (x ** 5.0d0)))) * (-0.5d0)
end function
public static double code(double x) {
return ((-2.0 * x) + ((-0.6666666666666666 * Math.pow(x, 3.0)) + (-0.4 * Math.pow(x, 5.0)))) * -0.5;
}
def code(x): return ((-2.0 * x) + ((-0.6666666666666666 * math.pow(x, 3.0)) + (-0.4 * math.pow(x, 5.0)))) * -0.5
function code(x) return Float64(Float64(Float64(-2.0 * x) + Float64(Float64(-0.6666666666666666 * (x ^ 3.0)) + Float64(-0.4 * (x ^ 5.0)))) * -0.5) end
function tmp = code(x) tmp = ((-2.0 * x) + ((-0.6666666666666666 * (x ^ 3.0)) + (-0.4 * (x ^ 5.0)))) * -0.5; end
code[x_] := N[(N[(N[(-2.0 * x), $MachinePrecision] + N[(N[(-0.6666666666666666 * N[Power[x, 3.0], $MachinePrecision]), $MachinePrecision] + N[(-0.4 * N[Power[x, 5.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * -0.5), $MachinePrecision]
\begin{array}{l}
\\
\left(-2 \cdot x + \left(-0.6666666666666666 \cdot {x}^{3} + -0.4 \cdot {x}^{5}\right)\right) \cdot -0.5
\end{array}
Initial program 6.7%
*-commutative6.7%
log-div6.7%
sub-neg6.7%
remove-double-neg6.7%
sub-neg6.7%
+-commutative6.7%
neg-sub06.7%
associate--r-6.7%
sub-neg6.7%
log-div6.7%
neg-sub06.7%
distribute-lft-neg-in6.7%
distribute-rgt-neg-in6.7%
Simplified100.0%
Taylor expanded in x around 0 100.0%
Final simplification100.0%
(FPCore (x) :precision binary64 (* -0.5 (- (log1p (- x)) (log1p x))))
double code(double x) {
return -0.5 * (log1p(-x) - log1p(x));
}
public static double code(double x) {
return -0.5 * (Math.log1p(-x) - Math.log1p(x));
}
def code(x): return -0.5 * (math.log1p(-x) - math.log1p(x))
function code(x) return Float64(-0.5 * Float64(log1p(Float64(-x)) - log1p(x))) end
code[x_] := N[(-0.5 * N[(N[Log[1 + (-x)], $MachinePrecision] - N[Log[1 + x], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
-0.5 \cdot \left(\mathsf{log1p}\left(-x\right) - \mathsf{log1p}\left(x\right)\right)
\end{array}
Initial program 6.7%
*-commutative6.7%
log-div6.7%
sub-neg6.7%
remove-double-neg6.7%
sub-neg6.7%
+-commutative6.7%
neg-sub06.7%
associate--r-6.7%
sub-neg6.7%
log-div6.7%
neg-sub06.7%
distribute-lft-neg-in6.7%
distribute-rgt-neg-in6.7%
Simplified100.0%
Final simplification100.0%
(FPCore (x) :precision binary64 (* -0.5 (+ (* -2.0 x) (* -0.6666666666666666 (pow x 3.0)))))
double code(double x) {
return -0.5 * ((-2.0 * x) + (-0.6666666666666666 * pow(x, 3.0)));
}
real(8) function code(x)
real(8), intent (in) :: x
code = (-0.5d0) * (((-2.0d0) * x) + ((-0.6666666666666666d0) * (x ** 3.0d0)))
end function
public static double code(double x) {
return -0.5 * ((-2.0 * x) + (-0.6666666666666666 * Math.pow(x, 3.0)));
}
def code(x): return -0.5 * ((-2.0 * x) + (-0.6666666666666666 * math.pow(x, 3.0)))
function code(x) return Float64(-0.5 * Float64(Float64(-2.0 * x) + Float64(-0.6666666666666666 * (x ^ 3.0)))) end
function tmp = code(x) tmp = -0.5 * ((-2.0 * x) + (-0.6666666666666666 * (x ^ 3.0))); end
code[x_] := N[(-0.5 * N[(N[(-2.0 * x), $MachinePrecision] + N[(-0.6666666666666666 * N[Power[x, 3.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
-0.5 \cdot \left(-2 \cdot x + -0.6666666666666666 \cdot {x}^{3}\right)
\end{array}
Initial program 6.7%
*-commutative6.7%
log-div6.7%
sub-neg6.7%
remove-double-neg6.7%
sub-neg6.7%
+-commutative6.7%
neg-sub06.7%
associate--r-6.7%
sub-neg6.7%
log-div6.7%
neg-sub06.7%
distribute-lft-neg-in6.7%
distribute-rgt-neg-in6.7%
Simplified100.0%
Taylor expanded in x around 0 99.9%
Final simplification99.9%
(FPCore (x) :precision binary64 (* (* -2.0 x) -0.5))
double code(double x) {
return (-2.0 * x) * -0.5;
}
real(8) function code(x)
real(8), intent (in) :: x
code = ((-2.0d0) * x) * (-0.5d0)
end function
public static double code(double x) {
return (-2.0 * x) * -0.5;
}
def code(x): return (-2.0 * x) * -0.5
function code(x) return Float64(Float64(-2.0 * x) * -0.5) end
function tmp = code(x) tmp = (-2.0 * x) * -0.5; end
code[x_] := N[(N[(-2.0 * x), $MachinePrecision] * -0.5), $MachinePrecision]
\begin{array}{l}
\\
\left(-2 \cdot x\right) \cdot -0.5
\end{array}
Initial program 6.7%
*-commutative6.7%
log-div6.7%
sub-neg6.7%
remove-double-neg6.7%
sub-neg6.7%
+-commutative6.7%
neg-sub06.7%
associate--r-6.7%
sub-neg6.7%
log-div6.7%
neg-sub06.7%
distribute-lft-neg-in6.7%
distribute-rgt-neg-in6.7%
Simplified100.0%
Taylor expanded in x around 0 99.6%
*-commutative99.6%
Simplified99.6%
Final simplification99.6%
(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 6.7%
metadata-eval6.7%
Simplified6.7%
log-div6.7%
log1p-udef19.6%
sub-neg19.6%
log1p-udef100.0%
*-un-lft-identity100.0%
cancel-sign-sub-inv100.0%
metadata-eval100.0%
add-sqr-sqrt52.5%
sqrt-unprod35.0%
sqr-neg35.0%
sqrt-unprod3.1%
add-sqr-sqrt5.6%
neg-mul-15.6%
Applied egg-rr5.6%
sub-neg5.6%
+-inverses5.6%
Simplified5.6%
Final simplification5.6%
herbie shell --seed 2023321
(FPCore (x)
:name "Hyperbolic arc-(co)tangent"
:precision binary64
(* (/ 1.0 2.0) (log (/ (+ 1.0 x) (- 1.0 x)))))