
(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 14 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
(let* ((t_1 (sqrt (+ t a))))
(if (<=
(+ (/ (* z t_1) t) (* (- b c) (- (/ 2.0 (* t 3.0)) (+ a (/ 5.0 6.0)))))
INFINITY)
(/
x
(+
x
(*
y
(exp
(*
2.0
(+
(* t_1 (/ z t))
(*
(- b c)
(+ (/ 0.6666666666666666 t) (- -0.8333333333333334 a)))))))))
(/
x
(+
x
(*
y
(exp
(*
2.0
(* b (- (/ 0.6666666666666666 t) (+ a 0.8333333333333334)))))))))))
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 ((((z * t_1) / t) + ((b - c) * ((2.0 / (t * 3.0)) - (a + (5.0 / 6.0))))) <= ((double) INFINITY)) {
tmp = x / (x + (y * exp((2.0 * ((t_1 * (z / t)) + ((b - c) * ((0.6666666666666666 / t) + (-0.8333333333333334 - a))))))));
} else {
tmp = x / (x + (y * exp((2.0 * (b * ((0.6666666666666666 / t) - (a + 0.8333333333333334)))))));
}
return tmp;
}
public static double code(double x, double y, double z, double t, double a, double b, double c) {
double t_1 = Math.sqrt((t + a));
double tmp;
if ((((z * t_1) / t) + ((b - c) * ((2.0 / (t * 3.0)) - (a + (5.0 / 6.0))))) <= Double.POSITIVE_INFINITY) {
tmp = x / (x + (y * Math.exp((2.0 * ((t_1 * (z / t)) + ((b - c) * ((0.6666666666666666 / t) + (-0.8333333333333334 - a))))))));
} else {
tmp = x / (x + (y * Math.exp((2.0 * (b * ((0.6666666666666666 / t) - (a + 0.8333333333333334)))))));
}
return tmp;
}
def code(x, y, z, t, a, b, c): t_1 = math.sqrt((t + a)) tmp = 0 if (((z * t_1) / t) + ((b - c) * ((2.0 / (t * 3.0)) - (a + (5.0 / 6.0))))) <= math.inf: tmp = x / (x + (y * math.exp((2.0 * ((t_1 * (z / t)) + ((b - c) * ((0.6666666666666666 / t) + (-0.8333333333333334 - a)))))))) else: tmp = x / (x + (y * math.exp((2.0 * (b * ((0.6666666666666666 / t) - (a + 0.8333333333333334))))))) return tmp
function code(x, y, z, t, a, b, c) t_1 = sqrt(Float64(t + a)) tmp = 0.0 if (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))))) <= Inf) tmp = Float64(x / Float64(x + Float64(y * exp(Float64(2.0 * Float64(Float64(t_1 * Float64(z / t)) + Float64(Float64(b - c) * Float64(Float64(0.6666666666666666 / t) + Float64(-0.8333333333333334 - a))))))))); else tmp = Float64(x / Float64(x + Float64(y * exp(Float64(2.0 * Float64(b * Float64(Float64(0.6666666666666666 / t) - Float64(a + 0.8333333333333334)))))))); end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c) t_1 = sqrt((t + a)); tmp = 0.0; if ((((z * t_1) / t) + ((b - c) * ((2.0 / (t * 3.0)) - (a + (5.0 / 6.0))))) <= Inf) tmp = x / (x + (y * exp((2.0 * ((t_1 * (z / t)) + ((b - c) * ((0.6666666666666666 / t) + (-0.8333333333333334 - a)))))))); else tmp = x / (x + (y * exp((2.0 * (b * ((0.6666666666666666 / t) - (a + 0.8333333333333334))))))); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_] := Block[{t$95$1 = N[Sqrt[N[(t + a), $MachinePrecision]], $MachinePrecision]}, If[LessEqual[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], Infinity], N[(x / N[(x + N[(y * N[Exp[N[(2.0 * N[(N[(t$95$1 * N[(z / t), $MachinePrecision]), $MachinePrecision] + N[(N[(b - c), $MachinePrecision] * N[(N[(0.6666666666666666 / t), $MachinePrecision] + N[(-0.8333333333333334 - a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x / N[(x + N[(y * N[Exp[N[(2.0 * N[(b * N[(N[(0.6666666666666666 / t), $MachinePrecision] - N[(a + 0.8333333333333334), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \sqrt{t + a}\\
\mathbf{if}\;\frac{z \cdot t\_1}{t} + \left(b - c\right) \cdot \left(\frac{2}{t \cdot 3} - \left(a + \frac{5}{6}\right)\right) \leq \infty:\\
\;\;\;\;\frac{x}{x + y \cdot e^{2 \cdot \left(t\_1 \cdot \frac{z}{t} + \left(b - c\right) \cdot \left(\frac{0.6666666666666666}{t} + \left(-0.8333333333333334 - a\right)\right)\right)}}\\
\mathbf{else}:\\
\;\;\;\;\frac{x}{x + y \cdot e^{2 \cdot \left(b \cdot \left(\frac{0.6666666666666666}{t} - \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)))))) < +inf.0Initial program 97.9%
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
count-2N/A
exp-lowering-exp.f64N/A
Simplified100.0%
if +inf.0 < (-.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 0.0%
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
count-2N/A
exp-lowering-exp.f64N/A
Simplified9.5%
Taylor expanded in b around inf
*-lowering-*.f64N/A
--lowering--.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f64N/A
+-lowering-+.f6463.1%
Simplified63.1%
Final simplification97.0%
(FPCore (x y z t a b c)
:precision binary64
(let* ((t_1
(/
x
(+
x
(*
y
(exp
(*
2.0
(*
b
(- (/ 0.6666666666666666 t) (+ a 0.8333333333333334))))))))))
(if (<= b -1.35e+39)
t_1
(if (<= b -2.45e-205)
(/ x (+ x (* y (exp (* 2.0 (/ (* c -0.6666666666666666) t))))))
(if (<= b 1.45e-98)
(/ x (+ x (* y (exp (* 2.0 (* c (+ a 0.8333333333333334)))))))
t_1)))))
double code(double x, double y, double z, double t, double a, double b, double c) {
double t_1 = x / (x + (y * exp((2.0 * (b * ((0.6666666666666666 / t) - (a + 0.8333333333333334)))))));
double tmp;
if (b <= -1.35e+39) {
tmp = t_1;
} else if (b <= -2.45e-205) {
tmp = x / (x + (y * exp((2.0 * ((c * -0.6666666666666666) / t)))));
} else if (b <= 1.45e-98) {
tmp = x / (x + (y * exp((2.0 * (c * (a + 0.8333333333333334))))));
} else {
tmp = t_1;
}
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) :: tmp
t_1 = x / (x + (y * exp((2.0d0 * (b * ((0.6666666666666666d0 / t) - (a + 0.8333333333333334d0)))))))
if (b <= (-1.35d+39)) then
tmp = t_1
else if (b <= (-2.45d-205)) then
tmp = x / (x + (y * exp((2.0d0 * ((c * (-0.6666666666666666d0)) / t)))))
else if (b <= 1.45d-98) then
tmp = x / (x + (y * exp((2.0d0 * (c * (a + 0.8333333333333334d0))))))
else
tmp = t_1
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 = x / (x + (y * Math.exp((2.0 * (b * ((0.6666666666666666 / t) - (a + 0.8333333333333334)))))));
double tmp;
if (b <= -1.35e+39) {
tmp = t_1;
} else if (b <= -2.45e-205) {
tmp = x / (x + (y * Math.exp((2.0 * ((c * -0.6666666666666666) / t)))));
} else if (b <= 1.45e-98) {
tmp = x / (x + (y * Math.exp((2.0 * (c * (a + 0.8333333333333334))))));
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b, c): t_1 = x / (x + (y * math.exp((2.0 * (b * ((0.6666666666666666 / t) - (a + 0.8333333333333334))))))) tmp = 0 if b <= -1.35e+39: tmp = t_1 elif b <= -2.45e-205: tmp = x / (x + (y * math.exp((2.0 * ((c * -0.6666666666666666) / t))))) elif b <= 1.45e-98: tmp = x / (x + (y * math.exp((2.0 * (c * (a + 0.8333333333333334)))))) else: tmp = t_1 return tmp
function code(x, y, z, t, a, b, c) t_1 = Float64(x / Float64(x + Float64(y * exp(Float64(2.0 * Float64(b * Float64(Float64(0.6666666666666666 / t) - Float64(a + 0.8333333333333334)))))))) tmp = 0.0 if (b <= -1.35e+39) tmp = t_1; elseif (b <= -2.45e-205) tmp = Float64(x / Float64(x + Float64(y * exp(Float64(2.0 * Float64(Float64(c * -0.6666666666666666) / t)))))); elseif (b <= 1.45e-98) tmp = Float64(x / Float64(x + Float64(y * exp(Float64(2.0 * Float64(c * Float64(a + 0.8333333333333334))))))); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c) t_1 = x / (x + (y * exp((2.0 * (b * ((0.6666666666666666 / t) - (a + 0.8333333333333334))))))); tmp = 0.0; if (b <= -1.35e+39) tmp = t_1; elseif (b <= -2.45e-205) tmp = x / (x + (y * exp((2.0 * ((c * -0.6666666666666666) / t))))); elseif (b <= 1.45e-98) tmp = x / (x + (y * exp((2.0 * (c * (a + 0.8333333333333334)))))); else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_] := Block[{t$95$1 = N[(x / N[(x + N[(y * N[Exp[N[(2.0 * N[(b * N[(N[(0.6666666666666666 / t), $MachinePrecision] - N[(a + 0.8333333333333334), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[b, -1.35e+39], t$95$1, If[LessEqual[b, -2.45e-205], N[(x / N[(x + N[(y * N[Exp[N[(2.0 * N[(N[(c * -0.6666666666666666), $MachinePrecision] / t), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[b, 1.45e-98], N[(x / N[(x + N[(y * N[Exp[N[(2.0 * N[(c * N[(a + 0.8333333333333334), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$1]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{x}{x + y \cdot e^{2 \cdot \left(b \cdot \left(\frac{0.6666666666666666}{t} - \left(a + 0.8333333333333334\right)\right)\right)}}\\
\mathbf{if}\;b \leq -1.35 \cdot 10^{+39}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;b \leq -2.45 \cdot 10^{-205}:\\
\;\;\;\;\frac{x}{x + y \cdot e^{2 \cdot \frac{c \cdot -0.6666666666666666}{t}}}\\
\mathbf{elif}\;b \leq 1.45 \cdot 10^{-98}:\\
\;\;\;\;\frac{x}{x + y \cdot e^{2 \cdot \left(c \cdot \left(a + 0.8333333333333334\right)\right)}}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if b < -1.35000000000000002e39 or 1.45e-98 < b Initial program 86.7%
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
count-2N/A
exp-lowering-exp.f64N/A
Simplified88.9%
Taylor expanded in b around inf
*-lowering-*.f64N/A
--lowering--.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f64N/A
+-lowering-+.f6486.4%
Simplified86.4%
if -1.35000000000000002e39 < b < -2.4499999999999999e-205Initial program 93.8%
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
count-2N/A
exp-lowering-exp.f64N/A
Simplified95.8%
Taylor expanded in c around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-lowering-*.f64N/A
--lowering--.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f64N/A
+-lowering-+.f6477.9%
Simplified77.9%
Taylor expanded in t around 0
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f64N/A
metadata-evalN/A
*-lowering-*.f6469.8%
Simplified69.8%
if -2.4499999999999999e-205 < b < 1.45e-98Initial program 93.3%
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
count-2N/A
exp-lowering-exp.f64N/A
Simplified97.2%
Taylor expanded in c around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-lowering-*.f64N/A
--lowering--.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f64N/A
+-lowering-+.f6479.7%
Simplified79.7%
Taylor expanded in t around inf
*-lowering-*.f64N/A
+-lowering-+.f6478.5%
Simplified78.5%
Final simplification81.0%
(FPCore (x y z t a b c)
:precision binary64
(let* ((t_1 (- (/ 0.6666666666666666 t) (+ a 0.8333333333333334)))
(t_2 (* t_1 t_1)))
(if (<= b -1.7e+42)
(/
x
(+
x
(*
y
(+
(*
b
(+
(* 2.0 t_1)
(* b (+ (* (* b 1.3333333333333333) (* t_1 t_2)) (* 2.0 t_2)))))
1.0))))
(if (<= b -1.12e-204)
(/ x (+ x (* y (exp (* 2.0 (/ (* c -0.6666666666666666) t))))))
(if (<= b 1.6e-56)
(/ x (+ x (* y (exp (* 2.0 (* c (+ a 0.8333333333333334)))))))
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 = t_1 * t_1;
double tmp;
if (b <= -1.7e+42) {
tmp = x / (x + (y * ((b * ((2.0 * t_1) + (b * (((b * 1.3333333333333333) * (t_1 * t_2)) + (2.0 * t_2))))) + 1.0)));
} else if (b <= -1.12e-204) {
tmp = x / (x + (y * exp((2.0 * ((c * -0.6666666666666666) / t)))));
} else if (b <= 1.6e-56) {
tmp = x / (x + (y * exp((2.0 * (c * (a + 0.8333333333333334))))));
} 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) :: t_1
real(8) :: t_2
real(8) :: tmp
t_1 = (0.6666666666666666d0 / t) - (a + 0.8333333333333334d0)
t_2 = t_1 * t_1
if (b <= (-1.7d+42)) then
tmp = x / (x + (y * ((b * ((2.0d0 * t_1) + (b * (((b * 1.3333333333333333d0) * (t_1 * t_2)) + (2.0d0 * t_2))))) + 1.0d0)))
else if (b <= (-1.12d-204)) then
tmp = x / (x + (y * exp((2.0d0 * ((c * (-0.6666666666666666d0)) / t)))))
else if (b <= 1.6d-56) then
tmp = x / (x + (y * exp((2.0d0 * (c * (a + 0.8333333333333334d0))))))
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 t_1 = (0.6666666666666666 / t) - (a + 0.8333333333333334);
double t_2 = t_1 * t_1;
double tmp;
if (b <= -1.7e+42) {
tmp = x / (x + (y * ((b * ((2.0 * t_1) + (b * (((b * 1.3333333333333333) * (t_1 * t_2)) + (2.0 * t_2))))) + 1.0)));
} else if (b <= -1.12e-204) {
tmp = x / (x + (y * Math.exp((2.0 * ((c * -0.6666666666666666) / t)))));
} else if (b <= 1.6e-56) {
tmp = x / (x + (y * Math.exp((2.0 * (c * (a + 0.8333333333333334))))));
} else {
tmp = 1.0;
}
return tmp;
}
def code(x, y, z, t, a, b, c): t_1 = (0.6666666666666666 / t) - (a + 0.8333333333333334) t_2 = t_1 * t_1 tmp = 0 if b <= -1.7e+42: tmp = x / (x + (y * ((b * ((2.0 * t_1) + (b * (((b * 1.3333333333333333) * (t_1 * t_2)) + (2.0 * t_2))))) + 1.0))) elif b <= -1.12e-204: tmp = x / (x + (y * math.exp((2.0 * ((c * -0.6666666666666666) / t))))) elif b <= 1.6e-56: tmp = x / (x + (y * math.exp((2.0 * (c * (a + 0.8333333333333334)))))) 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)) t_2 = Float64(t_1 * t_1) tmp = 0.0 if (b <= -1.7e+42) tmp = Float64(x / Float64(x + Float64(y * Float64(Float64(b * Float64(Float64(2.0 * t_1) + Float64(b * Float64(Float64(Float64(b * 1.3333333333333333) * Float64(t_1 * t_2)) + Float64(2.0 * t_2))))) + 1.0)))); elseif (b <= -1.12e-204) tmp = Float64(x / Float64(x + Float64(y * exp(Float64(2.0 * Float64(Float64(c * -0.6666666666666666) / t)))))); elseif (b <= 1.6e-56) tmp = Float64(x / Float64(x + Float64(y * exp(Float64(2.0 * Float64(c * Float64(a + 0.8333333333333334))))))); else tmp = 1.0; end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c) t_1 = (0.6666666666666666 / t) - (a + 0.8333333333333334); t_2 = t_1 * t_1; tmp = 0.0; if (b <= -1.7e+42) tmp = x / (x + (y * ((b * ((2.0 * t_1) + (b * (((b * 1.3333333333333333) * (t_1 * t_2)) + (2.0 * t_2))))) + 1.0))); elseif (b <= -1.12e-204) tmp = x / (x + (y * exp((2.0 * ((c * -0.6666666666666666) / t))))); elseif (b <= 1.6e-56) tmp = x / (x + (y * exp((2.0 * (c * (a + 0.8333333333333334)))))); else tmp = 1.0; end tmp_2 = 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[(t$95$1 * t$95$1), $MachinePrecision]}, If[LessEqual[b, -1.7e+42], N[(x / N[(x + N[(y * N[(N[(b * N[(N[(2.0 * t$95$1), $MachinePrecision] + N[(b * N[(N[(N[(b * 1.3333333333333333), $MachinePrecision] * N[(t$95$1 * t$95$2), $MachinePrecision]), $MachinePrecision] + N[(2.0 * t$95$2), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[b, -1.12e-204], N[(x / N[(x + N[(y * N[Exp[N[(2.0 * N[(N[(c * -0.6666666666666666), $MachinePrecision] / t), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[b, 1.6e-56], N[(x / N[(x + N[(y * N[Exp[N[(2.0 * N[(c * N[(a + 0.8333333333333334), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 1.0]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{0.6666666666666666}{t} - \left(a + 0.8333333333333334\right)\\
t_2 := t\_1 \cdot t\_1\\
\mathbf{if}\;b \leq -1.7 \cdot 10^{+42}:\\
\;\;\;\;\frac{x}{x + y \cdot \left(b \cdot \left(2 \cdot t\_1 + b \cdot \left(\left(b \cdot 1.3333333333333333\right) \cdot \left(t\_1 \cdot t\_2\right) + 2 \cdot t\_2\right)\right) + 1\right)}\\
\mathbf{elif}\;b \leq -1.12 \cdot 10^{-204}:\\
\;\;\;\;\frac{x}{x + y \cdot e^{2 \cdot \frac{c \cdot -0.6666666666666666}{t}}}\\
\mathbf{elif}\;b \leq 1.6 \cdot 10^{-56}:\\
\;\;\;\;\frac{x}{x + y \cdot e^{2 \cdot \left(c \cdot \left(a + 0.8333333333333334\right)\right)}}\\
\mathbf{else}:\\
\;\;\;\;1\\
\end{array}
\end{array}
if b < -1.69999999999999988e42Initial program 83.7%
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
count-2N/A
exp-lowering-exp.f64N/A
Simplified85.5%
Taylor expanded in b around inf
*-lowering-*.f64N/A
--lowering--.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f64N/A
+-lowering-+.f6485.9%
Simplified85.9%
Taylor expanded in b around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
Simplified68.0%
if -1.69999999999999988e42 < b < -1.11999999999999997e-204Initial program 93.8%
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
count-2N/A
exp-lowering-exp.f64N/A
Simplified95.8%
Taylor expanded in c around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-lowering-*.f64N/A
--lowering--.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f64N/A
+-lowering-+.f6477.9%
Simplified77.9%
Taylor expanded in t around 0
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f64N/A
metadata-evalN/A
*-lowering-*.f6469.8%
Simplified69.8%
if -1.11999999999999997e-204 < b < 1.59999999999999993e-56Initial program 94.0%
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
count-2N/A
exp-lowering-exp.f64N/A
Simplified97.5%
Taylor expanded in c around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-lowering-*.f64N/A
--lowering--.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f64N/A
+-lowering-+.f6479.6%
Simplified79.6%
Taylor expanded in t around inf
*-lowering-*.f64N/A
+-lowering-+.f6474.9%
Simplified74.9%
if 1.59999999999999993e-56 < b Initial program 87.4%
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
count-2N/A
exp-lowering-exp.f64N/A
Simplified90.1%
Taylor expanded in x around inf
Simplified61.8%
Final simplification68.8%
(FPCore (x y z t a b c)
:precision binary64
(let* ((t_1 (- (/ 0.6666666666666666 t) (+ a 0.8333333333333334)))
(t_2 (* t_1 t_1)))
(if (<= b -5.5e+61)
(/
x
(+
x
(*
y
(+
(*
b
(+
(* 2.0 t_1)
(* b (+ (* (* b 1.3333333333333333) (* t_1 t_2)) (* 2.0 t_2)))))
1.0))))
(if (<= b -2.7e-128)
(/
x
(+
x
(*
(* b b)
(+
(*
y
(+
(* -2.6666666666666665 (/ (+ a 0.8333333333333334) t))
(/ 0.8888888888888888 (* t t))))
(/ (- (/ y b) (/ (* y -1.3333333333333333) t)) b)))))
(if (<= b 2.35e-57)
(/ x (+ x (* y (exp (* 2.0 (* c (+ a 0.8333333333333334)))))))
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 = t_1 * t_1;
double tmp;
if (b <= -5.5e+61) {
tmp = x / (x + (y * ((b * ((2.0 * t_1) + (b * (((b * 1.3333333333333333) * (t_1 * t_2)) + (2.0 * t_2))))) + 1.0)));
} else if (b <= -2.7e-128) {
tmp = x / (x + ((b * b) * ((y * ((-2.6666666666666665 * ((a + 0.8333333333333334) / t)) + (0.8888888888888888 / (t * t)))) + (((y / b) - ((y * -1.3333333333333333) / t)) / b))));
} else if (b <= 2.35e-57) {
tmp = x / (x + (y * exp((2.0 * (c * (a + 0.8333333333333334))))));
} 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) :: t_1
real(8) :: t_2
real(8) :: tmp
t_1 = (0.6666666666666666d0 / t) - (a + 0.8333333333333334d0)
t_2 = t_1 * t_1
if (b <= (-5.5d+61)) then
tmp = x / (x + (y * ((b * ((2.0d0 * t_1) + (b * (((b * 1.3333333333333333d0) * (t_1 * t_2)) + (2.0d0 * t_2))))) + 1.0d0)))
else if (b <= (-2.7d-128)) then
tmp = x / (x + ((b * b) * ((y * (((-2.6666666666666665d0) * ((a + 0.8333333333333334d0) / t)) + (0.8888888888888888d0 / (t * t)))) + (((y / b) - ((y * (-1.3333333333333333d0)) / t)) / b))))
else if (b <= 2.35d-57) then
tmp = x / (x + (y * exp((2.0d0 * (c * (a + 0.8333333333333334d0))))))
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 t_1 = (0.6666666666666666 / t) - (a + 0.8333333333333334);
double t_2 = t_1 * t_1;
double tmp;
if (b <= -5.5e+61) {
tmp = x / (x + (y * ((b * ((2.0 * t_1) + (b * (((b * 1.3333333333333333) * (t_1 * t_2)) + (2.0 * t_2))))) + 1.0)));
} else if (b <= -2.7e-128) {
tmp = x / (x + ((b * b) * ((y * ((-2.6666666666666665 * ((a + 0.8333333333333334) / t)) + (0.8888888888888888 / (t * t)))) + (((y / b) - ((y * -1.3333333333333333) / t)) / b))));
} else if (b <= 2.35e-57) {
tmp = x / (x + (y * Math.exp((2.0 * (c * (a + 0.8333333333333334))))));
} else {
tmp = 1.0;
}
return tmp;
}
def code(x, y, z, t, a, b, c): t_1 = (0.6666666666666666 / t) - (a + 0.8333333333333334) t_2 = t_1 * t_1 tmp = 0 if b <= -5.5e+61: tmp = x / (x + (y * ((b * ((2.0 * t_1) + (b * (((b * 1.3333333333333333) * (t_1 * t_2)) + (2.0 * t_2))))) + 1.0))) elif b <= -2.7e-128: tmp = x / (x + ((b * b) * ((y * ((-2.6666666666666665 * ((a + 0.8333333333333334) / t)) + (0.8888888888888888 / (t * t)))) + (((y / b) - ((y * -1.3333333333333333) / t)) / b)))) elif b <= 2.35e-57: tmp = x / (x + (y * math.exp((2.0 * (c * (a + 0.8333333333333334)))))) 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)) t_2 = Float64(t_1 * t_1) tmp = 0.0 if (b <= -5.5e+61) tmp = Float64(x / Float64(x + Float64(y * Float64(Float64(b * Float64(Float64(2.0 * t_1) + Float64(b * Float64(Float64(Float64(b * 1.3333333333333333) * Float64(t_1 * t_2)) + Float64(2.0 * t_2))))) + 1.0)))); elseif (b <= -2.7e-128) tmp = Float64(x / Float64(x + Float64(Float64(b * b) * Float64(Float64(y * Float64(Float64(-2.6666666666666665 * Float64(Float64(a + 0.8333333333333334) / t)) + Float64(0.8888888888888888 / Float64(t * t)))) + Float64(Float64(Float64(y / b) - Float64(Float64(y * -1.3333333333333333) / t)) / b))))); elseif (b <= 2.35e-57) tmp = Float64(x / Float64(x + Float64(y * exp(Float64(2.0 * Float64(c * Float64(a + 0.8333333333333334))))))); else tmp = 1.0; end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c) t_1 = (0.6666666666666666 / t) - (a + 0.8333333333333334); t_2 = t_1 * t_1; tmp = 0.0; if (b <= -5.5e+61) tmp = x / (x + (y * ((b * ((2.0 * t_1) + (b * (((b * 1.3333333333333333) * (t_1 * t_2)) + (2.0 * t_2))))) + 1.0))); elseif (b <= -2.7e-128) tmp = x / (x + ((b * b) * ((y * ((-2.6666666666666665 * ((a + 0.8333333333333334) / t)) + (0.8888888888888888 / (t * t)))) + (((y / b) - ((y * -1.3333333333333333) / t)) / b)))); elseif (b <= 2.35e-57) tmp = x / (x + (y * exp((2.0 * (c * (a + 0.8333333333333334)))))); else tmp = 1.0; end tmp_2 = 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[(t$95$1 * t$95$1), $MachinePrecision]}, If[LessEqual[b, -5.5e+61], N[(x / N[(x + N[(y * N[(N[(b * N[(N[(2.0 * t$95$1), $MachinePrecision] + N[(b * N[(N[(N[(b * 1.3333333333333333), $MachinePrecision] * N[(t$95$1 * t$95$2), $MachinePrecision]), $MachinePrecision] + N[(2.0 * t$95$2), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[b, -2.7e-128], N[(x / N[(x + N[(N[(b * b), $MachinePrecision] * N[(N[(y * N[(N[(-2.6666666666666665 * N[(N[(a + 0.8333333333333334), $MachinePrecision] / t), $MachinePrecision]), $MachinePrecision] + N[(0.8888888888888888 / N[(t * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(N[(N[(y / b), $MachinePrecision] - N[(N[(y * -1.3333333333333333), $MachinePrecision] / t), $MachinePrecision]), $MachinePrecision] / b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[b, 2.35e-57], N[(x / N[(x + N[(y * N[Exp[N[(2.0 * N[(c * N[(a + 0.8333333333333334), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 1.0]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{0.6666666666666666}{t} - \left(a + 0.8333333333333334\right)\\
t_2 := t\_1 \cdot t\_1\\
\mathbf{if}\;b \leq -5.5 \cdot 10^{+61}:\\
\;\;\;\;\frac{x}{x + y \cdot \left(b \cdot \left(2 \cdot t\_1 + b \cdot \left(\left(b \cdot 1.3333333333333333\right) \cdot \left(t\_1 \cdot t\_2\right) + 2 \cdot t\_2\right)\right) + 1\right)}\\
\mathbf{elif}\;b \leq -2.7 \cdot 10^{-128}:\\
\;\;\;\;\frac{x}{x + \left(b \cdot b\right) \cdot \left(y \cdot \left(-2.6666666666666665 \cdot \frac{a + 0.8333333333333334}{t} + \frac{0.8888888888888888}{t \cdot t}\right) + \frac{\frac{y}{b} - \frac{y \cdot -1.3333333333333333}{t}}{b}\right)}\\
\mathbf{elif}\;b \leq 2.35 \cdot 10^{-57}:\\
\;\;\;\;\frac{x}{x + y \cdot e^{2 \cdot \left(c \cdot \left(a + 0.8333333333333334\right)\right)}}\\
\mathbf{else}:\\
\;\;\;\;1\\
\end{array}
\end{array}
if b < -5.50000000000000036e61Initial program 83.1%
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
count-2N/A
exp-lowering-exp.f64N/A
Simplified84.9%
Taylor expanded in b around inf
*-lowering-*.f64N/A
--lowering--.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f64N/A
+-lowering-+.f6485.4%
Simplified85.4%
Taylor expanded in b around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
Simplified68.6%
if -5.50000000000000036e61 < b < -2.70000000000000006e-128Initial program 94.9%
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
count-2N/A
exp-lowering-exp.f64N/A
Simplified97.4%
Taylor expanded in b around inf
*-lowering-*.f64N/A
--lowering--.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f64N/A
+-lowering-+.f6450.4%
Simplified50.4%
Taylor expanded in b around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
distribute-lft-outN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
Simplified52.9%
Taylor expanded in t around 0
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
unpow2N/A
*-lowering-*.f6450.2%
Simplified50.2%
Taylor expanded in b around -inf
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
Simplified62.3%
if -2.70000000000000006e-128 < b < 2.3499999999999999e-57Initial program 93.6%
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
count-2N/A
exp-lowering-exp.f64N/A
Simplified96.8%
Taylor expanded in c around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-lowering-*.f64N/A
--lowering--.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f64N/A
+-lowering-+.f6479.9%
Simplified79.9%
Taylor expanded in t around inf
*-lowering-*.f64N/A
+-lowering-+.f6472.7%
Simplified72.7%
if 2.3499999999999999e-57 < b Initial program 87.4%
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
count-2N/A
exp-lowering-exp.f64N/A
Simplified90.1%
Taylor expanded in x around inf
Simplified61.8%
Final simplification67.3%
(FPCore (x y z t a b c)
:precision binary64
(let* ((t_1
(/
x
(+
x
(*
y
(exp
(*
2.0
(*
b
(- (/ 0.6666666666666666 t) (+ a 0.8333333333333334))))))))))
(if (<= b -2e+39)
t_1
(if (<= b 4.2e-76)
(/
x
(+
x
(*
y
(exp
(*
2.0
(* c (- (+ a 0.8333333333333334) (/ 0.6666666666666666 t))))))))
t_1))))
double code(double x, double y, double z, double t, double a, double b, double c) {
double t_1 = x / (x + (y * exp((2.0 * (b * ((0.6666666666666666 / t) - (a + 0.8333333333333334)))))));
double tmp;
if (b <= -2e+39) {
tmp = t_1;
} else if (b <= 4.2e-76) {
tmp = x / (x + (y * exp((2.0 * (c * ((a + 0.8333333333333334) - (0.6666666666666666 / t)))))));
} else {
tmp = t_1;
}
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) :: tmp
t_1 = x / (x + (y * exp((2.0d0 * (b * ((0.6666666666666666d0 / t) - (a + 0.8333333333333334d0)))))))
if (b <= (-2d+39)) then
tmp = t_1
else if (b <= 4.2d-76) then
tmp = x / (x + (y * exp((2.0d0 * (c * ((a + 0.8333333333333334d0) - (0.6666666666666666d0 / t)))))))
else
tmp = t_1
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 = x / (x + (y * Math.exp((2.0 * (b * ((0.6666666666666666 / t) - (a + 0.8333333333333334)))))));
double tmp;
if (b <= -2e+39) {
tmp = t_1;
} else if (b <= 4.2e-76) {
tmp = x / (x + (y * Math.exp((2.0 * (c * ((a + 0.8333333333333334) - (0.6666666666666666 / t)))))));
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b, c): t_1 = x / (x + (y * math.exp((2.0 * (b * ((0.6666666666666666 / t) - (a + 0.8333333333333334))))))) tmp = 0 if b <= -2e+39: tmp = t_1 elif b <= 4.2e-76: tmp = x / (x + (y * math.exp((2.0 * (c * ((a + 0.8333333333333334) - (0.6666666666666666 / t))))))) else: tmp = t_1 return tmp
function code(x, y, z, t, a, b, c) t_1 = Float64(x / Float64(x + Float64(y * exp(Float64(2.0 * Float64(b * Float64(Float64(0.6666666666666666 / t) - Float64(a + 0.8333333333333334)))))))) tmp = 0.0 if (b <= -2e+39) tmp = t_1; elseif (b <= 4.2e-76) tmp = Float64(x / Float64(x + Float64(y * exp(Float64(2.0 * Float64(c * Float64(Float64(a + 0.8333333333333334) - Float64(0.6666666666666666 / t)))))))); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c) t_1 = x / (x + (y * exp((2.0 * (b * ((0.6666666666666666 / t) - (a + 0.8333333333333334))))))); tmp = 0.0; if (b <= -2e+39) tmp = t_1; elseif (b <= 4.2e-76) tmp = x / (x + (y * exp((2.0 * (c * ((a + 0.8333333333333334) - (0.6666666666666666 / t))))))); else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_] := Block[{t$95$1 = N[(x / N[(x + N[(y * N[Exp[N[(2.0 * N[(b * N[(N[(0.6666666666666666 / t), $MachinePrecision] - N[(a + 0.8333333333333334), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[b, -2e+39], t$95$1, If[LessEqual[b, 4.2e-76], N[(x / N[(x + N[(y * N[Exp[N[(2.0 * N[(c * N[(N[(a + 0.8333333333333334), $MachinePrecision] - N[(0.6666666666666666 / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{x}{x + y \cdot e^{2 \cdot \left(b \cdot \left(\frac{0.6666666666666666}{t} - \left(a + 0.8333333333333334\right)\right)\right)}}\\
\mathbf{if}\;b \leq -2 \cdot 10^{+39}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;b \leq 4.2 \cdot 10^{-76}:\\
\;\;\;\;\frac{x}{x + y \cdot e^{2 \cdot \left(c \cdot \left(\left(a + 0.8333333333333334\right) - \frac{0.6666666666666666}{t}\right)\right)}}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if b < -1.99999999999999988e39 or 4.19999999999999985e-76 < b Initial program 86.2%
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
count-2N/A
exp-lowering-exp.f64N/A
Simplified88.5%
Taylor expanded in b around inf
*-lowering-*.f64N/A
--lowering--.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f64N/A
+-lowering-+.f6487.3%
Simplified87.3%
if -1.99999999999999988e39 < b < 4.19999999999999985e-76Initial program 93.7%
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
count-2N/A
exp-lowering-exp.f64N/A
Simplified96.8%
Taylor expanded in c around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-lowering-*.f64N/A
--lowering--.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f64N/A
+-lowering-+.f6478.3%
Simplified78.3%
Final simplification82.9%
(FPCore (x y z t a b c)
:precision binary64
(let* ((t_1 (- (/ 0.6666666666666666 t) (+ a 0.8333333333333334)))
(t_2 (* t_1 t_1)))
(if (<= b -1.45e+64)
(/
x
(+
x
(*
y
(+
(*
b
(+
(* 2.0 t_1)
(* b (+ (* (* b 1.3333333333333333) (* t_1 t_2)) (* 2.0 t_2)))))
1.0))))
(if (<= b -2.7e-129)
(/
x
(+
x
(*
(* b b)
(+
(*
y
(+
(* -2.6666666666666665 (/ (+ a 0.8333333333333334) t))
(/ 0.8888888888888888 (* t t))))
(/ (- (/ y b) (/ (* y -1.3333333333333333) t)) b)))))
(if (<= b 5.4e-83) (/ x (+ x (* y (exp (* 2.0 (* a c)))))) 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 = t_1 * t_1;
double tmp;
if (b <= -1.45e+64) {
tmp = x / (x + (y * ((b * ((2.0 * t_1) + (b * (((b * 1.3333333333333333) * (t_1 * t_2)) + (2.0 * t_2))))) + 1.0)));
} else if (b <= -2.7e-129) {
tmp = x / (x + ((b * b) * ((y * ((-2.6666666666666665 * ((a + 0.8333333333333334) / t)) + (0.8888888888888888 / (t * t)))) + (((y / b) - ((y * -1.3333333333333333) / t)) / b))));
} else if (b <= 5.4e-83) {
tmp = x / (x + (y * exp((2.0 * (a * c)))));
} 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) :: t_1
real(8) :: t_2
real(8) :: tmp
t_1 = (0.6666666666666666d0 / t) - (a + 0.8333333333333334d0)
t_2 = t_1 * t_1
if (b <= (-1.45d+64)) then
tmp = x / (x + (y * ((b * ((2.0d0 * t_1) + (b * (((b * 1.3333333333333333d0) * (t_1 * t_2)) + (2.0d0 * t_2))))) + 1.0d0)))
else if (b <= (-2.7d-129)) then
tmp = x / (x + ((b * b) * ((y * (((-2.6666666666666665d0) * ((a + 0.8333333333333334d0) / t)) + (0.8888888888888888d0 / (t * t)))) + (((y / b) - ((y * (-1.3333333333333333d0)) / t)) / b))))
else if (b <= 5.4d-83) then
tmp = x / (x + (y * exp((2.0d0 * (a * c)))))
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 t_1 = (0.6666666666666666 / t) - (a + 0.8333333333333334);
double t_2 = t_1 * t_1;
double tmp;
if (b <= -1.45e+64) {
tmp = x / (x + (y * ((b * ((2.0 * t_1) + (b * (((b * 1.3333333333333333) * (t_1 * t_2)) + (2.0 * t_2))))) + 1.0)));
} else if (b <= -2.7e-129) {
tmp = x / (x + ((b * b) * ((y * ((-2.6666666666666665 * ((a + 0.8333333333333334) / t)) + (0.8888888888888888 / (t * t)))) + (((y / b) - ((y * -1.3333333333333333) / t)) / b))));
} else if (b <= 5.4e-83) {
tmp = x / (x + (y * Math.exp((2.0 * (a * c)))));
} else {
tmp = 1.0;
}
return tmp;
}
def code(x, y, z, t, a, b, c): t_1 = (0.6666666666666666 / t) - (a + 0.8333333333333334) t_2 = t_1 * t_1 tmp = 0 if b <= -1.45e+64: tmp = x / (x + (y * ((b * ((2.0 * t_1) + (b * (((b * 1.3333333333333333) * (t_1 * t_2)) + (2.0 * t_2))))) + 1.0))) elif b <= -2.7e-129: tmp = x / (x + ((b * b) * ((y * ((-2.6666666666666665 * ((a + 0.8333333333333334) / t)) + (0.8888888888888888 / (t * t)))) + (((y / b) - ((y * -1.3333333333333333) / t)) / b)))) elif b <= 5.4e-83: tmp = x / (x + (y * math.exp((2.0 * (a * c))))) 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)) t_2 = Float64(t_1 * t_1) tmp = 0.0 if (b <= -1.45e+64) tmp = Float64(x / Float64(x + Float64(y * Float64(Float64(b * Float64(Float64(2.0 * t_1) + Float64(b * Float64(Float64(Float64(b * 1.3333333333333333) * Float64(t_1 * t_2)) + Float64(2.0 * t_2))))) + 1.0)))); elseif (b <= -2.7e-129) tmp = Float64(x / Float64(x + Float64(Float64(b * b) * Float64(Float64(y * Float64(Float64(-2.6666666666666665 * Float64(Float64(a + 0.8333333333333334) / t)) + Float64(0.8888888888888888 / Float64(t * t)))) + Float64(Float64(Float64(y / b) - Float64(Float64(y * -1.3333333333333333) / t)) / b))))); elseif (b <= 5.4e-83) tmp = Float64(x / Float64(x + Float64(y * exp(Float64(2.0 * Float64(a * c)))))); else tmp = 1.0; end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c) t_1 = (0.6666666666666666 / t) - (a + 0.8333333333333334); t_2 = t_1 * t_1; tmp = 0.0; if (b <= -1.45e+64) tmp = x / (x + (y * ((b * ((2.0 * t_1) + (b * (((b * 1.3333333333333333) * (t_1 * t_2)) + (2.0 * t_2))))) + 1.0))); elseif (b <= -2.7e-129) tmp = x / (x + ((b * b) * ((y * ((-2.6666666666666665 * ((a + 0.8333333333333334) / t)) + (0.8888888888888888 / (t * t)))) + (((y / b) - ((y * -1.3333333333333333) / t)) / b)))); elseif (b <= 5.4e-83) tmp = x / (x + (y * exp((2.0 * (a * c))))); else tmp = 1.0; end tmp_2 = 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[(t$95$1 * t$95$1), $MachinePrecision]}, If[LessEqual[b, -1.45e+64], N[(x / N[(x + N[(y * N[(N[(b * N[(N[(2.0 * t$95$1), $MachinePrecision] + N[(b * N[(N[(N[(b * 1.3333333333333333), $MachinePrecision] * N[(t$95$1 * t$95$2), $MachinePrecision]), $MachinePrecision] + N[(2.0 * t$95$2), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[b, -2.7e-129], N[(x / N[(x + N[(N[(b * b), $MachinePrecision] * N[(N[(y * N[(N[(-2.6666666666666665 * N[(N[(a + 0.8333333333333334), $MachinePrecision] / t), $MachinePrecision]), $MachinePrecision] + N[(0.8888888888888888 / N[(t * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(N[(N[(y / b), $MachinePrecision] - N[(N[(y * -1.3333333333333333), $MachinePrecision] / t), $MachinePrecision]), $MachinePrecision] / b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[b, 5.4e-83], N[(x / N[(x + N[(y * N[Exp[N[(2.0 * N[(a * c), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 1.0]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{0.6666666666666666}{t} - \left(a + 0.8333333333333334\right)\\
t_2 := t\_1 \cdot t\_1\\
\mathbf{if}\;b \leq -1.45 \cdot 10^{+64}:\\
\;\;\;\;\frac{x}{x + y \cdot \left(b \cdot \left(2 \cdot t\_1 + b \cdot \left(\left(b \cdot 1.3333333333333333\right) \cdot \left(t\_1 \cdot t\_2\right) + 2 \cdot t\_2\right)\right) + 1\right)}\\
\mathbf{elif}\;b \leq -2.7 \cdot 10^{-129}:\\
\;\;\;\;\frac{x}{x + \left(b \cdot b\right) \cdot \left(y \cdot \left(-2.6666666666666665 \cdot \frac{a + 0.8333333333333334}{t} + \frac{0.8888888888888888}{t \cdot t}\right) + \frac{\frac{y}{b} - \frac{y \cdot -1.3333333333333333}{t}}{b}\right)}\\
\mathbf{elif}\;b \leq 5.4 \cdot 10^{-83}:\\
\;\;\;\;\frac{x}{x + y \cdot e^{2 \cdot \left(a \cdot c\right)}}\\
\mathbf{else}:\\
\;\;\;\;1\\
\end{array}
\end{array}
if b < -1.44999999999999997e64Initial program 83.1%
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
count-2N/A
exp-lowering-exp.f64N/A
Simplified84.9%
Taylor expanded in b around inf
*-lowering-*.f64N/A
--lowering--.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f64N/A
+-lowering-+.f6485.4%
Simplified85.4%
Taylor expanded in b around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
Simplified68.6%
if -1.44999999999999997e64 < b < -2.69999999999999999e-129Initial program 94.9%
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
count-2N/A
exp-lowering-exp.f64N/A
Simplified97.4%
Taylor expanded in b around inf
*-lowering-*.f64N/A
--lowering--.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f64N/A
+-lowering-+.f6450.4%
Simplified50.4%
Taylor expanded in b around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
distribute-lft-outN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
Simplified52.9%
Taylor expanded in t around 0
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
unpow2N/A
*-lowering-*.f6450.2%
Simplified50.2%
Taylor expanded in b around -inf
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
Simplified62.3%
if -2.69999999999999999e-129 < b < 5.39999999999999982e-83Initial program 93.3%
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
count-2N/A
exp-lowering-exp.f64N/A
Simplified96.6%
Taylor expanded in c around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-lowering-*.f64N/A
--lowering--.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f64N/A
+-lowering-+.f6478.8%
Simplified78.8%
Taylor expanded in a around inf
*-lowering-*.f6465.7%
Simplified65.7%
if 5.39999999999999982e-83 < b Initial program 88.2%
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
count-2N/A
exp-lowering-exp.f64N/A
Simplified90.8%
Taylor expanded in x around inf
Simplified61.8%
Final simplification64.6%
(FPCore (x y z t a b c)
:precision binary64
(if (<= t 2.25e-187)
(/ x (+ x (* y (exp (* 2.0 (/ (* b 0.6666666666666666) t))))))
(if (<= t 1.56e-21)
(/ x (+ x (* y (exp (* 2.0 (/ (* c -0.6666666666666666) t))))))
(/ x (+ x (* y (exp (* 2.0 (* c (+ a 0.8333333333333334))))))))))
double code(double x, double y, double z, double t, double a, double b, double c) {
double tmp;
if (t <= 2.25e-187) {
tmp = x / (x + (y * exp((2.0 * ((b * 0.6666666666666666) / t)))));
} else if (t <= 1.56e-21) {
tmp = x / (x + (y * exp((2.0 * ((c * -0.6666666666666666) / t)))));
} else {
tmp = x / (x + (y * exp((2.0 * (c * (a + 0.8333333333333334))))));
}
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 (t <= 2.25d-187) then
tmp = x / (x + (y * exp((2.0d0 * ((b * 0.6666666666666666d0) / t)))))
else if (t <= 1.56d-21) then
tmp = x / (x + (y * exp((2.0d0 * ((c * (-0.6666666666666666d0)) / t)))))
else
tmp = x / (x + (y * exp((2.0d0 * (c * (a + 0.8333333333333334d0))))))
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 (t <= 2.25e-187) {
tmp = x / (x + (y * Math.exp((2.0 * ((b * 0.6666666666666666) / t)))));
} else if (t <= 1.56e-21) {
tmp = x / (x + (y * Math.exp((2.0 * ((c * -0.6666666666666666) / t)))));
} else {
tmp = x / (x + (y * Math.exp((2.0 * (c * (a + 0.8333333333333334))))));
}
return tmp;
}
def code(x, y, z, t, a, b, c): tmp = 0 if t <= 2.25e-187: tmp = x / (x + (y * math.exp((2.0 * ((b * 0.6666666666666666) / t))))) elif t <= 1.56e-21: tmp = x / (x + (y * math.exp((2.0 * ((c * -0.6666666666666666) / t))))) else: tmp = x / (x + (y * math.exp((2.0 * (c * (a + 0.8333333333333334)))))) return tmp
function code(x, y, z, t, a, b, c) tmp = 0.0 if (t <= 2.25e-187) tmp = Float64(x / Float64(x + Float64(y * exp(Float64(2.0 * Float64(Float64(b * 0.6666666666666666) / t)))))); elseif (t <= 1.56e-21) tmp = Float64(x / Float64(x + Float64(y * exp(Float64(2.0 * Float64(Float64(c * -0.6666666666666666) / t)))))); else tmp = Float64(x / Float64(x + Float64(y * exp(Float64(2.0 * Float64(c * Float64(a + 0.8333333333333334))))))); end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c) tmp = 0.0; if (t <= 2.25e-187) tmp = x / (x + (y * exp((2.0 * ((b * 0.6666666666666666) / t))))); elseif (t <= 1.56e-21) tmp = x / (x + (y * exp((2.0 * ((c * -0.6666666666666666) / t))))); else tmp = x / (x + (y * exp((2.0 * (c * (a + 0.8333333333333334)))))); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_] := If[LessEqual[t, 2.25e-187], N[(x / N[(x + N[(y * N[Exp[N[(2.0 * N[(N[(b * 0.6666666666666666), $MachinePrecision] / t), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t, 1.56e-21], N[(x / N[(x + N[(y * N[Exp[N[(2.0 * N[(N[(c * -0.6666666666666666), $MachinePrecision] / t), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x / N[(x + N[(y * N[Exp[N[(2.0 * N[(c * N[(a + 0.8333333333333334), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;t \leq 2.25 \cdot 10^{-187}:\\
\;\;\;\;\frac{x}{x + y \cdot e^{2 \cdot \frac{b \cdot 0.6666666666666666}{t}}}\\
\mathbf{elif}\;t \leq 1.56 \cdot 10^{-21}:\\
\;\;\;\;\frac{x}{x + y \cdot e^{2 \cdot \frac{c \cdot -0.6666666666666666}{t}}}\\
\mathbf{else}:\\
\;\;\;\;\frac{x}{x + y \cdot e^{2 \cdot \left(c \cdot \left(a + 0.8333333333333334\right)\right)}}\\
\end{array}
\end{array}
if t < 2.2499999999999999e-187Initial program 83.8%
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
count-2N/A
exp-lowering-exp.f64N/A
Simplified84.8%
Taylor expanded in b around inf
*-lowering-*.f64N/A
--lowering--.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f64N/A
+-lowering-+.f6468.8%
Simplified68.8%
Taylor expanded in t around 0
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f6469.8%
Simplified69.8%
if 2.2499999999999999e-187 < t < 1.55999999999999999e-21Initial program 95.7%
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
count-2N/A
exp-lowering-exp.f64N/A
Simplified95.7%
Taylor expanded in c around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-lowering-*.f64N/A
--lowering--.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f64N/A
+-lowering-+.f6470.8%
Simplified70.8%
Taylor expanded in t around 0
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f64N/A
metadata-evalN/A
*-lowering-*.f6466.5%
Simplified66.5%
if 1.55999999999999999e-21 < t Initial program 93.4%
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
count-2N/A
exp-lowering-exp.f64N/A
Simplified99.0%
Taylor expanded in c around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-lowering-*.f64N/A
--lowering--.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f64N/A
+-lowering-+.f6468.6%
Simplified68.6%
Taylor expanded in t around inf
*-lowering-*.f64N/A
+-lowering-+.f6469.6%
Simplified69.6%
Final simplification69.1%
(FPCore (x y z t a b c)
:precision binary64
(let* ((t_1 (- (/ 0.6666666666666666 t) (+ a 0.8333333333333334)))
(t_2 (* t_1 t_1))
(t_3 (* t_1 t_2)))
(if (<= b -1.25e+66)
(/
x
(+
x
(*
y
(+
(*
b
(+
(* 2.0 t_1)
(* b (+ (* (* b 1.3333333333333333) t_3) (* 2.0 t_2)))))
1.0))))
(if (<= b -1.16e-46)
1.0
(if (<= b 4.3e-195)
(/
x
(+
(+ x y)
(*
c
(+
(* t_1 (* y -2.0))
(*
c
(+
(* t_2 (* 2.0 y))
(* -1.3333333333333333 (* y (* c t_3)))))))))
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 = t_1 * t_1;
double t_3 = t_1 * t_2;
double tmp;
if (b <= -1.25e+66) {
tmp = x / (x + (y * ((b * ((2.0 * t_1) + (b * (((b * 1.3333333333333333) * t_3) + (2.0 * t_2))))) + 1.0)));
} else if (b <= -1.16e-46) {
tmp = 1.0;
} else if (b <= 4.3e-195) {
tmp = x / ((x + y) + (c * ((t_1 * (y * -2.0)) + (c * ((t_2 * (2.0 * y)) + (-1.3333333333333333 * (y * (c * t_3))))))));
} 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) :: t_1
real(8) :: t_2
real(8) :: t_3
real(8) :: tmp
t_1 = (0.6666666666666666d0 / t) - (a + 0.8333333333333334d0)
t_2 = t_1 * t_1
t_3 = t_1 * t_2
if (b <= (-1.25d+66)) then
tmp = x / (x + (y * ((b * ((2.0d0 * t_1) + (b * (((b * 1.3333333333333333d0) * t_3) + (2.0d0 * t_2))))) + 1.0d0)))
else if (b <= (-1.16d-46)) then
tmp = 1.0d0
else if (b <= 4.3d-195) then
tmp = x / ((x + y) + (c * ((t_1 * (y * (-2.0d0))) + (c * ((t_2 * (2.0d0 * y)) + ((-1.3333333333333333d0) * (y * (c * t_3))))))))
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 t_1 = (0.6666666666666666 / t) - (a + 0.8333333333333334);
double t_2 = t_1 * t_1;
double t_3 = t_1 * t_2;
double tmp;
if (b <= -1.25e+66) {
tmp = x / (x + (y * ((b * ((2.0 * t_1) + (b * (((b * 1.3333333333333333) * t_3) + (2.0 * t_2))))) + 1.0)));
} else if (b <= -1.16e-46) {
tmp = 1.0;
} else if (b <= 4.3e-195) {
tmp = x / ((x + y) + (c * ((t_1 * (y * -2.0)) + (c * ((t_2 * (2.0 * y)) + (-1.3333333333333333 * (y * (c * t_3))))))));
} else {
tmp = 1.0;
}
return tmp;
}
def code(x, y, z, t, a, b, c): t_1 = (0.6666666666666666 / t) - (a + 0.8333333333333334) t_2 = t_1 * t_1 t_3 = t_1 * t_2 tmp = 0 if b <= -1.25e+66: tmp = x / (x + (y * ((b * ((2.0 * t_1) + (b * (((b * 1.3333333333333333) * t_3) + (2.0 * t_2))))) + 1.0))) elif b <= -1.16e-46: tmp = 1.0 elif b <= 4.3e-195: tmp = x / ((x + y) + (c * ((t_1 * (y * -2.0)) + (c * ((t_2 * (2.0 * y)) + (-1.3333333333333333 * (y * (c * t_3)))))))) 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)) t_2 = Float64(t_1 * t_1) t_3 = Float64(t_1 * t_2) tmp = 0.0 if (b <= -1.25e+66) tmp = Float64(x / Float64(x + Float64(y * Float64(Float64(b * Float64(Float64(2.0 * t_1) + Float64(b * Float64(Float64(Float64(b * 1.3333333333333333) * t_3) + Float64(2.0 * t_2))))) + 1.0)))); elseif (b <= -1.16e-46) tmp = 1.0; elseif (b <= 4.3e-195) tmp = Float64(x / Float64(Float64(x + y) + Float64(c * Float64(Float64(t_1 * Float64(y * -2.0)) + Float64(c * Float64(Float64(t_2 * Float64(2.0 * y)) + Float64(-1.3333333333333333 * Float64(y * Float64(c * t_3))))))))); else tmp = 1.0; end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c) t_1 = (0.6666666666666666 / t) - (a + 0.8333333333333334); t_2 = t_1 * t_1; t_3 = t_1 * t_2; tmp = 0.0; if (b <= -1.25e+66) tmp = x / (x + (y * ((b * ((2.0 * t_1) + (b * (((b * 1.3333333333333333) * t_3) + (2.0 * t_2))))) + 1.0))); elseif (b <= -1.16e-46) tmp = 1.0; elseif (b <= 4.3e-195) tmp = x / ((x + y) + (c * ((t_1 * (y * -2.0)) + (c * ((t_2 * (2.0 * y)) + (-1.3333333333333333 * (y * (c * t_3)))))))); else tmp = 1.0; end tmp_2 = 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[(t$95$1 * t$95$1), $MachinePrecision]}, Block[{t$95$3 = N[(t$95$1 * t$95$2), $MachinePrecision]}, If[LessEqual[b, -1.25e+66], N[(x / N[(x + N[(y * N[(N[(b * N[(N[(2.0 * t$95$1), $MachinePrecision] + N[(b * N[(N[(N[(b * 1.3333333333333333), $MachinePrecision] * t$95$3), $MachinePrecision] + N[(2.0 * t$95$2), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[b, -1.16e-46], 1.0, If[LessEqual[b, 4.3e-195], N[(x / N[(N[(x + y), $MachinePrecision] + N[(c * N[(N[(t$95$1 * N[(y * -2.0), $MachinePrecision]), $MachinePrecision] + N[(c * N[(N[(t$95$2 * N[(2.0 * y), $MachinePrecision]), $MachinePrecision] + N[(-1.3333333333333333 * N[(y * N[(c * t$95$3), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 1.0]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{0.6666666666666666}{t} - \left(a + 0.8333333333333334\right)\\
t_2 := t\_1 \cdot t\_1\\
t_3 := t\_1 \cdot t\_2\\
\mathbf{if}\;b \leq -1.25 \cdot 10^{+66}:\\
\;\;\;\;\frac{x}{x + y \cdot \left(b \cdot \left(2 \cdot t\_1 + b \cdot \left(\left(b \cdot 1.3333333333333333\right) \cdot t\_3 + 2 \cdot t\_2\right)\right) + 1\right)}\\
\mathbf{elif}\;b \leq -1.16 \cdot 10^{-46}:\\
\;\;\;\;1\\
\mathbf{elif}\;b \leq 4.3 \cdot 10^{-195}:\\
\;\;\;\;\frac{x}{\left(x + y\right) + c \cdot \left(t\_1 \cdot \left(y \cdot -2\right) + c \cdot \left(t\_2 \cdot \left(2 \cdot y\right) + -1.3333333333333333 \cdot \left(y \cdot \left(c \cdot t\_3\right)\right)\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;1\\
\end{array}
\end{array}
if b < -1.24999999999999998e66Initial program 83.1%
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
count-2N/A
exp-lowering-exp.f64N/A
Simplified84.9%
Taylor expanded in b around inf
*-lowering-*.f64N/A
--lowering--.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f64N/A
+-lowering-+.f6485.4%
Simplified85.4%
Taylor expanded in b around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
Simplified68.6%
if -1.24999999999999998e66 < b < -1.16e-46 or 4.3000000000000004e-195 < b Initial program 89.5%
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
count-2N/A
exp-lowering-exp.f64N/A
Simplified92.1%
Taylor expanded in x around inf
Simplified62.6%
if -1.16e-46 < b < 4.3000000000000004e-195Initial program 94.4%
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
count-2N/A
exp-lowering-exp.f64N/A
Simplified97.7%
Taylor expanded in c around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-lowering-*.f64N/A
--lowering--.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f64N/A
+-lowering-+.f6479.0%
Simplified79.0%
Taylor expanded in c around 0
associate-+r+N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
Simplified53.8%
*-commutativeN/A
associate-*r*N/A
*-lowering-*.f64N/A
Applied egg-rr58.4%
Final simplification62.4%
(FPCore (x y z t a b c)
:precision binary64
(let* ((t_1 (- (/ 0.6666666666666666 t) (+ a 0.8333333333333334))))
(if (<= b 1.9e-195)
(/ x (+ x (* y (+ (* b (* 2.0 (+ t_1 (* b (* 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 (b <= 1.9e-195) {
tmp = x / (x + (y * ((b * (2.0 * (t_1 + (b * (t_1 * t_1))))) + 1.0)));
} 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) :: t_1
real(8) :: tmp
t_1 = (0.6666666666666666d0 / t) - (a + 0.8333333333333334d0)
if (b <= 1.9d-195) then
tmp = x / (x + (y * ((b * (2.0d0 * (t_1 + (b * (t_1 * t_1))))) + 1.0d0)))
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 t_1 = (0.6666666666666666 / t) - (a + 0.8333333333333334);
double tmp;
if (b <= 1.9e-195) {
tmp = x / (x + (y * ((b * (2.0 * (t_1 + (b * (t_1 * t_1))))) + 1.0)));
} else {
tmp = 1.0;
}
return tmp;
}
def code(x, y, z, t, a, b, c): t_1 = (0.6666666666666666 / t) - (a + 0.8333333333333334) tmp = 0 if b <= 1.9e-195: tmp = x / (x + (y * ((b * (2.0 * (t_1 + (b * (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 (b <= 1.9e-195) tmp = Float64(x / Float64(x + Float64(y * Float64(Float64(b * Float64(2.0 * Float64(t_1 + Float64(b * Float64(t_1 * t_1))))) + 1.0)))); else tmp = 1.0; end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c) t_1 = (0.6666666666666666 / t) - (a + 0.8333333333333334); tmp = 0.0; if (b <= 1.9e-195) tmp = x / (x + (y * ((b * (2.0 * (t_1 + (b * (t_1 * t_1))))) + 1.0))); else tmp = 1.0; end tmp_2 = 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[b, 1.9e-195], N[(x / N[(x + N[(y * N[(N[(b * N[(2.0 * N[(t$95$1 + N[(b * N[(t$95$1 * t$95$1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $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}\;b \leq 1.9 \cdot 10^{-195}:\\
\;\;\;\;\frac{x}{x + y \cdot \left(b \cdot \left(2 \cdot \left(t\_1 + b \cdot \left(t\_1 \cdot t\_1\right)\right)\right) + 1\right)}\\
\mathbf{else}:\\
\;\;\;\;1\\
\end{array}
\end{array}
if b < 1.90000000000000006e-195Initial program 90.6%
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
count-2N/A
exp-lowering-exp.f64N/A
Simplified93.1%
Taylor expanded in b around inf
*-lowering-*.f64N/A
--lowering--.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f64N/A
+-lowering-+.f6463.6%
Simplified63.6%
Taylor expanded in b around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
distribute-lft-outN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
Simplified57.5%
if 1.90000000000000006e-195 < b Initial program 88.7%
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
count-2N/A
exp-lowering-exp.f64N/A
Simplified91.8%
Taylor expanded in x around inf
Simplified61.1%
Final simplification58.8%
(FPCore (x y z t a b c)
:precision binary64
(if (<= b -9.2e+62)
(/
x
(+
x
(*
y
(+
(*
b
(*
2.0
(-
(* b (* (+ a 0.8333333333333334) (+ a 0.8333333333333334)))
(+ a 0.8333333333333334))))
1.0))))
(if (<= b -6.2e-45)
1.0
(if (<= b 6.8e-197)
(/
x
(+
(+ x y)
(*
-2.0
(* (- (/ 0.6666666666666666 t) (+ a 0.8333333333333334)) (* c y)))))
1.0))))
double code(double x, double y, double z, double t, double a, double b, double c) {
double tmp;
if (b <= -9.2e+62) {
tmp = x / (x + (y * ((b * (2.0 * ((b * ((a + 0.8333333333333334) * (a + 0.8333333333333334))) - (a + 0.8333333333333334)))) + 1.0)));
} else if (b <= -6.2e-45) {
tmp = 1.0;
} else if (b <= 6.8e-197) {
tmp = x / ((x + y) + (-2.0 * (((0.6666666666666666 / t) - (a + 0.8333333333333334)) * (c * 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 (b <= (-9.2d+62)) then
tmp = x / (x + (y * ((b * (2.0d0 * ((b * ((a + 0.8333333333333334d0) * (a + 0.8333333333333334d0))) - (a + 0.8333333333333334d0)))) + 1.0d0)))
else if (b <= (-6.2d-45)) then
tmp = 1.0d0
else if (b <= 6.8d-197) then
tmp = x / ((x + y) + ((-2.0d0) * (((0.6666666666666666d0 / t) - (a + 0.8333333333333334d0)) * (c * 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 (b <= -9.2e+62) {
tmp = x / (x + (y * ((b * (2.0 * ((b * ((a + 0.8333333333333334) * (a + 0.8333333333333334))) - (a + 0.8333333333333334)))) + 1.0)));
} else if (b <= -6.2e-45) {
tmp = 1.0;
} else if (b <= 6.8e-197) {
tmp = x / ((x + y) + (-2.0 * (((0.6666666666666666 / t) - (a + 0.8333333333333334)) * (c * y))));
} else {
tmp = 1.0;
}
return tmp;
}
def code(x, y, z, t, a, b, c): tmp = 0 if b <= -9.2e+62: tmp = x / (x + (y * ((b * (2.0 * ((b * ((a + 0.8333333333333334) * (a + 0.8333333333333334))) - (a + 0.8333333333333334)))) + 1.0))) elif b <= -6.2e-45: tmp = 1.0 elif b <= 6.8e-197: tmp = x / ((x + y) + (-2.0 * (((0.6666666666666666 / t) - (a + 0.8333333333333334)) * (c * y)))) else: tmp = 1.0 return tmp
function code(x, y, z, t, a, b, c) tmp = 0.0 if (b <= -9.2e+62) tmp = Float64(x / Float64(x + Float64(y * Float64(Float64(b * Float64(2.0 * Float64(Float64(b * Float64(Float64(a + 0.8333333333333334) * Float64(a + 0.8333333333333334))) - Float64(a + 0.8333333333333334)))) + 1.0)))); elseif (b <= -6.2e-45) tmp = 1.0; elseif (b <= 6.8e-197) tmp = Float64(x / Float64(Float64(x + y) + Float64(-2.0 * Float64(Float64(Float64(0.6666666666666666 / t) - Float64(a + 0.8333333333333334)) * Float64(c * y))))); else tmp = 1.0; end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c) tmp = 0.0; if (b <= -9.2e+62) tmp = x / (x + (y * ((b * (2.0 * ((b * ((a + 0.8333333333333334) * (a + 0.8333333333333334))) - (a + 0.8333333333333334)))) + 1.0))); elseif (b <= -6.2e-45) tmp = 1.0; elseif (b <= 6.8e-197) tmp = x / ((x + y) + (-2.0 * (((0.6666666666666666 / t) - (a + 0.8333333333333334)) * (c * y)))); else tmp = 1.0; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_] := If[LessEqual[b, -9.2e+62], N[(x / N[(x + N[(y * N[(N[(b * N[(2.0 * N[(N[(b * N[(N[(a + 0.8333333333333334), $MachinePrecision] * N[(a + 0.8333333333333334), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - N[(a + 0.8333333333333334), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[b, -6.2e-45], 1.0, If[LessEqual[b, 6.8e-197], N[(x / N[(N[(x + y), $MachinePrecision] + N[(-2.0 * N[(N[(N[(0.6666666666666666 / t), $MachinePrecision] - N[(a + 0.8333333333333334), $MachinePrecision]), $MachinePrecision] * N[(c * y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 1.0]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;b \leq -9.2 \cdot 10^{+62}:\\
\;\;\;\;\frac{x}{x + y \cdot \left(b \cdot \left(2 \cdot \left(b \cdot \left(\left(a + 0.8333333333333334\right) \cdot \left(a + 0.8333333333333334\right)\right) - \left(a + 0.8333333333333334\right)\right)\right) + 1\right)}\\
\mathbf{elif}\;b \leq -6.2 \cdot 10^{-45}:\\
\;\;\;\;1\\
\mathbf{elif}\;b \leq 6.8 \cdot 10^{-197}:\\
\;\;\;\;\frac{x}{\left(x + y\right) + -2 \cdot \left(\left(\frac{0.6666666666666666}{t} - \left(a + 0.8333333333333334\right)\right) \cdot \left(c \cdot y\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;1\\
\end{array}
\end{array}
if b < -9.19999999999999936e62Initial program 83.1%
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
count-2N/A
exp-lowering-exp.f64N/A
Simplified84.9%
Taylor expanded in b around inf
*-lowering-*.f64N/A
--lowering--.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f64N/A
+-lowering-+.f6485.4%
Simplified85.4%
Taylor expanded in b around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
distribute-lft-outN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
Simplified65.3%
Taylor expanded in t around inf
--lowering--.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f6463.5%
Simplified63.5%
if -9.19999999999999936e62 < b < -6.2000000000000002e-45 or 6.7999999999999996e-197 < b Initial program 89.4%
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
count-2N/A
exp-lowering-exp.f64N/A
Simplified92.0%
Taylor expanded in x around inf
Simplified63.1%
if -6.2000000000000002e-45 < b < 6.7999999999999996e-197Initial program 94.5%
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
count-2N/A
exp-lowering-exp.f64N/A
Simplified97.8%
Taylor expanded in c around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-lowering-*.f64N/A
--lowering--.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f64N/A
+-lowering-+.f6478.2%
Simplified78.2%
Taylor expanded in c around 0
associate-+r+N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
--lowering--.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f64N/A
+-lowering-+.f6454.6%
Simplified54.6%
Final simplification60.2%
(FPCore (x y z t a b c)
:precision binary64
(if (<= b -1.2e+84)
(/ x (* -2.6666666666666665 (* a (/ (* y (* b b)) t))))
(if (<= b -1.25e-45)
1.0
(if (<= b 1.1e-195)
(/
x
(+
(+ x y)
(*
-2.0
(* (- (/ 0.6666666666666666 t) (+ a 0.8333333333333334)) (* c y)))))
1.0))))
double code(double x, double y, double z, double t, double a, double b, double c) {
double tmp;
if (b <= -1.2e+84) {
tmp = x / (-2.6666666666666665 * (a * ((y * (b * b)) / t)));
} else if (b <= -1.25e-45) {
tmp = 1.0;
} else if (b <= 1.1e-195) {
tmp = x / ((x + y) + (-2.0 * (((0.6666666666666666 / t) - (a + 0.8333333333333334)) * (c * 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 (b <= (-1.2d+84)) then
tmp = x / ((-2.6666666666666665d0) * (a * ((y * (b * b)) / t)))
else if (b <= (-1.25d-45)) then
tmp = 1.0d0
else if (b <= 1.1d-195) then
tmp = x / ((x + y) + ((-2.0d0) * (((0.6666666666666666d0 / t) - (a + 0.8333333333333334d0)) * (c * 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 (b <= -1.2e+84) {
tmp = x / (-2.6666666666666665 * (a * ((y * (b * b)) / t)));
} else if (b <= -1.25e-45) {
tmp = 1.0;
} else if (b <= 1.1e-195) {
tmp = x / ((x + y) + (-2.0 * (((0.6666666666666666 / t) - (a + 0.8333333333333334)) * (c * y))));
} else {
tmp = 1.0;
}
return tmp;
}
def code(x, y, z, t, a, b, c): tmp = 0 if b <= -1.2e+84: tmp = x / (-2.6666666666666665 * (a * ((y * (b * b)) / t))) elif b <= -1.25e-45: tmp = 1.0 elif b <= 1.1e-195: tmp = x / ((x + y) + (-2.0 * (((0.6666666666666666 / t) - (a + 0.8333333333333334)) * (c * y)))) else: tmp = 1.0 return tmp
function code(x, y, z, t, a, b, c) tmp = 0.0 if (b <= -1.2e+84) tmp = Float64(x / Float64(-2.6666666666666665 * Float64(a * Float64(Float64(y * Float64(b * b)) / t)))); elseif (b <= -1.25e-45) tmp = 1.0; elseif (b <= 1.1e-195) tmp = Float64(x / Float64(Float64(x + y) + Float64(-2.0 * Float64(Float64(Float64(0.6666666666666666 / t) - Float64(a + 0.8333333333333334)) * Float64(c * y))))); else tmp = 1.0; end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c) tmp = 0.0; if (b <= -1.2e+84) tmp = x / (-2.6666666666666665 * (a * ((y * (b * b)) / t))); elseif (b <= -1.25e-45) tmp = 1.0; elseif (b <= 1.1e-195) tmp = x / ((x + y) + (-2.0 * (((0.6666666666666666 / t) - (a + 0.8333333333333334)) * (c * y)))); else tmp = 1.0; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_] := If[LessEqual[b, -1.2e+84], N[(x / N[(-2.6666666666666665 * N[(a * N[(N[(y * N[(b * b), $MachinePrecision]), $MachinePrecision] / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[b, -1.25e-45], 1.0, If[LessEqual[b, 1.1e-195], N[(x / N[(N[(x + y), $MachinePrecision] + N[(-2.0 * N[(N[(N[(0.6666666666666666 / t), $MachinePrecision] - N[(a + 0.8333333333333334), $MachinePrecision]), $MachinePrecision] * N[(c * y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 1.0]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;b \leq -1.2 \cdot 10^{+84}:\\
\;\;\;\;\frac{x}{-2.6666666666666665 \cdot \left(a \cdot \frac{y \cdot \left(b \cdot b\right)}{t}\right)}\\
\mathbf{elif}\;b \leq -1.25 \cdot 10^{-45}:\\
\;\;\;\;1\\
\mathbf{elif}\;b \leq 1.1 \cdot 10^{-195}:\\
\;\;\;\;\frac{x}{\left(x + y\right) + -2 \cdot \left(\left(\frac{0.6666666666666666}{t} - \left(a + 0.8333333333333334\right)\right) \cdot \left(c \cdot y\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;1\\
\end{array}
\end{array}
if b < -1.2e84Initial program 85.1%
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
count-2N/A
exp-lowering-exp.f64N/A
Simplified85.1%
Taylor expanded in b around inf
*-lowering-*.f64N/A
--lowering--.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f64N/A
+-lowering-+.f6489.7%
Simplified89.7%
Taylor expanded in b around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
distribute-lft-outN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
Simplified67.0%
Taylor expanded in t around 0
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
unpow2N/A
*-lowering-*.f6440.9%
Simplified40.9%
Taylor expanded in a around inf
*-lowering-*.f64N/A
associate-/l*N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6467.0%
Simplified67.0%
if -1.2e84 < b < -1.24999999999999994e-45 or 1.10000000000000003e-195 < b Initial program 88.3%
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
count-2N/A
exp-lowering-exp.f64N/A
Simplified91.6%
Taylor expanded in x around inf
Simplified61.7%
if -1.24999999999999994e-45 < b < 1.10000000000000003e-195Initial program 94.5%
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
count-2N/A
exp-lowering-exp.f64N/A
Simplified97.8%
Taylor expanded in c around inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-lowering-*.f64N/A
--lowering--.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f64N/A
+-lowering-+.f6478.2%
Simplified78.2%
Taylor expanded in c around 0
associate-+r+N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
--lowering--.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f64N/A
+-lowering-+.f6454.6%
Simplified54.6%
Final simplification60.2%
(FPCore (x y z t a b c) :precision binary64 (if (<= t 3.1e-282) 1.0 (if (<= t 1.85e-188) (* 1.125 (/ (* x (* t t)) (* y (* b b)))) 1.0)))
double code(double x, double y, double z, double t, double a, double b, double c) {
double tmp;
if (t <= 3.1e-282) {
tmp = 1.0;
} else if (t <= 1.85e-188) {
tmp = 1.125 * ((x * (t * t)) / (y * (b * b)));
} 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 (t <= 3.1d-282) then
tmp = 1.0d0
else if (t <= 1.85d-188) then
tmp = 1.125d0 * ((x * (t * t)) / (y * (b * b)))
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 (t <= 3.1e-282) {
tmp = 1.0;
} else if (t <= 1.85e-188) {
tmp = 1.125 * ((x * (t * t)) / (y * (b * b)));
} else {
tmp = 1.0;
}
return tmp;
}
def code(x, y, z, t, a, b, c): tmp = 0 if t <= 3.1e-282: tmp = 1.0 elif t <= 1.85e-188: tmp = 1.125 * ((x * (t * t)) / (y * (b * b))) else: tmp = 1.0 return tmp
function code(x, y, z, t, a, b, c) tmp = 0.0 if (t <= 3.1e-282) tmp = 1.0; elseif (t <= 1.85e-188) tmp = Float64(1.125 * Float64(Float64(x * Float64(t * t)) / Float64(y * Float64(b * b)))); else tmp = 1.0; end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c) tmp = 0.0; if (t <= 3.1e-282) tmp = 1.0; elseif (t <= 1.85e-188) tmp = 1.125 * ((x * (t * t)) / (y * (b * b))); else tmp = 1.0; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_] := If[LessEqual[t, 3.1e-282], 1.0, If[LessEqual[t, 1.85e-188], N[(1.125 * N[(N[(x * N[(t * t), $MachinePrecision]), $MachinePrecision] / N[(y * N[(b * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 1.0]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;t \leq 3.1 \cdot 10^{-282}:\\
\;\;\;\;1\\
\mathbf{elif}\;t \leq 1.85 \cdot 10^{-188}:\\
\;\;\;\;1.125 \cdot \frac{x \cdot \left(t \cdot t\right)}{y \cdot \left(b \cdot b\right)}\\
\mathbf{else}:\\
\;\;\;\;1\\
\end{array}
\end{array}
if t < 3.10000000000000013e-282 or 1.84999999999999986e-188 < t Initial program 91.7%
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
count-2N/A
exp-lowering-exp.f64N/A
Simplified94.7%
Taylor expanded in x around inf
Simplified52.4%
if 3.10000000000000013e-282 < t < 1.84999999999999986e-188Initial program 75.0%
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
count-2N/A
exp-lowering-exp.f64N/A
Simplified75.0%
Taylor expanded in b around inf
*-lowering-*.f64N/A
--lowering--.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f64N/A
+-lowering-+.f6469.1%
Simplified69.1%
Taylor expanded in b around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
distribute-lft-outN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
Simplified82.7%
Taylor expanded in t around 0
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6471.7%
Simplified71.7%
Final simplification54.5%
(FPCore (x y z t a b c) :precision binary64 (if (<= b -2.15e+85) (/ x (* -2.6666666666666665 (* a (/ (* y (* b b)) t)))) 1.0))
double code(double x, double y, double z, double t, double a, double b, double c) {
double tmp;
if (b <= -2.15e+85) {
tmp = x / (-2.6666666666666665 * (a * ((y * (b * b)) / t)));
} 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 (b <= (-2.15d+85)) then
tmp = x / ((-2.6666666666666665d0) * (a * ((y * (b * b)) / t)))
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 (b <= -2.15e+85) {
tmp = x / (-2.6666666666666665 * (a * ((y * (b * b)) / t)));
} else {
tmp = 1.0;
}
return tmp;
}
def code(x, y, z, t, a, b, c): tmp = 0 if b <= -2.15e+85: tmp = x / (-2.6666666666666665 * (a * ((y * (b * b)) / t))) else: tmp = 1.0 return tmp
function code(x, y, z, t, a, b, c) tmp = 0.0 if (b <= -2.15e+85) tmp = Float64(x / Float64(-2.6666666666666665 * Float64(a * Float64(Float64(y * Float64(b * b)) / t)))); else tmp = 1.0; end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c) tmp = 0.0; if (b <= -2.15e+85) tmp = x / (-2.6666666666666665 * (a * ((y * (b * b)) / t))); else tmp = 1.0; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_] := If[LessEqual[b, -2.15e+85], N[(x / N[(-2.6666666666666665 * N[(a * N[(N[(y * N[(b * b), $MachinePrecision]), $MachinePrecision] / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 1.0]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;b \leq -2.15 \cdot 10^{+85}:\\
\;\;\;\;\frac{x}{-2.6666666666666665 \cdot \left(a \cdot \frac{y \cdot \left(b \cdot b\right)}{t}\right)}\\
\mathbf{else}:\\
\;\;\;\;1\\
\end{array}
\end{array}
if b < -2.15e85Initial program 85.1%
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
count-2N/A
exp-lowering-exp.f64N/A
Simplified85.1%
Taylor expanded in b around inf
*-lowering-*.f64N/A
--lowering--.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f64N/A
+-lowering-+.f6489.7%
Simplified89.7%
Taylor expanded in b around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
distribute-lft-outN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
Simplified67.0%
Taylor expanded in t around 0
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
unpow2N/A
*-lowering-*.f6440.9%
Simplified40.9%
Taylor expanded in a around inf
*-lowering-*.f64N/A
associate-/l*N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6467.0%
Simplified67.0%
if -2.15e85 < b Initial program 91.0%
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
count-2N/A
exp-lowering-exp.f64N/A
Simplified94.3%
Taylor expanded in x around inf
Simplified52.7%
Final simplification55.3%
(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 89.9%
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
count-2N/A
exp-lowering-exp.f64N/A
Simplified92.6%
Taylor expanded in x around inf
Simplified48.9%
(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 2024161
(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)))))))))))