
(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 10 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 (log (- (* x 2.0) (/ 0.5 x))))
double code(double x) {
return log(((x * 2.0) - (0.5 / x)));
}
real(8) function code(x)
real(8), intent (in) :: x
code = log(((x * 2.0d0) - (0.5d0 / x)))
end function
public static double code(double x) {
return Math.log(((x * 2.0) - (0.5 / x)));
}
def code(x): return math.log(((x * 2.0) - (0.5 / x)))
function code(x) return log(Float64(Float64(x * 2.0) - Float64(0.5 / x))) end
function tmp = code(x) tmp = log(((x * 2.0) - (0.5 / x))); end
code[x_] := N[Log[N[(N[(x * 2.0), $MachinePrecision] - N[(0.5 / x), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]
\begin{array}{l}
\\
\log \left(x \cdot 2 - \frac{0.5}{x}\right)
\end{array}
Initial program 49.7%
Taylor expanded in x around inf 99.1%
*-commutative99.1%
associate-*r/99.1%
metadata-eval99.1%
Simplified99.1%
Final simplification99.1%
(FPCore (x) :precision binary64 (log (+ x x)))
double code(double x) {
return log((x + x));
}
real(8) function code(x)
real(8), intent (in) :: x
code = log((x + x))
end function
public static double code(double x) {
return Math.log((x + x));
}
def code(x): return math.log((x + x))
function code(x) return log(Float64(x + x)) end
function tmp = code(x) tmp = log((x + x)); end
code[x_] := N[Log[N[(x + x), $MachinePrecision]], $MachinePrecision]
\begin{array}{l}
\\
\log \left(x + x\right)
\end{array}
Initial program 49.7%
Taylor expanded in x around inf 98.8%
Final simplification98.8%
(FPCore (x) :precision binary64 -1.0)
double code(double x) {
return -1.0;
}
real(8) function code(x)
real(8), intent (in) :: x
code = -1.0d0
end function
public static double code(double x) {
return -1.0;
}
def code(x): return -1.0
function code(x) return -1.0 end
function tmp = code(x) tmp = -1.0; end
code[x_] := -1.0
\begin{array}{l}
\\
-1
\end{array}
Initial program 49.7%
Taylor expanded in x around inf 98.8%
flip-+0.0%
difference-of-squares0.0%
count-20.0%
*-commutative0.0%
associate-*r/0.0%
+-inverses0.0%
+-inverses0.0%
flip-+10.6%
count-210.6%
*-commutative10.6%
sum-log18.7%
*-commutative18.7%
count-218.7%
flip-+0.0%
+-inverses0.0%
+-inverses0.0%
*-commutative0.0%
count-20.0%
flip-+0.0%
+-inverses0.0%
+-inverses0.0%
Applied egg-rr0.0%
Simplified1.6%
Final simplification1.6%
(FPCore (x) :precision binary64 1.5)
double code(double x) {
return 1.5;
}
real(8) function code(x)
real(8), intent (in) :: x
code = 1.5d0
end function
public static double code(double x) {
return 1.5;
}
def code(x): return 1.5
function code(x) return 1.5 end
function tmp = code(x) tmp = 1.5; end
code[x_] := 1.5
\begin{array}{l}
\\
1.5
\end{array}
Initial program 49.7%
Taylor expanded in x around inf 98.8%
flip-+0.0%
difference-of-squares0.0%
count-20.0%
*-commutative0.0%
associate-*r/0.0%
+-inverses0.0%
+-inverses0.0%
flip-+10.6%
count-210.6%
*-commutative10.6%
swap-sqr10.6%
rem-exp-log10.6%
rem-exp-log10.6%
log-prod10.6%
pow210.6%
rem-exp-log10.6%
rem-exp-log10.6%
metadata-eval10.6%
Applied egg-rr10.6%
Simplified13.7%
Applied egg-rr14.3%
Final simplification14.3%
(FPCore (x) :precision binary64 3.0)
double code(double x) {
return 3.0;
}
real(8) function code(x)
real(8), intent (in) :: x
code = 3.0d0
end function
public static double code(double x) {
return 3.0;
}
def code(x): return 3.0
function code(x) return 3.0 end
function tmp = code(x) tmp = 3.0; end
code[x_] := 3.0
\begin{array}{l}
\\
3
\end{array}
Initial program 49.7%
Taylor expanded in x around inf 98.8%
flip-+0.0%
difference-of-squares0.0%
count-20.0%
*-commutative0.0%
associate-*r/0.0%
+-inverses0.0%
+-inverses0.0%
flip-+10.6%
count-210.6%
*-commutative10.6%
swap-sqr10.6%
rem-exp-log10.6%
rem-exp-log10.6%
log-prod10.6%
pow210.6%
rem-exp-log10.6%
rem-exp-log10.6%
metadata-eval10.6%
Applied egg-rr10.6%
Simplified13.7%
Applied egg-rr14.6%
Final simplification14.6%
(FPCore (x) :precision binary64 6.0)
double code(double x) {
return 6.0;
}
real(8) function code(x)
real(8), intent (in) :: x
code = 6.0d0
end function
public static double code(double x) {
return 6.0;
}
def code(x): return 6.0
function code(x) return 6.0 end
function tmp = code(x) tmp = 6.0; end
code[x_] := 6.0
\begin{array}{l}
\\
6
\end{array}
Initial program 49.7%
Taylor expanded in x around inf 98.8%
flip-+0.0%
difference-of-squares0.0%
count-20.0%
*-commutative0.0%
associate-*r/0.0%
+-inverses0.0%
+-inverses0.0%
flip-+10.6%
count-210.6%
*-commutative10.6%
swap-sqr10.6%
rem-exp-log10.6%
rem-exp-log10.6%
log-prod10.6%
pow210.6%
rem-exp-log10.6%
rem-exp-log10.6%
metadata-eval10.6%
Applied egg-rr10.6%
Simplified13.7%
Applied egg-rr15.0%
Final simplification15.0%
(FPCore (x) :precision binary64 9.0)
double code(double x) {
return 9.0;
}
real(8) function code(x)
real(8), intent (in) :: x
code = 9.0d0
end function
public static double code(double x) {
return 9.0;
}
def code(x): return 9.0
function code(x) return 9.0 end
function tmp = code(x) tmp = 9.0; end
code[x_] := 9.0
\begin{array}{l}
\\
9
\end{array}
Initial program 49.7%
Taylor expanded in x around inf 98.8%
flip-+0.0%
difference-of-squares0.0%
count-20.0%
*-commutative0.0%
associate-*r/0.0%
+-inverses0.0%
+-inverses0.0%
flip-+10.6%
count-210.6%
*-commutative10.6%
swap-sqr10.6%
rem-exp-log10.6%
rem-exp-log10.6%
log-prod10.6%
pow210.6%
rem-exp-log10.6%
rem-exp-log10.6%
metadata-eval10.6%
Applied egg-rr10.6%
Simplified13.7%
Applied egg-rr15.3%
Final simplification15.3%
(FPCore (x) :precision binary64 16.0)
double code(double x) {
return 16.0;
}
real(8) function code(x)
real(8), intent (in) :: x
code = 16.0d0
end function
public static double code(double x) {
return 16.0;
}
def code(x): return 16.0
function code(x) return 16.0 end
function tmp = code(x) tmp = 16.0; end
code[x_] := 16.0
\begin{array}{l}
\\
16
\end{array}
Initial program 49.7%
Taylor expanded in x around inf 98.8%
flip-+0.0%
difference-of-squares0.0%
count-20.0%
*-commutative0.0%
associate-*r/0.0%
+-inverses0.0%
+-inverses0.0%
flip-+10.6%
count-210.6%
*-commutative10.6%
swap-sqr10.6%
rem-exp-log10.6%
rem-exp-log10.6%
log-prod10.6%
pow210.6%
rem-exp-log10.6%
rem-exp-log10.6%
metadata-eval10.6%
Applied egg-rr10.6%
Simplified13.7%
Applied egg-rr15.8%
Final simplification15.8%
(FPCore (x) :precision binary64 27.0)
double code(double x) {
return 27.0;
}
real(8) function code(x)
real(8), intent (in) :: x
code = 27.0d0
end function
public static double code(double x) {
return 27.0;
}
def code(x): return 27.0
function code(x) return 27.0 end
function tmp = code(x) tmp = 27.0; end
code[x_] := 27.0
\begin{array}{l}
\\
27
\end{array}
Initial program 49.7%
Taylor expanded in x around inf 98.8%
flip-+0.0%
difference-of-squares0.0%
count-20.0%
*-commutative0.0%
associate-*r/0.0%
+-inverses0.0%
+-inverses0.0%
flip-+10.6%
count-210.6%
*-commutative10.6%
swap-sqr10.6%
rem-exp-log10.6%
rem-exp-log10.6%
log-prod10.6%
pow210.6%
rem-exp-log10.6%
rem-exp-log10.6%
metadata-eval10.6%
Applied egg-rr10.6%
Simplified13.7%
Applied egg-rr16.3%
Final simplification16.3%
(FPCore (x) :precision binary64 64.0)
double code(double x) {
return 64.0;
}
real(8) function code(x)
real(8), intent (in) :: x
code = 64.0d0
end function
public static double code(double x) {
return 64.0;
}
def code(x): return 64.0
function code(x) return 64.0 end
function tmp = code(x) tmp = 64.0; end
code[x_] := 64.0
\begin{array}{l}
\\
64
\end{array}
Initial program 49.7%
Taylor expanded in x around inf 98.8%
flip-+0.0%
difference-of-squares0.0%
count-20.0%
*-commutative0.0%
associate-*r/0.0%
+-inverses0.0%
+-inverses0.0%
flip-+10.6%
count-210.6%
*-commutative10.6%
swap-sqr10.6%
rem-exp-log10.6%
rem-exp-log10.6%
log-prod10.6%
pow210.6%
rem-exp-log10.6%
rem-exp-log10.6%
metadata-eval10.6%
Applied egg-rr10.6%
Simplified13.7%
Applied egg-rr17.4%
Final simplification17.4%
herbie shell --seed 2023325
(FPCore (x)
:name "Hyperbolic arc-cosine"
:precision binary64
(log (+ x (sqrt (- (* x x) 1.0)))))