
(FPCore (x y z t a b c)
:precision binary64
(/
x
(+
x
(*
y
(exp
(*
2.0
(-
(/ (* z (sqrt (+ t a))) t)
(* (- b c) (- (+ a (/ 5.0 6.0)) (/ 2.0 (* t 3.0)))))))))))
double code(double x, double y, double z, double t, double a, double b, double c) {
return x / (x + (y * exp((2.0 * (((z * sqrt((t + a))) / t) - ((b - c) * ((a + (5.0 / 6.0)) - (2.0 / (t * 3.0)))))))));
}
real(8) function code(x, y, z, t, a, b, c)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8), intent (in) :: c
code = x / (x + (y * exp((2.0d0 * (((z * sqrt((t + a))) / t) - ((b - c) * ((a + (5.0d0 / 6.0d0)) - (2.0d0 / (t * 3.0d0)))))))))
end function
public static double code(double x, double y, double z, double t, double a, double b, double c) {
return x / (x + (y * Math.exp((2.0 * (((z * Math.sqrt((t + a))) / t) - ((b - c) * ((a + (5.0 / 6.0)) - (2.0 / (t * 3.0)))))))));
}
def code(x, y, z, t, a, b, c): return x / (x + (y * math.exp((2.0 * (((z * math.sqrt((t + a))) / t) - ((b - c) * ((a + (5.0 / 6.0)) - (2.0 / (t * 3.0)))))))))
function code(x, y, z, t, a, b, c) return Float64(x / Float64(x + Float64(y * exp(Float64(2.0 * Float64(Float64(Float64(z * sqrt(Float64(t + a))) / t) - Float64(Float64(b - c) * Float64(Float64(a + Float64(5.0 / 6.0)) - Float64(2.0 / Float64(t * 3.0)))))))))) end
function tmp = code(x, y, z, t, a, b, c) tmp = x / (x + (y * exp((2.0 * (((z * sqrt((t + a))) / t) - ((b - c) * ((a + (5.0 / 6.0)) - (2.0 / (t * 3.0))))))))); end
code[x_, y_, z_, t_, a_, b_, c_] := N[(x / N[(x + N[(y * N[Exp[N[(2.0 * N[(N[(N[(z * N[Sqrt[N[(t + a), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] / t), $MachinePrecision] - N[(N[(b - c), $MachinePrecision] * N[(N[(a + N[(5.0 / 6.0), $MachinePrecision]), $MachinePrecision] - N[(2.0 / N[(t * 3.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{x}{x + y \cdot e^{2 \cdot \left(\frac{z \cdot \sqrt{t + a}}{t} - \left(b - c\right) \cdot \left(\left(a + \frac{5}{6}\right) - \frac{2}{t \cdot 3}\right)\right)}}
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 23 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (x y z t a b c)
:precision binary64
(/
x
(+
x
(*
y
(exp
(*
2.0
(-
(/ (* z (sqrt (+ t a))) t)
(* (- b c) (- (+ a (/ 5.0 6.0)) (/ 2.0 (* t 3.0)))))))))))
double code(double x, double y, double z, double t, double a, double b, double c) {
return x / (x + (y * exp((2.0 * (((z * sqrt((t + a))) / t) - ((b - c) * ((a + (5.0 / 6.0)) - (2.0 / (t * 3.0)))))))));
}
real(8) function code(x, y, z, t, a, b, c)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8), intent (in) :: c
code = x / (x + (y * exp((2.0d0 * (((z * sqrt((t + a))) / t) - ((b - c) * ((a + (5.0d0 / 6.0d0)) - (2.0d0 / (t * 3.0d0)))))))))
end function
public static double code(double x, double y, double z, double t, double a, double b, double c) {
return x / (x + (y * Math.exp((2.0 * (((z * Math.sqrt((t + a))) / t) - ((b - c) * ((a + (5.0 / 6.0)) - (2.0 / (t * 3.0)))))))));
}
def code(x, y, z, t, a, b, c): return x / (x + (y * math.exp((2.0 * (((z * math.sqrt((t + a))) / t) - ((b - c) * ((a + (5.0 / 6.0)) - (2.0 / (t * 3.0)))))))))
function code(x, y, z, t, a, b, c) return Float64(x / Float64(x + Float64(y * exp(Float64(2.0 * Float64(Float64(Float64(z * sqrt(Float64(t + a))) / t) - Float64(Float64(b - c) * Float64(Float64(a + Float64(5.0 / 6.0)) - Float64(2.0 / Float64(t * 3.0)))))))))) end
function tmp = code(x, y, z, t, a, b, c) tmp = x / (x + (y * exp((2.0 * (((z * sqrt((t + a))) / t) - ((b - c) * ((a + (5.0 / 6.0)) - (2.0 / (t * 3.0))))))))); end
code[x_, y_, z_, t_, a_, b_, c_] := N[(x / N[(x + N[(y * N[Exp[N[(2.0 * N[(N[(N[(z * N[Sqrt[N[(t + a), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] / t), $MachinePrecision] - N[(N[(b - c), $MachinePrecision] * N[(N[(a + N[(5.0 / 6.0), $MachinePrecision]), $MachinePrecision] - N[(2.0 / N[(t * 3.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{x}{x + y \cdot e^{2 \cdot \left(\frac{z \cdot \sqrt{t + a}}{t} - \left(b - c\right) \cdot \left(\left(a + \frac{5}{6}\right) - \frac{2}{t \cdot 3}\right)\right)}}
\end{array}
(FPCore (x y z t a b c)
:precision binary64
(/
x
(+
x
(*
y
(exp
(*
2.0
(+
(/ (* z (sqrt (+ t a))) t)
(* (- b c) (- (/ 2.0 (* t 3.0)) (+ a (/ 5.0 6.0)))))))))))
double code(double x, double y, double z, double t, double a, double b, double c) {
return x / (x + (y * exp((2.0 * (((z * sqrt((t + a))) / t) + ((b - c) * ((2.0 / (t * 3.0)) - (a + (5.0 / 6.0)))))))));
}
real(8) function code(x, y, z, t, a, b, c)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8), intent (in) :: c
code = x / (x + (y * exp((2.0d0 * (((z * sqrt((t + a))) / t) + ((b - c) * ((2.0d0 / (t * 3.0d0)) - (a + (5.0d0 / 6.0d0)))))))))
end function
public static double code(double x, double y, double z, double t, double a, double b, double c) {
return x / (x + (y * Math.exp((2.0 * (((z * Math.sqrt((t + a))) / t) + ((b - c) * ((2.0 / (t * 3.0)) - (a + (5.0 / 6.0)))))))));
}
def code(x, y, z, t, a, b, c): return x / (x + (y * math.exp((2.0 * (((z * math.sqrt((t + a))) / t) + ((b - c) * ((2.0 / (t * 3.0)) - (a + (5.0 / 6.0)))))))))
function code(x, y, z, t, a, b, c) return Float64(x / Float64(x + Float64(y * exp(Float64(2.0 * Float64(Float64(Float64(z * sqrt(Float64(t + a))) / t) + Float64(Float64(b - c) * Float64(Float64(2.0 / Float64(t * 3.0)) - Float64(a + Float64(5.0 / 6.0)))))))))) end
function tmp = code(x, y, z, t, a, b, c) tmp = x / (x + (y * exp((2.0 * (((z * sqrt((t + a))) / t) + ((b - c) * ((2.0 / (t * 3.0)) - (a + (5.0 / 6.0))))))))); end
code[x_, y_, z_, t_, a_, b_, c_] := N[(x / N[(x + N[(y * N[Exp[N[(2.0 * N[(N[(N[(z * N[Sqrt[N[(t + a), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] / t), $MachinePrecision] + N[(N[(b - c), $MachinePrecision] * N[(N[(2.0 / N[(t * 3.0), $MachinePrecision]), $MachinePrecision] - N[(a + N[(5.0 / 6.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{x}{x + y \cdot e^{2 \cdot \left(\frac{z \cdot \sqrt{t + a}}{t} + \left(b - c\right) \cdot \left(\frac{2}{t \cdot 3} - \left(a + \frac{5}{6}\right)\right)\right)}}
\end{array}
Initial program 98.1%
Final simplification98.1%
(FPCore (x y z t a b c)
:precision binary64
(let* ((t_1 (sqrt (+ t a))))
(if (<=
(/
x
(+
x
(*
y
(exp
(*
2.0
(+
(/ (* z t_1) t)
(* (- b c) (- (/ 2.0 (* t 3.0)) (+ a (/ 5.0 6.0))))))))))
0.9999999999999916)
(/
x
(*
y
(+
(exp
(*
2.0
(fma
(/ z t)
t_1
(*
(+ a (+ 0.8333333333333334 (/ -0.6666666666666666 t)))
(- c b)))))
(/ x y))))
1.0)))
double code(double x, double y, double z, double t, double a, double b, double c) {
double t_1 = sqrt((t + a));
double tmp;
if ((x / (x + (y * exp((2.0 * (((z * t_1) / t) + ((b - c) * ((2.0 / (t * 3.0)) - (a + (5.0 / 6.0)))))))))) <= 0.9999999999999916) {
tmp = x / (y * (exp((2.0 * fma((z / t), t_1, ((a + (0.8333333333333334 + (-0.6666666666666666 / t))) * (c - b))))) + (x / y)));
} else {
tmp = 1.0;
}
return tmp;
}
function code(x, y, z, t, a, b, c) t_1 = sqrt(Float64(t + a)) tmp = 0.0 if (Float64(x / Float64(x + Float64(y * exp(Float64(2.0 * Float64(Float64(Float64(z * t_1) / t) + Float64(Float64(b - c) * Float64(Float64(2.0 / Float64(t * 3.0)) - Float64(a + Float64(5.0 / 6.0)))))))))) <= 0.9999999999999916) tmp = Float64(x / Float64(y * Float64(exp(Float64(2.0 * fma(Float64(z / t), t_1, Float64(Float64(a + Float64(0.8333333333333334 + Float64(-0.6666666666666666 / t))) * Float64(c - b))))) + Float64(x / y)))); else tmp = 1.0; end return tmp end
code[x_, y_, z_, t_, a_, b_, c_] := Block[{t$95$1 = N[Sqrt[N[(t + a), $MachinePrecision]], $MachinePrecision]}, If[LessEqual[N[(x / N[(x + N[(y * N[Exp[N[(2.0 * N[(N[(N[(z * t$95$1), $MachinePrecision] / t), $MachinePrecision] + N[(N[(b - c), $MachinePrecision] * N[(N[(2.0 / N[(t * 3.0), $MachinePrecision]), $MachinePrecision] - N[(a + N[(5.0 / 6.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 0.9999999999999916], N[(x / N[(y * N[(N[Exp[N[(2.0 * N[(N[(z / t), $MachinePrecision] * t$95$1 + N[(N[(a + N[(0.8333333333333334 + N[(-0.6666666666666666 / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(c - b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] + N[(x / y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 1.0]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \sqrt{t + a}\\
\mathbf{if}\;\frac{x}{x + y \cdot e^{2 \cdot \left(\frac{z \cdot t\_1}{t} + \left(b - c\right) \cdot \left(\frac{2}{t \cdot 3} - \left(a + \frac{5}{6}\right)\right)\right)}} \leq 0.9999999999999916:\\
\;\;\;\;\frac{x}{y \cdot \left(e^{2 \cdot \mathsf{fma}\left(\frac{z}{t}, t\_1, \left(a + \left(0.8333333333333334 + \frac{-0.6666666666666666}{t}\right)\right) \cdot \left(c - b\right)\right)} + \frac{x}{y}\right)}\\
\mathbf{else}:\\
\;\;\;\;1\\
\end{array}
\end{array}
if (/.f64 x (+.f64 x (*.f64 y (exp.f64 (*.f64 #s(literal 2 binary64) (-.f64 (/.f64 (*.f64 z (sqrt.f64 (+.f64 t a))) t) (*.f64 (-.f64 b c) (-.f64 (+.f64 a (/.f64 #s(literal 5 binary64) #s(literal 6 binary64))) (/.f64 #s(literal 2 binary64) (*.f64 t #s(literal 3 binary64))))))))))) < 0.99999999999999156Initial program 99.3%
Taylor expanded in y around inf
Simplified91.0%
if 0.99999999999999156 < (/.f64 x (+.f64 x (*.f64 y (exp.f64 (*.f64 #s(literal 2 binary64) (-.f64 (/.f64 (*.f64 z (sqrt.f64 (+.f64 t a))) t) (*.f64 (-.f64 b c) (-.f64 (+.f64 a (/.f64 #s(literal 5 binary64) #s(literal 6 binary64))) (/.f64 #s(literal 2 binary64) (*.f64 t #s(literal 3 binary64))))))))))) Initial program 96.4%
Taylor expanded in x around inf
Simplified98.3%
*-inverses98.3
Applied egg-rr98.3%
Final simplification94.2%
(FPCore (x y z t a b c)
:precision binary64
(let* ((t_1
(+
(/ (* z (sqrt (+ t a))) t)
(* (- b c) (- (/ 2.0 (* t 3.0)) (+ a (/ 5.0 6.0))))))
(t_2 (+ 0.8333333333333334 (- a (/ 0.6666666666666666 t))))
(t_3 (- (/ 0.6666666666666666 t) (+ a 0.8333333333333334))))
(if (<= t_1 -5e+126)
1.0
(if (<= t_1 2e+177)
(/ x (+ x (* y (exp (* 2.0 (* b t_3))))))
(if (<= t_1 4e+288)
(/ x (+ x (fma c (* 2.0 (fma c (* y (* t_2 t_2)) (* y t_2))) y)))
(/ x (+ x (* y (fma b (* 2.0 (fma b (* t_3 t_3) t_3)) 1.0)))))))))
double code(double x, double y, double z, double t, double a, double b, double c) {
double t_1 = ((z * sqrt((t + a))) / t) + ((b - c) * ((2.0 / (t * 3.0)) - (a + (5.0 / 6.0))));
double t_2 = 0.8333333333333334 + (a - (0.6666666666666666 / t));
double t_3 = (0.6666666666666666 / t) - (a + 0.8333333333333334);
double tmp;
if (t_1 <= -5e+126) {
tmp = 1.0;
} else if (t_1 <= 2e+177) {
tmp = x / (x + (y * exp((2.0 * (b * t_3)))));
} else if (t_1 <= 4e+288) {
tmp = x / (x + fma(c, (2.0 * fma(c, (y * (t_2 * t_2)), (y * t_2))), y));
} else {
tmp = x / (x + (y * fma(b, (2.0 * fma(b, (t_3 * t_3), t_3)), 1.0)));
}
return tmp;
}
function code(x, y, z, t, a, b, c) t_1 = Float64(Float64(Float64(z * sqrt(Float64(t + a))) / t) + Float64(Float64(b - c) * Float64(Float64(2.0 / Float64(t * 3.0)) - Float64(a + Float64(5.0 / 6.0))))) t_2 = Float64(0.8333333333333334 + Float64(a - Float64(0.6666666666666666 / t))) t_3 = Float64(Float64(0.6666666666666666 / t) - Float64(a + 0.8333333333333334)) tmp = 0.0 if (t_1 <= -5e+126) tmp = 1.0; elseif (t_1 <= 2e+177) tmp = Float64(x / Float64(x + Float64(y * exp(Float64(2.0 * Float64(b * t_3)))))); elseif (t_1 <= 4e+288) tmp = Float64(x / Float64(x + fma(c, Float64(2.0 * fma(c, Float64(y * Float64(t_2 * t_2)), Float64(y * t_2))), y))); else tmp = Float64(x / Float64(x + Float64(y * fma(b, Float64(2.0 * fma(b, Float64(t_3 * t_3), t_3)), 1.0)))); end return tmp end
code[x_, y_, z_, t_, a_, b_, c_] := Block[{t$95$1 = N[(N[(N[(z * N[Sqrt[N[(t + a), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] / t), $MachinePrecision] + N[(N[(b - c), $MachinePrecision] * N[(N[(2.0 / N[(t * 3.0), $MachinePrecision]), $MachinePrecision] - N[(a + N[(5.0 / 6.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(0.8333333333333334 + N[(a - N[(0.6666666666666666 / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$3 = N[(N[(0.6666666666666666 / t), $MachinePrecision] - N[(a + 0.8333333333333334), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$1, -5e+126], 1.0, If[LessEqual[t$95$1, 2e+177], N[(x / N[(x + N[(y * N[Exp[N[(2.0 * N[(b * t$95$3), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$1, 4e+288], N[(x / N[(x + N[(c * N[(2.0 * N[(c * N[(y * N[(t$95$2 * t$95$2), $MachinePrecision]), $MachinePrecision] + N[(y * t$95$2), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x / N[(x + N[(y * N[(b * N[(2.0 * N[(b * N[(t$95$3 * t$95$3), $MachinePrecision] + t$95$3), $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{z \cdot \sqrt{t + a}}{t} + \left(b - c\right) \cdot \left(\frac{2}{t \cdot 3} - \left(a + \frac{5}{6}\right)\right)\\
t_2 := 0.8333333333333334 + \left(a - \frac{0.6666666666666666}{t}\right)\\
t_3 := \frac{0.6666666666666666}{t} - \left(a + 0.8333333333333334\right)\\
\mathbf{if}\;t\_1 \leq -5 \cdot 10^{+126}:\\
\;\;\;\;1\\
\mathbf{elif}\;t\_1 \leq 2 \cdot 10^{+177}:\\
\;\;\;\;\frac{x}{x + y \cdot e^{2 \cdot \left(b \cdot t\_3\right)}}\\
\mathbf{elif}\;t\_1 \leq 4 \cdot 10^{+288}:\\
\;\;\;\;\frac{x}{x + \mathsf{fma}\left(c, 2 \cdot \mathsf{fma}\left(c, y \cdot \left(t\_2 \cdot t\_2\right), y \cdot t\_2\right), y\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{x}{x + y \cdot \mathsf{fma}\left(b, 2 \cdot \mathsf{fma}\left(b, t\_3 \cdot t\_3, t\_3\right), 1\right)}\\
\end{array}
\end{array}
if (-.f64 (/.f64 (*.f64 z (sqrt.f64 (+.f64 t a))) t) (*.f64 (-.f64 b c) (-.f64 (+.f64 a (/.f64 #s(literal 5 binary64) #s(literal 6 binary64))) (/.f64 #s(literal 2 binary64) (*.f64 t #s(literal 3 binary64)))))) < -4.99999999999999977e126Initial program 99.0%
Taylor expanded in x around inf
Simplified99.0%
*-inverses99.0
Applied egg-rr99.0%
if -4.99999999999999977e126 < (-.f64 (/.f64 (*.f64 z (sqrt.f64 (+.f64 t a))) t) (*.f64 (-.f64 b c) (-.f64 (+.f64 a (/.f64 #s(literal 5 binary64) #s(literal 6 binary64))) (/.f64 #s(literal 2 binary64) (*.f64 t #s(literal 3 binary64)))))) < 2e177Initial program 99.9%
Taylor expanded in b around inf
*-lowering-*.f64N/A
--lowering--.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f64N/A
+-commutativeN/A
+-lowering-+.f6476.7
Simplified76.7%
if 2e177 < (-.f64 (/.f64 (*.f64 z (sqrt.f64 (+.f64 t a))) t) (*.f64 (-.f64 b c) (-.f64 (+.f64 a (/.f64 #s(literal 5 binary64) #s(literal 6 binary64))) (/.f64 #s(literal 2 binary64) (*.f64 t #s(literal 3 binary64)))))) < 4e288Initial program 100.0%
Taylor expanded in c around inf
*-lowering-*.f64N/A
+-commutativeN/A
associate--l+N/A
+-lowering-+.f64N/A
sub-negN/A
+-lowering-+.f64N/A
associate-*r/N/A
metadata-evalN/A
distribute-neg-fracN/A
metadata-evalN/A
/-lowering-/.f6482.0
Simplified82.0%
Taylor expanded in c around 0
+-lowering-+.f64N/A
+-commutativeN/A
accelerator-lowering-fma.f64N/A
Simplified79.1%
if 4e288 < (-.f64 (/.f64 (*.f64 z (sqrt.f64 (+.f64 t a))) t) (*.f64 (-.f64 b c) (-.f64 (+.f64 a (/.f64 #s(literal 5 binary64) #s(literal 6 binary64))) (/.f64 #s(literal 2 binary64) (*.f64 t #s(literal 3 binary64)))))) Initial program 94.8%
Taylor expanded in b around inf
*-lowering-*.f64N/A
--lowering--.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f64N/A
+-commutativeN/A
+-lowering-+.f6468.3
Simplified68.3%
Taylor expanded in b around 0
+-commutativeN/A
accelerator-lowering-fma.f64N/A
Simplified83.6%
Final simplification87.3%
(FPCore (x y z t a b c)
:precision binary64
(let* ((t_1 (- (/ 0.6666666666666666 t) (+ a 0.8333333333333334)))
(t_2 (+ 0.8333333333333334 (- a (/ 0.6666666666666666 t))))
(t_3
(+
(/ (* z (sqrt (+ t a))) t)
(* (- b c) (- (/ 2.0 (* t 3.0)) (+ a (/ 5.0 6.0)))))))
(if (<= t_3 -200.0)
1.0
(if (<= t_3 4e+146)
(/ x (fma y (exp (* 2.0 (* c (+ a 0.8333333333333334)))) x))
(if (<= t_3 4e+288)
(/ x (+ x (fma c (* 2.0 (fma c (* y (* t_2 t_2)) (* y t_2))) y)))
(/ x (+ x (* y (fma b (* 2.0 (fma b (* t_1 t_1) t_1)) 1.0)))))))))
double code(double x, double y, double z, double t, double a, double b, double c) {
double t_1 = (0.6666666666666666 / t) - (a + 0.8333333333333334);
double t_2 = 0.8333333333333334 + (a - (0.6666666666666666 / t));
double t_3 = ((z * sqrt((t + a))) / t) + ((b - c) * ((2.0 / (t * 3.0)) - (a + (5.0 / 6.0))));
double tmp;
if (t_3 <= -200.0) {
tmp = 1.0;
} else if (t_3 <= 4e+146) {
tmp = x / fma(y, exp((2.0 * (c * (a + 0.8333333333333334)))), x);
} else if (t_3 <= 4e+288) {
tmp = x / (x + fma(c, (2.0 * fma(c, (y * (t_2 * t_2)), (y * t_2))), y));
} else {
tmp = x / (x + (y * fma(b, (2.0 * fma(b, (t_1 * t_1), t_1)), 1.0)));
}
return tmp;
}
function code(x, y, z, t, a, b, c) t_1 = Float64(Float64(0.6666666666666666 / t) - Float64(a + 0.8333333333333334)) t_2 = Float64(0.8333333333333334 + Float64(a - Float64(0.6666666666666666 / t))) t_3 = Float64(Float64(Float64(z * sqrt(Float64(t + a))) / t) + Float64(Float64(b - c) * Float64(Float64(2.0 / Float64(t * 3.0)) - Float64(a + Float64(5.0 / 6.0))))) tmp = 0.0 if (t_3 <= -200.0) tmp = 1.0; elseif (t_3 <= 4e+146) tmp = Float64(x / fma(y, exp(Float64(2.0 * Float64(c * Float64(a + 0.8333333333333334)))), x)); elseif (t_3 <= 4e+288) tmp = Float64(x / Float64(x + fma(c, Float64(2.0 * fma(c, Float64(y * Float64(t_2 * t_2)), Float64(y * t_2))), y))); else tmp = Float64(x / Float64(x + Float64(y * fma(b, Float64(2.0 * fma(b, Float64(t_1 * t_1), t_1)), 1.0)))); end return tmp end
code[x_, y_, z_, t_, a_, b_, c_] := Block[{t$95$1 = N[(N[(0.6666666666666666 / t), $MachinePrecision] - N[(a + 0.8333333333333334), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(0.8333333333333334 + N[(a - N[(0.6666666666666666 / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$3 = N[(N[(N[(z * N[Sqrt[N[(t + a), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] / t), $MachinePrecision] + N[(N[(b - c), $MachinePrecision] * N[(N[(2.0 / N[(t * 3.0), $MachinePrecision]), $MachinePrecision] - N[(a + N[(5.0 / 6.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$3, -200.0], 1.0, If[LessEqual[t$95$3, 4e+146], N[(x / N[(y * N[Exp[N[(2.0 * N[(c * N[(a + 0.8333333333333334), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] + x), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$3, 4e+288], N[(x / N[(x + N[(c * N[(2.0 * N[(c * N[(y * N[(t$95$2 * t$95$2), $MachinePrecision]), $MachinePrecision] + N[(y * t$95$2), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x / N[(x + N[(y * N[(b * N[(2.0 * N[(b * N[(t$95$1 * t$95$1), $MachinePrecision] + t$95$1), $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{0.6666666666666666}{t} - \left(a + 0.8333333333333334\right)\\
t_2 := 0.8333333333333334 + \left(a - \frac{0.6666666666666666}{t}\right)\\
t_3 := \frac{z \cdot \sqrt{t + a}}{t} + \left(b - c\right) \cdot \left(\frac{2}{t \cdot 3} - \left(a + \frac{5}{6}\right)\right)\\
\mathbf{if}\;t\_3 \leq -200:\\
\;\;\;\;1\\
\mathbf{elif}\;t\_3 \leq 4 \cdot 10^{+146}:\\
\;\;\;\;\frac{x}{\mathsf{fma}\left(y, e^{2 \cdot \left(c \cdot \left(a + 0.8333333333333334\right)\right)}, x\right)}\\
\mathbf{elif}\;t\_3 \leq 4 \cdot 10^{+288}:\\
\;\;\;\;\frac{x}{x + \mathsf{fma}\left(c, 2 \cdot \mathsf{fma}\left(c, y \cdot \left(t\_2 \cdot t\_2\right), y \cdot t\_2\right), y\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{x}{x + y \cdot \mathsf{fma}\left(b, 2 \cdot \mathsf{fma}\left(b, t\_1 \cdot t\_1, t\_1\right), 1\right)}\\
\end{array}
\end{array}
if (-.f64 (/.f64 (*.f64 z (sqrt.f64 (+.f64 t a))) t) (*.f64 (-.f64 b c) (-.f64 (+.f64 a (/.f64 #s(literal 5 binary64) #s(literal 6 binary64))) (/.f64 #s(literal 2 binary64) (*.f64 t #s(literal 3 binary64)))))) < -200Initial program 99.1%
Taylor expanded in x around inf
Simplified99.1%
*-inverses99.1
Applied egg-rr99.1%
if -200 < (-.f64 (/.f64 (*.f64 z (sqrt.f64 (+.f64 t a))) t) (*.f64 (-.f64 b c) (-.f64 (+.f64 a (/.f64 #s(literal 5 binary64) #s(literal 6 binary64))) (/.f64 #s(literal 2 binary64) (*.f64 t #s(literal 3 binary64)))))) < 3.99999999999999973e146Initial program 99.9%
Taylor expanded in c around inf
*-lowering-*.f64N/A
+-commutativeN/A
associate--l+N/A
+-lowering-+.f64N/A
sub-negN/A
+-lowering-+.f64N/A
associate-*r/N/A
metadata-evalN/A
distribute-neg-fracN/A
metadata-evalN/A
/-lowering-/.f6475.1
Simplified75.1%
Taylor expanded in t around inf
+-commutativeN/A
accelerator-lowering-fma.f64N/A
exp-lowering-exp.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f6470.8
Simplified70.8%
if 3.99999999999999973e146 < (-.f64 (/.f64 (*.f64 z (sqrt.f64 (+.f64 t a))) t) (*.f64 (-.f64 b c) (-.f64 (+.f64 a (/.f64 #s(literal 5 binary64) #s(literal 6 binary64))) (/.f64 #s(literal 2 binary64) (*.f64 t #s(literal 3 binary64)))))) < 4e288Initial program 100.0%
Taylor expanded in c around inf
*-lowering-*.f64N/A
+-commutativeN/A
associate--l+N/A
+-lowering-+.f64N/A
sub-negN/A
+-lowering-+.f64N/A
associate-*r/N/A
metadata-evalN/A
distribute-neg-fracN/A
metadata-evalN/A
/-lowering-/.f6473.0
Simplified73.0%
Taylor expanded in c around 0
+-lowering-+.f64N/A
+-commutativeN/A
accelerator-lowering-fma.f64N/A
Simplified75.6%
if 4e288 < (-.f64 (/.f64 (*.f64 z (sqrt.f64 (+.f64 t a))) t) (*.f64 (-.f64 b c) (-.f64 (+.f64 a (/.f64 #s(literal 5 binary64) #s(literal 6 binary64))) (/.f64 #s(literal 2 binary64) (*.f64 t #s(literal 3 binary64)))))) Initial program 94.8%
Taylor expanded in b around inf
*-lowering-*.f64N/A
--lowering--.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f64N/A
+-commutativeN/A
+-lowering-+.f6468.3
Simplified68.3%
Taylor expanded in b around 0
+-commutativeN/A
accelerator-lowering-fma.f64N/A
Simplified83.6%
Final simplification86.6%
(FPCore (x y z t a b c)
:precision binary64
(let* ((t_1 (- (/ 0.6666666666666666 t) (+ a 0.8333333333333334)))
(t_2 (+ 0.8333333333333334 (- a (/ 0.6666666666666666 t))))
(t_3
(+
(/ (* z (sqrt (+ t a))) t)
(* (- b c) (- (/ 2.0 (* t 3.0)) (+ a (/ 5.0 6.0)))))))
(if (<= t_3 -200.0)
1.0
(if (<= t_3 2e+138)
(/ x (fma y (exp (* -2.0 (* b (+ a 0.8333333333333334)))) x))
(if (<= t_3 4e+288)
(/ x (+ x (fma c (* 2.0 (fma c (* y (* t_2 t_2)) (* y t_2))) y)))
(/ x (+ x (* y (fma b (* 2.0 (fma b (* t_1 t_1) t_1)) 1.0)))))))))
double code(double x, double y, double z, double t, double a, double b, double c) {
double t_1 = (0.6666666666666666 / t) - (a + 0.8333333333333334);
double t_2 = 0.8333333333333334 + (a - (0.6666666666666666 / t));
double t_3 = ((z * sqrt((t + a))) / t) + ((b - c) * ((2.0 / (t * 3.0)) - (a + (5.0 / 6.0))));
double tmp;
if (t_3 <= -200.0) {
tmp = 1.0;
} else if (t_3 <= 2e+138) {
tmp = x / fma(y, exp((-2.0 * (b * (a + 0.8333333333333334)))), x);
} else if (t_3 <= 4e+288) {
tmp = x / (x + fma(c, (2.0 * fma(c, (y * (t_2 * t_2)), (y * t_2))), y));
} else {
tmp = x / (x + (y * fma(b, (2.0 * fma(b, (t_1 * t_1), t_1)), 1.0)));
}
return tmp;
}
function code(x, y, z, t, a, b, c) t_1 = Float64(Float64(0.6666666666666666 / t) - Float64(a + 0.8333333333333334)) t_2 = Float64(0.8333333333333334 + Float64(a - Float64(0.6666666666666666 / t))) t_3 = Float64(Float64(Float64(z * sqrt(Float64(t + a))) / t) + Float64(Float64(b - c) * Float64(Float64(2.0 / Float64(t * 3.0)) - Float64(a + Float64(5.0 / 6.0))))) tmp = 0.0 if (t_3 <= -200.0) tmp = 1.0; elseif (t_3 <= 2e+138) tmp = Float64(x / fma(y, exp(Float64(-2.0 * Float64(b * Float64(a + 0.8333333333333334)))), x)); elseif (t_3 <= 4e+288) tmp = Float64(x / Float64(x + fma(c, Float64(2.0 * fma(c, Float64(y * Float64(t_2 * t_2)), Float64(y * t_2))), y))); else tmp = Float64(x / Float64(x + Float64(y * fma(b, Float64(2.0 * fma(b, Float64(t_1 * t_1), t_1)), 1.0)))); end return tmp end
code[x_, y_, z_, t_, a_, b_, c_] := Block[{t$95$1 = N[(N[(0.6666666666666666 / t), $MachinePrecision] - N[(a + 0.8333333333333334), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(0.8333333333333334 + N[(a - N[(0.6666666666666666 / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$3 = N[(N[(N[(z * N[Sqrt[N[(t + a), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] / t), $MachinePrecision] + N[(N[(b - c), $MachinePrecision] * N[(N[(2.0 / N[(t * 3.0), $MachinePrecision]), $MachinePrecision] - N[(a + N[(5.0 / 6.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$3, -200.0], 1.0, If[LessEqual[t$95$3, 2e+138], N[(x / N[(y * N[Exp[N[(-2.0 * N[(b * N[(a + 0.8333333333333334), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] + x), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$3, 4e+288], N[(x / N[(x + N[(c * N[(2.0 * N[(c * N[(y * N[(t$95$2 * t$95$2), $MachinePrecision]), $MachinePrecision] + N[(y * t$95$2), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x / N[(x + N[(y * N[(b * N[(2.0 * N[(b * N[(t$95$1 * t$95$1), $MachinePrecision] + t$95$1), $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{0.6666666666666666}{t} - \left(a + 0.8333333333333334\right)\\
t_2 := 0.8333333333333334 + \left(a - \frac{0.6666666666666666}{t}\right)\\
t_3 := \frac{z \cdot \sqrt{t + a}}{t} + \left(b - c\right) \cdot \left(\frac{2}{t \cdot 3} - \left(a + \frac{5}{6}\right)\right)\\
\mathbf{if}\;t\_3 \leq -200:\\
\;\;\;\;1\\
\mathbf{elif}\;t\_3 \leq 2 \cdot 10^{+138}:\\
\;\;\;\;\frac{x}{\mathsf{fma}\left(y, e^{-2 \cdot \left(b \cdot \left(a + 0.8333333333333334\right)\right)}, x\right)}\\
\mathbf{elif}\;t\_3 \leq 4 \cdot 10^{+288}:\\
\;\;\;\;\frac{x}{x + \mathsf{fma}\left(c, 2 \cdot \mathsf{fma}\left(c, y \cdot \left(t\_2 \cdot t\_2\right), y \cdot t\_2\right), y\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{x}{x + y \cdot \mathsf{fma}\left(b, 2 \cdot \mathsf{fma}\left(b, t\_1 \cdot t\_1, t\_1\right), 1\right)}\\
\end{array}
\end{array}
if (-.f64 (/.f64 (*.f64 z (sqrt.f64 (+.f64 t a))) t) (*.f64 (-.f64 b c) (-.f64 (+.f64 a (/.f64 #s(literal 5 binary64) #s(literal 6 binary64))) (/.f64 #s(literal 2 binary64) (*.f64 t #s(literal 3 binary64)))))) < -200Initial program 99.1%
Taylor expanded in x around inf
Simplified99.1%
*-inverses99.1
Applied egg-rr99.1%
if -200 < (-.f64 (/.f64 (*.f64 z (sqrt.f64 (+.f64 t a))) t) (*.f64 (-.f64 b c) (-.f64 (+.f64 a (/.f64 #s(literal 5 binary64) #s(literal 6 binary64))) (/.f64 #s(literal 2 binary64) (*.f64 t #s(literal 3 binary64)))))) < 2.0000000000000001e138Initial program 99.9%
Taylor expanded in b around inf
*-lowering-*.f64N/A
--lowering--.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f64N/A
+-commutativeN/A
+-lowering-+.f6473.0
Simplified73.0%
Taylor expanded in t around inf
+-commutativeN/A
accelerator-lowering-fma.f64N/A
exp-lowering-exp.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f6470.3
Simplified70.3%
if 2.0000000000000001e138 < (-.f64 (/.f64 (*.f64 z (sqrt.f64 (+.f64 t a))) t) (*.f64 (-.f64 b c) (-.f64 (+.f64 a (/.f64 #s(literal 5 binary64) #s(literal 6 binary64))) (/.f64 #s(literal 2 binary64) (*.f64 t #s(literal 3 binary64)))))) < 4e288Initial program 100.0%
Taylor expanded in c around inf
*-lowering-*.f64N/A
+-commutativeN/A
associate--l+N/A
+-lowering-+.f64N/A
sub-negN/A
+-lowering-+.f64N/A
associate-*r/N/A
metadata-evalN/A
distribute-neg-fracN/A
metadata-evalN/A
/-lowering-/.f6473.3
Simplified73.3%
Taylor expanded in c around 0
+-lowering-+.f64N/A
+-commutativeN/A
accelerator-lowering-fma.f64N/A
Simplified73.5%
if 4e288 < (-.f64 (/.f64 (*.f64 z (sqrt.f64 (+.f64 t a))) t) (*.f64 (-.f64 b c) (-.f64 (+.f64 a (/.f64 #s(literal 5 binary64) #s(literal 6 binary64))) (/.f64 #s(literal 2 binary64) (*.f64 t #s(literal 3 binary64)))))) Initial program 94.8%
Taylor expanded in b around inf
*-lowering-*.f64N/A
--lowering--.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f64N/A
+-commutativeN/A
+-lowering-+.f6468.3
Simplified68.3%
Taylor expanded in b around 0
+-commutativeN/A
accelerator-lowering-fma.f64N/A
Simplified83.6%
Final simplification86.2%
(FPCore (x y z t a b c)
:precision binary64
(let* ((t_1 (- (/ 0.6666666666666666 t) (+ a 0.8333333333333334)))
(t_2 (+ 0.8333333333333334 (- a (/ 0.6666666666666666 t))))
(t_3
(+
(/ (* z (sqrt (+ t a))) t)
(* (- b c) (- (/ 2.0 (* t 3.0)) (+ a (/ 5.0 6.0)))))))
(if (<= t_3 -200.0)
1.0
(if (<= t_3 2e+104)
(/ x (fma y (exp (* c 1.6666666666666667)) x))
(if (<= t_3 4e+288)
(/ x (+ x (fma c (* 2.0 (fma c (* y (* t_2 t_2)) (* y t_2))) y)))
(/ x (+ x (* y (fma b (* 2.0 (fma b (* t_1 t_1) t_1)) 1.0)))))))))
double code(double x, double y, double z, double t, double a, double b, double c) {
double t_1 = (0.6666666666666666 / t) - (a + 0.8333333333333334);
double t_2 = 0.8333333333333334 + (a - (0.6666666666666666 / t));
double t_3 = ((z * sqrt((t + a))) / t) + ((b - c) * ((2.0 / (t * 3.0)) - (a + (5.0 / 6.0))));
double tmp;
if (t_3 <= -200.0) {
tmp = 1.0;
} else if (t_3 <= 2e+104) {
tmp = x / fma(y, exp((c * 1.6666666666666667)), x);
} else if (t_3 <= 4e+288) {
tmp = x / (x + fma(c, (2.0 * fma(c, (y * (t_2 * t_2)), (y * t_2))), y));
} else {
tmp = x / (x + (y * fma(b, (2.0 * fma(b, (t_1 * t_1), t_1)), 1.0)));
}
return tmp;
}
function code(x, y, z, t, a, b, c) t_1 = Float64(Float64(0.6666666666666666 / t) - Float64(a + 0.8333333333333334)) t_2 = Float64(0.8333333333333334 + Float64(a - Float64(0.6666666666666666 / t))) t_3 = Float64(Float64(Float64(z * sqrt(Float64(t + a))) / t) + Float64(Float64(b - c) * Float64(Float64(2.0 / Float64(t * 3.0)) - Float64(a + Float64(5.0 / 6.0))))) tmp = 0.0 if (t_3 <= -200.0) tmp = 1.0; elseif (t_3 <= 2e+104) tmp = Float64(x / fma(y, exp(Float64(c * 1.6666666666666667)), x)); elseif (t_3 <= 4e+288) tmp = Float64(x / Float64(x + fma(c, Float64(2.0 * fma(c, Float64(y * Float64(t_2 * t_2)), Float64(y * t_2))), y))); else tmp = Float64(x / Float64(x + Float64(y * fma(b, Float64(2.0 * fma(b, Float64(t_1 * t_1), t_1)), 1.0)))); end return tmp end
code[x_, y_, z_, t_, a_, b_, c_] := Block[{t$95$1 = N[(N[(0.6666666666666666 / t), $MachinePrecision] - N[(a + 0.8333333333333334), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(0.8333333333333334 + N[(a - N[(0.6666666666666666 / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$3 = N[(N[(N[(z * N[Sqrt[N[(t + a), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] / t), $MachinePrecision] + N[(N[(b - c), $MachinePrecision] * N[(N[(2.0 / N[(t * 3.0), $MachinePrecision]), $MachinePrecision] - N[(a + N[(5.0 / 6.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$3, -200.0], 1.0, If[LessEqual[t$95$3, 2e+104], N[(x / N[(y * N[Exp[N[(c * 1.6666666666666667), $MachinePrecision]], $MachinePrecision] + x), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$3, 4e+288], N[(x / N[(x + N[(c * N[(2.0 * N[(c * N[(y * N[(t$95$2 * t$95$2), $MachinePrecision]), $MachinePrecision] + N[(y * t$95$2), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x / N[(x + N[(y * N[(b * N[(2.0 * N[(b * N[(t$95$1 * t$95$1), $MachinePrecision] + t$95$1), $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{0.6666666666666666}{t} - \left(a + 0.8333333333333334\right)\\
t_2 := 0.8333333333333334 + \left(a - \frac{0.6666666666666666}{t}\right)\\
t_3 := \frac{z \cdot \sqrt{t + a}}{t} + \left(b - c\right) \cdot \left(\frac{2}{t \cdot 3} - \left(a + \frac{5}{6}\right)\right)\\
\mathbf{if}\;t\_3 \leq -200:\\
\;\;\;\;1\\
\mathbf{elif}\;t\_3 \leq 2 \cdot 10^{+104}:\\
\;\;\;\;\frac{x}{\mathsf{fma}\left(y, e^{c \cdot 1.6666666666666667}, x\right)}\\
\mathbf{elif}\;t\_3 \leq 4 \cdot 10^{+288}:\\
\;\;\;\;\frac{x}{x + \mathsf{fma}\left(c, 2 \cdot \mathsf{fma}\left(c, y \cdot \left(t\_2 \cdot t\_2\right), y \cdot t\_2\right), y\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{x}{x + y \cdot \mathsf{fma}\left(b, 2 \cdot \mathsf{fma}\left(b, t\_1 \cdot t\_1, t\_1\right), 1\right)}\\
\end{array}
\end{array}
if (-.f64 (/.f64 (*.f64 z (sqrt.f64 (+.f64 t a))) t) (*.f64 (-.f64 b c) (-.f64 (+.f64 a (/.f64 #s(literal 5 binary64) #s(literal 6 binary64))) (/.f64 #s(literal 2 binary64) (*.f64 t #s(literal 3 binary64)))))) < -200Initial program 99.1%
Taylor expanded in x around inf
Simplified99.1%
*-inverses99.1
Applied egg-rr99.1%
if -200 < (-.f64 (/.f64 (*.f64 z (sqrt.f64 (+.f64 t a))) t) (*.f64 (-.f64 b c) (-.f64 (+.f64 a (/.f64 #s(literal 5 binary64) #s(literal 6 binary64))) (/.f64 #s(literal 2 binary64) (*.f64 t #s(literal 3 binary64)))))) < 2e104Initial program 99.8%
Taylor expanded in c around inf
*-lowering-*.f64N/A
+-commutativeN/A
associate--l+N/A
+-lowering-+.f64N/A
sub-negN/A
+-lowering-+.f64N/A
associate-*r/N/A
metadata-evalN/A
distribute-neg-fracN/A
metadata-evalN/A
/-lowering-/.f6483.7
Simplified83.7%
Taylor expanded in t around inf
+-commutativeN/A
accelerator-lowering-fma.f64N/A
exp-lowering-exp.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f6478.1
Simplified78.1%
Taylor expanded in a around 0
*-lowering-*.f6478.1
Simplified78.1%
if 2e104 < (-.f64 (/.f64 (*.f64 z (sqrt.f64 (+.f64 t a))) t) (*.f64 (-.f64 b c) (-.f64 (+.f64 a (/.f64 #s(literal 5 binary64) #s(literal 6 binary64))) (/.f64 #s(literal 2 binary64) (*.f64 t #s(literal 3 binary64)))))) < 4e288Initial program 100.0%
Taylor expanded in c around inf
*-lowering-*.f64N/A
+-commutativeN/A
associate--l+N/A
+-lowering-+.f64N/A
sub-negN/A
+-lowering-+.f64N/A
associate-*r/N/A
metadata-evalN/A
distribute-neg-fracN/A
metadata-evalN/A
/-lowering-/.f6468.0
Simplified68.0%
Taylor expanded in c around 0
+-lowering-+.f64N/A
+-commutativeN/A
accelerator-lowering-fma.f64N/A
Simplified68.2%
if 4e288 < (-.f64 (/.f64 (*.f64 z (sqrt.f64 (+.f64 t a))) t) (*.f64 (-.f64 b c) (-.f64 (+.f64 a (/.f64 #s(literal 5 binary64) #s(literal 6 binary64))) (/.f64 #s(literal 2 binary64) (*.f64 t #s(literal 3 binary64)))))) Initial program 94.8%
Taylor expanded in b around inf
*-lowering-*.f64N/A
--lowering--.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f64N/A
+-commutativeN/A
+-lowering-+.f6468.3
Simplified68.3%
Taylor expanded in b around 0
+-commutativeN/A
accelerator-lowering-fma.f64N/A
Simplified83.6%
Final simplification86.2%
(FPCore (x y z t a b c)
:precision binary64
(let* ((t_1 (- (/ 0.6666666666666666 t) (+ a 0.8333333333333334)))
(t_2
(+
(/ (* z (sqrt (+ t a))) t)
(* (- b c) (- (/ 2.0 (* t 3.0)) (+ a (/ 5.0 6.0)))))))
(if (<= t_2 -200.0)
1.0
(if (<= t_2 5e+280)
(/
x
(+
x
(*
y
(exp
(*
2.0
(* c (+ a (+ 0.8333333333333334 (/ -0.6666666666666666 t)))))))))
(/ x (+ x (* y (fma b (* 2.0 (fma b (* t_1 t_1) t_1)) 1.0))))))))
double code(double x, double y, double z, double t, double a, double b, double c) {
double t_1 = (0.6666666666666666 / t) - (a + 0.8333333333333334);
double t_2 = ((z * sqrt((t + a))) / t) + ((b - c) * ((2.0 / (t * 3.0)) - (a + (5.0 / 6.0))));
double tmp;
if (t_2 <= -200.0) {
tmp = 1.0;
} else if (t_2 <= 5e+280) {
tmp = x / (x + (y * exp((2.0 * (c * (a + (0.8333333333333334 + (-0.6666666666666666 / t))))))));
} else {
tmp = x / (x + (y * fma(b, (2.0 * fma(b, (t_1 * t_1), t_1)), 1.0)));
}
return tmp;
}
function code(x, y, z, t, a, b, c) t_1 = Float64(Float64(0.6666666666666666 / t) - Float64(a + 0.8333333333333334)) t_2 = Float64(Float64(Float64(z * sqrt(Float64(t + a))) / t) + Float64(Float64(b - c) * Float64(Float64(2.0 / Float64(t * 3.0)) - Float64(a + Float64(5.0 / 6.0))))) tmp = 0.0 if (t_2 <= -200.0) tmp = 1.0; elseif (t_2 <= 5e+280) tmp = Float64(x / Float64(x + Float64(y * exp(Float64(2.0 * Float64(c * Float64(a + Float64(0.8333333333333334 + Float64(-0.6666666666666666 / t))))))))); else tmp = Float64(x / Float64(x + Float64(y * fma(b, Float64(2.0 * fma(b, Float64(t_1 * t_1), t_1)), 1.0)))); end return tmp end
code[x_, y_, z_, t_, a_, b_, c_] := Block[{t$95$1 = N[(N[(0.6666666666666666 / t), $MachinePrecision] - N[(a + 0.8333333333333334), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(N[(N[(z * N[Sqrt[N[(t + a), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] / t), $MachinePrecision] + N[(N[(b - c), $MachinePrecision] * N[(N[(2.0 / N[(t * 3.0), $MachinePrecision]), $MachinePrecision] - N[(a + N[(5.0 / 6.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$2, -200.0], 1.0, If[LessEqual[t$95$2, 5e+280], N[(x / N[(x + N[(y * N[Exp[N[(2.0 * N[(c * N[(a + N[(0.8333333333333334 + N[(-0.6666666666666666 / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x / N[(x + N[(y * N[(b * N[(2.0 * N[(b * N[(t$95$1 * t$95$1), $MachinePrecision] + t$95$1), $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{0.6666666666666666}{t} - \left(a + 0.8333333333333334\right)\\
t_2 := \frac{z \cdot \sqrt{t + a}}{t} + \left(b - c\right) \cdot \left(\frac{2}{t \cdot 3} - \left(a + \frac{5}{6}\right)\right)\\
\mathbf{if}\;t\_2 \leq -200:\\
\;\;\;\;1\\
\mathbf{elif}\;t\_2 \leq 5 \cdot 10^{+280}:\\
\;\;\;\;\frac{x}{x + y \cdot e^{2 \cdot \left(c \cdot \left(a + \left(0.8333333333333334 + \frac{-0.6666666666666666}{t}\right)\right)\right)}}\\
\mathbf{else}:\\
\;\;\;\;\frac{x}{x + y \cdot \mathsf{fma}\left(b, 2 \cdot \mathsf{fma}\left(b, t\_1 \cdot t\_1, t\_1\right), 1\right)}\\
\end{array}
\end{array}
if (-.f64 (/.f64 (*.f64 z (sqrt.f64 (+.f64 t a))) t) (*.f64 (-.f64 b c) (-.f64 (+.f64 a (/.f64 #s(literal 5 binary64) #s(literal 6 binary64))) (/.f64 #s(literal 2 binary64) (*.f64 t #s(literal 3 binary64)))))) < -200Initial program 99.1%
Taylor expanded in x around inf
Simplified99.1%
*-inverses99.1
Applied egg-rr99.1%
if -200 < (-.f64 (/.f64 (*.f64 z (sqrt.f64 (+.f64 t a))) t) (*.f64 (-.f64 b c) (-.f64 (+.f64 a (/.f64 #s(literal 5 binary64) #s(literal 6 binary64))) (/.f64 #s(literal 2 binary64) (*.f64 t #s(literal 3 binary64)))))) < 5.0000000000000002e280Initial program 99.9%
Taylor expanded in c around inf
*-lowering-*.f64N/A
+-commutativeN/A
associate--l+N/A
+-lowering-+.f64N/A
sub-negN/A
+-lowering-+.f64N/A
associate-*r/N/A
metadata-evalN/A
distribute-neg-fracN/A
metadata-evalN/A
/-lowering-/.f6473.9
Simplified73.9%
if 5.0000000000000002e280 < (-.f64 (/.f64 (*.f64 z (sqrt.f64 (+.f64 t a))) t) (*.f64 (-.f64 b c) (-.f64 (+.f64 a (/.f64 #s(literal 5 binary64) #s(literal 6 binary64))) (/.f64 #s(literal 2 binary64) (*.f64 t #s(literal 3 binary64)))))) Initial program 95.0%
Taylor expanded in b around inf
*-lowering-*.f64N/A
--lowering--.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f64N/A
+-commutativeN/A
+-lowering-+.f6468.7
Simplified68.7%
Taylor expanded in b around 0
+-commutativeN/A
accelerator-lowering-fma.f64N/A
Simplified83.2%
Final simplification86.8%
(FPCore (x y z t a b c)
:precision binary64
(let* ((t_1 (- (/ 0.6666666666666666 t) (+ a 0.8333333333333334))))
(if (<=
(/
x
(+
x
(*
y
(exp
(*
2.0
(+
(/ (* z (sqrt (+ t a))) t)
(* (- b c) (- (/ 2.0 (* t 3.0)) (+ a (/ 5.0 6.0))))))))))
0.9999999999999916)
(/ x (+ x (* y (fma b (* 2.0 (fma b (* t_1 t_1) t_1)) 1.0))))
1.0)))
double code(double x, double y, double z, double t, double a, double b, double c) {
double t_1 = (0.6666666666666666 / t) - (a + 0.8333333333333334);
double tmp;
if ((x / (x + (y * exp((2.0 * (((z * sqrt((t + a))) / t) + ((b - c) * ((2.0 / (t * 3.0)) - (a + (5.0 / 6.0)))))))))) <= 0.9999999999999916) {
tmp = x / (x + (y * fma(b, (2.0 * fma(b, (t_1 * t_1), t_1)), 1.0)));
} else {
tmp = 1.0;
}
return tmp;
}
function code(x, y, z, t, a, b, c) t_1 = Float64(Float64(0.6666666666666666 / t) - Float64(a + 0.8333333333333334)) tmp = 0.0 if (Float64(x / Float64(x + Float64(y * exp(Float64(2.0 * Float64(Float64(Float64(z * sqrt(Float64(t + a))) / t) + Float64(Float64(b - c) * Float64(Float64(2.0 / Float64(t * 3.0)) - Float64(a + Float64(5.0 / 6.0)))))))))) <= 0.9999999999999916) tmp = Float64(x / Float64(x + Float64(y * fma(b, Float64(2.0 * fma(b, Float64(t_1 * t_1), t_1)), 1.0)))); else tmp = 1.0; end return tmp end
code[x_, y_, z_, t_, a_, b_, c_] := Block[{t$95$1 = N[(N[(0.6666666666666666 / t), $MachinePrecision] - N[(a + 0.8333333333333334), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[N[(x / N[(x + N[(y * N[Exp[N[(2.0 * N[(N[(N[(z * N[Sqrt[N[(t + a), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] / t), $MachinePrecision] + N[(N[(b - c), $MachinePrecision] * N[(N[(2.0 / N[(t * 3.0), $MachinePrecision]), $MachinePrecision] - N[(a + N[(5.0 / 6.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 0.9999999999999916], N[(x / N[(x + N[(y * N[(b * N[(2.0 * N[(b * N[(t$95$1 * t$95$1), $MachinePrecision] + t$95$1), $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 1.0]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{0.6666666666666666}{t} - \left(a + 0.8333333333333334\right)\\
\mathbf{if}\;\frac{x}{x + y \cdot e^{2 \cdot \left(\frac{z \cdot \sqrt{t + a}}{t} + \left(b - c\right) \cdot \left(\frac{2}{t \cdot 3} - \left(a + \frac{5}{6}\right)\right)\right)}} \leq 0.9999999999999916:\\
\;\;\;\;\frac{x}{x + y \cdot \mathsf{fma}\left(b, 2 \cdot \mathsf{fma}\left(b, t\_1 \cdot t\_1, t\_1\right), 1\right)}\\
\mathbf{else}:\\
\;\;\;\;1\\
\end{array}
\end{array}
if (/.f64 x (+.f64 x (*.f64 y (exp.f64 (*.f64 #s(literal 2 binary64) (-.f64 (/.f64 (*.f64 z (sqrt.f64 (+.f64 t a))) t) (*.f64 (-.f64 b c) (-.f64 (+.f64 a (/.f64 #s(literal 5 binary64) #s(literal 6 binary64))) (/.f64 #s(literal 2 binary64) (*.f64 t #s(literal 3 binary64))))))))))) < 0.99999999999999156Initial program 99.3%
Taylor expanded in b around inf
*-lowering-*.f64N/A
--lowering--.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f64N/A
+-commutativeN/A
+-lowering-+.f6463.0
Simplified63.0%
Taylor expanded in b around 0
+-commutativeN/A
accelerator-lowering-fma.f64N/A
Simplified70.6%
if 0.99999999999999156 < (/.f64 x (+.f64 x (*.f64 y (exp.f64 (*.f64 #s(literal 2 binary64) (-.f64 (/.f64 (*.f64 z (sqrt.f64 (+.f64 t a))) t) (*.f64 (-.f64 b c) (-.f64 (+.f64 a (/.f64 #s(literal 5 binary64) #s(literal 6 binary64))) (/.f64 #s(literal 2 binary64) (*.f64 t #s(literal 3 binary64))))))))))) Initial program 96.4%
Taylor expanded in x around inf
Simplified98.3%
*-inverses98.3
Applied egg-rr98.3%
Final simplification82.6%
(FPCore (x y z t a b c)
:precision binary64
(let* ((t_1 (* (+ a 0.8333333333333334) (+ a 0.8333333333333334))))
(if (<=
(/
x
(+
x
(*
y
(exp
(*
2.0
(+
(/ (* z (sqrt (+ t a))) t)
(* (- b c) (- (/ 2.0 (* t 3.0)) (+ a (/ 5.0 6.0))))))))))
0.9999999999999916)
(/
x
(fma
y
(fma
c
(fma
c
(fma
(* c 1.3333333333333333)
(* (+ a 0.8333333333333334) t_1)
(* 2.0 t_1))
(+ 1.6666666666666667 (* 2.0 a)))
1.0)
x))
1.0)))
double code(double x, double y, double z, double t, double a, double b, double c) {
double t_1 = (a + 0.8333333333333334) * (a + 0.8333333333333334);
double tmp;
if ((x / (x + (y * exp((2.0 * (((z * sqrt((t + a))) / t) + ((b - c) * ((2.0 / (t * 3.0)) - (a + (5.0 / 6.0)))))))))) <= 0.9999999999999916) {
tmp = x / fma(y, fma(c, fma(c, fma((c * 1.3333333333333333), ((a + 0.8333333333333334) * t_1), (2.0 * t_1)), (1.6666666666666667 + (2.0 * a))), 1.0), x);
} else {
tmp = 1.0;
}
return tmp;
}
function code(x, y, z, t, a, b, c) t_1 = Float64(Float64(a + 0.8333333333333334) * Float64(a + 0.8333333333333334)) tmp = 0.0 if (Float64(x / Float64(x + Float64(y * exp(Float64(2.0 * Float64(Float64(Float64(z * sqrt(Float64(t + a))) / t) + Float64(Float64(b - c) * Float64(Float64(2.0 / Float64(t * 3.0)) - Float64(a + Float64(5.0 / 6.0)))))))))) <= 0.9999999999999916) tmp = Float64(x / fma(y, fma(c, fma(c, fma(Float64(c * 1.3333333333333333), Float64(Float64(a + 0.8333333333333334) * t_1), Float64(2.0 * t_1)), Float64(1.6666666666666667 + Float64(2.0 * a))), 1.0), x)); else tmp = 1.0; end return tmp end
code[x_, y_, z_, t_, a_, b_, c_] := Block[{t$95$1 = N[(N[(a + 0.8333333333333334), $MachinePrecision] * N[(a + 0.8333333333333334), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[N[(x / N[(x + N[(y * N[Exp[N[(2.0 * N[(N[(N[(z * N[Sqrt[N[(t + a), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] / t), $MachinePrecision] + N[(N[(b - c), $MachinePrecision] * N[(N[(2.0 / N[(t * 3.0), $MachinePrecision]), $MachinePrecision] - N[(a + N[(5.0 / 6.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 0.9999999999999916], N[(x / N[(y * N[(c * N[(c * N[(N[(c * 1.3333333333333333), $MachinePrecision] * N[(N[(a + 0.8333333333333334), $MachinePrecision] * t$95$1), $MachinePrecision] + N[(2.0 * t$95$1), $MachinePrecision]), $MachinePrecision] + N[(1.6666666666666667 + N[(2.0 * a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision] + x), $MachinePrecision]), $MachinePrecision], 1.0]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \left(a + 0.8333333333333334\right) \cdot \left(a + 0.8333333333333334\right)\\
\mathbf{if}\;\frac{x}{x + y \cdot e^{2 \cdot \left(\frac{z \cdot \sqrt{t + a}}{t} + \left(b - c\right) \cdot \left(\frac{2}{t \cdot 3} - \left(a + \frac{5}{6}\right)\right)\right)}} \leq 0.9999999999999916:\\
\;\;\;\;\frac{x}{\mathsf{fma}\left(y, \mathsf{fma}\left(c, \mathsf{fma}\left(c, \mathsf{fma}\left(c \cdot 1.3333333333333333, \left(a + 0.8333333333333334\right) \cdot t\_1, 2 \cdot t\_1\right), 1.6666666666666667 + 2 \cdot a\right), 1\right), x\right)}\\
\mathbf{else}:\\
\;\;\;\;1\\
\end{array}
\end{array}
if (/.f64 x (+.f64 x (*.f64 y (exp.f64 (*.f64 #s(literal 2 binary64) (-.f64 (/.f64 (*.f64 z (sqrt.f64 (+.f64 t a))) t) (*.f64 (-.f64 b c) (-.f64 (+.f64 a (/.f64 #s(literal 5 binary64) #s(literal 6 binary64))) (/.f64 #s(literal 2 binary64) (*.f64 t #s(literal 3 binary64))))))))))) < 0.99999999999999156Initial program 99.3%
Taylor expanded in c around inf
*-lowering-*.f64N/A
+-commutativeN/A
associate--l+N/A
+-lowering-+.f64N/A
sub-negN/A
+-lowering-+.f64N/A
associate-*r/N/A
metadata-evalN/A
distribute-neg-fracN/A
metadata-evalN/A
/-lowering-/.f6463.4
Simplified63.4%
Taylor expanded in t around inf
+-commutativeN/A
accelerator-lowering-fma.f64N/A
exp-lowering-exp.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f6454.3
Simplified54.3%
Taylor expanded in c around 0
+-commutativeN/A
accelerator-lowering-fma.f64N/A
Simplified60.8%
if 0.99999999999999156 < (/.f64 x (+.f64 x (*.f64 y (exp.f64 (*.f64 #s(literal 2 binary64) (-.f64 (/.f64 (*.f64 z (sqrt.f64 (+.f64 t a))) t) (*.f64 (-.f64 b c) (-.f64 (+.f64 a (/.f64 #s(literal 5 binary64) #s(literal 6 binary64))) (/.f64 #s(literal 2 binary64) (*.f64 t #s(literal 3 binary64))))))))))) Initial program 96.4%
Taylor expanded in x around inf
Simplified98.3%
*-inverses98.3
Applied egg-rr98.3%
Final simplification77.1%
(FPCore (x y z t a b c)
:precision binary64
(if (<=
(/
x
(+
x
(*
y
(exp
(*
2.0
(+
(/ (* z (sqrt (+ t a))) t)
(* (- b c) (- (/ 2.0 (* t 3.0)) (+ a (/ 5.0 6.0))))))))))
0.9999999999999916)
(/
x
(fma
y
(fma
c
(*
2.0
(fma
c
(* (+ a 0.8333333333333334) (+ a 0.8333333333333334))
(+ a 0.8333333333333334)))
1.0)
x))
1.0))
double code(double x, double y, double z, double t, double a, double b, double c) {
double tmp;
if ((x / (x + (y * exp((2.0 * (((z * sqrt((t + a))) / t) + ((b - c) * ((2.0 / (t * 3.0)) - (a + (5.0 / 6.0)))))))))) <= 0.9999999999999916) {
tmp = x / fma(y, fma(c, (2.0 * fma(c, ((a + 0.8333333333333334) * (a + 0.8333333333333334)), (a + 0.8333333333333334))), 1.0), x);
} else {
tmp = 1.0;
}
return tmp;
}
function code(x, y, z, t, a, b, c) tmp = 0.0 if (Float64(x / Float64(x + Float64(y * exp(Float64(2.0 * Float64(Float64(Float64(z * sqrt(Float64(t + a))) / t) + Float64(Float64(b - c) * Float64(Float64(2.0 / Float64(t * 3.0)) - Float64(a + Float64(5.0 / 6.0)))))))))) <= 0.9999999999999916) tmp = Float64(x / fma(y, fma(c, Float64(2.0 * fma(c, Float64(Float64(a + 0.8333333333333334) * Float64(a + 0.8333333333333334)), Float64(a + 0.8333333333333334))), 1.0), x)); else tmp = 1.0; end return tmp end
code[x_, y_, z_, t_, a_, b_, c_] := If[LessEqual[N[(x / N[(x + N[(y * N[Exp[N[(2.0 * N[(N[(N[(z * N[Sqrt[N[(t + a), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] / t), $MachinePrecision] + N[(N[(b - c), $MachinePrecision] * N[(N[(2.0 / N[(t * 3.0), $MachinePrecision]), $MachinePrecision] - N[(a + N[(5.0 / 6.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 0.9999999999999916], N[(x / N[(y * N[(c * N[(2.0 * N[(c * N[(N[(a + 0.8333333333333334), $MachinePrecision] * N[(a + 0.8333333333333334), $MachinePrecision]), $MachinePrecision] + N[(a + 0.8333333333333334), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision] + x), $MachinePrecision]), $MachinePrecision], 1.0]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\frac{x}{x + y \cdot e^{2 \cdot \left(\frac{z \cdot \sqrt{t + a}}{t} + \left(b - c\right) \cdot \left(\frac{2}{t \cdot 3} - \left(a + \frac{5}{6}\right)\right)\right)}} \leq 0.9999999999999916:\\
\;\;\;\;\frac{x}{\mathsf{fma}\left(y, \mathsf{fma}\left(c, 2 \cdot \mathsf{fma}\left(c, \left(a + 0.8333333333333334\right) \cdot \left(a + 0.8333333333333334\right), a + 0.8333333333333334\right), 1\right), x\right)}\\
\mathbf{else}:\\
\;\;\;\;1\\
\end{array}
\end{array}
if (/.f64 x (+.f64 x (*.f64 y (exp.f64 (*.f64 #s(literal 2 binary64) (-.f64 (/.f64 (*.f64 z (sqrt.f64 (+.f64 t a))) t) (*.f64 (-.f64 b c) (-.f64 (+.f64 a (/.f64 #s(literal 5 binary64) #s(literal 6 binary64))) (/.f64 #s(literal 2 binary64) (*.f64 t #s(literal 3 binary64))))))))))) < 0.99999999999999156Initial program 99.3%
Taylor expanded in c around inf
*-lowering-*.f64N/A
+-commutativeN/A
associate--l+N/A
+-lowering-+.f64N/A
sub-negN/A
+-lowering-+.f64N/A
associate-*r/N/A
metadata-evalN/A
distribute-neg-fracN/A
metadata-evalN/A
/-lowering-/.f6463.4
Simplified63.4%
Taylor expanded in t around inf
+-commutativeN/A
accelerator-lowering-fma.f64N/A
exp-lowering-exp.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f6454.3
Simplified54.3%
Taylor expanded in c around 0
+-commutativeN/A
accelerator-lowering-fma.f64N/A
distribute-lft-outN/A
*-lowering-*.f64N/A
accelerator-lowering-fma.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f6459.8
Simplified59.8%
if 0.99999999999999156 < (/.f64 x (+.f64 x (*.f64 y (exp.f64 (*.f64 #s(literal 2 binary64) (-.f64 (/.f64 (*.f64 z (sqrt.f64 (+.f64 t a))) t) (*.f64 (-.f64 b c) (-.f64 (+.f64 a (/.f64 #s(literal 5 binary64) #s(literal 6 binary64))) (/.f64 #s(literal 2 binary64) (*.f64 t #s(literal 3 binary64))))))))))) Initial program 96.4%
Taylor expanded in x around inf
Simplified98.3%
*-inverses98.3
Applied egg-rr98.3%
Final simplification76.4%
(FPCore (x y z t a b c)
:precision binary64
(if (<=
(/
x
(+
x
(*
y
(exp
(*
2.0
(+
(/ (* z (sqrt (+ t a))) t)
(* (- b c) (- (/ 2.0 (* t 3.0)) (+ a (/ 5.0 6.0))))))))))
5e-24)
(* (- y x) (/ x (* (+ x y) (- y x))))
1.0))
double code(double x, double y, double z, double t, double a, double b, double c) {
double tmp;
if ((x / (x + (y * exp((2.0 * (((z * sqrt((t + a))) / t) + ((b - c) * ((2.0 / (t * 3.0)) - (a + (5.0 / 6.0)))))))))) <= 5e-24) {
tmp = (y - x) * (x / ((x + y) * (y - x)));
} else {
tmp = 1.0;
}
return tmp;
}
real(8) function code(x, y, z, t, a, b, c)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8), intent (in) :: c
real(8) :: tmp
if ((x / (x + (y * exp((2.0d0 * (((z * sqrt((t + a))) / t) + ((b - c) * ((2.0d0 / (t * 3.0d0)) - (a + (5.0d0 / 6.0d0)))))))))) <= 5d-24) then
tmp = (y - x) * (x / ((x + y) * (y - x)))
else
tmp = 1.0d0
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b, double c) {
double tmp;
if ((x / (x + (y * Math.exp((2.0 * (((z * Math.sqrt((t + a))) / t) + ((b - c) * ((2.0 / (t * 3.0)) - (a + (5.0 / 6.0)))))))))) <= 5e-24) {
tmp = (y - x) * (x / ((x + y) * (y - x)));
} else {
tmp = 1.0;
}
return tmp;
}
def code(x, y, z, t, a, b, c): tmp = 0 if (x / (x + (y * math.exp((2.0 * (((z * math.sqrt((t + a))) / t) + ((b - c) * ((2.0 / (t * 3.0)) - (a + (5.0 / 6.0)))))))))) <= 5e-24: tmp = (y - x) * (x / ((x + y) * (y - x))) else: tmp = 1.0 return tmp
function code(x, y, z, t, a, b, c) tmp = 0.0 if (Float64(x / Float64(x + Float64(y * exp(Float64(2.0 * Float64(Float64(Float64(z * sqrt(Float64(t + a))) / t) + Float64(Float64(b - c) * Float64(Float64(2.0 / Float64(t * 3.0)) - Float64(a + Float64(5.0 / 6.0)))))))))) <= 5e-24) tmp = Float64(Float64(y - x) * Float64(x / Float64(Float64(x + y) * Float64(y - x)))); else tmp = 1.0; end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c) tmp = 0.0; if ((x / (x + (y * exp((2.0 * (((z * sqrt((t + a))) / t) + ((b - c) * ((2.0 / (t * 3.0)) - (a + (5.0 / 6.0)))))))))) <= 5e-24) tmp = (y - x) * (x / ((x + y) * (y - x))); else tmp = 1.0; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_] := If[LessEqual[N[(x / N[(x + N[(y * N[Exp[N[(2.0 * N[(N[(N[(z * N[Sqrt[N[(t + a), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] / t), $MachinePrecision] + N[(N[(b - c), $MachinePrecision] * N[(N[(2.0 / N[(t * 3.0), $MachinePrecision]), $MachinePrecision] - N[(a + N[(5.0 / 6.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 5e-24], N[(N[(y - x), $MachinePrecision] * N[(x / N[(N[(x + y), $MachinePrecision] * N[(y - x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 1.0]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\frac{x}{x + y \cdot e^{2 \cdot \left(\frac{z \cdot \sqrt{t + a}}{t} + \left(b - c\right) \cdot \left(\frac{2}{t \cdot 3} - \left(a + \frac{5}{6}\right)\right)\right)}} \leq 5 \cdot 10^{-24}:\\
\;\;\;\;\left(y - x\right) \cdot \frac{x}{\left(x + y\right) \cdot \left(y - x\right)}\\
\mathbf{else}:\\
\;\;\;\;1\\
\end{array}
\end{array}
if (/.f64 x (+.f64 x (*.f64 y (exp.f64 (*.f64 #s(literal 2 binary64) (-.f64 (/.f64 (*.f64 z (sqrt.f64 (+.f64 t a))) t) (*.f64 (-.f64 b c) (-.f64 (+.f64 a (/.f64 #s(literal 5 binary64) #s(literal 6 binary64))) (/.f64 #s(literal 2 binary64) (*.f64 t #s(literal 3 binary64))))))))))) < 4.9999999999999998e-24Initial program 99.3%
Taylor expanded in b around inf
*-lowering-*.f64N/A
--lowering--.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f64N/A
+-commutativeN/A
+-lowering-+.f6462.8
Simplified62.8%
clear-numN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
+-commutativeN/A
accelerator-lowering-fma.f64N/A
exp-lowering-exp.f64N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
+-commutativeN/A
associate--r+N/A
--lowering--.f64N/A
--lowering--.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f6462.3
Applied egg-rr62.3%
Taylor expanded in b around 0
/-lowering-/.f64N/A
+-commutativeN/A
+-lowering-+.f6418.1
Simplified18.1%
flip-+N/A
associate-/r/N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
difference-of-squaresN/A
*-lowering-*.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
--lowering--.f64N/A
--lowering--.f6449.4
Applied egg-rr49.4%
if 4.9999999999999998e-24 < (/.f64 x (+.f64 x (*.f64 y (exp.f64 (*.f64 #s(literal 2 binary64) (-.f64 (/.f64 (*.f64 z (sqrt.f64 (+.f64 t a))) t) (*.f64 (-.f64 b c) (-.f64 (+.f64 a (/.f64 #s(literal 5 binary64) #s(literal 6 binary64))) (/.f64 #s(literal 2 binary64) (*.f64 t #s(literal 3 binary64))))))))))) Initial program 96.4%
Taylor expanded in x around inf
Simplified98.2%
*-inverses98.2
Applied egg-rr98.2%
Final simplification70.8%
(FPCore (x y z t a b c)
:precision binary64
(if (<=
(/
x
(+
x
(*
y
(exp
(*
2.0
(+
(/ (* z (sqrt (+ t a))) t)
(* (- b c) (- (/ 2.0 (* t 3.0)) (+ a (/ 5.0 6.0))))))))))
0.9999999999999916)
(/ x (fma y (fma c (fma c 1.3888888888888888 1.6666666666666667) 1.0) x))
1.0))
double code(double x, double y, double z, double t, double a, double b, double c) {
double tmp;
if ((x / (x + (y * exp((2.0 * (((z * sqrt((t + a))) / t) + ((b - c) * ((2.0 / (t * 3.0)) - (a + (5.0 / 6.0)))))))))) <= 0.9999999999999916) {
tmp = x / fma(y, fma(c, fma(c, 1.3888888888888888, 1.6666666666666667), 1.0), x);
} else {
tmp = 1.0;
}
return tmp;
}
function code(x, y, z, t, a, b, c) tmp = 0.0 if (Float64(x / Float64(x + Float64(y * exp(Float64(2.0 * Float64(Float64(Float64(z * sqrt(Float64(t + a))) / t) + Float64(Float64(b - c) * Float64(Float64(2.0 / Float64(t * 3.0)) - Float64(a + Float64(5.0 / 6.0)))))))))) <= 0.9999999999999916) tmp = Float64(x / fma(y, fma(c, fma(c, 1.3888888888888888, 1.6666666666666667), 1.0), x)); else tmp = 1.0; end return tmp end
code[x_, y_, z_, t_, a_, b_, c_] := If[LessEqual[N[(x / N[(x + N[(y * N[Exp[N[(2.0 * N[(N[(N[(z * N[Sqrt[N[(t + a), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] / t), $MachinePrecision] + N[(N[(b - c), $MachinePrecision] * N[(N[(2.0 / N[(t * 3.0), $MachinePrecision]), $MachinePrecision] - N[(a + N[(5.0 / 6.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 0.9999999999999916], N[(x / N[(y * N[(c * N[(c * 1.3888888888888888 + 1.6666666666666667), $MachinePrecision] + 1.0), $MachinePrecision] + x), $MachinePrecision]), $MachinePrecision], 1.0]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\frac{x}{x + y \cdot e^{2 \cdot \left(\frac{z \cdot \sqrt{t + a}}{t} + \left(b - c\right) \cdot \left(\frac{2}{t \cdot 3} - \left(a + \frac{5}{6}\right)\right)\right)}} \leq 0.9999999999999916:\\
\;\;\;\;\frac{x}{\mathsf{fma}\left(y, \mathsf{fma}\left(c, \mathsf{fma}\left(c, 1.3888888888888888, 1.6666666666666667\right), 1\right), x\right)}\\
\mathbf{else}:\\
\;\;\;\;1\\
\end{array}
\end{array}
if (/.f64 x (+.f64 x (*.f64 y (exp.f64 (*.f64 #s(literal 2 binary64) (-.f64 (/.f64 (*.f64 z (sqrt.f64 (+.f64 t a))) t) (*.f64 (-.f64 b c) (-.f64 (+.f64 a (/.f64 #s(literal 5 binary64) #s(literal 6 binary64))) (/.f64 #s(literal 2 binary64) (*.f64 t #s(literal 3 binary64))))))))))) < 0.99999999999999156Initial program 99.3%
Taylor expanded in c around inf
*-lowering-*.f64N/A
+-commutativeN/A
associate--l+N/A
+-lowering-+.f64N/A
sub-negN/A
+-lowering-+.f64N/A
associate-*r/N/A
metadata-evalN/A
distribute-neg-fracN/A
metadata-evalN/A
/-lowering-/.f6463.4
Simplified63.4%
Taylor expanded in t around inf
+-commutativeN/A
accelerator-lowering-fma.f64N/A
exp-lowering-exp.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f6454.3
Simplified54.3%
Taylor expanded in a around 0
*-lowering-*.f6449.0
Simplified49.0%
Taylor expanded in c around 0
+-commutativeN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f6444.5
Simplified44.5%
if 0.99999999999999156 < (/.f64 x (+.f64 x (*.f64 y (exp.f64 (*.f64 #s(literal 2 binary64) (-.f64 (/.f64 (*.f64 z (sqrt.f64 (+.f64 t a))) t) (*.f64 (-.f64 b c) (-.f64 (+.f64 a (/.f64 #s(literal 5 binary64) #s(literal 6 binary64))) (/.f64 #s(literal 2 binary64) (*.f64 t #s(literal 3 binary64))))))))))) Initial program 96.4%
Taylor expanded in x around inf
Simplified98.3%
*-inverses98.3
Applied egg-rr98.3%
Final simplification67.8%
(FPCore (x y z t a b c)
:precision binary64
(if (<=
(/
x
(+
x
(*
y
(exp
(*
2.0
(+
(/ (* z (sqrt (+ t a))) t)
(* (- b c) (- (/ 2.0 (* t 3.0)) (+ a (/ 5.0 6.0))))))))))
0.9999999999999916)
(/ x (fma y (fma 1.6666666666666667 c 1.0) x))
1.0))
double code(double x, double y, double z, double t, double a, double b, double c) {
double tmp;
if ((x / (x + (y * exp((2.0 * (((z * sqrt((t + a))) / t) + ((b - c) * ((2.0 / (t * 3.0)) - (a + (5.0 / 6.0)))))))))) <= 0.9999999999999916) {
tmp = x / fma(y, fma(1.6666666666666667, c, 1.0), x);
} else {
tmp = 1.0;
}
return tmp;
}
function code(x, y, z, t, a, b, c) tmp = 0.0 if (Float64(x / Float64(x + Float64(y * exp(Float64(2.0 * Float64(Float64(Float64(z * sqrt(Float64(t + a))) / t) + Float64(Float64(b - c) * Float64(Float64(2.0 / Float64(t * 3.0)) - Float64(a + Float64(5.0 / 6.0)))))))))) <= 0.9999999999999916) tmp = Float64(x / fma(y, fma(1.6666666666666667, c, 1.0), x)); else tmp = 1.0; end return tmp end
code[x_, y_, z_, t_, a_, b_, c_] := If[LessEqual[N[(x / N[(x + N[(y * N[Exp[N[(2.0 * N[(N[(N[(z * N[Sqrt[N[(t + a), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] / t), $MachinePrecision] + N[(N[(b - c), $MachinePrecision] * N[(N[(2.0 / N[(t * 3.0), $MachinePrecision]), $MachinePrecision] - N[(a + N[(5.0 / 6.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 0.9999999999999916], N[(x / N[(y * N[(1.6666666666666667 * c + 1.0), $MachinePrecision] + x), $MachinePrecision]), $MachinePrecision], 1.0]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\frac{x}{x + y \cdot e^{2 \cdot \left(\frac{z \cdot \sqrt{t + a}}{t} + \left(b - c\right) \cdot \left(\frac{2}{t \cdot 3} - \left(a + \frac{5}{6}\right)\right)\right)}} \leq 0.9999999999999916:\\
\;\;\;\;\frac{x}{\mathsf{fma}\left(y, \mathsf{fma}\left(1.6666666666666667, c, 1\right), x\right)}\\
\mathbf{else}:\\
\;\;\;\;1\\
\end{array}
\end{array}
if (/.f64 x (+.f64 x (*.f64 y (exp.f64 (*.f64 #s(literal 2 binary64) (-.f64 (/.f64 (*.f64 z (sqrt.f64 (+.f64 t a))) t) (*.f64 (-.f64 b c) (-.f64 (+.f64 a (/.f64 #s(literal 5 binary64) #s(literal 6 binary64))) (/.f64 #s(literal 2 binary64) (*.f64 t #s(literal 3 binary64))))))))))) < 0.99999999999999156Initial program 99.3%
Taylor expanded in c around inf
*-lowering-*.f64N/A
+-commutativeN/A
associate--l+N/A
+-lowering-+.f64N/A
sub-negN/A
+-lowering-+.f64N/A
associate-*r/N/A
metadata-evalN/A
distribute-neg-fracN/A
metadata-evalN/A
/-lowering-/.f6463.4
Simplified63.4%
Taylor expanded in t around inf
+-commutativeN/A
accelerator-lowering-fma.f64N/A
exp-lowering-exp.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f6454.3
Simplified54.3%
Taylor expanded in a around 0
*-lowering-*.f6449.0
Simplified49.0%
Taylor expanded in c around 0
+-commutativeN/A
accelerator-lowering-fma.f6431.0
Simplified31.0%
if 0.99999999999999156 < (/.f64 x (+.f64 x (*.f64 y (exp.f64 (*.f64 #s(literal 2 binary64) (-.f64 (/.f64 (*.f64 z (sqrt.f64 (+.f64 t a))) t) (*.f64 (-.f64 b c) (-.f64 (+.f64 a (/.f64 #s(literal 5 binary64) #s(literal 6 binary64))) (/.f64 #s(literal 2 binary64) (*.f64 t #s(literal 3 binary64))))))))))) Initial program 96.4%
Taylor expanded in x around inf
Simplified98.3%
*-inverses98.3
Applied egg-rr98.3%
Final simplification60.1%
(FPCore (x y z t a b c)
:precision binary64
(if (<=
(/
x
(+
x
(*
y
(exp
(*
2.0
(+
(/ (* z (sqrt (+ t a))) t)
(* (- b c) (- (/ 2.0 (* t 3.0)) (+ a (/ 5.0 6.0))))))))))
0.9999999999999916)
(/ x (+ x y))
1.0))
double code(double x, double y, double z, double t, double a, double b, double c) {
double tmp;
if ((x / (x + (y * exp((2.0 * (((z * sqrt((t + a))) / t) + ((b - c) * ((2.0 / (t * 3.0)) - (a + (5.0 / 6.0)))))))))) <= 0.9999999999999916) {
tmp = x / (x + y);
} else {
tmp = 1.0;
}
return tmp;
}
real(8) function code(x, y, z, t, a, b, c)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8), intent (in) :: c
real(8) :: tmp
if ((x / (x + (y * exp((2.0d0 * (((z * sqrt((t + a))) / t) + ((b - c) * ((2.0d0 / (t * 3.0d0)) - (a + (5.0d0 / 6.0d0)))))))))) <= 0.9999999999999916d0) then
tmp = x / (x + y)
else
tmp = 1.0d0
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b, double c) {
double tmp;
if ((x / (x + (y * Math.exp((2.0 * (((z * Math.sqrt((t + a))) / t) + ((b - c) * ((2.0 / (t * 3.0)) - (a + (5.0 / 6.0)))))))))) <= 0.9999999999999916) {
tmp = x / (x + y);
} else {
tmp = 1.0;
}
return tmp;
}
def code(x, y, z, t, a, b, c): tmp = 0 if (x / (x + (y * math.exp((2.0 * (((z * math.sqrt((t + a))) / t) + ((b - c) * ((2.0 / (t * 3.0)) - (a + (5.0 / 6.0)))))))))) <= 0.9999999999999916: tmp = x / (x + y) else: tmp = 1.0 return tmp
function code(x, y, z, t, a, b, c) tmp = 0.0 if (Float64(x / Float64(x + Float64(y * exp(Float64(2.0 * Float64(Float64(Float64(z * sqrt(Float64(t + a))) / t) + Float64(Float64(b - c) * Float64(Float64(2.0 / Float64(t * 3.0)) - Float64(a + Float64(5.0 / 6.0)))))))))) <= 0.9999999999999916) tmp = Float64(x / Float64(x + y)); else tmp = 1.0; end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c) tmp = 0.0; if ((x / (x + (y * exp((2.0 * (((z * sqrt((t + a))) / t) + ((b - c) * ((2.0 / (t * 3.0)) - (a + (5.0 / 6.0)))))))))) <= 0.9999999999999916) tmp = x / (x + y); else tmp = 1.0; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_] := If[LessEqual[N[(x / N[(x + N[(y * N[Exp[N[(2.0 * N[(N[(N[(z * N[Sqrt[N[(t + a), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] / t), $MachinePrecision] + N[(N[(b - c), $MachinePrecision] * N[(N[(2.0 / N[(t * 3.0), $MachinePrecision]), $MachinePrecision] - N[(a + N[(5.0 / 6.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 0.9999999999999916], N[(x / N[(x + y), $MachinePrecision]), $MachinePrecision], 1.0]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\frac{x}{x + y \cdot e^{2 \cdot \left(\frac{z \cdot \sqrt{t + a}}{t} + \left(b - c\right) \cdot \left(\frac{2}{t \cdot 3} - \left(a + \frac{5}{6}\right)\right)\right)}} \leq 0.9999999999999916:\\
\;\;\;\;\frac{x}{x + y}\\
\mathbf{else}:\\
\;\;\;\;1\\
\end{array}
\end{array}
if (/.f64 x (+.f64 x (*.f64 y (exp.f64 (*.f64 #s(literal 2 binary64) (-.f64 (/.f64 (*.f64 z (sqrt.f64 (+.f64 t a))) t) (*.f64 (-.f64 b c) (-.f64 (+.f64 a (/.f64 #s(literal 5 binary64) #s(literal 6 binary64))) (/.f64 #s(literal 2 binary64) (*.f64 t #s(literal 3 binary64))))))))))) < 0.99999999999999156Initial program 99.3%
Taylor expanded in b around inf
*-lowering-*.f64N/A
--lowering--.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f64N/A
+-commutativeN/A
+-lowering-+.f6463.0
Simplified63.0%
Taylor expanded in b around 0
+-lowering-+.f6418.6
Simplified18.6%
if 0.99999999999999156 < (/.f64 x (+.f64 x (*.f64 y (exp.f64 (*.f64 #s(literal 2 binary64) (-.f64 (/.f64 (*.f64 z (sqrt.f64 (+.f64 t a))) t) (*.f64 (-.f64 b c) (-.f64 (+.f64 a (/.f64 #s(literal 5 binary64) #s(literal 6 binary64))) (/.f64 #s(literal 2 binary64) (*.f64 t #s(literal 3 binary64))))))))))) Initial program 96.4%
Taylor expanded in x around inf
Simplified98.3%
*-inverses98.3
Applied egg-rr98.3%
Final simplification53.2%
(FPCore (x y z t a b c)
:precision binary64
(if (<=
(/
x
(+
x
(*
y
(exp
(*
2.0
(+
(/ (* z (sqrt (+ t a))) t)
(* (- b c) (- (/ 2.0 (* t 3.0)) (+ a (/ 5.0 6.0))))))))))
5e-24)
(/ x y)
1.0))
double code(double x, double y, double z, double t, double a, double b, double c) {
double tmp;
if ((x / (x + (y * exp((2.0 * (((z * sqrt((t + a))) / t) + ((b - c) * ((2.0 / (t * 3.0)) - (a + (5.0 / 6.0)))))))))) <= 5e-24) {
tmp = x / y;
} else {
tmp = 1.0;
}
return tmp;
}
real(8) function code(x, y, z, t, a, b, c)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8), intent (in) :: c
real(8) :: tmp
if ((x / (x + (y * exp((2.0d0 * (((z * sqrt((t + a))) / t) + ((b - c) * ((2.0d0 / (t * 3.0d0)) - (a + (5.0d0 / 6.0d0)))))))))) <= 5d-24) then
tmp = x / y
else
tmp = 1.0d0
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b, double c) {
double tmp;
if ((x / (x + (y * Math.exp((2.0 * (((z * Math.sqrt((t + a))) / t) + ((b - c) * ((2.0 / (t * 3.0)) - (a + (5.0 / 6.0)))))))))) <= 5e-24) {
tmp = x / y;
} else {
tmp = 1.0;
}
return tmp;
}
def code(x, y, z, t, a, b, c): tmp = 0 if (x / (x + (y * math.exp((2.0 * (((z * math.sqrt((t + a))) / t) + ((b - c) * ((2.0 / (t * 3.0)) - (a + (5.0 / 6.0)))))))))) <= 5e-24: tmp = x / y else: tmp = 1.0 return tmp
function code(x, y, z, t, a, b, c) tmp = 0.0 if (Float64(x / Float64(x + Float64(y * exp(Float64(2.0 * Float64(Float64(Float64(z * sqrt(Float64(t + a))) / t) + Float64(Float64(b - c) * Float64(Float64(2.0 / Float64(t * 3.0)) - Float64(a + Float64(5.0 / 6.0)))))))))) <= 5e-24) tmp = Float64(x / y); else tmp = 1.0; end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c) tmp = 0.0; if ((x / (x + (y * exp((2.0 * (((z * sqrt((t + a))) / t) + ((b - c) * ((2.0 / (t * 3.0)) - (a + (5.0 / 6.0)))))))))) <= 5e-24) tmp = x / y; else tmp = 1.0; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_] := If[LessEqual[N[(x / N[(x + N[(y * N[Exp[N[(2.0 * N[(N[(N[(z * N[Sqrt[N[(t + a), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] / t), $MachinePrecision] + N[(N[(b - c), $MachinePrecision] * N[(N[(2.0 / N[(t * 3.0), $MachinePrecision]), $MachinePrecision] - N[(a + N[(5.0 / 6.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 5e-24], N[(x / y), $MachinePrecision], 1.0]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\frac{x}{x + y \cdot e^{2 \cdot \left(\frac{z \cdot \sqrt{t + a}}{t} + \left(b - c\right) \cdot \left(\frac{2}{t \cdot 3} - \left(a + \frac{5}{6}\right)\right)\right)}} \leq 5 \cdot 10^{-24}:\\
\;\;\;\;\frac{x}{y}\\
\mathbf{else}:\\
\;\;\;\;1\\
\end{array}
\end{array}
if (/.f64 x (+.f64 x (*.f64 y (exp.f64 (*.f64 #s(literal 2 binary64) (-.f64 (/.f64 (*.f64 z (sqrt.f64 (+.f64 t a))) t) (*.f64 (-.f64 b c) (-.f64 (+.f64 a (/.f64 #s(literal 5 binary64) #s(literal 6 binary64))) (/.f64 #s(literal 2 binary64) (*.f64 t #s(literal 3 binary64))))))))))) < 4.9999999999999998e-24Initial program 99.3%
Taylor expanded in b around inf
*-lowering-*.f64N/A
--lowering--.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f64N/A
+-commutativeN/A
+-lowering-+.f6462.8
Simplified62.8%
clear-numN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
+-commutativeN/A
accelerator-lowering-fma.f64N/A
exp-lowering-exp.f64N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
+-commutativeN/A
associate--r+N/A
--lowering--.f64N/A
--lowering--.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f6462.3
Applied egg-rr62.3%
Taylor expanded in b around 0
/-lowering-/.f64N/A
+-commutativeN/A
+-lowering-+.f6418.1
Simplified18.1%
Taylor expanded in x around 0
/-lowering-/.f6416.8
Simplified16.8%
if 4.9999999999999998e-24 < (/.f64 x (+.f64 x (*.f64 y (exp.f64 (*.f64 #s(literal 2 binary64) (-.f64 (/.f64 (*.f64 z (sqrt.f64 (+.f64 t a))) t) (*.f64 (-.f64 b c) (-.f64 (+.f64 a (/.f64 #s(literal 5 binary64) #s(literal 6 binary64))) (/.f64 #s(literal 2 binary64) (*.f64 t #s(literal 3 binary64))))))))))) Initial program 96.4%
Taylor expanded in x around inf
Simplified98.2%
*-inverses98.2
Applied egg-rr98.2%
Final simplification52.4%
(FPCore (x y z t a b c)
:precision binary64
(let* ((t_1
(+
(/ (* z (sqrt (+ t a))) t)
(* (- b c) (- (/ 2.0 (* t 3.0)) (+ a (/ 5.0 6.0)))))))
(if (<= t_1 -200.0)
1.0
(if (<= t_1 4e+288)
(/
x
(fma
y
(fma
c
(fma
c
(fma c 0.7716049382716049 1.3888888888888888)
1.6666666666666667)
1.0)
x))
(/
x
(+
x
(fma
(* 2.0 b)
(* (* (+ a 0.8333333333333334) (+ a 0.8333333333333334)) (* y b))
y)))))))
double code(double x, double y, double z, double t, double a, double b, double c) {
double t_1 = ((z * sqrt((t + a))) / t) + ((b - c) * ((2.0 / (t * 3.0)) - (a + (5.0 / 6.0))));
double tmp;
if (t_1 <= -200.0) {
tmp = 1.0;
} else if (t_1 <= 4e+288) {
tmp = x / fma(y, fma(c, fma(c, fma(c, 0.7716049382716049, 1.3888888888888888), 1.6666666666666667), 1.0), x);
} else {
tmp = x / (x + fma((2.0 * b), (((a + 0.8333333333333334) * (a + 0.8333333333333334)) * (y * b)), y));
}
return tmp;
}
function code(x, y, z, t, a, b, c) t_1 = Float64(Float64(Float64(z * sqrt(Float64(t + a))) / t) + Float64(Float64(b - c) * Float64(Float64(2.0 / Float64(t * 3.0)) - Float64(a + Float64(5.0 / 6.0))))) tmp = 0.0 if (t_1 <= -200.0) tmp = 1.0; elseif (t_1 <= 4e+288) tmp = Float64(x / fma(y, fma(c, fma(c, fma(c, 0.7716049382716049, 1.3888888888888888), 1.6666666666666667), 1.0), x)); else tmp = Float64(x / Float64(x + fma(Float64(2.0 * b), Float64(Float64(Float64(a + 0.8333333333333334) * Float64(a + 0.8333333333333334)) * Float64(y * b)), y))); end return tmp end
code[x_, y_, z_, t_, a_, b_, c_] := Block[{t$95$1 = N[(N[(N[(z * N[Sqrt[N[(t + a), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] / t), $MachinePrecision] + N[(N[(b - c), $MachinePrecision] * N[(N[(2.0 / N[(t * 3.0), $MachinePrecision]), $MachinePrecision] - N[(a + N[(5.0 / 6.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$1, -200.0], 1.0, If[LessEqual[t$95$1, 4e+288], N[(x / N[(y * N[(c * N[(c * N[(c * 0.7716049382716049 + 1.3888888888888888), $MachinePrecision] + 1.6666666666666667), $MachinePrecision] + 1.0), $MachinePrecision] + x), $MachinePrecision]), $MachinePrecision], N[(x / N[(x + N[(N[(2.0 * b), $MachinePrecision] * N[(N[(N[(a + 0.8333333333333334), $MachinePrecision] * N[(a + 0.8333333333333334), $MachinePrecision]), $MachinePrecision] * N[(y * b), $MachinePrecision]), $MachinePrecision] + y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{z \cdot \sqrt{t + a}}{t} + \left(b - c\right) \cdot \left(\frac{2}{t \cdot 3} - \left(a + \frac{5}{6}\right)\right)\\
\mathbf{if}\;t\_1 \leq -200:\\
\;\;\;\;1\\
\mathbf{elif}\;t\_1 \leq 4 \cdot 10^{+288}:\\
\;\;\;\;\frac{x}{\mathsf{fma}\left(y, \mathsf{fma}\left(c, \mathsf{fma}\left(c, \mathsf{fma}\left(c, 0.7716049382716049, 1.3888888888888888\right), 1.6666666666666667\right), 1\right), x\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{x}{x + \mathsf{fma}\left(2 \cdot b, \left(\left(a + 0.8333333333333334\right) \cdot \left(a + 0.8333333333333334\right)\right) \cdot \left(y \cdot b\right), y\right)}\\
\end{array}
\end{array}
if (-.f64 (/.f64 (*.f64 z (sqrt.f64 (+.f64 t a))) t) (*.f64 (-.f64 b c) (-.f64 (+.f64 a (/.f64 #s(literal 5 binary64) #s(literal 6 binary64))) (/.f64 #s(literal 2 binary64) (*.f64 t #s(literal 3 binary64)))))) < -200Initial program 99.1%
Taylor expanded in x around inf
Simplified99.1%
*-inverses99.1
Applied egg-rr99.1%
if -200 < (-.f64 (/.f64 (*.f64 z (sqrt.f64 (+.f64 t a))) t) (*.f64 (-.f64 b c) (-.f64 (+.f64 a (/.f64 #s(literal 5 binary64) #s(literal 6 binary64))) (/.f64 #s(literal 2 binary64) (*.f64 t #s(literal 3 binary64)))))) < 4e288Initial program 99.9%
Taylor expanded in c around inf
*-lowering-*.f64N/A
+-commutativeN/A
associate--l+N/A
+-lowering-+.f64N/A
sub-negN/A
+-lowering-+.f64N/A
associate-*r/N/A
metadata-evalN/A
distribute-neg-fracN/A
metadata-evalN/A
/-lowering-/.f6474.0
Simplified74.0%
Taylor expanded in t around inf
+-commutativeN/A
accelerator-lowering-fma.f64N/A
exp-lowering-exp.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f6464.4
Simplified64.4%
Taylor expanded in a around 0
*-lowering-*.f6458.2
Simplified58.2%
Taylor expanded in c around 0
+-commutativeN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f6459.6
Simplified59.6%
if 4e288 < (-.f64 (/.f64 (*.f64 z (sqrt.f64 (+.f64 t a))) t) (*.f64 (-.f64 b c) (-.f64 (+.f64 a (/.f64 #s(literal 5 binary64) #s(literal 6 binary64))) (/.f64 #s(literal 2 binary64) (*.f64 t #s(literal 3 binary64)))))) Initial program 94.8%
Taylor expanded in b around inf
*-lowering-*.f64N/A
--lowering--.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f64N/A
+-commutativeN/A
+-lowering-+.f6468.3
Simplified68.3%
Taylor expanded in b around 0
+-lowering-+.f64N/A
+-commutativeN/A
accelerator-lowering-fma.f64N/A
Simplified79.5%
Taylor expanded in t around inf
/-lowering-/.f64N/A
+-lowering-+.f64N/A
+-commutativeN/A
associate-*r*N/A
accelerator-lowering-fma.f64N/A
Simplified64.4%
Taylor expanded in b around inf
associate-*r*N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f6461.9
Simplified61.9%
Final simplification76.0%
(FPCore (x y z t a b c)
:precision binary64
(let* ((t_1
(+
(/ (* z (sqrt (+ t a))) t)
(* (- b c) (- (/ 2.0 (* t 3.0)) (+ a (/ 5.0 6.0)))))))
(if (<= t_1 -200.0)
1.0
(if (<= t_1 4e+288)
(/
x
(fma
y
(fma
c
(fma
c
(fma c 0.7716049382716049 1.3888888888888888)
1.6666666666666667)
1.0)
x))
(/ x (+ x (fma b (* (* b 0.8888888888888888) (/ y (* t t))) y)))))))
double code(double x, double y, double z, double t, double a, double b, double c) {
double t_1 = ((z * sqrt((t + a))) / t) + ((b - c) * ((2.0 / (t * 3.0)) - (a + (5.0 / 6.0))));
double tmp;
if (t_1 <= -200.0) {
tmp = 1.0;
} else if (t_1 <= 4e+288) {
tmp = x / fma(y, fma(c, fma(c, fma(c, 0.7716049382716049, 1.3888888888888888), 1.6666666666666667), 1.0), x);
} else {
tmp = x / (x + fma(b, ((b * 0.8888888888888888) * (y / (t * t))), y));
}
return tmp;
}
function code(x, y, z, t, a, b, c) t_1 = Float64(Float64(Float64(z * sqrt(Float64(t + a))) / t) + Float64(Float64(b - c) * Float64(Float64(2.0 / Float64(t * 3.0)) - Float64(a + Float64(5.0 / 6.0))))) tmp = 0.0 if (t_1 <= -200.0) tmp = 1.0; elseif (t_1 <= 4e+288) tmp = Float64(x / fma(y, fma(c, fma(c, fma(c, 0.7716049382716049, 1.3888888888888888), 1.6666666666666667), 1.0), x)); else tmp = Float64(x / Float64(x + fma(b, Float64(Float64(b * 0.8888888888888888) * Float64(y / Float64(t * t))), y))); end return tmp end
code[x_, y_, z_, t_, a_, b_, c_] := Block[{t$95$1 = N[(N[(N[(z * N[Sqrt[N[(t + a), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] / t), $MachinePrecision] + N[(N[(b - c), $MachinePrecision] * N[(N[(2.0 / N[(t * 3.0), $MachinePrecision]), $MachinePrecision] - N[(a + N[(5.0 / 6.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$1, -200.0], 1.0, If[LessEqual[t$95$1, 4e+288], N[(x / N[(y * N[(c * N[(c * N[(c * 0.7716049382716049 + 1.3888888888888888), $MachinePrecision] + 1.6666666666666667), $MachinePrecision] + 1.0), $MachinePrecision] + x), $MachinePrecision]), $MachinePrecision], N[(x / N[(x + N[(b * N[(N[(b * 0.8888888888888888), $MachinePrecision] * N[(y / N[(t * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{z \cdot \sqrt{t + a}}{t} + \left(b - c\right) \cdot \left(\frac{2}{t \cdot 3} - \left(a + \frac{5}{6}\right)\right)\\
\mathbf{if}\;t\_1 \leq -200:\\
\;\;\;\;1\\
\mathbf{elif}\;t\_1 \leq 4 \cdot 10^{+288}:\\
\;\;\;\;\frac{x}{\mathsf{fma}\left(y, \mathsf{fma}\left(c, \mathsf{fma}\left(c, \mathsf{fma}\left(c, 0.7716049382716049, 1.3888888888888888\right), 1.6666666666666667\right), 1\right), x\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{x}{x + \mathsf{fma}\left(b, \left(b \cdot 0.8888888888888888\right) \cdot \frac{y}{t \cdot t}, y\right)}\\
\end{array}
\end{array}
if (-.f64 (/.f64 (*.f64 z (sqrt.f64 (+.f64 t a))) t) (*.f64 (-.f64 b c) (-.f64 (+.f64 a (/.f64 #s(literal 5 binary64) #s(literal 6 binary64))) (/.f64 #s(literal 2 binary64) (*.f64 t #s(literal 3 binary64)))))) < -200Initial program 99.1%
Taylor expanded in x around inf
Simplified99.1%
*-inverses99.1
Applied egg-rr99.1%
if -200 < (-.f64 (/.f64 (*.f64 z (sqrt.f64 (+.f64 t a))) t) (*.f64 (-.f64 b c) (-.f64 (+.f64 a (/.f64 #s(literal 5 binary64) #s(literal 6 binary64))) (/.f64 #s(literal 2 binary64) (*.f64 t #s(literal 3 binary64)))))) < 4e288Initial program 99.9%
Taylor expanded in c around inf
*-lowering-*.f64N/A
+-commutativeN/A
associate--l+N/A
+-lowering-+.f64N/A
sub-negN/A
+-lowering-+.f64N/A
associate-*r/N/A
metadata-evalN/A
distribute-neg-fracN/A
metadata-evalN/A
/-lowering-/.f6474.0
Simplified74.0%
Taylor expanded in t around inf
+-commutativeN/A
accelerator-lowering-fma.f64N/A
exp-lowering-exp.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f6464.4
Simplified64.4%
Taylor expanded in a around 0
*-lowering-*.f6458.2
Simplified58.2%
Taylor expanded in c around 0
+-commutativeN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f6459.6
Simplified59.6%
if 4e288 < (-.f64 (/.f64 (*.f64 z (sqrt.f64 (+.f64 t a))) t) (*.f64 (-.f64 b c) (-.f64 (+.f64 a (/.f64 #s(literal 5 binary64) #s(literal 6 binary64))) (/.f64 #s(literal 2 binary64) (*.f64 t #s(literal 3 binary64)))))) Initial program 94.8%
Taylor expanded in b around inf
*-lowering-*.f64N/A
--lowering--.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f64N/A
+-commutativeN/A
+-lowering-+.f6468.3
Simplified68.3%
Taylor expanded in b around 0
+-lowering-+.f64N/A
+-commutativeN/A
accelerator-lowering-fma.f64N/A
Simplified79.5%
Taylor expanded in t around 0
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6454.1
Simplified54.1%
associate-*r*N/A
associate-/l*N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f6461.1
Applied egg-rr61.1%
Final simplification75.8%
(FPCore (x y z t a b c)
:precision binary64
(let* ((t_1
(+
(/ (* z (sqrt (+ t a))) t)
(* (- b c) (- (/ 2.0 (* t 3.0)) (+ a (/ 5.0 6.0)))))))
(if (<= t_1 -200.0)
1.0
(if (<= t_1 1e+301)
(/
x
(fma
y
(fma
c
(fma
c
(fma c 0.7716049382716049 1.3888888888888888)
1.6666666666666667)
1.0)
x))
(/
(* x 0.5)
(*
(* b b)
(* y (* (+ a 0.8333333333333334) (+ a 0.8333333333333334)))))))))
double code(double x, double y, double z, double t, double a, double b, double c) {
double t_1 = ((z * sqrt((t + a))) / t) + ((b - c) * ((2.0 / (t * 3.0)) - (a + (5.0 / 6.0))));
double tmp;
if (t_1 <= -200.0) {
tmp = 1.0;
} else if (t_1 <= 1e+301) {
tmp = x / fma(y, fma(c, fma(c, fma(c, 0.7716049382716049, 1.3888888888888888), 1.6666666666666667), 1.0), x);
} else {
tmp = (x * 0.5) / ((b * b) * (y * ((a + 0.8333333333333334) * (a + 0.8333333333333334))));
}
return tmp;
}
function code(x, y, z, t, a, b, c) t_1 = Float64(Float64(Float64(z * sqrt(Float64(t + a))) / t) + Float64(Float64(b - c) * Float64(Float64(2.0 / Float64(t * 3.0)) - Float64(a + Float64(5.0 / 6.0))))) tmp = 0.0 if (t_1 <= -200.0) tmp = 1.0; elseif (t_1 <= 1e+301) tmp = Float64(x / fma(y, fma(c, fma(c, fma(c, 0.7716049382716049, 1.3888888888888888), 1.6666666666666667), 1.0), x)); else tmp = Float64(Float64(x * 0.5) / Float64(Float64(b * b) * Float64(y * Float64(Float64(a + 0.8333333333333334) * Float64(a + 0.8333333333333334))))); end return tmp end
code[x_, y_, z_, t_, a_, b_, c_] := Block[{t$95$1 = N[(N[(N[(z * N[Sqrt[N[(t + a), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] / t), $MachinePrecision] + N[(N[(b - c), $MachinePrecision] * N[(N[(2.0 / N[(t * 3.0), $MachinePrecision]), $MachinePrecision] - N[(a + N[(5.0 / 6.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$1, -200.0], 1.0, If[LessEqual[t$95$1, 1e+301], N[(x / N[(y * N[(c * N[(c * N[(c * 0.7716049382716049 + 1.3888888888888888), $MachinePrecision] + 1.6666666666666667), $MachinePrecision] + 1.0), $MachinePrecision] + x), $MachinePrecision]), $MachinePrecision], N[(N[(x * 0.5), $MachinePrecision] / N[(N[(b * b), $MachinePrecision] * N[(y * N[(N[(a + 0.8333333333333334), $MachinePrecision] * N[(a + 0.8333333333333334), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{z \cdot \sqrt{t + a}}{t} + \left(b - c\right) \cdot \left(\frac{2}{t \cdot 3} - \left(a + \frac{5}{6}\right)\right)\\
\mathbf{if}\;t\_1 \leq -200:\\
\;\;\;\;1\\
\mathbf{elif}\;t\_1 \leq 10^{+301}:\\
\;\;\;\;\frac{x}{\mathsf{fma}\left(y, \mathsf{fma}\left(c, \mathsf{fma}\left(c, \mathsf{fma}\left(c, 0.7716049382716049, 1.3888888888888888\right), 1.6666666666666667\right), 1\right), x\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{x \cdot 0.5}{\left(b \cdot b\right) \cdot \left(y \cdot \left(\left(a + 0.8333333333333334\right) \cdot \left(a + 0.8333333333333334\right)\right)\right)}\\
\end{array}
\end{array}
if (-.f64 (/.f64 (*.f64 z (sqrt.f64 (+.f64 t a))) t) (*.f64 (-.f64 b c) (-.f64 (+.f64 a (/.f64 #s(literal 5 binary64) #s(literal 6 binary64))) (/.f64 #s(literal 2 binary64) (*.f64 t #s(literal 3 binary64)))))) < -200Initial program 99.1%
Taylor expanded in x around inf
Simplified99.1%
*-inverses99.1
Applied egg-rr99.1%
if -200 < (-.f64 (/.f64 (*.f64 z (sqrt.f64 (+.f64 t a))) t) (*.f64 (-.f64 b c) (-.f64 (+.f64 a (/.f64 #s(literal 5 binary64) #s(literal 6 binary64))) (/.f64 #s(literal 2 binary64) (*.f64 t #s(literal 3 binary64)))))) < 1.00000000000000005e301Initial program 99.9%
Taylor expanded in c around inf
*-lowering-*.f64N/A
+-commutativeN/A
associate--l+N/A
+-lowering-+.f64N/A
sub-negN/A
+-lowering-+.f64N/A
associate-*r/N/A
metadata-evalN/A
distribute-neg-fracN/A
metadata-evalN/A
/-lowering-/.f6472.6
Simplified72.6%
Taylor expanded in t around inf
+-commutativeN/A
accelerator-lowering-fma.f64N/A
exp-lowering-exp.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f6463.4
Simplified63.4%
Taylor expanded in a around 0
*-lowering-*.f6457.4
Simplified57.4%
Taylor expanded in c around 0
+-commutativeN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f6458.7
Simplified58.7%
if 1.00000000000000005e301 < (-.f64 (/.f64 (*.f64 z (sqrt.f64 (+.f64 t a))) t) (*.f64 (-.f64 b c) (-.f64 (+.f64 a (/.f64 #s(literal 5 binary64) #s(literal 6 binary64))) (/.f64 #s(literal 2 binary64) (*.f64 t #s(literal 3 binary64)))))) Initial program 94.6%
Taylor expanded in b around inf
*-lowering-*.f64N/A
--lowering--.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f64N/A
+-commutativeN/A
+-lowering-+.f6468.3
Simplified68.3%
Taylor expanded in b around 0
+-lowering-+.f64N/A
+-commutativeN/A
accelerator-lowering-fma.f64N/A
Simplified78.7%
Taylor expanded in t around inf
/-lowering-/.f64N/A
+-lowering-+.f64N/A
+-commutativeN/A
associate-*r*N/A
accelerator-lowering-fma.f64N/A
Simplified65.5%
Taylor expanded in b around inf
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f6461.2
Simplified61.2%
Final simplification75.5%
(FPCore (x y z t a b c)
:precision binary64
(let* ((t_1
(+
(/ (* z (sqrt (+ t a))) t)
(* (- b c) (- (/ 2.0 (* t 3.0)) (+ a (/ 5.0 6.0)))))))
(if (<= t_1 -200.0)
1.0
(if (<= t_1 2e+307)
(/
x
(fma
y
(fma
c
(fma
c
(fma c 0.7716049382716049 1.3888888888888888)
1.6666666666666667)
1.0)
x))
(/ x (+ x (fma (* 2.0 b) (* y (* b (* a a))) y)))))))
double code(double x, double y, double z, double t, double a, double b, double c) {
double t_1 = ((z * sqrt((t + a))) / t) + ((b - c) * ((2.0 / (t * 3.0)) - (a + (5.0 / 6.0))));
double tmp;
if (t_1 <= -200.0) {
tmp = 1.0;
} else if (t_1 <= 2e+307) {
tmp = x / fma(y, fma(c, fma(c, fma(c, 0.7716049382716049, 1.3888888888888888), 1.6666666666666667), 1.0), x);
} else {
tmp = x / (x + fma((2.0 * b), (y * (b * (a * a))), y));
}
return tmp;
}
function code(x, y, z, t, a, b, c) t_1 = Float64(Float64(Float64(z * sqrt(Float64(t + a))) / t) + Float64(Float64(b - c) * Float64(Float64(2.0 / Float64(t * 3.0)) - Float64(a + Float64(5.0 / 6.0))))) tmp = 0.0 if (t_1 <= -200.0) tmp = 1.0; elseif (t_1 <= 2e+307) tmp = Float64(x / fma(y, fma(c, fma(c, fma(c, 0.7716049382716049, 1.3888888888888888), 1.6666666666666667), 1.0), x)); else tmp = Float64(x / Float64(x + fma(Float64(2.0 * b), Float64(y * Float64(b * Float64(a * a))), y))); end return tmp end
code[x_, y_, z_, t_, a_, b_, c_] := Block[{t$95$1 = N[(N[(N[(z * N[Sqrt[N[(t + a), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] / t), $MachinePrecision] + N[(N[(b - c), $MachinePrecision] * N[(N[(2.0 / N[(t * 3.0), $MachinePrecision]), $MachinePrecision] - N[(a + N[(5.0 / 6.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$1, -200.0], 1.0, If[LessEqual[t$95$1, 2e+307], N[(x / N[(y * N[(c * N[(c * N[(c * 0.7716049382716049 + 1.3888888888888888), $MachinePrecision] + 1.6666666666666667), $MachinePrecision] + 1.0), $MachinePrecision] + x), $MachinePrecision]), $MachinePrecision], N[(x / N[(x + N[(N[(2.0 * b), $MachinePrecision] * N[(y * N[(b * N[(a * a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{z \cdot \sqrt{t + a}}{t} + \left(b - c\right) \cdot \left(\frac{2}{t \cdot 3} - \left(a + \frac{5}{6}\right)\right)\\
\mathbf{if}\;t\_1 \leq -200:\\
\;\;\;\;1\\
\mathbf{elif}\;t\_1 \leq 2 \cdot 10^{+307}:\\
\;\;\;\;\frac{x}{\mathsf{fma}\left(y, \mathsf{fma}\left(c, \mathsf{fma}\left(c, \mathsf{fma}\left(c, 0.7716049382716049, 1.3888888888888888\right), 1.6666666666666667\right), 1\right), x\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{x}{x + \mathsf{fma}\left(2 \cdot b, y \cdot \left(b \cdot \left(a \cdot a\right)\right), y\right)}\\
\end{array}
\end{array}
if (-.f64 (/.f64 (*.f64 z (sqrt.f64 (+.f64 t a))) t) (*.f64 (-.f64 b c) (-.f64 (+.f64 a (/.f64 #s(literal 5 binary64) #s(literal 6 binary64))) (/.f64 #s(literal 2 binary64) (*.f64 t #s(literal 3 binary64)))))) < -200Initial program 99.1%
Taylor expanded in x around inf
Simplified99.1%
*-inverses99.1
Applied egg-rr99.1%
if -200 < (-.f64 (/.f64 (*.f64 z (sqrt.f64 (+.f64 t a))) t) (*.f64 (-.f64 b c) (-.f64 (+.f64 a (/.f64 #s(literal 5 binary64) #s(literal 6 binary64))) (/.f64 #s(literal 2 binary64) (*.f64 t #s(literal 3 binary64)))))) < 1.99999999999999997e307Initial program 99.9%
Taylor expanded in c around inf
*-lowering-*.f64N/A
+-commutativeN/A
associate--l+N/A
+-lowering-+.f64N/A
sub-negN/A
+-lowering-+.f64N/A
associate-*r/N/A
metadata-evalN/A
distribute-neg-fracN/A
metadata-evalN/A
/-lowering-/.f6471.9
Simplified71.9%
Taylor expanded in t around inf
+-commutativeN/A
accelerator-lowering-fma.f64N/A
exp-lowering-exp.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f6462.1
Simplified62.1%
Taylor expanded in a around 0
*-lowering-*.f6456.5
Simplified56.5%
Taylor expanded in c around 0
+-commutativeN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f6457.8
Simplified57.8%
if 1.99999999999999997e307 < (-.f64 (/.f64 (*.f64 z (sqrt.f64 (+.f64 t a))) t) (*.f64 (-.f64 b c) (-.f64 (+.f64 a (/.f64 #s(literal 5 binary64) #s(literal 6 binary64))) (/.f64 #s(literal 2 binary64) (*.f64 t #s(literal 3 binary64)))))) Initial program 94.2%
Taylor expanded in b around inf
*-lowering-*.f64N/A
--lowering--.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f64N/A
+-commutativeN/A
+-lowering-+.f6468.8
Simplified68.8%
Taylor expanded in b around 0
+-lowering-+.f64N/A
+-commutativeN/A
accelerator-lowering-fma.f64N/A
Simplified78.5%
Taylor expanded in t around inf
/-lowering-/.f64N/A
+-lowering-+.f64N/A
+-commutativeN/A
associate-*r*N/A
accelerator-lowering-fma.f64N/A
Simplified64.4%
Taylor expanded in a around inf
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6457.7
Simplified57.7%
Final simplification74.2%
(FPCore (x y z t a b c)
:precision binary64
(let* ((t_1
(+
(/ (* z (sqrt (+ t a))) t)
(* (- b c) (- (/ 2.0 (* t 3.0)) (+ a (/ 5.0 6.0)))))))
(if (<= t_1 -200.0)
1.0
(if (<= t_1 5e+218)
(/ x (fma y (fma 1.6666666666666667 c 1.0) x))
(/ x (+ x (* 2.0 (* a (* y c)))))))))
double code(double x, double y, double z, double t, double a, double b, double c) {
double t_1 = ((z * sqrt((t + a))) / t) + ((b - c) * ((2.0 / (t * 3.0)) - (a + (5.0 / 6.0))));
double tmp;
if (t_1 <= -200.0) {
tmp = 1.0;
} else if (t_1 <= 5e+218) {
tmp = x / fma(y, fma(1.6666666666666667, c, 1.0), x);
} else {
tmp = x / (x + (2.0 * (a * (y * c))));
}
return tmp;
}
function code(x, y, z, t, a, b, c) t_1 = Float64(Float64(Float64(z * sqrt(Float64(t + a))) / t) + Float64(Float64(b - c) * Float64(Float64(2.0 / Float64(t * 3.0)) - Float64(a + Float64(5.0 / 6.0))))) tmp = 0.0 if (t_1 <= -200.0) tmp = 1.0; elseif (t_1 <= 5e+218) tmp = Float64(x / fma(y, fma(1.6666666666666667, c, 1.0), x)); else tmp = Float64(x / Float64(x + Float64(2.0 * Float64(a * Float64(y * c))))); end return tmp end
code[x_, y_, z_, t_, a_, b_, c_] := Block[{t$95$1 = N[(N[(N[(z * N[Sqrt[N[(t + a), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] / t), $MachinePrecision] + N[(N[(b - c), $MachinePrecision] * N[(N[(2.0 / N[(t * 3.0), $MachinePrecision]), $MachinePrecision] - N[(a + N[(5.0 / 6.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$1, -200.0], 1.0, If[LessEqual[t$95$1, 5e+218], N[(x / N[(y * N[(1.6666666666666667 * c + 1.0), $MachinePrecision] + x), $MachinePrecision]), $MachinePrecision], N[(x / N[(x + N[(2.0 * N[(a * N[(y * c), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{z \cdot \sqrt{t + a}}{t} + \left(b - c\right) \cdot \left(\frac{2}{t \cdot 3} - \left(a + \frac{5}{6}\right)\right)\\
\mathbf{if}\;t\_1 \leq -200:\\
\;\;\;\;1\\
\mathbf{elif}\;t\_1 \leq 5 \cdot 10^{+218}:\\
\;\;\;\;\frac{x}{\mathsf{fma}\left(y, \mathsf{fma}\left(1.6666666666666667, c, 1\right), x\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{x}{x + 2 \cdot \left(a \cdot \left(y \cdot c\right)\right)}\\
\end{array}
\end{array}
if (-.f64 (/.f64 (*.f64 z (sqrt.f64 (+.f64 t a))) t) (*.f64 (-.f64 b c) (-.f64 (+.f64 a (/.f64 #s(literal 5 binary64) #s(literal 6 binary64))) (/.f64 #s(literal 2 binary64) (*.f64 t #s(literal 3 binary64)))))) < -200Initial program 99.1%
Taylor expanded in x around inf
Simplified99.1%
*-inverses99.1
Applied egg-rr99.1%
if -200 < (-.f64 (/.f64 (*.f64 z (sqrt.f64 (+.f64 t a))) t) (*.f64 (-.f64 b c) (-.f64 (+.f64 a (/.f64 #s(literal 5 binary64) #s(literal 6 binary64))) (/.f64 #s(literal 2 binary64) (*.f64 t #s(literal 3 binary64)))))) < 4.99999999999999983e218Initial program 99.9%
Taylor expanded in c around inf
*-lowering-*.f64N/A
+-commutativeN/A
associate--l+N/A
+-lowering-+.f64N/A
sub-negN/A
+-lowering-+.f64N/A
associate-*r/N/A
metadata-evalN/A
distribute-neg-fracN/A
metadata-evalN/A
/-lowering-/.f6470.1
Simplified70.1%
Taylor expanded in t around inf
+-commutativeN/A
accelerator-lowering-fma.f64N/A
exp-lowering-exp.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f6463.9
Simplified63.9%
Taylor expanded in a around 0
*-lowering-*.f6462.2
Simplified62.2%
Taylor expanded in c around 0
+-commutativeN/A
accelerator-lowering-fma.f6447.3
Simplified47.3%
if 4.99999999999999983e218 < (-.f64 (/.f64 (*.f64 z (sqrt.f64 (+.f64 t a))) t) (*.f64 (-.f64 b c) (-.f64 (+.f64 a (/.f64 #s(literal 5 binary64) #s(literal 6 binary64))) (/.f64 #s(literal 2 binary64) (*.f64 t #s(literal 3 binary64)))))) Initial program 95.9%
Taylor expanded in c around inf
*-lowering-*.f64N/A
+-commutativeN/A
associate--l+N/A
+-lowering-+.f64N/A
sub-negN/A
+-lowering-+.f64N/A
associate-*r/N/A
metadata-evalN/A
distribute-neg-fracN/A
metadata-evalN/A
/-lowering-/.f6462.8
Simplified62.8%
Taylor expanded in c around 0
+-lowering-+.f64N/A
+-commutativeN/A
accelerator-lowering-fma.f64N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
associate--l+N/A
+-lowering-+.f64N/A
--lowering--.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f6443.8
Simplified43.8%
Taylor expanded in a around inf
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6438.2
Simplified38.2%
Final simplification64.5%
(FPCore (x y z t a b c)
:precision binary64
(let* ((t_1
(+
(/ (* z (sqrt (+ t a))) t)
(* (- b c) (- (/ 2.0 (* t 3.0)) (+ a (/ 5.0 6.0)))))))
(if (<= t_1 -200.0)
1.0
(if (<= t_1 1e+241)
(/ x (fma y (fma 1.6666666666666667 c 1.0) x))
(/ (* x 0.5) (* a (* y c)))))))
double code(double x, double y, double z, double t, double a, double b, double c) {
double t_1 = ((z * sqrt((t + a))) / t) + ((b - c) * ((2.0 / (t * 3.0)) - (a + (5.0 / 6.0))));
double tmp;
if (t_1 <= -200.0) {
tmp = 1.0;
} else if (t_1 <= 1e+241) {
tmp = x / fma(y, fma(1.6666666666666667, c, 1.0), x);
} else {
tmp = (x * 0.5) / (a * (y * c));
}
return tmp;
}
function code(x, y, z, t, a, b, c) t_1 = Float64(Float64(Float64(z * sqrt(Float64(t + a))) / t) + Float64(Float64(b - c) * Float64(Float64(2.0 / Float64(t * 3.0)) - Float64(a + Float64(5.0 / 6.0))))) tmp = 0.0 if (t_1 <= -200.0) tmp = 1.0; elseif (t_1 <= 1e+241) tmp = Float64(x / fma(y, fma(1.6666666666666667, c, 1.0), x)); else tmp = Float64(Float64(x * 0.5) / Float64(a * Float64(y * c))); end return tmp end
code[x_, y_, z_, t_, a_, b_, c_] := Block[{t$95$1 = N[(N[(N[(z * N[Sqrt[N[(t + a), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] / t), $MachinePrecision] + N[(N[(b - c), $MachinePrecision] * N[(N[(2.0 / N[(t * 3.0), $MachinePrecision]), $MachinePrecision] - N[(a + N[(5.0 / 6.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$1, -200.0], 1.0, If[LessEqual[t$95$1, 1e+241], N[(x / N[(y * N[(1.6666666666666667 * c + 1.0), $MachinePrecision] + x), $MachinePrecision]), $MachinePrecision], N[(N[(x * 0.5), $MachinePrecision] / N[(a * N[(y * c), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{z \cdot \sqrt{t + a}}{t} + \left(b - c\right) \cdot \left(\frac{2}{t \cdot 3} - \left(a + \frac{5}{6}\right)\right)\\
\mathbf{if}\;t\_1 \leq -200:\\
\;\;\;\;1\\
\mathbf{elif}\;t\_1 \leq 10^{+241}:\\
\;\;\;\;\frac{x}{\mathsf{fma}\left(y, \mathsf{fma}\left(1.6666666666666667, c, 1\right), x\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{x \cdot 0.5}{a \cdot \left(y \cdot c\right)}\\
\end{array}
\end{array}
if (-.f64 (/.f64 (*.f64 z (sqrt.f64 (+.f64 t a))) t) (*.f64 (-.f64 b c) (-.f64 (+.f64 a (/.f64 #s(literal 5 binary64) #s(literal 6 binary64))) (/.f64 #s(literal 2 binary64) (*.f64 t #s(literal 3 binary64)))))) < -200Initial program 99.1%
Taylor expanded in x around inf
Simplified99.1%
*-inverses99.1
Applied egg-rr99.1%
if -200 < (-.f64 (/.f64 (*.f64 z (sqrt.f64 (+.f64 t a))) t) (*.f64 (-.f64 b c) (-.f64 (+.f64 a (/.f64 #s(literal 5 binary64) #s(literal 6 binary64))) (/.f64 #s(literal 2 binary64) (*.f64 t #s(literal 3 binary64)))))) < 1.0000000000000001e241Initial program 99.9%
Taylor expanded in c around inf
*-lowering-*.f64N/A
+-commutativeN/A
associate--l+N/A
+-lowering-+.f64N/A
sub-negN/A
+-lowering-+.f64N/A
associate-*r/N/A
metadata-evalN/A
distribute-neg-fracN/A
metadata-evalN/A
/-lowering-/.f6470.0
Simplified70.0%
Taylor expanded in t around inf
+-commutativeN/A
accelerator-lowering-fma.f64N/A
exp-lowering-exp.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f6464.0
Simplified64.0%
Taylor expanded in a around 0
*-lowering-*.f6462.5
Simplified62.5%
Taylor expanded in c around 0
+-commutativeN/A
accelerator-lowering-fma.f6445.1
Simplified45.1%
if 1.0000000000000001e241 < (-.f64 (/.f64 (*.f64 z (sqrt.f64 (+.f64 t a))) t) (*.f64 (-.f64 b c) (-.f64 (+.f64 a (/.f64 #s(literal 5 binary64) #s(literal 6 binary64))) (/.f64 #s(literal 2 binary64) (*.f64 t #s(literal 3 binary64)))))) Initial program 95.7%
Taylor expanded in c around inf
*-lowering-*.f64N/A
+-commutativeN/A
associate--l+N/A
+-lowering-+.f64N/A
sub-negN/A
+-lowering-+.f64N/A
associate-*r/N/A
metadata-evalN/A
distribute-neg-fracN/A
metadata-evalN/A
/-lowering-/.f6462.6
Simplified62.6%
Taylor expanded in c around 0
+-lowering-+.f64N/A
+-commutativeN/A
accelerator-lowering-fma.f64N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
associate--l+N/A
+-lowering-+.f64N/A
--lowering--.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f6445.1
Simplified45.1%
Taylor expanded in a around inf
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6436.7
Simplified36.7%
Final simplification63.5%
(FPCore (x y z t a b c)
:precision binary64
(if (<=
(+
(/ (* z (sqrt (+ t a))) t)
(* (- b c) (- (/ 2.0 (* t 3.0)) (+ a (/ 5.0 6.0)))))
-200.0)
1.0
(/
x
(fma
y
(fma
c
(fma c (fma c 0.7716049382716049 1.3888888888888888) 1.6666666666666667)
1.0)
x))))
double code(double x, double y, double z, double t, double a, double b, double c) {
double tmp;
if ((((z * sqrt((t + a))) / t) + ((b - c) * ((2.0 / (t * 3.0)) - (a + (5.0 / 6.0))))) <= -200.0) {
tmp = 1.0;
} else {
tmp = x / fma(y, fma(c, fma(c, fma(c, 0.7716049382716049, 1.3888888888888888), 1.6666666666666667), 1.0), x);
}
return tmp;
}
function code(x, y, z, t, a, b, c) tmp = 0.0 if (Float64(Float64(Float64(z * sqrt(Float64(t + a))) / t) + Float64(Float64(b - c) * Float64(Float64(2.0 / Float64(t * 3.0)) - Float64(a + Float64(5.0 / 6.0))))) <= -200.0) tmp = 1.0; else tmp = Float64(x / fma(y, fma(c, fma(c, fma(c, 0.7716049382716049, 1.3888888888888888), 1.6666666666666667), 1.0), x)); end return tmp end
code[x_, y_, z_, t_, a_, b_, c_] := If[LessEqual[N[(N[(N[(z * N[Sqrt[N[(t + a), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] / t), $MachinePrecision] + N[(N[(b - c), $MachinePrecision] * N[(N[(2.0 / N[(t * 3.0), $MachinePrecision]), $MachinePrecision] - N[(a + N[(5.0 / 6.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], -200.0], 1.0, N[(x / N[(y * N[(c * N[(c * N[(c * 0.7716049382716049 + 1.3888888888888888), $MachinePrecision] + 1.6666666666666667), $MachinePrecision] + 1.0), $MachinePrecision] + x), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\frac{z \cdot \sqrt{t + a}}{t} + \left(b - c\right) \cdot \left(\frac{2}{t \cdot 3} - \left(a + \frac{5}{6}\right)\right) \leq -200:\\
\;\;\;\;1\\
\mathbf{else}:\\
\;\;\;\;\frac{x}{\mathsf{fma}\left(y, \mathsf{fma}\left(c, \mathsf{fma}\left(c, \mathsf{fma}\left(c, 0.7716049382716049, 1.3888888888888888\right), 1.6666666666666667\right), 1\right), x\right)}\\
\end{array}
\end{array}
if (-.f64 (/.f64 (*.f64 z (sqrt.f64 (+.f64 t a))) t) (*.f64 (-.f64 b c) (-.f64 (+.f64 a (/.f64 #s(literal 5 binary64) #s(literal 6 binary64))) (/.f64 #s(literal 2 binary64) (*.f64 t #s(literal 3 binary64)))))) < -200Initial program 99.1%
Taylor expanded in x around inf
Simplified99.1%
*-inverses99.1
Applied egg-rr99.1%
if -200 < (-.f64 (/.f64 (*.f64 z (sqrt.f64 (+.f64 t a))) t) (*.f64 (-.f64 b c) (-.f64 (+.f64 a (/.f64 #s(literal 5 binary64) #s(literal 6 binary64))) (/.f64 #s(literal 2 binary64) (*.f64 t #s(literal 3 binary64)))))) Initial program 97.4%
Taylor expanded in c around inf
*-lowering-*.f64N/A
+-commutativeN/A
associate--l+N/A
+-lowering-+.f64N/A
sub-negN/A
+-lowering-+.f64N/A
associate-*r/N/A
metadata-evalN/A
distribute-neg-fracN/A
metadata-evalN/A
/-lowering-/.f6465.5
Simplified65.5%
Taylor expanded in t around inf
+-commutativeN/A
accelerator-lowering-fma.f64N/A
exp-lowering-exp.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f6456.9
Simplified56.9%
Taylor expanded in a around 0
*-lowering-*.f6451.3
Simplified51.3%
Taylor expanded in c around 0
+-commutativeN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
accelerator-lowering-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
accelerator-lowering-fma.f6453.3
Simplified53.3%
Final simplification71.6%
(FPCore (x y z t a b c) :precision binary64 1.0)
double code(double x, double y, double z, double t, double a, double b, double c) {
return 1.0;
}
real(8) function code(x, y, z, t, a, b, c)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8), intent (in) :: c
code = 1.0d0
end function
public static double code(double x, double y, double z, double t, double a, double b, double c) {
return 1.0;
}
def code(x, y, z, t, a, b, c): return 1.0
function code(x, y, z, t, a, b, c) return 1.0 end
function tmp = code(x, y, z, t, a, b, c) tmp = 1.0; end
code[x_, y_, z_, t_, a_, b_, c_] := 1.0
\begin{array}{l}
\\
1
\end{array}
Initial program 98.1%
Taylor expanded in x around inf
Simplified45.1%
*-inverses45.1
Applied egg-rr45.1%
(FPCore (x y z t a b c)
:precision binary64
(let* ((t_1 (* z (sqrt (+ t a)))) (t_2 (- a (/ 5.0 6.0))))
(if (< t -2.118326644891581e-50)
(/
x
(+
x
(* y (exp (* 2.0 (- (+ (* a c) (* 0.8333333333333334 c)) (* a b)))))))
(if (< t 5.196588770651547e-123)
(/
x
(+
x
(*
y
(exp
(*
2.0
(/
(-
(* t_1 (* (* 3.0 t) t_2))
(*
(- (* (+ (/ 5.0 6.0) a) (* 3.0 t)) 2.0)
(* t_2 (* (- b c) t))))
(* (* (* t t) 3.0) t_2)))))))
(/
x
(+
x
(*
y
(exp
(*
2.0
(-
(/ t_1 t)
(* (- b c) (- (+ a (/ 5.0 6.0)) (/ 2.0 (* t 3.0))))))))))))))
double code(double x, double y, double z, double t, double a, double b, double c) {
double t_1 = z * sqrt((t + a));
double t_2 = a - (5.0 / 6.0);
double tmp;
if (t < -2.118326644891581e-50) {
tmp = x / (x + (y * exp((2.0 * (((a * c) + (0.8333333333333334 * c)) - (a * b))))));
} else if (t < 5.196588770651547e-123) {
tmp = x / (x + (y * exp((2.0 * (((t_1 * ((3.0 * t) * t_2)) - (((((5.0 / 6.0) + a) * (3.0 * t)) - 2.0) * (t_2 * ((b - c) * t)))) / (((t * t) * 3.0) * t_2))))));
} else {
tmp = x / (x + (y * exp((2.0 * ((t_1 / t) - ((b - c) * ((a + (5.0 / 6.0)) - (2.0 / (t * 3.0)))))))));
}
return tmp;
}
real(8) function code(x, y, z, t, a, b, c)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8), intent (in) :: c
real(8) :: t_1
real(8) :: t_2
real(8) :: tmp
t_1 = z * sqrt((t + a))
t_2 = a - (5.0d0 / 6.0d0)
if (t < (-2.118326644891581d-50)) then
tmp = x / (x + (y * exp((2.0d0 * (((a * c) + (0.8333333333333334d0 * c)) - (a * b))))))
else if (t < 5.196588770651547d-123) then
tmp = x / (x + (y * exp((2.0d0 * (((t_1 * ((3.0d0 * t) * t_2)) - (((((5.0d0 / 6.0d0) + a) * (3.0d0 * t)) - 2.0d0) * (t_2 * ((b - c) * t)))) / (((t * t) * 3.0d0) * t_2))))))
else
tmp = x / (x + (y * exp((2.0d0 * ((t_1 / t) - ((b - c) * ((a + (5.0d0 / 6.0d0)) - (2.0d0 / (t * 3.0d0)))))))))
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b, double c) {
double t_1 = z * Math.sqrt((t + a));
double t_2 = a - (5.0 / 6.0);
double tmp;
if (t < -2.118326644891581e-50) {
tmp = x / (x + (y * Math.exp((2.0 * (((a * c) + (0.8333333333333334 * c)) - (a * b))))));
} else if (t < 5.196588770651547e-123) {
tmp = x / (x + (y * Math.exp((2.0 * (((t_1 * ((3.0 * t) * t_2)) - (((((5.0 / 6.0) + a) * (3.0 * t)) - 2.0) * (t_2 * ((b - c) * t)))) / (((t * t) * 3.0) * t_2))))));
} else {
tmp = x / (x + (y * Math.exp((2.0 * ((t_1 / t) - ((b - c) * ((a + (5.0 / 6.0)) - (2.0 / (t * 3.0)))))))));
}
return tmp;
}
def code(x, y, z, t, a, b, c): t_1 = z * math.sqrt((t + a)) t_2 = a - (5.0 / 6.0) tmp = 0 if t < -2.118326644891581e-50: tmp = x / (x + (y * math.exp((2.0 * (((a * c) + (0.8333333333333334 * c)) - (a * b)))))) elif t < 5.196588770651547e-123: tmp = x / (x + (y * math.exp((2.0 * (((t_1 * ((3.0 * t) * t_2)) - (((((5.0 / 6.0) + a) * (3.0 * t)) - 2.0) * (t_2 * ((b - c) * t)))) / (((t * t) * 3.0) * t_2)))))) else: tmp = x / (x + (y * math.exp((2.0 * ((t_1 / t) - ((b - c) * ((a + (5.0 / 6.0)) - (2.0 / (t * 3.0))))))))) return tmp
function code(x, y, z, t, a, b, c) t_1 = Float64(z * sqrt(Float64(t + a))) t_2 = Float64(a - Float64(5.0 / 6.0)) tmp = 0.0 if (t < -2.118326644891581e-50) tmp = Float64(x / Float64(x + Float64(y * exp(Float64(2.0 * Float64(Float64(Float64(a * c) + Float64(0.8333333333333334 * c)) - Float64(a * b))))))); elseif (t < 5.196588770651547e-123) tmp = Float64(x / Float64(x + Float64(y * exp(Float64(2.0 * Float64(Float64(Float64(t_1 * Float64(Float64(3.0 * t) * t_2)) - Float64(Float64(Float64(Float64(Float64(5.0 / 6.0) + a) * Float64(3.0 * t)) - 2.0) * Float64(t_2 * Float64(Float64(b - c) * t)))) / Float64(Float64(Float64(t * t) * 3.0) * t_2))))))); else tmp = Float64(x / Float64(x + Float64(y * exp(Float64(2.0 * Float64(Float64(t_1 / t) - Float64(Float64(b - c) * Float64(Float64(a + Float64(5.0 / 6.0)) - Float64(2.0 / Float64(t * 3.0)))))))))); end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c) t_1 = z * sqrt((t + a)); t_2 = a - (5.0 / 6.0); tmp = 0.0; if (t < -2.118326644891581e-50) tmp = x / (x + (y * exp((2.0 * (((a * c) + (0.8333333333333334 * c)) - (a * b)))))); elseif (t < 5.196588770651547e-123) tmp = x / (x + (y * exp((2.0 * (((t_1 * ((3.0 * t) * t_2)) - (((((5.0 / 6.0) + a) * (3.0 * t)) - 2.0) * (t_2 * ((b - c) * t)))) / (((t * t) * 3.0) * t_2)))))); else tmp = x / (x + (y * exp((2.0 * ((t_1 / t) - ((b - c) * ((a + (5.0 / 6.0)) - (2.0 / (t * 3.0))))))))); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_] := Block[{t$95$1 = N[(z * N[Sqrt[N[(t + a), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(a - N[(5.0 / 6.0), $MachinePrecision]), $MachinePrecision]}, If[Less[t, -2.118326644891581e-50], N[(x / N[(x + N[(y * N[Exp[N[(2.0 * N[(N[(N[(a * c), $MachinePrecision] + N[(0.8333333333333334 * c), $MachinePrecision]), $MachinePrecision] - N[(a * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[Less[t, 5.196588770651547e-123], N[(x / N[(x + N[(y * N[Exp[N[(2.0 * N[(N[(N[(t$95$1 * N[(N[(3.0 * t), $MachinePrecision] * t$95$2), $MachinePrecision]), $MachinePrecision] - N[(N[(N[(N[(N[(5.0 / 6.0), $MachinePrecision] + a), $MachinePrecision] * N[(3.0 * t), $MachinePrecision]), $MachinePrecision] - 2.0), $MachinePrecision] * N[(t$95$2 * N[(N[(b - c), $MachinePrecision] * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(N[(N[(t * t), $MachinePrecision] * 3.0), $MachinePrecision] * t$95$2), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x / N[(x + N[(y * N[Exp[N[(2.0 * N[(N[(t$95$1 / t), $MachinePrecision] - N[(N[(b - c), $MachinePrecision] * N[(N[(a + N[(5.0 / 6.0), $MachinePrecision]), $MachinePrecision] - N[(2.0 / N[(t * 3.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := z \cdot \sqrt{t + a}\\
t_2 := a - \frac{5}{6}\\
\mathbf{if}\;t < -2.118326644891581 \cdot 10^{-50}:\\
\;\;\;\;\frac{x}{x + y \cdot e^{2 \cdot \left(\left(a \cdot c + 0.8333333333333334 \cdot c\right) - a \cdot b\right)}}\\
\mathbf{elif}\;t < 5.196588770651547 \cdot 10^{-123}:\\
\;\;\;\;\frac{x}{x + y \cdot e^{2 \cdot \frac{t\_1 \cdot \left(\left(3 \cdot t\right) \cdot t\_2\right) - \left(\left(\frac{5}{6} + a\right) \cdot \left(3 \cdot t\right) - 2\right) \cdot \left(t\_2 \cdot \left(\left(b - c\right) \cdot t\right)\right)}{\left(\left(t \cdot t\right) \cdot 3\right) \cdot t\_2}}}\\
\mathbf{else}:\\
\;\;\;\;\frac{x}{x + y \cdot e^{2 \cdot \left(\frac{t\_1}{t} - \left(b - c\right) \cdot \left(\left(a + \frac{5}{6}\right) - \frac{2}{t \cdot 3}\right)\right)}}\\
\end{array}
\end{array}
herbie shell --seed 2024195
(FPCore (x y z t a b c)
:name "Numeric.SpecFunctions:invIncompleteBetaWorker from math-functions-0.1.5.2, I"
:precision binary64
:alt
(! :herbie-platform default (if (< t -2118326644891581/100000000000000000000000000000000000000000000000000000000000000000) (/ x (+ x (* y (exp (* 2 (- (+ (* a c) (* 4166666666666667/5000000000000000 c)) (* a b))))))) (if (< t 5196588770651547/1000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000) (/ x (+ x (* y (exp (* 2 (/ (- (* (* z (sqrt (+ t a))) (* (* 3 t) (- a (/ 5 6)))) (* (- (* (+ (/ 5 6) a) (* 3 t)) 2) (* (- a (/ 5 6)) (* (- b c) t)))) (* (* (* t t) 3) (- a (/ 5 6))))))))) (/ x (+ x (* y (exp (* 2 (- (/ (* z (sqrt (+ t a))) t) (* (- b c) (- (+ a (/ 5 6)) (/ 2 (* t 3)))))))))))))
(/ x (+ x (* y (exp (* 2.0 (- (/ (* z (sqrt (+ t a))) t) (* (- b c) (- (+ a (/ 5.0 6.0)) (/ 2.0 (* t 3.0)))))))))))