
(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 21 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}
NOTE: alpha and beta should be sorted in increasing order before calling this function.
(FPCore (alpha beta)
:precision binary64
(if (<= beta 2e+153)
(/
(*
(pow (+ 2.0 (+ alpha beta)) -2.0)
(+ 1.0 (fma beta alpha (+ alpha beta))))
(+ 3.0 (+ alpha beta)))
(/
(/
(-
(+ (/ alpha beta) (+ (+ (/ 1.0 beta) alpha) 1.0))
(* (/ (fma 2.0 alpha 5.0) beta) (+ 1.0 alpha)))
beta)
(+ (+ 2.0 beta) alpha))))assert(alpha < beta);
double code(double alpha, double beta) {
double tmp;
if (beta <= 2e+153) {
tmp = (pow((2.0 + (alpha + beta)), -2.0) * (1.0 + fma(beta, alpha, (alpha + beta)))) / (3.0 + (alpha + beta));
} else {
tmp = ((((alpha / beta) + (((1.0 / beta) + alpha) + 1.0)) - ((fma(2.0, alpha, 5.0) / beta) * (1.0 + alpha))) / beta) / ((2.0 + beta) + alpha);
}
return tmp;
}
alpha, beta = sort([alpha, beta]) function code(alpha, beta) tmp = 0.0 if (beta <= 2e+153) tmp = Float64(Float64((Float64(2.0 + Float64(alpha + beta)) ^ -2.0) * Float64(1.0 + fma(beta, alpha, Float64(alpha + beta)))) / Float64(3.0 + Float64(alpha + beta))); else tmp = Float64(Float64(Float64(Float64(Float64(alpha / beta) + Float64(Float64(Float64(1.0 / beta) + alpha) + 1.0)) - Float64(Float64(fma(2.0, alpha, 5.0) / beta) * Float64(1.0 + alpha))) / beta) / Float64(Float64(2.0 + beta) + alpha)); end return tmp end
NOTE: alpha and beta should be sorted in increasing order before calling this function. code[alpha_, beta_] := If[LessEqual[beta, 2e+153], N[(N[(N[Power[N[(2.0 + N[(alpha + beta), $MachinePrecision]), $MachinePrecision], -2.0], $MachinePrecision] * N[(1.0 + N[(beta * alpha + N[(alpha + beta), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(3.0 + N[(alpha + beta), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(N[(N[(N[(alpha / beta), $MachinePrecision] + N[(N[(N[(1.0 / beta), $MachinePrecision] + alpha), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision] - N[(N[(N[(2.0 * alpha + 5.0), $MachinePrecision] / beta), $MachinePrecision] * N[(1.0 + alpha), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / beta), $MachinePrecision] / N[(N[(2.0 + beta), $MachinePrecision] + alpha), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
[alpha, beta] = \mathsf{sort}([alpha, beta])\\
\\
\begin{array}{l}
\mathbf{if}\;\beta \leq 2 \cdot 10^{+153}:\\
\;\;\;\;\frac{{\left(2 + \left(\alpha + \beta\right)\right)}^{-2} \cdot \left(1 + \mathsf{fma}\left(\beta, \alpha, \alpha + \beta\right)\right)}{3 + \left(\alpha + \beta\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{\left(\frac{\alpha}{\beta} + \left(\left(\frac{1}{\beta} + \alpha\right) + 1\right)\right) - \frac{\mathsf{fma}\left(2, \alpha, 5\right)}{\beta} \cdot \left(1 + \alpha\right)}{\beta}}{\left(2 + \beta\right) + \alpha}\\
\end{array}
\end{array}
if beta < 2e153Initial program 96.7%
Applied rewrites95.9%
if 2e153 < beta Initial program 72.7%
Taylor expanded in beta around -inf
mul-1-negN/A
lower-neg.f64N/A
sub-negN/A
mul-1-negN/A
distribute-neg-inN/A
+-commutativeN/A
distribute-neg-inN/A
metadata-evalN/A
unsub-negN/A
lower--.f6486.6
Applied rewrites86.6%
lift-/.f64N/A
lift-/.f64N/A
associate-/l/N/A
associate-/r*N/A
Applied rewrites86.6%
Taylor expanded in beta around inf
lower-/.f64N/A
lower--.f64N/A
associate-+r+N/A
associate-+r+N/A
lower-+.f64N/A
lower-+.f64N/A
lower-+.f64N/A
lower-/.f64N/A
lower-/.f64N/A
associate-/l*N/A
lower-*.f64N/A
lower-+.f64N/A
lower-/.f64N/A
+-commutativeN/A
lower-fma.f6486.4
Applied rewrites86.4%
Final simplification94.6%
NOTE: alpha and beta should be sorted in increasing order before calling this function.
(FPCore (alpha beta)
:precision binary64
(let* ((t_0 (+ 2.0 (+ alpha beta)))
(t_1 (/ (/ (+ (+ (* alpha beta) (+ alpha beta)) 1.0) t_0) t_0)))
(if (<= (/ t_1 (+ t_0 1.0)) 0.1)
(/ t_1 (+ (+ 1.0 (+ alpha beta)) 2.0))
(/ (/ (+ 1.0 alpha) (+ (+ 2.0 beta) alpha)) (+ 3.0 (+ alpha beta))))))assert(alpha < beta);
double code(double alpha, double beta) {
double t_0 = 2.0 + (alpha + beta);
double t_1 = ((((alpha * beta) + (alpha + beta)) + 1.0) / t_0) / t_0;
double tmp;
if ((t_1 / (t_0 + 1.0)) <= 0.1) {
tmp = t_1 / ((1.0 + (alpha + beta)) + 2.0);
} else {
tmp = ((1.0 + alpha) / ((2.0 + beta) + alpha)) / (3.0 + (alpha + beta));
}
return tmp;
}
NOTE: alpha and beta should be sorted in increasing order before calling this function.
real(8) function code(alpha, beta)
real(8), intent (in) :: alpha
real(8), intent (in) :: beta
real(8) :: t_0
real(8) :: t_1
real(8) :: tmp
t_0 = 2.0d0 + (alpha + beta)
t_1 = ((((alpha * beta) + (alpha + beta)) + 1.0d0) / t_0) / t_0
if ((t_1 / (t_0 + 1.0d0)) <= 0.1d0) then
tmp = t_1 / ((1.0d0 + (alpha + beta)) + 2.0d0)
else
tmp = ((1.0d0 + alpha) / ((2.0d0 + beta) + alpha)) / (3.0d0 + (alpha + beta))
end if
code = tmp
end function
assert alpha < beta;
public static double code(double alpha, double beta) {
double t_0 = 2.0 + (alpha + beta);
double t_1 = ((((alpha * beta) + (alpha + beta)) + 1.0) / t_0) / t_0;
double tmp;
if ((t_1 / (t_0 + 1.0)) <= 0.1) {
tmp = t_1 / ((1.0 + (alpha + beta)) + 2.0);
} else {
tmp = ((1.0 + alpha) / ((2.0 + beta) + alpha)) / (3.0 + (alpha + beta));
}
return tmp;
}
[alpha, beta] = sort([alpha, beta]) def code(alpha, beta): t_0 = 2.0 + (alpha + beta) t_1 = ((((alpha * beta) + (alpha + beta)) + 1.0) / t_0) / t_0 tmp = 0 if (t_1 / (t_0 + 1.0)) <= 0.1: tmp = t_1 / ((1.0 + (alpha + beta)) + 2.0) else: tmp = ((1.0 + alpha) / ((2.0 + beta) + alpha)) / (3.0 + (alpha + beta)) return tmp
alpha, beta = sort([alpha, beta]) function code(alpha, beta) t_0 = Float64(2.0 + Float64(alpha + beta)) t_1 = Float64(Float64(Float64(Float64(Float64(alpha * beta) + Float64(alpha + beta)) + 1.0) / t_0) / t_0) tmp = 0.0 if (Float64(t_1 / Float64(t_0 + 1.0)) <= 0.1) tmp = Float64(t_1 / Float64(Float64(1.0 + Float64(alpha + beta)) + 2.0)); else tmp = Float64(Float64(Float64(1.0 + alpha) / Float64(Float64(2.0 + beta) + alpha)) / Float64(3.0 + Float64(alpha + beta))); end return tmp end
alpha, beta = num2cell(sort([alpha, beta])){:}
function tmp_2 = code(alpha, beta)
t_0 = 2.0 + (alpha + beta);
t_1 = ((((alpha * beta) + (alpha + beta)) + 1.0) / t_0) / t_0;
tmp = 0.0;
if ((t_1 / (t_0 + 1.0)) <= 0.1)
tmp = t_1 / ((1.0 + (alpha + beta)) + 2.0);
else
tmp = ((1.0 + alpha) / ((2.0 + beta) + alpha)) / (3.0 + (alpha + beta));
end
tmp_2 = tmp;
end
NOTE: alpha and beta should be sorted in increasing order before calling this function.
code[alpha_, beta_] := Block[{t$95$0 = N[(2.0 + N[(alpha + beta), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(N[(N[(N[(N[(alpha * beta), $MachinePrecision] + N[(alpha + beta), $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision] / t$95$0), $MachinePrecision] / t$95$0), $MachinePrecision]}, If[LessEqual[N[(t$95$1 / N[(t$95$0 + 1.0), $MachinePrecision]), $MachinePrecision], 0.1], N[(t$95$1 / N[(N[(1.0 + N[(alpha + beta), $MachinePrecision]), $MachinePrecision] + 2.0), $MachinePrecision]), $MachinePrecision], N[(N[(N[(1.0 + alpha), $MachinePrecision] / N[(N[(2.0 + beta), $MachinePrecision] + alpha), $MachinePrecision]), $MachinePrecision] / N[(3.0 + N[(alpha + beta), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
[alpha, beta] = \mathsf{sort}([alpha, beta])\\
\\
\begin{array}{l}
t_0 := 2 + \left(\alpha + \beta\right)\\
t_1 := \frac{\frac{\left(\alpha \cdot \beta + \left(\alpha + \beta\right)\right) + 1}{t\_0}}{t\_0}\\
\mathbf{if}\;\frac{t\_1}{t\_0 + 1} \leq 0.1:\\
\;\;\;\;\frac{t\_1}{\left(1 + \left(\alpha + \beta\right)\right) + 2}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{1 + \alpha}{\left(2 + \beta\right) + \alpha}}{3 + \left(\alpha + \beta\right)}\\
\end{array}
\end{array}
if (/.f64 (/.f64 (/.f64 (+.f64 (+.f64 (+.f64 alpha beta) (*.f64 beta alpha)) #s(literal 1 binary64)) (+.f64 (+.f64 alpha beta) (*.f64 #s(literal 2 binary64) #s(literal 1 binary64)))) (+.f64 (+.f64 alpha beta) (*.f64 #s(literal 2 binary64) #s(literal 1 binary64)))) (+.f64 (+.f64 (+.f64 alpha beta) (*.f64 #s(literal 2 binary64) #s(literal 1 binary64))) #s(literal 1 binary64))) < 0.10000000000000001Initial program 99.8%
lift-+.f64N/A
+-commutativeN/A
lift-+.f64N/A
associate-+r+N/A
lower-+.f64N/A
lower-+.f6499.9
lift-+.f64N/A
+-commutativeN/A
lower-+.f6499.9
lift-*.f64N/A
metadata-eval99.9
Applied rewrites99.9%
if 0.10000000000000001 < (/.f64 (/.f64 (/.f64 (+.f64 (+.f64 (+.f64 alpha beta) (*.f64 beta alpha)) #s(literal 1 binary64)) (+.f64 (+.f64 alpha beta) (*.f64 #s(literal 2 binary64) #s(literal 1 binary64)))) (+.f64 (+.f64 alpha beta) (*.f64 #s(literal 2 binary64) #s(literal 1 binary64)))) (+.f64 (+.f64 (+.f64 alpha beta) (*.f64 #s(literal 2 binary64) #s(literal 1 binary64))) #s(literal 1 binary64))) Initial program 1.7%
Taylor expanded in beta around -inf
mul-1-negN/A
lower-neg.f64N/A
sub-negN/A
mul-1-negN/A
distribute-neg-inN/A
+-commutativeN/A
distribute-neg-inN/A
metadata-evalN/A
unsub-negN/A
lower--.f6434.0
Applied rewrites34.0%
Applied rewrites34.0%
Final simplification95.5%
NOTE: alpha and beta should be sorted in increasing order before calling this function.
(FPCore (alpha beta)
:precision binary64
(let* ((t_0 (+ 2.0 (+ alpha beta))))
(if (<= beta 2e+135)
(/
(/ (/ (+ (+ (* alpha beta) (+ alpha beta)) 1.0) t_0) t_0)
(+ (+ 1.0 (+ alpha beta)) 2.0))
(/
(/
(-
(+ (/ alpha beta) (+ (+ (/ 1.0 beta) alpha) 1.0))
(* (/ (fma 2.0 alpha 5.0) beta) (+ 1.0 alpha)))
beta)
(+ (+ 2.0 beta) alpha)))))assert(alpha < beta);
double code(double alpha, double beta) {
double t_0 = 2.0 + (alpha + beta);
double tmp;
if (beta <= 2e+135) {
tmp = (((((alpha * beta) + (alpha + beta)) + 1.0) / t_0) / t_0) / ((1.0 + (alpha + beta)) + 2.0);
} else {
tmp = ((((alpha / beta) + (((1.0 / beta) + alpha) + 1.0)) - ((fma(2.0, alpha, 5.0) / beta) * (1.0 + alpha))) / beta) / ((2.0 + beta) + alpha);
}
return tmp;
}
alpha, beta = sort([alpha, beta]) function code(alpha, beta) t_0 = Float64(2.0 + Float64(alpha + beta)) tmp = 0.0 if (beta <= 2e+135) tmp = Float64(Float64(Float64(Float64(Float64(Float64(alpha * beta) + Float64(alpha + beta)) + 1.0) / t_0) / t_0) / Float64(Float64(1.0 + Float64(alpha + beta)) + 2.0)); else tmp = Float64(Float64(Float64(Float64(Float64(alpha / beta) + Float64(Float64(Float64(1.0 / beta) + alpha) + 1.0)) - Float64(Float64(fma(2.0, alpha, 5.0) / beta) * Float64(1.0 + alpha))) / beta) / Float64(Float64(2.0 + beta) + alpha)); end return tmp end
NOTE: alpha and beta should be sorted in increasing order before calling this function.
code[alpha_, beta_] := Block[{t$95$0 = N[(2.0 + N[(alpha + beta), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[beta, 2e+135], N[(N[(N[(N[(N[(N[(alpha * beta), $MachinePrecision] + N[(alpha + beta), $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision] / t$95$0), $MachinePrecision] / t$95$0), $MachinePrecision] / N[(N[(1.0 + N[(alpha + beta), $MachinePrecision]), $MachinePrecision] + 2.0), $MachinePrecision]), $MachinePrecision], N[(N[(N[(N[(N[(alpha / beta), $MachinePrecision] + N[(N[(N[(1.0 / beta), $MachinePrecision] + alpha), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision] - N[(N[(N[(2.0 * alpha + 5.0), $MachinePrecision] / beta), $MachinePrecision] * N[(1.0 + alpha), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / beta), $MachinePrecision] / N[(N[(2.0 + beta), $MachinePrecision] + alpha), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
[alpha, beta] = \mathsf{sort}([alpha, beta])\\
\\
\begin{array}{l}
t_0 := 2 + \left(\alpha + \beta\right)\\
\mathbf{if}\;\beta \leq 2 \cdot 10^{+135}:\\
\;\;\;\;\frac{\frac{\frac{\left(\alpha \cdot \beta + \left(\alpha + \beta\right)\right) + 1}{t\_0}}{t\_0}}{\left(1 + \left(\alpha + \beta\right)\right) + 2}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{\left(\frac{\alpha}{\beta} + \left(\left(\frac{1}{\beta} + \alpha\right) + 1\right)\right) - \frac{\mathsf{fma}\left(2, \alpha, 5\right)}{\beta} \cdot \left(1 + \alpha\right)}{\beta}}{\left(2 + \beta\right) + \alpha}\\
\end{array}
\end{array}
if beta < 1.99999999999999992e135Initial program 97.1%
lift-+.f64N/A
+-commutativeN/A
lift-+.f64N/A
associate-+r+N/A
lower-+.f64N/A
lower-+.f6497.1
lift-+.f64N/A
+-commutativeN/A
lower-+.f6497.1
lift-*.f64N/A
metadata-eval97.1
Applied rewrites97.1%
if 1.99999999999999992e135 < beta Initial program 72.2%
Taylor expanded in beta around -inf
mul-1-negN/A
lower-neg.f64N/A
sub-negN/A
mul-1-negN/A
distribute-neg-inN/A
+-commutativeN/A
distribute-neg-inN/A
metadata-evalN/A
unsub-negN/A
lower--.f6485.1
Applied rewrites85.1%
lift-/.f64N/A
lift-/.f64N/A
associate-/l/N/A
associate-/r*N/A
Applied rewrites85.1%
Taylor expanded in beta around inf
lower-/.f64N/A
lower--.f64N/A
associate-+r+N/A
associate-+r+N/A
lower-+.f64N/A
lower-+.f64N/A
lower-+.f64N/A
lower-/.f64N/A
lower-/.f64N/A
associate-/l*N/A
lower-*.f64N/A
lower-+.f64N/A
lower-/.f64N/A
+-commutativeN/A
lower-fma.f6484.9
Applied rewrites84.9%
Final simplification95.3%
NOTE: alpha and beta should be sorted in increasing order before calling this function.
(FPCore (alpha beta)
:precision binary64
(if (<= beta 2e+144)
(/
(/ (* (+ 1.0 beta) (+ 1.0 alpha)) (+ 2.0 (+ alpha beta)))
(fma (+ (fma 2.0 alpha beta) 5.0) beta (* (+ 2.0 alpha) (+ 3.0 alpha))))
(/ (/ (+ 1.0 alpha) (+ (+ 2.0 beta) alpha)) (+ 3.0 (+ alpha beta)))))assert(alpha < beta);
double code(double alpha, double beta) {
double tmp;
if (beta <= 2e+144) {
tmp = (((1.0 + beta) * (1.0 + alpha)) / (2.0 + (alpha + beta))) / fma((fma(2.0, alpha, beta) + 5.0), beta, ((2.0 + alpha) * (3.0 + alpha)));
} else {
tmp = ((1.0 + alpha) / ((2.0 + beta) + alpha)) / (3.0 + (alpha + beta));
}
return tmp;
}
alpha, beta = sort([alpha, beta]) function code(alpha, beta) tmp = 0.0 if (beta <= 2e+144) tmp = Float64(Float64(Float64(Float64(1.0 + beta) * Float64(1.0 + alpha)) / Float64(2.0 + Float64(alpha + beta))) / fma(Float64(fma(2.0, alpha, beta) + 5.0), beta, Float64(Float64(2.0 + alpha) * Float64(3.0 + alpha)))); else tmp = Float64(Float64(Float64(1.0 + alpha) / Float64(Float64(2.0 + beta) + alpha)) / Float64(3.0 + Float64(alpha + beta))); end return tmp end
NOTE: alpha and beta should be sorted in increasing order before calling this function. code[alpha_, beta_] := If[LessEqual[beta, 2e+144], N[(N[(N[(N[(1.0 + beta), $MachinePrecision] * N[(1.0 + alpha), $MachinePrecision]), $MachinePrecision] / N[(2.0 + N[(alpha + beta), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(N[(N[(2.0 * alpha + beta), $MachinePrecision] + 5.0), $MachinePrecision] * beta + N[(N[(2.0 + alpha), $MachinePrecision] * N[(3.0 + alpha), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(N[(1.0 + alpha), $MachinePrecision] / N[(N[(2.0 + beta), $MachinePrecision] + alpha), $MachinePrecision]), $MachinePrecision] / N[(3.0 + N[(alpha + beta), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
[alpha, beta] = \mathsf{sort}([alpha, beta])\\
\\
\begin{array}{l}
\mathbf{if}\;\beta \leq 2 \cdot 10^{+144}:\\
\;\;\;\;\frac{\frac{\left(1 + \beta\right) \cdot \left(1 + \alpha\right)}{2 + \left(\alpha + \beta\right)}}{\mathsf{fma}\left(\mathsf{fma}\left(2, \alpha, \beta\right) + 5, \beta, \left(2 + \alpha\right) \cdot \left(3 + \alpha\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{1 + \alpha}{\left(2 + \beta\right) + \alpha}}{3 + \left(\alpha + \beta\right)}\\
\end{array}
\end{array}
if beta < 2.00000000000000005e144Initial program 97.1%
lift-/.f64N/A
clear-numN/A
lower-/.f64N/A
lower-/.f6497.1
lift-+.f64N/A
+-commutativeN/A
lower-+.f6497.1
lift-*.f64N/A
metadata-eval97.1
lift-+.f64N/A
+-commutativeN/A
lift-*.f64N/A
lower-fma.f6497.1
lift-+.f64N/A
+-commutativeN/A
lower-+.f6497.1
Applied rewrites97.1%
Applied rewrites88.9%
Taylor expanded in beta around 0
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-+.f64N/A
+-commutativeN/A
lower-fma.f64N/A
lower-*.f64N/A
lower-+.f64N/A
lower-+.f6489.0
Applied rewrites89.0%
Applied rewrites96.3%
if 2.00000000000000005e144 < beta Initial program 71.5%
Taylor expanded in beta around -inf
mul-1-negN/A
lower-neg.f64N/A
sub-negN/A
mul-1-negN/A
distribute-neg-inN/A
+-commutativeN/A
distribute-neg-inN/A
metadata-evalN/A
unsub-negN/A
lower--.f6484.8
Applied rewrites84.8%
Applied rewrites84.8%
Final simplification94.6%
NOTE: alpha and beta should be sorted in increasing order before calling this function.
(FPCore (alpha beta)
:precision binary64
(let* ((t_0 (+ 3.0 (+ alpha beta))) (t_1 (+ 2.0 (+ alpha beta))))
(if (<= beta 2e+135)
(/ (/ (/ (+ 1.0 (fma beta alpha (+ alpha beta))) t_1) t_0) t_1)
(/ (/ (+ 1.0 alpha) (+ (+ 2.0 beta) alpha)) t_0))))assert(alpha < beta);
double code(double alpha, double beta) {
double t_0 = 3.0 + (alpha + beta);
double t_1 = 2.0 + (alpha + beta);
double tmp;
if (beta <= 2e+135) {
tmp = (((1.0 + fma(beta, alpha, (alpha + beta))) / t_1) / t_0) / t_1;
} else {
tmp = ((1.0 + alpha) / ((2.0 + beta) + alpha)) / t_0;
}
return tmp;
}
alpha, beta = sort([alpha, beta]) function code(alpha, beta) t_0 = Float64(3.0 + Float64(alpha + beta)) t_1 = Float64(2.0 + Float64(alpha + beta)) tmp = 0.0 if (beta <= 2e+135) tmp = Float64(Float64(Float64(Float64(1.0 + fma(beta, alpha, Float64(alpha + beta))) / t_1) / t_0) / t_1); else tmp = Float64(Float64(Float64(1.0 + alpha) / Float64(Float64(2.0 + beta) + alpha)) / t_0); end return tmp end
NOTE: alpha and beta should be sorted in increasing order before calling this function.
code[alpha_, beta_] := Block[{t$95$0 = N[(3.0 + N[(alpha + beta), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(2.0 + N[(alpha + beta), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[beta, 2e+135], N[(N[(N[(N[(1.0 + N[(beta * alpha + N[(alpha + beta), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / t$95$1), $MachinePrecision] / t$95$0), $MachinePrecision] / t$95$1), $MachinePrecision], N[(N[(N[(1.0 + alpha), $MachinePrecision] / N[(N[(2.0 + beta), $MachinePrecision] + alpha), $MachinePrecision]), $MachinePrecision] / t$95$0), $MachinePrecision]]]]
\begin{array}{l}
[alpha, beta] = \mathsf{sort}([alpha, beta])\\
\\
\begin{array}{l}
t_0 := 3 + \left(\alpha + \beta\right)\\
t_1 := 2 + \left(\alpha + \beta\right)\\
\mathbf{if}\;\beta \leq 2 \cdot 10^{+135}:\\
\;\;\;\;\frac{\frac{\frac{1 + \mathsf{fma}\left(\beta, \alpha, \alpha + \beta\right)}{t\_1}}{t\_0}}{t\_1}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{1 + \alpha}{\left(2 + \beta\right) + \alpha}}{t\_0}\\
\end{array}
\end{array}
if beta < 1.99999999999999992e135Initial program 97.1%
lift-/.f64N/A
lift-/.f64N/A
associate-/l/N/A
associate-/r*N/A
lower-/.f64N/A
Applied rewrites97.1%
if 1.99999999999999992e135 < beta Initial program 72.2%
Taylor expanded in beta around -inf
mul-1-negN/A
lower-neg.f64N/A
sub-negN/A
mul-1-negN/A
distribute-neg-inN/A
+-commutativeN/A
distribute-neg-inN/A
metadata-evalN/A
unsub-negN/A
lower--.f6485.1
Applied rewrites85.1%
Applied rewrites85.1%
Final simplification95.3%
NOTE: alpha and beta should be sorted in increasing order before calling this function.
(FPCore (alpha beta)
:precision binary64
(let* ((t_0 (+ 3.0 (+ alpha beta))) (t_1 (+ 2.0 (+ alpha beta))))
(if (<= beta 5e+143)
(/ 1.0 (* (/ t_1 (+ 1.0 (fma beta alpha (+ alpha beta)))) (* t_0 t_1)))
(/ (/ (+ 1.0 alpha) (+ (+ 2.0 beta) alpha)) t_0))))assert(alpha < beta);
double code(double alpha, double beta) {
double t_0 = 3.0 + (alpha + beta);
double t_1 = 2.0 + (alpha + beta);
double tmp;
if (beta <= 5e+143) {
tmp = 1.0 / ((t_1 / (1.0 + fma(beta, alpha, (alpha + beta)))) * (t_0 * t_1));
} else {
tmp = ((1.0 + alpha) / ((2.0 + beta) + alpha)) / t_0;
}
return tmp;
}
alpha, beta = sort([alpha, beta]) function code(alpha, beta) t_0 = Float64(3.0 + Float64(alpha + beta)) t_1 = Float64(2.0 + Float64(alpha + beta)) tmp = 0.0 if (beta <= 5e+143) tmp = Float64(1.0 / Float64(Float64(t_1 / Float64(1.0 + fma(beta, alpha, Float64(alpha + beta)))) * Float64(t_0 * t_1))); else tmp = Float64(Float64(Float64(1.0 + alpha) / Float64(Float64(2.0 + beta) + alpha)) / t_0); end return tmp end
NOTE: alpha and beta should be sorted in increasing order before calling this function.
code[alpha_, beta_] := Block[{t$95$0 = N[(3.0 + N[(alpha + beta), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(2.0 + N[(alpha + beta), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[beta, 5e+143], N[(1.0 / N[(N[(t$95$1 / N[(1.0 + N[(beta * alpha + N[(alpha + beta), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(t$95$0 * t$95$1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(N[(1.0 + alpha), $MachinePrecision] / N[(N[(2.0 + beta), $MachinePrecision] + alpha), $MachinePrecision]), $MachinePrecision] / t$95$0), $MachinePrecision]]]]
\begin{array}{l}
[alpha, beta] = \mathsf{sort}([alpha, beta])\\
\\
\begin{array}{l}
t_0 := 3 + \left(\alpha + \beta\right)\\
t_1 := 2 + \left(\alpha + \beta\right)\\
\mathbf{if}\;\beta \leq 5 \cdot 10^{+143}:\\
\;\;\;\;\frac{1}{\frac{t\_1}{1 + \mathsf{fma}\left(\beta, \alpha, \alpha + \beta\right)} \cdot \left(t\_0 \cdot t\_1\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{1 + \alpha}{\left(2 + \beta\right) + \alpha}}{t\_0}\\
\end{array}
\end{array}
if beta < 5.00000000000000012e143Initial program 97.1%
lift-/.f64N/A
lift-/.f64N/A
associate-/l/N/A
lift-/.f64N/A
clear-numN/A
associate-/l/N/A
lower-/.f64N/A
lower-*.f64N/A
Applied rewrites96.2%
if 5.00000000000000012e143 < beta Initial program 71.5%
Taylor expanded in beta around -inf
mul-1-negN/A
lower-neg.f64N/A
sub-negN/A
mul-1-negN/A
distribute-neg-inN/A
+-commutativeN/A
distribute-neg-inN/A
metadata-evalN/A
unsub-negN/A
lower--.f6484.8
Applied rewrites84.8%
Applied rewrites84.8%
Final simplification94.5%
NOTE: alpha and beta should be sorted in increasing order before calling this function.
(FPCore (alpha beta)
:precision binary64
(let* ((t_0 (+ 3.0 (+ alpha beta))) (t_1 (+ (+ 2.0 beta) alpha)))
(if (<= beta 2e+144)
(/ (/ (* (+ 1.0 beta) (+ 1.0 alpha)) t_1) (* t_1 t_0))
(/ (/ (+ 1.0 alpha) t_1) t_0))))assert(alpha < beta);
double code(double alpha, double beta) {
double t_0 = 3.0 + (alpha + beta);
double t_1 = (2.0 + beta) + alpha;
double tmp;
if (beta <= 2e+144) {
tmp = (((1.0 + beta) * (1.0 + alpha)) / t_1) / (t_1 * t_0);
} else {
tmp = ((1.0 + alpha) / t_1) / t_0;
}
return tmp;
}
NOTE: alpha and beta should be sorted in increasing order before calling this function.
real(8) function code(alpha, beta)
real(8), intent (in) :: alpha
real(8), intent (in) :: beta
real(8) :: t_0
real(8) :: t_1
real(8) :: tmp
t_0 = 3.0d0 + (alpha + beta)
t_1 = (2.0d0 + beta) + alpha
if (beta <= 2d+144) then
tmp = (((1.0d0 + beta) * (1.0d0 + alpha)) / t_1) / (t_1 * t_0)
else
tmp = ((1.0d0 + alpha) / t_1) / t_0
end if
code = tmp
end function
assert alpha < beta;
public static double code(double alpha, double beta) {
double t_0 = 3.0 + (alpha + beta);
double t_1 = (2.0 + beta) + alpha;
double tmp;
if (beta <= 2e+144) {
tmp = (((1.0 + beta) * (1.0 + alpha)) / t_1) / (t_1 * t_0);
} else {
tmp = ((1.0 + alpha) / t_1) / t_0;
}
return tmp;
}
[alpha, beta] = sort([alpha, beta]) def code(alpha, beta): t_0 = 3.0 + (alpha + beta) t_1 = (2.0 + beta) + alpha tmp = 0 if beta <= 2e+144: tmp = (((1.0 + beta) * (1.0 + alpha)) / t_1) / (t_1 * t_0) else: tmp = ((1.0 + alpha) / t_1) / t_0 return tmp
alpha, beta = sort([alpha, beta]) function code(alpha, beta) t_0 = Float64(3.0 + Float64(alpha + beta)) t_1 = Float64(Float64(2.0 + beta) + alpha) tmp = 0.0 if (beta <= 2e+144) tmp = Float64(Float64(Float64(Float64(1.0 + beta) * Float64(1.0 + alpha)) / t_1) / Float64(t_1 * t_0)); else tmp = Float64(Float64(Float64(1.0 + alpha) / t_1) / t_0); end return tmp end
alpha, beta = num2cell(sort([alpha, beta])){:}
function tmp_2 = code(alpha, beta)
t_0 = 3.0 + (alpha + beta);
t_1 = (2.0 + beta) + alpha;
tmp = 0.0;
if (beta <= 2e+144)
tmp = (((1.0 + beta) * (1.0 + alpha)) / t_1) / (t_1 * t_0);
else
tmp = ((1.0 + alpha) / t_1) / t_0;
end
tmp_2 = tmp;
end
NOTE: alpha and beta should be sorted in increasing order before calling this function.
code[alpha_, beta_] := Block[{t$95$0 = N[(3.0 + N[(alpha + beta), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(N[(2.0 + beta), $MachinePrecision] + alpha), $MachinePrecision]}, If[LessEqual[beta, 2e+144], N[(N[(N[(N[(1.0 + beta), $MachinePrecision] * N[(1.0 + alpha), $MachinePrecision]), $MachinePrecision] / t$95$1), $MachinePrecision] / N[(t$95$1 * t$95$0), $MachinePrecision]), $MachinePrecision], N[(N[(N[(1.0 + alpha), $MachinePrecision] / t$95$1), $MachinePrecision] / t$95$0), $MachinePrecision]]]]
\begin{array}{l}
[alpha, beta] = \mathsf{sort}([alpha, beta])\\
\\
\begin{array}{l}
t_0 := 3 + \left(\alpha + \beta\right)\\
t_1 := \left(2 + \beta\right) + \alpha\\
\mathbf{if}\;\beta \leq 2 \cdot 10^{+144}:\\
\;\;\;\;\frac{\frac{\left(1 + \beta\right) \cdot \left(1 + \alpha\right)}{t\_1}}{t\_1 \cdot t\_0}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{1 + \alpha}{t\_1}}{t\_0}\\
\end{array}
\end{array}
if beta < 2.00000000000000005e144Initial program 97.1%
lift-+.f64N/A
+-commutativeN/A
lift-+.f64N/A
associate-+r+N/A
lower-+.f64N/A
lower-+.f6497.1
lift-+.f64N/A
+-commutativeN/A
lower-+.f6497.1
lift-*.f64N/A
metadata-eval97.1
Applied rewrites97.1%
Applied rewrites96.2%
if 2.00000000000000005e144 < beta Initial program 71.5%
Taylor expanded in beta around -inf
mul-1-negN/A
lower-neg.f64N/A
sub-negN/A
mul-1-negN/A
distribute-neg-inN/A
+-commutativeN/A
distribute-neg-inN/A
metadata-evalN/A
unsub-negN/A
lower--.f6484.8
Applied rewrites84.8%
Applied rewrites84.8%
Final simplification94.5%
NOTE: alpha and beta should be sorted in increasing order before calling this function.
(FPCore (alpha beta)
:precision binary64
(let* ((t_0 (+ 2.0 (+ alpha beta))) (t_1 (+ 3.0 (+ alpha beta))))
(if (<= beta 1e+41)
(/ (+ 1.0 (fma beta alpha (+ alpha beta))) (* (* t_1 t_0) t_0))
(/ (/ (+ 1.0 alpha) (+ (+ 2.0 beta) alpha)) t_1))))assert(alpha < beta);
double code(double alpha, double beta) {
double t_0 = 2.0 + (alpha + beta);
double t_1 = 3.0 + (alpha + beta);
double tmp;
if (beta <= 1e+41) {
tmp = (1.0 + fma(beta, alpha, (alpha + beta))) / ((t_1 * t_0) * t_0);
} else {
tmp = ((1.0 + alpha) / ((2.0 + beta) + alpha)) / t_1;
}
return tmp;
}
alpha, beta = sort([alpha, beta]) function code(alpha, beta) t_0 = Float64(2.0 + Float64(alpha + beta)) t_1 = Float64(3.0 + Float64(alpha + beta)) tmp = 0.0 if (beta <= 1e+41) tmp = Float64(Float64(1.0 + fma(beta, alpha, Float64(alpha + beta))) / Float64(Float64(t_1 * t_0) * t_0)); else tmp = Float64(Float64(Float64(1.0 + alpha) / Float64(Float64(2.0 + beta) + alpha)) / t_1); end return tmp end
NOTE: alpha and beta should be sorted in increasing order before calling this function.
code[alpha_, beta_] := Block[{t$95$0 = N[(2.0 + N[(alpha + beta), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(3.0 + N[(alpha + beta), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[beta, 1e+41], N[(N[(1.0 + N[(beta * alpha + N[(alpha + beta), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(N[(t$95$1 * t$95$0), $MachinePrecision] * t$95$0), $MachinePrecision]), $MachinePrecision], N[(N[(N[(1.0 + alpha), $MachinePrecision] / N[(N[(2.0 + beta), $MachinePrecision] + alpha), $MachinePrecision]), $MachinePrecision] / t$95$1), $MachinePrecision]]]]
\begin{array}{l}
[alpha, beta] = \mathsf{sort}([alpha, beta])\\
\\
\begin{array}{l}
t_0 := 2 + \left(\alpha + \beta\right)\\
t_1 := 3 + \left(\alpha + \beta\right)\\
\mathbf{if}\;\beta \leq 10^{+41}:\\
\;\;\;\;\frac{1 + \mathsf{fma}\left(\beta, \alpha, \alpha + \beta\right)}{\left(t\_1 \cdot t\_0\right) \cdot t\_0}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{1 + \alpha}{\left(2 + \beta\right) + \alpha}}{t\_1}\\
\end{array}
\end{array}
if beta < 1.00000000000000001e41Initial program 99.3%
lift-/.f64N/A
lift-/.f64N/A
associate-/l/N/A
lift-/.f64N/A
associate-/l/N/A
lower-/.f64N/A
Applied rewrites94.8%
if 1.00000000000000001e41 < beta Initial program 77.7%
Taylor expanded in beta around -inf
mul-1-negN/A
lower-neg.f64N/A
sub-negN/A
mul-1-negN/A
distribute-neg-inN/A
+-commutativeN/A
distribute-neg-inN/A
metadata-evalN/A
unsub-negN/A
lower--.f6479.3
Applied rewrites79.3%
Applied rewrites79.3%
Final simplification90.5%
NOTE: alpha and beta should be sorted in increasing order before calling this function.
(FPCore (alpha beta)
:precision binary64
(let* ((t_0 (+ (+ 2.0 beta) alpha)))
(if (<= beta 3.7e+28)
(/ (/ (+ 1.0 beta) (fma (+ 5.0 beta) beta 6.0)) t_0)
(/ (/ (+ 1.0 alpha) t_0) (+ 3.0 (+ alpha beta))))))assert(alpha < beta);
double code(double alpha, double beta) {
double t_0 = (2.0 + beta) + alpha;
double tmp;
if (beta <= 3.7e+28) {
tmp = ((1.0 + beta) / fma((5.0 + beta), beta, 6.0)) / t_0;
} else {
tmp = ((1.0 + alpha) / t_0) / (3.0 + (alpha + beta));
}
return tmp;
}
alpha, beta = sort([alpha, beta]) function code(alpha, beta) t_0 = Float64(Float64(2.0 + beta) + alpha) tmp = 0.0 if (beta <= 3.7e+28) tmp = Float64(Float64(Float64(1.0 + beta) / fma(Float64(5.0 + beta), beta, 6.0)) / t_0); else tmp = Float64(Float64(Float64(1.0 + alpha) / t_0) / Float64(3.0 + Float64(alpha + beta))); end return tmp end
NOTE: alpha and beta should be sorted in increasing order before calling this function.
code[alpha_, beta_] := Block[{t$95$0 = N[(N[(2.0 + beta), $MachinePrecision] + alpha), $MachinePrecision]}, If[LessEqual[beta, 3.7e+28], N[(N[(N[(1.0 + beta), $MachinePrecision] / N[(N[(5.0 + beta), $MachinePrecision] * beta + 6.0), $MachinePrecision]), $MachinePrecision] / t$95$0), $MachinePrecision], N[(N[(N[(1.0 + alpha), $MachinePrecision] / t$95$0), $MachinePrecision] / N[(3.0 + N[(alpha + beta), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
[alpha, beta] = \mathsf{sort}([alpha, beta])\\
\\
\begin{array}{l}
t_0 := \left(2 + \beta\right) + \alpha\\
\mathbf{if}\;\beta \leq 3.7 \cdot 10^{+28}:\\
\;\;\;\;\frac{\frac{1 + \beta}{\mathsf{fma}\left(5 + \beta, \beta, 6\right)}}{t\_0}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{1 + \alpha}{t\_0}}{3 + \left(\alpha + \beta\right)}\\
\end{array}
\end{array}
if beta < 3.6999999999999999e28Initial program 99.3%
Taylor expanded in beta around -inf
mul-1-negN/A
lower-neg.f64N/A
sub-negN/A
mul-1-negN/A
distribute-neg-inN/A
+-commutativeN/A
distribute-neg-inN/A
metadata-evalN/A
unsub-negN/A
lower--.f6415.4
Applied rewrites15.4%
lift-/.f64N/A
lift-/.f64N/A
associate-/l/N/A
associate-/r*N/A
Applied rewrites15.4%
Taylor expanded in alpha around 0
lower-/.f64N/A
lower-+.f64N/A
lower-*.f64N/A
lower-+.f64N/A
lower-+.f6466.2
Applied rewrites66.2%
Taylor expanded in beta around 0
Applied rewrites66.2%
if 3.6999999999999999e28 < beta Initial program 78.3%
Taylor expanded in beta around -inf
mul-1-negN/A
lower-neg.f64N/A
sub-negN/A
mul-1-negN/A
distribute-neg-inN/A
+-commutativeN/A
distribute-neg-inN/A
metadata-evalN/A
unsub-negN/A
lower--.f6479.8
Applied rewrites79.8%
Applied rewrites79.8%
Final simplification70.1%
NOTE: alpha and beta should be sorted in increasing order before calling this function.
(FPCore (alpha beta)
:precision binary64
(let* ((t_0 (+ (+ 2.0 beta) alpha)))
(if (<= beta 1.5)
(/
(fma
(fma
(fma 0.03780864197530864 beta -0.05092592592592592)
beta
0.027777777777777776)
beta
0.16666666666666666)
t_0)
(/ (/ (+ 1.0 alpha) t_0) (+ 3.0 (+ alpha beta))))))assert(alpha < beta);
double code(double alpha, double beta) {
double t_0 = (2.0 + beta) + alpha;
double tmp;
if (beta <= 1.5) {
tmp = fma(fma(fma(0.03780864197530864, beta, -0.05092592592592592), beta, 0.027777777777777776), beta, 0.16666666666666666) / t_0;
} else {
tmp = ((1.0 + alpha) / t_0) / (3.0 + (alpha + beta));
}
return tmp;
}
alpha, beta = sort([alpha, beta]) function code(alpha, beta) t_0 = Float64(Float64(2.0 + beta) + alpha) tmp = 0.0 if (beta <= 1.5) tmp = Float64(fma(fma(fma(0.03780864197530864, beta, -0.05092592592592592), beta, 0.027777777777777776), beta, 0.16666666666666666) / t_0); else tmp = Float64(Float64(Float64(1.0 + alpha) / t_0) / Float64(3.0 + Float64(alpha + beta))); end return tmp end
NOTE: alpha and beta should be sorted in increasing order before calling this function.
code[alpha_, beta_] := Block[{t$95$0 = N[(N[(2.0 + beta), $MachinePrecision] + alpha), $MachinePrecision]}, If[LessEqual[beta, 1.5], N[(N[(N[(N[(0.03780864197530864 * beta + -0.05092592592592592), $MachinePrecision] * beta + 0.027777777777777776), $MachinePrecision] * beta + 0.16666666666666666), $MachinePrecision] / t$95$0), $MachinePrecision], N[(N[(N[(1.0 + alpha), $MachinePrecision] / t$95$0), $MachinePrecision] / N[(3.0 + N[(alpha + beta), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
[alpha, beta] = \mathsf{sort}([alpha, beta])\\
\\
\begin{array}{l}
t_0 := \left(2 + \beta\right) + \alpha\\
\mathbf{if}\;\beta \leq 1.5:\\
\;\;\;\;\frac{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(0.03780864197530864, \beta, -0.05092592592592592\right), \beta, 0.027777777777777776\right), \beta, 0.16666666666666666\right)}{t\_0}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{1 + \alpha}{t\_0}}{3 + \left(\alpha + \beta\right)}\\
\end{array}
\end{array}
if beta < 1.5Initial program 99.9%
Taylor expanded in beta around -inf
mul-1-negN/A
lower-neg.f64N/A
sub-negN/A
mul-1-negN/A
distribute-neg-inN/A
+-commutativeN/A
distribute-neg-inN/A
metadata-evalN/A
unsub-negN/A
lower--.f6414.5
Applied rewrites14.5%
lift-/.f64N/A
lift-/.f64N/A
associate-/l/N/A
associate-/r*N/A
Applied rewrites14.5%
Taylor expanded in alpha around 0
lower-/.f64N/A
lower-+.f64N/A
lower-*.f64N/A
lower-+.f64N/A
lower-+.f6466.1
Applied rewrites66.1%
Taylor expanded in beta around 0
Applied rewrites65.7%
if 1.5 < beta Initial program 78.4%
Taylor expanded in beta around -inf
mul-1-negN/A
lower-neg.f64N/A
sub-negN/A
mul-1-negN/A
distribute-neg-inN/A
+-commutativeN/A
distribute-neg-inN/A
metadata-evalN/A
unsub-negN/A
lower--.f6477.9
Applied rewrites77.9%
Applied rewrites77.9%
Final simplification69.4%
NOTE: alpha and beta should be sorted in increasing order before calling this function.
(FPCore (alpha beta)
:precision binary64
(let* ((t_0 (+ (+ 2.0 beta) alpha)))
(if (<= beta 1.5)
(/
(fma
(fma
(fma 0.03780864197530864 beta -0.05092592592592592)
beta
0.027777777777777776)
beta
0.16666666666666666)
t_0)
(/ (/ (+ 1.0 alpha) (+ 3.0 (+ alpha beta))) t_0))))assert(alpha < beta);
double code(double alpha, double beta) {
double t_0 = (2.0 + beta) + alpha;
double tmp;
if (beta <= 1.5) {
tmp = fma(fma(fma(0.03780864197530864, beta, -0.05092592592592592), beta, 0.027777777777777776), beta, 0.16666666666666666) / t_0;
} else {
tmp = ((1.0 + alpha) / (3.0 + (alpha + beta))) / t_0;
}
return tmp;
}
alpha, beta = sort([alpha, beta]) function code(alpha, beta) t_0 = Float64(Float64(2.0 + beta) + alpha) tmp = 0.0 if (beta <= 1.5) tmp = Float64(fma(fma(fma(0.03780864197530864, beta, -0.05092592592592592), beta, 0.027777777777777776), beta, 0.16666666666666666) / t_0); else tmp = Float64(Float64(Float64(1.0 + alpha) / Float64(3.0 + Float64(alpha + beta))) / t_0); end return tmp end
NOTE: alpha and beta should be sorted in increasing order before calling this function.
code[alpha_, beta_] := Block[{t$95$0 = N[(N[(2.0 + beta), $MachinePrecision] + alpha), $MachinePrecision]}, If[LessEqual[beta, 1.5], N[(N[(N[(N[(0.03780864197530864 * beta + -0.05092592592592592), $MachinePrecision] * beta + 0.027777777777777776), $MachinePrecision] * beta + 0.16666666666666666), $MachinePrecision] / t$95$0), $MachinePrecision], N[(N[(N[(1.0 + alpha), $MachinePrecision] / N[(3.0 + N[(alpha + beta), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / t$95$0), $MachinePrecision]]]
\begin{array}{l}
[alpha, beta] = \mathsf{sort}([alpha, beta])\\
\\
\begin{array}{l}
t_0 := \left(2 + \beta\right) + \alpha\\
\mathbf{if}\;\beta \leq 1.5:\\
\;\;\;\;\frac{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(0.03780864197530864, \beta, -0.05092592592592592\right), \beta, 0.027777777777777776\right), \beta, 0.16666666666666666\right)}{t\_0}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{1 + \alpha}{3 + \left(\alpha + \beta\right)}}{t\_0}\\
\end{array}
\end{array}
if beta < 1.5Initial program 99.9%
Taylor expanded in beta around -inf
mul-1-negN/A
lower-neg.f64N/A
sub-negN/A
mul-1-negN/A
distribute-neg-inN/A
+-commutativeN/A
distribute-neg-inN/A
metadata-evalN/A
unsub-negN/A
lower--.f6414.5
Applied rewrites14.5%
lift-/.f64N/A
lift-/.f64N/A
associate-/l/N/A
associate-/r*N/A
Applied rewrites14.5%
Taylor expanded in alpha around 0
lower-/.f64N/A
lower-+.f64N/A
lower-*.f64N/A
lower-+.f64N/A
lower-+.f6466.1
Applied rewrites66.1%
Taylor expanded in beta around 0
Applied rewrites65.7%
if 1.5 < beta Initial program 78.4%
Taylor expanded in beta around -inf
mul-1-negN/A
lower-neg.f64N/A
sub-negN/A
mul-1-negN/A
distribute-neg-inN/A
+-commutativeN/A
distribute-neg-inN/A
metadata-evalN/A
unsub-negN/A
lower--.f6477.9
Applied rewrites77.9%
lift-/.f64N/A
lift-/.f64N/A
associate-/l/N/A
associate-/r*N/A
Applied rewrites77.9%
Final simplification69.4%
NOTE: alpha and beta should be sorted in increasing order before calling this function.
(FPCore (alpha beta)
:precision binary64
(let* ((t_0 (+ (+ 2.0 beta) alpha)))
(if (<= beta 2.4)
(/
(fma
(fma
(fma 0.03780864197530864 beta -0.05092592592592592)
beta
0.027777777777777776)
beta
0.16666666666666666)
t_0)
(/ (/ (+ 1.0 alpha) beta) t_0))))assert(alpha < beta);
double code(double alpha, double beta) {
double t_0 = (2.0 + beta) + alpha;
double tmp;
if (beta <= 2.4) {
tmp = fma(fma(fma(0.03780864197530864, beta, -0.05092592592592592), beta, 0.027777777777777776), beta, 0.16666666666666666) / t_0;
} else {
tmp = ((1.0 + alpha) / beta) / t_0;
}
return tmp;
}
alpha, beta = sort([alpha, beta]) function code(alpha, beta) t_0 = Float64(Float64(2.0 + beta) + alpha) tmp = 0.0 if (beta <= 2.4) tmp = Float64(fma(fma(fma(0.03780864197530864, beta, -0.05092592592592592), beta, 0.027777777777777776), beta, 0.16666666666666666) / t_0); else tmp = Float64(Float64(Float64(1.0 + alpha) / beta) / t_0); end return tmp end
NOTE: alpha and beta should be sorted in increasing order before calling this function.
code[alpha_, beta_] := Block[{t$95$0 = N[(N[(2.0 + beta), $MachinePrecision] + alpha), $MachinePrecision]}, If[LessEqual[beta, 2.4], N[(N[(N[(N[(0.03780864197530864 * beta + -0.05092592592592592), $MachinePrecision] * beta + 0.027777777777777776), $MachinePrecision] * beta + 0.16666666666666666), $MachinePrecision] / t$95$0), $MachinePrecision], N[(N[(N[(1.0 + alpha), $MachinePrecision] / beta), $MachinePrecision] / t$95$0), $MachinePrecision]]]
\begin{array}{l}
[alpha, beta] = \mathsf{sort}([alpha, beta])\\
\\
\begin{array}{l}
t_0 := \left(2 + \beta\right) + \alpha\\
\mathbf{if}\;\beta \leq 2.4:\\
\;\;\;\;\frac{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(0.03780864197530864, \beta, -0.05092592592592592\right), \beta, 0.027777777777777776\right), \beta, 0.16666666666666666\right)}{t\_0}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{1 + \alpha}{\beta}}{t\_0}\\
\end{array}
\end{array}
if beta < 2.39999999999999991Initial program 99.9%
Taylor expanded in beta around -inf
mul-1-negN/A
lower-neg.f64N/A
sub-negN/A
mul-1-negN/A
distribute-neg-inN/A
+-commutativeN/A
distribute-neg-inN/A
metadata-evalN/A
unsub-negN/A
lower--.f6414.5
Applied rewrites14.5%
lift-/.f64N/A
lift-/.f64N/A
associate-/l/N/A
associate-/r*N/A
Applied rewrites14.5%
Taylor expanded in alpha around 0
lower-/.f64N/A
lower-+.f64N/A
lower-*.f64N/A
lower-+.f64N/A
lower-+.f6466.1
Applied rewrites66.1%
Taylor expanded in beta around 0
Applied rewrites65.7%
if 2.39999999999999991 < beta Initial program 78.4%
Taylor expanded in beta around -inf
mul-1-negN/A
lower-neg.f64N/A
sub-negN/A
mul-1-negN/A
distribute-neg-inN/A
+-commutativeN/A
distribute-neg-inN/A
metadata-evalN/A
unsub-negN/A
lower--.f6477.9
Applied rewrites77.9%
lift-/.f64N/A
lift-/.f64N/A
associate-/l/N/A
associate-/r*N/A
Applied rewrites77.9%
Taylor expanded in beta around inf
lower-/.f64N/A
lower-+.f6477.1
Applied rewrites77.1%
NOTE: alpha and beta should be sorted in increasing order before calling this function.
(FPCore (alpha beta)
:precision binary64
(let* ((t_0 (+ (+ 2.0 beta) alpha)))
(if (<= beta 1.85)
(/
(fma
(fma -0.05092592592592592 beta 0.027777777777777776)
beta
0.16666666666666666)
t_0)
(/ (/ (+ 1.0 alpha) beta) t_0))))assert(alpha < beta);
double code(double alpha, double beta) {
double t_0 = (2.0 + beta) + alpha;
double tmp;
if (beta <= 1.85) {
tmp = fma(fma(-0.05092592592592592, beta, 0.027777777777777776), beta, 0.16666666666666666) / t_0;
} else {
tmp = ((1.0 + alpha) / beta) / t_0;
}
return tmp;
}
alpha, beta = sort([alpha, beta]) function code(alpha, beta) t_0 = Float64(Float64(2.0 + beta) + alpha) tmp = 0.0 if (beta <= 1.85) tmp = Float64(fma(fma(-0.05092592592592592, beta, 0.027777777777777776), beta, 0.16666666666666666) / t_0); else tmp = Float64(Float64(Float64(1.0 + alpha) / beta) / t_0); end return tmp end
NOTE: alpha and beta should be sorted in increasing order before calling this function.
code[alpha_, beta_] := Block[{t$95$0 = N[(N[(2.0 + beta), $MachinePrecision] + alpha), $MachinePrecision]}, If[LessEqual[beta, 1.85], N[(N[(N[(-0.05092592592592592 * beta + 0.027777777777777776), $MachinePrecision] * beta + 0.16666666666666666), $MachinePrecision] / t$95$0), $MachinePrecision], N[(N[(N[(1.0 + alpha), $MachinePrecision] / beta), $MachinePrecision] / t$95$0), $MachinePrecision]]]
\begin{array}{l}
[alpha, beta] = \mathsf{sort}([alpha, beta])\\
\\
\begin{array}{l}
t_0 := \left(2 + \beta\right) + \alpha\\
\mathbf{if}\;\beta \leq 1.85:\\
\;\;\;\;\frac{\mathsf{fma}\left(\mathsf{fma}\left(-0.05092592592592592, \beta, 0.027777777777777776\right), \beta, 0.16666666666666666\right)}{t\_0}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{1 + \alpha}{\beta}}{t\_0}\\
\end{array}
\end{array}
if beta < 1.8500000000000001Initial program 99.9%
Taylor expanded in beta around -inf
mul-1-negN/A
lower-neg.f64N/A
sub-negN/A
mul-1-negN/A
distribute-neg-inN/A
+-commutativeN/A
distribute-neg-inN/A
metadata-evalN/A
unsub-negN/A
lower--.f6414.5
Applied rewrites14.5%
lift-/.f64N/A
lift-/.f64N/A
associate-/l/N/A
associate-/r*N/A
Applied rewrites14.5%
Taylor expanded in alpha around 0
lower-/.f64N/A
lower-+.f64N/A
lower-*.f64N/A
lower-+.f64N/A
lower-+.f6466.1
Applied rewrites66.1%
Taylor expanded in beta around 0
Applied rewrites65.3%
if 1.8500000000000001 < beta Initial program 78.4%
Taylor expanded in beta around -inf
mul-1-negN/A
lower-neg.f64N/A
sub-negN/A
mul-1-negN/A
distribute-neg-inN/A
+-commutativeN/A
distribute-neg-inN/A
metadata-evalN/A
unsub-negN/A
lower--.f6477.9
Applied rewrites77.9%
lift-/.f64N/A
lift-/.f64N/A
associate-/l/N/A
associate-/r*N/A
Applied rewrites77.9%
Taylor expanded in beta around inf
lower-/.f64N/A
lower-+.f6477.1
Applied rewrites77.1%
NOTE: alpha and beta should be sorted in increasing order before calling this function.
(FPCore (alpha beta)
:precision binary64
(if (<= beta 2.0)
(/
(fma
(fma -0.05092592592592592 beta 0.027777777777777776)
beta
0.16666666666666666)
(+ (+ 2.0 beta) alpha))
(/ (/ (+ 1.0 alpha) beta) beta)))assert(alpha < beta);
double code(double alpha, double beta) {
double tmp;
if (beta <= 2.0) {
tmp = fma(fma(-0.05092592592592592, beta, 0.027777777777777776), beta, 0.16666666666666666) / ((2.0 + beta) + alpha);
} else {
tmp = ((1.0 + alpha) / beta) / beta;
}
return tmp;
}
alpha, beta = sort([alpha, beta]) function code(alpha, beta) tmp = 0.0 if (beta <= 2.0) tmp = Float64(fma(fma(-0.05092592592592592, beta, 0.027777777777777776), beta, 0.16666666666666666) / Float64(Float64(2.0 + beta) + alpha)); else tmp = Float64(Float64(Float64(1.0 + alpha) / beta) / beta); end return tmp end
NOTE: alpha and beta should be sorted in increasing order before calling this function. code[alpha_, beta_] := If[LessEqual[beta, 2.0], N[(N[(N[(-0.05092592592592592 * beta + 0.027777777777777776), $MachinePrecision] * beta + 0.16666666666666666), $MachinePrecision] / N[(N[(2.0 + beta), $MachinePrecision] + alpha), $MachinePrecision]), $MachinePrecision], N[(N[(N[(1.0 + alpha), $MachinePrecision] / beta), $MachinePrecision] / beta), $MachinePrecision]]
\begin{array}{l}
[alpha, beta] = \mathsf{sort}([alpha, beta])\\
\\
\begin{array}{l}
\mathbf{if}\;\beta \leq 2:\\
\;\;\;\;\frac{\mathsf{fma}\left(\mathsf{fma}\left(-0.05092592592592592, \beta, 0.027777777777777776\right), \beta, 0.16666666666666666\right)}{\left(2 + \beta\right) + \alpha}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{1 + \alpha}{\beta}}{\beta}\\
\end{array}
\end{array}
if beta < 2Initial program 99.9%
Taylor expanded in beta around -inf
mul-1-negN/A
lower-neg.f64N/A
sub-negN/A
mul-1-negN/A
distribute-neg-inN/A
+-commutativeN/A
distribute-neg-inN/A
metadata-evalN/A
unsub-negN/A
lower--.f6414.5
Applied rewrites14.5%
lift-/.f64N/A
lift-/.f64N/A
associate-/l/N/A
associate-/r*N/A
Applied rewrites14.5%
Taylor expanded in alpha around 0
lower-/.f64N/A
lower-+.f64N/A
lower-*.f64N/A
lower-+.f64N/A
lower-+.f6466.1
Applied rewrites66.1%
Taylor expanded in beta around 0
Applied rewrites65.3%
if 2 < beta Initial program 78.4%
Taylor expanded in beta around inf
lower-/.f64N/A
lower-+.f64N/A
unpow2N/A
lower-*.f6471.9
Applied rewrites71.9%
Applied rewrites76.8%
Final simplification68.8%
NOTE: alpha and beta should be sorted in increasing order before calling this function.
(FPCore (alpha beta)
:precision binary64
(if (<= beta 5.5)
(/
(fma 0.027777777777777776 beta 0.16666666666666666)
(+ (+ 2.0 beta) alpha))
(if (<= beta 5.2e+155)
(/ (+ 1.0 alpha) (* beta beta))
(/ (/ alpha beta) beta))))assert(alpha < beta);
double code(double alpha, double beta) {
double tmp;
if (beta <= 5.5) {
tmp = fma(0.027777777777777776, beta, 0.16666666666666666) / ((2.0 + beta) + alpha);
} else if (beta <= 5.2e+155) {
tmp = (1.0 + alpha) / (beta * beta);
} else {
tmp = (alpha / beta) / beta;
}
return tmp;
}
alpha, beta = sort([alpha, beta]) function code(alpha, beta) tmp = 0.0 if (beta <= 5.5) tmp = Float64(fma(0.027777777777777776, beta, 0.16666666666666666) / Float64(Float64(2.0 + beta) + alpha)); elseif (beta <= 5.2e+155) tmp = Float64(Float64(1.0 + alpha) / Float64(beta * beta)); else tmp = Float64(Float64(alpha / beta) / beta); end return tmp end
NOTE: alpha and beta should be sorted in increasing order before calling this function. code[alpha_, beta_] := If[LessEqual[beta, 5.5], N[(N[(0.027777777777777776 * beta + 0.16666666666666666), $MachinePrecision] / N[(N[(2.0 + beta), $MachinePrecision] + alpha), $MachinePrecision]), $MachinePrecision], If[LessEqual[beta, 5.2e+155], N[(N[(1.0 + alpha), $MachinePrecision] / N[(beta * beta), $MachinePrecision]), $MachinePrecision], N[(N[(alpha / beta), $MachinePrecision] / beta), $MachinePrecision]]]
\begin{array}{l}
[alpha, beta] = \mathsf{sort}([alpha, beta])\\
\\
\begin{array}{l}
\mathbf{if}\;\beta \leq 5.5:\\
\;\;\;\;\frac{\mathsf{fma}\left(0.027777777777777776, \beta, 0.16666666666666666\right)}{\left(2 + \beta\right) + \alpha}\\
\mathbf{elif}\;\beta \leq 5.2 \cdot 10^{+155}:\\
\;\;\;\;\frac{1 + \alpha}{\beta \cdot \beta}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{\alpha}{\beta}}{\beta}\\
\end{array}
\end{array}
if beta < 5.5Initial program 99.9%
Taylor expanded in beta around -inf
mul-1-negN/A
lower-neg.f64N/A
sub-negN/A
mul-1-negN/A
distribute-neg-inN/A
+-commutativeN/A
distribute-neg-inN/A
metadata-evalN/A
unsub-negN/A
lower--.f6414.5
Applied rewrites14.5%
lift-/.f64N/A
lift-/.f64N/A
associate-/l/N/A
associate-/r*N/A
Applied rewrites14.5%
Taylor expanded in alpha around 0
lower-/.f64N/A
lower-+.f64N/A
lower-*.f64N/A
lower-+.f64N/A
lower-+.f6466.1
Applied rewrites66.1%
Taylor expanded in beta around 0
Applied rewrites64.9%
if 5.5 < beta < 5.2000000000000004e155Initial program 84.0%
Taylor expanded in beta around inf
lower-/.f64N/A
lower-+.f64N/A
unpow2N/A
lower-*.f6467.2
Applied rewrites67.2%
if 5.2000000000000004e155 < beta Initial program 71.1%
Taylor expanded in beta around inf
lower-/.f64N/A
lower-+.f64N/A
unpow2N/A
lower-*.f6477.9
Applied rewrites77.9%
Taylor expanded in alpha around inf
Applied rewrites77.9%
Applied rewrites84.5%
NOTE: alpha and beta should be sorted in increasing order before calling this function.
(FPCore (alpha beta)
:precision binary64
(if (<= beta 8.5)
(/ 0.16666666666666666 (+ (+ 2.0 beta) alpha))
(if (<= beta 5.2e+155)
(/ (+ 1.0 alpha) (* beta beta))
(/ (/ alpha beta) beta))))assert(alpha < beta);
double code(double alpha, double beta) {
double tmp;
if (beta <= 8.5) {
tmp = 0.16666666666666666 / ((2.0 + beta) + alpha);
} else if (beta <= 5.2e+155) {
tmp = (1.0 + alpha) / (beta * beta);
} else {
tmp = (alpha / beta) / beta;
}
return tmp;
}
NOTE: alpha and beta should be sorted in increasing order before calling this function.
real(8) function code(alpha, beta)
real(8), intent (in) :: alpha
real(8), intent (in) :: beta
real(8) :: tmp
if (beta <= 8.5d0) then
tmp = 0.16666666666666666d0 / ((2.0d0 + beta) + alpha)
else if (beta <= 5.2d+155) then
tmp = (1.0d0 + alpha) / (beta * beta)
else
tmp = (alpha / beta) / beta
end if
code = tmp
end function
assert alpha < beta;
public static double code(double alpha, double beta) {
double tmp;
if (beta <= 8.5) {
tmp = 0.16666666666666666 / ((2.0 + beta) + alpha);
} else if (beta <= 5.2e+155) {
tmp = (1.0 + alpha) / (beta * beta);
} else {
tmp = (alpha / beta) / beta;
}
return tmp;
}
[alpha, beta] = sort([alpha, beta]) def code(alpha, beta): tmp = 0 if beta <= 8.5: tmp = 0.16666666666666666 / ((2.0 + beta) + alpha) elif beta <= 5.2e+155: tmp = (1.0 + alpha) / (beta * beta) else: tmp = (alpha / beta) / beta return tmp
alpha, beta = sort([alpha, beta]) function code(alpha, beta) tmp = 0.0 if (beta <= 8.5) tmp = Float64(0.16666666666666666 / Float64(Float64(2.0 + beta) + alpha)); elseif (beta <= 5.2e+155) tmp = Float64(Float64(1.0 + alpha) / Float64(beta * beta)); else tmp = Float64(Float64(alpha / beta) / beta); end return tmp end
alpha, beta = num2cell(sort([alpha, beta])){:}
function tmp_2 = code(alpha, beta)
tmp = 0.0;
if (beta <= 8.5)
tmp = 0.16666666666666666 / ((2.0 + beta) + alpha);
elseif (beta <= 5.2e+155)
tmp = (1.0 + alpha) / (beta * beta);
else
tmp = (alpha / beta) / beta;
end
tmp_2 = tmp;
end
NOTE: alpha and beta should be sorted in increasing order before calling this function. code[alpha_, beta_] := If[LessEqual[beta, 8.5], N[(0.16666666666666666 / N[(N[(2.0 + beta), $MachinePrecision] + alpha), $MachinePrecision]), $MachinePrecision], If[LessEqual[beta, 5.2e+155], N[(N[(1.0 + alpha), $MachinePrecision] / N[(beta * beta), $MachinePrecision]), $MachinePrecision], N[(N[(alpha / beta), $MachinePrecision] / beta), $MachinePrecision]]]
\begin{array}{l}
[alpha, beta] = \mathsf{sort}([alpha, beta])\\
\\
\begin{array}{l}
\mathbf{if}\;\beta \leq 8.5:\\
\;\;\;\;\frac{0.16666666666666666}{\left(2 + \beta\right) + \alpha}\\
\mathbf{elif}\;\beta \leq 5.2 \cdot 10^{+155}:\\
\;\;\;\;\frac{1 + \alpha}{\beta \cdot \beta}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{\alpha}{\beta}}{\beta}\\
\end{array}
\end{array}
if beta < 8.5Initial program 99.9%
Taylor expanded in beta around -inf
mul-1-negN/A
lower-neg.f64N/A
sub-negN/A
mul-1-negN/A
distribute-neg-inN/A
+-commutativeN/A
distribute-neg-inN/A
metadata-evalN/A
unsub-negN/A
lower--.f6414.5
Applied rewrites14.5%
lift-/.f64N/A
lift-/.f64N/A
associate-/l/N/A
associate-/r*N/A
Applied rewrites14.5%
Taylor expanded in alpha around 0
lower-/.f64N/A
lower-+.f64N/A
lower-*.f64N/A
lower-+.f64N/A
lower-+.f6466.1
Applied rewrites66.1%
Taylor expanded in beta around 0
Applied rewrites64.1%
if 8.5 < beta < 5.2000000000000004e155Initial program 84.0%
Taylor expanded in beta around inf
lower-/.f64N/A
lower-+.f64N/A
unpow2N/A
lower-*.f6467.2
Applied rewrites67.2%
if 5.2000000000000004e155 < beta Initial program 71.1%
Taylor expanded in beta around inf
lower-/.f64N/A
lower-+.f64N/A
unpow2N/A
lower-*.f6477.9
Applied rewrites77.9%
Taylor expanded in alpha around inf
Applied rewrites77.9%
Applied rewrites84.5%
NOTE: alpha and beta should be sorted in increasing order before calling this function.
(FPCore (alpha beta)
:precision binary64
(if (<= beta 5.5)
(/
(fma 0.027777777777777776 beta 0.16666666666666666)
(+ (+ 2.0 beta) alpha))
(/ (/ (+ 1.0 alpha) beta) beta)))assert(alpha < beta);
double code(double alpha, double beta) {
double tmp;
if (beta <= 5.5) {
tmp = fma(0.027777777777777776, beta, 0.16666666666666666) / ((2.0 + beta) + alpha);
} else {
tmp = ((1.0 + alpha) / beta) / beta;
}
return tmp;
}
alpha, beta = sort([alpha, beta]) function code(alpha, beta) tmp = 0.0 if (beta <= 5.5) tmp = Float64(fma(0.027777777777777776, beta, 0.16666666666666666) / Float64(Float64(2.0 + beta) + alpha)); else tmp = Float64(Float64(Float64(1.0 + alpha) / beta) / beta); end return tmp end
NOTE: alpha and beta should be sorted in increasing order before calling this function. code[alpha_, beta_] := If[LessEqual[beta, 5.5], N[(N[(0.027777777777777776 * beta + 0.16666666666666666), $MachinePrecision] / N[(N[(2.0 + beta), $MachinePrecision] + alpha), $MachinePrecision]), $MachinePrecision], N[(N[(N[(1.0 + alpha), $MachinePrecision] / beta), $MachinePrecision] / beta), $MachinePrecision]]
\begin{array}{l}
[alpha, beta] = \mathsf{sort}([alpha, beta])\\
\\
\begin{array}{l}
\mathbf{if}\;\beta \leq 5.5:\\
\;\;\;\;\frac{\mathsf{fma}\left(0.027777777777777776, \beta, 0.16666666666666666\right)}{\left(2 + \beta\right) + \alpha}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{1 + \alpha}{\beta}}{\beta}\\
\end{array}
\end{array}
if beta < 5.5Initial program 99.9%
Taylor expanded in beta around -inf
mul-1-negN/A
lower-neg.f64N/A
sub-negN/A
mul-1-negN/A
distribute-neg-inN/A
+-commutativeN/A
distribute-neg-inN/A
metadata-evalN/A
unsub-negN/A
lower--.f6414.5
Applied rewrites14.5%
lift-/.f64N/A
lift-/.f64N/A
associate-/l/N/A
associate-/r*N/A
Applied rewrites14.5%
Taylor expanded in alpha around 0
lower-/.f64N/A
lower-+.f64N/A
lower-*.f64N/A
lower-+.f64N/A
lower-+.f6466.1
Applied rewrites66.1%
Taylor expanded in beta around 0
Applied rewrites64.9%
if 5.5 < beta Initial program 78.4%
Taylor expanded in beta around inf
lower-/.f64N/A
lower-+.f64N/A
unpow2N/A
lower-*.f6471.9
Applied rewrites71.9%
Applied rewrites76.8%
Final simplification68.5%
NOTE: alpha and beta should be sorted in increasing order before calling this function. (FPCore (alpha beta) :precision binary64 (if (<= beta 8.5) (/ 0.16666666666666666 (+ (+ 2.0 beta) alpha)) (/ (+ 1.0 alpha) (* beta beta))))
assert(alpha < beta);
double code(double alpha, double beta) {
double tmp;
if (beta <= 8.5) {
tmp = 0.16666666666666666 / ((2.0 + beta) + alpha);
} else {
tmp = (1.0 + alpha) / (beta * beta);
}
return tmp;
}
NOTE: alpha and beta should be sorted in increasing order before calling this function.
real(8) function code(alpha, beta)
real(8), intent (in) :: alpha
real(8), intent (in) :: beta
real(8) :: tmp
if (beta <= 8.5d0) then
tmp = 0.16666666666666666d0 / ((2.0d0 + beta) + alpha)
else
tmp = (1.0d0 + alpha) / (beta * beta)
end if
code = tmp
end function
assert alpha < beta;
public static double code(double alpha, double beta) {
double tmp;
if (beta <= 8.5) {
tmp = 0.16666666666666666 / ((2.0 + beta) + alpha);
} else {
tmp = (1.0 + alpha) / (beta * beta);
}
return tmp;
}
[alpha, beta] = sort([alpha, beta]) def code(alpha, beta): tmp = 0 if beta <= 8.5: tmp = 0.16666666666666666 / ((2.0 + beta) + alpha) else: tmp = (1.0 + alpha) / (beta * beta) return tmp
alpha, beta = sort([alpha, beta]) function code(alpha, beta) tmp = 0.0 if (beta <= 8.5) tmp = Float64(0.16666666666666666 / Float64(Float64(2.0 + beta) + alpha)); else tmp = Float64(Float64(1.0 + alpha) / Float64(beta * beta)); end return tmp end
alpha, beta = num2cell(sort([alpha, beta])){:}
function tmp_2 = code(alpha, beta)
tmp = 0.0;
if (beta <= 8.5)
tmp = 0.16666666666666666 / ((2.0 + beta) + alpha);
else
tmp = (1.0 + alpha) / (beta * beta);
end
tmp_2 = tmp;
end
NOTE: alpha and beta should be sorted in increasing order before calling this function. code[alpha_, beta_] := If[LessEqual[beta, 8.5], N[(0.16666666666666666 / N[(N[(2.0 + beta), $MachinePrecision] + alpha), $MachinePrecision]), $MachinePrecision], N[(N[(1.0 + alpha), $MachinePrecision] / N[(beta * beta), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
[alpha, beta] = \mathsf{sort}([alpha, beta])\\
\\
\begin{array}{l}
\mathbf{if}\;\beta \leq 8.5:\\
\;\;\;\;\frac{0.16666666666666666}{\left(2 + \beta\right) + \alpha}\\
\mathbf{else}:\\
\;\;\;\;\frac{1 + \alpha}{\beta \cdot \beta}\\
\end{array}
\end{array}
if beta < 8.5Initial program 99.9%
Taylor expanded in beta around -inf
mul-1-negN/A
lower-neg.f64N/A
sub-negN/A
mul-1-negN/A
distribute-neg-inN/A
+-commutativeN/A
distribute-neg-inN/A
metadata-evalN/A
unsub-negN/A
lower--.f6414.5
Applied rewrites14.5%
lift-/.f64N/A
lift-/.f64N/A
associate-/l/N/A
associate-/r*N/A
Applied rewrites14.5%
Taylor expanded in alpha around 0
lower-/.f64N/A
lower-+.f64N/A
lower-*.f64N/A
lower-+.f64N/A
lower-+.f6466.1
Applied rewrites66.1%
Taylor expanded in beta around 0
Applied rewrites64.1%
if 8.5 < beta Initial program 78.4%
Taylor expanded in beta around inf
lower-/.f64N/A
lower-+.f64N/A
unpow2N/A
lower-*.f6471.9
Applied rewrites71.9%
NOTE: alpha and beta should be sorted in increasing order before calling this function. (FPCore (alpha beta) :precision binary64 (if (<= beta 8.0) (/ 0.16666666666666666 (+ (+ 2.0 beta) alpha)) (/ 1.0 (* beta beta))))
assert(alpha < beta);
double code(double alpha, double beta) {
double tmp;
if (beta <= 8.0) {
tmp = 0.16666666666666666 / ((2.0 + beta) + alpha);
} else {
tmp = 1.0 / (beta * beta);
}
return tmp;
}
NOTE: alpha and beta should be sorted in increasing order before calling this function.
real(8) function code(alpha, beta)
real(8), intent (in) :: alpha
real(8), intent (in) :: beta
real(8) :: tmp
if (beta <= 8.0d0) then
tmp = 0.16666666666666666d0 / ((2.0d0 + beta) + alpha)
else
tmp = 1.0d0 / (beta * beta)
end if
code = tmp
end function
assert alpha < beta;
public static double code(double alpha, double beta) {
double tmp;
if (beta <= 8.0) {
tmp = 0.16666666666666666 / ((2.0 + beta) + alpha);
} else {
tmp = 1.0 / (beta * beta);
}
return tmp;
}
[alpha, beta] = sort([alpha, beta]) def code(alpha, beta): tmp = 0 if beta <= 8.0: tmp = 0.16666666666666666 / ((2.0 + beta) + alpha) else: tmp = 1.0 / (beta * beta) return tmp
alpha, beta = sort([alpha, beta]) function code(alpha, beta) tmp = 0.0 if (beta <= 8.0) tmp = Float64(0.16666666666666666 / Float64(Float64(2.0 + beta) + alpha)); else tmp = Float64(1.0 / Float64(beta * beta)); end return tmp end
alpha, beta = num2cell(sort([alpha, beta])){:}
function tmp_2 = code(alpha, beta)
tmp = 0.0;
if (beta <= 8.0)
tmp = 0.16666666666666666 / ((2.0 + beta) + alpha);
else
tmp = 1.0 / (beta * beta);
end
tmp_2 = tmp;
end
NOTE: alpha and beta should be sorted in increasing order before calling this function. code[alpha_, beta_] := If[LessEqual[beta, 8.0], N[(0.16666666666666666 / N[(N[(2.0 + beta), $MachinePrecision] + alpha), $MachinePrecision]), $MachinePrecision], N[(1.0 / N[(beta * beta), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
[alpha, beta] = \mathsf{sort}([alpha, beta])\\
\\
\begin{array}{l}
\mathbf{if}\;\beta \leq 8:\\
\;\;\;\;\frac{0.16666666666666666}{\left(2 + \beta\right) + \alpha}\\
\mathbf{else}:\\
\;\;\;\;\frac{1}{\beta \cdot \beta}\\
\end{array}
\end{array}
if beta < 8Initial program 99.9%
Taylor expanded in beta around -inf
mul-1-negN/A
lower-neg.f64N/A
sub-negN/A
mul-1-negN/A
distribute-neg-inN/A
+-commutativeN/A
distribute-neg-inN/A
metadata-evalN/A
unsub-negN/A
lower--.f6414.5
Applied rewrites14.5%
lift-/.f64N/A
lift-/.f64N/A
associate-/l/N/A
associate-/r*N/A
Applied rewrites14.5%
Taylor expanded in alpha around 0
lower-/.f64N/A
lower-+.f64N/A
lower-*.f64N/A
lower-+.f64N/A
lower-+.f6466.1
Applied rewrites66.1%
Taylor expanded in beta around 0
Applied rewrites64.1%
if 8 < beta Initial program 78.4%
Taylor expanded in beta around inf
lower-/.f64N/A
lower-+.f64N/A
unpow2N/A
lower-*.f6471.9
Applied rewrites71.9%
Taylor expanded in alpha around 0
Applied rewrites68.4%
NOTE: alpha and beta should be sorted in increasing order before calling this function. (FPCore (alpha beta) :precision binary64 (/ 1.0 (* beta beta)))
assert(alpha < beta);
double code(double alpha, double beta) {
return 1.0 / (beta * beta);
}
NOTE: alpha and beta should be sorted in increasing order before calling this function.
real(8) function code(alpha, beta)
real(8), intent (in) :: alpha
real(8), intent (in) :: beta
code = 1.0d0 / (beta * beta)
end function
assert alpha < beta;
public static double code(double alpha, double beta) {
return 1.0 / (beta * beta);
}
[alpha, beta] = sort([alpha, beta]) def code(alpha, beta): return 1.0 / (beta * beta)
alpha, beta = sort([alpha, beta]) function code(alpha, beta) return Float64(1.0 / Float64(beta * beta)) end
alpha, beta = num2cell(sort([alpha, beta])){:}
function tmp = code(alpha, beta)
tmp = 1.0 / (beta * beta);
end
NOTE: alpha and beta should be sorted in increasing order before calling this function. code[alpha_, beta_] := N[(1.0 / N[(beta * beta), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
[alpha, beta] = \mathsf{sort}([alpha, beta])\\
\\
\frac{1}{\beta \cdot \beta}
\end{array}
Initial program 93.3%
Taylor expanded in beta around inf
lower-/.f64N/A
lower-+.f64N/A
unpow2N/A
lower-*.f6424.3
Applied rewrites24.3%
Taylor expanded in alpha around 0
Applied rewrites23.3%
NOTE: alpha and beta should be sorted in increasing order before calling this function. (FPCore (alpha beta) :precision binary64 (/ alpha (* beta beta)))
assert(alpha < beta);
double code(double alpha, double beta) {
return alpha / (beta * beta);
}
NOTE: alpha and beta should be sorted in increasing order before calling this function.
real(8) function code(alpha, beta)
real(8), intent (in) :: alpha
real(8), intent (in) :: beta
code = alpha / (beta * beta)
end function
assert alpha < beta;
public static double code(double alpha, double beta) {
return alpha / (beta * beta);
}
[alpha, beta] = sort([alpha, beta]) def code(alpha, beta): return alpha / (beta * beta)
alpha, beta = sort([alpha, beta]) function code(alpha, beta) return Float64(alpha / Float64(beta * beta)) end
alpha, beta = num2cell(sort([alpha, beta])){:}
function tmp = code(alpha, beta)
tmp = alpha / (beta * beta);
end
NOTE: alpha and beta should be sorted in increasing order before calling this function. code[alpha_, beta_] := N[(alpha / N[(beta * beta), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
[alpha, beta] = \mathsf{sort}([alpha, beta])\\
\\
\frac{\alpha}{\beta \cdot \beta}
\end{array}
Initial program 93.3%
Taylor expanded in beta around inf
lower-/.f64N/A
lower-+.f64N/A
unpow2N/A
lower-*.f6424.3
Applied rewrites24.3%
Taylor expanded in alpha around inf
Applied rewrites13.7%
herbie shell --seed 2024296
(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)))