
(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 22 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))))) 1.0)))
double code(double x, double y, double z, double t, double a, double b, double c) {
double t_1 = ((z * sqrt((t + a))) / t) + ((b - c) * ((2.0 / (t * 3.0)) - (a + 0.8333333333333334)));
double tmp;
if (t_1 <= ((double) INFINITY)) {
tmp = x / (x + (y * exp((2.0 * t_1))));
} else {
tmp = 1.0;
}
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 = 1.0;
}
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 = 1.0 return tmp
function code(x, y, z, t, a, b, c) t_1 = Float64(Float64(Float64(z * sqrt(Float64(t + a))) / t) + Float64(Float64(b - c) * Float64(Float64(2.0 / Float64(t * 3.0)) - Float64(a + 0.8333333333333334)))) tmp = 0.0 if (t_1 <= Inf) tmp = Float64(x / Float64(x + Float64(y * exp(Float64(2.0 * t_1))))); else tmp = 1.0; 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 = 1.0; 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], 1.0]]
\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}:\\
\;\;\;\;1\\
\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 98.8%
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%
Simplified60.3%
Taylor expanded in x around inf 70.9%
Final simplification97.7%
(FPCore (x y z t a b c)
:precision binary64
(/
x
(fma
y
(pow
(exp 2.0)
(fma
z
(/ (sqrt (+ t a)) t)
(* (+ a (- 0.8333333333333334 (/ 0.6666666666666666 t))) (- c b))))
x)))
double code(double x, double y, double z, double t, double a, double b, double c) {
return x / fma(y, pow(exp(2.0), fma(z, (sqrt((t + a)) / t), ((a + (0.8333333333333334 - (0.6666666666666666 / t))) * (c - b)))), x);
}
function code(x, y, z, t, a, b, c) return Float64(x / fma(y, (exp(2.0) ^ fma(z, Float64(sqrt(Float64(t + a)) / t), Float64(Float64(a + Float64(0.8333333333333334 - Float64(0.6666666666666666 / t))) * Float64(c - b)))), x)) end
code[x_, y_, z_, t_, a_, b_, c_] := N[(x / N[(y * N[Power[N[Exp[2.0], $MachinePrecision], N[(z * N[(N[Sqrt[N[(t + a), $MachinePrecision]], $MachinePrecision] / t), $MachinePrecision] + N[(N[(a + N[(0.8333333333333334 - N[(0.6666666666666666 / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(c - b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] + x), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{x}{\mathsf{fma}\left(y, {\left(e^{2}\right)}^{\left(\mathsf{fma}\left(z, \frac{\sqrt{t + a}}{t}, \left(a + \left(0.8333333333333334 - \frac{0.6666666666666666}{t}\right)\right) \cdot \left(c - b\right)\right)\right)}, x\right)}
\end{array}
Initial program 95.0%
Simplified97.3%
(FPCore (x y z t a b c)
:precision binary64
(if (or (<= c -6e+70) (not (<= c 5e+85)))
(/
x
(+
x
(*
y
(exp
(* 2.0 (* c (+ 0.8333333333333334 (- a (/ 0.6666666666666666 t)))))))))
(/
x
(+ x (* y (exp (* 2.0 (+ (/ (* z (sqrt (+ t a))) t) (* a (- c b))))))))))
double code(double x, double y, double z, double t, double a, double b, double c) {
double tmp;
if ((c <= -6e+70) || !(c <= 5e+85)) {
tmp = x / (x + (y * exp((2.0 * (c * (0.8333333333333334 + (a - (0.6666666666666666 / t))))))));
} else {
tmp = x / (x + (y * exp((2.0 * (((z * sqrt((t + a))) / t) + (a * (c - 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 ((c <= (-6d+70)) .or. (.not. (c <= 5d+85))) then
tmp = x / (x + (y * exp((2.0d0 * (c * (0.8333333333333334d0 + (a - (0.6666666666666666d0 / t))))))))
else
tmp = x / (x + (y * exp((2.0d0 * (((z * sqrt((t + a))) / t) + (a * (c - 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 ((c <= -6e+70) || !(c <= 5e+85)) {
tmp = x / (x + (y * Math.exp((2.0 * (c * (0.8333333333333334 + (a - (0.6666666666666666 / t))))))));
} else {
tmp = x / (x + (y * Math.exp((2.0 * (((z * Math.sqrt((t + a))) / t) + (a * (c - b)))))));
}
return tmp;
}
def code(x, y, z, t, a, b, c): tmp = 0 if (c <= -6e+70) or not (c <= 5e+85): tmp = x / (x + (y * math.exp((2.0 * (c * (0.8333333333333334 + (a - (0.6666666666666666 / t)))))))) else: tmp = x / (x + (y * math.exp((2.0 * (((z * math.sqrt((t + a))) / t) + (a * (c - b))))))) return tmp
function code(x, y, z, t, a, b, c) tmp = 0.0 if ((c <= -6e+70) || !(c <= 5e+85)) tmp = Float64(x / Float64(x + Float64(y * exp(Float64(2.0 * Float64(c * Float64(0.8333333333333334 + Float64(a - Float64(0.6666666666666666 / t))))))))); else tmp = Float64(x / Float64(x + Float64(y * exp(Float64(2.0 * Float64(Float64(Float64(z * sqrt(Float64(t + a))) / t) + Float64(a * Float64(c - b)))))))); end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c) tmp = 0.0; if ((c <= -6e+70) || ~((c <= 5e+85))) tmp = x / (x + (y * exp((2.0 * (c * (0.8333333333333334 + (a - (0.6666666666666666 / t)))))))); else tmp = x / (x + (y * exp((2.0 * (((z * sqrt((t + a))) / t) + (a * (c - b))))))); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_] := If[Or[LessEqual[c, -6e+70], N[Not[LessEqual[c, 5e+85]], $MachinePrecision]], N[(x / N[(x + N[(y * N[Exp[N[(2.0 * N[(c * N[(0.8333333333333334 + N[(a - N[(0.6666666666666666 / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 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[(a * N[(c - b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;c \leq -6 \cdot 10^{+70} \lor \neg \left(c \leq 5 \cdot 10^{+85}\right):\\
\;\;\;\;\frac{x}{x + y \cdot e^{2 \cdot \left(c \cdot \left(0.8333333333333334 + \left(a - \frac{0.6666666666666666}{t}\right)\right)\right)}}\\
\mathbf{else}:\\
\;\;\;\;\frac{x}{x + y \cdot e^{2 \cdot \left(\frac{z \cdot \sqrt{t + a}}{t} + a \cdot \left(c - b\right)\right)}}\\
\end{array}
\end{array}
if c < -5.99999999999999952e70 or 5.0000000000000001e85 < c Initial program 95.6%
Taylor expanded in c around inf 91.4%
associate--l+91.4%
associate-*r/91.4%
metadata-eval91.4%
Simplified91.4%
if -5.99999999999999952e70 < c < 5.0000000000000001e85Initial program 94.4%
Taylor expanded in a around inf 83.7%
Final simplification87.1%
(FPCore (x y z t a b c)
:precision binary64
(if (or (<= c -2.2e-32) (not (<= c 2.5e+87)))
(/
x
(+
x
(*
y
(exp
(* 2.0 (* c (+ 0.8333333333333334 (- a (/ 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 <= -2.2e-32) || !(c <= 2.5e+87)) {
tmp = x / (x + (y * exp((2.0 * (c * (0.8333333333333334 + (a - (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 <= (-2.2d-32)) .or. (.not. (c <= 2.5d+87))) then
tmp = x / (x + (y * exp((2.0d0 * (c * (0.8333333333333334d0 + (a - (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 <= -2.2e-32) || !(c <= 2.5e+87)) {
tmp = x / (x + (y * Math.exp((2.0 * (c * (0.8333333333333334 + (a - (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 <= -2.2e-32) or not (c <= 2.5e+87): tmp = x / (x + (y * math.exp((2.0 * (c * (0.8333333333333334 + (a - (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 <= -2.2e-32) || !(c <= 2.5e+87)) tmp = Float64(x / Float64(x + Float64(y * exp(Float64(2.0 * Float64(c * Float64(0.8333333333333334 + Float64(a - 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 <= -2.2e-32) || ~((c <= 2.5e+87))) tmp = x / (x + (y * exp((2.0 * (c * (0.8333333333333334 + (a - (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, -2.2e-32], N[Not[LessEqual[c, 2.5e+87]], $MachinePrecision]], N[(x / N[(x + N[(y * N[Exp[N[(2.0 * N[(c * N[(0.8333333333333334 + N[(a - N[(0.6666666666666666 / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x / N[(x + N[(y * N[Exp[N[(2.0 * N[(b * N[(N[(0.6666666666666666 / t), $MachinePrecision] - N[(a + 0.8333333333333334), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;c \leq -2.2 \cdot 10^{-32} \lor \neg \left(c \leq 2.5 \cdot 10^{+87}\right):\\
\;\;\;\;\frac{x}{x + y \cdot e^{2 \cdot \left(c \cdot \left(0.8333333333333334 + \left(a - \frac{0.6666666666666666}{t}\right)\right)\right)}}\\
\mathbf{else}:\\
\;\;\;\;\frac{x}{x + y \cdot e^{2 \cdot \left(b \cdot \left(\frac{0.6666666666666666}{t} - \left(a + 0.8333333333333334\right)\right)\right)}}\\
\end{array}
\end{array}
if c < -2.2e-32 or 2.4999999999999999e87 < c Initial program 96.2%
Taylor expanded in c around inf 89.6%
associate--l+89.6%
associate-*r/89.6%
metadata-eval89.6%
Simplified89.6%
if -2.2e-32 < c < 2.4999999999999999e87Initial program 93.7%
Taylor expanded in b around inf 77.6%
associate-*r/77.6%
metadata-eval77.6%
Simplified77.6%
Final simplification83.7%
(FPCore (x y z t a b c)
:precision binary64
(if (or (<= c -3.4e+118) (not (<= c 7.8e+95)))
(/ x (+ x (* y (exp (* 2.0 (* c (+ a 0.8333333333333334)))))))
(/
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 <= -3.4e+118) || !(c <= 7.8e+95)) {
tmp = x / (x + (y * exp((2.0 * (c * (a + 0.8333333333333334))))));
} 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 <= (-3.4d+118)) .or. (.not. (c <= 7.8d+95))) then
tmp = x / (x + (y * exp((2.0d0 * (c * (a + 0.8333333333333334d0))))))
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 <= -3.4e+118) || !(c <= 7.8e+95)) {
tmp = x / (x + (y * Math.exp((2.0 * (c * (a + 0.8333333333333334))))));
} 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 <= -3.4e+118) or not (c <= 7.8e+95): tmp = x / (x + (y * math.exp((2.0 * (c * (a + 0.8333333333333334)))))) 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 <= -3.4e+118) || !(c <= 7.8e+95)) 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(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 <= -3.4e+118) || ~((c <= 7.8e+95))) tmp = x / (x + (y * exp((2.0 * (c * (a + 0.8333333333333334)))))); 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, -3.4e+118], N[Not[LessEqual[c, 7.8e+95]], $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[(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 -3.4 \cdot 10^{+118} \lor \neg \left(c \leq 7.8 \cdot 10^{+95}\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(\frac{0.6666666666666666}{t} - \left(a + 0.8333333333333334\right)\right)\right)}}\\
\end{array}
\end{array}
if c < -3.39999999999999986e118 or 7.7999999999999994e95 < c Initial program 96.1%
Taylor expanded in c around inf 93.3%
associate--l+93.3%
associate-*r/93.3%
metadata-eval93.3%
Simplified93.3%
Taylor expanded in t around inf 77.0%
+-commutative77.0%
Simplified77.0%
if -3.39999999999999986e118 < c < 7.7999999999999994e95Initial program 94.2%
Taylor expanded in b around inf 75.6%
associate-*r/75.6%
metadata-eval75.6%
Simplified75.6%
Final simplification76.2%
(FPCore (x y z t a b c)
:precision binary64
(if (<= c -2.9e+119)
1.0
(if (<= c 7.5e-177)
(/ x (+ x (* y (exp (* (* a b) (- 2.0))))))
(if (<= c 1.8e+87)
(/ x (+ x (* y (exp (* b -1.6666666666666667)))))
(/ x (+ x (* y (exp (* 2.0 (* a c))))))))))
double code(double x, double y, double z, double t, double a, double b, double c) {
double tmp;
if (c <= -2.9e+119) {
tmp = 1.0;
} else if (c <= 7.5e-177) {
tmp = x / (x + (y * exp(((a * b) * -2.0))));
} else if (c <= 1.8e+87) {
tmp = x / (x + (y * exp((b * -1.6666666666666667))));
} else {
tmp = x / (x + (y * exp((2.0 * (a * c)))));
}
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 <= (-2.9d+119)) then
tmp = 1.0d0
else if (c <= 7.5d-177) then
tmp = x / (x + (y * exp(((a * b) * -2.0d0))))
else if (c <= 1.8d+87) then
tmp = x / (x + (y * exp((b * (-1.6666666666666667d0)))))
else
tmp = x / (x + (y * exp((2.0d0 * (a * c)))))
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 <= -2.9e+119) {
tmp = 1.0;
} else if (c <= 7.5e-177) {
tmp = x / (x + (y * Math.exp(((a * b) * -2.0))));
} else if (c <= 1.8e+87) {
tmp = x / (x + (y * Math.exp((b * -1.6666666666666667))));
} else {
tmp = x / (x + (y * Math.exp((2.0 * (a * c)))));
}
return tmp;
}
def code(x, y, z, t, a, b, c): tmp = 0 if c <= -2.9e+119: tmp = 1.0 elif c <= 7.5e-177: tmp = x / (x + (y * math.exp(((a * b) * -2.0)))) elif c <= 1.8e+87: tmp = x / (x + (y * math.exp((b * -1.6666666666666667)))) else: tmp = x / (x + (y * math.exp((2.0 * (a * c))))) return tmp
function code(x, y, z, t, a, b, c) tmp = 0.0 if (c <= -2.9e+119) tmp = 1.0; elseif (c <= 7.5e-177) tmp = Float64(x / Float64(x + Float64(y * exp(Float64(Float64(a * b) * Float64(-2.0)))))); elseif (c <= 1.8e+87) tmp = Float64(x / Float64(x + Float64(y * exp(Float64(b * -1.6666666666666667))))); else tmp = Float64(x / Float64(x + Float64(y * exp(Float64(2.0 * Float64(a * c)))))); end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c) tmp = 0.0; if (c <= -2.9e+119) tmp = 1.0; elseif (c <= 7.5e-177) tmp = x / (x + (y * exp(((a * b) * -2.0)))); elseif (c <= 1.8e+87) tmp = x / (x + (y * exp((b * -1.6666666666666667)))); else tmp = x / (x + (y * exp((2.0 * (a * c))))); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_] := If[LessEqual[c, -2.9e+119], 1.0, If[LessEqual[c, 7.5e-177], N[(x / N[(x + N[(y * N[Exp[N[(N[(a * b), $MachinePrecision] * (-2.0)), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[c, 1.8e+87], N[(x / N[(x + N[(y * N[Exp[N[(b * -1.6666666666666667), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x / N[(x + N[(y * N[Exp[N[(2.0 * N[(a * c), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;c \leq -2.9 \cdot 10^{+119}:\\
\;\;\;\;1\\
\mathbf{elif}\;c \leq 7.5 \cdot 10^{-177}:\\
\;\;\;\;\frac{x}{x + y \cdot e^{\left(a \cdot b\right) \cdot \left(-2\right)}}\\
\mathbf{elif}\;c \leq 1.8 \cdot 10^{+87}:\\
\;\;\;\;\frac{x}{x + y \cdot e^{b \cdot -1.6666666666666667}}\\
\mathbf{else}:\\
\;\;\;\;\frac{x}{x + y \cdot e^{2 \cdot \left(a \cdot c\right)}}\\
\end{array}
\end{array}
if c < -2.90000000000000007e119Initial program 100.0%
Simplified100.0%
Taylor expanded in x around inf 80.1%
if -2.90000000000000007e119 < c < 7.5e-177Initial program 96.2%
Taylor expanded in b around inf 76.2%
associate-*r/76.2%
metadata-eval76.2%
Simplified76.2%
Taylor expanded in a around inf 67.2%
associate-*r*67.2%
neg-mul-167.2%
Simplified67.2%
if 7.5e-177 < c < 1.79999999999999997e87Initial program 89.7%
Taylor expanded in b around inf 71.8%
associate-*r/71.8%
metadata-eval71.8%
Simplified71.8%
Taylor expanded in t around inf 63.9%
Taylor expanded in a around 0 63.7%
*-commutative63.7%
Simplified63.7%
if 1.79999999999999997e87 < c Initial program 93.7%
Taylor expanded in c around inf 92.3%
associate--l+92.3%
associate-*r/92.3%
metadata-eval92.3%
Simplified92.3%
Taylor expanded in a around inf 73.2%
Final simplification70.0%
(FPCore (x y z t a b c) :precision binary64 (if (or (<= c -1.1e+120) (not (<= c 6e+78))) (/ x (+ x (* y (exp (* 2.0 (* c (+ a 0.8333333333333334))))))) (/ x (+ x (* y (exp (* -2.0 (* b (+ a 0.8333333333333334)))))))))
double code(double x, double y, double z, double t, double a, double b, double c) {
double tmp;
if ((c <= -1.1e+120) || !(c <= 6e+78)) {
tmp = x / (x + (y * exp((2.0 * (c * (a + 0.8333333333333334))))));
} else {
tmp = x / (x + (y * exp((-2.0 * (b * (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 <= (-1.1d+120)) .or. (.not. (c <= 6d+78))) then
tmp = x / (x + (y * exp((2.0d0 * (c * (a + 0.8333333333333334d0))))))
else
tmp = x / (x + (y * exp(((-2.0d0) * (b * (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 <= -1.1e+120) || !(c <= 6e+78)) {
tmp = x / (x + (y * Math.exp((2.0 * (c * (a + 0.8333333333333334))))));
} else {
tmp = x / (x + (y * Math.exp((-2.0 * (b * (a + 0.8333333333333334))))));
}
return tmp;
}
def code(x, y, z, t, a, b, c): tmp = 0 if (c <= -1.1e+120) or not (c <= 6e+78): tmp = x / (x + (y * math.exp((2.0 * (c * (a + 0.8333333333333334)))))) else: tmp = x / (x + (y * math.exp((-2.0 * (b * (a + 0.8333333333333334)))))) return tmp
function code(x, y, z, t, a, b, c) tmp = 0.0 if ((c <= -1.1e+120) || !(c <= 6e+78)) 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(a + 0.8333333333333334))))))); end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c) tmp = 0.0; if ((c <= -1.1e+120) || ~((c <= 6e+78))) tmp = x / (x + (y * exp((2.0 * (c * (a + 0.8333333333333334)))))); else tmp = x / (x + (y * exp((-2.0 * (b * (a + 0.8333333333333334)))))); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_] := If[Or[LessEqual[c, -1.1e+120], N[Not[LessEqual[c, 6e+78]], $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[(a + 0.8333333333333334), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;c \leq -1.1 \cdot 10^{+120} \lor \neg \left(c \leq 6 \cdot 10^{+78}\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(a + 0.8333333333333334\right)\right)}}\\
\end{array}
\end{array}
if c < -1.1000000000000001e120 or 5.99999999999999964e78 < c Initial program 94.3%
Taylor expanded in c around inf 90.7%
associate--l+90.7%
associate-*r/90.7%
metadata-eval90.7%
Simplified90.7%
Taylor expanded in t around inf 75.8%
+-commutative75.8%
Simplified75.8%
if -1.1000000000000001e120 < c < 5.99999999999999964e78Initial program 95.4%
Taylor expanded in b around inf 75.8%
associate-*r/75.8%
metadata-eval75.8%
Simplified75.8%
Taylor expanded in t around inf 66.9%
Final simplification70.5%
(FPCore (x y z t a b c)
:precision binary64
(if (<= c -4.2e+119)
1.0
(if (<= c 2.1e+84)
(/ x (+ x (* y (exp (* -2.0 (* b (+ a 0.8333333333333334)))))))
(/ x (+ x (* y (exp (* 2.0 (* a c)))))))))
double code(double x, double y, double z, double t, double a, double b, double c) {
double tmp;
if (c <= -4.2e+119) {
tmp = 1.0;
} else if (c <= 2.1e+84) {
tmp = x / (x + (y * exp((-2.0 * (b * (a + 0.8333333333333334))))));
} else {
tmp = x / (x + (y * exp((2.0 * (a * c)))));
}
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 <= (-4.2d+119)) then
tmp = 1.0d0
else if (c <= 2.1d+84) then
tmp = x / (x + (y * exp(((-2.0d0) * (b * (a + 0.8333333333333334d0))))))
else
tmp = x / (x + (y * exp((2.0d0 * (a * c)))))
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 <= -4.2e+119) {
tmp = 1.0;
} else if (c <= 2.1e+84) {
tmp = x / (x + (y * Math.exp((-2.0 * (b * (a + 0.8333333333333334))))));
} else {
tmp = x / (x + (y * Math.exp((2.0 * (a * c)))));
}
return tmp;
}
def code(x, y, z, t, a, b, c): tmp = 0 if c <= -4.2e+119: tmp = 1.0 elif c <= 2.1e+84: tmp = x / (x + (y * math.exp((-2.0 * (b * (a + 0.8333333333333334)))))) else: tmp = x / (x + (y * math.exp((2.0 * (a * c))))) return tmp
function code(x, y, z, t, a, b, c) tmp = 0.0 if (c <= -4.2e+119) tmp = 1.0; elseif (c <= 2.1e+84) tmp = Float64(x / Float64(x + Float64(y * exp(Float64(-2.0 * Float64(b * Float64(a + 0.8333333333333334))))))); else tmp = Float64(x / Float64(x + Float64(y * exp(Float64(2.0 * Float64(a * c)))))); end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c) tmp = 0.0; if (c <= -4.2e+119) tmp = 1.0; elseif (c <= 2.1e+84) tmp = x / (x + (y * exp((-2.0 * (b * (a + 0.8333333333333334)))))); else tmp = x / (x + (y * exp((2.0 * (a * c))))); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_] := If[LessEqual[c, -4.2e+119], 1.0, If[LessEqual[c, 2.1e+84], N[(x / N[(x + N[(y * N[Exp[N[(-2.0 * N[(b * N[(a + 0.8333333333333334), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x / N[(x + N[(y * N[Exp[N[(2.0 * N[(a * c), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;c \leq -4.2 \cdot 10^{+119}:\\
\;\;\;\;1\\
\mathbf{elif}\;c \leq 2.1 \cdot 10^{+84}:\\
\;\;\;\;\frac{x}{x + y \cdot e^{-2 \cdot \left(b \cdot \left(a + 0.8333333333333334\right)\right)}}\\
\mathbf{else}:\\
\;\;\;\;\frac{x}{x + y \cdot e^{2 \cdot \left(a \cdot c\right)}}\\
\end{array}
\end{array}
if c < -4.19999999999999966e119Initial program 100.0%
Simplified100.0%
Taylor expanded in x around inf 80.1%
if -4.19999999999999966e119 < c < 2.10000000000000019e84Initial program 95.4%
Taylor expanded in b around inf 75.8%
associate-*r/75.8%
metadata-eval75.8%
Simplified75.8%
Taylor expanded in t around inf 66.9%
if 2.10000000000000019e84 < c Initial program 90.8%
Taylor expanded in c around inf 89.6%
associate--l+89.6%
associate-*r/89.6%
metadata-eval89.6%
Simplified89.6%
Taylor expanded in a around inf 71.0%
Final simplification70.0%
(FPCore (x y z t a b c)
:precision binary64
(if (<= c -6.2e+59)
1.0
(if (<= c 1.85e+87)
(/ x (+ x (* y (exp (* b -1.6666666666666667)))))
(/ x (+ x (* y (exp (* 2.0 (* a c)))))))))
double code(double x, double y, double z, double t, double a, double b, double c) {
double tmp;
if (c <= -6.2e+59) {
tmp = 1.0;
} else if (c <= 1.85e+87) {
tmp = x / (x + (y * exp((b * -1.6666666666666667))));
} else {
tmp = x / (x + (y * exp((2.0 * (a * c)))));
}
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 <= (-6.2d+59)) then
tmp = 1.0d0
else if (c <= 1.85d+87) then
tmp = x / (x + (y * exp((b * (-1.6666666666666667d0)))))
else
tmp = x / (x + (y * exp((2.0d0 * (a * c)))))
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 <= -6.2e+59) {
tmp = 1.0;
} else if (c <= 1.85e+87) {
tmp = x / (x + (y * Math.exp((b * -1.6666666666666667))));
} else {
tmp = x / (x + (y * Math.exp((2.0 * (a * c)))));
}
return tmp;
}
def code(x, y, z, t, a, b, c): tmp = 0 if c <= -6.2e+59: tmp = 1.0 elif c <= 1.85e+87: tmp = x / (x + (y * math.exp((b * -1.6666666666666667)))) else: tmp = x / (x + (y * math.exp((2.0 * (a * c))))) return tmp
function code(x, y, z, t, a, b, c) tmp = 0.0 if (c <= -6.2e+59) tmp = 1.0; elseif (c <= 1.85e+87) tmp = Float64(x / Float64(x + Float64(y * exp(Float64(b * -1.6666666666666667))))); else tmp = Float64(x / Float64(x + Float64(y * exp(Float64(2.0 * Float64(a * c)))))); end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c) tmp = 0.0; if (c <= -6.2e+59) tmp = 1.0; elseif (c <= 1.85e+87) tmp = x / (x + (y * exp((b * -1.6666666666666667)))); else tmp = x / (x + (y * exp((2.0 * (a * c))))); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_] := If[LessEqual[c, -6.2e+59], 1.0, If[LessEqual[c, 1.85e+87], N[(x / N[(x + N[(y * N[Exp[N[(b * -1.6666666666666667), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x / N[(x + N[(y * N[Exp[N[(2.0 * N[(a * c), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;c \leq -6.2 \cdot 10^{+59}:\\
\;\;\;\;1\\
\mathbf{elif}\;c \leq 1.85 \cdot 10^{+87}:\\
\;\;\;\;\frac{x}{x + y \cdot e^{b \cdot -1.6666666666666667}}\\
\mathbf{else}:\\
\;\;\;\;\frac{x}{x + y \cdot e^{2 \cdot \left(a \cdot c\right)}}\\
\end{array}
\end{array}
if c < -6.20000000000000029e59Initial program 100.0%
Simplified100.0%
Taylor expanded in x around inf 79.1%
if -6.20000000000000029e59 < c < 1.85000000000000001e87Initial program 93.7%
Taylor expanded in b around inf 74.7%
associate-*r/74.7%
metadata-eval74.7%
Simplified74.7%
Taylor expanded in t around inf 65.3%
Taylor expanded in a around 0 58.5%
*-commutative58.5%
Simplified58.5%
if 1.85000000000000001e87 < c Initial program 93.7%
Taylor expanded in c around inf 92.3%
associate--l+92.3%
associate-*r/92.3%
metadata-eval92.3%
Simplified92.3%
Taylor expanded in a around inf 73.2%
(FPCore (x y z t a b c)
:precision binary64
(if (<= c -4e+57)
1.0
(if (<= c 8.5e+72)
(/ x (+ x (* y (exp (* b -1.6666666666666667)))))
(/ 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 (c <= -4e+57) {
tmp = 1.0;
} else if (c <= 8.5e+72) {
tmp = x / (x + (y * exp((b * -1.6666666666666667))));
} 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 (c <= (-4d+57)) then
tmp = 1.0d0
else if (c <= 8.5d+72) then
tmp = x / (x + (y * exp((b * (-1.6666666666666667d0)))))
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 (c <= -4e+57) {
tmp = 1.0;
} else if (c <= 8.5e+72) {
tmp = x / (x + (y * Math.exp((b * -1.6666666666666667))));
} else {
tmp = x / (x + (y * Math.exp((c * 1.6666666666666667))));
}
return tmp;
}
def code(x, y, z, t, a, b, c): tmp = 0 if c <= -4e+57: tmp = 1.0 elif c <= 8.5e+72: tmp = x / (x + (y * math.exp((b * -1.6666666666666667)))) 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 (c <= -4e+57) tmp = 1.0; elseif (c <= 8.5e+72) tmp = Float64(x / Float64(x + Float64(y * exp(Float64(b * -1.6666666666666667))))); 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 (c <= -4e+57) tmp = 1.0; elseif (c <= 8.5e+72) tmp = x / (x + (y * exp((b * -1.6666666666666667)))); else tmp = x / (x + (y * exp((c * 1.6666666666666667)))); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_] := If[LessEqual[c, -4e+57], 1.0, If[LessEqual[c, 8.5e+72], N[(x / N[(x + N[(y * N[Exp[N[(b * -1.6666666666666667), $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}\;c \leq -4 \cdot 10^{+57}:\\
\;\;\;\;1\\
\mathbf{elif}\;c \leq 8.5 \cdot 10^{+72}:\\
\;\;\;\;\frac{x}{x + y \cdot e^{b \cdot -1.6666666666666667}}\\
\mathbf{else}:\\
\;\;\;\;\frac{x}{x + y \cdot e^{c \cdot 1.6666666666666667}}\\
\end{array}
\end{array}
if c < -4.00000000000000019e57Initial program 100.0%
Simplified100.0%
Taylor expanded in x around inf 79.1%
if -4.00000000000000019e57 < c < 8.5000000000000004e72Initial program 95.0%
Taylor expanded in b around inf 75.6%
associate-*r/75.6%
metadata-eval75.6%
Simplified75.6%
Taylor expanded in t around inf 65.9%
Taylor expanded in a around 0 59.7%
*-commutative59.7%
Simplified59.7%
if 8.5000000000000004e72 < c Initial program 91.0%
Taylor expanded in c around inf 88.3%
associate--l+88.3%
associate-*r/88.3%
metadata-eval88.3%
Simplified88.3%
Taylor expanded in t around inf 70.6%
+-commutative70.6%
Simplified70.6%
Taylor expanded in a around 0 67.7%
*-commutative67.7%
Simplified67.7%
(FPCore (x y z t a b c)
:precision binary64
(if (<= c -1e+59)
1.0
(if (<= c 11600000.0)
(/ x (+ x (* y (exp (* b -1.6666666666666667)))))
(if (<= c 6e+124)
1.0
(/
x
(*
y
(-
1.0
(*
2.0
(/
(* c (- 0.6666666666666666 (* t (+ a 0.8333333333333334))))
t)))))))))
double code(double x, double y, double z, double t, double a, double b, double c) {
double tmp;
if (c <= -1e+59) {
tmp = 1.0;
} else if (c <= 11600000.0) {
tmp = x / (x + (y * exp((b * -1.6666666666666667))));
} else if (c <= 6e+124) {
tmp = 1.0;
} else {
tmp = x / (y * (1.0 - (2.0 * ((c * (0.6666666666666666 - (t * (a + 0.8333333333333334)))) / t))));
}
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 <= (-1d+59)) then
tmp = 1.0d0
else if (c <= 11600000.0d0) then
tmp = x / (x + (y * exp((b * (-1.6666666666666667d0)))))
else if (c <= 6d+124) then
tmp = 1.0d0
else
tmp = x / (y * (1.0d0 - (2.0d0 * ((c * (0.6666666666666666d0 - (t * (a + 0.8333333333333334d0)))) / t))))
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 <= -1e+59) {
tmp = 1.0;
} else if (c <= 11600000.0) {
tmp = x / (x + (y * Math.exp((b * -1.6666666666666667))));
} else if (c <= 6e+124) {
tmp = 1.0;
} else {
tmp = x / (y * (1.0 - (2.0 * ((c * (0.6666666666666666 - (t * (a + 0.8333333333333334)))) / t))));
}
return tmp;
}
def code(x, y, z, t, a, b, c): tmp = 0 if c <= -1e+59: tmp = 1.0 elif c <= 11600000.0: tmp = x / (x + (y * math.exp((b * -1.6666666666666667)))) elif c <= 6e+124: tmp = 1.0 else: tmp = x / (y * (1.0 - (2.0 * ((c * (0.6666666666666666 - (t * (a + 0.8333333333333334)))) / t)))) return tmp
function code(x, y, z, t, a, b, c) tmp = 0.0 if (c <= -1e+59) tmp = 1.0; elseif (c <= 11600000.0) tmp = Float64(x / Float64(x + Float64(y * exp(Float64(b * -1.6666666666666667))))); elseif (c <= 6e+124) tmp = 1.0; else tmp = Float64(x / Float64(y * Float64(1.0 - Float64(2.0 * Float64(Float64(c * Float64(0.6666666666666666 - Float64(t * Float64(a + 0.8333333333333334)))) / t))))); end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c) tmp = 0.0; if (c <= -1e+59) tmp = 1.0; elseif (c <= 11600000.0) tmp = x / (x + (y * exp((b * -1.6666666666666667)))); elseif (c <= 6e+124) tmp = 1.0; else tmp = x / (y * (1.0 - (2.0 * ((c * (0.6666666666666666 - (t * (a + 0.8333333333333334)))) / t)))); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_] := If[LessEqual[c, -1e+59], 1.0, If[LessEqual[c, 11600000.0], N[(x / N[(x + N[(y * N[Exp[N[(b * -1.6666666666666667), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[c, 6e+124], 1.0, N[(x / N[(y * N[(1.0 - N[(2.0 * N[(N[(c * N[(0.6666666666666666 - N[(t * N[(a + 0.8333333333333334), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;c \leq -1 \cdot 10^{+59}:\\
\;\;\;\;1\\
\mathbf{elif}\;c \leq 11600000:\\
\;\;\;\;\frac{x}{x + y \cdot e^{b \cdot -1.6666666666666667}}\\
\mathbf{elif}\;c \leq 6 \cdot 10^{+124}:\\
\;\;\;\;1\\
\mathbf{else}:\\
\;\;\;\;\frac{x}{y \cdot \left(1 - 2 \cdot \frac{c \cdot \left(0.6666666666666666 - t \cdot \left(a + 0.8333333333333334\right)\right)}{t}\right)}\\
\end{array}
\end{array}
if c < -9.99999999999999972e58 or 1.16e7 < c < 5.9999999999999999e124Initial program 95.2%
Simplified95.2%
Taylor expanded in x around inf 74.3%
if -9.99999999999999972e58 < c < 1.16e7Initial program 95.3%
Taylor expanded in b around inf 75.4%
associate-*r/75.4%
metadata-eval75.4%
Simplified75.4%
Taylor expanded in t around inf 65.5%
Taylor expanded in a around 0 58.6%
*-commutative58.6%
Simplified58.6%
if 5.9999999999999999e124 < c Initial program 93.6%
Taylor expanded in c around inf 93.8%
associate--l+93.8%
associate-*r/93.8%
metadata-eval93.8%
Simplified93.8%
Taylor expanded in c around 0 55.9%
associate-*r*55.9%
associate-*r/55.9%
metadata-eval55.9%
+-commutative55.9%
associate--l+55.9%
metadata-eval55.9%
associate-*r/55.9%
sub-neg55.9%
associate-*r/55.9%
metadata-eval55.9%
distribute-neg-frac55.9%
metadata-eval55.9%
Simplified55.9%
Taylor expanded in t around 0 53.8%
Taylor expanded in y around inf 71.3%
Final simplification66.0%
(FPCore (x y z t a b c)
:precision binary64
(let* ((t_1 (- 0.6666666666666666 (* t (+ a 0.8333333333333334)))))
(if (<= c -1.1e-103)
1.0
(if (<= c 520000000.0)
(/ x (- x (- (* (* 2.0 c) (* y (/ t_1 t))) y)))
(if (<= c 1e+125) 1.0 (/ x (* y (- 1.0 (* 2.0 (/ (* c t_1) t))))))))))
double code(double x, double y, double z, double t, double a, double b, double c) {
double t_1 = 0.6666666666666666 - (t * (a + 0.8333333333333334));
double tmp;
if (c <= -1.1e-103) {
tmp = 1.0;
} else if (c <= 520000000.0) {
tmp = x / (x - (((2.0 * c) * (y * (t_1 / t))) - y));
} else if (c <= 1e+125) {
tmp = 1.0;
} else {
tmp = x / (y * (1.0 - (2.0 * ((c * t_1) / t))));
}
return tmp;
}
real(8) function code(x, y, z, t, a, b, c)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8), intent (in) :: c
real(8) :: t_1
real(8) :: tmp
t_1 = 0.6666666666666666d0 - (t * (a + 0.8333333333333334d0))
if (c <= (-1.1d-103)) then
tmp = 1.0d0
else if (c <= 520000000.0d0) then
tmp = x / (x - (((2.0d0 * c) * (y * (t_1 / t))) - y))
else if (c <= 1d+125) then
tmp = 1.0d0
else
tmp = x / (y * (1.0d0 - (2.0d0 * ((c * t_1) / t))))
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b, double c) {
double t_1 = 0.6666666666666666 - (t * (a + 0.8333333333333334));
double tmp;
if (c <= -1.1e-103) {
tmp = 1.0;
} else if (c <= 520000000.0) {
tmp = x / (x - (((2.0 * c) * (y * (t_1 / t))) - y));
} else if (c <= 1e+125) {
tmp = 1.0;
} else {
tmp = x / (y * (1.0 - (2.0 * ((c * t_1) / t))));
}
return tmp;
}
def code(x, y, z, t, a, b, c): t_1 = 0.6666666666666666 - (t * (a + 0.8333333333333334)) tmp = 0 if c <= -1.1e-103: tmp = 1.0 elif c <= 520000000.0: tmp = x / (x - (((2.0 * c) * (y * (t_1 / t))) - y)) elif c <= 1e+125: tmp = 1.0 else: tmp = x / (y * (1.0 - (2.0 * ((c * t_1) / t)))) return tmp
function code(x, y, z, t, a, b, c) t_1 = Float64(0.6666666666666666 - Float64(t * Float64(a + 0.8333333333333334))) tmp = 0.0 if (c <= -1.1e-103) tmp = 1.0; elseif (c <= 520000000.0) tmp = Float64(x / Float64(x - Float64(Float64(Float64(2.0 * c) * Float64(y * Float64(t_1 / t))) - y))); elseif (c <= 1e+125) tmp = 1.0; else tmp = Float64(x / Float64(y * Float64(1.0 - Float64(2.0 * Float64(Float64(c * t_1) / t))))); end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c) t_1 = 0.6666666666666666 - (t * (a + 0.8333333333333334)); tmp = 0.0; if (c <= -1.1e-103) tmp = 1.0; elseif (c <= 520000000.0) tmp = x / (x - (((2.0 * c) * (y * (t_1 / t))) - y)); elseif (c <= 1e+125) tmp = 1.0; else tmp = x / (y * (1.0 - (2.0 * ((c * t_1) / t)))); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_] := Block[{t$95$1 = N[(0.6666666666666666 - N[(t * N[(a + 0.8333333333333334), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[c, -1.1e-103], 1.0, If[LessEqual[c, 520000000.0], N[(x / N[(x - N[(N[(N[(2.0 * c), $MachinePrecision] * N[(y * N[(t$95$1 / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[c, 1e+125], 1.0, N[(x / N[(y * N[(1.0 - N[(2.0 * N[(N[(c * t$95$1), $MachinePrecision] / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := 0.6666666666666666 - t \cdot \left(a + 0.8333333333333334\right)\\
\mathbf{if}\;c \leq -1.1 \cdot 10^{-103}:\\
\;\;\;\;1\\
\mathbf{elif}\;c \leq 520000000:\\
\;\;\;\;\frac{x}{x - \left(\left(2 \cdot c\right) \cdot \left(y \cdot \frac{t\_1}{t}\right) - y\right)}\\
\mathbf{elif}\;c \leq 10^{+125}:\\
\;\;\;\;1\\
\mathbf{else}:\\
\;\;\;\;\frac{x}{y \cdot \left(1 - 2 \cdot \frac{c \cdot t\_1}{t}\right)}\\
\end{array}
\end{array}
if c < -1.1e-103 or 5.2e8 < c < 9.9999999999999992e124Initial program 95.7%
Simplified94.8%
Taylor expanded in x around inf 68.0%
if -1.1e-103 < c < 5.2e8Initial program 94.7%
Taylor expanded in c around inf 54.7%
associate--l+54.7%
associate-*r/54.7%
metadata-eval54.7%
Simplified54.7%
Taylor expanded in c around 0 50.7%
associate-*r*50.7%
associate-*r/50.7%
metadata-eval50.7%
+-commutative50.7%
associate--l+50.7%
metadata-eval50.7%
associate-*r/50.7%
sub-neg50.7%
associate-*r/50.7%
metadata-eval50.7%
distribute-neg-frac50.7%
metadata-eval50.7%
Simplified50.7%
Taylor expanded in t around 0 54.8%
if 9.9999999999999992e124 < c Initial program 93.6%
Taylor expanded in c around inf 93.8%
associate--l+93.8%
associate-*r/93.8%
metadata-eval93.8%
Simplified93.8%
Taylor expanded in c around 0 55.9%
associate-*r*55.9%
associate-*r/55.9%
metadata-eval55.9%
+-commutative55.9%
associate--l+55.9%
metadata-eval55.9%
associate-*r/55.9%
sub-neg55.9%
associate-*r/55.9%
metadata-eval55.9%
distribute-neg-frac55.9%
metadata-eval55.9%
Simplified55.9%
Taylor expanded in t around 0 53.8%
Taylor expanded in y around inf 71.3%
Final simplification63.7%
(FPCore (x y z t a b c)
:precision binary64
(if (<= c -3.5e-125)
1.0
(if (<= c 160000000.0)
(/ x (+ x (* y (+ (* -2.0 (* a b)) 1.0))))
(if (<= c 1.85e+125)
1.0
(/
x
(*
y
(-
1.0
(*
2.0
(/
(* c (- 0.6666666666666666 (* t (+ a 0.8333333333333334))))
t)))))))))
double code(double x, double y, double z, double t, double a, double b, double c) {
double tmp;
if (c <= -3.5e-125) {
tmp = 1.0;
} else if (c <= 160000000.0) {
tmp = x / (x + (y * ((-2.0 * (a * b)) + 1.0)));
} else if (c <= 1.85e+125) {
tmp = 1.0;
} else {
tmp = x / (y * (1.0 - (2.0 * ((c * (0.6666666666666666 - (t * (a + 0.8333333333333334)))) / t))));
}
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 <= (-3.5d-125)) then
tmp = 1.0d0
else if (c <= 160000000.0d0) then
tmp = x / (x + (y * (((-2.0d0) * (a * b)) + 1.0d0)))
else if (c <= 1.85d+125) then
tmp = 1.0d0
else
tmp = x / (y * (1.0d0 - (2.0d0 * ((c * (0.6666666666666666d0 - (t * (a + 0.8333333333333334d0)))) / t))))
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 <= -3.5e-125) {
tmp = 1.0;
} else if (c <= 160000000.0) {
tmp = x / (x + (y * ((-2.0 * (a * b)) + 1.0)));
} else if (c <= 1.85e+125) {
tmp = 1.0;
} else {
tmp = x / (y * (1.0 - (2.0 * ((c * (0.6666666666666666 - (t * (a + 0.8333333333333334)))) / t))));
}
return tmp;
}
def code(x, y, z, t, a, b, c): tmp = 0 if c <= -3.5e-125: tmp = 1.0 elif c <= 160000000.0: tmp = x / (x + (y * ((-2.0 * (a * b)) + 1.0))) elif c <= 1.85e+125: tmp = 1.0 else: tmp = x / (y * (1.0 - (2.0 * ((c * (0.6666666666666666 - (t * (a + 0.8333333333333334)))) / t)))) return tmp
function code(x, y, z, t, a, b, c) tmp = 0.0 if (c <= -3.5e-125) tmp = 1.0; elseif (c <= 160000000.0) tmp = Float64(x / Float64(x + Float64(y * Float64(Float64(-2.0 * Float64(a * b)) + 1.0)))); elseif (c <= 1.85e+125) tmp = 1.0; else tmp = Float64(x / Float64(y * Float64(1.0 - Float64(2.0 * Float64(Float64(c * Float64(0.6666666666666666 - Float64(t * Float64(a + 0.8333333333333334)))) / t))))); end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c) tmp = 0.0; if (c <= -3.5e-125) tmp = 1.0; elseif (c <= 160000000.0) tmp = x / (x + (y * ((-2.0 * (a * b)) + 1.0))); elseif (c <= 1.85e+125) tmp = 1.0; else tmp = x / (y * (1.0 - (2.0 * ((c * (0.6666666666666666 - (t * (a + 0.8333333333333334)))) / t)))); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_] := If[LessEqual[c, -3.5e-125], 1.0, If[LessEqual[c, 160000000.0], N[(x / N[(x + N[(y * N[(N[(-2.0 * N[(a * b), $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[c, 1.85e+125], 1.0, N[(x / N[(y * N[(1.0 - N[(2.0 * N[(N[(c * N[(0.6666666666666666 - N[(t * N[(a + 0.8333333333333334), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;c \leq -3.5 \cdot 10^{-125}:\\
\;\;\;\;1\\
\mathbf{elif}\;c \leq 160000000:\\
\;\;\;\;\frac{x}{x + y \cdot \left(-2 \cdot \left(a \cdot b\right) + 1\right)}\\
\mathbf{elif}\;c \leq 1.85 \cdot 10^{+125}:\\
\;\;\;\;1\\
\mathbf{else}:\\
\;\;\;\;\frac{x}{y \cdot \left(1 - 2 \cdot \frac{c \cdot \left(0.6666666666666666 - t \cdot \left(a + 0.8333333333333334\right)\right)}{t}\right)}\\
\end{array}
\end{array}
if c < -3.49999999999999998e-125 or 1.6e8 < c < 1.8499999999999999e125Initial program 95.8%
Simplified95.0%
Taylor expanded in x around inf 68.3%
if -3.49999999999999998e-125 < c < 1.6e8Initial program 94.5%
Taylor expanded in b around inf 77.3%
associate-*r/77.3%
metadata-eval77.3%
Simplified77.3%
Taylor expanded in a around inf 59.2%
associate-*r*59.2%
neg-mul-159.2%
Simplified59.2%
Taylor expanded in a around 0 53.7%
if 1.8499999999999999e125 < c Initial program 93.6%
Taylor expanded in c around inf 93.8%
associate--l+93.8%
associate-*r/93.8%
metadata-eval93.8%
Simplified93.8%
Taylor expanded in c around 0 55.9%
associate-*r*55.9%
associate-*r/55.9%
metadata-eval55.9%
+-commutative55.9%
associate--l+55.9%
metadata-eval55.9%
associate-*r/55.9%
sub-neg55.9%
associate-*r/55.9%
metadata-eval55.9%
distribute-neg-frac55.9%
metadata-eval55.9%
Simplified55.9%
Taylor expanded in t around 0 53.8%
Taylor expanded in y around inf 71.3%
Final simplification63.7%
(FPCore (x y z t a b c)
:precision binary64
(if (<= c -7e-124)
1.0
(if (<= c 1.8e+18)
(/ x (+ x (* y (+ (* -2.0 (* a b)) 1.0))))
(if (<= c 1.26e+126)
1.0
(/
-1.0
(*
y
(/
(-
-1.0
(*
(* 2.0 c)
(+ 0.8333333333333334 (+ a (/ -0.6666666666666666 t)))))
x)))))))
double code(double x, double y, double z, double t, double a, double b, double c) {
double tmp;
if (c <= -7e-124) {
tmp = 1.0;
} else if (c <= 1.8e+18) {
tmp = x / (x + (y * ((-2.0 * (a * b)) + 1.0)));
} else if (c <= 1.26e+126) {
tmp = 1.0;
} else {
tmp = -1.0 / (y * ((-1.0 - ((2.0 * c) * (0.8333333333333334 + (a + (-0.6666666666666666 / t))))) / x));
}
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 <= (-7d-124)) then
tmp = 1.0d0
else if (c <= 1.8d+18) then
tmp = x / (x + (y * (((-2.0d0) * (a * b)) + 1.0d0)))
else if (c <= 1.26d+126) then
tmp = 1.0d0
else
tmp = (-1.0d0) / (y * (((-1.0d0) - ((2.0d0 * c) * (0.8333333333333334d0 + (a + ((-0.6666666666666666d0) / t))))) / x))
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 <= -7e-124) {
tmp = 1.0;
} else if (c <= 1.8e+18) {
tmp = x / (x + (y * ((-2.0 * (a * b)) + 1.0)));
} else if (c <= 1.26e+126) {
tmp = 1.0;
} else {
tmp = -1.0 / (y * ((-1.0 - ((2.0 * c) * (0.8333333333333334 + (a + (-0.6666666666666666 / t))))) / x));
}
return tmp;
}
def code(x, y, z, t, a, b, c): tmp = 0 if c <= -7e-124: tmp = 1.0 elif c <= 1.8e+18: tmp = x / (x + (y * ((-2.0 * (a * b)) + 1.0))) elif c <= 1.26e+126: tmp = 1.0 else: tmp = -1.0 / (y * ((-1.0 - ((2.0 * c) * (0.8333333333333334 + (a + (-0.6666666666666666 / t))))) / x)) return tmp
function code(x, y, z, t, a, b, c) tmp = 0.0 if (c <= -7e-124) tmp = 1.0; elseif (c <= 1.8e+18) tmp = Float64(x / Float64(x + Float64(y * Float64(Float64(-2.0 * Float64(a * b)) + 1.0)))); elseif (c <= 1.26e+126) tmp = 1.0; else tmp = Float64(-1.0 / Float64(y * Float64(Float64(-1.0 - Float64(Float64(2.0 * c) * Float64(0.8333333333333334 + Float64(a + Float64(-0.6666666666666666 / t))))) / x))); end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c) tmp = 0.0; if (c <= -7e-124) tmp = 1.0; elseif (c <= 1.8e+18) tmp = x / (x + (y * ((-2.0 * (a * b)) + 1.0))); elseif (c <= 1.26e+126) tmp = 1.0; else tmp = -1.0 / (y * ((-1.0 - ((2.0 * c) * (0.8333333333333334 + (a + (-0.6666666666666666 / t))))) / x)); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_] := If[LessEqual[c, -7e-124], 1.0, If[LessEqual[c, 1.8e+18], N[(x / N[(x + N[(y * N[(N[(-2.0 * N[(a * b), $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[c, 1.26e+126], 1.0, N[(-1.0 / N[(y * N[(N[(-1.0 - N[(N[(2.0 * c), $MachinePrecision] * N[(0.8333333333333334 + N[(a + N[(-0.6666666666666666 / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;c \leq -7 \cdot 10^{-124}:\\
\;\;\;\;1\\
\mathbf{elif}\;c \leq 1.8 \cdot 10^{+18}:\\
\;\;\;\;\frac{x}{x + y \cdot \left(-2 \cdot \left(a \cdot b\right) + 1\right)}\\
\mathbf{elif}\;c \leq 1.26 \cdot 10^{+126}:\\
\;\;\;\;1\\
\mathbf{else}:\\
\;\;\;\;\frac{-1}{y \cdot \frac{-1 - \left(2 \cdot c\right) \cdot \left(0.8333333333333334 + \left(a + \frac{-0.6666666666666666}{t}\right)\right)}{x}}\\
\end{array}
\end{array}
if c < -6.9999999999999997e-124 or 1.8e18 < c < 1.26000000000000004e126Initial program 95.9%
Simplified95.0%
Taylor expanded in x around inf 67.7%
if -6.9999999999999997e-124 < c < 1.8e18Initial program 94.5%
Taylor expanded in b around inf 77.3%
associate-*r/77.3%
metadata-eval77.3%
Simplified77.3%
Taylor expanded in a around inf 59.2%
associate-*r*59.2%
neg-mul-159.2%
Simplified59.2%
Taylor expanded in a around 0 53.7%
if 1.26000000000000004e126 < c Initial program 93.5%
Taylor expanded in c around inf 93.7%
associate--l+93.7%
associate-*r/93.7%
metadata-eval93.7%
Simplified93.7%
Taylor expanded in c around 0 57.0%
associate-*r*57.0%
associate-*r/57.0%
metadata-eval57.0%
+-commutative57.0%
associate--l+57.0%
metadata-eval57.0%
associate-*r/57.0%
sub-neg57.0%
associate-*r/57.0%
metadata-eval57.0%
distribute-neg-frac57.0%
metadata-eval57.0%
Simplified57.0%
clear-num58.2%
inv-pow58.2%
+-commutative58.2%
associate-*r*58.2%
fma-define58.2%
*-commutative58.2%
associate-+r+58.2%
Applied egg-rr58.2%
unpow-158.2%
associate-*l*58.2%
associate-+l+58.2%
metadata-eval58.2%
associate-*r/58.2%
metadata-eval58.2%
cancel-sign-sub-inv58.2%
associate-*r/58.2%
metadata-eval58.2%
Simplified58.2%
Taylor expanded in y around inf 62.4%
associate-/l*64.3%
associate-*r*64.3%
associate--l+64.3%
associate-*r/64.3%
metadata-eval64.3%
sub-neg64.3%
distribute-neg-frac64.3%
metadata-eval64.3%
Simplified64.3%
Final simplification62.2%
(FPCore (x y z t a b c)
:precision binary64
(if (<= c -2.5e-127)
1.0
(if (<= c 420000000.0)
(/ x (+ x (* y (+ (* -2.0 (* a b)) 1.0))))
(if (<= c 1.15e+126)
1.0
(/
(* x 0.5)
(* c (* y (+ a (- 0.8333333333333334 (/ 0.6666666666666666 t))))))))))
double code(double x, double y, double z, double t, double a, double b, double c) {
double tmp;
if (c <= -2.5e-127) {
tmp = 1.0;
} else if (c <= 420000000.0) {
tmp = x / (x + (y * ((-2.0 * (a * b)) + 1.0)));
} else if (c <= 1.15e+126) {
tmp = 1.0;
} else {
tmp = (x * 0.5) / (c * (y * (a + (0.8333333333333334 - (0.6666666666666666 / t)))));
}
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 <= (-2.5d-127)) then
tmp = 1.0d0
else if (c <= 420000000.0d0) then
tmp = x / (x + (y * (((-2.0d0) * (a * b)) + 1.0d0)))
else if (c <= 1.15d+126) then
tmp = 1.0d0
else
tmp = (x * 0.5d0) / (c * (y * (a + (0.8333333333333334d0 - (0.6666666666666666d0 / t)))))
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 <= -2.5e-127) {
tmp = 1.0;
} else if (c <= 420000000.0) {
tmp = x / (x + (y * ((-2.0 * (a * b)) + 1.0)));
} else if (c <= 1.15e+126) {
tmp = 1.0;
} else {
tmp = (x * 0.5) / (c * (y * (a + (0.8333333333333334 - (0.6666666666666666 / t)))));
}
return tmp;
}
def code(x, y, z, t, a, b, c): tmp = 0 if c <= -2.5e-127: tmp = 1.0 elif c <= 420000000.0: tmp = x / (x + (y * ((-2.0 * (a * b)) + 1.0))) elif c <= 1.15e+126: tmp = 1.0 else: tmp = (x * 0.5) / (c * (y * (a + (0.8333333333333334 - (0.6666666666666666 / t))))) return tmp
function code(x, y, z, t, a, b, c) tmp = 0.0 if (c <= -2.5e-127) tmp = 1.0; elseif (c <= 420000000.0) tmp = Float64(x / Float64(x + Float64(y * Float64(Float64(-2.0 * Float64(a * b)) + 1.0)))); elseif (c <= 1.15e+126) tmp = 1.0; else tmp = Float64(Float64(x * 0.5) / Float64(c * Float64(y * Float64(a + Float64(0.8333333333333334 - Float64(0.6666666666666666 / t)))))); end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c) tmp = 0.0; if (c <= -2.5e-127) tmp = 1.0; elseif (c <= 420000000.0) tmp = x / (x + (y * ((-2.0 * (a * b)) + 1.0))); elseif (c <= 1.15e+126) tmp = 1.0; else tmp = (x * 0.5) / (c * (y * (a + (0.8333333333333334 - (0.6666666666666666 / t))))); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_] := If[LessEqual[c, -2.5e-127], 1.0, If[LessEqual[c, 420000000.0], N[(x / N[(x + N[(y * N[(N[(-2.0 * N[(a * b), $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[c, 1.15e+126], 1.0, N[(N[(x * 0.5), $MachinePrecision] / N[(c * N[(y * N[(a + N[(0.8333333333333334 - N[(0.6666666666666666 / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;c \leq -2.5 \cdot 10^{-127}:\\
\;\;\;\;1\\
\mathbf{elif}\;c \leq 420000000:\\
\;\;\;\;\frac{x}{x + y \cdot \left(-2 \cdot \left(a \cdot b\right) + 1\right)}\\
\mathbf{elif}\;c \leq 1.15 \cdot 10^{+126}:\\
\;\;\;\;1\\
\mathbf{else}:\\
\;\;\;\;\frac{x \cdot 0.5}{c \cdot \left(y \cdot \left(a + \left(0.8333333333333334 - \frac{0.6666666666666666}{t}\right)\right)\right)}\\
\end{array}
\end{array}
if c < -2.4999999999999999e-127 or 4.2e8 < c < 1.15e126Initial program 95.9%
Simplified95.0%
Taylor expanded in x around inf 67.7%
if -2.4999999999999999e-127 < c < 4.2e8Initial program 94.5%
Taylor expanded in b around inf 77.3%
associate-*r/77.3%
metadata-eval77.3%
Simplified77.3%
Taylor expanded in a around inf 59.2%
associate-*r*59.2%
neg-mul-159.2%
Simplified59.2%
Taylor expanded in a around 0 53.7%
if 1.15e126 < c Initial program 93.5%
Taylor expanded in c around inf 93.7%
associate--l+93.7%
associate-*r/93.7%
metadata-eval93.7%
Simplified93.7%
Taylor expanded in c around 0 57.0%
associate-*r*57.0%
associate-*r/57.0%
metadata-eval57.0%
+-commutative57.0%
associate--l+57.0%
metadata-eval57.0%
associate-*r/57.0%
sub-neg57.0%
associate-*r/57.0%
metadata-eval57.0%
distribute-neg-frac57.0%
metadata-eval57.0%
Simplified57.0%
Taylor expanded in c around inf 54.9%
associate-*r/54.9%
cancel-sign-sub-inv54.9%
+-commutative54.9%
metadata-eval54.9%
associate-*r/54.9%
metadata-eval54.9%
associate-+l+54.9%
metadata-eval54.9%
associate-*r/54.9%
metadata-eval54.9%
cancel-sign-sub-inv54.9%
associate-*r/54.9%
metadata-eval54.9%
Simplified54.9%
Final simplification60.5%
(FPCore (x y z t a b c)
:precision binary64
(if (<= c -9.5e-124)
1.0
(if (<= c 1250000000000.0)
(/ x (+ x (* y (+ (* -2.0 (* a b)) 1.0))))
(if (<= c 1.95e+87) 1.0 (/ x (+ x (+ y (* (* 2.0 c) (* y a)))))))))
double code(double x, double y, double z, double t, double a, double b, double c) {
double tmp;
if (c <= -9.5e-124) {
tmp = 1.0;
} else if (c <= 1250000000000.0) {
tmp = x / (x + (y * ((-2.0 * (a * b)) + 1.0)));
} else if (c <= 1.95e+87) {
tmp = 1.0;
} else {
tmp = x / (x + (y + ((2.0 * c) * (y * 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 (c <= (-9.5d-124)) then
tmp = 1.0d0
else if (c <= 1250000000000.0d0) then
tmp = x / (x + (y * (((-2.0d0) * (a * b)) + 1.0d0)))
else if (c <= 1.95d+87) then
tmp = 1.0d0
else
tmp = x / (x + (y + ((2.0d0 * c) * (y * 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 (c <= -9.5e-124) {
tmp = 1.0;
} else if (c <= 1250000000000.0) {
tmp = x / (x + (y * ((-2.0 * (a * b)) + 1.0)));
} else if (c <= 1.95e+87) {
tmp = 1.0;
} else {
tmp = x / (x + (y + ((2.0 * c) * (y * a))));
}
return tmp;
}
def code(x, y, z, t, a, b, c): tmp = 0 if c <= -9.5e-124: tmp = 1.0 elif c <= 1250000000000.0: tmp = x / (x + (y * ((-2.0 * (a * b)) + 1.0))) elif c <= 1.95e+87: tmp = 1.0 else: tmp = x / (x + (y + ((2.0 * c) * (y * a)))) return tmp
function code(x, y, z, t, a, b, c) tmp = 0.0 if (c <= -9.5e-124) tmp = 1.0; elseif (c <= 1250000000000.0) tmp = Float64(x / Float64(x + Float64(y * Float64(Float64(-2.0 * Float64(a * b)) + 1.0)))); elseif (c <= 1.95e+87) tmp = 1.0; else tmp = Float64(x / Float64(x + Float64(y + Float64(Float64(2.0 * c) * Float64(y * a))))); end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c) tmp = 0.0; if (c <= -9.5e-124) tmp = 1.0; elseif (c <= 1250000000000.0) tmp = x / (x + (y * ((-2.0 * (a * b)) + 1.0))); elseif (c <= 1.95e+87) tmp = 1.0; else tmp = x / (x + (y + ((2.0 * c) * (y * a)))); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_] := If[LessEqual[c, -9.5e-124], 1.0, If[LessEqual[c, 1250000000000.0], N[(x / N[(x + N[(y * N[(N[(-2.0 * N[(a * b), $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[c, 1.95e+87], 1.0, N[(x / N[(x + N[(y + N[(N[(2.0 * c), $MachinePrecision] * N[(y * a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;c \leq -9.5 \cdot 10^{-124}:\\
\;\;\;\;1\\
\mathbf{elif}\;c \leq 1250000000000:\\
\;\;\;\;\frac{x}{x + y \cdot \left(-2 \cdot \left(a \cdot b\right) + 1\right)}\\
\mathbf{elif}\;c \leq 1.95 \cdot 10^{+87}:\\
\;\;\;\;1\\
\mathbf{else}:\\
\;\;\;\;\frac{x}{x + \left(y + \left(2 \cdot c\right) \cdot \left(y \cdot a\right)\right)}\\
\end{array}
\end{array}
if c < -9.49999999999999989e-124 or 1.25e12 < c < 1.9500000000000001e87Initial program 96.1%
Simplified96.1%
Taylor expanded in x around inf 69.9%
if -9.49999999999999989e-124 < c < 1.25e12Initial program 94.5%
Taylor expanded in b around inf 77.3%
associate-*r/77.3%
metadata-eval77.3%
Simplified77.3%
Taylor expanded in a around inf 59.2%
associate-*r*59.2%
neg-mul-159.2%
Simplified59.2%
Taylor expanded in a around 0 53.7%
if 1.9500000000000001e87 < c Initial program 93.7%
Taylor expanded in c around inf 92.3%
associate--l+92.3%
associate-*r/92.3%
metadata-eval92.3%
Simplified92.3%
Taylor expanded in c around 0 52.5%
associate-*r*52.5%
associate-*r/52.5%
metadata-eval52.5%
+-commutative52.5%
associate--l+52.5%
metadata-eval52.5%
associate-*r/52.5%
sub-neg52.5%
associate-*r/52.5%
metadata-eval52.5%
distribute-neg-frac52.5%
metadata-eval52.5%
Simplified52.5%
Taylor expanded in a around inf 54.7%
Final simplification60.5%
(FPCore (x y z t a b c)
:precision binary64
(if (<= c -5.4e-127)
1.0
(if (<= c 3.5e+14)
(/ x (+ x (* y (+ (* -2.0 (* a b)) 1.0))))
(if (<= c 1.75e+126) 1.0 (/ 1.0 (/ (* (* 2.0 a) (* y c)) x))))))
double code(double x, double y, double z, double t, double a, double b, double c) {
double tmp;
if (c <= -5.4e-127) {
tmp = 1.0;
} else if (c <= 3.5e+14) {
tmp = x / (x + (y * ((-2.0 * (a * b)) + 1.0)));
} else if (c <= 1.75e+126) {
tmp = 1.0;
} else {
tmp = 1.0 / (((2.0 * a) * (y * c)) / x);
}
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 <= (-5.4d-127)) then
tmp = 1.0d0
else if (c <= 3.5d+14) then
tmp = x / (x + (y * (((-2.0d0) * (a * b)) + 1.0d0)))
else if (c <= 1.75d+126) then
tmp = 1.0d0
else
tmp = 1.0d0 / (((2.0d0 * a) * (y * c)) / x)
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 <= -5.4e-127) {
tmp = 1.0;
} else if (c <= 3.5e+14) {
tmp = x / (x + (y * ((-2.0 * (a * b)) + 1.0)));
} else if (c <= 1.75e+126) {
tmp = 1.0;
} else {
tmp = 1.0 / (((2.0 * a) * (y * c)) / x);
}
return tmp;
}
def code(x, y, z, t, a, b, c): tmp = 0 if c <= -5.4e-127: tmp = 1.0 elif c <= 3.5e+14: tmp = x / (x + (y * ((-2.0 * (a * b)) + 1.0))) elif c <= 1.75e+126: tmp = 1.0 else: tmp = 1.0 / (((2.0 * a) * (y * c)) / x) return tmp
function code(x, y, z, t, a, b, c) tmp = 0.0 if (c <= -5.4e-127) tmp = 1.0; elseif (c <= 3.5e+14) tmp = Float64(x / Float64(x + Float64(y * Float64(Float64(-2.0 * Float64(a * b)) + 1.0)))); elseif (c <= 1.75e+126) tmp = 1.0; else tmp = Float64(1.0 / Float64(Float64(Float64(2.0 * a) * Float64(y * c)) / x)); end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c) tmp = 0.0; if (c <= -5.4e-127) tmp = 1.0; elseif (c <= 3.5e+14) tmp = x / (x + (y * ((-2.0 * (a * b)) + 1.0))); elseif (c <= 1.75e+126) tmp = 1.0; else tmp = 1.0 / (((2.0 * a) * (y * c)) / x); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_] := If[LessEqual[c, -5.4e-127], 1.0, If[LessEqual[c, 3.5e+14], N[(x / N[(x + N[(y * N[(N[(-2.0 * N[(a * b), $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[c, 1.75e+126], 1.0, N[(1.0 / N[(N[(N[(2.0 * a), $MachinePrecision] * N[(y * c), $MachinePrecision]), $MachinePrecision] / x), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;c \leq -5.4 \cdot 10^{-127}:\\
\;\;\;\;1\\
\mathbf{elif}\;c \leq 3.5 \cdot 10^{+14}:\\
\;\;\;\;\frac{x}{x + y \cdot \left(-2 \cdot \left(a \cdot b\right) + 1\right)}\\
\mathbf{elif}\;c \leq 1.75 \cdot 10^{+126}:\\
\;\;\;\;1\\
\mathbf{else}:\\
\;\;\;\;\frac{1}{\frac{\left(2 \cdot a\right) \cdot \left(y \cdot c\right)}{x}}\\
\end{array}
\end{array}
if c < -5.3999999999999999e-127 or 3.5e14 < c < 1.7500000000000001e126Initial program 95.9%
Simplified95.0%
Taylor expanded in x around inf 67.7%
if -5.3999999999999999e-127 < c < 3.5e14Initial program 94.5%
Taylor expanded in b around inf 77.3%
associate-*r/77.3%
metadata-eval77.3%
Simplified77.3%
Taylor expanded in a around inf 59.2%
associate-*r*59.2%
neg-mul-159.2%
Simplified59.2%
Taylor expanded in a around 0 53.7%
if 1.7500000000000001e126 < c Initial program 93.5%
Taylor expanded in c around inf 93.7%
associate--l+93.7%
associate-*r/93.7%
metadata-eval93.7%
Simplified93.7%
Taylor expanded in c around 0 57.0%
associate-*r*57.0%
associate-*r/57.0%
metadata-eval57.0%
+-commutative57.0%
associate--l+57.0%
metadata-eval57.0%
associate-*r/57.0%
sub-neg57.0%
associate-*r/57.0%
metadata-eval57.0%
distribute-neg-frac57.0%
metadata-eval57.0%
Simplified57.0%
clear-num58.2%
inv-pow58.2%
+-commutative58.2%
associate-*r*58.2%
fma-define58.2%
*-commutative58.2%
associate-+r+58.2%
Applied egg-rr58.2%
unpow-158.2%
associate-*l*58.2%
associate-+l+58.2%
metadata-eval58.2%
associate-*r/58.2%
metadata-eval58.2%
cancel-sign-sub-inv58.2%
associate-*r/58.2%
metadata-eval58.2%
Simplified58.2%
Taylor expanded in a around inf 49.8%
associate-*r/49.8%
associate-*r*49.8%
*-commutative49.8%
*-commutative49.8%
Simplified49.8%
Final simplification59.6%
(FPCore (x y z t a b c) :precision binary64 (if (<= c 2.2e+126) 1.0 (/ 1.0 (/ (* (* 2.0 a) (* y c)) x))))
double code(double x, double y, double z, double t, double a, double b, double c) {
double tmp;
if (c <= 2.2e+126) {
tmp = 1.0;
} else {
tmp = 1.0 / (((2.0 * a) * (y * c)) / x);
}
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 <= 2.2d+126) then
tmp = 1.0d0
else
tmp = 1.0d0 / (((2.0d0 * a) * (y * c)) / x)
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 <= 2.2e+126) {
tmp = 1.0;
} else {
tmp = 1.0 / (((2.0 * a) * (y * c)) / x);
}
return tmp;
}
def code(x, y, z, t, a, b, c): tmp = 0 if c <= 2.2e+126: tmp = 1.0 else: tmp = 1.0 / (((2.0 * a) * (y * c)) / x) return tmp
function code(x, y, z, t, a, b, c) tmp = 0.0 if (c <= 2.2e+126) tmp = 1.0; else tmp = Float64(1.0 / Float64(Float64(Float64(2.0 * a) * Float64(y * c)) / x)); end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c) tmp = 0.0; if (c <= 2.2e+126) tmp = 1.0; else tmp = 1.0 / (((2.0 * a) * (y * c)) / x); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_] := If[LessEqual[c, 2.2e+126], 1.0, N[(1.0 / N[(N[(N[(2.0 * a), $MachinePrecision] * N[(y * c), $MachinePrecision]), $MachinePrecision] / x), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;c \leq 2.2 \cdot 10^{+126}:\\
\;\;\;\;1\\
\mathbf{else}:\\
\;\;\;\;\frac{1}{\frac{\left(2 \cdot a\right) \cdot \left(y \cdot c\right)}{x}}\\
\end{array}
\end{array}
if c < 2.19999999999999999e126Initial program 95.3%
Simplified97.2%
Taylor expanded in x around inf 57.2%
if 2.19999999999999999e126 < c Initial program 93.5%
Taylor expanded in c around inf 93.7%
associate--l+93.7%
associate-*r/93.7%
metadata-eval93.7%
Simplified93.7%
Taylor expanded in c around 0 57.0%
associate-*r*57.0%
associate-*r/57.0%
metadata-eval57.0%
+-commutative57.0%
associate--l+57.0%
metadata-eval57.0%
associate-*r/57.0%
sub-neg57.0%
associate-*r/57.0%
metadata-eval57.0%
distribute-neg-frac57.0%
metadata-eval57.0%
Simplified57.0%
clear-num58.2%
inv-pow58.2%
+-commutative58.2%
associate-*r*58.2%
fma-define58.2%
*-commutative58.2%
associate-+r+58.2%
Applied egg-rr58.2%
unpow-158.2%
associate-*l*58.2%
associate-+l+58.2%
metadata-eval58.2%
associate-*r/58.2%
metadata-eval58.2%
cancel-sign-sub-inv58.2%
associate-*r/58.2%
metadata-eval58.2%
Simplified58.2%
Taylor expanded in a around inf 49.8%
associate-*r/49.8%
associate-*r*49.8%
*-commutative49.8%
*-commutative49.8%
Simplified49.8%
Final simplification55.9%
(FPCore (x y z t a b c) :precision binary64 (if (<= c 8.5e+125) 1.0 (* 0.5 (/ x (* a (* y c))))))
double code(double x, double y, double z, double t, double a, double b, double c) {
double tmp;
if (c <= 8.5e+125) {
tmp = 1.0;
} else {
tmp = 0.5 * (x / (a * (y * c)));
}
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.5d+125) then
tmp = 1.0d0
else
tmp = 0.5d0 * (x / (a * (y * c)))
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.5e+125) {
tmp = 1.0;
} else {
tmp = 0.5 * (x / (a * (y * c)));
}
return tmp;
}
def code(x, y, z, t, a, b, c): tmp = 0 if c <= 8.5e+125: tmp = 1.0 else: tmp = 0.5 * (x / (a * (y * c))) return tmp
function code(x, y, z, t, a, b, c) tmp = 0.0 if (c <= 8.5e+125) tmp = 1.0; else tmp = Float64(0.5 * Float64(x / Float64(a * Float64(y * c)))); end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c) tmp = 0.0; if (c <= 8.5e+125) tmp = 1.0; else tmp = 0.5 * (x / (a * (y * c))); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_] := If[LessEqual[c, 8.5e+125], 1.0, N[(0.5 * N[(x / N[(a * N[(y * c), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;c \leq 8.5 \cdot 10^{+125}:\\
\;\;\;\;1\\
\mathbf{else}:\\
\;\;\;\;0.5 \cdot \frac{x}{a \cdot \left(y \cdot c\right)}\\
\end{array}
\end{array}
if c < 8.49999999999999974e125Initial program 95.3%
Simplified97.2%
Taylor expanded in x around inf 57.2%
if 8.49999999999999974e125 < c Initial program 93.5%
Taylor expanded in c around inf 93.7%
associate--l+93.7%
associate-*r/93.7%
metadata-eval93.7%
Simplified93.7%
Taylor expanded in c around 0 57.0%
associate-*r*57.0%
associate-*r/57.0%
metadata-eval57.0%
+-commutative57.0%
associate--l+57.0%
metadata-eval57.0%
associate-*r/57.0%
sub-neg57.0%
associate-*r/57.0%
metadata-eval57.0%
distribute-neg-frac57.0%
metadata-eval57.0%
Simplified57.0%
Taylor expanded in a around inf 48.6%
Final simplification55.6%
(FPCore (x y z t a b c) :precision binary64 (if (<= b -3.8e+193) (* -0.5 (/ x (* a (* y b)))) 1.0))
double code(double x, double y, double z, double t, double a, double b, double c) {
double tmp;
if (b <= -3.8e+193) {
tmp = -0.5 * (x / (a * (y * 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 (b <= (-3.8d+193)) then
tmp = (-0.5d0) * (x / (a * (y * 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 (b <= -3.8e+193) {
tmp = -0.5 * (x / (a * (y * b)));
} else {
tmp = 1.0;
}
return tmp;
}
def code(x, y, z, t, a, b, c): tmp = 0 if b <= -3.8e+193: tmp = -0.5 * (x / (a * (y * b))) else: tmp = 1.0 return tmp
function code(x, y, z, t, a, b, c) tmp = 0.0 if (b <= -3.8e+193) tmp = Float64(-0.5 * Float64(x / Float64(a * Float64(y * b)))); else tmp = 1.0; end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c) tmp = 0.0; if (b <= -3.8e+193) tmp = -0.5 * (x / (a * (y * b))); else tmp = 1.0; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_] := If[LessEqual[b, -3.8e+193], N[(-0.5 * N[(x / N[(a * N[(y * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 1.0]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;b \leq -3.8 \cdot 10^{+193}:\\
\;\;\;\;-0.5 \cdot \frac{x}{a \cdot \left(y \cdot b\right)}\\
\mathbf{else}:\\
\;\;\;\;1\\
\end{array}
\end{array}
if b < -3.79999999999999972e193Initial program 93.1%
Taylor expanded in b around inf 100.0%
associate-*r/100.0%
metadata-eval100.0%
Simplified100.0%
Taylor expanded in a around inf 72.4%
associate-*r*72.4%
neg-mul-172.4%
Simplified72.4%
Taylor expanded in a around 0 65.4%
Taylor expanded in a around inf 51.7%
if -3.79999999999999972e193 < b Initial program 95.1%
Simplified97.1%
Taylor expanded in x around inf 53.6%
Final simplification53.5%
(FPCore (x y z t a b c) :precision binary64 (if (<= c 1e+126) 1.0 (* -0.75 (/ (* x t) (* y c)))))
double code(double x, double y, double z, double t, double a, double b, double c) {
double tmp;
if (c <= 1e+126) {
tmp = 1.0;
} else {
tmp = -0.75 * ((x * t) / (y * c));
}
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 <= 1d+126) then
tmp = 1.0d0
else
tmp = (-0.75d0) * ((x * t) / (y * c))
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 <= 1e+126) {
tmp = 1.0;
} else {
tmp = -0.75 * ((x * t) / (y * c));
}
return tmp;
}
def code(x, y, z, t, a, b, c): tmp = 0 if c <= 1e+126: tmp = 1.0 else: tmp = -0.75 * ((x * t) / (y * c)) return tmp
function code(x, y, z, t, a, b, c) tmp = 0.0 if (c <= 1e+126) tmp = 1.0; else tmp = Float64(-0.75 * Float64(Float64(x * t) / Float64(y * c))); end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c) tmp = 0.0; if (c <= 1e+126) tmp = 1.0; else tmp = -0.75 * ((x * t) / (y * c)); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_] := If[LessEqual[c, 1e+126], 1.0, N[(-0.75 * N[(N[(x * t), $MachinePrecision] / N[(y * c), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;c \leq 10^{+126}:\\
\;\;\;\;1\\
\mathbf{else}:\\
\;\;\;\;-0.75 \cdot \frac{x \cdot t}{y \cdot c}\\
\end{array}
\end{array}
if c < 9.99999999999999925e125Initial program 95.3%
Simplified97.2%
Taylor expanded in x around inf 57.2%
if 9.99999999999999925e125 < c Initial program 93.5%
Taylor expanded in c around inf 93.7%
associate--l+93.7%
associate-*r/93.7%
metadata-eval93.7%
Simplified93.7%
Taylor expanded in c around 0 57.0%
associate-*r*57.0%
associate-*r/57.0%
metadata-eval57.0%
+-commutative57.0%
associate--l+57.0%
metadata-eval57.0%
associate-*r/57.0%
sub-neg57.0%
associate-*r/57.0%
metadata-eval57.0%
distribute-neg-frac57.0%
metadata-eval57.0%
Simplified57.0%
Taylor expanded in t around 0 36.7%
Final simplification53.5%
(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 95.0%
Simplified97.3%
Taylor expanded in x around inf 51.6%
(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 2024180
(FPCore (x y z t a b c)
:name "Numeric.SpecFunctions:invIncompleteBetaWorker from math-functions-0.1.5.2, I"
:precision binary64
:alt
(! :herbie-platform default (if (< t -2118326644891581/100000000000000000000000000000000000000000000000000000000000000000) (/ x (+ x (* y (exp (* 2 (- (+ (* a c) (* 4166666666666667/5000000000000000 c)) (* a b))))))) (if (< t 5196588770651547/1000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000) (/ x (+ x (* y (exp (* 2 (/ (- (* (* z (sqrt (+ t a))) (* (* 3 t) (- a (/ 5 6)))) (* (- (* (+ (/ 5 6) a) (* 3 t)) 2) (* (- a (/ 5 6)) (* (- b c) t)))) (* (* (* t t) 3) (- a (/ 5 6))))))))) (/ x (+ x (* y (exp (* 2 (- (/ (* z (sqrt (+ t a))) t) (* (- b c) (- (+ a (/ 5 6)) (/ 2 (* t 3)))))))))))))
(/ x (+ x (* y (exp (* 2.0 (- (/ (* z (sqrt (+ t a))) t) (* (- b c) (- (+ a (/ 5.0 6.0)) (/ 2.0 (* t 3.0)))))))))))