
(FPCore (x y) :precision binary64 (let* ((t_0 (/ x (* y 2.0)))) (/ (tan t_0) (sin t_0))))
double code(double x, double y) {
double t_0 = x / (y * 2.0);
return tan(t_0) / sin(t_0);
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8) :: t_0
t_0 = x / (y * 2.0d0)
code = tan(t_0) / sin(t_0)
end function
public static double code(double x, double y) {
double t_0 = x / (y * 2.0);
return Math.tan(t_0) / Math.sin(t_0);
}
def code(x, y): t_0 = x / (y * 2.0) return math.tan(t_0) / math.sin(t_0)
function code(x, y) t_0 = Float64(x / Float64(y * 2.0)) return Float64(tan(t_0) / sin(t_0)) end
function tmp = code(x, y) t_0 = x / (y * 2.0); tmp = tan(t_0) / sin(t_0); end
code[x_, y_] := Block[{t$95$0 = N[(x / N[(y * 2.0), $MachinePrecision]), $MachinePrecision]}, N[(N[Tan[t$95$0], $MachinePrecision] / N[Sin[t$95$0], $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \frac{x}{y \cdot 2}\\
\frac{\tan t\_0}{\sin t\_0}
\end{array}
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 6 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (x y) :precision binary64 (let* ((t_0 (/ x (* y 2.0)))) (/ (tan t_0) (sin t_0))))
double code(double x, double y) {
double t_0 = x / (y * 2.0);
return tan(t_0) / sin(t_0);
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8) :: t_0
t_0 = x / (y * 2.0d0)
code = tan(t_0) / sin(t_0)
end function
public static double code(double x, double y) {
double t_0 = x / (y * 2.0);
return Math.tan(t_0) / Math.sin(t_0);
}
def code(x, y): t_0 = x / (y * 2.0) return math.tan(t_0) / math.sin(t_0)
function code(x, y) t_0 = Float64(x / Float64(y * 2.0)) return Float64(tan(t_0) / sin(t_0)) end
function tmp = code(x, y) t_0 = x / (y * 2.0); tmp = tan(t_0) / sin(t_0); end
code[x_, y_] := Block[{t$95$0 = N[(x / N[(y * 2.0), $MachinePrecision]), $MachinePrecision]}, N[(N[Tan[t$95$0], $MachinePrecision] / N[Sin[t$95$0], $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \frac{x}{y \cdot 2}\\
\frac{\tan t\_0}{\sin t\_0}
\end{array}
\end{array}
y_m = (fabs.f64 y) (FPCore (x y_m) :precision binary64 (/ 1.0 (cbrt (pow (cos (/ (pow y_m -0.5) (* (sqrt y_m) (/ -2.0 x)))) 3.0))))
y_m = fabs(y);
double code(double x, double y_m) {
return 1.0 / cbrt(pow(cos((pow(y_m, -0.5) / (sqrt(y_m) * (-2.0 / x)))), 3.0));
}
y_m = Math.abs(y);
public static double code(double x, double y_m) {
return 1.0 / Math.cbrt(Math.pow(Math.cos((Math.pow(y_m, -0.5) / (Math.sqrt(y_m) * (-2.0 / x)))), 3.0));
}
y_m = abs(y) function code(x, y_m) return Float64(1.0 / cbrt((cos(Float64((y_m ^ -0.5) / Float64(sqrt(y_m) * Float64(-2.0 / x)))) ^ 3.0))) end
y_m = N[Abs[y], $MachinePrecision] code[x_, y$95$m_] := N[(1.0 / N[Power[N[Power[N[Cos[N[(N[Power[y$95$m, -0.5], $MachinePrecision] / N[(N[Sqrt[y$95$m], $MachinePrecision] * N[(-2.0 / x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision], 3.0], $MachinePrecision], 1/3], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
y_m = \left|y\right|
\\
\frac{1}{\sqrt[3]{{\cos \left(\frac{{y\_m}^{-0.5}}{\sqrt{y\_m} \cdot \frac{-2}{x}}\right)}^{3}}}
\end{array}
Initial program 40.4%
remove-double-neg40.4%
distribute-frac-neg40.4%
tan-neg40.4%
distribute-frac-neg240.4%
distribute-lft-neg-out40.4%
distribute-frac-neg240.4%
distribute-lft-neg-out40.4%
distribute-frac-neg240.4%
distribute-frac-neg40.4%
neg-mul-140.4%
*-commutative40.4%
associate-/l*40.5%
*-commutative40.5%
associate-/r*40.5%
metadata-eval40.5%
sin-neg40.5%
distribute-frac-neg40.5%
Simplified40.6%
Taylor expanded in x around inf 54.5%
associate-*r/54.5%
Simplified54.5%
Taylor expanded in x around inf 54.5%
associate-*r/54.5%
*-commutative54.5%
associate-/l*54.7%
Simplified54.7%
add-cbrt-cube54.7%
pow354.7%
Applied egg-rr54.7%
add-sqr-sqrt27.9%
sqrt-unprod50.8%
frac-times51.0%
metadata-eval51.0%
metadata-eval51.0%
frac-times50.8%
sqrt-unprod26.8%
add-sqr-sqrt54.7%
associate-*r/54.5%
add-sqr-sqrt28.3%
associate-/l/28.4%
div-inv28.3%
clear-num28.3%
associate-*l/28.3%
*-un-lft-identity28.3%
pow1/228.3%
pow-flip28.6%
metadata-eval28.6%
div-inv28.6%
*-commutative28.6%
associate-/r*28.6%
metadata-eval28.6%
Applied egg-rr28.6%
Final simplification28.6%
y_m = (fabs.f64 y) (FPCore (x y_m) :precision binary64 (let* ((t_0 (cbrt (/ x y_m)))) (/ 1.0 (cos (* (pow t_0 2.0) (* -0.5 t_0))))))
y_m = fabs(y);
double code(double x, double y_m) {
double t_0 = cbrt((x / y_m));
return 1.0 / cos((pow(t_0, 2.0) * (-0.5 * t_0)));
}
y_m = Math.abs(y);
public static double code(double x, double y_m) {
double t_0 = Math.cbrt((x / y_m));
return 1.0 / Math.cos((Math.pow(t_0, 2.0) * (-0.5 * t_0)));
}
y_m = abs(y) function code(x, y_m) t_0 = cbrt(Float64(x / y_m)) return Float64(1.0 / cos(Float64((t_0 ^ 2.0) * Float64(-0.5 * t_0)))) end
y_m = N[Abs[y], $MachinePrecision]
code[x_, y$95$m_] := Block[{t$95$0 = N[Power[N[(x / y$95$m), $MachinePrecision], 1/3], $MachinePrecision]}, N[(1.0 / N[Cos[N[(N[Power[t$95$0, 2.0], $MachinePrecision] * N[(-0.5 * t$95$0), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
y_m = \left|y\right|
\\
\begin{array}{l}
t_0 := \sqrt[3]{\frac{x}{y\_m}}\\
\frac{1}{\cos \left({t\_0}^{2} \cdot \left(-0.5 \cdot t\_0\right)\right)}
\end{array}
\end{array}
Initial program 40.4%
remove-double-neg40.4%
distribute-frac-neg40.4%
tan-neg40.4%
distribute-frac-neg240.4%
distribute-lft-neg-out40.4%
distribute-frac-neg240.4%
distribute-lft-neg-out40.4%
distribute-frac-neg240.4%
distribute-frac-neg40.4%
neg-mul-140.4%
*-commutative40.4%
associate-/l*40.5%
*-commutative40.5%
associate-/r*40.5%
metadata-eval40.5%
sin-neg40.5%
distribute-frac-neg40.5%
Simplified40.6%
Taylor expanded in x around inf 54.5%
associate-*r/54.5%
Simplified54.5%
Taylor expanded in x around inf 54.5%
associate-*r/54.5%
*-commutative54.5%
associate-/l*54.7%
Simplified54.7%
associate-*r/54.5%
add-sqr-sqrt28.3%
associate-/r*28.4%
Applied egg-rr28.4%
associate-/l/28.3%
add-sqr-sqrt54.5%
associate-*l/54.5%
add-cube-cbrt55.3%
associate-*l*55.3%
pow255.3%
Applied egg-rr55.3%
Final simplification55.3%
y_m = (fabs.f64 y) (FPCore (x y_m) :precision binary64 (/ 1.0 (cos (/ (/ (* -0.5 x) (sqrt y_m)) (sqrt y_m)))))
y_m = fabs(y);
double code(double x, double y_m) {
return 1.0 / cos((((-0.5 * x) / sqrt(y_m)) / sqrt(y_m)));
}
y_m = abs(y)
real(8) function code(x, y_m)
real(8), intent (in) :: x
real(8), intent (in) :: y_m
code = 1.0d0 / cos(((((-0.5d0) * x) / sqrt(y_m)) / sqrt(y_m)))
end function
y_m = Math.abs(y);
public static double code(double x, double y_m) {
return 1.0 / Math.cos((((-0.5 * x) / Math.sqrt(y_m)) / Math.sqrt(y_m)));
}
y_m = math.fabs(y) def code(x, y_m): return 1.0 / math.cos((((-0.5 * x) / math.sqrt(y_m)) / math.sqrt(y_m)))
y_m = abs(y) function code(x, y_m) return Float64(1.0 / cos(Float64(Float64(Float64(-0.5 * x) / sqrt(y_m)) / sqrt(y_m)))) end
y_m = abs(y); function tmp = code(x, y_m) tmp = 1.0 / cos((((-0.5 * x) / sqrt(y_m)) / sqrt(y_m))); end
y_m = N[Abs[y], $MachinePrecision] code[x_, y$95$m_] := N[(1.0 / N[Cos[N[(N[(N[(-0.5 * x), $MachinePrecision] / N[Sqrt[y$95$m], $MachinePrecision]), $MachinePrecision] / N[Sqrt[y$95$m], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
y_m = \left|y\right|
\\
\frac{1}{\cos \left(\frac{\frac{-0.5 \cdot x}{\sqrt{y\_m}}}{\sqrt{y\_m}}\right)}
\end{array}
Initial program 40.4%
remove-double-neg40.4%
distribute-frac-neg40.4%
tan-neg40.4%
distribute-frac-neg240.4%
distribute-lft-neg-out40.4%
distribute-frac-neg240.4%
distribute-lft-neg-out40.4%
distribute-frac-neg240.4%
distribute-frac-neg40.4%
neg-mul-140.4%
*-commutative40.4%
associate-/l*40.5%
*-commutative40.5%
associate-/r*40.5%
metadata-eval40.5%
sin-neg40.5%
distribute-frac-neg40.5%
Simplified40.6%
Taylor expanded in x around inf 54.5%
associate-*r/54.5%
Simplified54.5%
Taylor expanded in x around inf 54.5%
associate-*r/54.5%
*-commutative54.5%
associate-/l*54.7%
Simplified54.7%
associate-*r/54.5%
add-sqr-sqrt28.3%
associate-/r*28.4%
Applied egg-rr28.4%
Final simplification28.4%
y_m = (fabs.f64 y) (FPCore (x y_m) :precision binary64 (/ 1.0 (cos (* (/ 1.0 y_m) (/ 1.0 (/ -2.0 x))))))
y_m = fabs(y);
double code(double x, double y_m) {
return 1.0 / cos(((1.0 / y_m) * (1.0 / (-2.0 / x))));
}
y_m = abs(y)
real(8) function code(x, y_m)
real(8), intent (in) :: x
real(8), intent (in) :: y_m
code = 1.0d0 / cos(((1.0d0 / y_m) * (1.0d0 / ((-2.0d0) / x))))
end function
y_m = Math.abs(y);
public static double code(double x, double y_m) {
return 1.0 / Math.cos(((1.0 / y_m) * (1.0 / (-2.0 / x))));
}
y_m = math.fabs(y) def code(x, y_m): return 1.0 / math.cos(((1.0 / y_m) * (1.0 / (-2.0 / x))))
y_m = abs(y) function code(x, y_m) return Float64(1.0 / cos(Float64(Float64(1.0 / y_m) * Float64(1.0 / Float64(-2.0 / x))))) end
y_m = abs(y); function tmp = code(x, y_m) tmp = 1.0 / cos(((1.0 / y_m) * (1.0 / (-2.0 / x)))); end
y_m = N[Abs[y], $MachinePrecision] code[x_, y$95$m_] := N[(1.0 / N[Cos[N[(N[(1.0 / y$95$m), $MachinePrecision] * N[(1.0 / N[(-2.0 / x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
y_m = \left|y\right|
\\
\frac{1}{\cos \left(\frac{1}{y\_m} \cdot \frac{1}{\frac{-2}{x}}\right)}
\end{array}
Initial program 40.4%
remove-double-neg40.4%
distribute-frac-neg40.4%
tan-neg40.4%
distribute-frac-neg240.4%
distribute-lft-neg-out40.4%
distribute-frac-neg240.4%
distribute-lft-neg-out40.4%
distribute-frac-neg240.4%
distribute-frac-neg40.4%
neg-mul-140.4%
*-commutative40.4%
associate-/l*40.5%
*-commutative40.5%
associate-/r*40.5%
metadata-eval40.5%
sin-neg40.5%
distribute-frac-neg40.5%
Simplified40.6%
Taylor expanded in x around inf 54.5%
associate-*r/54.5%
Simplified54.5%
Taylor expanded in x around inf 54.5%
associate-*r/54.5%
*-commutative54.5%
associate-/l*54.7%
Simplified54.7%
associate-*r/54.5%
add-sqr-sqrt28.3%
associate-/r*28.4%
Applied egg-rr28.4%
associate-/l/28.3%
add-sqr-sqrt54.5%
clear-num54.4%
inv-pow54.4%
div-inv54.6%
unpow-prod-down55.0%
inv-pow55.0%
*-commutative55.0%
associate-/r*55.0%
metadata-eval55.0%
Applied egg-rr55.0%
unpow-155.0%
Simplified55.0%
Final simplification55.0%
y_m = (fabs.f64 y) (FPCore (x y_m) :precision binary64 (/ 1.0 (cos (* x (/ -0.5 y_m)))))
y_m = fabs(y);
double code(double x, double y_m) {
return 1.0 / cos((x * (-0.5 / y_m)));
}
y_m = abs(y)
real(8) function code(x, y_m)
real(8), intent (in) :: x
real(8), intent (in) :: y_m
code = 1.0d0 / cos((x * ((-0.5d0) / y_m)))
end function
y_m = Math.abs(y);
public static double code(double x, double y_m) {
return 1.0 / Math.cos((x * (-0.5 / y_m)));
}
y_m = math.fabs(y) def code(x, y_m): return 1.0 / math.cos((x * (-0.5 / y_m)))
y_m = abs(y) function code(x, y_m) return Float64(1.0 / cos(Float64(x * Float64(-0.5 / y_m)))) end
y_m = abs(y); function tmp = code(x, y_m) tmp = 1.0 / cos((x * (-0.5 / y_m))); end
y_m = N[Abs[y], $MachinePrecision] code[x_, y$95$m_] := N[(1.0 / N[Cos[N[(x * N[(-0.5 / y$95$m), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
y_m = \left|y\right|
\\
\frac{1}{\cos \left(x \cdot \frac{-0.5}{y\_m}\right)}
\end{array}
Initial program 40.4%
remove-double-neg40.4%
distribute-frac-neg40.4%
tan-neg40.4%
distribute-frac-neg240.4%
distribute-lft-neg-out40.4%
distribute-frac-neg240.4%
distribute-lft-neg-out40.4%
distribute-frac-neg240.4%
distribute-frac-neg40.4%
neg-mul-140.4%
*-commutative40.4%
associate-/l*40.5%
*-commutative40.5%
associate-/r*40.5%
metadata-eval40.5%
sin-neg40.5%
distribute-frac-neg40.5%
Simplified40.6%
Taylor expanded in x around inf 54.5%
associate-*r/54.5%
Simplified54.5%
Taylor expanded in x around inf 54.5%
associate-*r/54.5%
*-commutative54.5%
associate-/l*54.7%
Simplified54.7%
Final simplification54.7%
y_m = (fabs.f64 y) (FPCore (x y_m) :precision binary64 1.0)
y_m = fabs(y);
double code(double x, double y_m) {
return 1.0;
}
y_m = abs(y)
real(8) function code(x, y_m)
real(8), intent (in) :: x
real(8), intent (in) :: y_m
code = 1.0d0
end function
y_m = Math.abs(y);
public static double code(double x, double y_m) {
return 1.0;
}
y_m = math.fabs(y) def code(x, y_m): return 1.0
y_m = abs(y) function code(x, y_m) return 1.0 end
y_m = abs(y); function tmp = code(x, y_m) tmp = 1.0; end
y_m = N[Abs[y], $MachinePrecision] code[x_, y$95$m_] := 1.0
\begin{array}{l}
y_m = \left|y\right|
\\
1
\end{array}
Initial program 40.4%
remove-double-neg40.4%
distribute-frac-neg40.4%
tan-neg40.4%
distribute-frac-neg240.4%
distribute-lft-neg-out40.4%
distribute-frac-neg240.4%
distribute-lft-neg-out40.4%
distribute-frac-neg240.4%
distribute-frac-neg40.4%
neg-mul-140.4%
*-commutative40.4%
associate-/l*40.5%
*-commutative40.5%
associate-/r*40.5%
metadata-eval40.5%
sin-neg40.5%
distribute-frac-neg40.5%
Simplified40.6%
Taylor expanded in x around 0 54.6%
Final simplification54.6%
(FPCore (x y)
:precision binary64
(let* ((t_0 (/ x (* y 2.0))) (t_1 (sin t_0)))
(if (< y -1.2303690911306994e+114)
1.0
(if (< y -9.102852406811914e-222)
(/ t_1 (* t_1 (log (exp (cos t_0)))))
1.0))))
double code(double x, double y) {
double t_0 = x / (y * 2.0);
double t_1 = sin(t_0);
double tmp;
if (y < -1.2303690911306994e+114) {
tmp = 1.0;
} else if (y < -9.102852406811914e-222) {
tmp = t_1 / (t_1 * log(exp(cos(t_0))));
} else {
tmp = 1.0;
}
return tmp;
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8) :: t_0
real(8) :: t_1
real(8) :: tmp
t_0 = x / (y * 2.0d0)
t_1 = sin(t_0)
if (y < (-1.2303690911306994d+114)) then
tmp = 1.0d0
else if (y < (-9.102852406811914d-222)) then
tmp = t_1 / (t_1 * log(exp(cos(t_0))))
else
tmp = 1.0d0
end if
code = tmp
end function
public static double code(double x, double y) {
double t_0 = x / (y * 2.0);
double t_1 = Math.sin(t_0);
double tmp;
if (y < -1.2303690911306994e+114) {
tmp = 1.0;
} else if (y < -9.102852406811914e-222) {
tmp = t_1 / (t_1 * Math.log(Math.exp(Math.cos(t_0))));
} else {
tmp = 1.0;
}
return tmp;
}
def code(x, y): t_0 = x / (y * 2.0) t_1 = math.sin(t_0) tmp = 0 if y < -1.2303690911306994e+114: tmp = 1.0 elif y < -9.102852406811914e-222: tmp = t_1 / (t_1 * math.log(math.exp(math.cos(t_0)))) else: tmp = 1.0 return tmp
function code(x, y) t_0 = Float64(x / Float64(y * 2.0)) t_1 = sin(t_0) tmp = 0.0 if (y < -1.2303690911306994e+114) tmp = 1.0; elseif (y < -9.102852406811914e-222) tmp = Float64(t_1 / Float64(t_1 * log(exp(cos(t_0))))); else tmp = 1.0; end return tmp end
function tmp_2 = code(x, y) t_0 = x / (y * 2.0); t_1 = sin(t_0); tmp = 0.0; if (y < -1.2303690911306994e+114) tmp = 1.0; elseif (y < -9.102852406811914e-222) tmp = t_1 / (t_1 * log(exp(cos(t_0)))); else tmp = 1.0; end tmp_2 = tmp; end
code[x_, y_] := Block[{t$95$0 = N[(x / N[(y * 2.0), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[Sin[t$95$0], $MachinePrecision]}, If[Less[y, -1.2303690911306994e+114], 1.0, If[Less[y, -9.102852406811914e-222], N[(t$95$1 / N[(t$95$1 * N[Log[N[Exp[N[Cos[t$95$0], $MachinePrecision]], $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 1.0]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \frac{x}{y \cdot 2}\\
t_1 := \sin t\_0\\
\mathbf{if}\;y < -1.2303690911306994 \cdot 10^{+114}:\\
\;\;\;\;1\\
\mathbf{elif}\;y < -9.102852406811914 \cdot 10^{-222}:\\
\;\;\;\;\frac{t\_1}{t\_1 \cdot \log \left(e^{\cos t\_0}\right)}\\
\mathbf{else}:\\
\;\;\;\;1\\
\end{array}
\end{array}
herbie shell --seed 2024066
(FPCore (x y)
:name "Diagrams.TwoD.Layout.CirclePacking:approxRadius from diagrams-contrib-1.3.0.5"
:precision binary64
:alt
(if (< y -1.2303690911306994e+114) 1.0 (if (< y -9.102852406811914e-222) (/ (sin (/ x (* y 2.0))) (* (sin (/ x (* y 2.0))) (log (exp (cos (/ x (* y 2.0))))))) 1.0))
(/ (tan (/ x (* y 2.0))) (sin (/ x (* y 2.0)))))