
(FPCore (alpha beta) :precision binary64 (let* ((t_0 (+ (+ alpha beta) (* 2.0 1.0)))) (/ (/ (/ (+ (+ (+ alpha beta) (* beta alpha)) 1.0) t_0) t_0) (+ t_0 1.0))))
double code(double alpha, double beta) {
double t_0 = (alpha + beta) + (2.0 * 1.0);
return (((((alpha + beta) + (beta * alpha)) + 1.0) / t_0) / t_0) / (t_0 + 1.0);
}
real(8) function code(alpha, beta)
real(8), intent (in) :: alpha
real(8), intent (in) :: beta
real(8) :: t_0
t_0 = (alpha + beta) + (2.0d0 * 1.0d0)
code = (((((alpha + beta) + (beta * alpha)) + 1.0d0) / t_0) / t_0) / (t_0 + 1.0d0)
end function
public static double code(double alpha, double beta) {
double t_0 = (alpha + beta) + (2.0 * 1.0);
return (((((alpha + beta) + (beta * alpha)) + 1.0) / t_0) / t_0) / (t_0 + 1.0);
}
def code(alpha, beta): t_0 = (alpha + beta) + (2.0 * 1.0) return (((((alpha + beta) + (beta * alpha)) + 1.0) / t_0) / t_0) / (t_0 + 1.0)
function code(alpha, beta) t_0 = Float64(Float64(alpha + beta) + Float64(2.0 * 1.0)) return Float64(Float64(Float64(Float64(Float64(Float64(alpha + beta) + Float64(beta * alpha)) + 1.0) / t_0) / t_0) / Float64(t_0 + 1.0)) end
function tmp = code(alpha, beta) t_0 = (alpha + beta) + (2.0 * 1.0); tmp = (((((alpha + beta) + (beta * alpha)) + 1.0) / t_0) / t_0) / (t_0 + 1.0); end
code[alpha_, beta_] := Block[{t$95$0 = N[(N[(alpha + beta), $MachinePrecision] + N[(2.0 * 1.0), $MachinePrecision]), $MachinePrecision]}, N[(N[(N[(N[(N[(N[(alpha + beta), $MachinePrecision] + N[(beta * alpha), $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision] / t$95$0), $MachinePrecision] / t$95$0), $MachinePrecision] / N[(t$95$0 + 1.0), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \left(\alpha + \beta\right) + 2 \cdot 1\\
\frac{\frac{\frac{\left(\left(\alpha + \beta\right) + \beta \cdot \alpha\right) + 1}{t\_0}}{t\_0}}{t\_0 + 1}
\end{array}
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 16 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (alpha beta) :precision binary64 (let* ((t_0 (+ (+ alpha beta) (* 2.0 1.0)))) (/ (/ (/ (+ (+ (+ alpha beta) (* beta alpha)) 1.0) t_0) t_0) (+ t_0 1.0))))
double code(double alpha, double beta) {
double t_0 = (alpha + beta) + (2.0 * 1.0);
return (((((alpha + beta) + (beta * alpha)) + 1.0) / t_0) / t_0) / (t_0 + 1.0);
}
real(8) function code(alpha, beta)
real(8), intent (in) :: alpha
real(8), intent (in) :: beta
real(8) :: t_0
t_0 = (alpha + beta) + (2.0d0 * 1.0d0)
code = (((((alpha + beta) + (beta * alpha)) + 1.0d0) / t_0) / t_0) / (t_0 + 1.0d0)
end function
public static double code(double alpha, double beta) {
double t_0 = (alpha + beta) + (2.0 * 1.0);
return (((((alpha + beta) + (beta * alpha)) + 1.0) / t_0) / t_0) / (t_0 + 1.0);
}
def code(alpha, beta): t_0 = (alpha + beta) + (2.0 * 1.0) return (((((alpha + beta) + (beta * alpha)) + 1.0) / t_0) / t_0) / (t_0 + 1.0)
function code(alpha, beta) t_0 = Float64(Float64(alpha + beta) + Float64(2.0 * 1.0)) return Float64(Float64(Float64(Float64(Float64(Float64(alpha + beta) + Float64(beta * alpha)) + 1.0) / t_0) / t_0) / Float64(t_0 + 1.0)) end
function tmp = code(alpha, beta) t_0 = (alpha + beta) + (2.0 * 1.0); tmp = (((((alpha + beta) + (beta * alpha)) + 1.0) / t_0) / t_0) / (t_0 + 1.0); end
code[alpha_, beta_] := Block[{t$95$0 = N[(N[(alpha + beta), $MachinePrecision] + N[(2.0 * 1.0), $MachinePrecision]), $MachinePrecision]}, N[(N[(N[(N[(N[(N[(alpha + beta), $MachinePrecision] + N[(beta * alpha), $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision] / t$95$0), $MachinePrecision] / t$95$0), $MachinePrecision] / N[(t$95$0 + 1.0), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \left(\alpha + \beta\right) + 2 \cdot 1\\
\frac{\frac{\frac{\left(\left(\alpha + \beta\right) + \beta \cdot \alpha\right) + 1}{t\_0}}{t\_0}}{t\_0 + 1}
\end{array}
\end{array}
(FPCore (alpha beta)
:precision binary64
(let* ((t_0 (+ (+ beta 2.0) alpha)))
(/
1.0
(*
t_0
(* (/ t_0 (+ 1.0 beta)) (/ (+ alpha (+ beta 3.0)) (+ 1.0 alpha)))))))
double code(double alpha, double beta) {
double t_0 = (beta + 2.0) + alpha;
return 1.0 / (t_0 * ((t_0 / (1.0 + beta)) * ((alpha + (beta + 3.0)) / (1.0 + alpha))));
}
real(8) function code(alpha, beta)
real(8), intent (in) :: alpha
real(8), intent (in) :: beta
real(8) :: t_0
t_0 = (beta + 2.0d0) + alpha
code = 1.0d0 / (t_0 * ((t_0 / (1.0d0 + beta)) * ((alpha + (beta + 3.0d0)) / (1.0d0 + alpha))))
end function
public static double code(double alpha, double beta) {
double t_0 = (beta + 2.0) + alpha;
return 1.0 / (t_0 * ((t_0 / (1.0 + beta)) * ((alpha + (beta + 3.0)) / (1.0 + alpha))));
}
def code(alpha, beta): t_0 = (beta + 2.0) + alpha return 1.0 / (t_0 * ((t_0 / (1.0 + beta)) * ((alpha + (beta + 3.0)) / (1.0 + alpha))))
function code(alpha, beta) t_0 = Float64(Float64(beta + 2.0) + alpha) return Float64(1.0 / Float64(t_0 * Float64(Float64(t_0 / Float64(1.0 + beta)) * Float64(Float64(alpha + Float64(beta + 3.0)) / Float64(1.0 + alpha))))) end
function tmp = code(alpha, beta) t_0 = (beta + 2.0) + alpha; tmp = 1.0 / (t_0 * ((t_0 / (1.0 + beta)) * ((alpha + (beta + 3.0)) / (1.0 + alpha)))); end
code[alpha_, beta_] := Block[{t$95$0 = N[(N[(beta + 2.0), $MachinePrecision] + alpha), $MachinePrecision]}, N[(1.0 / N[(t$95$0 * N[(N[(t$95$0 / N[(1.0 + beta), $MachinePrecision]), $MachinePrecision] * N[(N[(alpha + N[(beta + 3.0), $MachinePrecision]), $MachinePrecision] / N[(1.0 + alpha), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \left(\beta + 2\right) + \alpha\\
\frac{1}{t\_0 \cdot \left(\frac{t\_0}{1 + \beta} \cdot \frac{\alpha + \left(\beta + 3\right)}{1 + \alpha}\right)}
\end{array}
\end{array}
Initial program 95.6%
Simplified83.0%
associate-+r+83.0%
fma-undefine83.0%
*-commutative83.0%
associate-+l+83.0%
+-commutative83.0%
associate-+l+83.0%
*-commutative83.0%
associate-*r*83.0%
associate-+r+83.0%
+-commutative83.0%
associate-/l/93.8%
clear-num93.8%
inv-pow93.8%
Applied egg-rr93.8%
unpow-193.8%
associate-/r/93.8%
Simplified99.8%
Final simplification99.8%
(FPCore (alpha beta) :precision binary64 (let* ((t_0 (+ (+ beta 2.0) alpha))) (* (/ (+ 1.0 alpha) t_0) (/ (/ (+ 1.0 beta) t_0) (+ alpha (+ beta 3.0))))))
double code(double alpha, double beta) {
double t_0 = (beta + 2.0) + alpha;
return ((1.0 + alpha) / t_0) * (((1.0 + beta) / t_0) / (alpha + (beta + 3.0)));
}
real(8) function code(alpha, beta)
real(8), intent (in) :: alpha
real(8), intent (in) :: beta
real(8) :: t_0
t_0 = (beta + 2.0d0) + alpha
code = ((1.0d0 + alpha) / t_0) * (((1.0d0 + beta) / t_0) / (alpha + (beta + 3.0d0)))
end function
public static double code(double alpha, double beta) {
double t_0 = (beta + 2.0) + alpha;
return ((1.0 + alpha) / t_0) * (((1.0 + beta) / t_0) / (alpha + (beta + 3.0)));
}
def code(alpha, beta): t_0 = (beta + 2.0) + alpha return ((1.0 + alpha) / t_0) * (((1.0 + beta) / t_0) / (alpha + (beta + 3.0)))
function code(alpha, beta) t_0 = Float64(Float64(beta + 2.0) + alpha) return Float64(Float64(Float64(1.0 + alpha) / t_0) * Float64(Float64(Float64(1.0 + beta) / t_0) / Float64(alpha + Float64(beta + 3.0)))) end
function tmp = code(alpha, beta) t_0 = (beta + 2.0) + alpha; tmp = ((1.0 + alpha) / t_0) * (((1.0 + beta) / t_0) / (alpha + (beta + 3.0))); end
code[alpha_, beta_] := Block[{t$95$0 = N[(N[(beta + 2.0), $MachinePrecision] + alpha), $MachinePrecision]}, N[(N[(N[(1.0 + alpha), $MachinePrecision] / t$95$0), $MachinePrecision] * N[(N[(N[(1.0 + beta), $MachinePrecision] / t$95$0), $MachinePrecision] / N[(alpha + N[(beta + 3.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \left(\beta + 2\right) + \alpha\\
\frac{1 + \alpha}{t\_0} \cdot \frac{\frac{1 + \beta}{t\_0}}{\alpha + \left(\beta + 3\right)}
\end{array}
\end{array}
Initial program 95.6%
Simplified83.0%
times-frac97.6%
+-commutative97.6%
Applied egg-rr97.6%
+-commutative97.6%
+-commutative97.6%
associate-/r*99.8%
+-commutative99.8%
+-commutative99.8%
+-commutative99.8%
+-commutative99.8%
Simplified99.8%
Final simplification99.8%
(FPCore (alpha beta)
:precision binary64
(if (<= beta 6.8e+38)
(/
1.0
(* (+ (+ beta 2.0) alpha) (/ (* (+ beta 2.0) (+ beta 3.0)) (+ 1.0 beta))))
(/
(* (/ (+ 1.0 alpha) (+ beta (+ 2.0 alpha))) (- 1.0 (* 2.0 (/ alpha beta))))
beta)))
double code(double alpha, double beta) {
double tmp;
if (beta <= 6.8e+38) {
tmp = 1.0 / (((beta + 2.0) + alpha) * (((beta + 2.0) * (beta + 3.0)) / (1.0 + beta)));
} else {
tmp = (((1.0 + alpha) / (beta + (2.0 + alpha))) * (1.0 - (2.0 * (alpha / beta)))) / beta;
}
return tmp;
}
real(8) function code(alpha, beta)
real(8), intent (in) :: alpha
real(8), intent (in) :: beta
real(8) :: tmp
if (beta <= 6.8d+38) then
tmp = 1.0d0 / (((beta + 2.0d0) + alpha) * (((beta + 2.0d0) * (beta + 3.0d0)) / (1.0d0 + beta)))
else
tmp = (((1.0d0 + alpha) / (beta + (2.0d0 + alpha))) * (1.0d0 - (2.0d0 * (alpha / beta)))) / beta
end if
code = tmp
end function
public static double code(double alpha, double beta) {
double tmp;
if (beta <= 6.8e+38) {
tmp = 1.0 / (((beta + 2.0) + alpha) * (((beta + 2.0) * (beta + 3.0)) / (1.0 + beta)));
} else {
tmp = (((1.0 + alpha) / (beta + (2.0 + alpha))) * (1.0 - (2.0 * (alpha / beta)))) / beta;
}
return tmp;
}
def code(alpha, beta): tmp = 0 if beta <= 6.8e+38: tmp = 1.0 / (((beta + 2.0) + alpha) * (((beta + 2.0) * (beta + 3.0)) / (1.0 + beta))) else: tmp = (((1.0 + alpha) / (beta + (2.0 + alpha))) * (1.0 - (2.0 * (alpha / beta)))) / beta return tmp
function code(alpha, beta) tmp = 0.0 if (beta <= 6.8e+38) tmp = Float64(1.0 / Float64(Float64(Float64(beta + 2.0) + alpha) * Float64(Float64(Float64(beta + 2.0) * Float64(beta + 3.0)) / Float64(1.0 + beta)))); else tmp = Float64(Float64(Float64(Float64(1.0 + alpha) / Float64(beta + Float64(2.0 + alpha))) * Float64(1.0 - Float64(2.0 * Float64(alpha / beta)))) / beta); end return tmp end
function tmp_2 = code(alpha, beta) tmp = 0.0; if (beta <= 6.8e+38) tmp = 1.0 / (((beta + 2.0) + alpha) * (((beta + 2.0) * (beta + 3.0)) / (1.0 + beta))); else tmp = (((1.0 + alpha) / (beta + (2.0 + alpha))) * (1.0 - (2.0 * (alpha / beta)))) / beta; end tmp_2 = tmp; end
code[alpha_, beta_] := If[LessEqual[beta, 6.8e+38], N[(1.0 / N[(N[(N[(beta + 2.0), $MachinePrecision] + alpha), $MachinePrecision] * N[(N[(N[(beta + 2.0), $MachinePrecision] * N[(beta + 3.0), $MachinePrecision]), $MachinePrecision] / N[(1.0 + beta), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(N[(N[(1.0 + alpha), $MachinePrecision] / N[(beta + N[(2.0 + alpha), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(1.0 - N[(2.0 * N[(alpha / beta), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / beta), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\beta \leq 6.8 \cdot 10^{+38}:\\
\;\;\;\;\frac{1}{\left(\left(\beta + 2\right) + \alpha\right) \cdot \frac{\left(\beta + 2\right) \cdot \left(\beta + 3\right)}{1 + \beta}}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{1 + \alpha}{\beta + \left(2 + \alpha\right)} \cdot \left(1 - 2 \cdot \frac{\alpha}{\beta}\right)}{\beta}\\
\end{array}
\end{array}
if beta < 6.79999999999999992e38Initial program 99.8%
Simplified92.5%
associate-+r+92.5%
fma-undefine92.5%
*-commutative92.5%
associate-+l+92.5%
+-commutative92.5%
associate-+l+92.5%
*-commutative92.5%
associate-*r*92.5%
associate-+r+92.5%
+-commutative92.5%
associate-/l/99.9%
clear-num99.8%
inv-pow99.8%
Applied egg-rr99.8%
unpow-199.8%
associate-/r/99.9%
Simplified99.9%
Taylor expanded in alpha around 0 69.0%
if 6.79999999999999992e38 < beta Initial program 88.1%
Simplified66.2%
times-frac93.8%
+-commutative93.8%
Applied egg-rr93.8%
Taylor expanded in beta around inf 84.0%
mul-1-neg84.0%
*-commutative84.0%
Simplified84.0%
Taylor expanded in alpha around inf 84.0%
associate-*r/84.0%
Simplified84.0%
associate-*r/84.1%
+-commutative84.1%
+-commutative84.1%
associate-+l+84.1%
+-commutative84.1%
unsub-neg84.1%
associate-/l*84.1%
Applied egg-rr84.1%
Final simplification74.5%
(FPCore (alpha beta)
:precision binary64
(if (<= beta 5e+15)
(/
1.0
(* (+ (+ beta 2.0) alpha) (/ (* (+ beta 2.0) (+ beta 3.0)) (+ 1.0 beta))))
(/ 1.0 (* beta (/ (+ alpha (+ beta 3.0)) (+ 1.0 alpha))))))
double code(double alpha, double beta) {
double tmp;
if (beta <= 5e+15) {
tmp = 1.0 / (((beta + 2.0) + alpha) * (((beta + 2.0) * (beta + 3.0)) / (1.0 + beta)));
} else {
tmp = 1.0 / (beta * ((alpha + (beta + 3.0)) / (1.0 + alpha)));
}
return tmp;
}
real(8) function code(alpha, beta)
real(8), intent (in) :: alpha
real(8), intent (in) :: beta
real(8) :: tmp
if (beta <= 5d+15) then
tmp = 1.0d0 / (((beta + 2.0d0) + alpha) * (((beta + 2.0d0) * (beta + 3.0d0)) / (1.0d0 + beta)))
else
tmp = 1.0d0 / (beta * ((alpha + (beta + 3.0d0)) / (1.0d0 + alpha)))
end if
code = tmp
end function
public static double code(double alpha, double beta) {
double tmp;
if (beta <= 5e+15) {
tmp = 1.0 / (((beta + 2.0) + alpha) * (((beta + 2.0) * (beta + 3.0)) / (1.0 + beta)));
} else {
tmp = 1.0 / (beta * ((alpha + (beta + 3.0)) / (1.0 + alpha)));
}
return tmp;
}
def code(alpha, beta): tmp = 0 if beta <= 5e+15: tmp = 1.0 / (((beta + 2.0) + alpha) * (((beta + 2.0) * (beta + 3.0)) / (1.0 + beta))) else: tmp = 1.0 / (beta * ((alpha + (beta + 3.0)) / (1.0 + alpha))) return tmp
function code(alpha, beta) tmp = 0.0 if (beta <= 5e+15) tmp = Float64(1.0 / Float64(Float64(Float64(beta + 2.0) + alpha) * Float64(Float64(Float64(beta + 2.0) * Float64(beta + 3.0)) / Float64(1.0 + beta)))); else tmp = Float64(1.0 / Float64(beta * Float64(Float64(alpha + Float64(beta + 3.0)) / Float64(1.0 + alpha)))); end return tmp end
function tmp_2 = code(alpha, beta) tmp = 0.0; if (beta <= 5e+15) tmp = 1.0 / (((beta + 2.0) + alpha) * (((beta + 2.0) * (beta + 3.0)) / (1.0 + beta))); else tmp = 1.0 / (beta * ((alpha + (beta + 3.0)) / (1.0 + alpha))); end tmp_2 = tmp; end
code[alpha_, beta_] := If[LessEqual[beta, 5e+15], N[(1.0 / N[(N[(N[(beta + 2.0), $MachinePrecision] + alpha), $MachinePrecision] * N[(N[(N[(beta + 2.0), $MachinePrecision] * N[(beta + 3.0), $MachinePrecision]), $MachinePrecision] / N[(1.0 + beta), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(1.0 / N[(beta * N[(N[(alpha + N[(beta + 3.0), $MachinePrecision]), $MachinePrecision] / N[(1.0 + alpha), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\beta \leq 5 \cdot 10^{+15}:\\
\;\;\;\;\frac{1}{\left(\left(\beta + 2\right) + \alpha\right) \cdot \frac{\left(\beta + 2\right) \cdot \left(\beta + 3\right)}{1 + \beta}}\\
\mathbf{else}:\\
\;\;\;\;\frac{1}{\beta \cdot \frac{\alpha + \left(\beta + 3\right)}{1 + \alpha}}\\
\end{array}
\end{array}
if beta < 5e15Initial program 99.8%
Simplified92.3%
associate-+r+92.3%
fma-undefine92.3%
*-commutative92.3%
associate-+l+92.3%
+-commutative92.3%
associate-+l+92.3%
*-commutative92.3%
associate-*r*92.3%
associate-+r+92.3%
+-commutative92.3%
associate-/l/99.9%
clear-num99.8%
inv-pow99.8%
Applied egg-rr99.8%
unpow-199.8%
associate-/r/99.9%
Simplified99.9%
Taylor expanded in alpha around 0 68.0%
if 5e15 < beta Initial program 88.7%
Taylor expanded in beta around inf 85.3%
clear-num85.3%
inv-pow85.3%
metadata-eval85.3%
associate-+l+85.3%
metadata-eval85.3%
associate-+r+85.3%
+-commutative85.3%
associate-+l+85.3%
Applied egg-rr85.3%
unpow-185.3%
associate-/r/85.3%
associate-+r+85.3%
+-commutative85.3%
Simplified85.3%
Final simplification74.7%
(FPCore (alpha beta) :precision binary64 (if (<= beta 7.1e+16) (/ (/ (+ 1.0 beta) (+ beta 2.0)) (* (+ beta 2.0) (+ beta 3.0))) (/ 1.0 (* beta (/ (+ alpha (+ beta 3.0)) (+ 1.0 alpha))))))
double code(double alpha, double beta) {
double tmp;
if (beta <= 7.1e+16) {
tmp = ((1.0 + beta) / (beta + 2.0)) / ((beta + 2.0) * (beta + 3.0));
} else {
tmp = 1.0 / (beta * ((alpha + (beta + 3.0)) / (1.0 + alpha)));
}
return tmp;
}
real(8) function code(alpha, beta)
real(8), intent (in) :: alpha
real(8), intent (in) :: beta
real(8) :: tmp
if (beta <= 7.1d+16) then
tmp = ((1.0d0 + beta) / (beta + 2.0d0)) / ((beta + 2.0d0) * (beta + 3.0d0))
else
tmp = 1.0d0 / (beta * ((alpha + (beta + 3.0d0)) / (1.0d0 + alpha)))
end if
code = tmp
end function
public static double code(double alpha, double beta) {
double tmp;
if (beta <= 7.1e+16) {
tmp = ((1.0 + beta) / (beta + 2.0)) / ((beta + 2.0) * (beta + 3.0));
} else {
tmp = 1.0 / (beta * ((alpha + (beta + 3.0)) / (1.0 + alpha)));
}
return tmp;
}
def code(alpha, beta): tmp = 0 if beta <= 7.1e+16: tmp = ((1.0 + beta) / (beta + 2.0)) / ((beta + 2.0) * (beta + 3.0)) else: tmp = 1.0 / (beta * ((alpha + (beta + 3.0)) / (1.0 + alpha))) return tmp
function code(alpha, beta) tmp = 0.0 if (beta <= 7.1e+16) tmp = Float64(Float64(Float64(1.0 + beta) / Float64(beta + 2.0)) / Float64(Float64(beta + 2.0) * Float64(beta + 3.0))); else tmp = Float64(1.0 / Float64(beta * Float64(Float64(alpha + Float64(beta + 3.0)) / Float64(1.0 + alpha)))); end return tmp end
function tmp_2 = code(alpha, beta) tmp = 0.0; if (beta <= 7.1e+16) tmp = ((1.0 + beta) / (beta + 2.0)) / ((beta + 2.0) * (beta + 3.0)); else tmp = 1.0 / (beta * ((alpha + (beta + 3.0)) / (1.0 + alpha))); end tmp_2 = tmp; end
code[alpha_, beta_] := If[LessEqual[beta, 7.1e+16], N[(N[(N[(1.0 + beta), $MachinePrecision] / N[(beta + 2.0), $MachinePrecision]), $MachinePrecision] / N[(N[(beta + 2.0), $MachinePrecision] * N[(beta + 3.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(1.0 / N[(beta * N[(N[(alpha + N[(beta + 3.0), $MachinePrecision]), $MachinePrecision] / N[(1.0 + alpha), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\beta \leq 7.1 \cdot 10^{+16}:\\
\;\;\;\;\frac{\frac{1 + \beta}{\beta + 2}}{\left(\beta + 2\right) \cdot \left(\beta + 3\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{1}{\beta \cdot \frac{\alpha + \left(\beta + 3\right)}{1 + \alpha}}\\
\end{array}
\end{array}
if beta < 7.1e16Initial program 99.8%
associate-/l/99.9%
+-commutative99.9%
associate-+l+99.9%
*-commutative99.9%
metadata-eval99.9%
associate-+l+99.9%
metadata-eval99.9%
associate-+l+99.9%
metadata-eval99.9%
metadata-eval99.9%
associate-+l+99.9%
Simplified99.9%
Taylor expanded in alpha around 0 82.1%
+-commutative82.1%
Simplified82.1%
Taylor expanded in alpha around 0 67.2%
if 7.1e16 < beta Initial program 88.7%
Taylor expanded in beta around inf 85.3%
clear-num85.3%
inv-pow85.3%
metadata-eval85.3%
associate-+l+85.3%
metadata-eval85.3%
associate-+r+85.3%
+-commutative85.3%
associate-+l+85.3%
Applied egg-rr85.3%
unpow-185.3%
associate-/r/85.3%
associate-+r+85.3%
+-commutative85.3%
Simplified85.3%
Final simplification74.1%
(FPCore (alpha beta) :precision binary64 (if (<= beta 2.5) (/ 1.0 (* (+ beta 2.0) (+ 6.0 (* beta (+ (* beta 2.0) -1.0))))) (/ 1.0 (* beta (/ (+ alpha (+ beta 3.0)) (+ 1.0 alpha))))))
double code(double alpha, double beta) {
double tmp;
if (beta <= 2.5) {
tmp = 1.0 / ((beta + 2.0) * (6.0 + (beta * ((beta * 2.0) + -1.0))));
} else {
tmp = 1.0 / (beta * ((alpha + (beta + 3.0)) / (1.0 + alpha)));
}
return tmp;
}
real(8) function code(alpha, beta)
real(8), intent (in) :: alpha
real(8), intent (in) :: beta
real(8) :: tmp
if (beta <= 2.5d0) then
tmp = 1.0d0 / ((beta + 2.0d0) * (6.0d0 + (beta * ((beta * 2.0d0) + (-1.0d0)))))
else
tmp = 1.0d0 / (beta * ((alpha + (beta + 3.0d0)) / (1.0d0 + alpha)))
end if
code = tmp
end function
public static double code(double alpha, double beta) {
double tmp;
if (beta <= 2.5) {
tmp = 1.0 / ((beta + 2.0) * (6.0 + (beta * ((beta * 2.0) + -1.0))));
} else {
tmp = 1.0 / (beta * ((alpha + (beta + 3.0)) / (1.0 + alpha)));
}
return tmp;
}
def code(alpha, beta): tmp = 0 if beta <= 2.5: tmp = 1.0 / ((beta + 2.0) * (6.0 + (beta * ((beta * 2.0) + -1.0)))) else: tmp = 1.0 / (beta * ((alpha + (beta + 3.0)) / (1.0 + alpha))) return tmp
function code(alpha, beta) tmp = 0.0 if (beta <= 2.5) tmp = Float64(1.0 / Float64(Float64(beta + 2.0) * Float64(6.0 + Float64(beta * Float64(Float64(beta * 2.0) + -1.0))))); else tmp = Float64(1.0 / Float64(beta * Float64(Float64(alpha + Float64(beta + 3.0)) / Float64(1.0 + alpha)))); end return tmp end
function tmp_2 = code(alpha, beta) tmp = 0.0; if (beta <= 2.5) tmp = 1.0 / ((beta + 2.0) * (6.0 + (beta * ((beta * 2.0) + -1.0)))); else tmp = 1.0 / (beta * ((alpha + (beta + 3.0)) / (1.0 + alpha))); end tmp_2 = tmp; end
code[alpha_, beta_] := If[LessEqual[beta, 2.5], N[(1.0 / N[(N[(beta + 2.0), $MachinePrecision] * N[(6.0 + N[(beta * N[(N[(beta * 2.0), $MachinePrecision] + -1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(1.0 / N[(beta * N[(N[(alpha + N[(beta + 3.0), $MachinePrecision]), $MachinePrecision] / N[(1.0 + alpha), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\beta \leq 2.5:\\
\;\;\;\;\frac{1}{\left(\beta + 2\right) \cdot \left(6 + \beta \cdot \left(\beta \cdot 2 + -1\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{1}{\beta \cdot \frac{\alpha + \left(\beta + 3\right)}{1 + \alpha}}\\
\end{array}
\end{array}
if beta < 2.5Initial program 99.8%
Simplified93.4%
associate-+r+93.4%
fma-undefine93.4%
*-commutative93.4%
associate-+l+93.4%
+-commutative93.4%
associate-+l+93.4%
*-commutative93.4%
associate-*r*93.4%
associate-+r+93.4%
+-commutative93.4%
associate-/l/99.9%
clear-num99.9%
inv-pow99.9%
Applied egg-rr99.9%
unpow-199.9%
associate-/r/99.9%
Simplified99.9%
Taylor expanded in alpha around 0 68.8%
Taylor expanded in beta around 0 68.9%
Taylor expanded in alpha around 0 68.0%
if 2.5 < beta Initial program 88.9%
Taylor expanded in beta around inf 83.6%
clear-num83.6%
inv-pow83.6%
metadata-eval83.6%
associate-+l+83.6%
metadata-eval83.6%
associate-+r+83.6%
+-commutative83.6%
associate-+l+83.6%
Applied egg-rr83.6%
unpow-183.6%
associate-/r/83.7%
associate-+r+83.7%
+-commutative83.7%
Simplified83.7%
Final simplification74.1%
(FPCore (alpha beta) :precision binary64 (if (<= beta 2.8) (/ 1.0 (* (+ (+ beta 2.0) alpha) (- 6.0 beta))) (/ 1.0 (* beta (/ (+ alpha (+ beta 3.0)) (+ 1.0 alpha))))))
double code(double alpha, double beta) {
double tmp;
if (beta <= 2.8) {
tmp = 1.0 / (((beta + 2.0) + alpha) * (6.0 - beta));
} else {
tmp = 1.0 / (beta * ((alpha + (beta + 3.0)) / (1.0 + alpha)));
}
return tmp;
}
real(8) function code(alpha, beta)
real(8), intent (in) :: alpha
real(8), intent (in) :: beta
real(8) :: tmp
if (beta <= 2.8d0) then
tmp = 1.0d0 / (((beta + 2.0d0) + alpha) * (6.0d0 - beta))
else
tmp = 1.0d0 / (beta * ((alpha + (beta + 3.0d0)) / (1.0d0 + alpha)))
end if
code = tmp
end function
public static double code(double alpha, double beta) {
double tmp;
if (beta <= 2.8) {
tmp = 1.0 / (((beta + 2.0) + alpha) * (6.0 - beta));
} else {
tmp = 1.0 / (beta * ((alpha + (beta + 3.0)) / (1.0 + alpha)));
}
return tmp;
}
def code(alpha, beta): tmp = 0 if beta <= 2.8: tmp = 1.0 / (((beta + 2.0) + alpha) * (6.0 - beta)) else: tmp = 1.0 / (beta * ((alpha + (beta + 3.0)) / (1.0 + alpha))) return tmp
function code(alpha, beta) tmp = 0.0 if (beta <= 2.8) tmp = Float64(1.0 / Float64(Float64(Float64(beta + 2.0) + alpha) * Float64(6.0 - beta))); else tmp = Float64(1.0 / Float64(beta * Float64(Float64(alpha + Float64(beta + 3.0)) / Float64(1.0 + alpha)))); end return tmp end
function tmp_2 = code(alpha, beta) tmp = 0.0; if (beta <= 2.8) tmp = 1.0 / (((beta + 2.0) + alpha) * (6.0 - beta)); else tmp = 1.0 / (beta * ((alpha + (beta + 3.0)) / (1.0 + alpha))); end tmp_2 = tmp; end
code[alpha_, beta_] := If[LessEqual[beta, 2.8], N[(1.0 / N[(N[(N[(beta + 2.0), $MachinePrecision] + alpha), $MachinePrecision] * N[(6.0 - beta), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(1.0 / N[(beta * N[(N[(alpha + N[(beta + 3.0), $MachinePrecision]), $MachinePrecision] / N[(1.0 + alpha), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\beta \leq 2.8:\\
\;\;\;\;\frac{1}{\left(\left(\beta + 2\right) + \alpha\right) \cdot \left(6 - \beta\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{1}{\beta \cdot \frac{\alpha + \left(\beta + 3\right)}{1 + \alpha}}\\
\end{array}
\end{array}
if beta < 2.7999999999999998Initial program 99.8%
Simplified93.4%
associate-+r+93.4%
fma-undefine93.4%
*-commutative93.4%
associate-+l+93.4%
+-commutative93.4%
associate-+l+93.4%
*-commutative93.4%
associate-*r*93.4%
associate-+r+93.4%
+-commutative93.4%
associate-/l/99.9%
clear-num99.9%
inv-pow99.9%
Applied egg-rr99.9%
unpow-199.9%
associate-/r/99.9%
Simplified99.9%
Taylor expanded in alpha around 0 68.8%
Taylor expanded in beta around 0 68.6%
mul-1-neg68.6%
unsub-neg68.6%
Simplified68.6%
if 2.7999999999999998 < beta Initial program 88.9%
Taylor expanded in beta around inf 83.6%
clear-num83.6%
inv-pow83.6%
metadata-eval83.6%
associate-+l+83.6%
metadata-eval83.6%
associate-+r+83.6%
+-commutative83.6%
associate-+l+83.6%
Applied egg-rr83.6%
unpow-183.6%
associate-/r/83.7%
associate-+r+83.7%
+-commutative83.7%
Simplified83.7%
Final simplification74.5%
(FPCore (alpha beta) :precision binary64 (if (<= beta 2.9) (/ 1.0 (* (+ (+ beta 2.0) alpha) (- 6.0 beta))) (/ (/ (+ 1.0 alpha) beta) (+ alpha (+ beta 3.0)))))
double code(double alpha, double beta) {
double tmp;
if (beta <= 2.9) {
tmp = 1.0 / (((beta + 2.0) + alpha) * (6.0 - beta));
} else {
tmp = ((1.0 + alpha) / beta) / (alpha + (beta + 3.0));
}
return tmp;
}
real(8) function code(alpha, beta)
real(8), intent (in) :: alpha
real(8), intent (in) :: beta
real(8) :: tmp
if (beta <= 2.9d0) then
tmp = 1.0d0 / (((beta + 2.0d0) + alpha) * (6.0d0 - beta))
else
tmp = ((1.0d0 + alpha) / beta) / (alpha + (beta + 3.0d0))
end if
code = tmp
end function
public static double code(double alpha, double beta) {
double tmp;
if (beta <= 2.9) {
tmp = 1.0 / (((beta + 2.0) + alpha) * (6.0 - beta));
} else {
tmp = ((1.0 + alpha) / beta) / (alpha + (beta + 3.0));
}
return tmp;
}
def code(alpha, beta): tmp = 0 if beta <= 2.9: tmp = 1.0 / (((beta + 2.0) + alpha) * (6.0 - beta)) else: tmp = ((1.0 + alpha) / beta) / (alpha + (beta + 3.0)) return tmp
function code(alpha, beta) tmp = 0.0 if (beta <= 2.9) tmp = Float64(1.0 / Float64(Float64(Float64(beta + 2.0) + alpha) * Float64(6.0 - beta))); else tmp = Float64(Float64(Float64(1.0 + alpha) / beta) / Float64(alpha + Float64(beta + 3.0))); end return tmp end
function tmp_2 = code(alpha, beta) tmp = 0.0; if (beta <= 2.9) tmp = 1.0 / (((beta + 2.0) + alpha) * (6.0 - beta)); else tmp = ((1.0 + alpha) / beta) / (alpha + (beta + 3.0)); end tmp_2 = tmp; end
code[alpha_, beta_] := If[LessEqual[beta, 2.9], N[(1.0 / N[(N[(N[(beta + 2.0), $MachinePrecision] + alpha), $MachinePrecision] * N[(6.0 - beta), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(N[(1.0 + alpha), $MachinePrecision] / beta), $MachinePrecision] / N[(alpha + N[(beta + 3.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\beta \leq 2.9:\\
\;\;\;\;\frac{1}{\left(\left(\beta + 2\right) + \alpha\right) \cdot \left(6 - \beta\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{1 + \alpha}{\beta}}{\alpha + \left(\beta + 3\right)}\\
\end{array}
\end{array}
if beta < 2.89999999999999991Initial program 99.8%
Simplified93.4%
associate-+r+93.4%
fma-undefine93.4%
*-commutative93.4%
associate-+l+93.4%
+-commutative93.4%
associate-+l+93.4%
*-commutative93.4%
associate-*r*93.4%
associate-+r+93.4%
+-commutative93.4%
associate-/l/99.9%
clear-num99.9%
inv-pow99.9%
Applied egg-rr99.9%
unpow-199.9%
associate-/r/99.9%
Simplified99.9%
Taylor expanded in alpha around 0 68.8%
Taylor expanded in beta around 0 68.6%
mul-1-neg68.6%
unsub-neg68.6%
Simplified68.6%
if 2.89999999999999991 < beta Initial program 88.9%
Taylor expanded in beta around inf 83.6%
*-un-lft-identity83.6%
metadata-eval83.6%
associate-+l+83.6%
metadata-eval83.6%
associate-+r+83.6%
+-commutative83.6%
associate-+l+83.6%
Applied egg-rr83.6%
*-lft-identity83.6%
associate-+r+83.6%
+-commutative83.6%
Simplified83.6%
Final simplification74.5%
(FPCore (alpha beta) :precision binary64 (if (<= beta 2.9) (/ 1.0 (* (+ (+ beta 2.0) alpha) (- 6.0 beta))) (* (/ (+ 1.0 alpha) beta) (/ 1.0 (+ beta 3.0)))))
double code(double alpha, double beta) {
double tmp;
if (beta <= 2.9) {
tmp = 1.0 / (((beta + 2.0) + alpha) * (6.0 - beta));
} else {
tmp = ((1.0 + alpha) / beta) * (1.0 / (beta + 3.0));
}
return tmp;
}
real(8) function code(alpha, beta)
real(8), intent (in) :: alpha
real(8), intent (in) :: beta
real(8) :: tmp
if (beta <= 2.9d0) then
tmp = 1.0d0 / (((beta + 2.0d0) + alpha) * (6.0d0 - beta))
else
tmp = ((1.0d0 + alpha) / beta) * (1.0d0 / (beta + 3.0d0))
end if
code = tmp
end function
public static double code(double alpha, double beta) {
double tmp;
if (beta <= 2.9) {
tmp = 1.0 / (((beta + 2.0) + alpha) * (6.0 - beta));
} else {
tmp = ((1.0 + alpha) / beta) * (1.0 / (beta + 3.0));
}
return tmp;
}
def code(alpha, beta): tmp = 0 if beta <= 2.9: tmp = 1.0 / (((beta + 2.0) + alpha) * (6.0 - beta)) else: tmp = ((1.0 + alpha) / beta) * (1.0 / (beta + 3.0)) return tmp
function code(alpha, beta) tmp = 0.0 if (beta <= 2.9) tmp = Float64(1.0 / Float64(Float64(Float64(beta + 2.0) + alpha) * Float64(6.0 - beta))); else tmp = Float64(Float64(Float64(1.0 + alpha) / beta) * Float64(1.0 / Float64(beta + 3.0))); end return tmp end
function tmp_2 = code(alpha, beta) tmp = 0.0; if (beta <= 2.9) tmp = 1.0 / (((beta + 2.0) + alpha) * (6.0 - beta)); else tmp = ((1.0 + alpha) / beta) * (1.0 / (beta + 3.0)); end tmp_2 = tmp; end
code[alpha_, beta_] := If[LessEqual[beta, 2.9], N[(1.0 / N[(N[(N[(beta + 2.0), $MachinePrecision] + alpha), $MachinePrecision] * N[(6.0 - beta), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(N[(1.0 + alpha), $MachinePrecision] / beta), $MachinePrecision] * N[(1.0 / N[(beta + 3.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\beta \leq 2.9:\\
\;\;\;\;\frac{1}{\left(\left(\beta + 2\right) + \alpha\right) \cdot \left(6 - \beta\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{1 + \alpha}{\beta} \cdot \frac{1}{\beta + 3}\\
\end{array}
\end{array}
if beta < 2.89999999999999991Initial program 99.8%
Simplified93.4%
associate-+r+93.4%
fma-undefine93.4%
*-commutative93.4%
associate-+l+93.4%
+-commutative93.4%
associate-+l+93.4%
*-commutative93.4%
associate-*r*93.4%
associate-+r+93.4%
+-commutative93.4%
associate-/l/99.9%
clear-num99.9%
inv-pow99.9%
Applied egg-rr99.9%
unpow-199.9%
associate-/r/99.9%
Simplified99.9%
Taylor expanded in alpha around 0 68.8%
Taylor expanded in beta around 0 68.6%
mul-1-neg68.6%
unsub-neg68.6%
Simplified68.6%
if 2.89999999999999991 < beta Initial program 88.9%
Taylor expanded in beta around inf 83.6%
div-inv83.5%
metadata-eval83.5%
associate-+l+83.5%
metadata-eval83.5%
associate-+r+83.5%
+-commutative83.5%
associate-+l+83.5%
Applied egg-rr83.5%
Taylor expanded in alpha around 0 83.3%
+-commutative83.3%
Simplified83.3%
Final simplification74.4%
(FPCore (alpha beta) :precision binary64 (if (<= beta 5.3) (/ 1.0 (* (+ (+ beta 2.0) alpha) 6.0)) (* (/ (+ 1.0 alpha) beta) (/ 1.0 (+ beta 3.0)))))
double code(double alpha, double beta) {
double tmp;
if (beta <= 5.3) {
tmp = 1.0 / (((beta + 2.0) + alpha) * 6.0);
} else {
tmp = ((1.0 + alpha) / beta) * (1.0 / (beta + 3.0));
}
return tmp;
}
real(8) function code(alpha, beta)
real(8), intent (in) :: alpha
real(8), intent (in) :: beta
real(8) :: tmp
if (beta <= 5.3d0) then
tmp = 1.0d0 / (((beta + 2.0d0) + alpha) * 6.0d0)
else
tmp = ((1.0d0 + alpha) / beta) * (1.0d0 / (beta + 3.0d0))
end if
code = tmp
end function
public static double code(double alpha, double beta) {
double tmp;
if (beta <= 5.3) {
tmp = 1.0 / (((beta + 2.0) + alpha) * 6.0);
} else {
tmp = ((1.0 + alpha) / beta) * (1.0 / (beta + 3.0));
}
return tmp;
}
def code(alpha, beta): tmp = 0 if beta <= 5.3: tmp = 1.0 / (((beta + 2.0) + alpha) * 6.0) else: tmp = ((1.0 + alpha) / beta) * (1.0 / (beta + 3.0)) return tmp
function code(alpha, beta) tmp = 0.0 if (beta <= 5.3) tmp = Float64(1.0 / Float64(Float64(Float64(beta + 2.0) + alpha) * 6.0)); else tmp = Float64(Float64(Float64(1.0 + alpha) / beta) * Float64(1.0 / Float64(beta + 3.0))); end return tmp end
function tmp_2 = code(alpha, beta) tmp = 0.0; if (beta <= 5.3) tmp = 1.0 / (((beta + 2.0) + alpha) * 6.0); else tmp = ((1.0 + alpha) / beta) * (1.0 / (beta + 3.0)); end tmp_2 = tmp; end
code[alpha_, beta_] := If[LessEqual[beta, 5.3], N[(1.0 / N[(N[(N[(beta + 2.0), $MachinePrecision] + alpha), $MachinePrecision] * 6.0), $MachinePrecision]), $MachinePrecision], N[(N[(N[(1.0 + alpha), $MachinePrecision] / beta), $MachinePrecision] * N[(1.0 / N[(beta + 3.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\beta \leq 5.3:\\
\;\;\;\;\frac{1}{\left(\left(\beta + 2\right) + \alpha\right) \cdot 6}\\
\mathbf{else}:\\
\;\;\;\;\frac{1 + \alpha}{\beta} \cdot \frac{1}{\beta + 3}\\
\end{array}
\end{array}
if beta < 5.29999999999999982Initial program 99.8%
Simplified93.4%
associate-+r+93.4%
fma-undefine93.4%
*-commutative93.4%
associate-+l+93.4%
+-commutative93.4%
associate-+l+93.4%
*-commutative93.4%
associate-*r*93.4%
associate-+r+93.4%
+-commutative93.4%
associate-/l/99.9%
clear-num99.9%
inv-pow99.9%
Applied egg-rr99.9%
unpow-199.9%
associate-/r/99.9%
Simplified99.9%
Taylor expanded in alpha around 0 68.8%
Taylor expanded in beta around 0 68.0%
if 5.29999999999999982 < beta Initial program 88.9%
Taylor expanded in beta around inf 83.6%
div-inv83.5%
metadata-eval83.5%
associate-+l+83.5%
metadata-eval83.5%
associate-+r+83.5%
+-commutative83.5%
associate-+l+83.5%
Applied egg-rr83.5%
Taylor expanded in alpha around 0 83.3%
+-commutative83.3%
Simplified83.3%
Final simplification74.0%
(FPCore (alpha beta) :precision binary64 (if (<= beta 5.5) (/ 1.0 (* (+ (+ beta 2.0) alpha) 6.0)) (/ 1.0 (* beta (+ beta 3.0)))))
double code(double alpha, double beta) {
double tmp;
if (beta <= 5.5) {
tmp = 1.0 / (((beta + 2.0) + alpha) * 6.0);
} else {
tmp = 1.0 / (beta * (beta + 3.0));
}
return tmp;
}
real(8) function code(alpha, beta)
real(8), intent (in) :: alpha
real(8), intent (in) :: beta
real(8) :: tmp
if (beta <= 5.5d0) then
tmp = 1.0d0 / (((beta + 2.0d0) + alpha) * 6.0d0)
else
tmp = 1.0d0 / (beta * (beta + 3.0d0))
end if
code = tmp
end function
public static double code(double alpha, double beta) {
double tmp;
if (beta <= 5.5) {
tmp = 1.0 / (((beta + 2.0) + alpha) * 6.0);
} else {
tmp = 1.0 / (beta * (beta + 3.0));
}
return tmp;
}
def code(alpha, beta): tmp = 0 if beta <= 5.5: tmp = 1.0 / (((beta + 2.0) + alpha) * 6.0) else: tmp = 1.0 / (beta * (beta + 3.0)) return tmp
function code(alpha, beta) tmp = 0.0 if (beta <= 5.5) tmp = Float64(1.0 / Float64(Float64(Float64(beta + 2.0) + alpha) * 6.0)); else tmp = Float64(1.0 / Float64(beta * Float64(beta + 3.0))); end return tmp end
function tmp_2 = code(alpha, beta) tmp = 0.0; if (beta <= 5.5) tmp = 1.0 / (((beta + 2.0) + alpha) * 6.0); else tmp = 1.0 / (beta * (beta + 3.0)); end tmp_2 = tmp; end
code[alpha_, beta_] := If[LessEqual[beta, 5.5], N[(1.0 / N[(N[(N[(beta + 2.0), $MachinePrecision] + alpha), $MachinePrecision] * 6.0), $MachinePrecision]), $MachinePrecision], N[(1.0 / N[(beta * N[(beta + 3.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\beta \leq 5.5:\\
\;\;\;\;\frac{1}{\left(\left(\beta + 2\right) + \alpha\right) \cdot 6}\\
\mathbf{else}:\\
\;\;\;\;\frac{1}{\beta \cdot \left(\beta + 3\right)}\\
\end{array}
\end{array}
if beta < 5.5Initial program 99.8%
Simplified93.4%
associate-+r+93.4%
fma-undefine93.4%
*-commutative93.4%
associate-+l+93.4%
+-commutative93.4%
associate-+l+93.4%
*-commutative93.4%
associate-*r*93.4%
associate-+r+93.4%
+-commutative93.4%
associate-/l/99.9%
clear-num99.9%
inv-pow99.9%
Applied egg-rr99.9%
unpow-199.9%
associate-/r/99.9%
Simplified99.9%
Taylor expanded in alpha around 0 68.8%
Taylor expanded in beta around 0 68.0%
if 5.5 < beta Initial program 88.9%
Taylor expanded in beta around inf 83.6%
Taylor expanded in alpha around 0 79.4%
Final simplification72.5%
(FPCore (alpha beta) :precision binary64 (if (<= beta 1.96) (+ 0.08333333333333333 (* alpha -0.041666666666666664)) (/ 1.0 (* beta (+ beta 3.0)))))
double code(double alpha, double beta) {
double tmp;
if (beta <= 1.96) {
tmp = 0.08333333333333333 + (alpha * -0.041666666666666664);
} else {
tmp = 1.0 / (beta * (beta + 3.0));
}
return tmp;
}
real(8) function code(alpha, beta)
real(8), intent (in) :: alpha
real(8), intent (in) :: beta
real(8) :: tmp
if (beta <= 1.96d0) then
tmp = 0.08333333333333333d0 + (alpha * (-0.041666666666666664d0))
else
tmp = 1.0d0 / (beta * (beta + 3.0d0))
end if
code = tmp
end function
public static double code(double alpha, double beta) {
double tmp;
if (beta <= 1.96) {
tmp = 0.08333333333333333 + (alpha * -0.041666666666666664);
} else {
tmp = 1.0 / (beta * (beta + 3.0));
}
return tmp;
}
def code(alpha, beta): tmp = 0 if beta <= 1.96: tmp = 0.08333333333333333 + (alpha * -0.041666666666666664) else: tmp = 1.0 / (beta * (beta + 3.0)) return tmp
function code(alpha, beta) tmp = 0.0 if (beta <= 1.96) tmp = Float64(0.08333333333333333 + Float64(alpha * -0.041666666666666664)); else tmp = Float64(1.0 / Float64(beta * Float64(beta + 3.0))); end return tmp end
function tmp_2 = code(alpha, beta) tmp = 0.0; if (beta <= 1.96) tmp = 0.08333333333333333 + (alpha * -0.041666666666666664); else tmp = 1.0 / (beta * (beta + 3.0)); end tmp_2 = tmp; end
code[alpha_, beta_] := If[LessEqual[beta, 1.96], N[(0.08333333333333333 + N[(alpha * -0.041666666666666664), $MachinePrecision]), $MachinePrecision], N[(1.0 / N[(beta * N[(beta + 3.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\beta \leq 1.96:\\
\;\;\;\;0.08333333333333333 + \alpha \cdot -0.041666666666666664\\
\mathbf{else}:\\
\;\;\;\;\frac{1}{\beta \cdot \left(\beta + 3\right)}\\
\end{array}
\end{array}
if beta < 1.96Initial program 99.8%
Simplified93.4%
associate-+r+93.4%
fma-undefine93.4%
*-commutative93.4%
associate-+l+93.4%
+-commutative93.4%
associate-+l+93.4%
*-commutative93.4%
associate-*r*93.4%
associate-+r+93.4%
+-commutative93.4%
associate-/l/99.9%
clear-num99.9%
inv-pow99.9%
Applied egg-rr99.9%
unpow-199.9%
associate-/r/99.9%
Simplified99.9%
Taylor expanded in alpha around 0 68.8%
Taylor expanded in beta around 0 67.9%
Taylor expanded in alpha around 0 66.4%
*-commutative66.4%
Simplified66.4%
if 1.96 < beta Initial program 88.9%
Taylor expanded in beta around inf 83.6%
Taylor expanded in alpha around 0 79.4%
Final simplification71.5%
(FPCore (alpha beta) :precision binary64 (if (<= beta 6.6) (+ 0.08333333333333333 (* alpha -0.041666666666666664)) (/ 1.0 beta)))
double code(double alpha, double beta) {
double tmp;
if (beta <= 6.6) {
tmp = 0.08333333333333333 + (alpha * -0.041666666666666664);
} else {
tmp = 1.0 / beta;
}
return tmp;
}
real(8) function code(alpha, beta)
real(8), intent (in) :: alpha
real(8), intent (in) :: beta
real(8) :: tmp
if (beta <= 6.6d0) then
tmp = 0.08333333333333333d0 + (alpha * (-0.041666666666666664d0))
else
tmp = 1.0d0 / beta
end if
code = tmp
end function
public static double code(double alpha, double beta) {
double tmp;
if (beta <= 6.6) {
tmp = 0.08333333333333333 + (alpha * -0.041666666666666664);
} else {
tmp = 1.0 / beta;
}
return tmp;
}
def code(alpha, beta): tmp = 0 if beta <= 6.6: tmp = 0.08333333333333333 + (alpha * -0.041666666666666664) else: tmp = 1.0 / beta return tmp
function code(alpha, beta) tmp = 0.0 if (beta <= 6.6) tmp = Float64(0.08333333333333333 + Float64(alpha * -0.041666666666666664)); else tmp = Float64(1.0 / beta); end return tmp end
function tmp_2 = code(alpha, beta) tmp = 0.0; if (beta <= 6.6) tmp = 0.08333333333333333 + (alpha * -0.041666666666666664); else tmp = 1.0 / beta; end tmp_2 = tmp; end
code[alpha_, beta_] := If[LessEqual[beta, 6.6], N[(0.08333333333333333 + N[(alpha * -0.041666666666666664), $MachinePrecision]), $MachinePrecision], N[(1.0 / beta), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\beta \leq 6.6:\\
\;\;\;\;0.08333333333333333 + \alpha \cdot -0.041666666666666664\\
\mathbf{else}:\\
\;\;\;\;\frac{1}{\beta}\\
\end{array}
\end{array}
if beta < 6.5999999999999996Initial program 99.8%
Simplified93.4%
associate-+r+93.4%
fma-undefine93.4%
*-commutative93.4%
associate-+l+93.4%
+-commutative93.4%
associate-+l+93.4%
*-commutative93.4%
associate-*r*93.4%
associate-+r+93.4%
+-commutative93.4%
associate-/l/99.9%
clear-num99.9%
inv-pow99.9%
Applied egg-rr99.9%
unpow-199.9%
associate-/r/99.9%
Simplified99.9%
Taylor expanded in alpha around 0 68.8%
Taylor expanded in beta around 0 67.9%
Taylor expanded in alpha around 0 66.4%
*-commutative66.4%
Simplified66.4%
if 6.5999999999999996 < beta Initial program 88.9%
Taylor expanded in beta around inf 83.6%
Taylor expanded in alpha around inf 6.4%
(FPCore (alpha beta) :precision binary64 (if (<= beta 12.0) 0.08333333333333333 (/ 1.0 beta)))
double code(double alpha, double beta) {
double tmp;
if (beta <= 12.0) {
tmp = 0.08333333333333333;
} else {
tmp = 1.0 / beta;
}
return tmp;
}
real(8) function code(alpha, beta)
real(8), intent (in) :: alpha
real(8), intent (in) :: beta
real(8) :: tmp
if (beta <= 12.0d0) then
tmp = 0.08333333333333333d0
else
tmp = 1.0d0 / beta
end if
code = tmp
end function
public static double code(double alpha, double beta) {
double tmp;
if (beta <= 12.0) {
tmp = 0.08333333333333333;
} else {
tmp = 1.0 / beta;
}
return tmp;
}
def code(alpha, beta): tmp = 0 if beta <= 12.0: tmp = 0.08333333333333333 else: tmp = 1.0 / beta return tmp
function code(alpha, beta) tmp = 0.0 if (beta <= 12.0) tmp = 0.08333333333333333; else tmp = Float64(1.0 / beta); end return tmp end
function tmp_2 = code(alpha, beta) tmp = 0.0; if (beta <= 12.0) tmp = 0.08333333333333333; else tmp = 1.0 / beta; end tmp_2 = tmp; end
code[alpha_, beta_] := If[LessEqual[beta, 12.0], 0.08333333333333333, N[(1.0 / beta), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\beta \leq 12:\\
\;\;\;\;0.08333333333333333\\
\mathbf{else}:\\
\;\;\;\;\frac{1}{\beta}\\
\end{array}
\end{array}
if beta < 12Initial program 99.8%
Simplified93.4%
associate-+r+93.4%
fma-undefine93.4%
*-commutative93.4%
associate-+l+93.4%
+-commutative93.4%
associate-+l+93.4%
*-commutative93.4%
associate-*r*93.4%
associate-+r+93.4%
+-commutative93.4%
associate-/l/99.9%
clear-num99.9%
inv-pow99.9%
Applied egg-rr99.9%
unpow-199.9%
associate-/r/99.9%
Simplified99.9%
Taylor expanded in alpha around 0 68.8%
Taylor expanded in beta around 0 67.9%
Taylor expanded in alpha around 0 67.1%
if 12 < beta Initial program 88.9%
Taylor expanded in beta around inf 83.6%
Taylor expanded in alpha around inf 6.4%
(FPCore (alpha beta) :precision binary64 (/ 0.16666666666666666 (+ 2.0 alpha)))
double code(double alpha, double beta) {
return 0.16666666666666666 / (2.0 + alpha);
}
real(8) function code(alpha, beta)
real(8), intent (in) :: alpha
real(8), intent (in) :: beta
code = 0.16666666666666666d0 / (2.0d0 + alpha)
end function
public static double code(double alpha, double beta) {
return 0.16666666666666666 / (2.0 + alpha);
}
def code(alpha, beta): return 0.16666666666666666 / (2.0 + alpha)
function code(alpha, beta) return Float64(0.16666666666666666 / Float64(2.0 + alpha)) end
function tmp = code(alpha, beta) tmp = 0.16666666666666666 / (2.0 + alpha); end
code[alpha_, beta_] := N[(0.16666666666666666 / N[(2.0 + alpha), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{0.16666666666666666}{2 + \alpha}
\end{array}
Initial program 95.6%
Simplified83.0%
associate-+r+83.0%
fma-undefine83.0%
*-commutative83.0%
associate-+l+83.0%
+-commutative83.0%
associate-+l+83.0%
*-commutative83.0%
associate-*r*83.0%
associate-+r+83.0%
+-commutative83.0%
associate-/l/93.8%
clear-num93.8%
inv-pow93.8%
Applied egg-rr93.8%
unpow-193.8%
associate-/r/93.8%
Simplified99.8%
Taylor expanded in alpha around 0 75.6%
Taylor expanded in beta around 0 43.2%
(FPCore (alpha beta) :precision binary64 0.08333333333333333)
double code(double alpha, double beta) {
return 0.08333333333333333;
}
real(8) function code(alpha, beta)
real(8), intent (in) :: alpha
real(8), intent (in) :: beta
code = 0.08333333333333333d0
end function
public static double code(double alpha, double beta) {
return 0.08333333333333333;
}
def code(alpha, beta): return 0.08333333333333333
function code(alpha, beta) return 0.08333333333333333 end
function tmp = code(alpha, beta) tmp = 0.08333333333333333; end
code[alpha_, beta_] := 0.08333333333333333
\begin{array}{l}
\\
0.08333333333333333
\end{array}
Initial program 95.6%
Simplified83.0%
associate-+r+83.0%
fma-undefine83.0%
*-commutative83.0%
associate-+l+83.0%
+-commutative83.0%
associate-+l+83.0%
*-commutative83.0%
associate-*r*83.0%
associate-+r+83.0%
+-commutative83.0%
associate-/l/93.8%
clear-num93.8%
inv-pow93.8%
Applied egg-rr93.8%
unpow-193.8%
associate-/r/93.8%
Simplified99.8%
Taylor expanded in alpha around 0 75.6%
Taylor expanded in beta around 0 43.2%
Taylor expanded in alpha around 0 42.4%
herbie shell --seed 2024103
(FPCore (alpha beta)
:name "Octave 3.8, jcobi/3"
:precision binary64
:pre (and (> alpha -1.0) (> beta -1.0))
(/ (/ (/ (+ (+ (+ alpha beta) (* beta alpha)) 1.0) (+ (+ alpha beta) (* 2.0 1.0))) (+ (+ alpha beta) (* 2.0 1.0))) (+ (+ (+ alpha beta) (* 2.0 1.0)) 1.0)))