
(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 (- x (/ 0.5 x)))))
double code(double x) {
return log((x + (x - (0.5 / x))));
}
real(8) function code(x)
real(8), intent (in) :: x
code = log((x + (x - (0.5d0 / x))))
end function
public static double code(double x) {
return Math.log((x + (x - (0.5 / x))));
}
def code(x): return math.log((x + (x - (0.5 / x))))
function code(x) return log(Float64(x + Float64(x - Float64(0.5 / x)))) end
function tmp = code(x) tmp = log((x + (x - (0.5 / x)))); end
code[x_] := N[Log[N[(x + N[(x - N[(0.5 / x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]
\begin{array}{l}
\\
\log \left(x + \left(x - \frac{0.5}{x}\right)\right)
\end{array}
Initial program 47.0%
Taylor expanded in x around inf 99.4%
associate-*r/99.4%
metadata-eval99.4%
Simplified99.4%
Final simplification99.4%
(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 47.0%
Taylor expanded in x around inf 99.2%
Final simplification99.2%
(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 47.0%
Taylor expanded in x around inf 99.2%
Applied egg-rr0.0%
Simplified1.6%
Final simplification1.6%
(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 47.0%
Taylor expanded in x around inf 99.2%
flip-+0.0%
difference-of-squares0.0%
associate-*r/0.0%
+-inverses0.0%
+-inverses0.0%
flip-+10.1%
count-210.1%
count-210.1%
swap-sqr10.1%
log-prod10.2%
metadata-eval10.2%
Applied egg-rr10.2%
Simplified13.7%
Applied egg-rr14.6%
Final simplification14.6%
(FPCore (x) :precision binary64 8.0)
double code(double x) {
return 8.0;
}
real(8) function code(x)
real(8), intent (in) :: x
code = 8.0d0
end function
public static double code(double x) {
return 8.0;
}
def code(x): return 8.0
function code(x) return 8.0 end
function tmp = code(x) tmp = 8.0; end
code[x_] := 8.0
\begin{array}{l}
\\
8
\end{array}
Initial program 47.0%
Taylor expanded in x around inf 99.2%
flip-+0.0%
difference-of-squares0.0%
associate-*r/0.0%
+-inverses0.0%
+-inverses0.0%
flip-+10.1%
count-210.1%
count-210.1%
swap-sqr10.1%
log-prod10.2%
metadata-eval10.2%
Applied egg-rr10.2%
Simplified13.7%
Applied egg-rr15.2%
Final simplification15.2%
(FPCore (x) :precision binary64 17.0)
double code(double x) {
return 17.0;
}
real(8) function code(x)
real(8), intent (in) :: x
code = 17.0d0
end function
public static double code(double x) {
return 17.0;
}
def code(x): return 17.0
function code(x) return 17.0 end
function tmp = code(x) tmp = 17.0; end
code[x_] := 17.0
\begin{array}{l}
\\
17
\end{array}
Initial program 47.0%
Taylor expanded in x around inf 99.2%
flip-+0.0%
difference-of-squares0.0%
associate-*r/0.0%
+-inverses0.0%
+-inverses0.0%
flip-+10.1%
count-210.1%
count-210.1%
swap-sqr10.1%
log-prod10.2%
metadata-eval10.2%
Applied egg-rr10.2%
Simplified13.7%
Applied egg-rr15.9%
Final simplification15.9%
(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 47.0%
Taylor expanded in x around inf 99.2%
flip-+0.0%
difference-of-squares0.0%
associate-*r/0.0%
+-inverses0.0%
+-inverses0.0%
flip-+10.1%
count-210.1%
count-210.1%
swap-sqr10.1%
log-prod10.2%
metadata-eval10.2%
Applied egg-rr10.2%
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 47.0%
Taylor expanded in x around inf 99.2%
flip-+0.0%
difference-of-squares0.0%
associate-*r/0.0%
+-inverses0.0%
+-inverses0.0%
flip-+10.1%
count-210.1%
count-210.1%
swap-sqr10.1%
log-prod10.2%
metadata-eval10.2%
Applied egg-rr10.2%
Simplified13.7%
Applied egg-rr17.3%
Final simplification17.3%
(FPCore (x) :precision binary64 65.0)
double code(double x) {
return 65.0;
}
real(8) function code(x)
real(8), intent (in) :: x
code = 65.0d0
end function
public static double code(double x) {
return 65.0;
}
def code(x): return 65.0
function code(x) return 65.0 end
function tmp = code(x) tmp = 65.0; end
code[x_] := 65.0
\begin{array}{l}
\\
65
\end{array}
Initial program 47.0%
Taylor expanded in x around inf 99.2%
flip-+0.0%
difference-of-squares0.0%
associate-*r/0.0%
+-inverses0.0%
+-inverses0.0%
flip-+10.1%
count-210.1%
count-210.1%
swap-sqr10.1%
log-prod10.2%
metadata-eval10.2%
Applied egg-rr10.2%
Simplified13.7%
Applied egg-rr17.3%
Final simplification17.3%
(FPCore (x) :precision binary64 256.0)
double code(double x) {
return 256.0;
}
real(8) function code(x)
real(8), intent (in) :: x
code = 256.0d0
end function
public static double code(double x) {
return 256.0;
}
def code(x): return 256.0
function code(x) return 256.0 end
function tmp = code(x) tmp = 256.0; end
code[x_] := 256.0
\begin{array}{l}
\\
256
\end{array}
Initial program 47.0%
Taylor expanded in x around inf 99.2%
flip-+0.0%
difference-of-squares0.0%
associate-*r/0.0%
+-inverses0.0%
+-inverses0.0%
flip-+10.1%
count-210.1%
count-210.1%
swap-sqr10.1%
log-prod10.2%
metadata-eval10.2%
Applied egg-rr10.2%
Simplified13.7%
Applied egg-rr19.3%
Final simplification19.3%
herbie shell --seed 2023279
(FPCore (x)
:name "Hyperbolic arc-cosine"
:precision binary64
(log (+ x (sqrt (- (* x x) 1.0)))))