
(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 25 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
(+
(/ (* z (sqrt (+ t a))) t)
(* (- b c) (- (/ 2.0 (* t 3.0)) (+ a 0.8333333333333334))))))
(if (<= t_1 INFINITY)
(/ x (+ x (* y (exp (* 2.0 t_1)))))
(/
x
(+
x
(*
y
(exp
(*
2.0
(/ (- (* z (sqrt a)) (* (- b c) -0.6666666666666666)) t)))))))))
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 + 0.8333333333333334)));
double tmp;
if (t_1 <= ((double) INFINITY)) {
tmp = x / (x + (y * exp((2.0 * t_1))));
} else {
tmp = x / (x + (y * exp((2.0 * (((z * sqrt(a)) - ((b - c) * -0.6666666666666666)) / t)))));
}
return tmp;
}
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))) / t) + ((b - c) * ((2.0 / (t * 3.0)) - (a + 0.8333333333333334)));
double tmp;
if (t_1 <= Double.POSITIVE_INFINITY) {
tmp = x / (x + (y * Math.exp((2.0 * t_1))));
} else {
tmp = x / (x + (y * Math.exp((2.0 * (((z * Math.sqrt(a)) - ((b - c) * -0.6666666666666666)) / t)))));
}
return tmp;
}
def code(x, y, z, t, a, b, c): t_1 = ((z * math.sqrt((t + a))) / t) + ((b - c) * ((2.0 / (t * 3.0)) - (a + 0.8333333333333334))) tmp = 0 if t_1 <= math.inf: tmp = x / (x + (y * math.exp((2.0 * t_1)))) else: tmp = x / (x + (y * math.exp((2.0 * (((z * math.sqrt(a)) - ((b - c) * -0.6666666666666666)) / t))))) 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 + 0.8333333333333334)))) tmp = 0.0 if (t_1 <= Inf) tmp = Float64(x / Float64(x + Float64(y * exp(Float64(2.0 * t_1))))); else tmp = Float64(x / Float64(x + Float64(y * exp(Float64(2.0 * Float64(Float64(Float64(z * sqrt(a)) - Float64(Float64(b - c) * -0.6666666666666666)) / t)))))); end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c) t_1 = ((z * sqrt((t + a))) / t) + ((b - c) * ((2.0 / (t * 3.0)) - (a + 0.8333333333333334))); tmp = 0.0; if (t_1 <= Inf) tmp = x / (x + (y * exp((2.0 * t_1)))); else tmp = x / (x + (y * exp((2.0 * (((z * sqrt(a)) - ((b - c) * -0.6666666666666666)) / t))))); end tmp_2 = 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 + 0.8333333333333334), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$1, Infinity], N[(x / N[(x + N[(y * N[Exp[N[(2.0 * t$95$1), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x / N[(x + N[(y * N[Exp[N[(2.0 * N[(N[(N[(z * N[Sqrt[a], $MachinePrecision]), $MachinePrecision] - N[(N[(b - c), $MachinePrecision] * -0.6666666666666666), $MachinePrecision]), $MachinePrecision] / t), $MachinePrecision]), $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 + 0.8333333333333334\right)\right)\\
\mathbf{if}\;t\_1 \leq \infty:\\
\;\;\;\;\frac{x}{x + y \cdot e^{2 \cdot t\_1}}\\
\mathbf{else}:\\
\;\;\;\;\frac{x}{x + y \cdot e^{2 \cdot \frac{z \cdot \sqrt{a} - \left(b - c\right) \cdot -0.6666666666666666}{t}}}\\
\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 99.6%
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%
Taylor expanded in t around 0 77.4%
Final simplification98.5%
(FPCore (x y z t a b c)
:precision binary64
(if (<= t 9.6e-272)
(/
x
(+
x
(*
y
(exp (* 2.0 (/ (- (* z (sqrt a)) (* (- b c) -0.6666666666666666)) t))))))
(if (<= t 6.2e-12)
(/
x
(+
x
(*
y
(exp
(*
2.0
(* b (- (/ 0.6666666666666666 t) (+ a 0.8333333333333334))))))))
(/
x
(+
x
(*
y
(exp
(*
2.0
(-
(* z (sqrt (/ 1.0 t)))
(* (- b c) (+ a 0.8333333333333334)))))))))))
double code(double x, double y, double z, double t, double a, double b, double c) {
double tmp;
if (t <= 9.6e-272) {
tmp = x / (x + (y * exp((2.0 * (((z * sqrt(a)) - ((b - c) * -0.6666666666666666)) / t)))));
} else if (t <= 6.2e-12) {
tmp = x / (x + (y * exp((2.0 * (b * ((0.6666666666666666 / t) - (a + 0.8333333333333334)))))));
} else {
tmp = x / (x + (y * exp((2.0 * ((z * sqrt((1.0 / t))) - ((b - 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 <= 9.6d-272) then
tmp = x / (x + (y * exp((2.0d0 * (((z * sqrt(a)) - ((b - c) * (-0.6666666666666666d0))) / t)))))
else if (t <= 6.2d-12) then
tmp = x / (x + (y * exp((2.0d0 * (b * ((0.6666666666666666d0 / t) - (a + 0.8333333333333334d0)))))))
else
tmp = x / (x + (y * exp((2.0d0 * ((z * sqrt((1.0d0 / t))) - ((b - 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 <= 9.6e-272) {
tmp = x / (x + (y * Math.exp((2.0 * (((z * Math.sqrt(a)) - ((b - c) * -0.6666666666666666)) / t)))));
} else if (t <= 6.2e-12) {
tmp = x / (x + (y * Math.exp((2.0 * (b * ((0.6666666666666666 / t) - (a + 0.8333333333333334)))))));
} else {
tmp = x / (x + (y * Math.exp((2.0 * ((z * Math.sqrt((1.0 / t))) - ((b - c) * (a + 0.8333333333333334)))))));
}
return tmp;
}
def code(x, y, z, t, a, b, c): tmp = 0 if t <= 9.6e-272: tmp = x / (x + (y * math.exp((2.0 * (((z * math.sqrt(a)) - ((b - c) * -0.6666666666666666)) / t))))) elif t <= 6.2e-12: tmp = x / (x + (y * math.exp((2.0 * (b * ((0.6666666666666666 / t) - (a + 0.8333333333333334))))))) else: tmp = x / (x + (y * math.exp((2.0 * ((z * math.sqrt((1.0 / t))) - ((b - c) * (a + 0.8333333333333334))))))) return tmp
function code(x, y, z, t, a, b, c) tmp = 0.0 if (t <= 9.6e-272) tmp = Float64(x / Float64(x + Float64(y * exp(Float64(2.0 * Float64(Float64(Float64(z * sqrt(a)) - Float64(Float64(b - c) * -0.6666666666666666)) / t)))))); elseif (t <= 6.2e-12) tmp = Float64(x / Float64(x + Float64(y * exp(Float64(2.0 * Float64(b * Float64(Float64(0.6666666666666666 / t) - Float64(a + 0.8333333333333334)))))))); else tmp = Float64(x / Float64(x + Float64(y * exp(Float64(2.0 * Float64(Float64(z * sqrt(Float64(1.0 / t))) - Float64(Float64(b - 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 <= 9.6e-272) tmp = x / (x + (y * exp((2.0 * (((z * sqrt(a)) - ((b - c) * -0.6666666666666666)) / t))))); elseif (t <= 6.2e-12) tmp = x / (x + (y * exp((2.0 * (b * ((0.6666666666666666 / t) - (a + 0.8333333333333334))))))); else tmp = x / (x + (y * exp((2.0 * ((z * sqrt((1.0 / t))) - ((b - c) * (a + 0.8333333333333334))))))); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_] := If[LessEqual[t, 9.6e-272], N[(x / N[(x + N[(y * N[Exp[N[(2.0 * N[(N[(N[(z * N[Sqrt[a], $MachinePrecision]), $MachinePrecision] - N[(N[(b - c), $MachinePrecision] * -0.6666666666666666), $MachinePrecision]), $MachinePrecision] / t), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t, 6.2e-12], 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], N[(x / N[(x + N[(y * N[Exp[N[(2.0 * N[(N[(z * N[Sqrt[N[(1.0 / t), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] - N[(N[(b - c), $MachinePrecision] * N[(a + 0.8333333333333334), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;t \leq 9.6 \cdot 10^{-272}:\\
\;\;\;\;\frac{x}{x + y \cdot e^{2 \cdot \frac{z \cdot \sqrt{a} - \left(b - c\right) \cdot -0.6666666666666666}{t}}}\\
\mathbf{elif}\;t \leq 6.2 \cdot 10^{-12}:\\
\;\;\;\;\frac{x}{x + y \cdot e^{2 \cdot \left(b \cdot \left(\frac{0.6666666666666666}{t} - \left(a + 0.8333333333333334\right)\right)\right)}}\\
\mathbf{else}:\\
\;\;\;\;\frac{x}{x + y \cdot e^{2 \cdot \left(z \cdot \sqrt{\frac{1}{t}} - \left(b - c\right) \cdot \left(a + 0.8333333333333334\right)\right)}}\\
\end{array}
\end{array}
if t < 9.59999999999999959e-272Initial program 88.7%
Taylor expanded in t around 0 95.9%
if 9.59999999999999959e-272 < t < 6.2000000000000002e-12Initial program 95.5%
Taylor expanded in b around inf 74.1%
associate-*r/74.1%
metadata-eval74.1%
+-commutative74.1%
Simplified74.1%
if 6.2000000000000002e-12 < t Initial program 97.4%
Taylor expanded in t around inf 98.3%
Final simplification91.3%
(FPCore (x y z t a b c)
:precision binary64
(let* ((t_1 (/ x (+ x (* y (exp (* 2.0 (* a (- b)))))))))
(if (<= a -54000000.0)
t_1
(if (<= a -5.2e-166)
(/
x
(+
x
(*
b
(/
(+
(* y 1.3333333333333333)
(* t (+ (* -2.0 (* y (+ a 0.8333333333333334))) (/ y b))))
t))))
(if (<= a 7.2e-268)
(/ x (+ x (* y (exp (* -1.3333333333333333 (/ c t))))))
(if (<= a 1.05e-41)
(/ x (+ x (* y (exp (* 1.3333333333333333 (/ b t))))))
(if (or (<= a 1.4e+179) (not (<= a 4.6e+280)))
t_1
(/ x (+ x (* y (exp (* c 1.6666666666666667))))))))))))
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 * (a * -b)))));
double tmp;
if (a <= -54000000.0) {
tmp = t_1;
} else if (a <= -5.2e-166) {
tmp = x / (x + (b * (((y * 1.3333333333333333) + (t * ((-2.0 * (y * (a + 0.8333333333333334))) + (y / b)))) / t)));
} else if (a <= 7.2e-268) {
tmp = x / (x + (y * exp((-1.3333333333333333 * (c / t)))));
} else if (a <= 1.05e-41) {
tmp = x / (x + (y * exp((1.3333333333333333 * (b / t)))));
} else if ((a <= 1.4e+179) || !(a <= 4.6e+280)) {
tmp = t_1;
} else {
tmp = x / (x + (y * exp((c * 1.6666666666666667))));
}
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 * (a * -b)))))
if (a <= (-54000000.0d0)) then
tmp = t_1
else if (a <= (-5.2d-166)) then
tmp = x / (x + (b * (((y * 1.3333333333333333d0) + (t * (((-2.0d0) * (y * (a + 0.8333333333333334d0))) + (y / b)))) / t)))
else if (a <= 7.2d-268) then
tmp = x / (x + (y * exp(((-1.3333333333333333d0) * (c / t)))))
else if (a <= 1.05d-41) then
tmp = x / (x + (y * exp((1.3333333333333333d0 * (b / t)))))
else if ((a <= 1.4d+179) .or. (.not. (a <= 4.6d+280))) then
tmp = t_1
else
tmp = x / (x + (y * exp((c * 1.6666666666666667d0))))
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 * (a * -b)))));
double tmp;
if (a <= -54000000.0) {
tmp = t_1;
} else if (a <= -5.2e-166) {
tmp = x / (x + (b * (((y * 1.3333333333333333) + (t * ((-2.0 * (y * (a + 0.8333333333333334))) + (y / b)))) / t)));
} else if (a <= 7.2e-268) {
tmp = x / (x + (y * Math.exp((-1.3333333333333333 * (c / t)))));
} else if (a <= 1.05e-41) {
tmp = x / (x + (y * Math.exp((1.3333333333333333 * (b / t)))));
} else if ((a <= 1.4e+179) || !(a <= 4.6e+280)) {
tmp = t_1;
} else {
tmp = x / (x + (y * Math.exp((c * 1.6666666666666667))));
}
return tmp;
}
def code(x, y, z, t, a, b, c): t_1 = x / (x + (y * math.exp((2.0 * (a * -b))))) tmp = 0 if a <= -54000000.0: tmp = t_1 elif a <= -5.2e-166: tmp = x / (x + (b * (((y * 1.3333333333333333) + (t * ((-2.0 * (y * (a + 0.8333333333333334))) + (y / b)))) / t))) elif a <= 7.2e-268: tmp = x / (x + (y * math.exp((-1.3333333333333333 * (c / t))))) elif a <= 1.05e-41: tmp = x / (x + (y * math.exp((1.3333333333333333 * (b / t))))) elif (a <= 1.4e+179) or not (a <= 4.6e+280): tmp = t_1 else: tmp = x / (x + (y * math.exp((c * 1.6666666666666667)))) return tmp
function code(x, y, z, t, a, b, c) t_1 = Float64(x / Float64(x + Float64(y * exp(Float64(2.0 * Float64(a * Float64(-b))))))) tmp = 0.0 if (a <= -54000000.0) tmp = t_1; elseif (a <= -5.2e-166) tmp = Float64(x / Float64(x + Float64(b * Float64(Float64(Float64(y * 1.3333333333333333) + Float64(t * Float64(Float64(-2.0 * Float64(y * Float64(a + 0.8333333333333334))) + Float64(y / b)))) / t)))); elseif (a <= 7.2e-268) tmp = Float64(x / Float64(x + Float64(y * exp(Float64(-1.3333333333333333 * Float64(c / t)))))); elseif (a <= 1.05e-41) tmp = Float64(x / Float64(x + Float64(y * exp(Float64(1.3333333333333333 * Float64(b / t)))))); elseif ((a <= 1.4e+179) || !(a <= 4.6e+280)) tmp = t_1; else tmp = Float64(x / Float64(x + Float64(y * exp(Float64(c * 1.6666666666666667))))); end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c) t_1 = x / (x + (y * exp((2.0 * (a * -b))))); tmp = 0.0; if (a <= -54000000.0) tmp = t_1; elseif (a <= -5.2e-166) tmp = x / (x + (b * (((y * 1.3333333333333333) + (t * ((-2.0 * (y * (a + 0.8333333333333334))) + (y / b)))) / t))); elseif (a <= 7.2e-268) tmp = x / (x + (y * exp((-1.3333333333333333 * (c / t))))); elseif (a <= 1.05e-41) tmp = x / (x + (y * exp((1.3333333333333333 * (b / t))))); elseif ((a <= 1.4e+179) || ~((a <= 4.6e+280))) tmp = t_1; else tmp = x / (x + (y * exp((c * 1.6666666666666667)))); 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[(a * (-b)), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[a, -54000000.0], t$95$1, If[LessEqual[a, -5.2e-166], N[(x / N[(x + N[(b * N[(N[(N[(y * 1.3333333333333333), $MachinePrecision] + N[(t * N[(N[(-2.0 * N[(y * N[(a + 0.8333333333333334), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(y / b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[a, 7.2e-268], N[(x / N[(x + N[(y * N[Exp[N[(-1.3333333333333333 * N[(c / t), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[a, 1.05e-41], N[(x / N[(x + N[(y * N[Exp[N[(1.3333333333333333 * N[(b / t), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[Or[LessEqual[a, 1.4e+179], N[Not[LessEqual[a, 4.6e+280]], $MachinePrecision]], t$95$1, N[(x / N[(x + N[(y * N[Exp[N[(c * 1.6666666666666667), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{x}{x + y \cdot e^{2 \cdot \left(a \cdot \left(-b\right)\right)}}\\
\mathbf{if}\;a \leq -54000000:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;a \leq -5.2 \cdot 10^{-166}:\\
\;\;\;\;\frac{x}{x + b \cdot \frac{y \cdot 1.3333333333333333 + t \cdot \left(-2 \cdot \left(y \cdot \left(a + 0.8333333333333334\right)\right) + \frac{y}{b}\right)}{t}}\\
\mathbf{elif}\;a \leq 7.2 \cdot 10^{-268}:\\
\;\;\;\;\frac{x}{x + y \cdot e^{-1.3333333333333333 \cdot \frac{c}{t}}}\\
\mathbf{elif}\;a \leq 1.05 \cdot 10^{-41}:\\
\;\;\;\;\frac{x}{x + y \cdot e^{1.3333333333333333 \cdot \frac{b}{t}}}\\
\mathbf{elif}\;a \leq 1.4 \cdot 10^{+179} \lor \neg \left(a \leq 4.6 \cdot 10^{+280}\right):\\
\;\;\;\;t\_1\\
\mathbf{else}:\\
\;\;\;\;\frac{x}{x + y \cdot e^{c \cdot 1.6666666666666667}}\\
\end{array}
\end{array}
if a < -5.4e7 or 1.05000000000000006e-41 < a < 1.4e179 or 4.59999999999999999e280 < a Initial program 95.2%
Taylor expanded in b around inf 77.8%
associate-*r/77.8%
metadata-eval77.8%
+-commutative77.8%
Simplified77.8%
Taylor expanded in a around inf 76.0%
associate-*r*76.0%
mul-1-neg76.0%
Simplified76.0%
if -5.4e7 < a < -5.19999999999999979e-166Initial program 100.0%
Taylor expanded in b around inf 66.4%
associate-*r/66.4%
metadata-eval66.4%
+-commutative66.4%
Simplified66.4%
Taylor expanded in b around 0 62.2%
Taylor expanded in b around inf 70.6%
Taylor expanded in t around 0 78.9%
if -5.19999999999999979e-166 < a < 7.2000000000000002e-268Initial program 94.8%
Taylor expanded in t around 0 13.4%
Taylor expanded in z around 0 67.1%
Taylor expanded in b around 0 62.0%
if 7.2000000000000002e-268 < a < 1.05000000000000006e-41Initial program 95.1%
Taylor expanded in t around 0 81.5%
Taylor expanded in z around 0 70.5%
Taylor expanded in c around 0 69.0%
if 1.4e179 < a < 4.59999999999999999e280Initial program 85.7%
Taylor expanded in c around inf 62.1%
+-commutative62.1%
associate-*r/62.1%
metadata-eval62.1%
Simplified62.1%
Taylor expanded in t around inf 69.0%
Taylor expanded in a around 0 72.8%
Final simplification72.1%
(FPCore (x y z t a b c)
:precision binary64
(if (<= y -2.06e+40)
(/
x
(+
x
(*
y
(+
(*
2.0
(/
(- (* b 0.6666666666666666) (* (+ a 0.8333333333333334) (* t b)))
t))
1.0))))
(if (<= y -2.2e-12)
1.0
(if (<= y -1.5e-216)
(/
x
(+
x
(*
y
(+
(*
2.0
(*
a
(-
(* b (/ (+ (/ 0.6666666666666666 t) -0.8333333333333334) a))
b)))
1.0))))
(if (<= y 1.86e-144)
(/ x (+ x (* y (exp (* b -1.6666666666666667)))))
1.0)))))
double code(double x, double y, double z, double t, double a, double b, double c) {
double tmp;
if (y <= -2.06e+40) {
tmp = x / (x + (y * ((2.0 * (((b * 0.6666666666666666) - ((a + 0.8333333333333334) * (t * b))) / t)) + 1.0)));
} else if (y <= -2.2e-12) {
tmp = 1.0;
} else if (y <= -1.5e-216) {
tmp = x / (x + (y * ((2.0 * (a * ((b * (((0.6666666666666666 / t) + -0.8333333333333334) / a)) - b))) + 1.0)));
} else if (y <= 1.86e-144) {
tmp = x / (x + (y * exp((b * -1.6666666666666667))));
} 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 (y <= (-2.06d+40)) then
tmp = x / (x + (y * ((2.0d0 * (((b * 0.6666666666666666d0) - ((a + 0.8333333333333334d0) * (t * b))) / t)) + 1.0d0)))
else if (y <= (-2.2d-12)) then
tmp = 1.0d0
else if (y <= (-1.5d-216)) then
tmp = x / (x + (y * ((2.0d0 * (a * ((b * (((0.6666666666666666d0 / t) + (-0.8333333333333334d0)) / a)) - b))) + 1.0d0)))
else if (y <= 1.86d-144) then
tmp = x / (x + (y * exp((b * (-1.6666666666666667d0)))))
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 (y <= -2.06e+40) {
tmp = x / (x + (y * ((2.0 * (((b * 0.6666666666666666) - ((a + 0.8333333333333334) * (t * b))) / t)) + 1.0)));
} else if (y <= -2.2e-12) {
tmp = 1.0;
} else if (y <= -1.5e-216) {
tmp = x / (x + (y * ((2.0 * (a * ((b * (((0.6666666666666666 / t) + -0.8333333333333334) / a)) - b))) + 1.0)));
} else if (y <= 1.86e-144) {
tmp = x / (x + (y * Math.exp((b * -1.6666666666666667))));
} else {
tmp = 1.0;
}
return tmp;
}
def code(x, y, z, t, a, b, c): tmp = 0 if y <= -2.06e+40: tmp = x / (x + (y * ((2.0 * (((b * 0.6666666666666666) - ((a + 0.8333333333333334) * (t * b))) / t)) + 1.0))) elif y <= -2.2e-12: tmp = 1.0 elif y <= -1.5e-216: tmp = x / (x + (y * ((2.0 * (a * ((b * (((0.6666666666666666 / t) + -0.8333333333333334) / a)) - b))) + 1.0))) elif y <= 1.86e-144: tmp = x / (x + (y * math.exp((b * -1.6666666666666667)))) else: tmp = 1.0 return tmp
function code(x, y, z, t, a, b, c) tmp = 0.0 if (y <= -2.06e+40) tmp = Float64(x / Float64(x + Float64(y * Float64(Float64(2.0 * Float64(Float64(Float64(b * 0.6666666666666666) - Float64(Float64(a + 0.8333333333333334) * Float64(t * b))) / t)) + 1.0)))); elseif (y <= -2.2e-12) tmp = 1.0; elseif (y <= -1.5e-216) tmp = Float64(x / Float64(x + Float64(y * Float64(Float64(2.0 * Float64(a * Float64(Float64(b * Float64(Float64(Float64(0.6666666666666666 / t) + -0.8333333333333334) / a)) - b))) + 1.0)))); elseif (y <= 1.86e-144) tmp = Float64(x / Float64(x + Float64(y * exp(Float64(b * -1.6666666666666667))))); else tmp = 1.0; end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c) tmp = 0.0; if (y <= -2.06e+40) tmp = x / (x + (y * ((2.0 * (((b * 0.6666666666666666) - ((a + 0.8333333333333334) * (t * b))) / t)) + 1.0))); elseif (y <= -2.2e-12) tmp = 1.0; elseif (y <= -1.5e-216) tmp = x / (x + (y * ((2.0 * (a * ((b * (((0.6666666666666666 / t) + -0.8333333333333334) / a)) - b))) + 1.0))); elseif (y <= 1.86e-144) tmp = x / (x + (y * exp((b * -1.6666666666666667)))); else tmp = 1.0; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_] := If[LessEqual[y, -2.06e+40], N[(x / N[(x + N[(y * N[(N[(2.0 * N[(N[(N[(b * 0.6666666666666666), $MachinePrecision] - N[(N[(a + 0.8333333333333334), $MachinePrecision] * N[(t * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / t), $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, -2.2e-12], 1.0, If[LessEqual[y, -1.5e-216], N[(x / N[(x + N[(y * N[(N[(2.0 * N[(a * N[(N[(b * N[(N[(N[(0.6666666666666666 / t), $MachinePrecision] + -0.8333333333333334), $MachinePrecision] / a), $MachinePrecision]), $MachinePrecision] - b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 1.86e-144], N[(x / N[(x + N[(y * N[Exp[N[(b * -1.6666666666666667), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 1.0]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -2.06 \cdot 10^{+40}:\\
\;\;\;\;\frac{x}{x + y \cdot \left(2 \cdot \frac{b \cdot 0.6666666666666666 - \left(a + 0.8333333333333334\right) \cdot \left(t \cdot b\right)}{t} + 1\right)}\\
\mathbf{elif}\;y \leq -2.2 \cdot 10^{-12}:\\
\;\;\;\;1\\
\mathbf{elif}\;y \leq -1.5 \cdot 10^{-216}:\\
\;\;\;\;\frac{x}{x + y \cdot \left(2 \cdot \left(a \cdot \left(b \cdot \frac{\frac{0.6666666666666666}{t} + -0.8333333333333334}{a} - b\right)\right) + 1\right)}\\
\mathbf{elif}\;y \leq 1.86 \cdot 10^{-144}:\\
\;\;\;\;\frac{x}{x + y \cdot e^{b \cdot -1.6666666666666667}}\\
\mathbf{else}:\\
\;\;\;\;1\\
\end{array}
\end{array}
if y < -2.05999999999999999e40Initial program 93.7%
Taylor expanded in b around inf 75.7%
associate-*r/75.7%
metadata-eval75.7%
+-commutative75.7%
Simplified75.7%
Taylor expanded in b around 0 67.8%
Taylor expanded in t around 0 71.7%
+-commutative71.7%
mul-1-neg71.7%
unsub-neg71.7%
*-commutative71.7%
associate-*r*71.7%
Simplified71.7%
if -2.05999999999999999e40 < y < -2.19999999999999992e-12 or 1.8599999999999999e-144 < y Initial program 94.7%
Taylor expanded in c around inf 71.7%
+-commutative71.7%
associate-*r/71.7%
metadata-eval71.7%
Simplified71.7%
Taylor expanded in t around inf 58.8%
Taylor expanded in x around inf 61.4%
if -2.19999999999999992e-12 < y < -1.50000000000000006e-216Initial program 91.4%
Taylor expanded in b around inf 64.4%
associate-*r/64.4%
metadata-eval64.4%
+-commutative64.4%
Simplified64.4%
Taylor expanded in b around 0 45.8%
Taylor expanded in a around inf 60.3%
+-commutative60.3%
mul-1-neg60.3%
unsub-neg60.3%
associate-/l*64.3%
sub-neg64.3%
associate-*r/64.3%
metadata-eval64.3%
metadata-eval64.3%
Simplified64.3%
if -1.50000000000000006e-216 < y < 1.8599999999999999e-144Initial program 97.9%
Taylor expanded in b around inf 78.5%
associate-*r/78.5%
metadata-eval78.5%
+-commutative78.5%
Simplified78.5%
Taylor expanded in t around inf 70.6%
mul-1-neg70.6%
+-commutative70.6%
distribute-rgt-neg-in70.6%
+-commutative70.6%
mul-1-neg70.6%
distribute-lft-in70.6%
metadata-eval70.6%
mul-1-neg70.6%
unsub-neg70.6%
Simplified70.6%
Taylor expanded in a around 0 68.6%
Final simplification65.2%
(FPCore (x y z t a b c)
:precision binary64
(let* ((t_1 (/ x (+ x (* y (exp (* b -1.6666666666666667)))))))
(if (<= t -4.9e-213)
t_1
(if (<= t 4.5e-64)
(/ x (+ x (* y (exp (* -1.3333333333333333 (/ c t))))))
(if (<= t 0.8)
1.0
(if (<= t 5.5e+174)
t_1
(/ x (+ x (* y (exp (* c 1.6666666666666667)))))))))))
double code(double x, double y, double z, double t, double a, double b, double c) {
double t_1 = x / (x + (y * exp((b * -1.6666666666666667))));
double tmp;
if (t <= -4.9e-213) {
tmp = t_1;
} else if (t <= 4.5e-64) {
tmp = x / (x + (y * exp((-1.3333333333333333 * (c / t)))));
} else if (t <= 0.8) {
tmp = 1.0;
} else if (t <= 5.5e+174) {
tmp = t_1;
} else {
tmp = x / (x + (y * exp((c * 1.6666666666666667))));
}
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((b * (-1.6666666666666667d0)))))
if (t <= (-4.9d-213)) then
tmp = t_1
else if (t <= 4.5d-64) then
tmp = x / (x + (y * exp(((-1.3333333333333333d0) * (c / t)))))
else if (t <= 0.8d0) then
tmp = 1.0d0
else if (t <= 5.5d+174) then
tmp = t_1
else
tmp = x / (x + (y * exp((c * 1.6666666666666667d0))))
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((b * -1.6666666666666667))));
double tmp;
if (t <= -4.9e-213) {
tmp = t_1;
} else if (t <= 4.5e-64) {
tmp = x / (x + (y * Math.exp((-1.3333333333333333 * (c / t)))));
} else if (t <= 0.8) {
tmp = 1.0;
} else if (t <= 5.5e+174) {
tmp = t_1;
} else {
tmp = x / (x + (y * Math.exp((c * 1.6666666666666667))));
}
return tmp;
}
def code(x, y, z, t, a, b, c): t_1 = x / (x + (y * math.exp((b * -1.6666666666666667)))) tmp = 0 if t <= -4.9e-213: tmp = t_1 elif t <= 4.5e-64: tmp = x / (x + (y * math.exp((-1.3333333333333333 * (c / t))))) elif t <= 0.8: tmp = 1.0 elif t <= 5.5e+174: tmp = t_1 else: tmp = x / (x + (y * math.exp((c * 1.6666666666666667)))) return tmp
function code(x, y, z, t, a, b, c) t_1 = Float64(x / Float64(x + Float64(y * exp(Float64(b * -1.6666666666666667))))) tmp = 0.0 if (t <= -4.9e-213) tmp = t_1; elseif (t <= 4.5e-64) tmp = Float64(x / Float64(x + Float64(y * exp(Float64(-1.3333333333333333 * Float64(c / t)))))); elseif (t <= 0.8) tmp = 1.0; elseif (t <= 5.5e+174) tmp = t_1; else tmp = Float64(x / Float64(x + Float64(y * exp(Float64(c * 1.6666666666666667))))); end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c) t_1 = x / (x + (y * exp((b * -1.6666666666666667)))); tmp = 0.0; if (t <= -4.9e-213) tmp = t_1; elseif (t <= 4.5e-64) tmp = x / (x + (y * exp((-1.3333333333333333 * (c / t))))); elseif (t <= 0.8) tmp = 1.0; elseif (t <= 5.5e+174) tmp = t_1; else tmp = x / (x + (y * exp((c * 1.6666666666666667)))); 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[(b * -1.6666666666666667), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t, -4.9e-213], t$95$1, If[LessEqual[t, 4.5e-64], N[(x / N[(x + N[(y * N[Exp[N[(-1.3333333333333333 * N[(c / t), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t, 0.8], 1.0, If[LessEqual[t, 5.5e+174], t$95$1, N[(x / N[(x + N[(y * N[Exp[N[(c * 1.6666666666666667), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{x}{x + y \cdot e^{b \cdot -1.6666666666666667}}\\
\mathbf{if}\;t \leq -4.9 \cdot 10^{-213}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;t \leq 4.5 \cdot 10^{-64}:\\
\;\;\;\;\frac{x}{x + y \cdot e^{-1.3333333333333333 \cdot \frac{c}{t}}}\\
\mathbf{elif}\;t \leq 0.8:\\
\;\;\;\;1\\
\mathbf{elif}\;t \leq 5.5 \cdot 10^{+174}:\\
\;\;\;\;t\_1\\
\mathbf{else}:\\
\;\;\;\;\frac{x}{x + y \cdot e^{c \cdot 1.6666666666666667}}\\
\end{array}
\end{array}
if t < -4.8999999999999998e-213 or 0.80000000000000004 < t < 5.4999999999999998e174Initial program 97.1%
Taylor expanded in b around inf 69.3%
associate-*r/69.3%
metadata-eval69.3%
+-commutative69.3%
Simplified69.3%
Taylor expanded in t around inf 69.3%
mul-1-neg69.3%
+-commutative69.3%
distribute-rgt-neg-in69.3%
+-commutative69.3%
mul-1-neg69.3%
distribute-lft-in69.3%
metadata-eval69.3%
mul-1-neg69.3%
unsub-neg69.3%
Simplified69.3%
Taylor expanded in a around 0 68.3%
if -4.8999999999999998e-213 < t < 4.5000000000000001e-64Initial program 86.1%
Taylor expanded in t around 0 76.6%
Taylor expanded in z around 0 75.8%
Taylor expanded in b around 0 67.8%
if 4.5000000000000001e-64 < t < 0.80000000000000004Initial program 100.0%
Taylor expanded in c around inf 68.9%
+-commutative68.9%
associate-*r/68.9%
metadata-eval68.9%
Simplified68.9%
Taylor expanded in t around inf 51.6%
Taylor expanded in x around inf 72.3%
if 5.4999999999999998e174 < t Initial program 98.0%
Taylor expanded in c around inf 76.5%
+-commutative76.5%
associate-*r/76.5%
metadata-eval76.5%
Simplified76.5%
Taylor expanded in t around inf 76.5%
Taylor expanded in a around 0 66.5%
Final simplification68.2%
(FPCore (x y z t a b c)
:precision binary64
(if (<= t -5e-213)
(/ x (+ x (* y (exp (* 1.3333333333333333 (/ b t))))))
(if (<= t 3.2e-151)
(/ x (+ x (* y (exp (* -1.3333333333333333 (/ c t))))))
(if (<= t 2.6e-126)
(/
x
(-
x
(*
y
(+
-1.0
(*
2.0
(*
a
(-
b
(*
b
(/ (+ (/ 0.6666666666666666 t) -0.8333333333333334) a)))))))))
(if (<= t 1.65e-10)
1.0
(/ x (+ x (* y (exp (* c 1.6666666666666667))))))))))
double code(double x, double y, double z, double t, double a, double b, double c) {
double tmp;
if (t <= -5e-213) {
tmp = x / (x + (y * exp((1.3333333333333333 * (b / t)))));
} else if (t <= 3.2e-151) {
tmp = x / (x + (y * exp((-1.3333333333333333 * (c / t)))));
} else if (t <= 2.6e-126) {
tmp = x / (x - (y * (-1.0 + (2.0 * (a * (b - (b * (((0.6666666666666666 / t) + -0.8333333333333334) / a))))))));
} else if (t <= 1.65e-10) {
tmp = 1.0;
} else {
tmp = x / (x + (y * exp((c * 1.6666666666666667))));
}
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 <= (-5d-213)) then
tmp = x / (x + (y * exp((1.3333333333333333d0 * (b / t)))))
else if (t <= 3.2d-151) then
tmp = x / (x + (y * exp(((-1.3333333333333333d0) * (c / t)))))
else if (t <= 2.6d-126) then
tmp = x / (x - (y * ((-1.0d0) + (2.0d0 * (a * (b - (b * (((0.6666666666666666d0 / t) + (-0.8333333333333334d0)) / a))))))))
else if (t <= 1.65d-10) then
tmp = 1.0d0
else
tmp = x / (x + (y * exp((c * 1.6666666666666667d0))))
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 <= -5e-213) {
tmp = x / (x + (y * Math.exp((1.3333333333333333 * (b / t)))));
} else if (t <= 3.2e-151) {
tmp = x / (x + (y * Math.exp((-1.3333333333333333 * (c / t)))));
} else if (t <= 2.6e-126) {
tmp = x / (x - (y * (-1.0 + (2.0 * (a * (b - (b * (((0.6666666666666666 / t) + -0.8333333333333334) / a))))))));
} else if (t <= 1.65e-10) {
tmp = 1.0;
} else {
tmp = x / (x + (y * Math.exp((c * 1.6666666666666667))));
}
return tmp;
}
def code(x, y, z, t, a, b, c): tmp = 0 if t <= -5e-213: tmp = x / (x + (y * math.exp((1.3333333333333333 * (b / t))))) elif t <= 3.2e-151: tmp = x / (x + (y * math.exp((-1.3333333333333333 * (c / t))))) elif t <= 2.6e-126: tmp = x / (x - (y * (-1.0 + (2.0 * (a * (b - (b * (((0.6666666666666666 / t) + -0.8333333333333334) / a)))))))) elif t <= 1.65e-10: tmp = 1.0 else: tmp = x / (x + (y * math.exp((c * 1.6666666666666667)))) return tmp
function code(x, y, z, t, a, b, c) tmp = 0.0 if (t <= -5e-213) tmp = Float64(x / Float64(x + Float64(y * exp(Float64(1.3333333333333333 * Float64(b / t)))))); elseif (t <= 3.2e-151) tmp = Float64(x / Float64(x + Float64(y * exp(Float64(-1.3333333333333333 * Float64(c / t)))))); elseif (t <= 2.6e-126) tmp = Float64(x / Float64(x - Float64(y * Float64(-1.0 + Float64(2.0 * Float64(a * Float64(b - Float64(b * Float64(Float64(Float64(0.6666666666666666 / t) + -0.8333333333333334) / a))))))))); elseif (t <= 1.65e-10) tmp = 1.0; else tmp = Float64(x / Float64(x + Float64(y * exp(Float64(c * 1.6666666666666667))))); end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c) tmp = 0.0; if (t <= -5e-213) tmp = x / (x + (y * exp((1.3333333333333333 * (b / t))))); elseif (t <= 3.2e-151) tmp = x / (x + (y * exp((-1.3333333333333333 * (c / t))))); elseif (t <= 2.6e-126) tmp = x / (x - (y * (-1.0 + (2.0 * (a * (b - (b * (((0.6666666666666666 / t) + -0.8333333333333334) / a)))))))); elseif (t <= 1.65e-10) tmp = 1.0; else tmp = x / (x + (y * exp((c * 1.6666666666666667)))); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_] := If[LessEqual[t, -5e-213], N[(x / N[(x + N[(y * N[Exp[N[(1.3333333333333333 * N[(b / t), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t, 3.2e-151], N[(x / N[(x + N[(y * N[Exp[N[(-1.3333333333333333 * N[(c / t), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t, 2.6e-126], N[(x / N[(x - N[(y * N[(-1.0 + N[(2.0 * N[(a * N[(b - N[(b * N[(N[(N[(0.6666666666666666 / t), $MachinePrecision] + -0.8333333333333334), $MachinePrecision] / a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t, 1.65e-10], 1.0, N[(x / N[(x + N[(y * N[Exp[N[(c * 1.6666666666666667), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;t \leq -5 \cdot 10^{-213}:\\
\;\;\;\;\frac{x}{x + y \cdot e^{1.3333333333333333 \cdot \frac{b}{t}}}\\
\mathbf{elif}\;t \leq 3.2 \cdot 10^{-151}:\\
\;\;\;\;\frac{x}{x + y \cdot e^{-1.3333333333333333 \cdot \frac{c}{t}}}\\
\mathbf{elif}\;t \leq 2.6 \cdot 10^{-126}:\\
\;\;\;\;\frac{x}{x - y \cdot \left(-1 + 2 \cdot \left(a \cdot \left(b - b \cdot \frac{\frac{0.6666666666666666}{t} + -0.8333333333333334}{a}\right)\right)\right)}\\
\mathbf{elif}\;t \leq 1.65 \cdot 10^{-10}:\\
\;\;\;\;1\\
\mathbf{else}:\\
\;\;\;\;\frac{x}{x + y \cdot e^{c \cdot 1.6666666666666667}}\\
\end{array}
\end{array}
if t < -4.99999999999999977e-213Initial program 97.6%
Taylor expanded in t around 0 93.2%
Taylor expanded in z around 0 88.7%
Taylor expanded in c around 0 77.7%
if -4.99999999999999977e-213 < t < 3.20000000000000021e-151Initial program 82.4%
Taylor expanded in t around 0 88.4%
Taylor expanded in z around 0 73.4%
Taylor expanded in b around 0 69.7%
if 3.20000000000000021e-151 < t < 2.59999999999999999e-126Initial program 100.0%
Taylor expanded in b around inf 89.2%
associate-*r/89.2%
metadata-eval89.2%
+-commutative89.2%
Simplified89.2%
Taylor expanded in b around 0 57.2%
Taylor expanded in a around inf 78.6%
+-commutative78.6%
mul-1-neg78.6%
unsub-neg78.6%
associate-/l*89.2%
sub-neg89.2%
associate-*r/89.2%
metadata-eval89.2%
metadata-eval89.2%
Simplified89.2%
if 2.59999999999999999e-126 < t < 1.65e-10Initial program 97.2%
Taylor expanded in c around inf 73.1%
+-commutative73.1%
associate-*r/73.1%
metadata-eval73.1%
Simplified73.1%
Taylor expanded in t around inf 43.5%
Taylor expanded in x around inf 70.4%
if 1.65e-10 < t Initial program 97.4%
Taylor expanded in c around inf 71.0%
+-commutative71.0%
associate-*r/71.0%
metadata-eval71.0%
Simplified71.0%
Taylor expanded in t around inf 72.6%
Taylor expanded in a around 0 64.3%
Final simplification69.3%
(FPCore (x y z t a b c)
:precision binary64
(if (<= t 0.04)
(/ x (+ x (* y (exp (* 1.3333333333333333 (/ (- b c) t))))))
(if (or (<= t 1.55e+190) (not (<= t 3.4e+248)))
(/ x (+ x (* y (exp (* 2.0 (* b (- -0.8333333333333334 a)))))))
(/ x (+ x (* y (exp (* c 1.6666666666666667))))))))
double code(double x, double y, double z, double t, double a, double b, double c) {
double tmp;
if (t <= 0.04) {
tmp = x / (x + (y * exp((1.3333333333333333 * ((b - c) / t)))));
} else if ((t <= 1.55e+190) || !(t <= 3.4e+248)) {
tmp = x / (x + (y * exp((2.0 * (b * (-0.8333333333333334 - a))))));
} else {
tmp = x / (x + (y * exp((c * 1.6666666666666667))));
}
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 <= 0.04d0) then
tmp = x / (x + (y * exp((1.3333333333333333d0 * ((b - c) / t)))))
else if ((t <= 1.55d+190) .or. (.not. (t <= 3.4d+248))) then
tmp = x / (x + (y * exp((2.0d0 * (b * ((-0.8333333333333334d0) - a))))))
else
tmp = x / (x + (y * exp((c * 1.6666666666666667d0))))
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 <= 0.04) {
tmp = x / (x + (y * Math.exp((1.3333333333333333 * ((b - c) / t)))));
} else if ((t <= 1.55e+190) || !(t <= 3.4e+248)) {
tmp = x / (x + (y * Math.exp((2.0 * (b * (-0.8333333333333334 - a))))));
} else {
tmp = x / (x + (y * Math.exp((c * 1.6666666666666667))));
}
return tmp;
}
def code(x, y, z, t, a, b, c): tmp = 0 if t <= 0.04: tmp = x / (x + (y * math.exp((1.3333333333333333 * ((b - c) / t))))) elif (t <= 1.55e+190) or not (t <= 3.4e+248): tmp = x / (x + (y * math.exp((2.0 * (b * (-0.8333333333333334 - a)))))) else: tmp = x / (x + (y * math.exp((c * 1.6666666666666667)))) return tmp
function code(x, y, z, t, a, b, c) tmp = 0.0 if (t <= 0.04) tmp = Float64(x / Float64(x + Float64(y * exp(Float64(1.3333333333333333 * Float64(Float64(b - c) / t)))))); elseif ((t <= 1.55e+190) || !(t <= 3.4e+248)) tmp = Float64(x / Float64(x + Float64(y * exp(Float64(2.0 * Float64(b * Float64(-0.8333333333333334 - a))))))); else tmp = Float64(x / Float64(x + Float64(y * exp(Float64(c * 1.6666666666666667))))); end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c) tmp = 0.0; if (t <= 0.04) tmp = x / (x + (y * exp((1.3333333333333333 * ((b - c) / t))))); elseif ((t <= 1.55e+190) || ~((t <= 3.4e+248))) tmp = x / (x + (y * exp((2.0 * (b * (-0.8333333333333334 - a)))))); else tmp = x / (x + (y * exp((c * 1.6666666666666667)))); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_] := If[LessEqual[t, 0.04], N[(x / N[(x + N[(y * N[Exp[N[(1.3333333333333333 * N[(N[(b - c), $MachinePrecision] / t), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[Or[LessEqual[t, 1.55e+190], N[Not[LessEqual[t, 3.4e+248]], $MachinePrecision]], N[(x / N[(x + N[(y * N[Exp[N[(2.0 * N[(b * N[(-0.8333333333333334 - a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x / N[(x + N[(y * N[Exp[N[(c * 1.6666666666666667), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;t \leq 0.04:\\
\;\;\;\;\frac{x}{x + y \cdot e^{1.3333333333333333 \cdot \frac{b - c}{t}}}\\
\mathbf{elif}\;t \leq 1.55 \cdot 10^{+190} \lor \neg \left(t \leq 3.4 \cdot 10^{+248}\right):\\
\;\;\;\;\frac{x}{x + y \cdot e^{2 \cdot \left(b \cdot \left(-0.8333333333333334 - a\right)\right)}}\\
\mathbf{else}:\\
\;\;\;\;\frac{x}{x + y \cdot e^{c \cdot 1.6666666666666667}}\\
\end{array}
\end{array}
if t < 0.0400000000000000008Initial program 92.3%
Taylor expanded in t around 0 75.2%
Taylor expanded in z around 0 77.0%
if 0.0400000000000000008 < t < 1.5500000000000001e190 or 3.3999999999999999e248 < t Initial program 97.7%
Taylor expanded in b around inf 71.1%
associate-*r/71.1%
metadata-eval71.1%
+-commutative71.1%
Simplified71.1%
Taylor expanded in t around inf 71.1%
mul-1-neg71.1%
+-commutative71.1%
distribute-rgt-neg-in71.1%
+-commutative71.1%
mul-1-neg71.1%
distribute-lft-in71.1%
metadata-eval71.1%
mul-1-neg71.1%
unsub-neg71.1%
Simplified71.1%
if 1.5500000000000001e190 < t < 3.3999999999999999e248Initial program 95.8%
Taylor expanded in c around inf 83.3%
+-commutative83.3%
associate-*r/83.3%
metadata-eval83.3%
Simplified83.3%
Taylor expanded in t around inf 83.3%
Taylor expanded in a around 0 83.3%
Final simplification75.5%
(FPCore (x y z t a b c)
:precision binary64
(if (<= t 8e-13)
(/ x (+ x (* y (exp (* 1.3333333333333333 (/ (- b c) t))))))
(if (or (<= t 2.3e+131) (not (<= t 1.32e+188)))
(/ x (+ x (* y (exp (* 2.0 (* c (+ a 0.8333333333333334)))))))
(/ x (+ x (* y (exp (* 2.0 (* b (- -0.8333333333333334 a))))))))))
double code(double x, double y, double z, double t, double a, double b, double c) {
double tmp;
if (t <= 8e-13) {
tmp = x / (x + (y * exp((1.3333333333333333 * ((b - c) / t)))));
} else if ((t <= 2.3e+131) || !(t <= 1.32e+188)) {
tmp = x / (x + (y * exp((2.0 * (c * (a + 0.8333333333333334))))));
} else {
tmp = x / (x + (y * exp((2.0 * (b * (-0.8333333333333334 - a))))));
}
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 <= 8d-13) then
tmp = x / (x + (y * exp((1.3333333333333333d0 * ((b - c) / t)))))
else if ((t <= 2.3d+131) .or. (.not. (t <= 1.32d+188))) then
tmp = x / (x + (y * exp((2.0d0 * (c * (a + 0.8333333333333334d0))))))
else
tmp = x / (x + (y * exp((2.0d0 * (b * ((-0.8333333333333334d0) - a))))))
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 <= 8e-13) {
tmp = x / (x + (y * Math.exp((1.3333333333333333 * ((b - c) / t)))));
} else if ((t <= 2.3e+131) || !(t <= 1.32e+188)) {
tmp = x / (x + (y * Math.exp((2.0 * (c * (a + 0.8333333333333334))))));
} else {
tmp = x / (x + (y * Math.exp((2.0 * (b * (-0.8333333333333334 - a))))));
}
return tmp;
}
def code(x, y, z, t, a, b, c): tmp = 0 if t <= 8e-13: tmp = x / (x + (y * math.exp((1.3333333333333333 * ((b - c) / t))))) elif (t <= 2.3e+131) or not (t <= 1.32e+188): tmp = x / (x + (y * math.exp((2.0 * (c * (a + 0.8333333333333334)))))) else: tmp = x / (x + (y * math.exp((2.0 * (b * (-0.8333333333333334 - a)))))) return tmp
function code(x, y, z, t, a, b, c) tmp = 0.0 if (t <= 8e-13) tmp = Float64(x / Float64(x + Float64(y * exp(Float64(1.3333333333333333 * Float64(Float64(b - c) / t)))))); elseif ((t <= 2.3e+131) || !(t <= 1.32e+188)) tmp = Float64(x / Float64(x + Float64(y * exp(Float64(2.0 * Float64(c * Float64(a + 0.8333333333333334))))))); else tmp = Float64(x / Float64(x + Float64(y * exp(Float64(2.0 * Float64(b * Float64(-0.8333333333333334 - a))))))); end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c) tmp = 0.0; if (t <= 8e-13) tmp = x / (x + (y * exp((1.3333333333333333 * ((b - c) / t))))); elseif ((t <= 2.3e+131) || ~((t <= 1.32e+188))) tmp = x / (x + (y * exp((2.0 * (c * (a + 0.8333333333333334)))))); else tmp = x / (x + (y * exp((2.0 * (b * (-0.8333333333333334 - a)))))); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_] := If[LessEqual[t, 8e-13], N[(x / N[(x + N[(y * N[Exp[N[(1.3333333333333333 * N[(N[(b - c), $MachinePrecision] / t), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[Or[LessEqual[t, 2.3e+131], N[Not[LessEqual[t, 1.32e+188]], $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], N[(x / N[(x + N[(y * N[Exp[N[(2.0 * N[(b * N[(-0.8333333333333334 - a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;t \leq 8 \cdot 10^{-13}:\\
\;\;\;\;\frac{x}{x + y \cdot e^{1.3333333333333333 \cdot \frac{b - c}{t}}}\\
\mathbf{elif}\;t \leq 2.3 \cdot 10^{+131} \lor \neg \left(t \leq 1.32 \cdot 10^{+188}\right):\\
\;\;\;\;\frac{x}{x + y \cdot e^{2 \cdot \left(c \cdot \left(a + 0.8333333333333334\right)\right)}}\\
\mathbf{else}:\\
\;\;\;\;\frac{x}{x + y \cdot e^{2 \cdot \left(b \cdot \left(-0.8333333333333334 - a\right)\right)}}\\
\end{array}
\end{array}
if t < 8.0000000000000002e-13Initial program 92.0%
Taylor expanded in t around 0 76.4%
Taylor expanded in z around 0 78.2%
if 8.0000000000000002e-13 < t < 2.29999999999999992e131 or 1.3200000000000001e188 < t Initial program 97.9%
Taylor expanded in c around inf 75.0%
+-commutative75.0%
associate-*r/75.0%
metadata-eval75.0%
Simplified75.0%
Taylor expanded in t around inf 77.0%
if 2.29999999999999992e131 < t < 1.3200000000000001e188Initial program 95.4%
Taylor expanded in b around inf 77.7%
associate-*r/77.7%
metadata-eval77.7%
+-commutative77.7%
Simplified77.7%
Taylor expanded in t around inf 77.7%
mul-1-neg77.7%
+-commutative77.7%
distribute-rgt-neg-in77.7%
+-commutative77.7%
mul-1-neg77.7%
distribute-lft-in77.7%
metadata-eval77.7%
mul-1-neg77.7%
unsub-neg77.7%
Simplified77.7%
Final simplification77.7%
(FPCore (x y z t a b c)
:precision binary64
(if (<= t 2.3e-15)
(/ x (+ x (* y (exp (* 1.3333333333333333 (/ (- b c) t))))))
(if (or (<= t 3.3e+131) (not (<= t 1.22e+190)))
(/ x (+ x (* y (exp (* c 1.6666666666666667)))))
(/ x (+ x (* y (exp (* 2.0 (* a (- b))))))))))
double code(double x, double y, double z, double t, double a, double b, double c) {
double tmp;
if (t <= 2.3e-15) {
tmp = x / (x + (y * exp((1.3333333333333333 * ((b - c) / t)))));
} else if ((t <= 3.3e+131) || !(t <= 1.22e+190)) {
tmp = x / (x + (y * exp((c * 1.6666666666666667))));
} else {
tmp = x / (x + (y * exp((2.0 * (a * -b)))));
}
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.3d-15) then
tmp = x / (x + (y * exp((1.3333333333333333d0 * ((b - c) / t)))))
else if ((t <= 3.3d+131) .or. (.not. (t <= 1.22d+190))) then
tmp = x / (x + (y * exp((c * 1.6666666666666667d0))))
else
tmp = x / (x + (y * exp((2.0d0 * (a * -b)))))
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.3e-15) {
tmp = x / (x + (y * Math.exp((1.3333333333333333 * ((b - c) / t)))));
} else if ((t <= 3.3e+131) || !(t <= 1.22e+190)) {
tmp = x / (x + (y * Math.exp((c * 1.6666666666666667))));
} else {
tmp = x / (x + (y * Math.exp((2.0 * (a * -b)))));
}
return tmp;
}
def code(x, y, z, t, a, b, c): tmp = 0 if t <= 2.3e-15: tmp = x / (x + (y * math.exp((1.3333333333333333 * ((b - c) / t))))) elif (t <= 3.3e+131) or not (t <= 1.22e+190): tmp = x / (x + (y * math.exp((c * 1.6666666666666667)))) else: tmp = x / (x + (y * math.exp((2.0 * (a * -b))))) return tmp
function code(x, y, z, t, a, b, c) tmp = 0.0 if (t <= 2.3e-15) tmp = Float64(x / Float64(x + Float64(y * exp(Float64(1.3333333333333333 * Float64(Float64(b - c) / t)))))); elseif ((t <= 3.3e+131) || !(t <= 1.22e+190)) tmp = Float64(x / Float64(x + Float64(y * exp(Float64(c * 1.6666666666666667))))); else tmp = Float64(x / Float64(x + Float64(y * exp(Float64(2.0 * Float64(a * Float64(-b))))))); end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c) tmp = 0.0; if (t <= 2.3e-15) tmp = x / (x + (y * exp((1.3333333333333333 * ((b - c) / t))))); elseif ((t <= 3.3e+131) || ~((t <= 1.22e+190))) tmp = x / (x + (y * exp((c * 1.6666666666666667)))); else tmp = x / (x + (y * exp((2.0 * (a * -b))))); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_] := If[LessEqual[t, 2.3e-15], N[(x / N[(x + N[(y * N[Exp[N[(1.3333333333333333 * N[(N[(b - c), $MachinePrecision] / t), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[Or[LessEqual[t, 3.3e+131], N[Not[LessEqual[t, 1.22e+190]], $MachinePrecision]], N[(x / N[(x + N[(y * N[Exp[N[(c * 1.6666666666666667), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x / N[(x + N[(y * N[Exp[N[(2.0 * N[(a * (-b)), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;t \leq 2.3 \cdot 10^{-15}:\\
\;\;\;\;\frac{x}{x + y \cdot e^{1.3333333333333333 \cdot \frac{b - c}{t}}}\\
\mathbf{elif}\;t \leq 3.3 \cdot 10^{+131} \lor \neg \left(t \leq 1.22 \cdot 10^{+190}\right):\\
\;\;\;\;\frac{x}{x + y \cdot e^{c \cdot 1.6666666666666667}}\\
\mathbf{else}:\\
\;\;\;\;\frac{x}{x + y \cdot e^{2 \cdot \left(a \cdot \left(-b\right)\right)}}\\
\end{array}
\end{array}
if t < 2.2999999999999999e-15Initial program 91.9%
Taylor expanded in t around 0 76.8%
Taylor expanded in z around 0 78.6%
if 2.2999999999999999e-15 < t < 3.2999999999999998e131 or 1.21999999999999995e190 < t Initial program 98.0%
Taylor expanded in c around inf 73.5%
+-commutative73.5%
associate-*r/73.5%
metadata-eval73.5%
Simplified73.5%
Taylor expanded in t around inf 75.5%
Taylor expanded in a around 0 66.6%
if 3.2999999999999998e131 < t < 1.21999999999999995e190Initial program 95.4%
Taylor expanded in b around inf 77.7%
associate-*r/77.7%
metadata-eval77.7%
+-commutative77.7%
Simplified77.7%
Taylor expanded in a around inf 73.2%
associate-*r*73.2%
mul-1-neg73.2%
Simplified73.2%
Final simplification73.5%
(FPCore (x y z t a b c)
:precision binary64
(if (or (<= b -4.1e+19) (not (<= b 7.2e-152)))
(/
x
(+
x
(*
y
(exp
(* 2.0 (* b (- (/ 0.6666666666666666 t) (+ a 0.8333333333333334))))))))
(/ 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 ((b <= -4.1e+19) || !(b <= 7.2e-152)) {
tmp = x / (x + (y * exp((2.0 * (b * ((0.6666666666666666 / t) - (a + 0.8333333333333334)))))));
} 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 ((b <= (-4.1d+19)) .or. (.not. (b <= 7.2d-152))) then
tmp = x / (x + (y * exp((2.0d0 * (b * ((0.6666666666666666d0 / t) - (a + 0.8333333333333334d0)))))))
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 ((b <= -4.1e+19) || !(b <= 7.2e-152)) {
tmp = x / (x + (y * Math.exp((2.0 * (b * ((0.6666666666666666 / t) - (a + 0.8333333333333334)))))));
} 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 (b <= -4.1e+19) or not (b <= 7.2e-152): tmp = x / (x + (y * math.exp((2.0 * (b * ((0.6666666666666666 / t) - (a + 0.8333333333333334))))))) 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 ((b <= -4.1e+19) || !(b <= 7.2e-152)) tmp = Float64(x / Float64(x + Float64(y * exp(Float64(2.0 * Float64(b * Float64(Float64(0.6666666666666666 / t) - Float64(a + 0.8333333333333334)))))))); 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 ((b <= -4.1e+19) || ~((b <= 7.2e-152))) tmp = x / (x + (y * exp((2.0 * (b * ((0.6666666666666666 / t) - (a + 0.8333333333333334))))))); 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[Or[LessEqual[b, -4.1e+19], N[Not[LessEqual[b, 7.2e-152]], $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], 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}\;b \leq -4.1 \cdot 10^{+19} \lor \neg \left(b \leq 7.2 \cdot 10^{-152}\right):\\
\;\;\;\;\frac{x}{x + y \cdot e^{2 \cdot \left(b \cdot \left(\frac{0.6666666666666666}{t} - \left(a + 0.8333333333333334\right)\right)\right)}}\\
\mathbf{else}:\\
\;\;\;\;\frac{x}{x + y \cdot e^{2 \cdot \left(c \cdot \left(a + 0.8333333333333334\right)\right)}}\\
\end{array}
\end{array}
if b < -4.1e19 or 7.2e-152 < b Initial program 94.2%
Taylor expanded in b around inf 81.5%
associate-*r/81.5%
metadata-eval81.5%
+-commutative81.5%
Simplified81.5%
if -4.1e19 < b < 7.2e-152Initial program 95.0%
Taylor expanded in c around inf 80.6%
+-commutative80.6%
associate-*r/80.6%
metadata-eval80.6%
Simplified80.6%
Taylor expanded in t around inf 67.1%
Final simplification75.9%
(FPCore (x y z t a b c)
:precision binary64
(if (or (<= c -8.6e+114) (not (<= c 2.4e+112)))
(/
x
(+
x
(*
y
(exp
(* 2.0 (* c (- (+ a 0.8333333333333334) (/ 0.6666666666666666 t))))))))
(/
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 tmp;
if ((c <= -8.6e+114) || !(c <= 2.4e+112)) {
tmp = x / (x + (y * exp((2.0 * (c * ((a + 0.8333333333333334) - (0.6666666666666666 / t)))))));
} else {
tmp = x / (x + (y * exp((2.0 * (b * ((0.6666666666666666 / t) - (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 ((c <= (-8.6d+114)) .or. (.not. (c <= 2.4d+112))) then
tmp = x / (x + (y * exp((2.0d0 * (c * ((a + 0.8333333333333334d0) - (0.6666666666666666d0 / t)))))))
else
tmp = x / (x + (y * exp((2.0d0 * (b * ((0.6666666666666666d0 / t) - (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 ((c <= -8.6e+114) || !(c <= 2.4e+112)) {
tmp = x / (x + (y * Math.exp((2.0 * (c * ((a + 0.8333333333333334) - (0.6666666666666666 / t)))))));
} 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): tmp = 0 if (c <= -8.6e+114) or not (c <= 2.4e+112): tmp = x / (x + (y * math.exp((2.0 * (c * ((a + 0.8333333333333334) - (0.6666666666666666 / t))))))) 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) tmp = 0.0 if ((c <= -8.6e+114) || !(c <= 2.4e+112)) tmp = Float64(x / Float64(x + Float64(y * exp(Float64(2.0 * Float64(c * Float64(Float64(a + 0.8333333333333334) - Float64(0.6666666666666666 / t)))))))); 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) tmp = 0.0; if ((c <= -8.6e+114) || ~((c <= 2.4e+112))) tmp = x / (x + (y * exp((2.0 * (c * ((a + 0.8333333333333334) - (0.6666666666666666 / t))))))); 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_] := If[Or[LessEqual[c, -8.6e+114], N[Not[LessEqual[c, 2.4e+112]], $MachinePrecision]], 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], 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}
\mathbf{if}\;c \leq -8.6 \cdot 10^{+114} \lor \neg \left(c \leq 2.4 \cdot 10^{+112}\right):\\
\;\;\;\;\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}:\\
\;\;\;\;\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 c < -8.6000000000000001e114 or 2.4e112 < c Initial program 90.5%
Taylor expanded in c around inf 91.9%
+-commutative91.9%
associate-*r/91.9%
metadata-eval91.9%
Simplified91.9%
if -8.6000000000000001e114 < c < 2.4e112Initial program 96.5%
Taylor expanded in b around inf 77.5%
associate-*r/77.5%
metadata-eval77.5%
+-commutative77.5%
Simplified77.5%
Final simplification82.2%
(FPCore (x y z t a b c)
:precision binary64
(if (<= t -8e-222)
(/ x (+ x (* y (exp (* b -1.6666666666666667)))))
(if (<= t 3.4e-126)
(/
x
(-
x
(*
y
(+
-1.0
(*
2.0
(*
a
(-
b
(* b (/ (+ (/ 0.6666666666666666 t) -0.8333333333333334) a)))))))))
(if (<= t 9.4e-11)
1.0
(/ x (+ x (* y (exp (* c 1.6666666666666667)))))))))
double code(double x, double y, double z, double t, double a, double b, double c) {
double tmp;
if (t <= -8e-222) {
tmp = x / (x + (y * exp((b * -1.6666666666666667))));
} else if (t <= 3.4e-126) {
tmp = x / (x - (y * (-1.0 + (2.0 * (a * (b - (b * (((0.6666666666666666 / t) + -0.8333333333333334) / a))))))));
} else if (t <= 9.4e-11) {
tmp = 1.0;
} else {
tmp = x / (x + (y * exp((c * 1.6666666666666667))));
}
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 <= (-8d-222)) then
tmp = x / (x + (y * exp((b * (-1.6666666666666667d0)))))
else if (t <= 3.4d-126) then
tmp = x / (x - (y * ((-1.0d0) + (2.0d0 * (a * (b - (b * (((0.6666666666666666d0 / t) + (-0.8333333333333334d0)) / a))))))))
else if (t <= 9.4d-11) then
tmp = 1.0d0
else
tmp = x / (x + (y * exp((c * 1.6666666666666667d0))))
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 <= -8e-222) {
tmp = x / (x + (y * Math.exp((b * -1.6666666666666667))));
} else if (t <= 3.4e-126) {
tmp = x / (x - (y * (-1.0 + (2.0 * (a * (b - (b * (((0.6666666666666666 / t) + -0.8333333333333334) / a))))))));
} else if (t <= 9.4e-11) {
tmp = 1.0;
} else {
tmp = x / (x + (y * Math.exp((c * 1.6666666666666667))));
}
return tmp;
}
def code(x, y, z, t, a, b, c): tmp = 0 if t <= -8e-222: tmp = x / (x + (y * math.exp((b * -1.6666666666666667)))) elif t <= 3.4e-126: tmp = x / (x - (y * (-1.0 + (2.0 * (a * (b - (b * (((0.6666666666666666 / t) + -0.8333333333333334) / a)))))))) elif t <= 9.4e-11: tmp = 1.0 else: tmp = x / (x + (y * math.exp((c * 1.6666666666666667)))) return tmp
function code(x, y, z, t, a, b, c) tmp = 0.0 if (t <= -8e-222) tmp = Float64(x / Float64(x + Float64(y * exp(Float64(b * -1.6666666666666667))))); elseif (t <= 3.4e-126) tmp = Float64(x / Float64(x - Float64(y * Float64(-1.0 + Float64(2.0 * Float64(a * Float64(b - Float64(b * Float64(Float64(Float64(0.6666666666666666 / t) + -0.8333333333333334) / a))))))))); elseif (t <= 9.4e-11) tmp = 1.0; else tmp = Float64(x / Float64(x + Float64(y * exp(Float64(c * 1.6666666666666667))))); end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c) tmp = 0.0; if (t <= -8e-222) tmp = x / (x + (y * exp((b * -1.6666666666666667)))); elseif (t <= 3.4e-126) tmp = x / (x - (y * (-1.0 + (2.0 * (a * (b - (b * (((0.6666666666666666 / t) + -0.8333333333333334) / a)))))))); elseif (t <= 9.4e-11) tmp = 1.0; else tmp = x / (x + (y * exp((c * 1.6666666666666667)))); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_] := If[LessEqual[t, -8e-222], N[(x / N[(x + N[(y * N[Exp[N[(b * -1.6666666666666667), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t, 3.4e-126], N[(x / N[(x - N[(y * N[(-1.0 + N[(2.0 * N[(a * N[(b - N[(b * N[(N[(N[(0.6666666666666666 / t), $MachinePrecision] + -0.8333333333333334), $MachinePrecision] / a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t, 9.4e-11], 1.0, N[(x / N[(x + N[(y * N[Exp[N[(c * 1.6666666666666667), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;t \leq -8 \cdot 10^{-222}:\\
\;\;\;\;\frac{x}{x + y \cdot e^{b \cdot -1.6666666666666667}}\\
\mathbf{elif}\;t \leq 3.4 \cdot 10^{-126}:\\
\;\;\;\;\frac{x}{x - y \cdot \left(-1 + 2 \cdot \left(a \cdot \left(b - b \cdot \frac{\frac{0.6666666666666666}{t} + -0.8333333333333334}{a}\right)\right)\right)}\\
\mathbf{elif}\;t \leq 9.4 \cdot 10^{-11}:\\
\;\;\;\;1\\
\mathbf{else}:\\
\;\;\;\;\frac{x}{x + y \cdot e^{c \cdot 1.6666666666666667}}\\
\end{array}
\end{array}
if t < -8.00000000000000038e-222Initial program 95.6%
Taylor expanded in b around inf 68.4%
associate-*r/68.4%
metadata-eval68.4%
+-commutative68.4%
Simplified68.4%
Taylor expanded in t around inf 72.6%
mul-1-neg72.6%
+-commutative72.6%
distribute-rgt-neg-in72.6%
+-commutative72.6%
mul-1-neg72.6%
distribute-lft-in72.6%
metadata-eval72.6%
mul-1-neg72.6%
unsub-neg72.6%
Simplified72.6%
Taylor expanded in a around 0 74.9%
if -8.00000000000000038e-222 < t < 3.4e-126Initial program 86.0%
Taylor expanded in b around inf 72.8%
associate-*r/72.8%
metadata-eval72.8%
+-commutative72.8%
Simplified72.8%
Taylor expanded in b around 0 47.6%
Taylor expanded in a around inf 54.4%
+-commutative54.4%
mul-1-neg54.4%
unsub-neg54.4%
associate-/l*57.6%
sub-neg57.6%
associate-*r/57.6%
metadata-eval57.6%
metadata-eval57.6%
Simplified57.6%
if 3.4e-126 < t < 9.39999999999999985e-11Initial program 97.2%
Taylor expanded in c around inf 73.1%
+-commutative73.1%
associate-*r/73.1%
metadata-eval73.1%
Simplified73.1%
Taylor expanded in t around inf 43.5%
Taylor expanded in x around inf 70.4%
if 9.39999999999999985e-11 < t Initial program 97.4%
Taylor expanded in c around inf 71.0%
+-commutative71.0%
associate-*r/71.0%
metadata-eval71.0%
Simplified71.0%
Taylor expanded in t around inf 72.6%
Taylor expanded in a around 0 64.3%
Final simplification65.6%
(FPCore (x y z t a b c)
:precision binary64
(if (<= y -8.5e+40)
(/
x
(+
x
(*
y
(+
(*
2.0
(/
(- (* b 0.6666666666666666) (* (+ a 0.8333333333333334) (* t b)))
t))
1.0))))
(if (<= y -2.7e-12)
1.0
(if (<= y -1.2e-271)
(/
x
(+
x
(*
y
(+
(*
2.0
(*
a
(-
(* b (/ (+ (/ 0.6666666666666666 t) -0.8333333333333334) a))
b)))
1.0))))
(if (<= y 2e-18)
1.0
(if (<= y 2.85e+34)
(/
x
(+
x
(*
b
(*
a
(+
(* y -2.0)
(+
(*
2.0
(/
(*
y
(- (* 0.6666666666666666 (/ 1.0 t)) 0.8333333333333334))
a))
(/ y (* a b))))))))
1.0))))))
double code(double x, double y, double z, double t, double a, double b, double c) {
double tmp;
if (y <= -8.5e+40) {
tmp = x / (x + (y * ((2.0 * (((b * 0.6666666666666666) - ((a + 0.8333333333333334) * (t * b))) / t)) + 1.0)));
} else if (y <= -2.7e-12) {
tmp = 1.0;
} else if (y <= -1.2e-271) {
tmp = x / (x + (y * ((2.0 * (a * ((b * (((0.6666666666666666 / t) + -0.8333333333333334) / a)) - b))) + 1.0)));
} else if (y <= 2e-18) {
tmp = 1.0;
} else if (y <= 2.85e+34) {
tmp = x / (x + (b * (a * ((y * -2.0) + ((2.0 * ((y * ((0.6666666666666666 * (1.0 / t)) - 0.8333333333333334)) / a)) + (y / (a * 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 (y <= (-8.5d+40)) then
tmp = x / (x + (y * ((2.0d0 * (((b * 0.6666666666666666d0) - ((a + 0.8333333333333334d0) * (t * b))) / t)) + 1.0d0)))
else if (y <= (-2.7d-12)) then
tmp = 1.0d0
else if (y <= (-1.2d-271)) then
tmp = x / (x + (y * ((2.0d0 * (a * ((b * (((0.6666666666666666d0 / t) + (-0.8333333333333334d0)) / a)) - b))) + 1.0d0)))
else if (y <= 2d-18) then
tmp = 1.0d0
else if (y <= 2.85d+34) then
tmp = x / (x + (b * (a * ((y * (-2.0d0)) + ((2.0d0 * ((y * ((0.6666666666666666d0 * (1.0d0 / t)) - 0.8333333333333334d0)) / a)) + (y / (a * 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 (y <= -8.5e+40) {
tmp = x / (x + (y * ((2.0 * (((b * 0.6666666666666666) - ((a + 0.8333333333333334) * (t * b))) / t)) + 1.0)));
} else if (y <= -2.7e-12) {
tmp = 1.0;
} else if (y <= -1.2e-271) {
tmp = x / (x + (y * ((2.0 * (a * ((b * (((0.6666666666666666 / t) + -0.8333333333333334) / a)) - b))) + 1.0)));
} else if (y <= 2e-18) {
tmp = 1.0;
} else if (y <= 2.85e+34) {
tmp = x / (x + (b * (a * ((y * -2.0) + ((2.0 * ((y * ((0.6666666666666666 * (1.0 / t)) - 0.8333333333333334)) / a)) + (y / (a * b)))))));
} else {
tmp = 1.0;
}
return tmp;
}
def code(x, y, z, t, a, b, c): tmp = 0 if y <= -8.5e+40: tmp = x / (x + (y * ((2.0 * (((b * 0.6666666666666666) - ((a + 0.8333333333333334) * (t * b))) / t)) + 1.0))) elif y <= -2.7e-12: tmp = 1.0 elif y <= -1.2e-271: tmp = x / (x + (y * ((2.0 * (a * ((b * (((0.6666666666666666 / t) + -0.8333333333333334) / a)) - b))) + 1.0))) elif y <= 2e-18: tmp = 1.0 elif y <= 2.85e+34: tmp = x / (x + (b * (a * ((y * -2.0) + ((2.0 * ((y * ((0.6666666666666666 * (1.0 / t)) - 0.8333333333333334)) / a)) + (y / (a * b))))))) else: tmp = 1.0 return tmp
function code(x, y, z, t, a, b, c) tmp = 0.0 if (y <= -8.5e+40) tmp = Float64(x / Float64(x + Float64(y * Float64(Float64(2.0 * Float64(Float64(Float64(b * 0.6666666666666666) - Float64(Float64(a + 0.8333333333333334) * Float64(t * b))) / t)) + 1.0)))); elseif (y <= -2.7e-12) tmp = 1.0; elseif (y <= -1.2e-271) tmp = Float64(x / Float64(x + Float64(y * Float64(Float64(2.0 * Float64(a * Float64(Float64(b * Float64(Float64(Float64(0.6666666666666666 / t) + -0.8333333333333334) / a)) - b))) + 1.0)))); elseif (y <= 2e-18) tmp = 1.0; elseif (y <= 2.85e+34) tmp = Float64(x / Float64(x + Float64(b * Float64(a * Float64(Float64(y * -2.0) + Float64(Float64(2.0 * Float64(Float64(y * Float64(Float64(0.6666666666666666 * Float64(1.0 / t)) - 0.8333333333333334)) / a)) + Float64(y / Float64(a * b)))))))); else tmp = 1.0; end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c) tmp = 0.0; if (y <= -8.5e+40) tmp = x / (x + (y * ((2.0 * (((b * 0.6666666666666666) - ((a + 0.8333333333333334) * (t * b))) / t)) + 1.0))); elseif (y <= -2.7e-12) tmp = 1.0; elseif (y <= -1.2e-271) tmp = x / (x + (y * ((2.0 * (a * ((b * (((0.6666666666666666 / t) + -0.8333333333333334) / a)) - b))) + 1.0))); elseif (y <= 2e-18) tmp = 1.0; elseif (y <= 2.85e+34) tmp = x / (x + (b * (a * ((y * -2.0) + ((2.0 * ((y * ((0.6666666666666666 * (1.0 / t)) - 0.8333333333333334)) / a)) + (y / (a * b))))))); else tmp = 1.0; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_] := If[LessEqual[y, -8.5e+40], N[(x / N[(x + N[(y * N[(N[(2.0 * N[(N[(N[(b * 0.6666666666666666), $MachinePrecision] - N[(N[(a + 0.8333333333333334), $MachinePrecision] * N[(t * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / t), $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, -2.7e-12], 1.0, If[LessEqual[y, -1.2e-271], N[(x / N[(x + N[(y * N[(N[(2.0 * N[(a * N[(N[(b * N[(N[(N[(0.6666666666666666 / t), $MachinePrecision] + -0.8333333333333334), $MachinePrecision] / a), $MachinePrecision]), $MachinePrecision] - b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 2e-18], 1.0, If[LessEqual[y, 2.85e+34], N[(x / N[(x + N[(b * N[(a * N[(N[(y * -2.0), $MachinePrecision] + N[(N[(2.0 * N[(N[(y * N[(N[(0.6666666666666666 * N[(1.0 / t), $MachinePrecision]), $MachinePrecision] - 0.8333333333333334), $MachinePrecision]), $MachinePrecision] / a), $MachinePrecision]), $MachinePrecision] + N[(y / N[(a * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 1.0]]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -8.5 \cdot 10^{+40}:\\
\;\;\;\;\frac{x}{x + y \cdot \left(2 \cdot \frac{b \cdot 0.6666666666666666 - \left(a + 0.8333333333333334\right) \cdot \left(t \cdot b\right)}{t} + 1\right)}\\
\mathbf{elif}\;y \leq -2.7 \cdot 10^{-12}:\\
\;\;\;\;1\\
\mathbf{elif}\;y \leq -1.2 \cdot 10^{-271}:\\
\;\;\;\;\frac{x}{x + y \cdot \left(2 \cdot \left(a \cdot \left(b \cdot \frac{\frac{0.6666666666666666}{t} + -0.8333333333333334}{a} - b\right)\right) + 1\right)}\\
\mathbf{elif}\;y \leq 2 \cdot 10^{-18}:\\
\;\;\;\;1\\
\mathbf{elif}\;y \leq 2.85 \cdot 10^{+34}:\\
\;\;\;\;\frac{x}{x + b \cdot \left(a \cdot \left(y \cdot -2 + \left(2 \cdot \frac{y \cdot \left(0.6666666666666666 \cdot \frac{1}{t} - 0.8333333333333334\right)}{a} + \frac{y}{a \cdot b}\right)\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;1\\
\end{array}
\end{array}
if y < -8.49999999999999996e40Initial program 93.7%
Taylor expanded in b around inf 75.7%
associate-*r/75.7%
metadata-eval75.7%
+-commutative75.7%
Simplified75.7%
Taylor expanded in b around 0 67.8%
Taylor expanded in t around 0 71.7%
+-commutative71.7%
mul-1-neg71.7%
unsub-neg71.7%
*-commutative71.7%
associate-*r*71.7%
Simplified71.7%
if -8.49999999999999996e40 < y < -2.6999999999999998e-12 or -1.2000000000000001e-271 < y < 2.0000000000000001e-18 or 2.84999999999999987e34 < y Initial program 96.4%
Taylor expanded in c around inf 68.3%
+-commutative68.3%
associate-*r/68.3%
metadata-eval68.3%
Simplified68.3%
Taylor expanded in t around inf 55.6%
Taylor expanded in x around inf 62.1%
if -2.6999999999999998e-12 < y < -1.2000000000000001e-271Initial program 91.4%
Taylor expanded in b around inf 71.7%
associate-*r/71.7%
metadata-eval71.7%
+-commutative71.7%
Simplified71.7%
Taylor expanded in b around 0 45.4%
Taylor expanded in a around inf 58.6%
+-commutative58.6%
mul-1-neg58.6%
unsub-neg58.6%
associate-/l*61.7%
sub-neg61.7%
associate-*r/61.7%
metadata-eval61.7%
metadata-eval61.7%
Simplified61.7%
if 2.0000000000000001e-18 < y < 2.84999999999999987e34Initial program 91.5%
Taylor expanded in b around inf 59.6%
associate-*r/59.6%
metadata-eval59.6%
+-commutative59.6%
Simplified59.6%
Taylor expanded in b around 0 46.7%
Taylor expanded in b around inf 54.8%
Taylor expanded in a around inf 78.9%
Final simplification64.6%
(FPCore (x y z t a b c)
:precision binary64
(if (<= (- b c) -2e+236)
(/
x
(+
x
(*
y
(+
(* 2.0 (* c (+ (+ a (/ -0.6666666666666666 t)) 0.8333333333333334)))
1.0))))
(if (<= (- b c) -2e+77)
(/
x
(+
x
(*
y
(+
(*
2.0
(/
(- (* b 0.6666666666666666) (* (+ a 0.8333333333333334) (* t b)))
t))
1.0))))
(if (<= (- b c) -1e-5)
1.0
(if (<= (- b c) 40.0)
(/ x (+ x (* y (+ (* -2.0 (* b (+ a 0.8333333333333334))) 1.0))))
1.0)))))
double code(double x, double y, double z, double t, double a, double b, double c) {
double tmp;
if ((b - c) <= -2e+236) {
tmp = x / (x + (y * ((2.0 * (c * ((a + (-0.6666666666666666 / t)) + 0.8333333333333334))) + 1.0)));
} else if ((b - c) <= -2e+77) {
tmp = x / (x + (y * ((2.0 * (((b * 0.6666666666666666) - ((a + 0.8333333333333334) * (t * b))) / t)) + 1.0)));
} else if ((b - c) <= -1e-5) {
tmp = 1.0;
} else if ((b - c) <= 40.0) {
tmp = x / (x + (y * ((-2.0 * (b * (a + 0.8333333333333334))) + 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) :: tmp
if ((b - c) <= (-2d+236)) then
tmp = x / (x + (y * ((2.0d0 * (c * ((a + ((-0.6666666666666666d0) / t)) + 0.8333333333333334d0))) + 1.0d0)))
else if ((b - c) <= (-2d+77)) then
tmp = x / (x + (y * ((2.0d0 * (((b * 0.6666666666666666d0) - ((a + 0.8333333333333334d0) * (t * b))) / t)) + 1.0d0)))
else if ((b - c) <= (-1d-5)) then
tmp = 1.0d0
else if ((b - c) <= 40.0d0) then
tmp = x / (x + (y * (((-2.0d0) * (b * (a + 0.8333333333333334d0))) + 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 tmp;
if ((b - c) <= -2e+236) {
tmp = x / (x + (y * ((2.0 * (c * ((a + (-0.6666666666666666 / t)) + 0.8333333333333334))) + 1.0)));
} else if ((b - c) <= -2e+77) {
tmp = x / (x + (y * ((2.0 * (((b * 0.6666666666666666) - ((a + 0.8333333333333334) * (t * b))) / t)) + 1.0)));
} else if ((b - c) <= -1e-5) {
tmp = 1.0;
} else if ((b - c) <= 40.0) {
tmp = x / (x + (y * ((-2.0 * (b * (a + 0.8333333333333334))) + 1.0)));
} else {
tmp = 1.0;
}
return tmp;
}
def code(x, y, z, t, a, b, c): tmp = 0 if (b - c) <= -2e+236: tmp = x / (x + (y * ((2.0 * (c * ((a + (-0.6666666666666666 / t)) + 0.8333333333333334))) + 1.0))) elif (b - c) <= -2e+77: tmp = x / (x + (y * ((2.0 * (((b * 0.6666666666666666) - ((a + 0.8333333333333334) * (t * b))) / t)) + 1.0))) elif (b - c) <= -1e-5: tmp = 1.0 elif (b - c) <= 40.0: tmp = x / (x + (y * ((-2.0 * (b * (a + 0.8333333333333334))) + 1.0))) else: tmp = 1.0 return tmp
function code(x, y, z, t, a, b, c) tmp = 0.0 if (Float64(b - c) <= -2e+236) tmp = Float64(x / Float64(x + Float64(y * Float64(Float64(2.0 * Float64(c * Float64(Float64(a + Float64(-0.6666666666666666 / t)) + 0.8333333333333334))) + 1.0)))); elseif (Float64(b - c) <= -2e+77) tmp = Float64(x / Float64(x + Float64(y * Float64(Float64(2.0 * Float64(Float64(Float64(b * 0.6666666666666666) - Float64(Float64(a + 0.8333333333333334) * Float64(t * b))) / t)) + 1.0)))); elseif (Float64(b - c) <= -1e-5) tmp = 1.0; elseif (Float64(b - c) <= 40.0) tmp = Float64(x / Float64(x + Float64(y * Float64(Float64(-2.0 * Float64(b * Float64(a + 0.8333333333333334))) + 1.0)))); else tmp = 1.0; end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c) tmp = 0.0; if ((b - c) <= -2e+236) tmp = x / (x + (y * ((2.0 * (c * ((a + (-0.6666666666666666 / t)) + 0.8333333333333334))) + 1.0))); elseif ((b - c) <= -2e+77) tmp = x / (x + (y * ((2.0 * (((b * 0.6666666666666666) - ((a + 0.8333333333333334) * (t * b))) / t)) + 1.0))); elseif ((b - c) <= -1e-5) tmp = 1.0; elseif ((b - c) <= 40.0) tmp = x / (x + (y * ((-2.0 * (b * (a + 0.8333333333333334))) + 1.0))); else tmp = 1.0; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_] := If[LessEqual[N[(b - c), $MachinePrecision], -2e+236], N[(x / N[(x + N[(y * N[(N[(2.0 * N[(c * N[(N[(a + N[(-0.6666666666666666 / t), $MachinePrecision]), $MachinePrecision] + 0.8333333333333334), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[N[(b - c), $MachinePrecision], -2e+77], N[(x / N[(x + N[(y * N[(N[(2.0 * N[(N[(N[(b * 0.6666666666666666), $MachinePrecision] - N[(N[(a + 0.8333333333333334), $MachinePrecision] * N[(t * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / t), $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[N[(b - c), $MachinePrecision], -1e-5], 1.0, If[LessEqual[N[(b - c), $MachinePrecision], 40.0], N[(x / N[(x + N[(y * N[(N[(-2.0 * N[(b * N[(a + 0.8333333333333334), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 1.0]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;b - c \leq -2 \cdot 10^{+236}:\\
\;\;\;\;\frac{x}{x + y \cdot \left(2 \cdot \left(c \cdot \left(\left(a + \frac{-0.6666666666666666}{t}\right) + 0.8333333333333334\right)\right) + 1\right)}\\
\mathbf{elif}\;b - c \leq -2 \cdot 10^{+77}:\\
\;\;\;\;\frac{x}{x + y \cdot \left(2 \cdot \frac{b \cdot 0.6666666666666666 - \left(a + 0.8333333333333334\right) \cdot \left(t \cdot b\right)}{t} + 1\right)}\\
\mathbf{elif}\;b - c \leq -1 \cdot 10^{-5}:\\
\;\;\;\;1\\
\mathbf{elif}\;b - c \leq 40:\\
\;\;\;\;\frac{x}{x + y \cdot \left(-2 \cdot \left(b \cdot \left(a + 0.8333333333333334\right)\right) + 1\right)}\\
\mathbf{else}:\\
\;\;\;\;1\\
\end{array}
\end{array}
if (-.f64 b c) < -2.00000000000000011e236Initial program 85.7%
Taylor expanded in c around inf 75.2%
+-commutative75.2%
associate-*r/75.2%
metadata-eval75.2%
Simplified75.2%
Taylor expanded in c around 0 58.7%
*-lft-identity58.7%
metadata-eval58.7%
cancel-sign-sub-inv58.7%
associate--r+58.7%
associate--r+58.7%
cancel-sign-sub-inv58.7%
metadata-eval58.7%
*-lft-identity58.7%
associate--l+58.7%
cancel-sign-sub-inv58.7%
metadata-eval58.7%
associate-*r/58.7%
metadata-eval58.7%
Simplified58.7%
if -2.00000000000000011e236 < (-.f64 b c) < -1.99999999999999997e77Initial program 86.6%
Taylor expanded in b around inf 68.4%
associate-*r/68.4%
metadata-eval68.4%
+-commutative68.4%
Simplified68.4%
Taylor expanded in b around 0 51.7%
Taylor expanded in t around 0 57.2%
+-commutative57.2%
mul-1-neg57.2%
unsub-neg57.2%
*-commutative57.2%
associate-*r*57.3%
Simplified57.3%
if -1.99999999999999997e77 < (-.f64 b c) < -1.00000000000000008e-5 or 40 < (-.f64 b c) Initial program 98.0%
Taylor expanded in c around inf 67.7%
+-commutative67.7%
associate-*r/67.7%
metadata-eval67.7%
Simplified67.7%
Taylor expanded in t around inf 57.3%
Taylor expanded in x around inf 64.9%
if -1.00000000000000008e-5 < (-.f64 b c) < 40Initial program 99.9%
Taylor expanded in b around inf 68.3%
associate-*r/68.3%
metadata-eval68.3%
+-commutative68.3%
Simplified68.3%
Taylor expanded in b around 0 63.1%
Taylor expanded in t around inf 66.1%
*-commutative66.1%
*-commutative66.1%
Simplified66.1%
Final simplification62.8%
(FPCore (x y z t a b c)
:precision binary64
(if (or (<= y -5.5e+41) (and (not (<= y -3.3e-12)) (<= y -1.5e-264)))
(/
x
(-
x
(*
y
(+
-1.0
(*
2.0
(*
a
(-
b
(* b (/ (+ (/ 0.6666666666666666 t) -0.8333333333333334) a)))))))))
1.0))
double code(double x, double y, double z, double t, double a, double b, double c) {
double tmp;
if ((y <= -5.5e+41) || (!(y <= -3.3e-12) && (y <= -1.5e-264))) {
tmp = x / (x - (y * (-1.0 + (2.0 * (a * (b - (b * (((0.6666666666666666 / t) + -0.8333333333333334) / a))))))));
} 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 ((y <= (-5.5d+41)) .or. (.not. (y <= (-3.3d-12))) .and. (y <= (-1.5d-264))) then
tmp = x / (x - (y * ((-1.0d0) + (2.0d0 * (a * (b - (b * (((0.6666666666666666d0 / t) + (-0.8333333333333334d0)) / a))))))))
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 ((y <= -5.5e+41) || (!(y <= -3.3e-12) && (y <= -1.5e-264))) {
tmp = x / (x - (y * (-1.0 + (2.0 * (a * (b - (b * (((0.6666666666666666 / t) + -0.8333333333333334) / a))))))));
} else {
tmp = 1.0;
}
return tmp;
}
def code(x, y, z, t, a, b, c): tmp = 0 if (y <= -5.5e+41) or (not (y <= -3.3e-12) and (y <= -1.5e-264)): tmp = x / (x - (y * (-1.0 + (2.0 * (a * (b - (b * (((0.6666666666666666 / t) + -0.8333333333333334) / a)))))))) else: tmp = 1.0 return tmp
function code(x, y, z, t, a, b, c) tmp = 0.0 if ((y <= -5.5e+41) || (!(y <= -3.3e-12) && (y <= -1.5e-264))) tmp = Float64(x / Float64(x - Float64(y * Float64(-1.0 + Float64(2.0 * Float64(a * Float64(b - Float64(b * Float64(Float64(Float64(0.6666666666666666 / t) + -0.8333333333333334) / a))))))))); else tmp = 1.0; end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c) tmp = 0.0; if ((y <= -5.5e+41) || (~((y <= -3.3e-12)) && (y <= -1.5e-264))) tmp = x / (x - (y * (-1.0 + (2.0 * (a * (b - (b * (((0.6666666666666666 / t) + -0.8333333333333334) / a)))))))); else tmp = 1.0; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_] := If[Or[LessEqual[y, -5.5e+41], And[N[Not[LessEqual[y, -3.3e-12]], $MachinePrecision], LessEqual[y, -1.5e-264]]], N[(x / N[(x - N[(y * N[(-1.0 + N[(2.0 * N[(a * N[(b - N[(b * N[(N[(N[(0.6666666666666666 / t), $MachinePrecision] + -0.8333333333333334), $MachinePrecision] / a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 1.0]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -5.5 \cdot 10^{+41} \lor \neg \left(y \leq -3.3 \cdot 10^{-12}\right) \land y \leq -1.5 \cdot 10^{-264}:\\
\;\;\;\;\frac{x}{x - y \cdot \left(-1 + 2 \cdot \left(a \cdot \left(b - b \cdot \frac{\frac{0.6666666666666666}{t} + -0.8333333333333334}{a}\right)\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;1\\
\end{array}
\end{array}
if y < -5.5000000000000003e41 or -3.3000000000000001e-12 < y < -1.5e-264Initial program 92.4%
Taylor expanded in b around inf 73.3%
associate-*r/73.3%
metadata-eval73.3%
+-commutative73.3%
Simplified73.3%
Taylor expanded in b around 0 56.0%
Taylor expanded in a around inf 64.2%
+-commutative64.2%
mul-1-neg64.2%
unsub-neg64.2%
associate-/l*65.9%
sub-neg65.9%
associate-*r/65.9%
metadata-eval65.9%
metadata-eval65.9%
Simplified65.9%
if -5.5000000000000003e41 < y < -3.3000000000000001e-12 or -1.5e-264 < y Initial program 96.0%
Taylor expanded in c around inf 67.3%
+-commutative67.3%
associate-*r/67.3%
metadata-eval67.3%
Simplified67.3%
Taylor expanded in t around inf 55.1%
Taylor expanded in x around inf 58.9%
Final simplification61.8%
(FPCore (x y z t a b c)
:precision binary64
(if (<= y -1.65e+41)
(/ x (- x (* y (- -1.0 (* -2.0 (* b (+ a 0.8333333333333334)))))))
(if (<= y 1.1e-25)
1.0
(if (<= y 2.8e+32)
(/
x
(+
x
(*
y
(+
(* 2.0 (* c (+ (+ a (/ -0.6666666666666666 t)) 0.8333333333333334)))
1.0))))
1.0))))
double code(double x, double y, double z, double t, double a, double b, double c) {
double tmp;
if (y <= -1.65e+41) {
tmp = x / (x - (y * (-1.0 - (-2.0 * (b * (a + 0.8333333333333334))))));
} else if (y <= 1.1e-25) {
tmp = 1.0;
} else if (y <= 2.8e+32) {
tmp = x / (x + (y * ((2.0 * (c * ((a + (-0.6666666666666666 / t)) + 0.8333333333333334))) + 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) :: tmp
if (y <= (-1.65d+41)) then
tmp = x / (x - (y * ((-1.0d0) - ((-2.0d0) * (b * (a + 0.8333333333333334d0))))))
else if (y <= 1.1d-25) then
tmp = 1.0d0
else if (y <= 2.8d+32) then
tmp = x / (x + (y * ((2.0d0 * (c * ((a + ((-0.6666666666666666d0) / t)) + 0.8333333333333334d0))) + 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 tmp;
if (y <= -1.65e+41) {
tmp = x / (x - (y * (-1.0 - (-2.0 * (b * (a + 0.8333333333333334))))));
} else if (y <= 1.1e-25) {
tmp = 1.0;
} else if (y <= 2.8e+32) {
tmp = x / (x + (y * ((2.0 * (c * ((a + (-0.6666666666666666 / t)) + 0.8333333333333334))) + 1.0)));
} else {
tmp = 1.0;
}
return tmp;
}
def code(x, y, z, t, a, b, c): tmp = 0 if y <= -1.65e+41: tmp = x / (x - (y * (-1.0 - (-2.0 * (b * (a + 0.8333333333333334)))))) elif y <= 1.1e-25: tmp = 1.0 elif y <= 2.8e+32: tmp = x / (x + (y * ((2.0 * (c * ((a + (-0.6666666666666666 / t)) + 0.8333333333333334))) + 1.0))) else: tmp = 1.0 return tmp
function code(x, y, z, t, a, b, c) tmp = 0.0 if (y <= -1.65e+41) tmp = Float64(x / Float64(x - Float64(y * Float64(-1.0 - Float64(-2.0 * Float64(b * Float64(a + 0.8333333333333334))))))); elseif (y <= 1.1e-25) tmp = 1.0; elseif (y <= 2.8e+32) tmp = Float64(x / Float64(x + Float64(y * Float64(Float64(2.0 * Float64(c * Float64(Float64(a + Float64(-0.6666666666666666 / t)) + 0.8333333333333334))) + 1.0)))); else tmp = 1.0; end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c) tmp = 0.0; if (y <= -1.65e+41) tmp = x / (x - (y * (-1.0 - (-2.0 * (b * (a + 0.8333333333333334)))))); elseif (y <= 1.1e-25) tmp = 1.0; elseif (y <= 2.8e+32) tmp = x / (x + (y * ((2.0 * (c * ((a + (-0.6666666666666666 / t)) + 0.8333333333333334))) + 1.0))); else tmp = 1.0; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_] := If[LessEqual[y, -1.65e+41], N[(x / N[(x - N[(y * N[(-1.0 - N[(-2.0 * N[(b * N[(a + 0.8333333333333334), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 1.1e-25], 1.0, If[LessEqual[y, 2.8e+32], N[(x / N[(x + N[(y * N[(N[(2.0 * N[(c * N[(N[(a + N[(-0.6666666666666666 / t), $MachinePrecision]), $MachinePrecision] + 0.8333333333333334), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 1.0]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -1.65 \cdot 10^{+41}:\\
\;\;\;\;\frac{x}{x - y \cdot \left(-1 - -2 \cdot \left(b \cdot \left(a + 0.8333333333333334\right)\right)\right)}\\
\mathbf{elif}\;y \leq 1.1 \cdot 10^{-25}:\\
\;\;\;\;1\\
\mathbf{elif}\;y \leq 2.8 \cdot 10^{+32}:\\
\;\;\;\;\frac{x}{x + y \cdot \left(2 \cdot \left(c \cdot \left(\left(a + \frac{-0.6666666666666666}{t}\right) + 0.8333333333333334\right)\right) + 1\right)}\\
\mathbf{else}:\\
\;\;\;\;1\\
\end{array}
\end{array}
if y < -1.65e41Initial program 93.7%
Taylor expanded in b around inf 75.7%
associate-*r/75.7%
metadata-eval75.7%
+-commutative75.7%
Simplified75.7%
Taylor expanded in b around 0 67.8%
Taylor expanded in t around inf 65.9%
*-commutative65.9%
*-commutative65.9%
Simplified65.9%
if -1.65e41 < y < 1.1000000000000001e-25 or 2.8e32 < y Initial program 95.4%
Taylor expanded in c around inf 68.4%
+-commutative68.4%
associate-*r/68.4%
metadata-eval68.4%
Simplified68.4%
Taylor expanded in t around inf 55.0%
Taylor expanded in x around inf 58.5%
if 1.1000000000000001e-25 < y < 2.8e32Initial program 85.6%
Taylor expanded in c around inf 60.3%
+-commutative60.3%
associate-*r/60.3%
metadata-eval60.3%
Simplified60.3%
Taylor expanded in c around 0 74.1%
*-lft-identity74.1%
metadata-eval74.1%
cancel-sign-sub-inv74.1%
associate--r+74.1%
associate--r+74.1%
cancel-sign-sub-inv74.1%
metadata-eval74.1%
*-lft-identity74.1%
associate--l+74.1%
cancel-sign-sub-inv74.1%
metadata-eval74.1%
associate-*r/74.1%
metadata-eval74.1%
Simplified74.1%
Final simplification60.7%
(FPCore (x y z t a b c)
:precision binary64
(if (<= y -1.55e+41)
(/
x
(+
x
(*
y
(+
(* (- (+ (/ 0.6666666666666666 t) -0.8333333333333334) a) (* b 2.0))
1.0))))
(if (<= y 1.75e-25)
1.0
(if (<= y 5.4e+33)
(/
x
(+
x
(*
y
(+
(* 2.0 (* c (+ (+ a (/ -0.6666666666666666 t)) 0.8333333333333334)))
1.0))))
1.0))))
double code(double x, double y, double z, double t, double a, double b, double c) {
double tmp;
if (y <= -1.55e+41) {
tmp = x / (x + (y * (((((0.6666666666666666 / t) + -0.8333333333333334) - a) * (b * 2.0)) + 1.0)));
} else if (y <= 1.75e-25) {
tmp = 1.0;
} else if (y <= 5.4e+33) {
tmp = x / (x + (y * ((2.0 * (c * ((a + (-0.6666666666666666 / t)) + 0.8333333333333334))) + 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) :: tmp
if (y <= (-1.55d+41)) then
tmp = x / (x + (y * (((((0.6666666666666666d0 / t) + (-0.8333333333333334d0)) - a) * (b * 2.0d0)) + 1.0d0)))
else if (y <= 1.75d-25) then
tmp = 1.0d0
else if (y <= 5.4d+33) then
tmp = x / (x + (y * ((2.0d0 * (c * ((a + ((-0.6666666666666666d0) / t)) + 0.8333333333333334d0))) + 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 tmp;
if (y <= -1.55e+41) {
tmp = x / (x + (y * (((((0.6666666666666666 / t) + -0.8333333333333334) - a) * (b * 2.0)) + 1.0)));
} else if (y <= 1.75e-25) {
tmp = 1.0;
} else if (y <= 5.4e+33) {
tmp = x / (x + (y * ((2.0 * (c * ((a + (-0.6666666666666666 / t)) + 0.8333333333333334))) + 1.0)));
} else {
tmp = 1.0;
}
return tmp;
}
def code(x, y, z, t, a, b, c): tmp = 0 if y <= -1.55e+41: tmp = x / (x + (y * (((((0.6666666666666666 / t) + -0.8333333333333334) - a) * (b * 2.0)) + 1.0))) elif y <= 1.75e-25: tmp = 1.0 elif y <= 5.4e+33: tmp = x / (x + (y * ((2.0 * (c * ((a + (-0.6666666666666666 / t)) + 0.8333333333333334))) + 1.0))) else: tmp = 1.0 return tmp
function code(x, y, z, t, a, b, c) tmp = 0.0 if (y <= -1.55e+41) tmp = Float64(x / Float64(x + Float64(y * Float64(Float64(Float64(Float64(Float64(0.6666666666666666 / t) + -0.8333333333333334) - a) * Float64(b * 2.0)) + 1.0)))); elseif (y <= 1.75e-25) tmp = 1.0; elseif (y <= 5.4e+33) tmp = Float64(x / Float64(x + Float64(y * Float64(Float64(2.0 * Float64(c * Float64(Float64(a + Float64(-0.6666666666666666 / t)) + 0.8333333333333334))) + 1.0)))); else tmp = 1.0; end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c) tmp = 0.0; if (y <= -1.55e+41) tmp = x / (x + (y * (((((0.6666666666666666 / t) + -0.8333333333333334) - a) * (b * 2.0)) + 1.0))); elseif (y <= 1.75e-25) tmp = 1.0; elseif (y <= 5.4e+33) tmp = x / (x + (y * ((2.0 * (c * ((a + (-0.6666666666666666 / t)) + 0.8333333333333334))) + 1.0))); else tmp = 1.0; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_] := If[LessEqual[y, -1.55e+41], N[(x / N[(x + N[(y * N[(N[(N[(N[(N[(0.6666666666666666 / t), $MachinePrecision] + -0.8333333333333334), $MachinePrecision] - a), $MachinePrecision] * N[(b * 2.0), $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 1.75e-25], 1.0, If[LessEqual[y, 5.4e+33], N[(x / N[(x + N[(y * N[(N[(2.0 * N[(c * N[(N[(a + N[(-0.6666666666666666 / t), $MachinePrecision]), $MachinePrecision] + 0.8333333333333334), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 1.0]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -1.55 \cdot 10^{+41}:\\
\;\;\;\;\frac{x}{x + y \cdot \left(\left(\left(\frac{0.6666666666666666}{t} + -0.8333333333333334\right) - a\right) \cdot \left(b \cdot 2\right) + 1\right)}\\
\mathbf{elif}\;y \leq 1.75 \cdot 10^{-25}:\\
\;\;\;\;1\\
\mathbf{elif}\;y \leq 5.4 \cdot 10^{+33}:\\
\;\;\;\;\frac{x}{x + y \cdot \left(2 \cdot \left(c \cdot \left(\left(a + \frac{-0.6666666666666666}{t}\right) + 0.8333333333333334\right)\right) + 1\right)}\\
\mathbf{else}:\\
\;\;\;\;1\\
\end{array}
\end{array}
if y < -1.55e41Initial program 93.7%
Taylor expanded in b around inf 75.7%
associate-*r/75.7%
metadata-eval75.7%
+-commutative75.7%
Simplified75.7%
Taylor expanded in b around 0 67.8%
associate-*r*67.8%
associate-*r/67.8%
metadata-eval67.8%
+-commutative67.8%
*-commutative67.8%
+-commutative67.8%
associate--r+67.8%
metadata-eval67.8%
associate-*r/67.8%
sub-neg67.8%
associate-*r/67.8%
metadata-eval67.8%
metadata-eval67.8%
*-commutative67.8%
Simplified67.8%
if -1.55e41 < y < 1.7500000000000001e-25 or 5.39999999999999982e33 < y Initial program 95.4%
Taylor expanded in c around inf 68.4%
+-commutative68.4%
associate-*r/68.4%
metadata-eval68.4%
Simplified68.4%
Taylor expanded in t around inf 55.0%
Taylor expanded in x around inf 58.5%
if 1.7500000000000001e-25 < y < 5.39999999999999982e33Initial program 85.6%
Taylor expanded in c around inf 60.3%
+-commutative60.3%
associate-*r/60.3%
metadata-eval60.3%
Simplified60.3%
Taylor expanded in c around 0 74.1%
*-lft-identity74.1%
metadata-eval74.1%
cancel-sign-sub-inv74.1%
associate--r+74.1%
associate--r+74.1%
cancel-sign-sub-inv74.1%
metadata-eval74.1%
*-lft-identity74.1%
associate--l+74.1%
cancel-sign-sub-inv74.1%
metadata-eval74.1%
associate-*r/74.1%
metadata-eval74.1%
Simplified74.1%
Final simplification61.1%
(FPCore (x y z t a b c)
:precision binary64
(if (<= y -2.8e+41)
(/ x (- x (* y (- -1.0 (* -2.0 (* b (+ a 0.8333333333333334)))))))
(if (<= y 1.6e-23)
1.0
(if (<= y 1e+36)
(/ x (+ x (+ y (* 2.0 (* c (* y (+ a 0.8333333333333334)))))))
1.0))))
double code(double x, double y, double z, double t, double a, double b, double c) {
double tmp;
if (y <= -2.8e+41) {
tmp = x / (x - (y * (-1.0 - (-2.0 * (b * (a + 0.8333333333333334))))));
} else if (y <= 1.6e-23) {
tmp = 1.0;
} else if (y <= 1e+36) {
tmp = x / (x + (y + (2.0 * (c * (y * (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) :: tmp
if (y <= (-2.8d+41)) then
tmp = x / (x - (y * ((-1.0d0) - ((-2.0d0) * (b * (a + 0.8333333333333334d0))))))
else if (y <= 1.6d-23) then
tmp = 1.0d0
else if (y <= 1d+36) then
tmp = x / (x + (y + (2.0d0 * (c * (y * (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 tmp;
if (y <= -2.8e+41) {
tmp = x / (x - (y * (-1.0 - (-2.0 * (b * (a + 0.8333333333333334))))));
} else if (y <= 1.6e-23) {
tmp = 1.0;
} else if (y <= 1e+36) {
tmp = x / (x + (y + (2.0 * (c * (y * (a + 0.8333333333333334))))));
} else {
tmp = 1.0;
}
return tmp;
}
def code(x, y, z, t, a, b, c): tmp = 0 if y <= -2.8e+41: tmp = x / (x - (y * (-1.0 - (-2.0 * (b * (a + 0.8333333333333334)))))) elif y <= 1.6e-23: tmp = 1.0 elif y <= 1e+36: tmp = x / (x + (y + (2.0 * (c * (y * (a + 0.8333333333333334)))))) else: tmp = 1.0 return tmp
function code(x, y, z, t, a, b, c) tmp = 0.0 if (y <= -2.8e+41) tmp = Float64(x / Float64(x - Float64(y * Float64(-1.0 - Float64(-2.0 * Float64(b * Float64(a + 0.8333333333333334))))))); elseif (y <= 1.6e-23) tmp = 1.0; elseif (y <= 1e+36) tmp = Float64(x / Float64(x + Float64(y + Float64(2.0 * Float64(c * Float64(y * Float64(a + 0.8333333333333334))))))); else tmp = 1.0; end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c) tmp = 0.0; if (y <= -2.8e+41) tmp = x / (x - (y * (-1.0 - (-2.0 * (b * (a + 0.8333333333333334)))))); elseif (y <= 1.6e-23) tmp = 1.0; elseif (y <= 1e+36) tmp = x / (x + (y + (2.0 * (c * (y * (a + 0.8333333333333334)))))); else tmp = 1.0; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_] := If[LessEqual[y, -2.8e+41], N[(x / N[(x - N[(y * N[(-1.0 - N[(-2.0 * N[(b * N[(a + 0.8333333333333334), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 1.6e-23], 1.0, If[LessEqual[y, 1e+36], N[(x / N[(x + N[(y + N[(2.0 * N[(c * N[(y * N[(a + 0.8333333333333334), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 1.0]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -2.8 \cdot 10^{+41}:\\
\;\;\;\;\frac{x}{x - y \cdot \left(-1 - -2 \cdot \left(b \cdot \left(a + 0.8333333333333334\right)\right)\right)}\\
\mathbf{elif}\;y \leq 1.6 \cdot 10^{-23}:\\
\;\;\;\;1\\
\mathbf{elif}\;y \leq 10^{+36}:\\
\;\;\;\;\frac{x}{x + \left(y + 2 \cdot \left(c \cdot \left(y \cdot \left(a + 0.8333333333333334\right)\right)\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;1\\
\end{array}
\end{array}
if y < -2.7999999999999999e41Initial program 93.7%
Taylor expanded in b around inf 75.7%
associate-*r/75.7%
metadata-eval75.7%
+-commutative75.7%
Simplified75.7%
Taylor expanded in b around 0 67.8%
Taylor expanded in t around inf 65.9%
*-commutative65.9%
*-commutative65.9%
Simplified65.9%
if -2.7999999999999999e41 < y < 1.59999999999999988e-23 or 1.00000000000000004e36 < y Initial program 95.4%
Taylor expanded in c around inf 68.4%
+-commutative68.4%
associate-*r/68.4%
metadata-eval68.4%
Simplified68.4%
Taylor expanded in t around inf 55.0%
Taylor expanded in x around inf 58.5%
if 1.59999999999999988e-23 < y < 1.00000000000000004e36Initial program 85.6%
Taylor expanded in c around inf 60.3%
+-commutative60.3%
associate-*r/60.3%
metadata-eval60.3%
Simplified60.3%
Taylor expanded in t around inf 53.3%
Taylor expanded in c around 0 67.2%
Final simplification60.3%
(FPCore (x y z t a b c) :precision binary64 (if (<= y -7.2e+216) (/ (* x -0.5) (* a (* b y))) (if (<= y -4.7e+114) 1.0 (if (<= y -6e+41) (/ 1.0 (/ (+ x y) x)) 1.0))))
double code(double x, double y, double z, double t, double a, double b, double c) {
double tmp;
if (y <= -7.2e+216) {
tmp = (x * -0.5) / (a * (b * y));
} else if (y <= -4.7e+114) {
tmp = 1.0;
} else if (y <= -6e+41) {
tmp = 1.0 / ((x + 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 (y <= (-7.2d+216)) then
tmp = (x * (-0.5d0)) / (a * (b * y))
else if (y <= (-4.7d+114)) then
tmp = 1.0d0
else if (y <= (-6d+41)) then
tmp = 1.0d0 / ((x + 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 (y <= -7.2e+216) {
tmp = (x * -0.5) / (a * (b * y));
} else if (y <= -4.7e+114) {
tmp = 1.0;
} else if (y <= -6e+41) {
tmp = 1.0 / ((x + y) / x);
} else {
tmp = 1.0;
}
return tmp;
}
def code(x, y, z, t, a, b, c): tmp = 0 if y <= -7.2e+216: tmp = (x * -0.5) / (a * (b * y)) elif y <= -4.7e+114: tmp = 1.0 elif y <= -6e+41: tmp = 1.0 / ((x + y) / x) else: tmp = 1.0 return tmp
function code(x, y, z, t, a, b, c) tmp = 0.0 if (y <= -7.2e+216) tmp = Float64(Float64(x * -0.5) / Float64(a * Float64(b * y))); elseif (y <= -4.7e+114) tmp = 1.0; elseif (y <= -6e+41) tmp = Float64(1.0 / Float64(Float64(x + 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 (y <= -7.2e+216) tmp = (x * -0.5) / (a * (b * y)); elseif (y <= -4.7e+114) tmp = 1.0; elseif (y <= -6e+41) tmp = 1.0 / ((x + y) / x); else tmp = 1.0; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_] := If[LessEqual[y, -7.2e+216], N[(N[(x * -0.5), $MachinePrecision] / N[(a * N[(b * y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, -4.7e+114], 1.0, If[LessEqual[y, -6e+41], N[(1.0 / N[(N[(x + y), $MachinePrecision] / x), $MachinePrecision]), $MachinePrecision], 1.0]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -7.2 \cdot 10^{+216}:\\
\;\;\;\;\frac{x \cdot -0.5}{a \cdot \left(b \cdot y\right)}\\
\mathbf{elif}\;y \leq -4.7 \cdot 10^{+114}:\\
\;\;\;\;1\\
\mathbf{elif}\;y \leq -6 \cdot 10^{+41}:\\
\;\;\;\;\frac{1}{\frac{x + y}{x}}\\
\mathbf{else}:\\
\;\;\;\;1\\
\end{array}
\end{array}
if y < -7.2000000000000004e216Initial program 100.0%
Taylor expanded in b around inf 75.8%
associate-*r/75.8%
metadata-eval75.8%
+-commutative75.8%
Simplified75.8%
Taylor expanded in b around 0 81.8%
Taylor expanded in a around inf 70.1%
associate-*r/70.1%
Simplified70.1%
if -7.2000000000000004e216 < y < -4.7000000000000001e114 or -5.9999999999999997e41 < y Initial program 94.6%
Taylor expanded in c around inf 67.5%
+-commutative67.5%
associate-*r/67.5%
metadata-eval67.5%
Simplified67.5%
Taylor expanded in t around inf 55.7%
Taylor expanded in x around inf 57.4%
if -4.7000000000000001e114 < y < -5.9999999999999997e41Initial program 88.1%
Taylor expanded in t around 0 47.9%
Taylor expanded in t around inf 49.4%
clear-num49.4%
inv-pow49.4%
Applied egg-rr49.4%
unpow-149.4%
Simplified49.4%
Final simplification57.6%
(FPCore (x y z t a b c) :precision binary64 (if (<= y -5.6e+41) (/ x (- x (* y (- -1.0 (* -2.0 (* b (+ a 0.8333333333333334))))))) 1.0))
double code(double x, double y, double z, double t, double a, double b, double c) {
double tmp;
if (y <= -5.6e+41) {
tmp = x / (x - (y * (-1.0 - (-2.0 * (b * (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) :: tmp
if (y <= (-5.6d+41)) then
tmp = x / (x - (y * ((-1.0d0) - ((-2.0d0) * (b * (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 tmp;
if (y <= -5.6e+41) {
tmp = x / (x - (y * (-1.0 - (-2.0 * (b * (a + 0.8333333333333334))))));
} else {
tmp = 1.0;
}
return tmp;
}
def code(x, y, z, t, a, b, c): tmp = 0 if y <= -5.6e+41: tmp = x / (x - (y * (-1.0 - (-2.0 * (b * (a + 0.8333333333333334)))))) else: tmp = 1.0 return tmp
function code(x, y, z, t, a, b, c) tmp = 0.0 if (y <= -5.6e+41) tmp = Float64(x / Float64(x - Float64(y * Float64(-1.0 - Float64(-2.0 * Float64(b * Float64(a + 0.8333333333333334))))))); else tmp = 1.0; end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c) tmp = 0.0; if (y <= -5.6e+41) tmp = x / (x - (y * (-1.0 - (-2.0 * (b * (a + 0.8333333333333334)))))); else tmp = 1.0; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_] := If[LessEqual[y, -5.6e+41], N[(x / N[(x - N[(y * N[(-1.0 - N[(-2.0 * N[(b * N[(a + 0.8333333333333334), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 1.0]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -5.6 \cdot 10^{+41}:\\
\;\;\;\;\frac{x}{x - y \cdot \left(-1 - -2 \cdot \left(b \cdot \left(a + 0.8333333333333334\right)\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;1\\
\end{array}
\end{array}
if y < -5.5999999999999999e41Initial program 93.7%
Taylor expanded in b around inf 75.7%
associate-*r/75.7%
metadata-eval75.7%
+-commutative75.7%
Simplified75.7%
Taylor expanded in b around 0 67.8%
Taylor expanded in t around inf 65.9%
*-commutative65.9%
*-commutative65.9%
Simplified65.9%
if -5.5999999999999999e41 < y Initial program 94.7%
Taylor expanded in c around inf 67.9%
+-commutative67.9%
associate-*r/67.9%
metadata-eval67.9%
Simplified67.9%
Taylor expanded in t around inf 54.9%
Taylor expanded in x around inf 57.1%
Final simplification58.7%
(FPCore (x y z t a b c) :precision binary64 (if (<= y -5.9e+41) (/ x (+ x (* -2.0 (* a (* b y))))) 1.0))
double code(double x, double y, double z, double t, double a, double b, double c) {
double tmp;
if (y <= -5.9e+41) {
tmp = x / (x + (-2.0 * (a * (b * 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 (y <= (-5.9d+41)) then
tmp = x / (x + ((-2.0d0) * (a * (b * 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 (y <= -5.9e+41) {
tmp = x / (x + (-2.0 * (a * (b * y))));
} else {
tmp = 1.0;
}
return tmp;
}
def code(x, y, z, t, a, b, c): tmp = 0 if y <= -5.9e+41: tmp = x / (x + (-2.0 * (a * (b * y)))) else: tmp = 1.0 return tmp
function code(x, y, z, t, a, b, c) tmp = 0.0 if (y <= -5.9e+41) tmp = Float64(x / Float64(x + Float64(-2.0 * Float64(a * Float64(b * y))))); else tmp = 1.0; end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c) tmp = 0.0; if (y <= -5.9e+41) tmp = x / (x + (-2.0 * (a * (b * y)))); else tmp = 1.0; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_] := If[LessEqual[y, -5.9e+41], N[(x / N[(x + N[(-2.0 * N[(a * N[(b * y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 1.0]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -5.9 \cdot 10^{+41}:\\
\;\;\;\;\frac{x}{x + -2 \cdot \left(a \cdot \left(b \cdot y\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;1\\
\end{array}
\end{array}
if y < -5.9000000000000001e41Initial program 93.7%
Taylor expanded in b around inf 75.7%
associate-*r/75.7%
metadata-eval75.7%
+-commutative75.7%
Simplified75.7%
Taylor expanded in b around 0 67.8%
Taylor expanded in a around inf 55.7%
if -5.9000000000000001e41 < y Initial program 94.7%
Taylor expanded in c around inf 67.9%
+-commutative67.9%
associate-*r/67.9%
metadata-eval67.9%
Simplified67.9%
Taylor expanded in t around inf 54.9%
Taylor expanded in x around inf 57.1%
Final simplification56.8%
(FPCore (x y z t a b c) :precision binary64 (if (<= y -4e+41) (/ x (+ x (* b (* y (* a -2.0))))) 1.0))
double code(double x, double y, double z, double t, double a, double b, double c) {
double tmp;
if (y <= -4e+41) {
tmp = x / (x + (b * (y * (a * -2.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) :: tmp
if (y <= (-4d+41)) then
tmp = x / (x + (b * (y * (a * (-2.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 tmp;
if (y <= -4e+41) {
tmp = x / (x + (b * (y * (a * -2.0))));
} else {
tmp = 1.0;
}
return tmp;
}
def code(x, y, z, t, a, b, c): tmp = 0 if y <= -4e+41: tmp = x / (x + (b * (y * (a * -2.0)))) else: tmp = 1.0 return tmp
function code(x, y, z, t, a, b, c) tmp = 0.0 if (y <= -4e+41) tmp = Float64(x / Float64(x + Float64(b * Float64(y * Float64(a * -2.0))))); else tmp = 1.0; end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c) tmp = 0.0; if (y <= -4e+41) tmp = x / (x + (b * (y * (a * -2.0)))); else tmp = 1.0; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_] := If[LessEqual[y, -4e+41], N[(x / N[(x + N[(b * N[(y * N[(a * -2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 1.0]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -4 \cdot 10^{+41}:\\
\;\;\;\;\frac{x}{x + b \cdot \left(y \cdot \left(a \cdot -2\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;1\\
\end{array}
\end{array}
if y < -4.00000000000000002e41Initial program 93.7%
Taylor expanded in b around inf 75.7%
associate-*r/75.7%
metadata-eval75.7%
+-commutative75.7%
Simplified75.7%
Taylor expanded in b around 0 67.8%
Taylor expanded in b around inf 59.6%
Taylor expanded in a around inf 56.0%
associate-*r*56.0%
Simplified56.0%
if -4.00000000000000002e41 < y Initial program 94.7%
Taylor expanded in c around inf 67.9%
+-commutative67.9%
associate-*r/67.9%
metadata-eval67.9%
Simplified67.9%
Taylor expanded in t around inf 54.9%
Taylor expanded in x around inf 57.1%
Final simplification56.9%
(FPCore (x y z t a b c) :precision binary64 (if (<= y -6e+215) (* -0.75 (* t (/ x (* c y)))) 1.0))
double code(double x, double y, double z, double t, double a, double b, double c) {
double tmp;
if (y <= -6e+215) {
tmp = -0.75 * (t * (x / (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 (y <= (-6d+215)) then
tmp = (-0.75d0) * (t * (x / (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 (y <= -6e+215) {
tmp = -0.75 * (t * (x / (c * y)));
} else {
tmp = 1.0;
}
return tmp;
}
def code(x, y, z, t, a, b, c): tmp = 0 if y <= -6e+215: tmp = -0.75 * (t * (x / (c * y))) else: tmp = 1.0 return tmp
function code(x, y, z, t, a, b, c) tmp = 0.0 if (y <= -6e+215) tmp = Float64(-0.75 * Float64(t * Float64(x / 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 (y <= -6e+215) tmp = -0.75 * (t * (x / (c * y))); else tmp = 1.0; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_] := If[LessEqual[y, -6e+215], N[(-0.75 * N[(t * N[(x / N[(c * y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 1.0]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -6 \cdot 10^{+215}:\\
\;\;\;\;-0.75 \cdot \left(t \cdot \frac{x}{c \cdot y}\right)\\
\mathbf{else}:\\
\;\;\;\;1\\
\end{array}
\end{array}
if y < -5.9999999999999998e215Initial program 100.0%
Taylor expanded in t around 0 56.8%
Taylor expanded in t around inf 56.8%
Taylor expanded in c around inf 51.7%
associate-/l*57.6%
*-commutative57.6%
Simplified57.6%
if -5.9999999999999998e215 < y Initial program 94.2%
Taylor expanded in c around inf 67.8%
+-commutative67.8%
associate-*r/67.8%
metadata-eval67.8%
Simplified67.8%
Taylor expanded in t around inf 57.3%
Taylor expanded in x around inf 54.8%
Final simplification54.9%
(FPCore (x y z t a b c) :precision binary64 (if (<= y -1.65e+217) (* -0.5 (/ (/ x a) (* b y))) 1.0))
double code(double x, double y, double z, double t, double a, double b, double c) {
double tmp;
if (y <= -1.65e+217) {
tmp = -0.5 * ((x / a) / (b * 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 (y <= (-1.65d+217)) then
tmp = (-0.5d0) * ((x / a) / (b * 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 (y <= -1.65e+217) {
tmp = -0.5 * ((x / a) / (b * y));
} else {
tmp = 1.0;
}
return tmp;
}
def code(x, y, z, t, a, b, c): tmp = 0 if y <= -1.65e+217: tmp = -0.5 * ((x / a) / (b * y)) else: tmp = 1.0 return tmp
function code(x, y, z, t, a, b, c) tmp = 0.0 if (y <= -1.65e+217) tmp = Float64(-0.5 * Float64(Float64(x / a) / Float64(b * y))); else tmp = 1.0; end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c) tmp = 0.0; if (y <= -1.65e+217) tmp = -0.5 * ((x / a) / (b * y)); else tmp = 1.0; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_] := If[LessEqual[y, -1.65e+217], N[(-0.5 * N[(N[(x / a), $MachinePrecision] / N[(b * y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 1.0]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -1.65 \cdot 10^{+217}:\\
\;\;\;\;-0.5 \cdot \frac{\frac{x}{a}}{b \cdot y}\\
\mathbf{else}:\\
\;\;\;\;1\\
\end{array}
\end{array}
if y < -1.65e217Initial program 100.0%
Taylor expanded in b around inf 75.8%
associate-*r/75.8%
metadata-eval75.8%
+-commutative75.8%
Simplified75.8%
Taylor expanded in b around 0 81.8%
Taylor expanded in a around inf 70.1%
associate-/r*58.4%
Simplified58.4%
if -1.65e217 < y Initial program 94.2%
Taylor expanded in c around inf 67.8%
+-commutative67.8%
associate-*r/67.8%
metadata-eval67.8%
Simplified67.8%
Taylor expanded in t around inf 57.3%
Taylor expanded in x around inf 54.8%
Final simplification55.0%
(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 94.5%
Taylor expanded in c around inf 67.9%
+-commutative67.9%
associate-*r/67.9%
metadata-eval67.9%
Simplified67.9%
Taylor expanded in t around inf 57.7%
Taylor expanded in x around inf 52.7%
Final simplification52.7%
(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 2024078
(FPCore (x y z t a b c)
:name "Numeric.SpecFunctions:invIncompleteBetaWorker from math-functions-0.1.5.2, I"
:precision binary64
:alt
(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 (/ (- (* (* z (sqrt (+ t a))) (* (* 3.0 t) (- a (/ 5.0 6.0)))) (* (- (* (+ (/ 5.0 6.0) a) (* 3.0 t)) 2.0) (* (- a (/ 5.0 6.0)) (* (- b c) t)))) (* (* (* t t) 3.0) (- a (/ 5.0 6.0))))))))) (/ x (+ x (* y (exp (* 2.0 (- (/ (* z (sqrt (+ t a))) t) (* (- b c) (- (+ a (/ 5.0 6.0)) (/ 2.0 (* t 3.0))))))))))))
(/ x (+ x (* y (exp (* 2.0 (- (/ (* z (sqrt (+ t a))) t) (* (- b c) (- (+ a (/ 5.0 6.0)) (/ 2.0 (* t 3.0)))))))))))