
(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 7 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}
x_m = (fabs.f64 x) y_m = (fabs.f64 y) (FPCore (x_m y_m) :precision binary64 (if (<= (/ x_m (* y_m 2.0)) 2e+64) (/ 1.0 (cos (/ (* x_m (/ -0.5 (cbrt y_m))) (pow (cbrt y_m) 2.0)))) (fabs (* 2.0 (pow (cbrt 0.5) 3.0)))))
x_m = fabs(x);
y_m = fabs(y);
double code(double x_m, double y_m) {
double tmp;
if ((x_m / (y_m * 2.0)) <= 2e+64) {
tmp = 1.0 / cos(((x_m * (-0.5 / cbrt(y_m))) / pow(cbrt(y_m), 2.0)));
} else {
tmp = fabs((2.0 * pow(cbrt(0.5), 3.0)));
}
return tmp;
}
x_m = Math.abs(x);
y_m = Math.abs(y);
public static double code(double x_m, double y_m) {
double tmp;
if ((x_m / (y_m * 2.0)) <= 2e+64) {
tmp = 1.0 / Math.cos(((x_m * (-0.5 / Math.cbrt(y_m))) / Math.pow(Math.cbrt(y_m), 2.0)));
} else {
tmp = Math.abs((2.0 * Math.pow(Math.cbrt(0.5), 3.0)));
}
return tmp;
}
x_m = abs(x) y_m = abs(y) function code(x_m, y_m) tmp = 0.0 if (Float64(x_m / Float64(y_m * 2.0)) <= 2e+64) tmp = Float64(1.0 / cos(Float64(Float64(x_m * Float64(-0.5 / cbrt(y_m))) / (cbrt(y_m) ^ 2.0)))); else tmp = abs(Float64(2.0 * (cbrt(0.5) ^ 3.0))); end return tmp end
x_m = N[Abs[x], $MachinePrecision] y_m = N[Abs[y], $MachinePrecision] code[x$95$m_, y$95$m_] := If[LessEqual[N[(x$95$m / N[(y$95$m * 2.0), $MachinePrecision]), $MachinePrecision], 2e+64], N[(1.0 / N[Cos[N[(N[(x$95$m * N[(-0.5 / N[Power[y$95$m, 1/3], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[Power[N[Power[y$95$m, 1/3], $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], N[Abs[N[(2.0 * N[Power[N[Power[0.5, 1/3], $MachinePrecision], 3.0], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]]
\begin{array}{l}
x_m = \left|x\right|
\\
y_m = \left|y\right|
\\
\begin{array}{l}
\mathbf{if}\;\frac{x\_m}{y\_m \cdot 2} \leq 2 \cdot 10^{+64}:\\
\;\;\;\;\frac{1}{\cos \left(\frac{x\_m \cdot \frac{-0.5}{\sqrt[3]{y\_m}}}{{\left(\sqrt[3]{y\_m}\right)}^{2}}\right)}\\
\mathbf{else}:\\
\;\;\;\;\left|2 \cdot {\left(\sqrt[3]{0.5}\right)}^{3}\right|\\
\end{array}
\end{array}
if (/.f64 x (*.f64 y #s(literal 2 binary64))) < 2.00000000000000004e64Initial program 52.2%
expm1-log1p-u52.2%
expm1-undefine6.6%
*-un-lft-identity6.6%
*-commutative6.6%
times-frac6.6%
metadata-eval6.6%
Applied egg-rr6.6%
expm1-define52.1%
*-commutative52.1%
associate-*l/52.2%
Simplified52.2%
Taylor expanded in x around inf 68.3%
associate-*r/68.3%
*-commutative68.3%
associate-*r/68.4%
Simplified68.4%
clear-num68.4%
div-inv68.4%
metadata-eval68.4%
div-inv68.3%
associate-/r*68.3%
clear-num68.5%
Applied egg-rr68.5%
add-sqr-sqrt37.2%
sqrt-unprod67.8%
associate-/r/67.7%
metadata-eval67.7%
associate-/r/67.6%
metadata-eval67.6%
swap-sqr67.6%
metadata-eval67.6%
metadata-eval67.6%
swap-sqr67.6%
associate-*r/67.6%
*-commutative67.6%
add-cube-cbrt67.7%
unpow267.7%
associate-/l/67.7%
rem-cube-cbrt67.4%
associate-*r/67.4%
*-commutative67.4%
add-cube-cbrt67.5%
unpow267.5%
associate-/l/67.5%
Applied egg-rr69.2%
if 2.00000000000000004e64 < (/.f64 x (*.f64 y #s(literal 2 binary64))) Initial program 8.1%
remove-double-neg8.1%
distribute-frac-neg8.1%
tan-neg8.1%
distribute-frac-neg28.1%
distribute-lft-neg-out8.1%
distribute-frac-neg28.1%
distribute-lft-neg-out8.1%
distribute-frac-neg28.1%
distribute-frac-neg8.1%
neg-mul-18.1%
*-commutative8.1%
associate-/l*8.0%
*-commutative8.0%
associate-/r*8.0%
metadata-eval8.0%
sin-neg8.0%
distribute-frac-neg8.0%
Simplified8.2%
add-sqr-sqrt5.0%
sqrt-unprod8.8%
pow28.8%
Applied egg-rr8.8%
unpow28.8%
rem-sqrt-square8.8%
remove-double-neg8.8%
distribute-frac-neg8.8%
distribute-neg-frac28.8%
tan-neg8.8%
associate-*r/8.6%
distribute-frac-neg8.6%
distribute-rgt-neg-in8.6%
metadata-eval8.6%
*-commutative8.6%
associate-*r/8.6%
sin-neg8.6%
associate-*r/8.8%
distribute-frac-neg8.8%
Simplified8.8%
associate-*r/8.6%
metadata-eval8.6%
distribute-rgt-neg-in8.6%
*-commutative8.6%
distribute-neg-frac8.6%
neg-sub08.6%
associate-*r/8.6%
clear-num7.9%
un-div-inv7.9%
Applied egg-rr7.9%
neg-sub07.9%
distribute-neg-frac7.9%
metadata-eval7.9%
Simplified7.9%
add-cube-cbrt7.9%
pow37.9%
associate-/r/7.6%
*-commutative7.6%
associate-*r/7.7%
Applied egg-rr7.7%
Taylor expanded in x around 0 11.9%
Final simplification54.9%
x_m = (fabs.f64 x) y_m = (fabs.f64 y) (FPCore (x_m y_m) :precision binary64 (if (<= (/ x_m (* y_m 2.0)) 100000000000.0) (/ 1.0 (cbrt (pow (cos (* -0.5 (/ x_m y_m))) 3.0))) (fabs (* 2.0 (pow (cbrt 0.5) 3.0)))))
x_m = fabs(x);
y_m = fabs(y);
double code(double x_m, double y_m) {
double tmp;
if ((x_m / (y_m * 2.0)) <= 100000000000.0) {
tmp = 1.0 / cbrt(pow(cos((-0.5 * (x_m / y_m))), 3.0));
} else {
tmp = fabs((2.0 * pow(cbrt(0.5), 3.0)));
}
return tmp;
}
x_m = Math.abs(x);
y_m = Math.abs(y);
public static double code(double x_m, double y_m) {
double tmp;
if ((x_m / (y_m * 2.0)) <= 100000000000.0) {
tmp = 1.0 / Math.cbrt(Math.pow(Math.cos((-0.5 * (x_m / y_m))), 3.0));
} else {
tmp = Math.abs((2.0 * Math.pow(Math.cbrt(0.5), 3.0)));
}
return tmp;
}
x_m = abs(x) y_m = abs(y) function code(x_m, y_m) tmp = 0.0 if (Float64(x_m / Float64(y_m * 2.0)) <= 100000000000.0) tmp = Float64(1.0 / cbrt((cos(Float64(-0.5 * Float64(x_m / y_m))) ^ 3.0))); else tmp = abs(Float64(2.0 * (cbrt(0.5) ^ 3.0))); end return tmp end
x_m = N[Abs[x], $MachinePrecision] y_m = N[Abs[y], $MachinePrecision] code[x$95$m_, y$95$m_] := If[LessEqual[N[(x$95$m / N[(y$95$m * 2.0), $MachinePrecision]), $MachinePrecision], 100000000000.0], N[(1.0 / N[Power[N[Power[N[Cos[N[(-0.5 * N[(x$95$m / y$95$m), $MachinePrecision]), $MachinePrecision]], $MachinePrecision], 3.0], $MachinePrecision], 1/3], $MachinePrecision]), $MachinePrecision], N[Abs[N[(2.0 * N[Power[N[Power[0.5, 1/3], $MachinePrecision], 3.0], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]]
\begin{array}{l}
x_m = \left|x\right|
\\
y_m = \left|y\right|
\\
\begin{array}{l}
\mathbf{if}\;\frac{x\_m}{y\_m \cdot 2} \leq 100000000000:\\
\;\;\;\;\frac{1}{\sqrt[3]{{\cos \left(-0.5 \cdot \frac{x\_m}{y\_m}\right)}^{3}}}\\
\mathbf{else}:\\
\;\;\;\;\left|2 \cdot {\left(\sqrt[3]{0.5}\right)}^{3}\right|\\
\end{array}
\end{array}
if (/.f64 x (*.f64 y #s(literal 2 binary64))) < 1e11Initial program 54.9%
remove-double-neg54.9%
distribute-frac-neg54.9%
tan-neg54.9%
distribute-frac-neg254.9%
distribute-lft-neg-out54.9%
distribute-frac-neg254.9%
distribute-lft-neg-out54.9%
distribute-frac-neg254.9%
distribute-frac-neg54.9%
neg-mul-154.9%
*-commutative54.9%
associate-/l*54.8%
*-commutative54.8%
associate-/r*54.8%
metadata-eval54.8%
sin-neg54.8%
distribute-frac-neg54.8%
Simplified54.9%
Taylor expanded in x around inf 72.1%
associate-*r/72.1%
Simplified72.1%
add-cbrt-cube72.1%
pow372.1%
associate-*r/72.1%
Applied egg-rr72.1%
if 1e11 < (/.f64 x (*.f64 y #s(literal 2 binary64))) Initial program 8.8%
remove-double-neg8.8%
distribute-frac-neg8.8%
tan-neg8.8%
distribute-frac-neg28.8%
distribute-lft-neg-out8.8%
distribute-frac-neg28.8%
distribute-lft-neg-out8.8%
distribute-frac-neg28.8%
distribute-frac-neg8.8%
neg-mul-18.8%
*-commutative8.8%
associate-/l*8.3%
*-commutative8.3%
associate-/r*8.3%
metadata-eval8.3%
sin-neg8.3%
distribute-frac-neg8.3%
Simplified8.9%
add-sqr-sqrt5.1%
sqrt-unprod9.1%
pow29.1%
Applied egg-rr9.1%
unpow29.1%
rem-sqrt-square9.1%
remove-double-neg9.1%
distribute-frac-neg9.1%
distribute-neg-frac29.1%
tan-neg9.1%
associate-*r/8.8%
distribute-frac-neg8.8%
distribute-rgt-neg-in8.8%
metadata-eval8.8%
*-commutative8.8%
associate-*r/8.8%
sin-neg8.8%
associate-*r/9.0%
distribute-frac-neg9.0%
Simplified9.1%
associate-*r/8.8%
metadata-eval8.8%
distribute-rgt-neg-in8.8%
*-commutative8.8%
distribute-neg-frac8.8%
neg-sub08.8%
associate-*r/8.8%
clear-num8.4%
un-div-inv8.4%
Applied egg-rr8.4%
neg-sub08.4%
distribute-neg-frac8.4%
metadata-eval8.4%
Simplified8.4%
add-cube-cbrt8.3%
pow38.3%
associate-/r/8.1%
*-commutative8.1%
associate-*r/8.1%
Applied egg-rr8.1%
Taylor expanded in x around 0 11.9%
x_m = (fabs.f64 x) y_m = (fabs.f64 y) (FPCore (x_m y_m) :precision binary64 (if (<= (/ x_m (* y_m 2.0)) 4000000000.0) (log1p (expm1 (/ 1.0 (cos (* x_m (/ -0.5 y_m)))))) (fabs (* 2.0 (pow (cbrt 0.5) 3.0)))))
x_m = fabs(x);
y_m = fabs(y);
double code(double x_m, double y_m) {
double tmp;
if ((x_m / (y_m * 2.0)) <= 4000000000.0) {
tmp = log1p(expm1((1.0 / cos((x_m * (-0.5 / y_m))))));
} else {
tmp = fabs((2.0 * pow(cbrt(0.5), 3.0)));
}
return tmp;
}
x_m = Math.abs(x);
y_m = Math.abs(y);
public static double code(double x_m, double y_m) {
double tmp;
if ((x_m / (y_m * 2.0)) <= 4000000000.0) {
tmp = Math.log1p(Math.expm1((1.0 / Math.cos((x_m * (-0.5 / y_m))))));
} else {
tmp = Math.abs((2.0 * Math.pow(Math.cbrt(0.5), 3.0)));
}
return tmp;
}
x_m = abs(x) y_m = abs(y) function code(x_m, y_m) tmp = 0.0 if (Float64(x_m / Float64(y_m * 2.0)) <= 4000000000.0) tmp = log1p(expm1(Float64(1.0 / cos(Float64(x_m * Float64(-0.5 / y_m)))))); else tmp = abs(Float64(2.0 * (cbrt(0.5) ^ 3.0))); end return tmp end
x_m = N[Abs[x], $MachinePrecision] y_m = N[Abs[y], $MachinePrecision] code[x$95$m_, y$95$m_] := If[LessEqual[N[(x$95$m / N[(y$95$m * 2.0), $MachinePrecision]), $MachinePrecision], 4000000000.0], N[Log[1 + N[(Exp[N[(1.0 / N[Cos[N[(x$95$m * N[(-0.5 / y$95$m), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]] - 1), $MachinePrecision]], $MachinePrecision], N[Abs[N[(2.0 * N[Power[N[Power[0.5, 1/3], $MachinePrecision], 3.0], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]]
\begin{array}{l}
x_m = \left|x\right|
\\
y_m = \left|y\right|
\\
\begin{array}{l}
\mathbf{if}\;\frac{x\_m}{y\_m \cdot 2} \leq 4000000000:\\
\;\;\;\;\mathsf{log1p}\left(\mathsf{expm1}\left(\frac{1}{\cos \left(x\_m \cdot \frac{-0.5}{y\_m}\right)}\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\left|2 \cdot {\left(\sqrt[3]{0.5}\right)}^{3}\right|\\
\end{array}
\end{array}
if (/.f64 x (*.f64 y #s(literal 2 binary64))) < 4e9Initial program 54.9%
expm1-log1p-u54.9%
expm1-undefine6.2%
*-un-lft-identity6.2%
*-commutative6.2%
times-frac6.2%
metadata-eval6.2%
Applied egg-rr6.2%
expm1-define54.7%
*-commutative54.7%
associate-*l/54.9%
Simplified54.9%
Taylor expanded in x around inf 72.1%
associate-*r/72.1%
*-commutative72.1%
associate-*r/72.1%
Simplified72.1%
clear-num72.1%
div-inv72.1%
metadata-eval72.1%
div-inv72.1%
associate-/r*72.1%
clear-num72.2%
Applied egg-rr72.2%
Applied egg-rr72.1%
if 4e9 < (/.f64 x (*.f64 y #s(literal 2 binary64))) Initial program 8.8%
remove-double-neg8.8%
distribute-frac-neg8.8%
tan-neg8.8%
distribute-frac-neg28.8%
distribute-lft-neg-out8.8%
distribute-frac-neg28.8%
distribute-lft-neg-out8.8%
distribute-frac-neg28.8%
distribute-frac-neg8.8%
neg-mul-18.8%
*-commutative8.8%
associate-/l*8.3%
*-commutative8.3%
associate-/r*8.3%
metadata-eval8.3%
sin-neg8.3%
distribute-frac-neg8.3%
Simplified8.9%
add-sqr-sqrt5.1%
sqrt-unprod9.1%
pow29.1%
Applied egg-rr9.1%
unpow29.1%
rem-sqrt-square9.1%
remove-double-neg9.1%
distribute-frac-neg9.1%
distribute-neg-frac29.1%
tan-neg9.1%
associate-*r/8.8%
distribute-frac-neg8.8%
distribute-rgt-neg-in8.8%
metadata-eval8.8%
*-commutative8.8%
associate-*r/8.8%
sin-neg8.8%
associate-*r/9.0%
distribute-frac-neg9.0%
Simplified9.1%
associate-*r/8.8%
metadata-eval8.8%
distribute-rgt-neg-in8.8%
*-commutative8.8%
distribute-neg-frac8.8%
neg-sub08.8%
associate-*r/8.8%
clear-num8.4%
un-div-inv8.4%
Applied egg-rr8.4%
neg-sub08.4%
distribute-neg-frac8.4%
metadata-eval8.4%
Simplified8.4%
add-cube-cbrt8.3%
pow38.3%
associate-/r/8.1%
*-commutative8.1%
associate-*r/8.1%
Applied egg-rr8.1%
Taylor expanded in x around 0 11.9%
x_m = (fabs.f64 x) y_m = (fabs.f64 y) (FPCore (x_m y_m) :precision binary64 (if (<= (/ x_m (* y_m 2.0)) 4000000000.0) (/ 1.0 (cos (* x_m (/ 0.5 y_m)))) (fabs (* 2.0 (pow (cbrt 0.5) 3.0)))))
x_m = fabs(x);
y_m = fabs(y);
double code(double x_m, double y_m) {
double tmp;
if ((x_m / (y_m * 2.0)) <= 4000000000.0) {
tmp = 1.0 / cos((x_m * (0.5 / y_m)));
} else {
tmp = fabs((2.0 * pow(cbrt(0.5), 3.0)));
}
return tmp;
}
x_m = Math.abs(x);
y_m = Math.abs(y);
public static double code(double x_m, double y_m) {
double tmp;
if ((x_m / (y_m * 2.0)) <= 4000000000.0) {
tmp = 1.0 / Math.cos((x_m * (0.5 / y_m)));
} else {
tmp = Math.abs((2.0 * Math.pow(Math.cbrt(0.5), 3.0)));
}
return tmp;
}
x_m = abs(x) y_m = abs(y) function code(x_m, y_m) tmp = 0.0 if (Float64(x_m / Float64(y_m * 2.0)) <= 4000000000.0) tmp = Float64(1.0 / cos(Float64(x_m * Float64(0.5 / y_m)))); else tmp = abs(Float64(2.0 * (cbrt(0.5) ^ 3.0))); end return tmp end
x_m = N[Abs[x], $MachinePrecision] y_m = N[Abs[y], $MachinePrecision] code[x$95$m_, y$95$m_] := If[LessEqual[N[(x$95$m / N[(y$95$m * 2.0), $MachinePrecision]), $MachinePrecision], 4000000000.0], N[(1.0 / N[Cos[N[(x$95$m * N[(0.5 / y$95$m), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], N[Abs[N[(2.0 * N[Power[N[Power[0.5, 1/3], $MachinePrecision], 3.0], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]]
\begin{array}{l}
x_m = \left|x\right|
\\
y_m = \left|y\right|
\\
\begin{array}{l}
\mathbf{if}\;\frac{x\_m}{y\_m \cdot 2} \leq 4000000000:\\
\;\;\;\;\frac{1}{\cos \left(x\_m \cdot \frac{0.5}{y\_m}\right)}\\
\mathbf{else}:\\
\;\;\;\;\left|2 \cdot {\left(\sqrt[3]{0.5}\right)}^{3}\right|\\
\end{array}
\end{array}
if (/.f64 x (*.f64 y #s(literal 2 binary64))) < 4e9Initial program 54.9%
expm1-log1p-u54.9%
expm1-undefine6.2%
*-un-lft-identity6.2%
*-commutative6.2%
times-frac6.2%
metadata-eval6.2%
Applied egg-rr6.2%
expm1-define54.7%
*-commutative54.7%
associate-*l/54.9%
Simplified54.9%
Taylor expanded in x around inf 72.1%
associate-*r/72.1%
*-commutative72.1%
associate-*r/72.1%
Simplified72.1%
if 4e9 < (/.f64 x (*.f64 y #s(literal 2 binary64))) Initial program 8.8%
remove-double-neg8.8%
distribute-frac-neg8.8%
tan-neg8.8%
distribute-frac-neg28.8%
distribute-lft-neg-out8.8%
distribute-frac-neg28.8%
distribute-lft-neg-out8.8%
distribute-frac-neg28.8%
distribute-frac-neg8.8%
neg-mul-18.8%
*-commutative8.8%
associate-/l*8.3%
*-commutative8.3%
associate-/r*8.3%
metadata-eval8.3%
sin-neg8.3%
distribute-frac-neg8.3%
Simplified8.9%
add-sqr-sqrt5.1%
sqrt-unprod9.1%
pow29.1%
Applied egg-rr9.1%
unpow29.1%
rem-sqrt-square9.1%
remove-double-neg9.1%
distribute-frac-neg9.1%
distribute-neg-frac29.1%
tan-neg9.1%
associate-*r/8.8%
distribute-frac-neg8.8%
distribute-rgt-neg-in8.8%
metadata-eval8.8%
*-commutative8.8%
associate-*r/8.8%
sin-neg8.8%
associate-*r/9.0%
distribute-frac-neg9.0%
Simplified9.1%
associate-*r/8.8%
metadata-eval8.8%
distribute-rgt-neg-in8.8%
*-commutative8.8%
distribute-neg-frac8.8%
neg-sub08.8%
associate-*r/8.8%
clear-num8.4%
un-div-inv8.4%
Applied egg-rr8.4%
neg-sub08.4%
distribute-neg-frac8.4%
metadata-eval8.4%
Simplified8.4%
add-cube-cbrt8.3%
pow38.3%
associate-/r/8.1%
*-commutative8.1%
associate-*r/8.1%
Applied egg-rr8.1%
Taylor expanded in x around 0 11.9%
x_m = (fabs.f64 x) y_m = (fabs.f64 y) (FPCore (x_m y_m) :precision binary64 (if (<= (/ x_m (* y_m 2.0)) 4000000000.0) (/ 1.0 (cos (* x_m (/ 0.5 y_m)))) 1.0))
x_m = fabs(x);
y_m = fabs(y);
double code(double x_m, double y_m) {
double tmp;
if ((x_m / (y_m * 2.0)) <= 4000000000.0) {
tmp = 1.0 / cos((x_m * (0.5 / y_m)));
} else {
tmp = 1.0;
}
return tmp;
}
x_m = abs(x)
y_m = abs(y)
real(8) function code(x_m, y_m)
real(8), intent (in) :: x_m
real(8), intent (in) :: y_m
real(8) :: tmp
if ((x_m / (y_m * 2.0d0)) <= 4000000000.0d0) then
tmp = 1.0d0 / cos((x_m * (0.5d0 / y_m)))
else
tmp = 1.0d0
end if
code = tmp
end function
x_m = Math.abs(x);
y_m = Math.abs(y);
public static double code(double x_m, double y_m) {
double tmp;
if ((x_m / (y_m * 2.0)) <= 4000000000.0) {
tmp = 1.0 / Math.cos((x_m * (0.5 / y_m)));
} else {
tmp = 1.0;
}
return tmp;
}
x_m = math.fabs(x) y_m = math.fabs(y) def code(x_m, y_m): tmp = 0 if (x_m / (y_m * 2.0)) <= 4000000000.0: tmp = 1.0 / math.cos((x_m * (0.5 / y_m))) else: tmp = 1.0 return tmp
x_m = abs(x) y_m = abs(y) function code(x_m, y_m) tmp = 0.0 if (Float64(x_m / Float64(y_m * 2.0)) <= 4000000000.0) tmp = Float64(1.0 / cos(Float64(x_m * Float64(0.5 / y_m)))); else tmp = 1.0; end return tmp end
x_m = abs(x); y_m = abs(y); function tmp_2 = code(x_m, y_m) tmp = 0.0; if ((x_m / (y_m * 2.0)) <= 4000000000.0) tmp = 1.0 / cos((x_m * (0.5 / y_m))); else tmp = 1.0; end tmp_2 = tmp; end
x_m = N[Abs[x], $MachinePrecision] y_m = N[Abs[y], $MachinePrecision] code[x$95$m_, y$95$m_] := If[LessEqual[N[(x$95$m / N[(y$95$m * 2.0), $MachinePrecision]), $MachinePrecision], 4000000000.0], N[(1.0 / N[Cos[N[(x$95$m * N[(0.5 / y$95$m), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], 1.0]
\begin{array}{l}
x_m = \left|x\right|
\\
y_m = \left|y\right|
\\
\begin{array}{l}
\mathbf{if}\;\frac{x\_m}{y\_m \cdot 2} \leq 4000000000:\\
\;\;\;\;\frac{1}{\cos \left(x\_m \cdot \frac{0.5}{y\_m}\right)}\\
\mathbf{else}:\\
\;\;\;\;1\\
\end{array}
\end{array}
if (/.f64 x (*.f64 y #s(literal 2 binary64))) < 4e9Initial program 54.9%
expm1-log1p-u54.9%
expm1-undefine6.2%
*-un-lft-identity6.2%
*-commutative6.2%
times-frac6.2%
metadata-eval6.2%
Applied egg-rr6.2%
expm1-define54.7%
*-commutative54.7%
associate-*l/54.9%
Simplified54.9%
Taylor expanded in x around inf 72.1%
associate-*r/72.1%
*-commutative72.1%
associate-*r/72.1%
Simplified72.1%
if 4e9 < (/.f64 x (*.f64 y #s(literal 2 binary64))) Initial program 8.8%
remove-double-neg8.8%
distribute-frac-neg8.8%
tan-neg8.8%
distribute-frac-neg28.8%
distribute-lft-neg-out8.8%
distribute-frac-neg28.8%
distribute-lft-neg-out8.8%
distribute-frac-neg28.8%
distribute-frac-neg8.8%
neg-mul-18.8%
*-commutative8.8%
associate-/l*8.3%
*-commutative8.3%
associate-/r*8.3%
metadata-eval8.3%
sin-neg8.3%
distribute-frac-neg8.3%
Simplified8.9%
Taylor expanded in x around 0 11.9%
x_m = (fabs.f64 x) y_m = (fabs.f64 y) (FPCore (x_m y_m) :precision binary64 1.0)
x_m = fabs(x);
y_m = fabs(y);
double code(double x_m, double y_m) {
return 1.0;
}
x_m = abs(x)
y_m = abs(y)
real(8) function code(x_m, y_m)
real(8), intent (in) :: x_m
real(8), intent (in) :: y_m
code = 1.0d0
end function
x_m = Math.abs(x);
y_m = Math.abs(y);
public static double code(double x_m, double y_m) {
return 1.0;
}
x_m = math.fabs(x) y_m = math.fabs(y) def code(x_m, y_m): return 1.0
x_m = abs(x) y_m = abs(y) function code(x_m, y_m) return 1.0 end
x_m = abs(x); y_m = abs(y); function tmp = code(x_m, y_m) tmp = 1.0; end
x_m = N[Abs[x], $MachinePrecision] y_m = N[Abs[y], $MachinePrecision] code[x$95$m_, y$95$m_] := 1.0
\begin{array}{l}
x_m = \left|x\right|
\\
y_m = \left|y\right|
\\
1
\end{array}
Initial program 41.2%
remove-double-neg41.2%
distribute-frac-neg41.2%
tan-neg41.2%
distribute-frac-neg241.2%
distribute-lft-neg-out41.2%
distribute-frac-neg241.2%
distribute-lft-neg-out41.2%
distribute-frac-neg241.2%
distribute-frac-neg41.2%
neg-mul-141.2%
*-commutative41.2%
associate-/l*41.0%
*-commutative41.0%
associate-/r*41.0%
metadata-eval41.0%
sin-neg41.0%
distribute-frac-neg41.0%
Simplified41.2%
Taylor expanded in x around 0 53.5%
x_m = (fabs.f64 x) y_m = (fabs.f64 y) (FPCore (x_m y_m) :precision binary64 -1.0)
x_m = fabs(x);
y_m = fabs(y);
double code(double x_m, double y_m) {
return -1.0;
}
x_m = abs(x)
y_m = abs(y)
real(8) function code(x_m, y_m)
real(8), intent (in) :: x_m
real(8), intent (in) :: y_m
code = -1.0d0
end function
x_m = Math.abs(x);
y_m = Math.abs(y);
public static double code(double x_m, double y_m) {
return -1.0;
}
x_m = math.fabs(x) y_m = math.fabs(y) def code(x_m, y_m): return -1.0
x_m = abs(x) y_m = abs(y) function code(x_m, y_m) return -1.0 end
x_m = abs(x); y_m = abs(y); function tmp = code(x_m, y_m) tmp = -1.0; end
x_m = N[Abs[x], $MachinePrecision] y_m = N[Abs[y], $MachinePrecision] code[x$95$m_, y$95$m_] := -1.0
\begin{array}{l}
x_m = \left|x\right|
\\
y_m = \left|y\right|
\\
-1
\end{array}
Initial program 41.2%
expm1-log1p-u41.2%
expm1-undefine7.0%
*-un-lft-identity7.0%
*-commutative7.0%
times-frac7.0%
metadata-eval7.0%
Applied egg-rr7.0%
expm1-define41.1%
*-commutative41.1%
associate-*l/41.2%
Simplified41.2%
*-commutative41.2%
associate-*r/41.1%
add-sqr-sqrt16.6%
metadata-eval16.6%
fabs-sqr16.6%
add-sqr-sqrt19.3%
fabs-mul19.3%
log1p-undefine5.2%
add-sqr-sqrt1.4%
fabs-sqr1.4%
add-sqr-sqrt4.3%
Applied egg-rr4.3%
Taylor expanded in x around 0 7.0%
(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 2024107
(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)))))