
(FPCore (x) :precision binary64 (log (+ x (sqrt (+ (* x x) 1.0)))))
double code(double x) {
return log((x + sqrt(((x * x) + 1.0))));
}
real(8) function code(x)
real(8), intent (in) :: x
code = log((x + sqrt(((x * x) + 1.0d0))))
end function
public static double code(double x) {
return Math.log((x + Math.sqrt(((x * x) + 1.0))));
}
def code(x): return math.log((x + math.sqrt(((x * x) + 1.0))))
function code(x) return log(Float64(x + sqrt(Float64(Float64(x * x) + 1.0)))) end
function tmp = code(x) tmp = log((x + sqrt(((x * x) + 1.0)))); end
code[x_] := N[Log[N[(x + N[Sqrt[N[(N[(x * x), $MachinePrecision] + 1.0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]
\begin{array}{l}
\\
\log \left(x + \sqrt{x \cdot x + 1}\right)
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 6 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (x) :precision binary64 (log (+ x (sqrt (+ (* x x) 1.0)))))
double code(double x) {
return log((x + sqrt(((x * x) + 1.0))));
}
real(8) function code(x)
real(8), intent (in) :: x
code = log((x + sqrt(((x * x) + 1.0d0))))
end function
public static double code(double x) {
return Math.log((x + Math.sqrt(((x * x) + 1.0))));
}
def code(x): return math.log((x + math.sqrt(((x * x) + 1.0))))
function code(x) return log(Float64(x + sqrt(Float64(Float64(x * x) + 1.0)))) end
function tmp = code(x) tmp = log((x + sqrt(((x * x) + 1.0)))); end
code[x_] := N[Log[N[(x + N[Sqrt[N[(N[(x * x), $MachinePrecision] + 1.0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]
\begin{array}{l}
\\
\log \left(x + \sqrt{x \cdot x + 1}\right)
\end{array}
(FPCore (x)
:precision binary64
(let* ((t_0 (* x (* x (* x x)))))
(if (<= x -1.25)
(log (/ (+ (/ 0.125 (* x x)) (+ -0.5 (/ -0.0625 t_0))) x))
(if (<= x 1.25)
(*
x
(/
(- 1.0 (* t_0 (+ 0.027777777777777776 (* (* x x) -0.025))))
(-
1.0
(*
(* x x)
(+
-0.16666666666666666
(* (* x x) (+ 0.075 (* x (* x -0.044642857142857144)))))))))
(log (+ (* x 2.0) (/ 0.5 x)))))))
double code(double x) {
double t_0 = x * (x * (x * x));
double tmp;
if (x <= -1.25) {
tmp = log((((0.125 / (x * x)) + (-0.5 + (-0.0625 / t_0))) / x));
} else if (x <= 1.25) {
tmp = x * ((1.0 - (t_0 * (0.027777777777777776 + ((x * x) * -0.025)))) / (1.0 - ((x * x) * (-0.16666666666666666 + ((x * x) * (0.075 + (x * (x * -0.044642857142857144))))))));
} else {
tmp = log(((x * 2.0) + (0.5 / x)));
}
return tmp;
}
real(8) function code(x)
real(8), intent (in) :: x
real(8) :: t_0
real(8) :: tmp
t_0 = x * (x * (x * x))
if (x <= (-1.25d0)) then
tmp = log((((0.125d0 / (x * x)) + ((-0.5d0) + ((-0.0625d0) / t_0))) / x))
else if (x <= 1.25d0) then
tmp = x * ((1.0d0 - (t_0 * (0.027777777777777776d0 + ((x * x) * (-0.025d0))))) / (1.0d0 - ((x * x) * ((-0.16666666666666666d0) + ((x * x) * (0.075d0 + (x * (x * (-0.044642857142857144d0)))))))))
else
tmp = log(((x * 2.0d0) + (0.5d0 / x)))
end if
code = tmp
end function
public static double code(double x) {
double t_0 = x * (x * (x * x));
double tmp;
if (x <= -1.25) {
tmp = Math.log((((0.125 / (x * x)) + (-0.5 + (-0.0625 / t_0))) / x));
} else if (x <= 1.25) {
tmp = x * ((1.0 - (t_0 * (0.027777777777777776 + ((x * x) * -0.025)))) / (1.0 - ((x * x) * (-0.16666666666666666 + ((x * x) * (0.075 + (x * (x * -0.044642857142857144))))))));
} else {
tmp = Math.log(((x * 2.0) + (0.5 / x)));
}
return tmp;
}
def code(x): t_0 = x * (x * (x * x)) tmp = 0 if x <= -1.25: tmp = math.log((((0.125 / (x * x)) + (-0.5 + (-0.0625 / t_0))) / x)) elif x <= 1.25: tmp = x * ((1.0 - (t_0 * (0.027777777777777776 + ((x * x) * -0.025)))) / (1.0 - ((x * x) * (-0.16666666666666666 + ((x * x) * (0.075 + (x * (x * -0.044642857142857144)))))))) else: tmp = math.log(((x * 2.0) + (0.5 / x))) return tmp
function code(x) t_0 = Float64(x * Float64(x * Float64(x * x))) tmp = 0.0 if (x <= -1.25) tmp = log(Float64(Float64(Float64(0.125 / Float64(x * x)) + Float64(-0.5 + Float64(-0.0625 / t_0))) / x)); elseif (x <= 1.25) tmp = Float64(x * Float64(Float64(1.0 - Float64(t_0 * Float64(0.027777777777777776 + Float64(Float64(x * x) * -0.025)))) / Float64(1.0 - Float64(Float64(x * x) * Float64(-0.16666666666666666 + Float64(Float64(x * x) * Float64(0.075 + Float64(x * Float64(x * -0.044642857142857144))))))))); else tmp = log(Float64(Float64(x * 2.0) + Float64(0.5 / x))); end return tmp end
function tmp_2 = code(x) t_0 = x * (x * (x * x)); tmp = 0.0; if (x <= -1.25) tmp = log((((0.125 / (x * x)) + (-0.5 + (-0.0625 / t_0))) / x)); elseif (x <= 1.25) tmp = x * ((1.0 - (t_0 * (0.027777777777777776 + ((x * x) * -0.025)))) / (1.0 - ((x * x) * (-0.16666666666666666 + ((x * x) * (0.075 + (x * (x * -0.044642857142857144)))))))); else tmp = log(((x * 2.0) + (0.5 / x))); end tmp_2 = tmp; end
code[x_] := Block[{t$95$0 = N[(x * N[(x * N[(x * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[x, -1.25], N[Log[N[(N[(N[(0.125 / N[(x * x), $MachinePrecision]), $MachinePrecision] + N[(-0.5 + N[(-0.0625 / t$95$0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / x), $MachinePrecision]], $MachinePrecision], If[LessEqual[x, 1.25], N[(x * N[(N[(1.0 - N[(t$95$0 * N[(0.027777777777777776 + N[(N[(x * x), $MachinePrecision] * -0.025), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(1.0 - N[(N[(x * x), $MachinePrecision] * N[(-0.16666666666666666 + N[(N[(x * x), $MachinePrecision] * N[(0.075 + N[(x * N[(x * -0.044642857142857144), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[Log[N[(N[(x * 2.0), $MachinePrecision] + N[(0.5 / x), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := x \cdot \left(x \cdot \left(x \cdot x\right)\right)\\
\mathbf{if}\;x \leq -1.25:\\
\;\;\;\;\log \left(\frac{\frac{0.125}{x \cdot x} + \left(-0.5 + \frac{-0.0625}{t\_0}\right)}{x}\right)\\
\mathbf{elif}\;x \leq 1.25:\\
\;\;\;\;x \cdot \frac{1 - t\_0 \cdot \left(0.027777777777777776 + \left(x \cdot x\right) \cdot -0.025\right)}{1 - \left(x \cdot x\right) \cdot \left(-0.16666666666666666 + \left(x \cdot x\right) \cdot \left(0.075 + x \cdot \left(x \cdot -0.044642857142857144\right)\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;\log \left(x \cdot 2 + \frac{0.5}{x}\right)\\
\end{array}
\end{array}
if x < -1.25Initial program 3.4%
log-lowering-log.f64N/A
+-lowering-+.f64N/A
+-commutativeN/A
hypot-1-defN/A
hypot-lowering-hypot.f644.7%
Simplified4.7%
Taylor expanded in x around -inf
associate-*r/N/A
/-lowering-/.f64N/A
Simplified100.0%
if -1.25 < x < 1.25Initial program 6.5%
log-lowering-log.f64N/A
+-lowering-+.f64N/A
+-commutativeN/A
hypot-1-defN/A
hypot-lowering-hypot.f646.5%
Simplified6.5%
Taylor expanded in x around 0
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6499.9%
Simplified99.9%
flip-+N/A
/-lowering-/.f64N/A
Applied egg-rr99.9%
Taylor expanded in x around 0
*-lowering-*.f64N/A
metadata-evalN/A
pow-sqrN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6499.9%
Simplified99.9%
Taylor expanded in x around 0
*-lowering-*.f64N/A
metadata-evalN/A
pow-sqrN/A
unpow2N/A
associate-*l*N/A
unpow2N/A
cube-multN/A
*-lowering-*.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6499.9%
Simplified99.9%
if 1.25 < x Initial program 53.0%
log-lowering-log.f64N/A
+-lowering-+.f64N/A
+-commutativeN/A
hypot-1-defN/A
hypot-lowering-hypot.f64100.0%
Simplified100.0%
Taylor expanded in x around inf
distribute-rgt-inN/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
associate-*r/N/A
metadata-evalN/A
associate-*l/N/A
unpow2N/A
associate-/r*N/A
*-lft-identityN/A
times-fracN/A
metadata-evalN/A
*-inversesN/A
metadata-evalN/A
/-lowering-/.f64100.0%
Simplified100.0%
(FPCore (x)
:precision binary64
(if (<= x -1.35)
(log (/ (+ (/ 0.125 (* x x)) -0.5) x))
(if (<= x 1.25)
(*
x
(/
(-
1.0
(* (* x (* x (* x x))) (+ 0.027777777777777776 (* (* x x) -0.025))))
(-
1.0
(*
(* x x)
(+
-0.16666666666666666
(* (* x x) (+ 0.075 (* x (* x -0.044642857142857144)))))))))
(log (+ (* x 2.0) (/ 0.5 x))))))
double code(double x) {
double tmp;
if (x <= -1.35) {
tmp = log((((0.125 / (x * x)) + -0.5) / x));
} else if (x <= 1.25) {
tmp = x * ((1.0 - ((x * (x * (x * x))) * (0.027777777777777776 + ((x * x) * -0.025)))) / (1.0 - ((x * x) * (-0.16666666666666666 + ((x * x) * (0.075 + (x * (x * -0.044642857142857144))))))));
} else {
tmp = log(((x * 2.0) + (0.5 / x)));
}
return tmp;
}
real(8) function code(x)
real(8), intent (in) :: x
real(8) :: tmp
if (x <= (-1.35d0)) then
tmp = log((((0.125d0 / (x * x)) + (-0.5d0)) / x))
else if (x <= 1.25d0) then
tmp = x * ((1.0d0 - ((x * (x * (x * x))) * (0.027777777777777776d0 + ((x * x) * (-0.025d0))))) / (1.0d0 - ((x * x) * ((-0.16666666666666666d0) + ((x * x) * (0.075d0 + (x * (x * (-0.044642857142857144d0)))))))))
else
tmp = log(((x * 2.0d0) + (0.5d0 / x)))
end if
code = tmp
end function
public static double code(double x) {
double tmp;
if (x <= -1.35) {
tmp = Math.log((((0.125 / (x * x)) + -0.5) / x));
} else if (x <= 1.25) {
tmp = x * ((1.0 - ((x * (x * (x * x))) * (0.027777777777777776 + ((x * x) * -0.025)))) / (1.0 - ((x * x) * (-0.16666666666666666 + ((x * x) * (0.075 + (x * (x * -0.044642857142857144))))))));
} else {
tmp = Math.log(((x * 2.0) + (0.5 / x)));
}
return tmp;
}
def code(x): tmp = 0 if x <= -1.35: tmp = math.log((((0.125 / (x * x)) + -0.5) / x)) elif x <= 1.25: tmp = x * ((1.0 - ((x * (x * (x * x))) * (0.027777777777777776 + ((x * x) * -0.025)))) / (1.0 - ((x * x) * (-0.16666666666666666 + ((x * x) * (0.075 + (x * (x * -0.044642857142857144)))))))) else: tmp = math.log(((x * 2.0) + (0.5 / x))) return tmp
function code(x) tmp = 0.0 if (x <= -1.35) tmp = log(Float64(Float64(Float64(0.125 / Float64(x * x)) + -0.5) / x)); elseif (x <= 1.25) tmp = Float64(x * Float64(Float64(1.0 - Float64(Float64(x * Float64(x * Float64(x * x))) * Float64(0.027777777777777776 + Float64(Float64(x * x) * -0.025)))) / Float64(1.0 - Float64(Float64(x * x) * Float64(-0.16666666666666666 + Float64(Float64(x * x) * Float64(0.075 + Float64(x * Float64(x * -0.044642857142857144))))))))); else tmp = log(Float64(Float64(x * 2.0) + Float64(0.5 / x))); end return tmp end
function tmp_2 = code(x) tmp = 0.0; if (x <= -1.35) tmp = log((((0.125 / (x * x)) + -0.5) / x)); elseif (x <= 1.25) tmp = x * ((1.0 - ((x * (x * (x * x))) * (0.027777777777777776 + ((x * x) * -0.025)))) / (1.0 - ((x * x) * (-0.16666666666666666 + ((x * x) * (0.075 + (x * (x * -0.044642857142857144)))))))); else tmp = log(((x * 2.0) + (0.5 / x))); end tmp_2 = tmp; end
code[x_] := If[LessEqual[x, -1.35], N[Log[N[(N[(N[(0.125 / N[(x * x), $MachinePrecision]), $MachinePrecision] + -0.5), $MachinePrecision] / x), $MachinePrecision]], $MachinePrecision], If[LessEqual[x, 1.25], N[(x * N[(N[(1.0 - N[(N[(x * N[(x * N[(x * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(0.027777777777777776 + N[(N[(x * x), $MachinePrecision] * -0.025), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(1.0 - N[(N[(x * x), $MachinePrecision] * N[(-0.16666666666666666 + N[(N[(x * x), $MachinePrecision] * N[(0.075 + N[(x * N[(x * -0.044642857142857144), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[Log[N[(N[(x * 2.0), $MachinePrecision] + N[(0.5 / x), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq -1.35:\\
\;\;\;\;\log \left(\frac{\frac{0.125}{x \cdot x} + -0.5}{x}\right)\\
\mathbf{elif}\;x \leq 1.25:\\
\;\;\;\;x \cdot \frac{1 - \left(x \cdot \left(x \cdot \left(x \cdot x\right)\right)\right) \cdot \left(0.027777777777777776 + \left(x \cdot x\right) \cdot -0.025\right)}{1 - \left(x \cdot x\right) \cdot \left(-0.16666666666666666 + \left(x \cdot x\right) \cdot \left(0.075 + x \cdot \left(x \cdot -0.044642857142857144\right)\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;\log \left(x \cdot 2 + \frac{0.5}{x}\right)\\
\end{array}
\end{array}
if x < -1.3500000000000001Initial program 3.4%
log-lowering-log.f64N/A
+-lowering-+.f64N/A
+-commutativeN/A
hypot-1-defN/A
hypot-lowering-hypot.f644.7%
Simplified4.7%
Taylor expanded in x around -inf
associate-*r/N/A
mul-1-negN/A
sub-negN/A
+-commutativeN/A
distribute-neg-inN/A
remove-double-negN/A
sub-negN/A
/-lowering-/.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f6499.7%
Simplified99.7%
if -1.3500000000000001 < x < 1.25Initial program 6.5%
log-lowering-log.f64N/A
+-lowering-+.f64N/A
+-commutativeN/A
hypot-1-defN/A
hypot-lowering-hypot.f646.5%
Simplified6.5%
Taylor expanded in x around 0
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6499.9%
Simplified99.9%
flip-+N/A
/-lowering-/.f64N/A
Applied egg-rr99.9%
Taylor expanded in x around 0
*-lowering-*.f64N/A
metadata-evalN/A
pow-sqrN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6499.9%
Simplified99.9%
Taylor expanded in x around 0
*-lowering-*.f64N/A
metadata-evalN/A
pow-sqrN/A
unpow2N/A
associate-*l*N/A
unpow2N/A
cube-multN/A
*-lowering-*.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6499.9%
Simplified99.9%
if 1.25 < x Initial program 53.0%
log-lowering-log.f64N/A
+-lowering-+.f64N/A
+-commutativeN/A
hypot-1-defN/A
hypot-lowering-hypot.f64100.0%
Simplified100.0%
Taylor expanded in x around inf
distribute-rgt-inN/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
associate-*r/N/A
metadata-evalN/A
associate-*l/N/A
unpow2N/A
associate-/r*N/A
*-lft-identityN/A
times-fracN/A
metadata-evalN/A
*-inversesN/A
metadata-evalN/A
/-lowering-/.f64100.0%
Simplified100.0%
Final simplification99.9%
(FPCore (x)
:precision binary64
(if (<= x -1.65)
(log (/ -0.5 x))
(if (<= x 1.25)
(*
x
(/
(-
1.0
(* (* x (* x (* x x))) (+ 0.027777777777777776 (* (* x x) -0.025))))
(-
1.0
(*
(* x x)
(+
-0.16666666666666666
(* (* x x) (+ 0.075 (* x (* x -0.044642857142857144)))))))))
(log (+ (* x 2.0) (/ 0.5 x))))))
double code(double x) {
double tmp;
if (x <= -1.65) {
tmp = log((-0.5 / x));
} else if (x <= 1.25) {
tmp = x * ((1.0 - ((x * (x * (x * x))) * (0.027777777777777776 + ((x * x) * -0.025)))) / (1.0 - ((x * x) * (-0.16666666666666666 + ((x * x) * (0.075 + (x * (x * -0.044642857142857144))))))));
} else {
tmp = log(((x * 2.0) + (0.5 / x)));
}
return tmp;
}
real(8) function code(x)
real(8), intent (in) :: x
real(8) :: tmp
if (x <= (-1.65d0)) then
tmp = log(((-0.5d0) / x))
else if (x <= 1.25d0) then
tmp = x * ((1.0d0 - ((x * (x * (x * x))) * (0.027777777777777776d0 + ((x * x) * (-0.025d0))))) / (1.0d0 - ((x * x) * ((-0.16666666666666666d0) + ((x * x) * (0.075d0 + (x * (x * (-0.044642857142857144d0)))))))))
else
tmp = log(((x * 2.0d0) + (0.5d0 / x)))
end if
code = tmp
end function
public static double code(double x) {
double tmp;
if (x <= -1.65) {
tmp = Math.log((-0.5 / x));
} else if (x <= 1.25) {
tmp = x * ((1.0 - ((x * (x * (x * x))) * (0.027777777777777776 + ((x * x) * -0.025)))) / (1.0 - ((x * x) * (-0.16666666666666666 + ((x * x) * (0.075 + (x * (x * -0.044642857142857144))))))));
} else {
tmp = Math.log(((x * 2.0) + (0.5 / x)));
}
return tmp;
}
def code(x): tmp = 0 if x <= -1.65: tmp = math.log((-0.5 / x)) elif x <= 1.25: tmp = x * ((1.0 - ((x * (x * (x * x))) * (0.027777777777777776 + ((x * x) * -0.025)))) / (1.0 - ((x * x) * (-0.16666666666666666 + ((x * x) * (0.075 + (x * (x * -0.044642857142857144)))))))) else: tmp = math.log(((x * 2.0) + (0.5 / x))) return tmp
function code(x) tmp = 0.0 if (x <= -1.65) tmp = log(Float64(-0.5 / x)); elseif (x <= 1.25) tmp = Float64(x * Float64(Float64(1.0 - Float64(Float64(x * Float64(x * Float64(x * x))) * Float64(0.027777777777777776 + Float64(Float64(x * x) * -0.025)))) / Float64(1.0 - Float64(Float64(x * x) * Float64(-0.16666666666666666 + Float64(Float64(x * x) * Float64(0.075 + Float64(x * Float64(x * -0.044642857142857144))))))))); else tmp = log(Float64(Float64(x * 2.0) + Float64(0.5 / x))); end return tmp end
function tmp_2 = code(x) tmp = 0.0; if (x <= -1.65) tmp = log((-0.5 / x)); elseif (x <= 1.25) tmp = x * ((1.0 - ((x * (x * (x * x))) * (0.027777777777777776 + ((x * x) * -0.025)))) / (1.0 - ((x * x) * (-0.16666666666666666 + ((x * x) * (0.075 + (x * (x * -0.044642857142857144)))))))); else tmp = log(((x * 2.0) + (0.5 / x))); end tmp_2 = tmp; end
code[x_] := If[LessEqual[x, -1.65], N[Log[N[(-0.5 / x), $MachinePrecision]], $MachinePrecision], If[LessEqual[x, 1.25], N[(x * N[(N[(1.0 - N[(N[(x * N[(x * N[(x * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(0.027777777777777776 + N[(N[(x * x), $MachinePrecision] * -0.025), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(1.0 - N[(N[(x * x), $MachinePrecision] * N[(-0.16666666666666666 + N[(N[(x * x), $MachinePrecision] * N[(0.075 + N[(x * N[(x * -0.044642857142857144), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[Log[N[(N[(x * 2.0), $MachinePrecision] + N[(0.5 / x), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq -1.65:\\
\;\;\;\;\log \left(\frac{-0.5}{x}\right)\\
\mathbf{elif}\;x \leq 1.25:\\
\;\;\;\;x \cdot \frac{1 - \left(x \cdot \left(x \cdot \left(x \cdot x\right)\right)\right) \cdot \left(0.027777777777777776 + \left(x \cdot x\right) \cdot -0.025\right)}{1 - \left(x \cdot x\right) \cdot \left(-0.16666666666666666 + \left(x \cdot x\right) \cdot \left(0.075 + x \cdot \left(x \cdot -0.044642857142857144\right)\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;\log \left(x \cdot 2 + \frac{0.5}{x}\right)\\
\end{array}
\end{array}
if x < -1.6499999999999999Initial program 3.4%
log-lowering-log.f64N/A
+-lowering-+.f64N/A
+-commutativeN/A
hypot-1-defN/A
hypot-lowering-hypot.f644.7%
Simplified4.7%
Taylor expanded in x around -inf
/-lowering-/.f6499.3%
Simplified99.3%
if -1.6499999999999999 < x < 1.25Initial program 6.5%
log-lowering-log.f64N/A
+-lowering-+.f64N/A
+-commutativeN/A
hypot-1-defN/A
hypot-lowering-hypot.f646.5%
Simplified6.5%
Taylor expanded in x around 0
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6499.9%
Simplified99.9%
flip-+N/A
/-lowering-/.f64N/A
Applied egg-rr99.9%
Taylor expanded in x around 0
*-lowering-*.f64N/A
metadata-evalN/A
pow-sqrN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6499.9%
Simplified99.9%
Taylor expanded in x around 0
*-lowering-*.f64N/A
metadata-evalN/A
pow-sqrN/A
unpow2N/A
associate-*l*N/A
unpow2N/A
cube-multN/A
*-lowering-*.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6499.9%
Simplified99.9%
if 1.25 < x Initial program 53.0%
log-lowering-log.f64N/A
+-lowering-+.f64N/A
+-commutativeN/A
hypot-1-defN/A
hypot-lowering-hypot.f64100.0%
Simplified100.0%
Taylor expanded in x around inf
distribute-rgt-inN/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
associate-*r/N/A
metadata-evalN/A
associate-*l/N/A
unpow2N/A
associate-/r*N/A
*-lft-identityN/A
times-fracN/A
metadata-evalN/A
*-inversesN/A
metadata-evalN/A
/-lowering-/.f64100.0%
Simplified100.0%
(FPCore (x)
:precision binary64
(if (<= x -1.65)
(log (/ -0.5 x))
(if (<= x 1.65)
(*
x
(/
(-
1.0
(* (* x (* x (* x x))) (+ 0.027777777777777776 (* (* x x) -0.025))))
(-
1.0
(*
(* x x)
(+
-0.16666666666666666
(* (* x x) (+ 0.075 (* x (* x -0.044642857142857144)))))))))
(log (+ x x)))))
double code(double x) {
double tmp;
if (x <= -1.65) {
tmp = log((-0.5 / x));
} else if (x <= 1.65) {
tmp = x * ((1.0 - ((x * (x * (x * x))) * (0.027777777777777776 + ((x * x) * -0.025)))) / (1.0 - ((x * x) * (-0.16666666666666666 + ((x * x) * (0.075 + (x * (x * -0.044642857142857144))))))));
} else {
tmp = log((x + x));
}
return tmp;
}
real(8) function code(x)
real(8), intent (in) :: x
real(8) :: tmp
if (x <= (-1.65d0)) then
tmp = log(((-0.5d0) / x))
else if (x <= 1.65d0) then
tmp = x * ((1.0d0 - ((x * (x * (x * x))) * (0.027777777777777776d0 + ((x * x) * (-0.025d0))))) / (1.0d0 - ((x * x) * ((-0.16666666666666666d0) + ((x * x) * (0.075d0 + (x * (x * (-0.044642857142857144d0)))))))))
else
tmp = log((x + x))
end if
code = tmp
end function
public static double code(double x) {
double tmp;
if (x <= -1.65) {
tmp = Math.log((-0.5 / x));
} else if (x <= 1.65) {
tmp = x * ((1.0 - ((x * (x * (x * x))) * (0.027777777777777776 + ((x * x) * -0.025)))) / (1.0 - ((x * x) * (-0.16666666666666666 + ((x * x) * (0.075 + (x * (x * -0.044642857142857144))))))));
} else {
tmp = Math.log((x + x));
}
return tmp;
}
def code(x): tmp = 0 if x <= -1.65: tmp = math.log((-0.5 / x)) elif x <= 1.65: tmp = x * ((1.0 - ((x * (x * (x * x))) * (0.027777777777777776 + ((x * x) * -0.025)))) / (1.0 - ((x * x) * (-0.16666666666666666 + ((x * x) * (0.075 + (x * (x * -0.044642857142857144)))))))) else: tmp = math.log((x + x)) return tmp
function code(x) tmp = 0.0 if (x <= -1.65) tmp = log(Float64(-0.5 / x)); elseif (x <= 1.65) tmp = Float64(x * Float64(Float64(1.0 - Float64(Float64(x * Float64(x * Float64(x * x))) * Float64(0.027777777777777776 + Float64(Float64(x * x) * -0.025)))) / Float64(1.0 - Float64(Float64(x * x) * Float64(-0.16666666666666666 + Float64(Float64(x * x) * Float64(0.075 + Float64(x * Float64(x * -0.044642857142857144))))))))); else tmp = log(Float64(x + x)); end return tmp end
function tmp_2 = code(x) tmp = 0.0; if (x <= -1.65) tmp = log((-0.5 / x)); elseif (x <= 1.65) tmp = x * ((1.0 - ((x * (x * (x * x))) * (0.027777777777777776 + ((x * x) * -0.025)))) / (1.0 - ((x * x) * (-0.16666666666666666 + ((x * x) * (0.075 + (x * (x * -0.044642857142857144)))))))); else tmp = log((x + x)); end tmp_2 = tmp; end
code[x_] := If[LessEqual[x, -1.65], N[Log[N[(-0.5 / x), $MachinePrecision]], $MachinePrecision], If[LessEqual[x, 1.65], N[(x * N[(N[(1.0 - N[(N[(x * N[(x * N[(x * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(0.027777777777777776 + N[(N[(x * x), $MachinePrecision] * -0.025), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(1.0 - N[(N[(x * x), $MachinePrecision] * N[(-0.16666666666666666 + N[(N[(x * x), $MachinePrecision] * N[(0.075 + N[(x * N[(x * -0.044642857142857144), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[Log[N[(x + x), $MachinePrecision]], $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq -1.65:\\
\;\;\;\;\log \left(\frac{-0.5}{x}\right)\\
\mathbf{elif}\;x \leq 1.65:\\
\;\;\;\;x \cdot \frac{1 - \left(x \cdot \left(x \cdot \left(x \cdot x\right)\right)\right) \cdot \left(0.027777777777777776 + \left(x \cdot x\right) \cdot -0.025\right)}{1 - \left(x \cdot x\right) \cdot \left(-0.16666666666666666 + \left(x \cdot x\right) \cdot \left(0.075 + x \cdot \left(x \cdot -0.044642857142857144\right)\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;\log \left(x + x\right)\\
\end{array}
\end{array}
if x < -1.6499999999999999Initial program 3.4%
log-lowering-log.f64N/A
+-lowering-+.f64N/A
+-commutativeN/A
hypot-1-defN/A
hypot-lowering-hypot.f644.7%
Simplified4.7%
Taylor expanded in x around -inf
/-lowering-/.f6499.3%
Simplified99.3%
if -1.6499999999999999 < x < 1.6499999999999999Initial program 6.5%
log-lowering-log.f64N/A
+-lowering-+.f64N/A
+-commutativeN/A
hypot-1-defN/A
hypot-lowering-hypot.f646.5%
Simplified6.5%
Taylor expanded in x around 0
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6499.9%
Simplified99.9%
flip-+N/A
/-lowering-/.f64N/A
Applied egg-rr99.9%
Taylor expanded in x around 0
*-lowering-*.f64N/A
metadata-evalN/A
pow-sqrN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6499.9%
Simplified99.9%
Taylor expanded in x around 0
*-lowering-*.f64N/A
metadata-evalN/A
pow-sqrN/A
unpow2N/A
associate-*l*N/A
unpow2N/A
cube-multN/A
*-lowering-*.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6499.9%
Simplified99.9%
if 1.6499999999999999 < x Initial program 53.0%
log-lowering-log.f64N/A
+-lowering-+.f64N/A
+-commutativeN/A
hypot-1-defN/A
hypot-lowering-hypot.f64100.0%
Simplified100.0%
Taylor expanded in x around inf
Simplified99.8%
(FPCore (x) :precision binary64 (if (<= x 1.26) x (log (+ x x))))
double code(double x) {
double tmp;
if (x <= 1.26) {
tmp = x;
} else {
tmp = log((x + x));
}
return tmp;
}
real(8) function code(x)
real(8), intent (in) :: x
real(8) :: tmp
if (x <= 1.26d0) then
tmp = x
else
tmp = log((x + x))
end if
code = tmp
end function
public static double code(double x) {
double tmp;
if (x <= 1.26) {
tmp = x;
} else {
tmp = Math.log((x + x));
}
return tmp;
}
def code(x): tmp = 0 if x <= 1.26: tmp = x else: tmp = math.log((x + x)) return tmp
function code(x) tmp = 0.0 if (x <= 1.26) tmp = x; else tmp = log(Float64(x + x)); end return tmp end
function tmp_2 = code(x) tmp = 0.0; if (x <= 1.26) tmp = x; else tmp = log((x + x)); end tmp_2 = tmp; end
code[x_] := If[LessEqual[x, 1.26], x, N[Log[N[(x + x), $MachinePrecision]], $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq 1.26:\\
\;\;\;\;x\\
\mathbf{else}:\\
\;\;\;\;\log \left(x + x\right)\\
\end{array}
\end{array}
if x < 1.26000000000000001Initial program 5.4%
log-lowering-log.f64N/A
+-lowering-+.f64N/A
+-commutativeN/A
hypot-1-defN/A
hypot-lowering-hypot.f645.8%
Simplified5.8%
Taylor expanded in x around 0
Simplified64.4%
if 1.26000000000000001 < x Initial program 53.0%
log-lowering-log.f64N/A
+-lowering-+.f64N/A
+-commutativeN/A
hypot-1-defN/A
hypot-lowering-hypot.f64100.0%
Simplified100.0%
Taylor expanded in x around inf
Simplified99.8%
(FPCore (x) :precision binary64 x)
double code(double x) {
return x;
}
real(8) function code(x)
real(8), intent (in) :: x
code = x
end function
public static double code(double x) {
return x;
}
def code(x): return x
function code(x) return x end
function tmp = code(x) tmp = x; end
code[x_] := x
\begin{array}{l}
\\
x
\end{array}
Initial program 17.6%
log-lowering-log.f64N/A
+-lowering-+.f64N/A
+-commutativeN/A
hypot-1-defN/A
hypot-lowering-hypot.f6430.1%
Simplified30.1%
Taylor expanded in x around 0
Simplified49.2%
(FPCore (x) :precision binary64 (let* ((t_0 (sqrt (+ (* x x) 1.0)))) (if (< x 0.0) (log (/ -1.0 (- x t_0))) (log (+ x t_0)))))
double code(double x) {
double t_0 = sqrt(((x * x) + 1.0));
double tmp;
if (x < 0.0) {
tmp = log((-1.0 / (x - t_0)));
} else {
tmp = log((x + t_0));
}
return tmp;
}
real(8) function code(x)
real(8), intent (in) :: x
real(8) :: t_0
real(8) :: tmp
t_0 = sqrt(((x * x) + 1.0d0))
if (x < 0.0d0) then
tmp = log(((-1.0d0) / (x - t_0)))
else
tmp = log((x + t_0))
end if
code = tmp
end function
public static double code(double x) {
double t_0 = Math.sqrt(((x * x) + 1.0));
double tmp;
if (x < 0.0) {
tmp = Math.log((-1.0 / (x - t_0)));
} else {
tmp = Math.log((x + t_0));
}
return tmp;
}
def code(x): t_0 = math.sqrt(((x * x) + 1.0)) tmp = 0 if x < 0.0: tmp = math.log((-1.0 / (x - t_0))) else: tmp = math.log((x + t_0)) return tmp
function code(x) t_0 = sqrt(Float64(Float64(x * x) + 1.0)) tmp = 0.0 if (x < 0.0) tmp = log(Float64(-1.0 / Float64(x - t_0))); else tmp = log(Float64(x + t_0)); end return tmp end
function tmp_2 = code(x) t_0 = sqrt(((x * x) + 1.0)); tmp = 0.0; if (x < 0.0) tmp = log((-1.0 / (x - t_0))); else tmp = log((x + t_0)); end tmp_2 = tmp; end
code[x_] := Block[{t$95$0 = N[Sqrt[N[(N[(x * x), $MachinePrecision] + 1.0), $MachinePrecision]], $MachinePrecision]}, If[Less[x, 0.0], N[Log[N[(-1.0 / N[(x - t$95$0), $MachinePrecision]), $MachinePrecision]], $MachinePrecision], N[Log[N[(x + t$95$0), $MachinePrecision]], $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \sqrt{x \cdot x + 1}\\
\mathbf{if}\;x < 0:\\
\;\;\;\;\log \left(\frac{-1}{x - t\_0}\right)\\
\mathbf{else}:\\
\;\;\;\;\log \left(x + t\_0\right)\\
\end{array}
\end{array}
herbie shell --seed 2024152
(FPCore (x)
:name "Hyperbolic arcsine"
:precision binary64
:alt
(! :herbie-platform default (if (< x 0) (log (/ -1 (- x (sqrt (+ (* x x) 1))))) (log (+ x (sqrt (+ (* x x) 1))))))
(log (+ x (sqrt (+ (* x x) 1.0)))))