
(FPCore (x) :precision binary64 (acosh x))
double code(double x) {
return acosh(x);
}
def code(x): return math.acosh(x)
function code(x) return acosh(x) end
function tmp = code(x) tmp = acosh(x); end
code[x_] := N[ArcCosh[x], $MachinePrecision]
\begin{array}{l}
\\
\cosh^{-1} x
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 12 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 59.1%
Taylor expanded in x around inf 99.3%
associate-*r/99.3%
metadata-eval99.3%
Simplified99.3%
Final simplification99.3%
(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 59.1%
Taylor expanded in x around inf 99.0%
Final simplification99.0%
(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 59.1%
Taylor expanded in x around inf 99.0%
Applied egg-rr0.0%
Simplified1.6%
Final simplification1.6%
(FPCore (x) :precision binary64 0.6666666666666666)
double code(double x) {
return 0.6666666666666666;
}
real(8) function code(x)
real(8), intent (in) :: x
code = 0.6666666666666666d0
end function
public static double code(double x) {
return 0.6666666666666666;
}
def code(x): return 0.6666666666666666
function code(x) return 0.6666666666666666 end
function tmp = code(x) tmp = 0.6666666666666666; end
code[x_] := 0.6666666666666666
\begin{array}{l}
\\
0.6666666666666666
\end{array}
Initial program 59.1%
Taylor expanded in x around inf 99.0%
flip-+0.0%
difference-of-squares0.0%
associate-*r/0.0%
+-inverses0.0%
+-inverses0.0%
flip-+12.1%
count-212.1%
count-212.1%
swap-sqr12.1%
log-prod12.2%
metadata-eval12.2%
Applied egg-rr12.2%
Simplified13.8%
Applied egg-rr14.0%
Final simplification14.0%
(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 59.1%
Taylor expanded in x around inf 99.0%
flip-+0.0%
difference-of-squares0.0%
associate-*r/0.0%
+-inverses0.0%
+-inverses0.0%
flip-+12.1%
count-212.1%
count-212.1%
swap-sqr12.1%
log-prod12.2%
metadata-eval12.2%
Applied egg-rr12.2%
Simplified13.8%
Applied egg-rr15.1%
Final simplification15.1%
(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 59.1%
Taylor expanded in x around inf 99.0%
flip-+0.0%
difference-of-squares0.0%
associate-*r/0.0%
+-inverses0.0%
+-inverses0.0%
flip-+12.1%
count-212.1%
count-212.1%
swap-sqr12.1%
log-prod12.2%
metadata-eval12.2%
Applied egg-rr12.2%
Simplified13.8%
Applied egg-rr15.4%
Final simplification15.4%
(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 59.1%
Taylor expanded in x around inf 99.0%
flip-+0.0%
difference-of-squares0.0%
associate-*r/0.0%
+-inverses0.0%
+-inverses0.0%
flip-+12.1%
count-212.1%
count-212.1%
swap-sqr12.1%
log-prod12.2%
metadata-eval12.2%
Applied egg-rr12.2%
Simplified13.8%
Applied egg-rr15.9%
Final simplification15.9%
(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 59.1%
Taylor expanded in x around inf 99.0%
flip-+0.0%
difference-of-squares0.0%
associate-*r/0.0%
+-inverses0.0%
+-inverses0.0%
flip-+12.1%
count-212.1%
count-212.1%
swap-sqr12.1%
log-prod12.2%
metadata-eval12.2%
Applied egg-rr12.2%
Simplified13.8%
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 59.1%
Taylor expanded in x around inf 99.0%
flip-+0.0%
difference-of-squares0.0%
associate-*r/0.0%
+-inverses0.0%
+-inverses0.0%
flip-+12.1%
count-212.1%
count-212.1%
swap-sqr12.1%
log-prod12.2%
metadata-eval12.2%
Applied egg-rr12.2%
Simplified13.8%
Applied egg-rr16.4%
Final simplification16.4%
(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 59.1%
Taylor expanded in x around inf 99.0%
flip-+0.0%
difference-of-squares0.0%
associate-*r/0.0%
+-inverses0.0%
+-inverses0.0%
flip-+12.1%
count-212.1%
count-212.1%
swap-sqr12.1%
log-prod12.2%
metadata-eval12.2%
Applied egg-rr12.2%
Simplified13.8%
Applied egg-rr17.5%
Final simplification17.5%
(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 59.1%
Taylor expanded in x around inf 99.0%
flip-+0.0%
difference-of-squares0.0%
associate-*r/0.0%
+-inverses0.0%
+-inverses0.0%
flip-+12.1%
count-212.1%
count-212.1%
swap-sqr12.1%
log-prod12.2%
metadata-eval12.2%
Applied egg-rr12.2%
Simplified13.8%
Applied egg-rr17.6%
Final simplification17.6%
(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 59.1%
Taylor expanded in x around inf 99.0%
flip-+0.0%
difference-of-squares0.0%
associate-*r/0.0%
+-inverses0.0%
+-inverses0.0%
flip-+12.1%
count-212.1%
count-212.1%
swap-sqr12.1%
log-prod12.2%
metadata-eval12.2%
Applied egg-rr12.2%
Simplified13.8%
Applied egg-rr19.6%
Final simplification19.6%
(FPCore (x) :precision binary64 (log (+ x (* (sqrt (- x 1.0)) (sqrt (+ x 1.0))))))
double code(double x) {
return log((x + (sqrt((x - 1.0)) * sqrt((x + 1.0)))));
}
real(8) function code(x)
real(8), intent (in) :: x
code = log((x + (sqrt((x - 1.0d0)) * sqrt((x + 1.0d0)))))
end function
public static double code(double x) {
return Math.log((x + (Math.sqrt((x - 1.0)) * Math.sqrt((x + 1.0)))));
}
def code(x): return math.log((x + (math.sqrt((x - 1.0)) * math.sqrt((x + 1.0)))))
function code(x) return log(Float64(x + Float64(sqrt(Float64(x - 1.0)) * sqrt(Float64(x + 1.0))))) end
function tmp = code(x) tmp = log((x + (sqrt((x - 1.0)) * sqrt((x + 1.0))))); end
code[x_] := N[Log[N[(x + N[(N[Sqrt[N[(x - 1.0), $MachinePrecision]], $MachinePrecision] * N[Sqrt[N[(x + 1.0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]
\begin{array}{l}
\\
\log \left(x + \sqrt{x - 1} \cdot \sqrt{x + 1}\right)
\end{array}
herbie shell --seed 2023275
(FPCore (x)
:name "Rust f64::acosh"
:precision binary64
:pre (>= x 1.0)
:herbie-target
(log (+ x (* (sqrt (- x 1.0)) (sqrt (+ x 1.0)))))
(log (+ x (sqrt (- (* x x) 1.0)))))