
(FPCore (f) :precision binary64 (let* ((t_0 (* (/ PI 4.0) f)) (t_1 (exp t_0)) (t_2 (exp (- t_0)))) (- (* (/ 1.0 (/ PI 4.0)) (log (/ (+ t_1 t_2) (- t_1 t_2)))))))
double code(double f) {
double t_0 = (((double) M_PI) / 4.0) * f;
double t_1 = exp(t_0);
double t_2 = exp(-t_0);
return -((1.0 / (((double) M_PI) / 4.0)) * log(((t_1 + t_2) / (t_1 - t_2))));
}
public static double code(double f) {
double t_0 = (Math.PI / 4.0) * f;
double t_1 = Math.exp(t_0);
double t_2 = Math.exp(-t_0);
return -((1.0 / (Math.PI / 4.0)) * Math.log(((t_1 + t_2) / (t_1 - t_2))));
}
def code(f): t_0 = (math.pi / 4.0) * f t_1 = math.exp(t_0) t_2 = math.exp(-t_0) return -((1.0 / (math.pi / 4.0)) * math.log(((t_1 + t_2) / (t_1 - t_2))))
function code(f) t_0 = Float64(Float64(pi / 4.0) * f) t_1 = exp(t_0) t_2 = exp(Float64(-t_0)) return Float64(-Float64(Float64(1.0 / Float64(pi / 4.0)) * log(Float64(Float64(t_1 + t_2) / Float64(t_1 - t_2))))) end
function tmp = code(f) t_0 = (pi / 4.0) * f; t_1 = exp(t_0); t_2 = exp(-t_0); tmp = -((1.0 / (pi / 4.0)) * log(((t_1 + t_2) / (t_1 - t_2)))); end
code[f_] := Block[{t$95$0 = N[(N[(Pi / 4.0), $MachinePrecision] * f), $MachinePrecision]}, Block[{t$95$1 = N[Exp[t$95$0], $MachinePrecision]}, Block[{t$95$2 = N[Exp[(-t$95$0)], $MachinePrecision]}, (-N[(N[(1.0 / N[(Pi / 4.0), $MachinePrecision]), $MachinePrecision] * N[Log[N[(N[(t$95$1 + t$95$2), $MachinePrecision] / N[(t$95$1 - t$95$2), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision])]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \frac{\pi}{4} \cdot f\\
t_1 := e^{t\_0}\\
t_2 := e^{-t\_0}\\
-\frac{1}{\frac{\pi}{4}} \cdot \log \left(\frac{t\_1 + t\_2}{t\_1 - t\_2}\right)
\end{array}
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 5 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (f) :precision binary64 (let* ((t_0 (* (/ PI 4.0) f)) (t_1 (exp t_0)) (t_2 (exp (- t_0)))) (- (* (/ 1.0 (/ PI 4.0)) (log (/ (+ t_1 t_2) (- t_1 t_2)))))))
double code(double f) {
double t_0 = (((double) M_PI) / 4.0) * f;
double t_1 = exp(t_0);
double t_2 = exp(-t_0);
return -((1.0 / (((double) M_PI) / 4.0)) * log(((t_1 + t_2) / (t_1 - t_2))));
}
public static double code(double f) {
double t_0 = (Math.PI / 4.0) * f;
double t_1 = Math.exp(t_0);
double t_2 = Math.exp(-t_0);
return -((1.0 / (Math.PI / 4.0)) * Math.log(((t_1 + t_2) / (t_1 - t_2))));
}
def code(f): t_0 = (math.pi / 4.0) * f t_1 = math.exp(t_0) t_2 = math.exp(-t_0) return -((1.0 / (math.pi / 4.0)) * math.log(((t_1 + t_2) / (t_1 - t_2))))
function code(f) t_0 = Float64(Float64(pi / 4.0) * f) t_1 = exp(t_0) t_2 = exp(Float64(-t_0)) return Float64(-Float64(Float64(1.0 / Float64(pi / 4.0)) * log(Float64(Float64(t_1 + t_2) / Float64(t_1 - t_2))))) end
function tmp = code(f) t_0 = (pi / 4.0) * f; t_1 = exp(t_0); t_2 = exp(-t_0); tmp = -((1.0 / (pi / 4.0)) * log(((t_1 + t_2) / (t_1 - t_2)))); end
code[f_] := Block[{t$95$0 = N[(N[(Pi / 4.0), $MachinePrecision] * f), $MachinePrecision]}, Block[{t$95$1 = N[Exp[t$95$0], $MachinePrecision]}, Block[{t$95$2 = N[Exp[(-t$95$0)], $MachinePrecision]}, (-N[(N[(1.0 / N[(Pi / 4.0), $MachinePrecision]), $MachinePrecision] * N[Log[N[(N[(t$95$1 + t$95$2), $MachinePrecision] / N[(t$95$1 - t$95$2), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision])]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \frac{\pi}{4} \cdot f\\
t_1 := e^{t\_0}\\
t_2 := e^{-t\_0}\\
-\frac{1}{\frac{\pi}{4}} \cdot \log \left(\frac{t\_1 + t\_2}{t\_1 - t\_2}\right)
\end{array}
\end{array}
(FPCore (f) :precision binary64 (/ (log (tanh (/ f (/ 4.0 PI)))) (/ PI 4.0)))
double code(double f) {
return log(tanh((f / (4.0 / ((double) M_PI))))) / (((double) M_PI) / 4.0);
}
public static double code(double f) {
return Math.log(Math.tanh((f / (4.0 / Math.PI)))) / (Math.PI / 4.0);
}
def code(f): return math.log(math.tanh((f / (4.0 / math.pi)))) / (math.pi / 4.0)
function code(f) return Float64(log(tanh(Float64(f / Float64(4.0 / pi)))) / Float64(pi / 4.0)) end
function tmp = code(f) tmp = log(tanh((f / (4.0 / pi)))) / (pi / 4.0); end
code[f_] := N[(N[Log[N[Tanh[N[(f / N[(4.0 / Pi), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]], $MachinePrecision] / N[(Pi / 4.0), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{\log \tanh \left(\frac{f}{\frac{4}{\pi}}\right)}{\frac{\pi}{4}}
\end{array}
Initial program 5.6%
distribute-rgt-neg-inN/A
clear-numN/A
distribute-rgt-neg-inN/A
*-commutativeN/A
clear-numN/A
Applied egg-rr99.5%
(FPCore (f) :precision binary64 (/ (/ (log (/ (+ (/ 4.0 PI) (* (* f 0.08333333333333333) (* f PI))) f)) PI) -0.25))
double code(double f) {
return (log((((4.0 / ((double) M_PI)) + ((f * 0.08333333333333333) * (f * ((double) M_PI)))) / f)) / ((double) M_PI)) / -0.25;
}
public static double code(double f) {
return (Math.log((((4.0 / Math.PI) + ((f * 0.08333333333333333) * (f * Math.PI))) / f)) / Math.PI) / -0.25;
}
def code(f): return (math.log((((4.0 / math.pi) + ((f * 0.08333333333333333) * (f * math.pi))) / f)) / math.pi) / -0.25
function code(f) return Float64(Float64(log(Float64(Float64(Float64(4.0 / pi) + Float64(Float64(f * 0.08333333333333333) * Float64(f * pi))) / f)) / pi) / -0.25) end
function tmp = code(f) tmp = (log((((4.0 / pi) + ((f * 0.08333333333333333) * (f * pi))) / f)) / pi) / -0.25; end
code[f_] := N[(N[(N[Log[N[(N[(N[(4.0 / Pi), $MachinePrecision] + N[(N[(f * 0.08333333333333333), $MachinePrecision] * N[(f * Pi), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / f), $MachinePrecision]], $MachinePrecision] / Pi), $MachinePrecision] / -0.25), $MachinePrecision]
\begin{array}{l}
\\
\frac{\frac{\log \left(\frac{\frac{4}{\pi} + \left(f \cdot 0.08333333333333333\right) \cdot \left(f \cdot \pi\right)}{f}\right)}{\pi}}{-0.25}
\end{array}
Initial program 5.6%
distribute-lft-neg-inN/A
distribute-neg-frac2N/A
associate-*l/N/A
*-lft-identityN/A
/-lowering-/.f64N/A
Simplified5.6%
Taylor expanded in f around 0
Simplified97.0%
Taylor expanded in f around 0
distribute-rgt-outN/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
PI-lowering-PI.f64N/A
metadata-eval97.0%
Simplified97.0%
div-invN/A
metadata-evalN/A
associate-/r*N/A
/-lowering-/.f64N/A
Applied egg-rr97.0%
Final simplification97.0%
(FPCore (f) :precision binary64 (/ (/ (log (/ 4.0 (* f PI))) PI) -0.25))
double code(double f) {
return (log((4.0 / (f * ((double) M_PI)))) / ((double) M_PI)) / -0.25;
}
public static double code(double f) {
return (Math.log((4.0 / (f * Math.PI))) / Math.PI) / -0.25;
}
def code(f): return (math.log((4.0 / (f * math.pi))) / math.pi) / -0.25
function code(f) return Float64(Float64(log(Float64(4.0 / Float64(f * pi))) / pi) / -0.25) end
function tmp = code(f) tmp = (log((4.0 / (f * pi))) / pi) / -0.25; end
code[f_] := N[(N[(N[Log[N[(4.0 / N[(f * Pi), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] / Pi), $MachinePrecision] / -0.25), $MachinePrecision]
\begin{array}{l}
\\
\frac{\frac{\log \left(\frac{4}{f \cdot \pi}\right)}{\pi}}{-0.25}
\end{array}
Initial program 5.6%
distribute-lft-neg-inN/A
distribute-neg-frac2N/A
associate-*l/N/A
*-lft-identityN/A
/-lowering-/.f64N/A
Simplified5.6%
Taylor expanded in f around 0
/-lowering-/.f64N/A
*-commutativeN/A
distribute-rgt-out--N/A
metadata-evalN/A
associate-*l*N/A
metadata-evalN/A
*-inversesN/A
metadata-evalN/A
times-fracN/A
*-commutativeN/A
metadata-evalN/A
distribute-rgt-out--N/A
associate-*r/N/A
*-lowering-*.f64N/A
PI-lowering-PI.f64N/A
Simplified96.6%
div-invN/A
metadata-evalN/A
associate-/r*N/A
/-lowering-/.f64N/A
Applied egg-rr96.6%
(FPCore (f) :precision binary64 (/ -4.0 (/ PI (log (/ (/ 4.0 PI) f)))))
double code(double f) {
return -4.0 / (((double) M_PI) / log(((4.0 / ((double) M_PI)) / f)));
}
public static double code(double f) {
return -4.0 / (Math.PI / Math.log(((4.0 / Math.PI) / f)));
}
def code(f): return -4.0 / (math.pi / math.log(((4.0 / math.pi) / f)))
function code(f) return Float64(-4.0 / Float64(pi / log(Float64(Float64(4.0 / pi) / f)))) end
function tmp = code(f) tmp = -4.0 / (pi / log(((4.0 / pi) / f))); end
code[f_] := N[(-4.0 / N[(Pi / N[Log[N[(N[(4.0 / Pi), $MachinePrecision] / f), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{-4}{\frac{\pi}{\log \left(\frac{\frac{4}{\pi}}{f}\right)}}
\end{array}
Initial program 5.6%
distribute-lft-neg-inN/A
distribute-neg-frac2N/A
associate-*l/N/A
*-lft-identityN/A
/-lowering-/.f64N/A
Simplified5.6%
Taylor expanded in f around 0
/-lowering-/.f64N/A
*-commutativeN/A
distribute-rgt-out--N/A
metadata-evalN/A
associate-*l*N/A
metadata-evalN/A
*-inversesN/A
metadata-evalN/A
times-fracN/A
*-commutativeN/A
metadata-evalN/A
distribute-rgt-out--N/A
associate-*r/N/A
*-lowering-*.f64N/A
PI-lowering-PI.f64N/A
Simplified96.6%
div-invN/A
metadata-evalN/A
associate-/r*N/A
/-lowering-/.f64N/A
Applied egg-rr96.6%
div-invN/A
clear-numN/A
metadata-evalN/A
associate-*l/N/A
metadata-evalN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
PI-lowering-PI.f64N/A
log-lowering-log.f64N/A
*-commutativeN/A
associate-/r*N/A
metadata-evalN/A
associate-/r*N/A
*-commutativeN/A
/-lowering-/.f64N/A
*-commutativeN/A
associate-/r*N/A
metadata-evalN/A
/-lowering-/.f64N/A
PI-lowering-PI.f6496.4%
Applied egg-rr96.4%
(FPCore (f) :precision binary64 (* (log (/ 4.0 (* f PI))) (/ -4.0 PI)))
double code(double f) {
return log((4.0 / (f * ((double) M_PI)))) * (-4.0 / ((double) M_PI));
}
public static double code(double f) {
return Math.log((4.0 / (f * Math.PI))) * (-4.0 / Math.PI);
}
def code(f): return math.log((4.0 / (f * math.pi))) * (-4.0 / math.pi)
function code(f) return Float64(log(Float64(4.0 / Float64(f * pi))) * Float64(-4.0 / pi)) end
function tmp = code(f) tmp = log((4.0 / (f * pi))) * (-4.0 / pi); end
code[f_] := N[(N[Log[N[(4.0 / N[(f * Pi), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * N[(-4.0 / Pi), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\log \left(\frac{4}{f \cdot \pi}\right) \cdot \frac{-4}{\pi}
\end{array}
Initial program 5.6%
distribute-lft-neg-inN/A
distribute-neg-frac2N/A
associate-*l/N/A
*-lft-identityN/A
/-lowering-/.f64N/A
Simplified5.6%
Taylor expanded in f around 0
/-lowering-/.f64N/A
*-commutativeN/A
distribute-rgt-out--N/A
metadata-evalN/A
associate-*l*N/A
metadata-evalN/A
*-inversesN/A
metadata-evalN/A
times-fracN/A
*-commutativeN/A
metadata-evalN/A
distribute-rgt-out--N/A
associate-*r/N/A
*-lowering-*.f64N/A
PI-lowering-PI.f64N/A
Simplified96.6%
div-invN/A
clear-numN/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
PI-lowering-PI.f64N/A
associate-/r*N/A
div-invN/A
*-commutativeN/A
times-fracN/A
metadata-evalN/A
associate-/l*N/A
div-invN/A
log-lowering-log.f64N/A
associate-/l/N/A
*-commutativeN/A
/-lowering-/.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
PI-lowering-PI.f6496.4%
Applied egg-rr96.4%
Final simplification96.4%
herbie shell --seed 2024170
(FPCore (f)
:name "VandenBroeck and Keller, Equation (20)"
:precision binary64
(- (* (/ 1.0 (/ PI 4.0)) (log (/ (+ (exp (* (/ PI 4.0) f)) (exp (- (* (/ PI 4.0) f)))) (- (exp (* (/ PI 4.0) f)) (exp (- (* (/ PI 4.0) f)))))))))