
(FPCore (alpha beta i) :precision binary64 (let* ((t_0 (+ (+ alpha beta) (* 2.0 i)))) (/ (+ (/ (/ (* (+ alpha beta) (- beta alpha)) t_0) (+ t_0 2.0)) 1.0) 2.0)))
double code(double alpha, double beta, double i) {
double t_0 = (alpha + beta) + (2.0 * i);
return (((((alpha + beta) * (beta - alpha)) / t_0) / (t_0 + 2.0)) + 1.0) / 2.0;
}
real(8) function code(alpha, beta, i)
real(8), intent (in) :: alpha
real(8), intent (in) :: beta
real(8), intent (in) :: i
real(8) :: t_0
t_0 = (alpha + beta) + (2.0d0 * i)
code = (((((alpha + beta) * (beta - alpha)) / t_0) / (t_0 + 2.0d0)) + 1.0d0) / 2.0d0
end function
public static double code(double alpha, double beta, double i) {
double t_0 = (alpha + beta) + (2.0 * i);
return (((((alpha + beta) * (beta - alpha)) / t_0) / (t_0 + 2.0)) + 1.0) / 2.0;
}
def code(alpha, beta, i): t_0 = (alpha + beta) + (2.0 * i) return (((((alpha + beta) * (beta - alpha)) / t_0) / (t_0 + 2.0)) + 1.0) / 2.0
function code(alpha, beta, i) t_0 = Float64(Float64(alpha + beta) + Float64(2.0 * i)) return Float64(Float64(Float64(Float64(Float64(Float64(alpha + beta) * Float64(beta - alpha)) / t_0) / Float64(t_0 + 2.0)) + 1.0) / 2.0) end
function tmp = code(alpha, beta, i) t_0 = (alpha + beta) + (2.0 * i); tmp = (((((alpha + beta) * (beta - alpha)) / t_0) / (t_0 + 2.0)) + 1.0) / 2.0; end
code[alpha_, beta_, i_] := Block[{t$95$0 = N[(N[(alpha + beta), $MachinePrecision] + N[(2.0 * i), $MachinePrecision]), $MachinePrecision]}, N[(N[(N[(N[(N[(N[(alpha + beta), $MachinePrecision] * N[(beta - alpha), $MachinePrecision]), $MachinePrecision] / t$95$0), $MachinePrecision] / N[(t$95$0 + 2.0), $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision] / 2.0), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \left(\alpha + \beta\right) + 2 \cdot i\\
\frac{\frac{\frac{\left(\alpha + \beta\right) \cdot \left(\beta - \alpha\right)}{t\_0}}{t\_0 + 2} + 1}{2}
\end{array}
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 19 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (alpha beta i) :precision binary64 (let* ((t_0 (+ (+ alpha beta) (* 2.0 i)))) (/ (+ (/ (/ (* (+ alpha beta) (- beta alpha)) t_0) (+ t_0 2.0)) 1.0) 2.0)))
double code(double alpha, double beta, double i) {
double t_0 = (alpha + beta) + (2.0 * i);
return (((((alpha + beta) * (beta - alpha)) / t_0) / (t_0 + 2.0)) + 1.0) / 2.0;
}
real(8) function code(alpha, beta, i)
real(8), intent (in) :: alpha
real(8), intent (in) :: beta
real(8), intent (in) :: i
real(8) :: t_0
t_0 = (alpha + beta) + (2.0d0 * i)
code = (((((alpha + beta) * (beta - alpha)) / t_0) / (t_0 + 2.0d0)) + 1.0d0) / 2.0d0
end function
public static double code(double alpha, double beta, double i) {
double t_0 = (alpha + beta) + (2.0 * i);
return (((((alpha + beta) * (beta - alpha)) / t_0) / (t_0 + 2.0)) + 1.0) / 2.0;
}
def code(alpha, beta, i): t_0 = (alpha + beta) + (2.0 * i) return (((((alpha + beta) * (beta - alpha)) / t_0) / (t_0 + 2.0)) + 1.0) / 2.0
function code(alpha, beta, i) t_0 = Float64(Float64(alpha + beta) + Float64(2.0 * i)) return Float64(Float64(Float64(Float64(Float64(Float64(alpha + beta) * Float64(beta - alpha)) / t_0) / Float64(t_0 + 2.0)) + 1.0) / 2.0) end
function tmp = code(alpha, beta, i) t_0 = (alpha + beta) + (2.0 * i); tmp = (((((alpha + beta) * (beta - alpha)) / t_0) / (t_0 + 2.0)) + 1.0) / 2.0; end
code[alpha_, beta_, i_] := Block[{t$95$0 = N[(N[(alpha + beta), $MachinePrecision] + N[(2.0 * i), $MachinePrecision]), $MachinePrecision]}, N[(N[(N[(N[(N[(N[(alpha + beta), $MachinePrecision] * N[(beta - alpha), $MachinePrecision]), $MachinePrecision] / t$95$0), $MachinePrecision] / N[(t$95$0 + 2.0), $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision] / 2.0), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \left(\alpha + \beta\right) + 2 \cdot i\\
\frac{\frac{\frac{\left(\alpha + \beta\right) \cdot \left(\beta - \alpha\right)}{t\_0}}{t\_0 + 2} + 1}{2}
\end{array}
\end{array}
(FPCore (alpha beta i)
:precision binary64
(let* ((t_0 (+ (* i 2.0) (+ beta alpha))) (t_1 (fma i 2.0 (+ beta alpha))))
(if (<= (/ (/ (* (- beta alpha) (+ beta alpha)) t_0) (+ t_0 2.0)) -0.99995)
(fma (/ (fma beta 2.0 2.0) alpha) 0.5 (* (/ i alpha) 2.0))
(fma (/ (* (/ (+ beta alpha) t_1) (- beta alpha)) (+ t_1 2.0)) 0.5 0.5))))
double code(double alpha, double beta, double i) {
double t_0 = (i * 2.0) + (beta + alpha);
double t_1 = fma(i, 2.0, (beta + alpha));
double tmp;
if (((((beta - alpha) * (beta + alpha)) / t_0) / (t_0 + 2.0)) <= -0.99995) {
tmp = fma((fma(beta, 2.0, 2.0) / alpha), 0.5, ((i / alpha) * 2.0));
} else {
tmp = fma(((((beta + alpha) / t_1) * (beta - alpha)) / (t_1 + 2.0)), 0.5, 0.5);
}
return tmp;
}
function code(alpha, beta, i) t_0 = Float64(Float64(i * 2.0) + Float64(beta + alpha)) t_1 = fma(i, 2.0, Float64(beta + alpha)) tmp = 0.0 if (Float64(Float64(Float64(Float64(beta - alpha) * Float64(beta + alpha)) / t_0) / Float64(t_0 + 2.0)) <= -0.99995) tmp = fma(Float64(fma(beta, 2.0, 2.0) / alpha), 0.5, Float64(Float64(i / alpha) * 2.0)); else tmp = fma(Float64(Float64(Float64(Float64(beta + alpha) / t_1) * Float64(beta - alpha)) / Float64(t_1 + 2.0)), 0.5, 0.5); end return tmp end
code[alpha_, beta_, i_] := Block[{t$95$0 = N[(N[(i * 2.0), $MachinePrecision] + N[(beta + alpha), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(i * 2.0 + N[(beta + alpha), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[N[(N[(N[(N[(beta - alpha), $MachinePrecision] * N[(beta + alpha), $MachinePrecision]), $MachinePrecision] / t$95$0), $MachinePrecision] / N[(t$95$0 + 2.0), $MachinePrecision]), $MachinePrecision], -0.99995], N[(N[(N[(beta * 2.0 + 2.0), $MachinePrecision] / alpha), $MachinePrecision] * 0.5 + N[(N[(i / alpha), $MachinePrecision] * 2.0), $MachinePrecision]), $MachinePrecision], N[(N[(N[(N[(N[(beta + alpha), $MachinePrecision] / t$95$1), $MachinePrecision] * N[(beta - alpha), $MachinePrecision]), $MachinePrecision] / N[(t$95$1 + 2.0), $MachinePrecision]), $MachinePrecision] * 0.5 + 0.5), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := i \cdot 2 + \left(\beta + \alpha\right)\\
t_1 := \mathsf{fma}\left(i, 2, \beta + \alpha\right)\\
\mathbf{if}\;\frac{\frac{\left(\beta - \alpha\right) \cdot \left(\beta + \alpha\right)}{t\_0}}{t\_0 + 2} \leq -0.99995:\\
\;\;\;\;\mathsf{fma}\left(\frac{\mathsf{fma}\left(\beta, 2, 2\right)}{\alpha}, 0.5, \frac{i}{\alpha} \cdot 2\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(\frac{\frac{\beta + \alpha}{t\_1} \cdot \left(\beta - \alpha\right)}{t\_1 + 2}, 0.5, 0.5\right)\\
\end{array}
\end{array}
if (/.f64 (/.f64 (*.f64 (+.f64 alpha beta) (-.f64 beta alpha)) (+.f64 (+.f64 alpha beta) (*.f64 #s(literal 2 binary64) i))) (+.f64 (+.f64 (+.f64 alpha beta) (*.f64 #s(literal 2 binary64) i)) #s(literal 2 binary64))) < -0.999950000000000006Initial program 3.5%
Taylor expanded in alpha around inf
*-commutativeN/A
distribute-rgt1-inN/A
metadata-evalN/A
mul0-lftN/A
neg-sub0N/A
mul-1-negN/A
remove-double-negN/A
lower-*.f64N/A
lower-/.f64N/A
+-commutativeN/A
lower-+.f64N/A
+-commutativeN/A
lower-fma.f64N/A
*-commutativeN/A
lower-*.f6494.0
Applied rewrites94.0%
Taylor expanded in i around 0
Applied rewrites94.0%
if -0.999950000000000006 < (/.f64 (/.f64 (*.f64 (+.f64 alpha beta) (-.f64 beta alpha)) (+.f64 (+.f64 alpha beta) (*.f64 #s(literal 2 binary64) i))) (+.f64 (+.f64 (+.f64 alpha beta) (*.f64 #s(literal 2 binary64) i)) #s(literal 2 binary64))) Initial program 75.1%
lift-/.f64N/A
clear-numN/A
lower-/.f64N/A
lower-/.f6475.1
lift-+.f64N/A
Applied rewrites99.8%
lift-/.f64N/A
lift-/.f64N/A
associate-/r/N/A
metadata-evalN/A
lift-fma.f64N/A
distribute-rgt-inN/A
metadata-evalN/A
lower-fma.f64N/A
Applied rewrites99.8%
Final simplification98.6%
(FPCore (alpha beta i)
:precision binary64
(let* ((t_0 (fma i 2.0 (+ beta alpha)))
(t_1 (+ (* i 2.0) (+ beta alpha)))
(t_2 (/ (/ (* (- beta alpha) (+ beta alpha)) t_1) (+ t_1 2.0))))
(if (<= t_2 -0.99995)
(fma (/ (fma beta 2.0 2.0) alpha) 0.5 (* (/ i alpha) 2.0))
(if (<= t_2 0.99999999)
(/ (fma (- beta alpha) (/ (+ beta alpha) (* (+ t_0 2.0) t_0)) 1.0) 2.0)
(fma
(* (/ beta (+ (fma i 2.0 beta) 2.0)) (/ beta (fma i 2.0 beta)))
0.5
0.5)))))
double code(double alpha, double beta, double i) {
double t_0 = fma(i, 2.0, (beta + alpha));
double t_1 = (i * 2.0) + (beta + alpha);
double t_2 = (((beta - alpha) * (beta + alpha)) / t_1) / (t_1 + 2.0);
double tmp;
if (t_2 <= -0.99995) {
tmp = fma((fma(beta, 2.0, 2.0) / alpha), 0.5, ((i / alpha) * 2.0));
} else if (t_2 <= 0.99999999) {
tmp = fma((beta - alpha), ((beta + alpha) / ((t_0 + 2.0) * t_0)), 1.0) / 2.0;
} else {
tmp = fma(((beta / (fma(i, 2.0, beta) + 2.0)) * (beta / fma(i, 2.0, beta))), 0.5, 0.5);
}
return tmp;
}
function code(alpha, beta, i) t_0 = fma(i, 2.0, Float64(beta + alpha)) t_1 = Float64(Float64(i * 2.0) + Float64(beta + alpha)) t_2 = Float64(Float64(Float64(Float64(beta - alpha) * Float64(beta + alpha)) / t_1) / Float64(t_1 + 2.0)) tmp = 0.0 if (t_2 <= -0.99995) tmp = fma(Float64(fma(beta, 2.0, 2.0) / alpha), 0.5, Float64(Float64(i / alpha) * 2.0)); elseif (t_2 <= 0.99999999) tmp = Float64(fma(Float64(beta - alpha), Float64(Float64(beta + alpha) / Float64(Float64(t_0 + 2.0) * t_0)), 1.0) / 2.0); else tmp = fma(Float64(Float64(beta / Float64(fma(i, 2.0, beta) + 2.0)) * Float64(beta / fma(i, 2.0, beta))), 0.5, 0.5); end return tmp end
code[alpha_, beta_, i_] := Block[{t$95$0 = N[(i * 2.0 + N[(beta + alpha), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(N[(i * 2.0), $MachinePrecision] + N[(beta + alpha), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(N[(N[(N[(beta - alpha), $MachinePrecision] * N[(beta + alpha), $MachinePrecision]), $MachinePrecision] / t$95$1), $MachinePrecision] / N[(t$95$1 + 2.0), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$2, -0.99995], N[(N[(N[(beta * 2.0 + 2.0), $MachinePrecision] / alpha), $MachinePrecision] * 0.5 + N[(N[(i / alpha), $MachinePrecision] * 2.0), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$2, 0.99999999], N[(N[(N[(beta - alpha), $MachinePrecision] * N[(N[(beta + alpha), $MachinePrecision] / N[(N[(t$95$0 + 2.0), $MachinePrecision] * t$95$0), $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision] / 2.0), $MachinePrecision], N[(N[(N[(beta / N[(N[(i * 2.0 + beta), $MachinePrecision] + 2.0), $MachinePrecision]), $MachinePrecision] * N[(beta / N[(i * 2.0 + beta), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * 0.5 + 0.5), $MachinePrecision]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \mathsf{fma}\left(i, 2, \beta + \alpha\right)\\
t_1 := i \cdot 2 + \left(\beta + \alpha\right)\\
t_2 := \frac{\frac{\left(\beta - \alpha\right) \cdot \left(\beta + \alpha\right)}{t\_1}}{t\_1 + 2}\\
\mathbf{if}\;t\_2 \leq -0.99995:\\
\;\;\;\;\mathsf{fma}\left(\frac{\mathsf{fma}\left(\beta, 2, 2\right)}{\alpha}, 0.5, \frac{i}{\alpha} \cdot 2\right)\\
\mathbf{elif}\;t\_2 \leq 0.99999999:\\
\;\;\;\;\frac{\mathsf{fma}\left(\beta - \alpha, \frac{\beta + \alpha}{\left(t\_0 + 2\right) \cdot t\_0}, 1\right)}{2}\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(\frac{\beta}{\mathsf{fma}\left(i, 2, \beta\right) + 2} \cdot \frac{\beta}{\mathsf{fma}\left(i, 2, \beta\right)}, 0.5, 0.5\right)\\
\end{array}
\end{array}
if (/.f64 (/.f64 (*.f64 (+.f64 alpha beta) (-.f64 beta alpha)) (+.f64 (+.f64 alpha beta) (*.f64 #s(literal 2 binary64) i))) (+.f64 (+.f64 (+.f64 alpha beta) (*.f64 #s(literal 2 binary64) i)) #s(literal 2 binary64))) < -0.999950000000000006Initial program 3.5%
Taylor expanded in alpha around inf
*-commutativeN/A
distribute-rgt1-inN/A
metadata-evalN/A
mul0-lftN/A
neg-sub0N/A
mul-1-negN/A
remove-double-negN/A
lower-*.f64N/A
lower-/.f64N/A
+-commutativeN/A
lower-+.f64N/A
+-commutativeN/A
lower-fma.f64N/A
*-commutativeN/A
lower-*.f6494.0
Applied rewrites94.0%
Taylor expanded in i around 0
Applied rewrites94.0%
if -0.999950000000000006 < (/.f64 (/.f64 (*.f64 (+.f64 alpha beta) (-.f64 beta alpha)) (+.f64 (+.f64 alpha beta) (*.f64 #s(literal 2 binary64) i))) (+.f64 (+.f64 (+.f64 alpha beta) (*.f64 #s(literal 2 binary64) i)) #s(literal 2 binary64))) < 0.99999998999999995Initial program 99.7%
lift-+.f64N/A
lift-/.f64N/A
lift-/.f64N/A
associate-/l/N/A
lift-*.f64N/A
*-commutativeN/A
associate-/l*N/A
lower-fma.f64N/A
Applied rewrites99.7%
if 0.99999998999999995 < (/.f64 (/.f64 (*.f64 (+.f64 alpha beta) (-.f64 beta alpha)) (+.f64 (+.f64 alpha beta) (*.f64 #s(literal 2 binary64) i))) (+.f64 (+.f64 (+.f64 alpha beta) (*.f64 #s(literal 2 binary64) i)) #s(literal 2 binary64))) Initial program 28.0%
Taylor expanded in alpha around 0
*-commutativeN/A
lower-*.f64N/A
+-commutativeN/A
unpow2N/A
times-fracN/A
lower-fma.f64N/A
lower-/.f64N/A
+-commutativeN/A
lower-+.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower-/.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f6499.9
Applied rewrites99.9%
Applied rewrites100.0%
Final simplification98.6%
(FPCore (alpha beta i)
:precision binary64
(let* ((t_0 (+ (* i 2.0) (+ beta alpha)))
(t_1 (/ (/ (* (- beta alpha) (+ beta alpha)) t_0) (+ t_0 2.0))))
(if (<= t_1 -0.99995)
(fma (/ (fma beta 2.0 2.0) alpha) 0.5 (* (/ i alpha) 2.0))
(if (<= t_1 0.002)
(fma
-0.5
(/ (* (/ alpha (fma 2.0 i alpha)) alpha) (+ (fma 2.0 i alpha) 2.0))
0.5)
(fma (/ (- -2.0 (fma 4.0 i (* 2.0 alpha))) beta) 0.5 1.0)))))
double code(double alpha, double beta, double i) {
double t_0 = (i * 2.0) + (beta + alpha);
double t_1 = (((beta - alpha) * (beta + alpha)) / t_0) / (t_0 + 2.0);
double tmp;
if (t_1 <= -0.99995) {
tmp = fma((fma(beta, 2.0, 2.0) / alpha), 0.5, ((i / alpha) * 2.0));
} else if (t_1 <= 0.002) {
tmp = fma(-0.5, (((alpha / fma(2.0, i, alpha)) * alpha) / (fma(2.0, i, alpha) + 2.0)), 0.5);
} else {
tmp = fma(((-2.0 - fma(4.0, i, (2.0 * alpha))) / beta), 0.5, 1.0);
}
return tmp;
}
function code(alpha, beta, i) t_0 = Float64(Float64(i * 2.0) + Float64(beta + alpha)) t_1 = Float64(Float64(Float64(Float64(beta - alpha) * Float64(beta + alpha)) / t_0) / Float64(t_0 + 2.0)) tmp = 0.0 if (t_1 <= -0.99995) tmp = fma(Float64(fma(beta, 2.0, 2.0) / alpha), 0.5, Float64(Float64(i / alpha) * 2.0)); elseif (t_1 <= 0.002) tmp = fma(-0.5, Float64(Float64(Float64(alpha / fma(2.0, i, alpha)) * alpha) / Float64(fma(2.0, i, alpha) + 2.0)), 0.5); else tmp = fma(Float64(Float64(-2.0 - fma(4.0, i, Float64(2.0 * alpha))) / beta), 0.5, 1.0); end return tmp end
code[alpha_, beta_, i_] := Block[{t$95$0 = N[(N[(i * 2.0), $MachinePrecision] + N[(beta + alpha), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(N[(N[(N[(beta - alpha), $MachinePrecision] * N[(beta + alpha), $MachinePrecision]), $MachinePrecision] / t$95$0), $MachinePrecision] / N[(t$95$0 + 2.0), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$1, -0.99995], N[(N[(N[(beta * 2.0 + 2.0), $MachinePrecision] / alpha), $MachinePrecision] * 0.5 + N[(N[(i / alpha), $MachinePrecision] * 2.0), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$1, 0.002], N[(-0.5 * N[(N[(N[(alpha / N[(2.0 * i + alpha), $MachinePrecision]), $MachinePrecision] * alpha), $MachinePrecision] / N[(N[(2.0 * i + alpha), $MachinePrecision] + 2.0), $MachinePrecision]), $MachinePrecision] + 0.5), $MachinePrecision], N[(N[(N[(-2.0 - N[(4.0 * i + N[(2.0 * alpha), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / beta), $MachinePrecision] * 0.5 + 1.0), $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := i \cdot 2 + \left(\beta + \alpha\right)\\
t_1 := \frac{\frac{\left(\beta - \alpha\right) \cdot \left(\beta + \alpha\right)}{t\_0}}{t\_0 + 2}\\
\mathbf{if}\;t\_1 \leq -0.99995:\\
\;\;\;\;\mathsf{fma}\left(\frac{\mathsf{fma}\left(\beta, 2, 2\right)}{\alpha}, 0.5, \frac{i}{\alpha} \cdot 2\right)\\
\mathbf{elif}\;t\_1 \leq 0.002:\\
\;\;\;\;\mathsf{fma}\left(-0.5, \frac{\frac{\alpha}{\mathsf{fma}\left(2, i, \alpha\right)} \cdot \alpha}{\mathsf{fma}\left(2, i, \alpha\right) + 2}, 0.5\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(\frac{-2 - \mathsf{fma}\left(4, i, 2 \cdot \alpha\right)}{\beta}, 0.5, 1\right)\\
\end{array}
\end{array}
if (/.f64 (/.f64 (*.f64 (+.f64 alpha beta) (-.f64 beta alpha)) (+.f64 (+.f64 alpha beta) (*.f64 #s(literal 2 binary64) i))) (+.f64 (+.f64 (+.f64 alpha beta) (*.f64 #s(literal 2 binary64) i)) #s(literal 2 binary64))) < -0.999950000000000006Initial program 3.5%
Taylor expanded in alpha around inf
*-commutativeN/A
distribute-rgt1-inN/A
metadata-evalN/A
mul0-lftN/A
neg-sub0N/A
mul-1-negN/A
remove-double-negN/A
lower-*.f64N/A
lower-/.f64N/A
+-commutativeN/A
lower-+.f64N/A
+-commutativeN/A
lower-fma.f64N/A
*-commutativeN/A
lower-*.f6494.0
Applied rewrites94.0%
Taylor expanded in i around 0
Applied rewrites94.0%
if -0.999950000000000006 < (/.f64 (/.f64 (*.f64 (+.f64 alpha beta) (-.f64 beta alpha)) (+.f64 (+.f64 alpha beta) (*.f64 #s(literal 2 binary64) i))) (+.f64 (+.f64 (+.f64 alpha beta) (*.f64 #s(literal 2 binary64) i)) #s(literal 2 binary64))) < 2e-3Initial program 99.7%
lift-/.f64N/A
clear-numN/A
lower-/.f64N/A
lower-/.f6499.7
lift-+.f64N/A
Applied rewrites99.7%
lift-/.f64N/A
lift-/.f64N/A
associate-/r/N/A
metadata-evalN/A
lift-fma.f64N/A
distribute-rgt-inN/A
metadata-evalN/A
lower-fma.f64N/A
Applied rewrites99.7%
Taylor expanded in beta around 0
+-commutativeN/A
associate-*r/N/A
times-fracN/A
lower-fma.f64N/A
lower-/.f64N/A
+-commutativeN/A
lower-+.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower-/.f64N/A
unpow2N/A
lower-*.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f6499.2
Applied rewrites99.2%
Applied rewrites99.2%
if 2e-3 < (/.f64 (/.f64 (*.f64 (+.f64 alpha beta) (-.f64 beta alpha)) (+.f64 (+.f64 alpha beta) (*.f64 #s(literal 2 binary64) i))) (+.f64 (+.f64 (+.f64 alpha beta) (*.f64 #s(literal 2 binary64) i)) #s(literal 2 binary64))) Initial program 30.0%
Taylor expanded in beta around inf
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower-/.f64N/A
associate--r+N/A
distribute-rgt1-inN/A
metadata-evalN/A
mul0-lftN/A
metadata-evalN/A
lower--.f64N/A
+-commutativeN/A
lower-fma.f64N/A
lower-*.f6494.4
Applied rewrites94.4%
Final simplification96.8%
(FPCore (alpha beta i)
:precision binary64
(let* ((t_0 (+ (* i 2.0) (+ beta alpha)))
(t_1 (/ (/ (* (- beta alpha) (+ beta alpha)) t_0) (+ t_0 2.0))))
(if (<= t_1 -0.99995)
(fma (/ (fma beta 2.0 2.0) alpha) 0.5 (* (/ i alpha) 2.0))
(if (<= t_1 0.002)
(fma
0.5
(/ (* alpha alpha) (* (- -2.0 (fma i 2.0 alpha)) (fma i 2.0 alpha)))
0.5)
(fma (/ (- -2.0 (fma 4.0 i (* 2.0 alpha))) beta) 0.5 1.0)))))
double code(double alpha, double beta, double i) {
double t_0 = (i * 2.0) + (beta + alpha);
double t_1 = (((beta - alpha) * (beta + alpha)) / t_0) / (t_0 + 2.0);
double tmp;
if (t_1 <= -0.99995) {
tmp = fma((fma(beta, 2.0, 2.0) / alpha), 0.5, ((i / alpha) * 2.0));
} else if (t_1 <= 0.002) {
tmp = fma(0.5, ((alpha * alpha) / ((-2.0 - fma(i, 2.0, alpha)) * fma(i, 2.0, alpha))), 0.5);
} else {
tmp = fma(((-2.0 - fma(4.0, i, (2.0 * alpha))) / beta), 0.5, 1.0);
}
return tmp;
}
function code(alpha, beta, i) t_0 = Float64(Float64(i * 2.0) + Float64(beta + alpha)) t_1 = Float64(Float64(Float64(Float64(beta - alpha) * Float64(beta + alpha)) / t_0) / Float64(t_0 + 2.0)) tmp = 0.0 if (t_1 <= -0.99995) tmp = fma(Float64(fma(beta, 2.0, 2.0) / alpha), 0.5, Float64(Float64(i / alpha) * 2.0)); elseif (t_1 <= 0.002) tmp = fma(0.5, Float64(Float64(alpha * alpha) / Float64(Float64(-2.0 - fma(i, 2.0, alpha)) * fma(i, 2.0, alpha))), 0.5); else tmp = fma(Float64(Float64(-2.0 - fma(4.0, i, Float64(2.0 * alpha))) / beta), 0.5, 1.0); end return tmp end
code[alpha_, beta_, i_] := Block[{t$95$0 = N[(N[(i * 2.0), $MachinePrecision] + N[(beta + alpha), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(N[(N[(N[(beta - alpha), $MachinePrecision] * N[(beta + alpha), $MachinePrecision]), $MachinePrecision] / t$95$0), $MachinePrecision] / N[(t$95$0 + 2.0), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$1, -0.99995], N[(N[(N[(beta * 2.0 + 2.0), $MachinePrecision] / alpha), $MachinePrecision] * 0.5 + N[(N[(i / alpha), $MachinePrecision] * 2.0), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$1, 0.002], N[(0.5 * N[(N[(alpha * alpha), $MachinePrecision] / N[(N[(-2.0 - N[(i * 2.0 + alpha), $MachinePrecision]), $MachinePrecision] * N[(i * 2.0 + alpha), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + 0.5), $MachinePrecision], N[(N[(N[(-2.0 - N[(4.0 * i + N[(2.0 * alpha), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / beta), $MachinePrecision] * 0.5 + 1.0), $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := i \cdot 2 + \left(\beta + \alpha\right)\\
t_1 := \frac{\frac{\left(\beta - \alpha\right) \cdot \left(\beta + \alpha\right)}{t\_0}}{t\_0 + 2}\\
\mathbf{if}\;t\_1 \leq -0.99995:\\
\;\;\;\;\mathsf{fma}\left(\frac{\mathsf{fma}\left(\beta, 2, 2\right)}{\alpha}, 0.5, \frac{i}{\alpha} \cdot 2\right)\\
\mathbf{elif}\;t\_1 \leq 0.002:\\
\;\;\;\;\mathsf{fma}\left(0.5, \frac{\alpha \cdot \alpha}{\left(-2 - \mathsf{fma}\left(i, 2, \alpha\right)\right) \cdot \mathsf{fma}\left(i, 2, \alpha\right)}, 0.5\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(\frac{-2 - \mathsf{fma}\left(4, i, 2 \cdot \alpha\right)}{\beta}, 0.5, 1\right)\\
\end{array}
\end{array}
if (/.f64 (/.f64 (*.f64 (+.f64 alpha beta) (-.f64 beta alpha)) (+.f64 (+.f64 alpha beta) (*.f64 #s(literal 2 binary64) i))) (+.f64 (+.f64 (+.f64 alpha beta) (*.f64 #s(literal 2 binary64) i)) #s(literal 2 binary64))) < -0.999950000000000006Initial program 3.5%
Taylor expanded in alpha around inf
*-commutativeN/A
distribute-rgt1-inN/A
metadata-evalN/A
mul0-lftN/A
neg-sub0N/A
mul-1-negN/A
remove-double-negN/A
lower-*.f64N/A
lower-/.f64N/A
+-commutativeN/A
lower-+.f64N/A
+-commutativeN/A
lower-fma.f64N/A
*-commutativeN/A
lower-*.f6494.0
Applied rewrites94.0%
Taylor expanded in i around 0
Applied rewrites94.0%
if -0.999950000000000006 < (/.f64 (/.f64 (*.f64 (+.f64 alpha beta) (-.f64 beta alpha)) (+.f64 (+.f64 alpha beta) (*.f64 #s(literal 2 binary64) i))) (+.f64 (+.f64 (+.f64 alpha beta) (*.f64 #s(literal 2 binary64) i)) #s(literal 2 binary64))) < 2e-3Initial program 99.7%
Taylor expanded in beta around 0
+-commutativeN/A
distribute-lft-inN/A
metadata-evalN/A
lower-fma.f64N/A
Applied rewrites99.2%
if 2e-3 < (/.f64 (/.f64 (*.f64 (+.f64 alpha beta) (-.f64 beta alpha)) (+.f64 (+.f64 alpha beta) (*.f64 #s(literal 2 binary64) i))) (+.f64 (+.f64 (+.f64 alpha beta) (*.f64 #s(literal 2 binary64) i)) #s(literal 2 binary64))) Initial program 30.0%
Taylor expanded in beta around inf
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower-/.f64N/A
associate--r+N/A
distribute-rgt1-inN/A
metadata-evalN/A
mul0-lftN/A
metadata-evalN/A
lower--.f64N/A
+-commutativeN/A
lower-fma.f64N/A
lower-*.f6494.4
Applied rewrites94.4%
Final simplification96.8%
(FPCore (alpha beta i)
:precision binary64
(let* ((t_0 (+ (* i 2.0) (+ beta alpha)))
(t_1 (/ (/ (* (- beta alpha) (+ beta alpha)) t_0) (+ t_0 2.0))))
(if (<= t_1 -0.5)
(fma (/ (fma beta 2.0 2.0) alpha) 0.5 (* (/ i alpha) 2.0))
(if (<= t_1 5e-84)
0.5
(fma (/ (- beta alpha) (+ 2.0 (+ beta alpha))) 0.5 0.5)))))
double code(double alpha, double beta, double i) {
double t_0 = (i * 2.0) + (beta + alpha);
double t_1 = (((beta - alpha) * (beta + alpha)) / t_0) / (t_0 + 2.0);
double tmp;
if (t_1 <= -0.5) {
tmp = fma((fma(beta, 2.0, 2.0) / alpha), 0.5, ((i / alpha) * 2.0));
} else if (t_1 <= 5e-84) {
tmp = 0.5;
} else {
tmp = fma(((beta - alpha) / (2.0 + (beta + alpha))), 0.5, 0.5);
}
return tmp;
}
function code(alpha, beta, i) t_0 = Float64(Float64(i * 2.0) + Float64(beta + alpha)) t_1 = Float64(Float64(Float64(Float64(beta - alpha) * Float64(beta + alpha)) / t_0) / Float64(t_0 + 2.0)) tmp = 0.0 if (t_1 <= -0.5) tmp = fma(Float64(fma(beta, 2.0, 2.0) / alpha), 0.5, Float64(Float64(i / alpha) * 2.0)); elseif (t_1 <= 5e-84) tmp = 0.5; else tmp = fma(Float64(Float64(beta - alpha) / Float64(2.0 + Float64(beta + alpha))), 0.5, 0.5); end return tmp end
code[alpha_, beta_, i_] := Block[{t$95$0 = N[(N[(i * 2.0), $MachinePrecision] + N[(beta + alpha), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(N[(N[(N[(beta - alpha), $MachinePrecision] * N[(beta + alpha), $MachinePrecision]), $MachinePrecision] / t$95$0), $MachinePrecision] / N[(t$95$0 + 2.0), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$1, -0.5], N[(N[(N[(beta * 2.0 + 2.0), $MachinePrecision] / alpha), $MachinePrecision] * 0.5 + N[(N[(i / alpha), $MachinePrecision] * 2.0), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$1, 5e-84], 0.5, N[(N[(N[(beta - alpha), $MachinePrecision] / N[(2.0 + N[(beta + alpha), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * 0.5 + 0.5), $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := i \cdot 2 + \left(\beta + \alpha\right)\\
t_1 := \frac{\frac{\left(\beta - \alpha\right) \cdot \left(\beta + \alpha\right)}{t\_0}}{t\_0 + 2}\\
\mathbf{if}\;t\_1 \leq -0.5:\\
\;\;\;\;\mathsf{fma}\left(\frac{\mathsf{fma}\left(\beta, 2, 2\right)}{\alpha}, 0.5, \frac{i}{\alpha} \cdot 2\right)\\
\mathbf{elif}\;t\_1 \leq 5 \cdot 10^{-84}:\\
\;\;\;\;0.5\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(\frac{\beta - \alpha}{2 + \left(\beta + \alpha\right)}, 0.5, 0.5\right)\\
\end{array}
\end{array}
if (/.f64 (/.f64 (*.f64 (+.f64 alpha beta) (-.f64 beta alpha)) (+.f64 (+.f64 alpha beta) (*.f64 #s(literal 2 binary64) i))) (+.f64 (+.f64 (+.f64 alpha beta) (*.f64 #s(literal 2 binary64) i)) #s(literal 2 binary64))) < -0.5Initial program 6.4%
Taylor expanded in alpha around inf
*-commutativeN/A
distribute-rgt1-inN/A
metadata-evalN/A
mul0-lftN/A
neg-sub0N/A
mul-1-negN/A
remove-double-negN/A
lower-*.f64N/A
lower-/.f64N/A
+-commutativeN/A
lower-+.f64N/A
+-commutativeN/A
lower-fma.f64N/A
*-commutativeN/A
lower-*.f6492.0
Applied rewrites92.0%
Taylor expanded in i around 0
Applied rewrites92.0%
if -0.5 < (/.f64 (/.f64 (*.f64 (+.f64 alpha beta) (-.f64 beta alpha)) (+.f64 (+.f64 alpha beta) (*.f64 #s(literal 2 binary64) i))) (+.f64 (+.f64 (+.f64 alpha beta) (*.f64 #s(literal 2 binary64) i)) #s(literal 2 binary64))) < 5.0000000000000002e-84Initial program 100.0%
Taylor expanded in i around inf
Applied rewrites97.6%
if 5.0000000000000002e-84 < (/.f64 (/.f64 (*.f64 (+.f64 alpha beta) (-.f64 beta alpha)) (+.f64 (+.f64 alpha beta) (*.f64 #s(literal 2 binary64) i))) (+.f64 (+.f64 (+.f64 alpha beta) (*.f64 #s(literal 2 binary64) i)) #s(literal 2 binary64))) Initial program 42.1%
lift-/.f64N/A
clear-numN/A
lower-/.f64N/A
lower-/.f6442.1
lift-+.f64N/A
Applied rewrites99.9%
lift-/.f64N/A
lift-/.f64N/A
associate-/r/N/A
metadata-evalN/A
lift-fma.f64N/A
distribute-rgt-inN/A
metadata-evalN/A
lower-fma.f64N/A
Applied rewrites99.9%
Taylor expanded in i around 0
lower-/.f64N/A
lower--.f64N/A
+-commutativeN/A
lower-+.f64N/A
+-commutativeN/A
lower-+.f6494.9
Applied rewrites94.9%
Final simplification95.5%
(FPCore (alpha beta i)
:precision binary64
(let* ((t_0 (+ (* i 2.0) (+ beta alpha)))
(t_1 (/ (/ (* (- beta alpha) (+ beta alpha)) t_0) (+ t_0 2.0))))
(if (<= t_1 -0.5)
(/ (fma (fma 4.0 i (* 2.0 beta)) 0.5 1.0) alpha)
(if (<= t_1 5e-84)
0.5
(fma (/ (- beta alpha) (+ 2.0 (+ beta alpha))) 0.5 0.5)))))
double code(double alpha, double beta, double i) {
double t_0 = (i * 2.0) + (beta + alpha);
double t_1 = (((beta - alpha) * (beta + alpha)) / t_0) / (t_0 + 2.0);
double tmp;
if (t_1 <= -0.5) {
tmp = fma(fma(4.0, i, (2.0 * beta)), 0.5, 1.0) / alpha;
} else if (t_1 <= 5e-84) {
tmp = 0.5;
} else {
tmp = fma(((beta - alpha) / (2.0 + (beta + alpha))), 0.5, 0.5);
}
return tmp;
}
function code(alpha, beta, i) t_0 = Float64(Float64(i * 2.0) + Float64(beta + alpha)) t_1 = Float64(Float64(Float64(Float64(beta - alpha) * Float64(beta + alpha)) / t_0) / Float64(t_0 + 2.0)) tmp = 0.0 if (t_1 <= -0.5) tmp = Float64(fma(fma(4.0, i, Float64(2.0 * beta)), 0.5, 1.0) / alpha); elseif (t_1 <= 5e-84) tmp = 0.5; else tmp = fma(Float64(Float64(beta - alpha) / Float64(2.0 + Float64(beta + alpha))), 0.5, 0.5); end return tmp end
code[alpha_, beta_, i_] := Block[{t$95$0 = N[(N[(i * 2.0), $MachinePrecision] + N[(beta + alpha), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(N[(N[(N[(beta - alpha), $MachinePrecision] * N[(beta + alpha), $MachinePrecision]), $MachinePrecision] / t$95$0), $MachinePrecision] / N[(t$95$0 + 2.0), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$1, -0.5], N[(N[(N[(4.0 * i + N[(2.0 * beta), $MachinePrecision]), $MachinePrecision] * 0.5 + 1.0), $MachinePrecision] / alpha), $MachinePrecision], If[LessEqual[t$95$1, 5e-84], 0.5, N[(N[(N[(beta - alpha), $MachinePrecision] / N[(2.0 + N[(beta + alpha), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * 0.5 + 0.5), $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := i \cdot 2 + \left(\beta + \alpha\right)\\
t_1 := \frac{\frac{\left(\beta - \alpha\right) \cdot \left(\beta + \alpha\right)}{t\_0}}{t\_0 + 2}\\
\mathbf{if}\;t\_1 \leq -0.5:\\
\;\;\;\;\frac{\mathsf{fma}\left(\mathsf{fma}\left(4, i, 2 \cdot \beta\right), 0.5, 1\right)}{\alpha}\\
\mathbf{elif}\;t\_1 \leq 5 \cdot 10^{-84}:\\
\;\;\;\;0.5\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(\frac{\beta - \alpha}{2 + \left(\beta + \alpha\right)}, 0.5, 0.5\right)\\
\end{array}
\end{array}
if (/.f64 (/.f64 (*.f64 (+.f64 alpha beta) (-.f64 beta alpha)) (+.f64 (+.f64 alpha beta) (*.f64 #s(literal 2 binary64) i))) (+.f64 (+.f64 (+.f64 alpha beta) (*.f64 #s(literal 2 binary64) i)) #s(literal 2 binary64))) < -0.5Initial program 6.4%
Taylor expanded in alpha around inf
*-commutativeN/A
distribute-rgt1-inN/A
metadata-evalN/A
mul0-lftN/A
neg-sub0N/A
mul-1-negN/A
remove-double-negN/A
lower-*.f64N/A
lower-/.f64N/A
+-commutativeN/A
lower-+.f64N/A
+-commutativeN/A
lower-fma.f64N/A
*-commutativeN/A
lower-*.f6492.0
Applied rewrites92.0%
Applied rewrites92.0%
if -0.5 < (/.f64 (/.f64 (*.f64 (+.f64 alpha beta) (-.f64 beta alpha)) (+.f64 (+.f64 alpha beta) (*.f64 #s(literal 2 binary64) i))) (+.f64 (+.f64 (+.f64 alpha beta) (*.f64 #s(literal 2 binary64) i)) #s(literal 2 binary64))) < 5.0000000000000002e-84Initial program 100.0%
Taylor expanded in i around inf
Applied rewrites97.6%
if 5.0000000000000002e-84 < (/.f64 (/.f64 (*.f64 (+.f64 alpha beta) (-.f64 beta alpha)) (+.f64 (+.f64 alpha beta) (*.f64 #s(literal 2 binary64) i))) (+.f64 (+.f64 (+.f64 alpha beta) (*.f64 #s(literal 2 binary64) i)) #s(literal 2 binary64))) Initial program 42.1%
lift-/.f64N/A
clear-numN/A
lower-/.f64N/A
lower-/.f6442.1
lift-+.f64N/A
Applied rewrites99.9%
lift-/.f64N/A
lift-/.f64N/A
associate-/r/N/A
metadata-evalN/A
lift-fma.f64N/A
distribute-rgt-inN/A
metadata-evalN/A
lower-fma.f64N/A
Applied rewrites99.9%
Taylor expanded in i around 0
lower-/.f64N/A
lower--.f64N/A
+-commutativeN/A
lower-+.f64N/A
+-commutativeN/A
lower-+.f6494.9
Applied rewrites94.9%
Final simplification95.5%
(FPCore (alpha beta i)
:precision binary64
(let* ((t_0 (+ (* i 2.0) (+ beta alpha)))
(t_1 (/ (/ (* (- beta alpha) (+ beta alpha)) t_0) (+ t_0 2.0))))
(if (<= t_1 -0.5)
(* (/ 0.5 alpha) (fma 4.0 i (fma beta 2.0 2.0)))
(if (<= t_1 5e-84)
0.5
(fma (/ (- beta alpha) (+ 2.0 (+ beta alpha))) 0.5 0.5)))))
double code(double alpha, double beta, double i) {
double t_0 = (i * 2.0) + (beta + alpha);
double t_1 = (((beta - alpha) * (beta + alpha)) / t_0) / (t_0 + 2.0);
double tmp;
if (t_1 <= -0.5) {
tmp = (0.5 / alpha) * fma(4.0, i, fma(beta, 2.0, 2.0));
} else if (t_1 <= 5e-84) {
tmp = 0.5;
} else {
tmp = fma(((beta - alpha) / (2.0 + (beta + alpha))), 0.5, 0.5);
}
return tmp;
}
function code(alpha, beta, i) t_0 = Float64(Float64(i * 2.0) + Float64(beta + alpha)) t_1 = Float64(Float64(Float64(Float64(beta - alpha) * Float64(beta + alpha)) / t_0) / Float64(t_0 + 2.0)) tmp = 0.0 if (t_1 <= -0.5) tmp = Float64(Float64(0.5 / alpha) * fma(4.0, i, fma(beta, 2.0, 2.0))); elseif (t_1 <= 5e-84) tmp = 0.5; else tmp = fma(Float64(Float64(beta - alpha) / Float64(2.0 + Float64(beta + alpha))), 0.5, 0.5); end return tmp end
code[alpha_, beta_, i_] := Block[{t$95$0 = N[(N[(i * 2.0), $MachinePrecision] + N[(beta + alpha), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(N[(N[(N[(beta - alpha), $MachinePrecision] * N[(beta + alpha), $MachinePrecision]), $MachinePrecision] / t$95$0), $MachinePrecision] / N[(t$95$0 + 2.0), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$1, -0.5], N[(N[(0.5 / alpha), $MachinePrecision] * N[(4.0 * i + N[(beta * 2.0 + 2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$1, 5e-84], 0.5, N[(N[(N[(beta - alpha), $MachinePrecision] / N[(2.0 + N[(beta + alpha), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * 0.5 + 0.5), $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := i \cdot 2 + \left(\beta + \alpha\right)\\
t_1 := \frac{\frac{\left(\beta - \alpha\right) \cdot \left(\beta + \alpha\right)}{t\_0}}{t\_0 + 2}\\
\mathbf{if}\;t\_1 \leq -0.5:\\
\;\;\;\;\frac{0.5}{\alpha} \cdot \mathsf{fma}\left(4, i, \mathsf{fma}\left(\beta, 2, 2\right)\right)\\
\mathbf{elif}\;t\_1 \leq 5 \cdot 10^{-84}:\\
\;\;\;\;0.5\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(\frac{\beta - \alpha}{2 + \left(\beta + \alpha\right)}, 0.5, 0.5\right)\\
\end{array}
\end{array}
if (/.f64 (/.f64 (*.f64 (+.f64 alpha beta) (-.f64 beta alpha)) (+.f64 (+.f64 alpha beta) (*.f64 #s(literal 2 binary64) i))) (+.f64 (+.f64 (+.f64 alpha beta) (*.f64 #s(literal 2 binary64) i)) #s(literal 2 binary64))) < -0.5Initial program 6.4%
Taylor expanded in alpha around inf
*-commutativeN/A
distribute-rgt1-inN/A
metadata-evalN/A
mul0-lftN/A
neg-sub0N/A
mul-1-negN/A
remove-double-negN/A
lower-*.f64N/A
lower-/.f64N/A
+-commutativeN/A
lower-+.f64N/A
+-commutativeN/A
lower-fma.f64N/A
*-commutativeN/A
lower-*.f6492.0
Applied rewrites92.0%
Applied rewrites91.8%
if -0.5 < (/.f64 (/.f64 (*.f64 (+.f64 alpha beta) (-.f64 beta alpha)) (+.f64 (+.f64 alpha beta) (*.f64 #s(literal 2 binary64) i))) (+.f64 (+.f64 (+.f64 alpha beta) (*.f64 #s(literal 2 binary64) i)) #s(literal 2 binary64))) < 5.0000000000000002e-84Initial program 100.0%
Taylor expanded in i around inf
Applied rewrites97.6%
if 5.0000000000000002e-84 < (/.f64 (/.f64 (*.f64 (+.f64 alpha beta) (-.f64 beta alpha)) (+.f64 (+.f64 alpha beta) (*.f64 #s(literal 2 binary64) i))) (+.f64 (+.f64 (+.f64 alpha beta) (*.f64 #s(literal 2 binary64) i)) #s(literal 2 binary64))) Initial program 42.1%
lift-/.f64N/A
clear-numN/A
lower-/.f64N/A
lower-/.f6442.1
lift-+.f64N/A
Applied rewrites99.9%
lift-/.f64N/A
lift-/.f64N/A
associate-/r/N/A
metadata-evalN/A
lift-fma.f64N/A
distribute-rgt-inN/A
metadata-evalN/A
lower-fma.f64N/A
Applied rewrites99.9%
Taylor expanded in i around 0
lower-/.f64N/A
lower--.f64N/A
+-commutativeN/A
lower-+.f64N/A
+-commutativeN/A
lower-+.f6494.9
Applied rewrites94.9%
Final simplification95.5%
(FPCore (alpha beta i)
:precision binary64
(let* ((t_0 (+ (* i 2.0) (+ beta alpha)))
(t_1 (/ (/ (* (- beta alpha) (+ beta alpha)) t_0) (+ t_0 2.0))))
(if (<= t_1 -0.99992)
(* (/ (fma beta 2.0 2.0) alpha) 0.5)
(if (<= t_1 5e-84)
0.5
(fma (/ (- beta alpha) (+ 2.0 (+ beta alpha))) 0.5 0.5)))))
double code(double alpha, double beta, double i) {
double t_0 = (i * 2.0) + (beta + alpha);
double t_1 = (((beta - alpha) * (beta + alpha)) / t_0) / (t_0 + 2.0);
double tmp;
if (t_1 <= -0.99992) {
tmp = (fma(beta, 2.0, 2.0) / alpha) * 0.5;
} else if (t_1 <= 5e-84) {
tmp = 0.5;
} else {
tmp = fma(((beta - alpha) / (2.0 + (beta + alpha))), 0.5, 0.5);
}
return tmp;
}
function code(alpha, beta, i) t_0 = Float64(Float64(i * 2.0) + Float64(beta + alpha)) t_1 = Float64(Float64(Float64(Float64(beta - alpha) * Float64(beta + alpha)) / t_0) / Float64(t_0 + 2.0)) tmp = 0.0 if (t_1 <= -0.99992) tmp = Float64(Float64(fma(beta, 2.0, 2.0) / alpha) * 0.5); elseif (t_1 <= 5e-84) tmp = 0.5; else tmp = fma(Float64(Float64(beta - alpha) / Float64(2.0 + Float64(beta + alpha))), 0.5, 0.5); end return tmp end
code[alpha_, beta_, i_] := Block[{t$95$0 = N[(N[(i * 2.0), $MachinePrecision] + N[(beta + alpha), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(N[(N[(N[(beta - alpha), $MachinePrecision] * N[(beta + alpha), $MachinePrecision]), $MachinePrecision] / t$95$0), $MachinePrecision] / N[(t$95$0 + 2.0), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$1, -0.99992], N[(N[(N[(beta * 2.0 + 2.0), $MachinePrecision] / alpha), $MachinePrecision] * 0.5), $MachinePrecision], If[LessEqual[t$95$1, 5e-84], 0.5, N[(N[(N[(beta - alpha), $MachinePrecision] / N[(2.0 + N[(beta + alpha), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * 0.5 + 0.5), $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := i \cdot 2 + \left(\beta + \alpha\right)\\
t_1 := \frac{\frac{\left(\beta - \alpha\right) \cdot \left(\beta + \alpha\right)}{t\_0}}{t\_0 + 2}\\
\mathbf{if}\;t\_1 \leq -0.99992:\\
\;\;\;\;\frac{\mathsf{fma}\left(\beta, 2, 2\right)}{\alpha} \cdot 0.5\\
\mathbf{elif}\;t\_1 \leq 5 \cdot 10^{-84}:\\
\;\;\;\;0.5\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(\frac{\beta - \alpha}{2 + \left(\beta + \alpha\right)}, 0.5, 0.5\right)\\
\end{array}
\end{array}
if (/.f64 (/.f64 (*.f64 (+.f64 alpha beta) (-.f64 beta alpha)) (+.f64 (+.f64 alpha beta) (*.f64 #s(literal 2 binary64) i))) (+.f64 (+.f64 (+.f64 alpha beta) (*.f64 #s(literal 2 binary64) i)) #s(literal 2 binary64))) < -0.999920000000000031Initial program 5.0%
Taylor expanded in alpha around inf
*-commutativeN/A
distribute-rgt1-inN/A
metadata-evalN/A
mul0-lftN/A
neg-sub0N/A
mul-1-negN/A
remove-double-negN/A
lower-*.f64N/A
lower-/.f64N/A
+-commutativeN/A
lower-+.f64N/A
+-commutativeN/A
lower-fma.f64N/A
*-commutativeN/A
lower-*.f6493.0
Applied rewrites93.0%
Taylor expanded in i around 0
Applied rewrites78.0%
if -0.999920000000000031 < (/.f64 (/.f64 (*.f64 (+.f64 alpha beta) (-.f64 beta alpha)) (+.f64 (+.f64 alpha beta) (*.f64 #s(literal 2 binary64) i))) (+.f64 (+.f64 (+.f64 alpha beta) (*.f64 #s(literal 2 binary64) i)) #s(literal 2 binary64))) < 5.0000000000000002e-84Initial program 99.8%
Taylor expanded in i around inf
Applied rewrites96.9%
if 5.0000000000000002e-84 < (/.f64 (/.f64 (*.f64 (+.f64 alpha beta) (-.f64 beta alpha)) (+.f64 (+.f64 alpha beta) (*.f64 #s(literal 2 binary64) i))) (+.f64 (+.f64 (+.f64 alpha beta) (*.f64 #s(literal 2 binary64) i)) #s(literal 2 binary64))) Initial program 42.1%
lift-/.f64N/A
clear-numN/A
lower-/.f64N/A
lower-/.f6442.1
lift-+.f64N/A
Applied rewrites99.9%
lift-/.f64N/A
lift-/.f64N/A
associate-/r/N/A
metadata-evalN/A
lift-fma.f64N/A
distribute-rgt-inN/A
metadata-evalN/A
lower-fma.f64N/A
Applied rewrites99.9%
Taylor expanded in i around 0
lower-/.f64N/A
lower--.f64N/A
+-commutativeN/A
lower-+.f64N/A
+-commutativeN/A
lower-+.f6494.9
Applied rewrites94.9%
Final simplification92.3%
(FPCore (alpha beta i)
:precision binary64
(let* ((t_0 (+ (* i 2.0) (+ beta alpha)))
(t_1 (/ (/ (* (- beta alpha) (+ beta alpha)) t_0) (+ t_0 2.0))))
(if (<= t_1 -0.99992)
(* (/ (fma beta 2.0 2.0) alpha) 0.5)
(if (<= t_1 0.002) 0.5 (fma (/ (fma 4.0 i 2.0) beta) -0.5 1.0)))))
double code(double alpha, double beta, double i) {
double t_0 = (i * 2.0) + (beta + alpha);
double t_1 = (((beta - alpha) * (beta + alpha)) / t_0) / (t_0 + 2.0);
double tmp;
if (t_1 <= -0.99992) {
tmp = (fma(beta, 2.0, 2.0) / alpha) * 0.5;
} else if (t_1 <= 0.002) {
tmp = 0.5;
} else {
tmp = fma((fma(4.0, i, 2.0) / beta), -0.5, 1.0);
}
return tmp;
}
function code(alpha, beta, i) t_0 = Float64(Float64(i * 2.0) + Float64(beta + alpha)) t_1 = Float64(Float64(Float64(Float64(beta - alpha) * Float64(beta + alpha)) / t_0) / Float64(t_0 + 2.0)) tmp = 0.0 if (t_1 <= -0.99992) tmp = Float64(Float64(fma(beta, 2.0, 2.0) / alpha) * 0.5); elseif (t_1 <= 0.002) tmp = 0.5; else tmp = fma(Float64(fma(4.0, i, 2.0) / beta), -0.5, 1.0); end return tmp end
code[alpha_, beta_, i_] := Block[{t$95$0 = N[(N[(i * 2.0), $MachinePrecision] + N[(beta + alpha), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(N[(N[(N[(beta - alpha), $MachinePrecision] * N[(beta + alpha), $MachinePrecision]), $MachinePrecision] / t$95$0), $MachinePrecision] / N[(t$95$0 + 2.0), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$1, -0.99992], N[(N[(N[(beta * 2.0 + 2.0), $MachinePrecision] / alpha), $MachinePrecision] * 0.5), $MachinePrecision], If[LessEqual[t$95$1, 0.002], 0.5, N[(N[(N[(4.0 * i + 2.0), $MachinePrecision] / beta), $MachinePrecision] * -0.5 + 1.0), $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := i \cdot 2 + \left(\beta + \alpha\right)\\
t_1 := \frac{\frac{\left(\beta - \alpha\right) \cdot \left(\beta + \alpha\right)}{t\_0}}{t\_0 + 2}\\
\mathbf{if}\;t\_1 \leq -0.99992:\\
\;\;\;\;\frac{\mathsf{fma}\left(\beta, 2, 2\right)}{\alpha} \cdot 0.5\\
\mathbf{elif}\;t\_1 \leq 0.002:\\
\;\;\;\;0.5\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(\frac{\mathsf{fma}\left(4, i, 2\right)}{\beta}, -0.5, 1\right)\\
\end{array}
\end{array}
if (/.f64 (/.f64 (*.f64 (+.f64 alpha beta) (-.f64 beta alpha)) (+.f64 (+.f64 alpha beta) (*.f64 #s(literal 2 binary64) i))) (+.f64 (+.f64 (+.f64 alpha beta) (*.f64 #s(literal 2 binary64) i)) #s(literal 2 binary64))) < -0.999920000000000031Initial program 5.0%
Taylor expanded in alpha around inf
*-commutativeN/A
distribute-rgt1-inN/A
metadata-evalN/A
mul0-lftN/A
neg-sub0N/A
mul-1-negN/A
remove-double-negN/A
lower-*.f64N/A
lower-/.f64N/A
+-commutativeN/A
lower-+.f64N/A
+-commutativeN/A
lower-fma.f64N/A
*-commutativeN/A
lower-*.f6493.0
Applied rewrites93.0%
Taylor expanded in i around 0
Applied rewrites78.0%
if -0.999920000000000031 < (/.f64 (/.f64 (*.f64 (+.f64 alpha beta) (-.f64 beta alpha)) (+.f64 (+.f64 alpha beta) (*.f64 #s(literal 2 binary64) i))) (+.f64 (+.f64 (+.f64 alpha beta) (*.f64 #s(literal 2 binary64) i)) #s(literal 2 binary64))) < 2e-3Initial program 99.9%
Taylor expanded in i around inf
Applied rewrites96.3%
if 2e-3 < (/.f64 (/.f64 (*.f64 (+.f64 alpha beta) (-.f64 beta alpha)) (+.f64 (+.f64 alpha beta) (*.f64 #s(literal 2 binary64) i))) (+.f64 (+.f64 (+.f64 alpha beta) (*.f64 #s(literal 2 binary64) i)) #s(literal 2 binary64))) Initial program 30.0%
Taylor expanded in alpha around 0
*-commutativeN/A
lower-*.f64N/A
+-commutativeN/A
unpow2N/A
times-fracN/A
lower-fma.f64N/A
lower-/.f64N/A
+-commutativeN/A
lower-+.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower-/.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f6499.1
Applied rewrites99.1%
Taylor expanded in beta around inf
Applied rewrites94.1%
Final simplification91.9%
(FPCore (alpha beta i)
:precision binary64
(let* ((t_0 (+ (* i 2.0) (+ beta alpha)))
(t_1 (/ (/ (* (- beta alpha) (+ beta alpha)) t_0) (+ t_0 2.0))))
(if (<= t_1 -0.99992)
(* (/ (fma beta 2.0 2.0) alpha) 0.5)
(if (<= t_1 0.002) 0.5 (fma (/ (* 4.0 i) beta) -0.5 1.0)))))
double code(double alpha, double beta, double i) {
double t_0 = (i * 2.0) + (beta + alpha);
double t_1 = (((beta - alpha) * (beta + alpha)) / t_0) / (t_0 + 2.0);
double tmp;
if (t_1 <= -0.99992) {
tmp = (fma(beta, 2.0, 2.0) / alpha) * 0.5;
} else if (t_1 <= 0.002) {
tmp = 0.5;
} else {
tmp = fma(((4.0 * i) / beta), -0.5, 1.0);
}
return tmp;
}
function code(alpha, beta, i) t_0 = Float64(Float64(i * 2.0) + Float64(beta + alpha)) t_1 = Float64(Float64(Float64(Float64(beta - alpha) * Float64(beta + alpha)) / t_0) / Float64(t_0 + 2.0)) tmp = 0.0 if (t_1 <= -0.99992) tmp = Float64(Float64(fma(beta, 2.0, 2.0) / alpha) * 0.5); elseif (t_1 <= 0.002) tmp = 0.5; else tmp = fma(Float64(Float64(4.0 * i) / beta), -0.5, 1.0); end return tmp end
code[alpha_, beta_, i_] := Block[{t$95$0 = N[(N[(i * 2.0), $MachinePrecision] + N[(beta + alpha), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(N[(N[(N[(beta - alpha), $MachinePrecision] * N[(beta + alpha), $MachinePrecision]), $MachinePrecision] / t$95$0), $MachinePrecision] / N[(t$95$0 + 2.0), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$1, -0.99992], N[(N[(N[(beta * 2.0 + 2.0), $MachinePrecision] / alpha), $MachinePrecision] * 0.5), $MachinePrecision], If[LessEqual[t$95$1, 0.002], 0.5, N[(N[(N[(4.0 * i), $MachinePrecision] / beta), $MachinePrecision] * -0.5 + 1.0), $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := i \cdot 2 + \left(\beta + \alpha\right)\\
t_1 := \frac{\frac{\left(\beta - \alpha\right) \cdot \left(\beta + \alpha\right)}{t\_0}}{t\_0 + 2}\\
\mathbf{if}\;t\_1 \leq -0.99992:\\
\;\;\;\;\frac{\mathsf{fma}\left(\beta, 2, 2\right)}{\alpha} \cdot 0.5\\
\mathbf{elif}\;t\_1 \leq 0.002:\\
\;\;\;\;0.5\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(\frac{4 \cdot i}{\beta}, -0.5, 1\right)\\
\end{array}
\end{array}
if (/.f64 (/.f64 (*.f64 (+.f64 alpha beta) (-.f64 beta alpha)) (+.f64 (+.f64 alpha beta) (*.f64 #s(literal 2 binary64) i))) (+.f64 (+.f64 (+.f64 alpha beta) (*.f64 #s(literal 2 binary64) i)) #s(literal 2 binary64))) < -0.999920000000000031Initial program 5.0%
Taylor expanded in alpha around inf
*-commutativeN/A
distribute-rgt1-inN/A
metadata-evalN/A
mul0-lftN/A
neg-sub0N/A
mul-1-negN/A
remove-double-negN/A
lower-*.f64N/A
lower-/.f64N/A
+-commutativeN/A
lower-+.f64N/A
+-commutativeN/A
lower-fma.f64N/A
*-commutativeN/A
lower-*.f6493.0
Applied rewrites93.0%
Taylor expanded in i around 0
Applied rewrites78.0%
if -0.999920000000000031 < (/.f64 (/.f64 (*.f64 (+.f64 alpha beta) (-.f64 beta alpha)) (+.f64 (+.f64 alpha beta) (*.f64 #s(literal 2 binary64) i))) (+.f64 (+.f64 (+.f64 alpha beta) (*.f64 #s(literal 2 binary64) i)) #s(literal 2 binary64))) < 2e-3Initial program 99.9%
Taylor expanded in i around inf
Applied rewrites96.3%
if 2e-3 < (/.f64 (/.f64 (*.f64 (+.f64 alpha beta) (-.f64 beta alpha)) (+.f64 (+.f64 alpha beta) (*.f64 #s(literal 2 binary64) i))) (+.f64 (+.f64 (+.f64 alpha beta) (*.f64 #s(literal 2 binary64) i)) #s(literal 2 binary64))) Initial program 30.0%
Taylor expanded in alpha around 0
*-commutativeN/A
lower-*.f64N/A
+-commutativeN/A
unpow2N/A
times-fracN/A
lower-fma.f64N/A
lower-/.f64N/A
+-commutativeN/A
lower-+.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower-/.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f6499.1
Applied rewrites99.1%
Taylor expanded in beta around inf
Applied rewrites94.1%
Taylor expanded in i around inf
Applied rewrites94.1%
Final simplification91.9%
(FPCore (alpha beta i)
:precision binary64
(let* ((t_0 (+ (* i 2.0) (+ beta alpha)))
(t_1 (/ (/ (* (- beta alpha) (+ beta alpha)) t_0) (+ t_0 2.0))))
(if (<= t_1 -0.99992)
(* (fma 2.0 beta 2.0) (/ 0.5 alpha))
(if (<= t_1 0.002) 0.5 (fma (/ (* 4.0 i) beta) -0.5 1.0)))))
double code(double alpha, double beta, double i) {
double t_0 = (i * 2.0) + (beta + alpha);
double t_1 = (((beta - alpha) * (beta + alpha)) / t_0) / (t_0 + 2.0);
double tmp;
if (t_1 <= -0.99992) {
tmp = fma(2.0, beta, 2.0) * (0.5 / alpha);
} else if (t_1 <= 0.002) {
tmp = 0.5;
} else {
tmp = fma(((4.0 * i) / beta), -0.5, 1.0);
}
return tmp;
}
function code(alpha, beta, i) t_0 = Float64(Float64(i * 2.0) + Float64(beta + alpha)) t_1 = Float64(Float64(Float64(Float64(beta - alpha) * Float64(beta + alpha)) / t_0) / Float64(t_0 + 2.0)) tmp = 0.0 if (t_1 <= -0.99992) tmp = Float64(fma(2.0, beta, 2.0) * Float64(0.5 / alpha)); elseif (t_1 <= 0.002) tmp = 0.5; else tmp = fma(Float64(Float64(4.0 * i) / beta), -0.5, 1.0); end return tmp end
code[alpha_, beta_, i_] := Block[{t$95$0 = N[(N[(i * 2.0), $MachinePrecision] + N[(beta + alpha), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(N[(N[(N[(beta - alpha), $MachinePrecision] * N[(beta + alpha), $MachinePrecision]), $MachinePrecision] / t$95$0), $MachinePrecision] / N[(t$95$0 + 2.0), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$1, -0.99992], N[(N[(2.0 * beta + 2.0), $MachinePrecision] * N[(0.5 / alpha), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$1, 0.002], 0.5, N[(N[(N[(4.0 * i), $MachinePrecision] / beta), $MachinePrecision] * -0.5 + 1.0), $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := i \cdot 2 + \left(\beta + \alpha\right)\\
t_1 := \frac{\frac{\left(\beta - \alpha\right) \cdot \left(\beta + \alpha\right)}{t\_0}}{t\_0 + 2}\\
\mathbf{if}\;t\_1 \leq -0.99992:\\
\;\;\;\;\mathsf{fma}\left(2, \beta, 2\right) \cdot \frac{0.5}{\alpha}\\
\mathbf{elif}\;t\_1 \leq 0.002:\\
\;\;\;\;0.5\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(\frac{4 \cdot i}{\beta}, -0.5, 1\right)\\
\end{array}
\end{array}
if (/.f64 (/.f64 (*.f64 (+.f64 alpha beta) (-.f64 beta alpha)) (+.f64 (+.f64 alpha beta) (*.f64 #s(literal 2 binary64) i))) (+.f64 (+.f64 (+.f64 alpha beta) (*.f64 #s(literal 2 binary64) i)) #s(literal 2 binary64))) < -0.999920000000000031Initial program 5.0%
Taylor expanded in alpha around inf
*-commutativeN/A
distribute-rgt1-inN/A
metadata-evalN/A
mul0-lftN/A
neg-sub0N/A
mul-1-negN/A
remove-double-negN/A
lower-*.f64N/A
lower-/.f64N/A
+-commutativeN/A
lower-+.f64N/A
+-commutativeN/A
lower-fma.f64N/A
*-commutativeN/A
lower-*.f6493.0
Applied rewrites93.0%
Taylor expanded in i around 0
Applied rewrites78.0%
Applied rewrites78.0%
if -0.999920000000000031 < (/.f64 (/.f64 (*.f64 (+.f64 alpha beta) (-.f64 beta alpha)) (+.f64 (+.f64 alpha beta) (*.f64 #s(literal 2 binary64) i))) (+.f64 (+.f64 (+.f64 alpha beta) (*.f64 #s(literal 2 binary64) i)) #s(literal 2 binary64))) < 2e-3Initial program 99.9%
Taylor expanded in i around inf
Applied rewrites96.3%
if 2e-3 < (/.f64 (/.f64 (*.f64 (+.f64 alpha beta) (-.f64 beta alpha)) (+.f64 (+.f64 alpha beta) (*.f64 #s(literal 2 binary64) i))) (+.f64 (+.f64 (+.f64 alpha beta) (*.f64 #s(literal 2 binary64) i)) #s(literal 2 binary64))) Initial program 30.0%
Taylor expanded in alpha around 0
*-commutativeN/A
lower-*.f64N/A
+-commutativeN/A
unpow2N/A
times-fracN/A
lower-fma.f64N/A
lower-/.f64N/A
+-commutativeN/A
lower-+.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower-/.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f6499.1
Applied rewrites99.1%
Taylor expanded in beta around inf
Applied rewrites94.1%
Taylor expanded in i around inf
Applied rewrites94.1%
Final simplification91.9%
(FPCore (alpha beta i)
:precision binary64
(let* ((t_0 (+ (* i 2.0) (+ beta alpha)))
(t_1 (/ (/ (* (- beta alpha) (+ beta alpha)) t_0) (+ t_0 2.0))))
(if (<= t_1 -0.9999999999999)
(/ beta alpha)
(if (<= t_1 0.002) 0.5 (fma (/ (* 4.0 i) beta) -0.5 1.0)))))
double code(double alpha, double beta, double i) {
double t_0 = (i * 2.0) + (beta + alpha);
double t_1 = (((beta - alpha) * (beta + alpha)) / t_0) / (t_0 + 2.0);
double tmp;
if (t_1 <= -0.9999999999999) {
tmp = beta / alpha;
} else if (t_1 <= 0.002) {
tmp = 0.5;
} else {
tmp = fma(((4.0 * i) / beta), -0.5, 1.0);
}
return tmp;
}
function code(alpha, beta, i) t_0 = Float64(Float64(i * 2.0) + Float64(beta + alpha)) t_1 = Float64(Float64(Float64(Float64(beta - alpha) * Float64(beta + alpha)) / t_0) / Float64(t_0 + 2.0)) tmp = 0.0 if (t_1 <= -0.9999999999999) tmp = Float64(beta / alpha); elseif (t_1 <= 0.002) tmp = 0.5; else tmp = fma(Float64(Float64(4.0 * i) / beta), -0.5, 1.0); end return tmp end
code[alpha_, beta_, i_] := Block[{t$95$0 = N[(N[(i * 2.0), $MachinePrecision] + N[(beta + alpha), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(N[(N[(N[(beta - alpha), $MachinePrecision] * N[(beta + alpha), $MachinePrecision]), $MachinePrecision] / t$95$0), $MachinePrecision] / N[(t$95$0 + 2.0), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$1, -0.9999999999999], N[(beta / alpha), $MachinePrecision], If[LessEqual[t$95$1, 0.002], 0.5, N[(N[(N[(4.0 * i), $MachinePrecision] / beta), $MachinePrecision] * -0.5 + 1.0), $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := i \cdot 2 + \left(\beta + \alpha\right)\\
t_1 := \frac{\frac{\left(\beta - \alpha\right) \cdot \left(\beta + \alpha\right)}{t\_0}}{t\_0 + 2}\\
\mathbf{if}\;t\_1 \leq -0.9999999999999:\\
\;\;\;\;\frac{\beta}{\alpha}\\
\mathbf{elif}\;t\_1 \leq 0.002:\\
\;\;\;\;0.5\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(\frac{4 \cdot i}{\beta}, -0.5, 1\right)\\
\end{array}
\end{array}
if (/.f64 (/.f64 (*.f64 (+.f64 alpha beta) (-.f64 beta alpha)) (+.f64 (+.f64 alpha beta) (*.f64 #s(literal 2 binary64) i))) (+.f64 (+.f64 (+.f64 alpha beta) (*.f64 #s(literal 2 binary64) i)) #s(literal 2 binary64))) < -0.999999999999899969Initial program 2.8%
Taylor expanded in alpha around inf
*-commutativeN/A
distribute-rgt1-inN/A
metadata-evalN/A
mul0-lftN/A
neg-sub0N/A
mul-1-negN/A
remove-double-negN/A
lower-*.f64N/A
lower-/.f64N/A
+-commutativeN/A
lower-+.f64N/A
+-commutativeN/A
lower-fma.f64N/A
*-commutativeN/A
lower-*.f6494.1
Applied rewrites94.1%
Taylor expanded in beta around inf
Applied rewrites28.0%
if -0.999999999999899969 < (/.f64 (/.f64 (*.f64 (+.f64 alpha beta) (-.f64 beta alpha)) (+.f64 (+.f64 alpha beta) (*.f64 #s(literal 2 binary64) i))) (+.f64 (+.f64 (+.f64 alpha beta) (*.f64 #s(literal 2 binary64) i)) #s(literal 2 binary64))) < 2e-3Initial program 99.3%
Taylor expanded in i around inf
Applied rewrites95.0%
if 2e-3 < (/.f64 (/.f64 (*.f64 (+.f64 alpha beta) (-.f64 beta alpha)) (+.f64 (+.f64 alpha beta) (*.f64 #s(literal 2 binary64) i))) (+.f64 (+.f64 (+.f64 alpha beta) (*.f64 #s(literal 2 binary64) i)) #s(literal 2 binary64))) Initial program 30.0%
Taylor expanded in alpha around 0
*-commutativeN/A
lower-*.f64N/A
+-commutativeN/A
unpow2N/A
times-fracN/A
lower-fma.f64N/A
lower-/.f64N/A
+-commutativeN/A
lower-+.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower-/.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f6499.1
Applied rewrites99.1%
Taylor expanded in beta around inf
Applied rewrites94.1%
Taylor expanded in i around inf
Applied rewrites94.1%
Final simplification81.4%
(FPCore (alpha beta i)
:precision binary64
(let* ((t_0 (+ (* i 2.0) (+ beta alpha)))
(t_1 (/ (/ (* (- beta alpha) (+ beta alpha)) t_0) (+ t_0 2.0))))
(if (<= t_1 -0.9999999999999)
(/ beta alpha)
(if (<= t_1 5e-84) 0.5 (fma (/ beta (+ 2.0 beta)) 0.5 0.5)))))
double code(double alpha, double beta, double i) {
double t_0 = (i * 2.0) + (beta + alpha);
double t_1 = (((beta - alpha) * (beta + alpha)) / t_0) / (t_0 + 2.0);
double tmp;
if (t_1 <= -0.9999999999999) {
tmp = beta / alpha;
} else if (t_1 <= 5e-84) {
tmp = 0.5;
} else {
tmp = fma((beta / (2.0 + beta)), 0.5, 0.5);
}
return tmp;
}
function code(alpha, beta, i) t_0 = Float64(Float64(i * 2.0) + Float64(beta + alpha)) t_1 = Float64(Float64(Float64(Float64(beta - alpha) * Float64(beta + alpha)) / t_0) / Float64(t_0 + 2.0)) tmp = 0.0 if (t_1 <= -0.9999999999999) tmp = Float64(beta / alpha); elseif (t_1 <= 5e-84) tmp = 0.5; else tmp = fma(Float64(beta / Float64(2.0 + beta)), 0.5, 0.5); end return tmp end
code[alpha_, beta_, i_] := Block[{t$95$0 = N[(N[(i * 2.0), $MachinePrecision] + N[(beta + alpha), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(N[(N[(N[(beta - alpha), $MachinePrecision] * N[(beta + alpha), $MachinePrecision]), $MachinePrecision] / t$95$0), $MachinePrecision] / N[(t$95$0 + 2.0), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$1, -0.9999999999999], N[(beta / alpha), $MachinePrecision], If[LessEqual[t$95$1, 5e-84], 0.5, N[(N[(beta / N[(2.0 + beta), $MachinePrecision]), $MachinePrecision] * 0.5 + 0.5), $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := i \cdot 2 + \left(\beta + \alpha\right)\\
t_1 := \frac{\frac{\left(\beta - \alpha\right) \cdot \left(\beta + \alpha\right)}{t\_0}}{t\_0 + 2}\\
\mathbf{if}\;t\_1 \leq -0.9999999999999:\\
\;\;\;\;\frac{\beta}{\alpha}\\
\mathbf{elif}\;t\_1 \leq 5 \cdot 10^{-84}:\\
\;\;\;\;0.5\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(\frac{\beta}{2 + \beta}, 0.5, 0.5\right)\\
\end{array}
\end{array}
if (/.f64 (/.f64 (*.f64 (+.f64 alpha beta) (-.f64 beta alpha)) (+.f64 (+.f64 alpha beta) (*.f64 #s(literal 2 binary64) i))) (+.f64 (+.f64 (+.f64 alpha beta) (*.f64 #s(literal 2 binary64) i)) #s(literal 2 binary64))) < -0.999999999999899969Initial program 2.8%
Taylor expanded in alpha around inf
*-commutativeN/A
distribute-rgt1-inN/A
metadata-evalN/A
mul0-lftN/A
neg-sub0N/A
mul-1-negN/A
remove-double-negN/A
lower-*.f64N/A
lower-/.f64N/A
+-commutativeN/A
lower-+.f64N/A
+-commutativeN/A
lower-fma.f64N/A
*-commutativeN/A
lower-*.f6494.1
Applied rewrites94.1%
Taylor expanded in beta around inf
Applied rewrites28.0%
if -0.999999999999899969 < (/.f64 (/.f64 (*.f64 (+.f64 alpha beta) (-.f64 beta alpha)) (+.f64 (+.f64 alpha beta) (*.f64 #s(literal 2 binary64) i))) (+.f64 (+.f64 (+.f64 alpha beta) (*.f64 #s(literal 2 binary64) i)) #s(literal 2 binary64))) < 5.0000000000000002e-84Initial program 99.2%
Taylor expanded in i around inf
Applied rewrites95.4%
if 5.0000000000000002e-84 < (/.f64 (/.f64 (*.f64 (+.f64 alpha beta) (-.f64 beta alpha)) (+.f64 (+.f64 alpha beta) (*.f64 #s(literal 2 binary64) i))) (+.f64 (+.f64 (+.f64 alpha beta) (*.f64 #s(literal 2 binary64) i)) #s(literal 2 binary64))) Initial program 42.1%
lift-/.f64N/A
clear-numN/A
lower-/.f64N/A
lower-/.f6442.1
lift-+.f64N/A
Applied rewrites99.9%
lift-/.f64N/A
lift-/.f64N/A
associate-/r/N/A
metadata-evalN/A
lift-fma.f64N/A
distribute-rgt-inN/A
metadata-evalN/A
lower-fma.f64N/A
Applied rewrites99.9%
Taylor expanded in alpha around 0
unpow2N/A
times-fracN/A
lower-*.f64N/A
lower-/.f64N/A
+-commutativeN/A
lower-+.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower-/.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f6497.8
Applied rewrites97.8%
Taylor expanded in i around 0
Applied rewrites92.7%
Final simplification81.1%
(FPCore (alpha beta i)
:precision binary64
(let* ((t_0 (+ (* i 2.0) (+ beta alpha)))
(t_1 (/ (/ (* (- beta alpha) (+ beta alpha)) t_0) (+ t_0 2.0))))
(if (<= t_1 -0.9999999999999)
(/ beta alpha)
(if (<= t_1 0.002) 0.5 (- 1.0 (/ 1.0 beta))))))
double code(double alpha, double beta, double i) {
double t_0 = (i * 2.0) + (beta + alpha);
double t_1 = (((beta - alpha) * (beta + alpha)) / t_0) / (t_0 + 2.0);
double tmp;
if (t_1 <= -0.9999999999999) {
tmp = beta / alpha;
} else if (t_1 <= 0.002) {
tmp = 0.5;
} else {
tmp = 1.0 - (1.0 / beta);
}
return tmp;
}
real(8) function code(alpha, beta, i)
real(8), intent (in) :: alpha
real(8), intent (in) :: beta
real(8), intent (in) :: i
real(8) :: t_0
real(8) :: t_1
real(8) :: tmp
t_0 = (i * 2.0d0) + (beta + alpha)
t_1 = (((beta - alpha) * (beta + alpha)) / t_0) / (t_0 + 2.0d0)
if (t_1 <= (-0.9999999999999d0)) then
tmp = beta / alpha
else if (t_1 <= 0.002d0) then
tmp = 0.5d0
else
tmp = 1.0d0 - (1.0d0 / beta)
end if
code = tmp
end function
public static double code(double alpha, double beta, double i) {
double t_0 = (i * 2.0) + (beta + alpha);
double t_1 = (((beta - alpha) * (beta + alpha)) / t_0) / (t_0 + 2.0);
double tmp;
if (t_1 <= -0.9999999999999) {
tmp = beta / alpha;
} else if (t_1 <= 0.002) {
tmp = 0.5;
} else {
tmp = 1.0 - (1.0 / beta);
}
return tmp;
}
def code(alpha, beta, i): t_0 = (i * 2.0) + (beta + alpha) t_1 = (((beta - alpha) * (beta + alpha)) / t_0) / (t_0 + 2.0) tmp = 0 if t_1 <= -0.9999999999999: tmp = beta / alpha elif t_1 <= 0.002: tmp = 0.5 else: tmp = 1.0 - (1.0 / beta) return tmp
function code(alpha, beta, i) t_0 = Float64(Float64(i * 2.0) + Float64(beta + alpha)) t_1 = Float64(Float64(Float64(Float64(beta - alpha) * Float64(beta + alpha)) / t_0) / Float64(t_0 + 2.0)) tmp = 0.0 if (t_1 <= -0.9999999999999) tmp = Float64(beta / alpha); elseif (t_1 <= 0.002) tmp = 0.5; else tmp = Float64(1.0 - Float64(1.0 / beta)); end return tmp end
function tmp_2 = code(alpha, beta, i) t_0 = (i * 2.0) + (beta + alpha); t_1 = (((beta - alpha) * (beta + alpha)) / t_0) / (t_0 + 2.0); tmp = 0.0; if (t_1 <= -0.9999999999999) tmp = beta / alpha; elseif (t_1 <= 0.002) tmp = 0.5; else tmp = 1.0 - (1.0 / beta); end tmp_2 = tmp; end
code[alpha_, beta_, i_] := Block[{t$95$0 = N[(N[(i * 2.0), $MachinePrecision] + N[(beta + alpha), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(N[(N[(N[(beta - alpha), $MachinePrecision] * N[(beta + alpha), $MachinePrecision]), $MachinePrecision] / t$95$0), $MachinePrecision] / N[(t$95$0 + 2.0), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$1, -0.9999999999999], N[(beta / alpha), $MachinePrecision], If[LessEqual[t$95$1, 0.002], 0.5, N[(1.0 - N[(1.0 / beta), $MachinePrecision]), $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := i \cdot 2 + \left(\beta + \alpha\right)\\
t_1 := \frac{\frac{\left(\beta - \alpha\right) \cdot \left(\beta + \alpha\right)}{t\_0}}{t\_0 + 2}\\
\mathbf{if}\;t\_1 \leq -0.9999999999999:\\
\;\;\;\;\frac{\beta}{\alpha}\\
\mathbf{elif}\;t\_1 \leq 0.002:\\
\;\;\;\;0.5\\
\mathbf{else}:\\
\;\;\;\;1 - \frac{1}{\beta}\\
\end{array}
\end{array}
if (/.f64 (/.f64 (*.f64 (+.f64 alpha beta) (-.f64 beta alpha)) (+.f64 (+.f64 alpha beta) (*.f64 #s(literal 2 binary64) i))) (+.f64 (+.f64 (+.f64 alpha beta) (*.f64 #s(literal 2 binary64) i)) #s(literal 2 binary64))) < -0.999999999999899969Initial program 2.8%
Taylor expanded in alpha around inf
*-commutativeN/A
distribute-rgt1-inN/A
metadata-evalN/A
mul0-lftN/A
neg-sub0N/A
mul-1-negN/A
remove-double-negN/A
lower-*.f64N/A
lower-/.f64N/A
+-commutativeN/A
lower-+.f64N/A
+-commutativeN/A
lower-fma.f64N/A
*-commutativeN/A
lower-*.f6494.1
Applied rewrites94.1%
Taylor expanded in beta around inf
Applied rewrites28.0%
if -0.999999999999899969 < (/.f64 (/.f64 (*.f64 (+.f64 alpha beta) (-.f64 beta alpha)) (+.f64 (+.f64 alpha beta) (*.f64 #s(literal 2 binary64) i))) (+.f64 (+.f64 (+.f64 alpha beta) (*.f64 #s(literal 2 binary64) i)) #s(literal 2 binary64))) < 2e-3Initial program 99.3%
Taylor expanded in i around inf
Applied rewrites95.0%
if 2e-3 < (/.f64 (/.f64 (*.f64 (+.f64 alpha beta) (-.f64 beta alpha)) (+.f64 (+.f64 alpha beta) (*.f64 #s(literal 2 binary64) i))) (+.f64 (+.f64 (+.f64 alpha beta) (*.f64 #s(literal 2 binary64) i)) #s(literal 2 binary64))) Initial program 30.0%
Taylor expanded in alpha around 0
*-commutativeN/A
lower-*.f64N/A
+-commutativeN/A
unpow2N/A
times-fracN/A
lower-fma.f64N/A
lower-/.f64N/A
+-commutativeN/A
lower-+.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower-/.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f6499.1
Applied rewrites99.1%
Taylor expanded in beta around inf
Applied rewrites94.1%
Taylor expanded in i around 0
Applied rewrites92.9%
Final simplification81.1%
(FPCore (alpha beta i)
:precision binary64
(let* ((t_0 (+ (* i 2.0) (+ beta alpha)))
(t_1 (/ (/ (* (- beta alpha) (+ beta alpha)) t_0) (+ t_0 2.0))))
(if (<= t_1 -0.9999999999999) (/ beta alpha) (if (<= t_1 0.002) 0.5 1.0))))
double code(double alpha, double beta, double i) {
double t_0 = (i * 2.0) + (beta + alpha);
double t_1 = (((beta - alpha) * (beta + alpha)) / t_0) / (t_0 + 2.0);
double tmp;
if (t_1 <= -0.9999999999999) {
tmp = beta / alpha;
} else if (t_1 <= 0.002) {
tmp = 0.5;
} else {
tmp = 1.0;
}
return tmp;
}
real(8) function code(alpha, beta, i)
real(8), intent (in) :: alpha
real(8), intent (in) :: beta
real(8), intent (in) :: i
real(8) :: t_0
real(8) :: t_1
real(8) :: tmp
t_0 = (i * 2.0d0) + (beta + alpha)
t_1 = (((beta - alpha) * (beta + alpha)) / t_0) / (t_0 + 2.0d0)
if (t_1 <= (-0.9999999999999d0)) then
tmp = beta / alpha
else if (t_1 <= 0.002d0) then
tmp = 0.5d0
else
tmp = 1.0d0
end if
code = tmp
end function
public static double code(double alpha, double beta, double i) {
double t_0 = (i * 2.0) + (beta + alpha);
double t_1 = (((beta - alpha) * (beta + alpha)) / t_0) / (t_0 + 2.0);
double tmp;
if (t_1 <= -0.9999999999999) {
tmp = beta / alpha;
} else if (t_1 <= 0.002) {
tmp = 0.5;
} else {
tmp = 1.0;
}
return tmp;
}
def code(alpha, beta, i): t_0 = (i * 2.0) + (beta + alpha) t_1 = (((beta - alpha) * (beta + alpha)) / t_0) / (t_0 + 2.0) tmp = 0 if t_1 <= -0.9999999999999: tmp = beta / alpha elif t_1 <= 0.002: tmp = 0.5 else: tmp = 1.0 return tmp
function code(alpha, beta, i) t_0 = Float64(Float64(i * 2.0) + Float64(beta + alpha)) t_1 = Float64(Float64(Float64(Float64(beta - alpha) * Float64(beta + alpha)) / t_0) / Float64(t_0 + 2.0)) tmp = 0.0 if (t_1 <= -0.9999999999999) tmp = Float64(beta / alpha); elseif (t_1 <= 0.002) tmp = 0.5; else tmp = 1.0; end return tmp end
function tmp_2 = code(alpha, beta, i) t_0 = (i * 2.0) + (beta + alpha); t_1 = (((beta - alpha) * (beta + alpha)) / t_0) / (t_0 + 2.0); tmp = 0.0; if (t_1 <= -0.9999999999999) tmp = beta / alpha; elseif (t_1 <= 0.002) tmp = 0.5; else tmp = 1.0; end tmp_2 = tmp; end
code[alpha_, beta_, i_] := Block[{t$95$0 = N[(N[(i * 2.0), $MachinePrecision] + N[(beta + alpha), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(N[(N[(N[(beta - alpha), $MachinePrecision] * N[(beta + alpha), $MachinePrecision]), $MachinePrecision] / t$95$0), $MachinePrecision] / N[(t$95$0 + 2.0), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$1, -0.9999999999999], N[(beta / alpha), $MachinePrecision], If[LessEqual[t$95$1, 0.002], 0.5, 1.0]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := i \cdot 2 + \left(\beta + \alpha\right)\\
t_1 := \frac{\frac{\left(\beta - \alpha\right) \cdot \left(\beta + \alpha\right)}{t\_0}}{t\_0 + 2}\\
\mathbf{if}\;t\_1 \leq -0.9999999999999:\\
\;\;\;\;\frac{\beta}{\alpha}\\
\mathbf{elif}\;t\_1 \leq 0.002:\\
\;\;\;\;0.5\\
\mathbf{else}:\\
\;\;\;\;1\\
\end{array}
\end{array}
if (/.f64 (/.f64 (*.f64 (+.f64 alpha beta) (-.f64 beta alpha)) (+.f64 (+.f64 alpha beta) (*.f64 #s(literal 2 binary64) i))) (+.f64 (+.f64 (+.f64 alpha beta) (*.f64 #s(literal 2 binary64) i)) #s(literal 2 binary64))) < -0.999999999999899969Initial program 2.8%
Taylor expanded in alpha around inf
*-commutativeN/A
distribute-rgt1-inN/A
metadata-evalN/A
mul0-lftN/A
neg-sub0N/A
mul-1-negN/A
remove-double-negN/A
lower-*.f64N/A
lower-/.f64N/A
+-commutativeN/A
lower-+.f64N/A
+-commutativeN/A
lower-fma.f64N/A
*-commutativeN/A
lower-*.f6494.1
Applied rewrites94.1%
Taylor expanded in beta around inf
Applied rewrites28.0%
if -0.999999999999899969 < (/.f64 (/.f64 (*.f64 (+.f64 alpha beta) (-.f64 beta alpha)) (+.f64 (+.f64 alpha beta) (*.f64 #s(literal 2 binary64) i))) (+.f64 (+.f64 (+.f64 alpha beta) (*.f64 #s(literal 2 binary64) i)) #s(literal 2 binary64))) < 2e-3Initial program 99.3%
Taylor expanded in i around inf
Applied rewrites95.0%
if 2e-3 < (/.f64 (/.f64 (*.f64 (+.f64 alpha beta) (-.f64 beta alpha)) (+.f64 (+.f64 alpha beta) (*.f64 #s(literal 2 binary64) i))) (+.f64 (+.f64 (+.f64 alpha beta) (*.f64 #s(literal 2 binary64) i)) #s(literal 2 binary64))) Initial program 30.0%
Taylor expanded in beta around inf
Applied rewrites92.9%
Final simplification81.1%
(FPCore (alpha beta i)
:precision binary64
(let* ((t_0 (+ (* i 2.0) (+ beta alpha))))
(if (<= (/ (/ (* (- beta alpha) (+ beta alpha)) t_0) (+ t_0 2.0)) -0.5)
(fma (/ (fma beta 2.0 2.0) alpha) 0.5 (* (/ i alpha) 2.0))
(fma
(* (/ beta (+ (fma i 2.0 beta) 2.0)) (/ beta (fma i 2.0 beta)))
0.5
0.5))))
double code(double alpha, double beta, double i) {
double t_0 = (i * 2.0) + (beta + alpha);
double tmp;
if (((((beta - alpha) * (beta + alpha)) / t_0) / (t_0 + 2.0)) <= -0.5) {
tmp = fma((fma(beta, 2.0, 2.0) / alpha), 0.5, ((i / alpha) * 2.0));
} else {
tmp = fma(((beta / (fma(i, 2.0, beta) + 2.0)) * (beta / fma(i, 2.0, beta))), 0.5, 0.5);
}
return tmp;
}
function code(alpha, beta, i) t_0 = Float64(Float64(i * 2.0) + Float64(beta + alpha)) tmp = 0.0 if (Float64(Float64(Float64(Float64(beta - alpha) * Float64(beta + alpha)) / t_0) / Float64(t_0 + 2.0)) <= -0.5) tmp = fma(Float64(fma(beta, 2.0, 2.0) / alpha), 0.5, Float64(Float64(i / alpha) * 2.0)); else tmp = fma(Float64(Float64(beta / Float64(fma(i, 2.0, beta) + 2.0)) * Float64(beta / fma(i, 2.0, beta))), 0.5, 0.5); end return tmp end
code[alpha_, beta_, i_] := Block[{t$95$0 = N[(N[(i * 2.0), $MachinePrecision] + N[(beta + alpha), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[N[(N[(N[(N[(beta - alpha), $MachinePrecision] * N[(beta + alpha), $MachinePrecision]), $MachinePrecision] / t$95$0), $MachinePrecision] / N[(t$95$0 + 2.0), $MachinePrecision]), $MachinePrecision], -0.5], N[(N[(N[(beta * 2.0 + 2.0), $MachinePrecision] / alpha), $MachinePrecision] * 0.5 + N[(N[(i / alpha), $MachinePrecision] * 2.0), $MachinePrecision]), $MachinePrecision], N[(N[(N[(beta / N[(N[(i * 2.0 + beta), $MachinePrecision] + 2.0), $MachinePrecision]), $MachinePrecision] * N[(beta / N[(i * 2.0 + beta), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * 0.5 + 0.5), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := i \cdot 2 + \left(\beta + \alpha\right)\\
\mathbf{if}\;\frac{\frac{\left(\beta - \alpha\right) \cdot \left(\beta + \alpha\right)}{t\_0}}{t\_0 + 2} \leq -0.5:\\
\;\;\;\;\mathsf{fma}\left(\frac{\mathsf{fma}\left(\beta, 2, 2\right)}{\alpha}, 0.5, \frac{i}{\alpha} \cdot 2\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(\frac{\beta}{\mathsf{fma}\left(i, 2, \beta\right) + 2} \cdot \frac{\beta}{\mathsf{fma}\left(i, 2, \beta\right)}, 0.5, 0.5\right)\\
\end{array}
\end{array}
if (/.f64 (/.f64 (*.f64 (+.f64 alpha beta) (-.f64 beta alpha)) (+.f64 (+.f64 alpha beta) (*.f64 #s(literal 2 binary64) i))) (+.f64 (+.f64 (+.f64 alpha beta) (*.f64 #s(literal 2 binary64) i)) #s(literal 2 binary64))) < -0.5Initial program 6.4%
Taylor expanded in alpha around inf
*-commutativeN/A
distribute-rgt1-inN/A
metadata-evalN/A
mul0-lftN/A
neg-sub0N/A
mul-1-negN/A
remove-double-negN/A
lower-*.f64N/A
lower-/.f64N/A
+-commutativeN/A
lower-+.f64N/A
+-commutativeN/A
lower-fma.f64N/A
*-commutativeN/A
lower-*.f6492.0
Applied rewrites92.0%
Taylor expanded in i around 0
Applied rewrites92.0%
if -0.5 < (/.f64 (/.f64 (*.f64 (+.f64 alpha beta) (-.f64 beta alpha)) (+.f64 (+.f64 alpha beta) (*.f64 #s(literal 2 binary64) i))) (+.f64 (+.f64 (+.f64 alpha beta) (*.f64 #s(literal 2 binary64) i)) #s(literal 2 binary64))) Initial program 75.1%
Taylor expanded in alpha around 0
*-commutativeN/A
lower-*.f64N/A
+-commutativeN/A
unpow2N/A
times-fracN/A
lower-fma.f64N/A
lower-/.f64N/A
+-commutativeN/A
lower-+.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower-/.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f6498.0
Applied rewrites98.0%
Applied rewrites98.0%
Final simplification96.7%
(FPCore (alpha beta i)
:precision binary64
(let* ((t_0 (+ (* i 2.0) (+ beta alpha))))
(if (<= (/ (/ (* (- beta alpha) (+ beta alpha)) t_0) (+ t_0 2.0)) -0.5)
(fma (/ (fma beta 2.0 2.0) alpha) 0.5 (* (/ i alpha) 2.0))
(*
(fma (/ beta (+ (fma i 2.0 beta) 2.0)) (/ beta (fma i 2.0 beta)) 1.0)
0.5))))
double code(double alpha, double beta, double i) {
double t_0 = (i * 2.0) + (beta + alpha);
double tmp;
if (((((beta - alpha) * (beta + alpha)) / t_0) / (t_0 + 2.0)) <= -0.5) {
tmp = fma((fma(beta, 2.0, 2.0) / alpha), 0.5, ((i / alpha) * 2.0));
} else {
tmp = fma((beta / (fma(i, 2.0, beta) + 2.0)), (beta / fma(i, 2.0, beta)), 1.0) * 0.5;
}
return tmp;
}
function code(alpha, beta, i) t_0 = Float64(Float64(i * 2.0) + Float64(beta + alpha)) tmp = 0.0 if (Float64(Float64(Float64(Float64(beta - alpha) * Float64(beta + alpha)) / t_0) / Float64(t_0 + 2.0)) <= -0.5) tmp = fma(Float64(fma(beta, 2.0, 2.0) / alpha), 0.5, Float64(Float64(i / alpha) * 2.0)); else tmp = Float64(fma(Float64(beta / Float64(fma(i, 2.0, beta) + 2.0)), Float64(beta / fma(i, 2.0, beta)), 1.0) * 0.5); end return tmp end
code[alpha_, beta_, i_] := Block[{t$95$0 = N[(N[(i * 2.0), $MachinePrecision] + N[(beta + alpha), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[N[(N[(N[(N[(beta - alpha), $MachinePrecision] * N[(beta + alpha), $MachinePrecision]), $MachinePrecision] / t$95$0), $MachinePrecision] / N[(t$95$0 + 2.0), $MachinePrecision]), $MachinePrecision], -0.5], N[(N[(N[(beta * 2.0 + 2.0), $MachinePrecision] / alpha), $MachinePrecision] * 0.5 + N[(N[(i / alpha), $MachinePrecision] * 2.0), $MachinePrecision]), $MachinePrecision], N[(N[(N[(beta / N[(N[(i * 2.0 + beta), $MachinePrecision] + 2.0), $MachinePrecision]), $MachinePrecision] * N[(beta / N[(i * 2.0 + beta), $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision] * 0.5), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := i \cdot 2 + \left(\beta + \alpha\right)\\
\mathbf{if}\;\frac{\frac{\left(\beta - \alpha\right) \cdot \left(\beta + \alpha\right)}{t\_0}}{t\_0 + 2} \leq -0.5:\\
\;\;\;\;\mathsf{fma}\left(\frac{\mathsf{fma}\left(\beta, 2, 2\right)}{\alpha}, 0.5, \frac{i}{\alpha} \cdot 2\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(\frac{\beta}{\mathsf{fma}\left(i, 2, \beta\right) + 2}, \frac{\beta}{\mathsf{fma}\left(i, 2, \beta\right)}, 1\right) \cdot 0.5\\
\end{array}
\end{array}
if (/.f64 (/.f64 (*.f64 (+.f64 alpha beta) (-.f64 beta alpha)) (+.f64 (+.f64 alpha beta) (*.f64 #s(literal 2 binary64) i))) (+.f64 (+.f64 (+.f64 alpha beta) (*.f64 #s(literal 2 binary64) i)) #s(literal 2 binary64))) < -0.5Initial program 6.4%
Taylor expanded in alpha around inf
*-commutativeN/A
distribute-rgt1-inN/A
metadata-evalN/A
mul0-lftN/A
neg-sub0N/A
mul-1-negN/A
remove-double-negN/A
lower-*.f64N/A
lower-/.f64N/A
+-commutativeN/A
lower-+.f64N/A
+-commutativeN/A
lower-fma.f64N/A
*-commutativeN/A
lower-*.f6492.0
Applied rewrites92.0%
Taylor expanded in i around 0
Applied rewrites92.0%
if -0.5 < (/.f64 (/.f64 (*.f64 (+.f64 alpha beta) (-.f64 beta alpha)) (+.f64 (+.f64 alpha beta) (*.f64 #s(literal 2 binary64) i))) (+.f64 (+.f64 (+.f64 alpha beta) (*.f64 #s(literal 2 binary64) i)) #s(literal 2 binary64))) Initial program 75.1%
Taylor expanded in alpha around 0
*-commutativeN/A
lower-*.f64N/A
+-commutativeN/A
unpow2N/A
times-fracN/A
lower-fma.f64N/A
lower-/.f64N/A
+-commutativeN/A
lower-+.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower-/.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f6498.0
Applied rewrites98.0%
Final simplification96.7%
(FPCore (alpha beta i)
:precision binary64
(let* ((t_0 (+ (* i 2.0) (+ beta alpha))))
(if (<= (/ (/ (* (- beta alpha) (+ beta alpha)) t_0) (+ t_0 2.0)) 0.5)
0.5
1.0)))
double code(double alpha, double beta, double i) {
double t_0 = (i * 2.0) + (beta + alpha);
double tmp;
if (((((beta - alpha) * (beta + alpha)) / t_0) / (t_0 + 2.0)) <= 0.5) {
tmp = 0.5;
} else {
tmp = 1.0;
}
return tmp;
}
real(8) function code(alpha, beta, i)
real(8), intent (in) :: alpha
real(8), intent (in) :: beta
real(8), intent (in) :: i
real(8) :: t_0
real(8) :: tmp
t_0 = (i * 2.0d0) + (beta + alpha)
if (((((beta - alpha) * (beta + alpha)) / t_0) / (t_0 + 2.0d0)) <= 0.5d0) then
tmp = 0.5d0
else
tmp = 1.0d0
end if
code = tmp
end function
public static double code(double alpha, double beta, double i) {
double t_0 = (i * 2.0) + (beta + alpha);
double tmp;
if (((((beta - alpha) * (beta + alpha)) / t_0) / (t_0 + 2.0)) <= 0.5) {
tmp = 0.5;
} else {
tmp = 1.0;
}
return tmp;
}
def code(alpha, beta, i): t_0 = (i * 2.0) + (beta + alpha) tmp = 0 if ((((beta - alpha) * (beta + alpha)) / t_0) / (t_0 + 2.0)) <= 0.5: tmp = 0.5 else: tmp = 1.0 return tmp
function code(alpha, beta, i) t_0 = Float64(Float64(i * 2.0) + Float64(beta + alpha)) tmp = 0.0 if (Float64(Float64(Float64(Float64(beta - alpha) * Float64(beta + alpha)) / t_0) / Float64(t_0 + 2.0)) <= 0.5) tmp = 0.5; else tmp = 1.0; end return tmp end
function tmp_2 = code(alpha, beta, i) t_0 = (i * 2.0) + (beta + alpha); tmp = 0.0; if (((((beta - alpha) * (beta + alpha)) / t_0) / (t_0 + 2.0)) <= 0.5) tmp = 0.5; else tmp = 1.0; end tmp_2 = tmp; end
code[alpha_, beta_, i_] := Block[{t$95$0 = N[(N[(i * 2.0), $MachinePrecision] + N[(beta + alpha), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[N[(N[(N[(N[(beta - alpha), $MachinePrecision] * N[(beta + alpha), $MachinePrecision]), $MachinePrecision] / t$95$0), $MachinePrecision] / N[(t$95$0 + 2.0), $MachinePrecision]), $MachinePrecision], 0.5], 0.5, 1.0]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := i \cdot 2 + \left(\beta + \alpha\right)\\
\mathbf{if}\;\frac{\frac{\left(\beta - \alpha\right) \cdot \left(\beta + \alpha\right)}{t\_0}}{t\_0 + 2} \leq 0.5:\\
\;\;\;\;0.5\\
\mathbf{else}:\\
\;\;\;\;1\\
\end{array}
\end{array}
if (/.f64 (/.f64 (*.f64 (+.f64 alpha beta) (-.f64 beta alpha)) (+.f64 (+.f64 alpha beta) (*.f64 #s(literal 2 binary64) i))) (+.f64 (+.f64 (+.f64 alpha beta) (*.f64 #s(literal 2 binary64) i)) #s(literal 2 binary64))) < 0.5Initial program 72.5%
Taylor expanded in i around inf
Applied rewrites71.7%
if 0.5 < (/.f64 (/.f64 (*.f64 (+.f64 alpha beta) (-.f64 beta alpha)) (+.f64 (+.f64 alpha beta) (*.f64 #s(literal 2 binary64) i))) (+.f64 (+.f64 (+.f64 alpha beta) (*.f64 #s(literal 2 binary64) i)) #s(literal 2 binary64))) Initial program 30.0%
Taylor expanded in beta around inf
Applied rewrites92.9%
Final simplification77.7%
(FPCore (alpha beta i) :precision binary64 0.5)
double code(double alpha, double beta, double i) {
return 0.5;
}
real(8) function code(alpha, beta, i)
real(8), intent (in) :: alpha
real(8), intent (in) :: beta
real(8), intent (in) :: i
code = 0.5d0
end function
public static double code(double alpha, double beta, double i) {
return 0.5;
}
def code(alpha, beta, i): return 0.5
function code(alpha, beta, i) return 0.5 end
function tmp = code(alpha, beta, i) tmp = 0.5; end
code[alpha_, beta_, i_] := 0.5
\begin{array}{l}
\\
0.5
\end{array}
Initial program 60.6%
Taylor expanded in i around inf
Applied rewrites57.8%
herbie shell --seed 2024248
(FPCore (alpha beta i)
:name "Octave 3.8, jcobi/2"
:precision binary64
:pre (and (and (> alpha -1.0) (> beta -1.0)) (> i 0.0))
(/ (+ (/ (/ (* (+ alpha beta) (- beta alpha)) (+ (+ alpha beta) (* 2.0 i))) (+ (+ (+ alpha beta) (* 2.0 i)) 2.0)) 1.0) 2.0))