
(FPCore (x y z t a b c i) :precision binary64 (/ (+ (* (+ (* (+ (* (+ (* x y) z) y) 27464.7644705) y) 230661.510616) y) t) (+ (* (+ (* (+ (* (+ y a) y) b) y) c) y) i)))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
return ((((((((x * y) + z) * y) + 27464.7644705) * y) + 230661.510616) * y) + t) / (((((((y + a) * y) + b) * y) + c) * y) + i);
}
real(8) function code(x, y, z, t, a, b, c, i)
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), intent (in) :: i
code = ((((((((x * y) + z) * y) + 27464.7644705d0) * y) + 230661.510616d0) * y) + t) / (((((((y + a) * y) + b) * y) + c) * y) + i)
end function
public static double code(double x, double y, double z, double t, double a, double b, double c, double i) {
return ((((((((x * y) + z) * y) + 27464.7644705) * y) + 230661.510616) * y) + t) / (((((((y + a) * y) + b) * y) + c) * y) + i);
}
def code(x, y, z, t, a, b, c, i): return ((((((((x * y) + z) * y) + 27464.7644705) * y) + 230661.510616) * y) + t) / (((((((y + a) * y) + b) * y) + c) * y) + i)
function code(x, y, z, t, a, b, c, i) return Float64(Float64(Float64(Float64(Float64(Float64(Float64(Float64(Float64(x * y) + z) * y) + 27464.7644705) * y) + 230661.510616) * y) + t) / Float64(Float64(Float64(Float64(Float64(Float64(Float64(y + a) * y) + b) * y) + c) * y) + i)) end
function tmp = code(x, y, z, t, a, b, c, i) tmp = ((((((((x * y) + z) * y) + 27464.7644705) * y) + 230661.510616) * y) + t) / (((((((y + a) * y) + b) * y) + c) * y) + i); end
code[x_, y_, z_, t_, a_, b_, c_, i_] := N[(N[(N[(N[(N[(N[(N[(N[(N[(x * y), $MachinePrecision] + z), $MachinePrecision] * y), $MachinePrecision] + 27464.7644705), $MachinePrecision] * y), $MachinePrecision] + 230661.510616), $MachinePrecision] * y), $MachinePrecision] + t), $MachinePrecision] / N[(N[(N[(N[(N[(N[(N[(y + a), $MachinePrecision] * y), $MachinePrecision] + b), $MachinePrecision] * y), $MachinePrecision] + c), $MachinePrecision] * y), $MachinePrecision] + i), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{\left(\left(\left(x \cdot y + z\right) \cdot y + 27464.7644705\right) \cdot y + 230661.510616\right) \cdot y + t}{\left(\left(\left(y + a\right) \cdot y + b\right) \cdot y + c\right) \cdot y + i}
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 27 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (x y z t a b c i) :precision binary64 (/ (+ (* (+ (* (+ (* (+ (* x y) z) y) 27464.7644705) y) 230661.510616) y) t) (+ (* (+ (* (+ (* (+ y a) y) b) y) c) y) i)))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
return ((((((((x * y) + z) * y) + 27464.7644705) * y) + 230661.510616) * y) + t) / (((((((y + a) * y) + b) * y) + c) * y) + i);
}
real(8) function code(x, y, z, t, a, b, c, i)
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), intent (in) :: i
code = ((((((((x * y) + z) * y) + 27464.7644705d0) * y) + 230661.510616d0) * y) + t) / (((((((y + a) * y) + b) * y) + c) * y) + i)
end function
public static double code(double x, double y, double z, double t, double a, double b, double c, double i) {
return ((((((((x * y) + z) * y) + 27464.7644705) * y) + 230661.510616) * y) + t) / (((((((y + a) * y) + b) * y) + c) * y) + i);
}
def code(x, y, z, t, a, b, c, i): return ((((((((x * y) + z) * y) + 27464.7644705) * y) + 230661.510616) * y) + t) / (((((((y + a) * y) + b) * y) + c) * y) + i)
function code(x, y, z, t, a, b, c, i) return Float64(Float64(Float64(Float64(Float64(Float64(Float64(Float64(Float64(x * y) + z) * y) + 27464.7644705) * y) + 230661.510616) * y) + t) / Float64(Float64(Float64(Float64(Float64(Float64(Float64(y + a) * y) + b) * y) + c) * y) + i)) end
function tmp = code(x, y, z, t, a, b, c, i) tmp = ((((((((x * y) + z) * y) + 27464.7644705) * y) + 230661.510616) * y) + t) / (((((((y + a) * y) + b) * y) + c) * y) + i); end
code[x_, y_, z_, t_, a_, b_, c_, i_] := N[(N[(N[(N[(N[(N[(N[(N[(N[(x * y), $MachinePrecision] + z), $MachinePrecision] * y), $MachinePrecision] + 27464.7644705), $MachinePrecision] * y), $MachinePrecision] + 230661.510616), $MachinePrecision] * y), $MachinePrecision] + t), $MachinePrecision] / N[(N[(N[(N[(N[(N[(N[(y + a), $MachinePrecision] * y), $MachinePrecision] + b), $MachinePrecision] * y), $MachinePrecision] + c), $MachinePrecision] * y), $MachinePrecision] + i), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{\left(\left(\left(x \cdot y + z\right) \cdot y + 27464.7644705\right) \cdot y + 230661.510616\right) \cdot y + t}{\left(\left(\left(y + a\right) \cdot y + b\right) \cdot y + c\right) \cdot y + i}
\end{array}
(FPCore (x y z t a b c i)
:precision binary64
(let* ((t_1 (* y (+ (* x y) z)))
(t_2 (+ x (- (/ z y) (* a (/ x y)))))
(t_3 (+ (* y (+ y a)) b))
(t_4 (* y t_3))
(t_5 (+ i (* y (+ c t_4)))))
(if (<= y -7.8e+135)
t_2
(if (<= y -2.2e+64)
(+
(* 27464.7644705 (/ 1.0 t_3))
(+
(* 230661.510616 (/ 1.0 t_4))
(+ (/ t (+ i (* t_3 (pow y 2.0)))) (/ t_1 t_3))))
(if (<= y 1.05e+44)
(+
(/ t t_5)
(/ (* y (+ 230661.510616 (* y (+ 27464.7644705 t_1)))) t_5))
t_2)))))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double t_1 = y * ((x * y) + z);
double t_2 = x + ((z / y) - (a * (x / y)));
double t_3 = (y * (y + a)) + b;
double t_4 = y * t_3;
double t_5 = i + (y * (c + t_4));
double tmp;
if (y <= -7.8e+135) {
tmp = t_2;
} else if (y <= -2.2e+64) {
tmp = (27464.7644705 * (1.0 / t_3)) + ((230661.510616 * (1.0 / t_4)) + ((t / (i + (t_3 * pow(y, 2.0)))) + (t_1 / t_3)));
} else if (y <= 1.05e+44) {
tmp = (t / t_5) + ((y * (230661.510616 + (y * (27464.7644705 + t_1)))) / t_5);
} else {
tmp = t_2;
}
return tmp;
}
real(8) function code(x, y, z, t, a, b, c, i)
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), intent (in) :: i
real(8) :: t_1
real(8) :: t_2
real(8) :: t_3
real(8) :: t_4
real(8) :: t_5
real(8) :: tmp
t_1 = y * ((x * y) + z)
t_2 = x + ((z / y) - (a * (x / y)))
t_3 = (y * (y + a)) + b
t_4 = y * t_3
t_5 = i + (y * (c + t_4))
if (y <= (-7.8d+135)) then
tmp = t_2
else if (y <= (-2.2d+64)) then
tmp = (27464.7644705d0 * (1.0d0 / t_3)) + ((230661.510616d0 * (1.0d0 / t_4)) + ((t / (i + (t_3 * (y ** 2.0d0)))) + (t_1 / t_3)))
else if (y <= 1.05d+44) then
tmp = (t / t_5) + ((y * (230661.510616d0 + (y * (27464.7644705d0 + t_1)))) / t_5)
else
tmp = t_2
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 i) {
double t_1 = y * ((x * y) + z);
double t_2 = x + ((z / y) - (a * (x / y)));
double t_3 = (y * (y + a)) + b;
double t_4 = y * t_3;
double t_5 = i + (y * (c + t_4));
double tmp;
if (y <= -7.8e+135) {
tmp = t_2;
} else if (y <= -2.2e+64) {
tmp = (27464.7644705 * (1.0 / t_3)) + ((230661.510616 * (1.0 / t_4)) + ((t / (i + (t_3 * Math.pow(y, 2.0)))) + (t_1 / t_3)));
} else if (y <= 1.05e+44) {
tmp = (t / t_5) + ((y * (230661.510616 + (y * (27464.7644705 + t_1)))) / t_5);
} else {
tmp = t_2;
}
return tmp;
}
def code(x, y, z, t, a, b, c, i): t_1 = y * ((x * y) + z) t_2 = x + ((z / y) - (a * (x / y))) t_3 = (y * (y + a)) + b t_4 = y * t_3 t_5 = i + (y * (c + t_4)) tmp = 0 if y <= -7.8e+135: tmp = t_2 elif y <= -2.2e+64: tmp = (27464.7644705 * (1.0 / t_3)) + ((230661.510616 * (1.0 / t_4)) + ((t / (i + (t_3 * math.pow(y, 2.0)))) + (t_1 / t_3))) elif y <= 1.05e+44: tmp = (t / t_5) + ((y * (230661.510616 + (y * (27464.7644705 + t_1)))) / t_5) else: tmp = t_2 return tmp
function code(x, y, z, t, a, b, c, i) t_1 = Float64(y * Float64(Float64(x * y) + z)) t_2 = Float64(x + Float64(Float64(z / y) - Float64(a * Float64(x / y)))) t_3 = Float64(Float64(y * Float64(y + a)) + b) t_4 = Float64(y * t_3) t_5 = Float64(i + Float64(y * Float64(c + t_4))) tmp = 0.0 if (y <= -7.8e+135) tmp = t_2; elseif (y <= -2.2e+64) tmp = Float64(Float64(27464.7644705 * Float64(1.0 / t_3)) + Float64(Float64(230661.510616 * Float64(1.0 / t_4)) + Float64(Float64(t / Float64(i + Float64(t_3 * (y ^ 2.0)))) + Float64(t_1 / t_3)))); elseif (y <= 1.05e+44) tmp = Float64(Float64(t / t_5) + Float64(Float64(y * Float64(230661.510616 + Float64(y * Float64(27464.7644705 + t_1)))) / t_5)); else tmp = t_2; end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c, i) t_1 = y * ((x * y) + z); t_2 = x + ((z / y) - (a * (x / y))); t_3 = (y * (y + a)) + b; t_4 = y * t_3; t_5 = i + (y * (c + t_4)); tmp = 0.0; if (y <= -7.8e+135) tmp = t_2; elseif (y <= -2.2e+64) tmp = (27464.7644705 * (1.0 / t_3)) + ((230661.510616 * (1.0 / t_4)) + ((t / (i + (t_3 * (y ^ 2.0)))) + (t_1 / t_3))); elseif (y <= 1.05e+44) tmp = (t / t_5) + ((y * (230661.510616 + (y * (27464.7644705 + t_1)))) / t_5); else tmp = t_2; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_, i_] := Block[{t$95$1 = N[(y * N[(N[(x * y), $MachinePrecision] + z), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(x + N[(N[(z / y), $MachinePrecision] - N[(a * N[(x / y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$3 = N[(N[(y * N[(y + a), $MachinePrecision]), $MachinePrecision] + b), $MachinePrecision]}, Block[{t$95$4 = N[(y * t$95$3), $MachinePrecision]}, Block[{t$95$5 = N[(i + N[(y * N[(c + t$95$4), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -7.8e+135], t$95$2, If[LessEqual[y, -2.2e+64], N[(N[(27464.7644705 * N[(1.0 / t$95$3), $MachinePrecision]), $MachinePrecision] + N[(N[(230661.510616 * N[(1.0 / t$95$4), $MachinePrecision]), $MachinePrecision] + N[(N[(t / N[(i + N[(t$95$3 * N[Power[y, 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(t$95$1 / t$95$3), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 1.05e+44], N[(N[(t / t$95$5), $MachinePrecision] + N[(N[(y * N[(230661.510616 + N[(y * N[(27464.7644705 + t$95$1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / t$95$5), $MachinePrecision]), $MachinePrecision], t$95$2]]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := y \cdot \left(x \cdot y + z\right)\\
t_2 := x + \left(\frac{z}{y} - a \cdot \frac{x}{y}\right)\\
t_3 := y \cdot \left(y + a\right) + b\\
t_4 := y \cdot t\_3\\
t_5 := i + y \cdot \left(c + t\_4\right)\\
\mathbf{if}\;y \leq -7.8 \cdot 10^{+135}:\\
\;\;\;\;t\_2\\
\mathbf{elif}\;y \leq -2.2 \cdot 10^{+64}:\\
\;\;\;\;27464.7644705 \cdot \frac{1}{t\_3} + \left(230661.510616 \cdot \frac{1}{t\_4} + \left(\frac{t}{i + t\_3 \cdot {y}^{2}} + \frac{t\_1}{t\_3}\right)\right)\\
\mathbf{elif}\;y \leq 1.05 \cdot 10^{+44}:\\
\;\;\;\;\frac{t}{t\_5} + \frac{y \cdot \left(230661.510616 + y \cdot \left(27464.7644705 + t\_1\right)\right)}{t\_5}\\
\mathbf{else}:\\
\;\;\;\;t\_2\\
\end{array}
\end{array}
if y < -7.80000000000000064e135 or 1.04999999999999993e44 < y Initial program 3.9%
Taylor expanded in y around inf 69.6%
associate--l+69.6%
associate-/l*79.3%
Simplified79.3%
if -7.80000000000000064e135 < y < -2.20000000000000002e64Initial program 1.8%
Taylor expanded in t around 0 1.8%
Taylor expanded in i around 0 35.3%
Taylor expanded in c around 0 67.6%
if -2.20000000000000002e64 < y < 1.04999999999999993e44Initial program 97.2%
Taylor expanded in t around 0 97.2%
Final simplification88.8%
(FPCore (x y z t a b c i)
:precision binary64
(let* ((t_1
(* y (+ 230661.510616 (* y (+ 27464.7644705 (* y (+ (* x y) z)))))))
(t_2 (+ i (* y (+ c (* y (+ (* y (+ y a)) b))))))
(t_3 (/ (+ t_1 t) t_2))
(t_4 (* x t_2)))
(if (<= t_3 (- INFINITY))
(*
x
(+
(/ t t_4)
(+
(/ (* y (+ 230661.510616 (* y 27464.7644705))) t_4)
(/ (pow y 4.0) t_2))))
(if (<= t_3 5e+243)
(+ (/ t t_2) (/ t_1 t_2))
(+ x (- (/ z y) (* a (/ x y))))))))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double t_1 = y * (230661.510616 + (y * (27464.7644705 + (y * ((x * y) + z)))));
double t_2 = i + (y * (c + (y * ((y * (y + a)) + b))));
double t_3 = (t_1 + t) / t_2;
double t_4 = x * t_2;
double tmp;
if (t_3 <= -((double) INFINITY)) {
tmp = x * ((t / t_4) + (((y * (230661.510616 + (y * 27464.7644705))) / t_4) + (pow(y, 4.0) / t_2)));
} else if (t_3 <= 5e+243) {
tmp = (t / t_2) + (t_1 / t_2);
} else {
tmp = x + ((z / y) - (a * (x / y)));
}
return tmp;
}
public static double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double t_1 = y * (230661.510616 + (y * (27464.7644705 + (y * ((x * y) + z)))));
double t_2 = i + (y * (c + (y * ((y * (y + a)) + b))));
double t_3 = (t_1 + t) / t_2;
double t_4 = x * t_2;
double tmp;
if (t_3 <= -Double.POSITIVE_INFINITY) {
tmp = x * ((t / t_4) + (((y * (230661.510616 + (y * 27464.7644705))) / t_4) + (Math.pow(y, 4.0) / t_2)));
} else if (t_3 <= 5e+243) {
tmp = (t / t_2) + (t_1 / t_2);
} else {
tmp = x + ((z / y) - (a * (x / y)));
}
return tmp;
}
def code(x, y, z, t, a, b, c, i): t_1 = y * (230661.510616 + (y * (27464.7644705 + (y * ((x * y) + z))))) t_2 = i + (y * (c + (y * ((y * (y + a)) + b)))) t_3 = (t_1 + t) / t_2 t_4 = x * t_2 tmp = 0 if t_3 <= -math.inf: tmp = x * ((t / t_4) + (((y * (230661.510616 + (y * 27464.7644705))) / t_4) + (math.pow(y, 4.0) / t_2))) elif t_3 <= 5e+243: tmp = (t / t_2) + (t_1 / t_2) else: tmp = x + ((z / y) - (a * (x / y))) return tmp
function code(x, y, z, t, a, b, c, i) t_1 = Float64(y * Float64(230661.510616 + Float64(y * Float64(27464.7644705 + Float64(y * Float64(Float64(x * y) + z)))))) t_2 = Float64(i + Float64(y * Float64(c + Float64(y * Float64(Float64(y * Float64(y + a)) + b))))) t_3 = Float64(Float64(t_1 + t) / t_2) t_4 = Float64(x * t_2) tmp = 0.0 if (t_3 <= Float64(-Inf)) tmp = Float64(x * Float64(Float64(t / t_4) + Float64(Float64(Float64(y * Float64(230661.510616 + Float64(y * 27464.7644705))) / t_4) + Float64((y ^ 4.0) / t_2)))); elseif (t_3 <= 5e+243) tmp = Float64(Float64(t / t_2) + Float64(t_1 / t_2)); else tmp = Float64(x + Float64(Float64(z / y) - Float64(a * Float64(x / y)))); end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c, i) t_1 = y * (230661.510616 + (y * (27464.7644705 + (y * ((x * y) + z))))); t_2 = i + (y * (c + (y * ((y * (y + a)) + b)))); t_3 = (t_1 + t) / t_2; t_4 = x * t_2; tmp = 0.0; if (t_3 <= -Inf) tmp = x * ((t / t_4) + (((y * (230661.510616 + (y * 27464.7644705))) / t_4) + ((y ^ 4.0) / t_2))); elseif (t_3 <= 5e+243) tmp = (t / t_2) + (t_1 / t_2); else tmp = x + ((z / y) - (a * (x / y))); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_, i_] := Block[{t$95$1 = N[(y * N[(230661.510616 + N[(y * N[(27464.7644705 + N[(y * N[(N[(x * y), $MachinePrecision] + z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(i + N[(y * N[(c + N[(y * N[(N[(y * N[(y + a), $MachinePrecision]), $MachinePrecision] + b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$3 = N[(N[(t$95$1 + t), $MachinePrecision] / t$95$2), $MachinePrecision]}, Block[{t$95$4 = N[(x * t$95$2), $MachinePrecision]}, If[LessEqual[t$95$3, (-Infinity)], N[(x * N[(N[(t / t$95$4), $MachinePrecision] + N[(N[(N[(y * N[(230661.510616 + N[(y * 27464.7644705), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / t$95$4), $MachinePrecision] + N[(N[Power[y, 4.0], $MachinePrecision] / t$95$2), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$3, 5e+243], N[(N[(t / t$95$2), $MachinePrecision] + N[(t$95$1 / t$95$2), $MachinePrecision]), $MachinePrecision], N[(x + N[(N[(z / y), $MachinePrecision] - N[(a * N[(x / y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := y \cdot \left(230661.510616 + y \cdot \left(27464.7644705 + y \cdot \left(x \cdot y + z\right)\right)\right)\\
t_2 := i + y \cdot \left(c + y \cdot \left(y \cdot \left(y + a\right) + b\right)\right)\\
t_3 := \frac{t\_1 + t}{t\_2}\\
t_4 := x \cdot t\_2\\
\mathbf{if}\;t\_3 \leq -\infty:\\
\;\;\;\;x \cdot \left(\frac{t}{t\_4} + \left(\frac{y \cdot \left(230661.510616 + y \cdot 27464.7644705\right)}{t\_4} + \frac{{y}^{4}}{t\_2}\right)\right)\\
\mathbf{elif}\;t\_3 \leq 5 \cdot 10^{+243}:\\
\;\;\;\;\frac{t}{t\_2} + \frac{t\_1}{t\_2}\\
\mathbf{else}:\\
\;\;\;\;x + \left(\frac{z}{y} - a \cdot \frac{x}{y}\right)\\
\end{array}
\end{array}
if (/.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 x y) z) y) #s(literal 54929528941/2000000 binary64)) y) #s(literal 28832688827/125000 binary64)) y) t) (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 y a) y) b) y) c) y) i)) < -inf.0Initial program 41.7%
fma-define41.7%
fma-define41.7%
fma-define41.7%
fma-define41.7%
fma-define41.7%
fma-define41.7%
fma-define41.7%
Simplified41.7%
Taylor expanded in x around inf 41.1%
Taylor expanded in z around 0 41.7%
Taylor expanded in x around inf 92.0%
if -inf.0 < (/.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 x y) z) y) #s(literal 54929528941/2000000 binary64)) y) #s(literal 28832688827/125000 binary64)) y) t) (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 y a) y) b) y) c) y) i)) < 5.00000000000000037e243Initial program 94.5%
Taylor expanded in t around 0 94.6%
if 5.00000000000000037e243 < (/.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 x y) z) y) #s(literal 54929528941/2000000 binary64)) y) #s(literal 28832688827/125000 binary64)) y) t) (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 y a) y) b) y) c) y) i)) Initial program 4.4%
Taylor expanded in y around inf 67.7%
associate--l+67.7%
associate-/l*76.3%
Simplified76.3%
Final simplification87.7%
(FPCore (x y z t a b c i)
:precision binary64
(let* ((t_1 (+ i (* y (+ c (* y (+ (* y (+ y a)) b))))))
(t_2 (/ t t_1))
(t_3
(* y (+ 230661.510616 (* y (+ 27464.7644705 (* y (+ (* x y) z)))))))
(t_4 (/ (+ t_3 t) t_1)))
(if (<= t_4 (- INFINITY))
(+ t_2 (* y (/ x y)))
(if (<= t_4 5e+243)
(+ t_2 (/ t_3 t_1))
(+ x (- (/ z y) (* a (/ x y))))))))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double t_1 = i + (y * (c + (y * ((y * (y + a)) + b))));
double t_2 = t / t_1;
double t_3 = y * (230661.510616 + (y * (27464.7644705 + (y * ((x * y) + z)))));
double t_4 = (t_3 + t) / t_1;
double tmp;
if (t_4 <= -((double) INFINITY)) {
tmp = t_2 + (y * (x / y));
} else if (t_4 <= 5e+243) {
tmp = t_2 + (t_3 / t_1);
} else {
tmp = x + ((z / y) - (a * (x / y)));
}
return tmp;
}
public static double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double t_1 = i + (y * (c + (y * ((y * (y + a)) + b))));
double t_2 = t / t_1;
double t_3 = y * (230661.510616 + (y * (27464.7644705 + (y * ((x * y) + z)))));
double t_4 = (t_3 + t) / t_1;
double tmp;
if (t_4 <= -Double.POSITIVE_INFINITY) {
tmp = t_2 + (y * (x / y));
} else if (t_4 <= 5e+243) {
tmp = t_2 + (t_3 / t_1);
} else {
tmp = x + ((z / y) - (a * (x / y)));
}
return tmp;
}
def code(x, y, z, t, a, b, c, i): t_1 = i + (y * (c + (y * ((y * (y + a)) + b)))) t_2 = t / t_1 t_3 = y * (230661.510616 + (y * (27464.7644705 + (y * ((x * y) + z))))) t_4 = (t_3 + t) / t_1 tmp = 0 if t_4 <= -math.inf: tmp = t_2 + (y * (x / y)) elif t_4 <= 5e+243: tmp = t_2 + (t_3 / t_1) else: tmp = x + ((z / y) - (a * (x / y))) return tmp
function code(x, y, z, t, a, b, c, i) t_1 = Float64(i + Float64(y * Float64(c + Float64(y * Float64(Float64(y * Float64(y + a)) + b))))) t_2 = Float64(t / t_1) t_3 = Float64(y * Float64(230661.510616 + Float64(y * Float64(27464.7644705 + Float64(y * Float64(Float64(x * y) + z)))))) t_4 = Float64(Float64(t_3 + t) / t_1) tmp = 0.0 if (t_4 <= Float64(-Inf)) tmp = Float64(t_2 + Float64(y * Float64(x / y))); elseif (t_4 <= 5e+243) tmp = Float64(t_2 + Float64(t_3 / t_1)); else tmp = Float64(x + Float64(Float64(z / y) - Float64(a * Float64(x / y)))); end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c, i) t_1 = i + (y * (c + (y * ((y * (y + a)) + b)))); t_2 = t / t_1; t_3 = y * (230661.510616 + (y * (27464.7644705 + (y * ((x * y) + z))))); t_4 = (t_3 + t) / t_1; tmp = 0.0; if (t_4 <= -Inf) tmp = t_2 + (y * (x / y)); elseif (t_4 <= 5e+243) tmp = t_2 + (t_3 / t_1); else tmp = x + ((z / y) - (a * (x / y))); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_, i_] := Block[{t$95$1 = N[(i + N[(y * N[(c + N[(y * N[(N[(y * N[(y + a), $MachinePrecision]), $MachinePrecision] + b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(t / t$95$1), $MachinePrecision]}, Block[{t$95$3 = N[(y * N[(230661.510616 + N[(y * N[(27464.7644705 + N[(y * N[(N[(x * y), $MachinePrecision] + z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$4 = N[(N[(t$95$3 + t), $MachinePrecision] / t$95$1), $MachinePrecision]}, If[LessEqual[t$95$4, (-Infinity)], N[(t$95$2 + N[(y * N[(x / y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$4, 5e+243], N[(t$95$2 + N[(t$95$3 / t$95$1), $MachinePrecision]), $MachinePrecision], N[(x + N[(N[(z / y), $MachinePrecision] - N[(a * N[(x / y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := i + y \cdot \left(c + y \cdot \left(y \cdot \left(y + a\right) + b\right)\right)\\
t_2 := \frac{t}{t\_1}\\
t_3 := y \cdot \left(230661.510616 + y \cdot \left(27464.7644705 + y \cdot \left(x \cdot y + z\right)\right)\right)\\
t_4 := \frac{t\_3 + t}{t\_1}\\
\mathbf{if}\;t\_4 \leq -\infty:\\
\;\;\;\;t\_2 + y \cdot \frac{x}{y}\\
\mathbf{elif}\;t\_4 \leq 5 \cdot 10^{+243}:\\
\;\;\;\;t\_2 + \frac{t\_3}{t\_1}\\
\mathbf{else}:\\
\;\;\;\;x + \left(\frac{z}{y} - a \cdot \frac{x}{y}\right)\\
\end{array}
\end{array}
if (/.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 x y) z) y) #s(literal 54929528941/2000000 binary64)) y) #s(literal 28832688827/125000 binary64)) y) t) (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 y a) y) b) y) c) y) i)) < -inf.0Initial program 41.7%
Taylor expanded in t around 0 41.7%
Taylor expanded in a around 0 41.7%
associate-/l*56.5%
distribute-lft-in56.5%
unpow256.5%
cube-mult56.5%
Simplified56.5%
Taylor expanded in y around inf 79.5%
if -inf.0 < (/.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 x y) z) y) #s(literal 54929528941/2000000 binary64)) y) #s(literal 28832688827/125000 binary64)) y) t) (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 y a) y) b) y) c) y) i)) < 5.00000000000000037e243Initial program 94.5%
Taylor expanded in t around 0 94.6%
if 5.00000000000000037e243 < (/.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 x y) z) y) #s(literal 54929528941/2000000 binary64)) y) #s(literal 28832688827/125000 binary64)) y) t) (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 y a) y) b) y) c) y) i)) Initial program 4.4%
Taylor expanded in y around inf 67.7%
associate--l+67.7%
associate-/l*76.3%
Simplified76.3%
Final simplification87.0%
(FPCore (x y z t a b c i)
:precision binary64
(let* ((t_1 (+ i (* y (+ c (* y (+ (* y (+ y a)) b))))))
(t_2
(/
(+
(* y (+ 230661.510616 (* y (+ 27464.7644705 (* y (+ (* x y) z))))))
t)
t_1)))
(if (<= t_2 (- INFINITY))
(+ (/ t t_1) (* y (/ x y)))
(if (<= t_2 5e+243) t_2 (+ x (- (/ z y) (* a (/ x y))))))))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double t_1 = i + (y * (c + (y * ((y * (y + a)) + b))));
double t_2 = ((y * (230661.510616 + (y * (27464.7644705 + (y * ((x * y) + z)))))) + t) / t_1;
double tmp;
if (t_2 <= -((double) INFINITY)) {
tmp = (t / t_1) + (y * (x / y));
} else if (t_2 <= 5e+243) {
tmp = t_2;
} else {
tmp = x + ((z / y) - (a * (x / y)));
}
return tmp;
}
public static double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double t_1 = i + (y * (c + (y * ((y * (y + a)) + b))));
double t_2 = ((y * (230661.510616 + (y * (27464.7644705 + (y * ((x * y) + z)))))) + t) / t_1;
double tmp;
if (t_2 <= -Double.POSITIVE_INFINITY) {
tmp = (t / t_1) + (y * (x / y));
} else if (t_2 <= 5e+243) {
tmp = t_2;
} else {
tmp = x + ((z / y) - (a * (x / y)));
}
return tmp;
}
def code(x, y, z, t, a, b, c, i): t_1 = i + (y * (c + (y * ((y * (y + a)) + b)))) t_2 = ((y * (230661.510616 + (y * (27464.7644705 + (y * ((x * y) + z)))))) + t) / t_1 tmp = 0 if t_2 <= -math.inf: tmp = (t / t_1) + (y * (x / y)) elif t_2 <= 5e+243: tmp = t_2 else: tmp = x + ((z / y) - (a * (x / y))) return tmp
function code(x, y, z, t, a, b, c, i) t_1 = Float64(i + Float64(y * Float64(c + Float64(y * Float64(Float64(y * Float64(y + a)) + b))))) t_2 = Float64(Float64(Float64(y * Float64(230661.510616 + Float64(y * Float64(27464.7644705 + Float64(y * Float64(Float64(x * y) + z)))))) + t) / t_1) tmp = 0.0 if (t_2 <= Float64(-Inf)) tmp = Float64(Float64(t / t_1) + Float64(y * Float64(x / y))); elseif (t_2 <= 5e+243) tmp = t_2; else tmp = Float64(x + Float64(Float64(z / y) - Float64(a * Float64(x / y)))); end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c, i) t_1 = i + (y * (c + (y * ((y * (y + a)) + b)))); t_2 = ((y * (230661.510616 + (y * (27464.7644705 + (y * ((x * y) + z)))))) + t) / t_1; tmp = 0.0; if (t_2 <= -Inf) tmp = (t / t_1) + (y * (x / y)); elseif (t_2 <= 5e+243) tmp = t_2; else tmp = x + ((z / y) - (a * (x / y))); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_, i_] := Block[{t$95$1 = N[(i + N[(y * N[(c + N[(y * N[(N[(y * N[(y + a), $MachinePrecision]), $MachinePrecision] + b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(N[(N[(y * N[(230661.510616 + N[(y * N[(27464.7644705 + N[(y * N[(N[(x * y), $MachinePrecision] + z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + t), $MachinePrecision] / t$95$1), $MachinePrecision]}, If[LessEqual[t$95$2, (-Infinity)], N[(N[(t / t$95$1), $MachinePrecision] + N[(y * N[(x / y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$2, 5e+243], t$95$2, N[(x + N[(N[(z / y), $MachinePrecision] - N[(a * N[(x / y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := i + y \cdot \left(c + y \cdot \left(y \cdot \left(y + a\right) + b\right)\right)\\
t_2 := \frac{y \cdot \left(230661.510616 + y \cdot \left(27464.7644705 + y \cdot \left(x \cdot y + z\right)\right)\right) + t}{t\_1}\\
\mathbf{if}\;t\_2 \leq -\infty:\\
\;\;\;\;\frac{t}{t\_1} + y \cdot \frac{x}{y}\\
\mathbf{elif}\;t\_2 \leq 5 \cdot 10^{+243}:\\
\;\;\;\;t\_2\\
\mathbf{else}:\\
\;\;\;\;x + \left(\frac{z}{y} - a \cdot \frac{x}{y}\right)\\
\end{array}
\end{array}
if (/.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 x y) z) y) #s(literal 54929528941/2000000 binary64)) y) #s(literal 28832688827/125000 binary64)) y) t) (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 y a) y) b) y) c) y) i)) < -inf.0Initial program 41.7%
Taylor expanded in t around 0 41.7%
Taylor expanded in a around 0 41.7%
associate-/l*56.5%
distribute-lft-in56.5%
unpow256.5%
cube-mult56.5%
Simplified56.5%
Taylor expanded in y around inf 79.5%
if -inf.0 < (/.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 x y) z) y) #s(literal 54929528941/2000000 binary64)) y) #s(literal 28832688827/125000 binary64)) y) t) (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 y a) y) b) y) c) y) i)) < 5.00000000000000037e243Initial program 94.5%
if 5.00000000000000037e243 < (/.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 x y) z) y) #s(literal 54929528941/2000000 binary64)) y) #s(literal 28832688827/125000 binary64)) y) t) (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 y a) y) b) y) c) y) i)) Initial program 4.4%
Taylor expanded in y around inf 67.7%
associate--l+67.7%
associate-/l*76.3%
Simplified76.3%
Final simplification87.0%
(FPCore (x y z t a b c i)
:precision binary64
(let* ((t_1 (/ t (+ i (* y c))))
(t_2 (/ t (* y (+ c (* y b)))))
(t_3 (+ x (- (/ z y) (* a (/ x y))))))
(if (<= y -9e+25)
t_3
(if (<= y -1.6e-34)
t_2
(if (<= y -1.2e-148)
t_1
(if (<= y -1.18e-256)
(+ (/ t i) (* y (/ 230661.510616 i)))
(if (<= y 2.3e-229)
t_1
(if (<= y 2.2e-86)
(/ (+ t (* y 230661.510616)) i)
(if (<= y 1.06e-78) t_2 (if (<= y 4.2e-75) (/ t i) t_3))))))))))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double t_1 = t / (i + (y * c));
double t_2 = t / (y * (c + (y * b)));
double t_3 = x + ((z / y) - (a * (x / y)));
double tmp;
if (y <= -9e+25) {
tmp = t_3;
} else if (y <= -1.6e-34) {
tmp = t_2;
} else if (y <= -1.2e-148) {
tmp = t_1;
} else if (y <= -1.18e-256) {
tmp = (t / i) + (y * (230661.510616 / i));
} else if (y <= 2.3e-229) {
tmp = t_1;
} else if (y <= 2.2e-86) {
tmp = (t + (y * 230661.510616)) / i;
} else if (y <= 1.06e-78) {
tmp = t_2;
} else if (y <= 4.2e-75) {
tmp = t / i;
} else {
tmp = t_3;
}
return tmp;
}
real(8) function code(x, y, z, t, a, b, c, i)
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), intent (in) :: i
real(8) :: t_1
real(8) :: t_2
real(8) :: t_3
real(8) :: tmp
t_1 = t / (i + (y * c))
t_2 = t / (y * (c + (y * b)))
t_3 = x + ((z / y) - (a * (x / y)))
if (y <= (-9d+25)) then
tmp = t_3
else if (y <= (-1.6d-34)) then
tmp = t_2
else if (y <= (-1.2d-148)) then
tmp = t_1
else if (y <= (-1.18d-256)) then
tmp = (t / i) + (y * (230661.510616d0 / i))
else if (y <= 2.3d-229) then
tmp = t_1
else if (y <= 2.2d-86) then
tmp = (t + (y * 230661.510616d0)) / i
else if (y <= 1.06d-78) then
tmp = t_2
else if (y <= 4.2d-75) then
tmp = t / i
else
tmp = t_3
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 i) {
double t_1 = t / (i + (y * c));
double t_2 = t / (y * (c + (y * b)));
double t_3 = x + ((z / y) - (a * (x / y)));
double tmp;
if (y <= -9e+25) {
tmp = t_3;
} else if (y <= -1.6e-34) {
tmp = t_2;
} else if (y <= -1.2e-148) {
tmp = t_1;
} else if (y <= -1.18e-256) {
tmp = (t / i) + (y * (230661.510616 / i));
} else if (y <= 2.3e-229) {
tmp = t_1;
} else if (y <= 2.2e-86) {
tmp = (t + (y * 230661.510616)) / i;
} else if (y <= 1.06e-78) {
tmp = t_2;
} else if (y <= 4.2e-75) {
tmp = t / i;
} else {
tmp = t_3;
}
return tmp;
}
def code(x, y, z, t, a, b, c, i): t_1 = t / (i + (y * c)) t_2 = t / (y * (c + (y * b))) t_3 = x + ((z / y) - (a * (x / y))) tmp = 0 if y <= -9e+25: tmp = t_3 elif y <= -1.6e-34: tmp = t_2 elif y <= -1.2e-148: tmp = t_1 elif y <= -1.18e-256: tmp = (t / i) + (y * (230661.510616 / i)) elif y <= 2.3e-229: tmp = t_1 elif y <= 2.2e-86: tmp = (t + (y * 230661.510616)) / i elif y <= 1.06e-78: tmp = t_2 elif y <= 4.2e-75: tmp = t / i else: tmp = t_3 return tmp
function code(x, y, z, t, a, b, c, i) t_1 = Float64(t / Float64(i + Float64(y * c))) t_2 = Float64(t / Float64(y * Float64(c + Float64(y * b)))) t_3 = Float64(x + Float64(Float64(z / y) - Float64(a * Float64(x / y)))) tmp = 0.0 if (y <= -9e+25) tmp = t_3; elseif (y <= -1.6e-34) tmp = t_2; elseif (y <= -1.2e-148) tmp = t_1; elseif (y <= -1.18e-256) tmp = Float64(Float64(t / i) + Float64(y * Float64(230661.510616 / i))); elseif (y <= 2.3e-229) tmp = t_1; elseif (y <= 2.2e-86) tmp = Float64(Float64(t + Float64(y * 230661.510616)) / i); elseif (y <= 1.06e-78) tmp = t_2; elseif (y <= 4.2e-75) tmp = Float64(t / i); else tmp = t_3; end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c, i) t_1 = t / (i + (y * c)); t_2 = t / (y * (c + (y * b))); t_3 = x + ((z / y) - (a * (x / y))); tmp = 0.0; if (y <= -9e+25) tmp = t_3; elseif (y <= -1.6e-34) tmp = t_2; elseif (y <= -1.2e-148) tmp = t_1; elseif (y <= -1.18e-256) tmp = (t / i) + (y * (230661.510616 / i)); elseif (y <= 2.3e-229) tmp = t_1; elseif (y <= 2.2e-86) tmp = (t + (y * 230661.510616)) / i; elseif (y <= 1.06e-78) tmp = t_2; elseif (y <= 4.2e-75) tmp = t / i; else tmp = t_3; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_, i_] := Block[{t$95$1 = N[(t / N[(i + N[(y * c), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(t / N[(y * N[(c + N[(y * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$3 = N[(x + N[(N[(z / y), $MachinePrecision] - N[(a * N[(x / y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -9e+25], t$95$3, If[LessEqual[y, -1.6e-34], t$95$2, If[LessEqual[y, -1.2e-148], t$95$1, If[LessEqual[y, -1.18e-256], N[(N[(t / i), $MachinePrecision] + N[(y * N[(230661.510616 / i), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 2.3e-229], t$95$1, If[LessEqual[y, 2.2e-86], N[(N[(t + N[(y * 230661.510616), $MachinePrecision]), $MachinePrecision] / i), $MachinePrecision], If[LessEqual[y, 1.06e-78], t$95$2, If[LessEqual[y, 4.2e-75], N[(t / i), $MachinePrecision], t$95$3]]]]]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{t}{i + y \cdot c}\\
t_2 := \frac{t}{y \cdot \left(c + y \cdot b\right)}\\
t_3 := x + \left(\frac{z}{y} - a \cdot \frac{x}{y}\right)\\
\mathbf{if}\;y \leq -9 \cdot 10^{+25}:\\
\;\;\;\;t\_3\\
\mathbf{elif}\;y \leq -1.6 \cdot 10^{-34}:\\
\;\;\;\;t\_2\\
\mathbf{elif}\;y \leq -1.2 \cdot 10^{-148}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y \leq -1.18 \cdot 10^{-256}:\\
\;\;\;\;\frac{t}{i} + y \cdot \frac{230661.510616}{i}\\
\mathbf{elif}\;y \leq 2.3 \cdot 10^{-229}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y \leq 2.2 \cdot 10^{-86}:\\
\;\;\;\;\frac{t + y \cdot 230661.510616}{i}\\
\mathbf{elif}\;y \leq 1.06 \cdot 10^{-78}:\\
\;\;\;\;t\_2\\
\mathbf{elif}\;y \leq 4.2 \cdot 10^{-75}:\\
\;\;\;\;\frac{t}{i}\\
\mathbf{else}:\\
\;\;\;\;t\_3\\
\end{array}
\end{array}
if y < -9.0000000000000006e25 or 4.2000000000000002e-75 < y Initial program 24.3%
Taylor expanded in y around inf 54.0%
associate--l+54.0%
associate-/l*59.9%
Simplified59.9%
if -9.0000000000000006e25 < y < -1.60000000000000001e-34 or 2.2000000000000002e-86 < y < 1.06e-78Initial program 93.2%
Taylor expanded in t around inf 52.2%
Taylor expanded in y around 0 45.6%
*-commutative45.6%
Simplified45.6%
Taylor expanded in i around 0 45.0%
if -1.60000000000000001e-34 < y < -1.2000000000000001e-148 or -1.18e-256 < y < 2.29999999999999996e-229Initial program 99.7%
Taylor expanded in t around inf 78.5%
Taylor expanded in y around 0 76.4%
*-commutative76.4%
Simplified76.4%
if -1.2000000000000001e-148 < y < -1.18e-256Initial program 99.9%
Taylor expanded in t around 0 99.9%
Taylor expanded in a around 0 99.9%
associate-/l*99.9%
distribute-lft-in99.9%
unpow299.9%
cube-mult99.9%
Simplified99.9%
Taylor expanded in y around 0 84.9%
Taylor expanded in i around inf 80.8%
if 2.29999999999999996e-229 < y < 2.2000000000000002e-86Initial program 99.7%
Taylor expanded in y around 0 96.8%
*-commutative96.8%
Simplified96.8%
Taylor expanded in i around inf 73.9%
if 1.06e-78 < y < 4.2000000000000002e-75Initial program 100.0%
Taylor expanded in y around 0 68.6%
Final simplification65.7%
(FPCore (x y z t a b c i)
:precision binary64
(let* ((t_1 (* y (+ 230661.510616 (* y (+ 27464.7644705 (* y z))))))
(t_2 (/ (+ t t_1) (+ i (* y c))))
(t_3 (+ x (- (/ z y) (* a (/ x y)))))
(t_4 (* y (+ c (* y (+ (* y (+ y a)) b))))))
(if (<= y -2.3e+27)
t_3
(if (<= y -1.3e-35)
(/ (+ t (* y (+ 230661.510616 (* y 27464.7644705)))) t_4)
(if (<= y -1.5e-70)
(/ t_1 (+ i (* y (+ c (* y b)))))
(if (<= y -4.6e-125)
t_2
(if (<= y 6.5e-41)
(/ (+ t (* y 230661.510616)) (+ i t_4))
(if (<= y 1.8e-10) t_2 t_3))))))))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double t_1 = y * (230661.510616 + (y * (27464.7644705 + (y * z))));
double t_2 = (t + t_1) / (i + (y * c));
double t_3 = x + ((z / y) - (a * (x / y)));
double t_4 = y * (c + (y * ((y * (y + a)) + b)));
double tmp;
if (y <= -2.3e+27) {
tmp = t_3;
} else if (y <= -1.3e-35) {
tmp = (t + (y * (230661.510616 + (y * 27464.7644705)))) / t_4;
} else if (y <= -1.5e-70) {
tmp = t_1 / (i + (y * (c + (y * b))));
} else if (y <= -4.6e-125) {
tmp = t_2;
} else if (y <= 6.5e-41) {
tmp = (t + (y * 230661.510616)) / (i + t_4);
} else if (y <= 1.8e-10) {
tmp = t_2;
} else {
tmp = t_3;
}
return tmp;
}
real(8) function code(x, y, z, t, a, b, c, i)
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), intent (in) :: i
real(8) :: t_1
real(8) :: t_2
real(8) :: t_3
real(8) :: t_4
real(8) :: tmp
t_1 = y * (230661.510616d0 + (y * (27464.7644705d0 + (y * z))))
t_2 = (t + t_1) / (i + (y * c))
t_3 = x + ((z / y) - (a * (x / y)))
t_4 = y * (c + (y * ((y * (y + a)) + b)))
if (y <= (-2.3d+27)) then
tmp = t_3
else if (y <= (-1.3d-35)) then
tmp = (t + (y * (230661.510616d0 + (y * 27464.7644705d0)))) / t_4
else if (y <= (-1.5d-70)) then
tmp = t_1 / (i + (y * (c + (y * b))))
else if (y <= (-4.6d-125)) then
tmp = t_2
else if (y <= 6.5d-41) then
tmp = (t + (y * 230661.510616d0)) / (i + t_4)
else if (y <= 1.8d-10) then
tmp = t_2
else
tmp = t_3
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 i) {
double t_1 = y * (230661.510616 + (y * (27464.7644705 + (y * z))));
double t_2 = (t + t_1) / (i + (y * c));
double t_3 = x + ((z / y) - (a * (x / y)));
double t_4 = y * (c + (y * ((y * (y + a)) + b)));
double tmp;
if (y <= -2.3e+27) {
tmp = t_3;
} else if (y <= -1.3e-35) {
tmp = (t + (y * (230661.510616 + (y * 27464.7644705)))) / t_4;
} else if (y <= -1.5e-70) {
tmp = t_1 / (i + (y * (c + (y * b))));
} else if (y <= -4.6e-125) {
tmp = t_2;
} else if (y <= 6.5e-41) {
tmp = (t + (y * 230661.510616)) / (i + t_4);
} else if (y <= 1.8e-10) {
tmp = t_2;
} else {
tmp = t_3;
}
return tmp;
}
def code(x, y, z, t, a, b, c, i): t_1 = y * (230661.510616 + (y * (27464.7644705 + (y * z)))) t_2 = (t + t_1) / (i + (y * c)) t_3 = x + ((z / y) - (a * (x / y))) t_4 = y * (c + (y * ((y * (y + a)) + b))) tmp = 0 if y <= -2.3e+27: tmp = t_3 elif y <= -1.3e-35: tmp = (t + (y * (230661.510616 + (y * 27464.7644705)))) / t_4 elif y <= -1.5e-70: tmp = t_1 / (i + (y * (c + (y * b)))) elif y <= -4.6e-125: tmp = t_2 elif y <= 6.5e-41: tmp = (t + (y * 230661.510616)) / (i + t_4) elif y <= 1.8e-10: tmp = t_2 else: tmp = t_3 return tmp
function code(x, y, z, t, a, b, c, i) t_1 = Float64(y * Float64(230661.510616 + Float64(y * Float64(27464.7644705 + Float64(y * z))))) t_2 = Float64(Float64(t + t_1) / Float64(i + Float64(y * c))) t_3 = Float64(x + Float64(Float64(z / y) - Float64(a * Float64(x / y)))) t_4 = Float64(y * Float64(c + Float64(y * Float64(Float64(y * Float64(y + a)) + b)))) tmp = 0.0 if (y <= -2.3e+27) tmp = t_3; elseif (y <= -1.3e-35) tmp = Float64(Float64(t + Float64(y * Float64(230661.510616 + Float64(y * 27464.7644705)))) / t_4); elseif (y <= -1.5e-70) tmp = Float64(t_1 / Float64(i + Float64(y * Float64(c + Float64(y * b))))); elseif (y <= -4.6e-125) tmp = t_2; elseif (y <= 6.5e-41) tmp = Float64(Float64(t + Float64(y * 230661.510616)) / Float64(i + t_4)); elseif (y <= 1.8e-10) tmp = t_2; else tmp = t_3; end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c, i) t_1 = y * (230661.510616 + (y * (27464.7644705 + (y * z)))); t_2 = (t + t_1) / (i + (y * c)); t_3 = x + ((z / y) - (a * (x / y))); t_4 = y * (c + (y * ((y * (y + a)) + b))); tmp = 0.0; if (y <= -2.3e+27) tmp = t_3; elseif (y <= -1.3e-35) tmp = (t + (y * (230661.510616 + (y * 27464.7644705)))) / t_4; elseif (y <= -1.5e-70) tmp = t_1 / (i + (y * (c + (y * b)))); elseif (y <= -4.6e-125) tmp = t_2; elseif (y <= 6.5e-41) tmp = (t + (y * 230661.510616)) / (i + t_4); elseif (y <= 1.8e-10) tmp = t_2; else tmp = t_3; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_, i_] := Block[{t$95$1 = N[(y * N[(230661.510616 + N[(y * N[(27464.7644705 + N[(y * z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(N[(t + t$95$1), $MachinePrecision] / N[(i + N[(y * c), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$3 = N[(x + N[(N[(z / y), $MachinePrecision] - N[(a * N[(x / y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$4 = N[(y * N[(c + N[(y * N[(N[(y * N[(y + a), $MachinePrecision]), $MachinePrecision] + b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -2.3e+27], t$95$3, If[LessEqual[y, -1.3e-35], N[(N[(t + N[(y * N[(230661.510616 + N[(y * 27464.7644705), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / t$95$4), $MachinePrecision], If[LessEqual[y, -1.5e-70], N[(t$95$1 / N[(i + N[(y * N[(c + N[(y * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, -4.6e-125], t$95$2, If[LessEqual[y, 6.5e-41], N[(N[(t + N[(y * 230661.510616), $MachinePrecision]), $MachinePrecision] / N[(i + t$95$4), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 1.8e-10], t$95$2, t$95$3]]]]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := y \cdot \left(230661.510616 + y \cdot \left(27464.7644705 + y \cdot z\right)\right)\\
t_2 := \frac{t + t\_1}{i + y \cdot c}\\
t_3 := x + \left(\frac{z}{y} - a \cdot \frac{x}{y}\right)\\
t_4 := y \cdot \left(c + y \cdot \left(y \cdot \left(y + a\right) + b\right)\right)\\
\mathbf{if}\;y \leq -2.3 \cdot 10^{+27}:\\
\;\;\;\;t\_3\\
\mathbf{elif}\;y \leq -1.3 \cdot 10^{-35}:\\
\;\;\;\;\frac{t + y \cdot \left(230661.510616 + y \cdot 27464.7644705\right)}{t\_4}\\
\mathbf{elif}\;y \leq -1.5 \cdot 10^{-70}:\\
\;\;\;\;\frac{t\_1}{i + y \cdot \left(c + y \cdot b\right)}\\
\mathbf{elif}\;y \leq -4.6 \cdot 10^{-125}:\\
\;\;\;\;t\_2\\
\mathbf{elif}\;y \leq 6.5 \cdot 10^{-41}:\\
\;\;\;\;\frac{t + y \cdot 230661.510616}{i + t\_4}\\
\mathbf{elif}\;y \leq 1.8 \cdot 10^{-10}:\\
\;\;\;\;t\_2\\
\mathbf{else}:\\
\;\;\;\;t\_3\\
\end{array}
\end{array}
if y < -2.3000000000000001e27 or 1.8e-10 < y Initial program 10.1%
Taylor expanded in y around inf 62.7%
associate--l+62.7%
associate-/l*69.8%
Simplified69.8%
if -2.3000000000000001e27 < y < -1.30000000000000002e-35Initial program 92.3%
Taylor expanded in i around 0 85.2%
Taylor expanded in y around 0 74.6%
if -1.30000000000000002e-35 < y < -1.5000000000000001e-70Initial program 98.9%
Taylor expanded in x around 0 98.9%
Taylor expanded in y around 0 83.3%
Taylor expanded in t around 0 66.6%
if -1.5000000000000001e-70 < y < -4.5999999999999998e-125 or 6.5000000000000004e-41 < y < 1.8e-10Initial program 99.7%
Taylor expanded in x around 0 91.3%
Taylor expanded in y around 0 88.9%
Taylor expanded in b around 0 80.6%
*-commutative80.6%
Simplified80.6%
if -4.5999999999999998e-125 < y < 6.5000000000000004e-41Initial program 99.8%
Taylor expanded in y around 0 97.0%
*-commutative97.0%
Simplified97.0%
Final simplification81.2%
(FPCore (x y z t a b c i)
:precision binary64
(let* ((t_1 (/ t (+ i (* y c)))) (t_2 (+ x (/ (- z (* x a)) y))))
(if (<= y -3.2e+25)
t_2
(if (<= y -5.4e-95)
t_1
(if (<= y -3.4e-125)
(/ (+ t (* y 230661.510616)) i)
(if (<= y -1.4e-235)
t_1
(if (<= y -5.7e-258)
(+ (/ t i) (* y (/ 230661.510616 i)))
(if (<= y 0.00036) t_1 (if (<= y 1e+265) t_2 x)))))))))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double t_1 = t / (i + (y * c));
double t_2 = x + ((z - (x * a)) / y);
double tmp;
if (y <= -3.2e+25) {
tmp = t_2;
} else if (y <= -5.4e-95) {
tmp = t_1;
} else if (y <= -3.4e-125) {
tmp = (t + (y * 230661.510616)) / i;
} else if (y <= -1.4e-235) {
tmp = t_1;
} else if (y <= -5.7e-258) {
tmp = (t / i) + (y * (230661.510616 / i));
} else if (y <= 0.00036) {
tmp = t_1;
} else if (y <= 1e+265) {
tmp = t_2;
} else {
tmp = x;
}
return tmp;
}
real(8) function code(x, y, z, t, a, b, c, i)
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), intent (in) :: i
real(8) :: t_1
real(8) :: t_2
real(8) :: tmp
t_1 = t / (i + (y * c))
t_2 = x + ((z - (x * a)) / y)
if (y <= (-3.2d+25)) then
tmp = t_2
else if (y <= (-5.4d-95)) then
tmp = t_1
else if (y <= (-3.4d-125)) then
tmp = (t + (y * 230661.510616d0)) / i
else if (y <= (-1.4d-235)) then
tmp = t_1
else if (y <= (-5.7d-258)) then
tmp = (t / i) + (y * (230661.510616d0 / i))
else if (y <= 0.00036d0) then
tmp = t_1
else if (y <= 1d+265) then
tmp = t_2
else
tmp = 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 i) {
double t_1 = t / (i + (y * c));
double t_2 = x + ((z - (x * a)) / y);
double tmp;
if (y <= -3.2e+25) {
tmp = t_2;
} else if (y <= -5.4e-95) {
tmp = t_1;
} else if (y <= -3.4e-125) {
tmp = (t + (y * 230661.510616)) / i;
} else if (y <= -1.4e-235) {
tmp = t_1;
} else if (y <= -5.7e-258) {
tmp = (t / i) + (y * (230661.510616 / i));
} else if (y <= 0.00036) {
tmp = t_1;
} else if (y <= 1e+265) {
tmp = t_2;
} else {
tmp = x;
}
return tmp;
}
def code(x, y, z, t, a, b, c, i): t_1 = t / (i + (y * c)) t_2 = x + ((z - (x * a)) / y) tmp = 0 if y <= -3.2e+25: tmp = t_2 elif y <= -5.4e-95: tmp = t_1 elif y <= -3.4e-125: tmp = (t + (y * 230661.510616)) / i elif y <= -1.4e-235: tmp = t_1 elif y <= -5.7e-258: tmp = (t / i) + (y * (230661.510616 / i)) elif y <= 0.00036: tmp = t_1 elif y <= 1e+265: tmp = t_2 else: tmp = x return tmp
function code(x, y, z, t, a, b, c, i) t_1 = Float64(t / Float64(i + Float64(y * c))) t_2 = Float64(x + Float64(Float64(z - Float64(x * a)) / y)) tmp = 0.0 if (y <= -3.2e+25) tmp = t_2; elseif (y <= -5.4e-95) tmp = t_1; elseif (y <= -3.4e-125) tmp = Float64(Float64(t + Float64(y * 230661.510616)) / i); elseif (y <= -1.4e-235) tmp = t_1; elseif (y <= -5.7e-258) tmp = Float64(Float64(t / i) + Float64(y * Float64(230661.510616 / i))); elseif (y <= 0.00036) tmp = t_1; elseif (y <= 1e+265) tmp = t_2; else tmp = x; end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c, i) t_1 = t / (i + (y * c)); t_2 = x + ((z - (x * a)) / y); tmp = 0.0; if (y <= -3.2e+25) tmp = t_2; elseif (y <= -5.4e-95) tmp = t_1; elseif (y <= -3.4e-125) tmp = (t + (y * 230661.510616)) / i; elseif (y <= -1.4e-235) tmp = t_1; elseif (y <= -5.7e-258) tmp = (t / i) + (y * (230661.510616 / i)); elseif (y <= 0.00036) tmp = t_1; elseif (y <= 1e+265) tmp = t_2; else tmp = x; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_, i_] := Block[{t$95$1 = N[(t / N[(i + N[(y * c), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(x + N[(N[(z - N[(x * a), $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -3.2e+25], t$95$2, If[LessEqual[y, -5.4e-95], t$95$1, If[LessEqual[y, -3.4e-125], N[(N[(t + N[(y * 230661.510616), $MachinePrecision]), $MachinePrecision] / i), $MachinePrecision], If[LessEqual[y, -1.4e-235], t$95$1, If[LessEqual[y, -5.7e-258], N[(N[(t / i), $MachinePrecision] + N[(y * N[(230661.510616 / i), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 0.00036], t$95$1, If[LessEqual[y, 1e+265], t$95$2, x]]]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{t}{i + y \cdot c}\\
t_2 := x + \frac{z - x \cdot a}{y}\\
\mathbf{if}\;y \leq -3.2 \cdot 10^{+25}:\\
\;\;\;\;t\_2\\
\mathbf{elif}\;y \leq -5.4 \cdot 10^{-95}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y \leq -3.4 \cdot 10^{-125}:\\
\;\;\;\;\frac{t + y \cdot 230661.510616}{i}\\
\mathbf{elif}\;y \leq -1.4 \cdot 10^{-235}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y \leq -5.7 \cdot 10^{-258}:\\
\;\;\;\;\frac{t}{i} + y \cdot \frac{230661.510616}{i}\\
\mathbf{elif}\;y \leq 0.00036:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y \leq 10^{+265}:\\
\;\;\;\;t\_2\\
\mathbf{else}:\\
\;\;\;\;x\\
\end{array}
\end{array}
if y < -3.1999999999999999e25 or 3.60000000000000023e-4 < y < 1.00000000000000007e265Initial program 9.3%
Taylor expanded in i around 0 6.6%
Taylor expanded in y around inf 62.5%
associate--l+62.5%
div-sub62.5%
Simplified62.5%
if -3.1999999999999999e25 < y < -5.4e-95 or -3.39999999999999975e-125 < y < -1.39999999999999998e-235 or -5.7000000000000002e-258 < y < 3.60000000000000023e-4Initial program 98.9%
Taylor expanded in t around inf 67.9%
Taylor expanded in y around 0 55.9%
*-commutative55.9%
Simplified55.9%
if -5.4e-95 < y < -3.39999999999999975e-125Initial program 99.6%
Taylor expanded in y around 0 77.0%
*-commutative77.0%
Simplified77.0%
Taylor expanded in i around inf 65.4%
if -1.39999999999999998e-235 < y < -5.7000000000000002e-258Initial program 99.7%
Taylor expanded in t around 0 99.7%
Taylor expanded in a around 0 99.7%
associate-/l*99.8%
distribute-lft-in99.8%
unpow299.8%
cube-mult99.8%
Simplified99.8%
Taylor expanded in y around 0 99.8%
Taylor expanded in i around inf 98.0%
if 1.00000000000000007e265 < y Initial program 0.0%
Taylor expanded in y around inf 99.0%
Final simplification62.0%
(FPCore (x y z t a b c i)
:precision binary64
(let* ((t_1 (/ t (+ i (* y c))))
(t_2 (+ x (/ (- z (* x a)) y)))
(t_3 (/ (+ t (* y 230661.510616)) i)))
(if (<= y -2.7e+25)
t_2
(if (<= y -1.75e-91)
t_1
(if (<= y -3.1e-125)
t_3
(if (<= y -9.8e-236)
t_1
(if (<= y -5.1e-258)
t_3
(if (<= y 0.000115) t_1 (if (<= y 5.8e+264) t_2 x)))))))))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double t_1 = t / (i + (y * c));
double t_2 = x + ((z - (x * a)) / y);
double t_3 = (t + (y * 230661.510616)) / i;
double tmp;
if (y <= -2.7e+25) {
tmp = t_2;
} else if (y <= -1.75e-91) {
tmp = t_1;
} else if (y <= -3.1e-125) {
tmp = t_3;
} else if (y <= -9.8e-236) {
tmp = t_1;
} else if (y <= -5.1e-258) {
tmp = t_3;
} else if (y <= 0.000115) {
tmp = t_1;
} else if (y <= 5.8e+264) {
tmp = t_2;
} else {
tmp = x;
}
return tmp;
}
real(8) function code(x, y, z, t, a, b, c, i)
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), intent (in) :: i
real(8) :: t_1
real(8) :: t_2
real(8) :: t_3
real(8) :: tmp
t_1 = t / (i + (y * c))
t_2 = x + ((z - (x * a)) / y)
t_3 = (t + (y * 230661.510616d0)) / i
if (y <= (-2.7d+25)) then
tmp = t_2
else if (y <= (-1.75d-91)) then
tmp = t_1
else if (y <= (-3.1d-125)) then
tmp = t_3
else if (y <= (-9.8d-236)) then
tmp = t_1
else if (y <= (-5.1d-258)) then
tmp = t_3
else if (y <= 0.000115d0) then
tmp = t_1
else if (y <= 5.8d+264) then
tmp = t_2
else
tmp = 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 i) {
double t_1 = t / (i + (y * c));
double t_2 = x + ((z - (x * a)) / y);
double t_3 = (t + (y * 230661.510616)) / i;
double tmp;
if (y <= -2.7e+25) {
tmp = t_2;
} else if (y <= -1.75e-91) {
tmp = t_1;
} else if (y <= -3.1e-125) {
tmp = t_3;
} else if (y <= -9.8e-236) {
tmp = t_1;
} else if (y <= -5.1e-258) {
tmp = t_3;
} else if (y <= 0.000115) {
tmp = t_1;
} else if (y <= 5.8e+264) {
tmp = t_2;
} else {
tmp = x;
}
return tmp;
}
def code(x, y, z, t, a, b, c, i): t_1 = t / (i + (y * c)) t_2 = x + ((z - (x * a)) / y) t_3 = (t + (y * 230661.510616)) / i tmp = 0 if y <= -2.7e+25: tmp = t_2 elif y <= -1.75e-91: tmp = t_1 elif y <= -3.1e-125: tmp = t_3 elif y <= -9.8e-236: tmp = t_1 elif y <= -5.1e-258: tmp = t_3 elif y <= 0.000115: tmp = t_1 elif y <= 5.8e+264: tmp = t_2 else: tmp = x return tmp
function code(x, y, z, t, a, b, c, i) t_1 = Float64(t / Float64(i + Float64(y * c))) t_2 = Float64(x + Float64(Float64(z - Float64(x * a)) / y)) t_3 = Float64(Float64(t + Float64(y * 230661.510616)) / i) tmp = 0.0 if (y <= -2.7e+25) tmp = t_2; elseif (y <= -1.75e-91) tmp = t_1; elseif (y <= -3.1e-125) tmp = t_3; elseif (y <= -9.8e-236) tmp = t_1; elseif (y <= -5.1e-258) tmp = t_3; elseif (y <= 0.000115) tmp = t_1; elseif (y <= 5.8e+264) tmp = t_2; else tmp = x; end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c, i) t_1 = t / (i + (y * c)); t_2 = x + ((z - (x * a)) / y); t_3 = (t + (y * 230661.510616)) / i; tmp = 0.0; if (y <= -2.7e+25) tmp = t_2; elseif (y <= -1.75e-91) tmp = t_1; elseif (y <= -3.1e-125) tmp = t_3; elseif (y <= -9.8e-236) tmp = t_1; elseif (y <= -5.1e-258) tmp = t_3; elseif (y <= 0.000115) tmp = t_1; elseif (y <= 5.8e+264) tmp = t_2; else tmp = x; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_, i_] := Block[{t$95$1 = N[(t / N[(i + N[(y * c), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(x + N[(N[(z - N[(x * a), $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$3 = N[(N[(t + N[(y * 230661.510616), $MachinePrecision]), $MachinePrecision] / i), $MachinePrecision]}, If[LessEqual[y, -2.7e+25], t$95$2, If[LessEqual[y, -1.75e-91], t$95$1, If[LessEqual[y, -3.1e-125], t$95$3, If[LessEqual[y, -9.8e-236], t$95$1, If[LessEqual[y, -5.1e-258], t$95$3, If[LessEqual[y, 0.000115], t$95$1, If[LessEqual[y, 5.8e+264], t$95$2, x]]]]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{t}{i + y \cdot c}\\
t_2 := x + \frac{z - x \cdot a}{y}\\
t_3 := \frac{t + y \cdot 230661.510616}{i}\\
\mathbf{if}\;y \leq -2.7 \cdot 10^{+25}:\\
\;\;\;\;t\_2\\
\mathbf{elif}\;y \leq -1.75 \cdot 10^{-91}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y \leq -3.1 \cdot 10^{-125}:\\
\;\;\;\;t\_3\\
\mathbf{elif}\;y \leq -9.8 \cdot 10^{-236}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y \leq -5.1 \cdot 10^{-258}:\\
\;\;\;\;t\_3\\
\mathbf{elif}\;y \leq 0.000115:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y \leq 5.8 \cdot 10^{+264}:\\
\;\;\;\;t\_2\\
\mathbf{else}:\\
\;\;\;\;x\\
\end{array}
\end{array}
if y < -2.7e25 or 1.15e-4 < y < 5.7999999999999996e264Initial program 9.3%
Taylor expanded in i around 0 6.6%
Taylor expanded in y around inf 62.5%
associate--l+62.5%
div-sub62.5%
Simplified62.5%
if -2.7e25 < y < -1.7499999999999999e-91 or -3.10000000000000013e-125 < y < -9.7999999999999993e-236 or -5.0999999999999997e-258 < y < 1.15e-4Initial program 98.9%
Taylor expanded in t around inf 67.9%
Taylor expanded in y around 0 55.9%
*-commutative55.9%
Simplified55.9%
if -1.7499999999999999e-91 < y < -3.10000000000000013e-125 or -9.7999999999999993e-236 < y < -5.0999999999999997e-258Initial program 99.6%
Taylor expanded in y around 0 89.0%
*-commutative89.0%
Simplified89.0%
Taylor expanded in i around inf 82.6%
if 5.7999999999999996e264 < y Initial program 0.0%
Taylor expanded in y around inf 99.0%
Final simplification62.0%
(FPCore (x y z t a b c i)
:precision binary64
(let* ((t_1 (+ x (- (/ z y) (* a (/ x y))))))
(if (<= y -7e+26)
t_1
(if (<= y -1.2e-21)
(/ t (* y (+ c (* y (+ (* y (+ y a)) b)))))
(if (<= y 1.6e-173)
(/ (* t (+ 1.0 (* 230661.510616 (/ y t)))) (+ i (* y c)))
(if (<= y 3.9e-23)
(/ t (+ i (* y (+ c (* y (+ b (* y a)))))))
(if (<= y 1.7e+69)
(/
(+ 230661.510616 (* y (+ 27464.7644705 (* y (+ (* x y) z)))))
c)
t_1)))))))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double t_1 = x + ((z / y) - (a * (x / y)));
double tmp;
if (y <= -7e+26) {
tmp = t_1;
} else if (y <= -1.2e-21) {
tmp = t / (y * (c + (y * ((y * (y + a)) + b))));
} else if (y <= 1.6e-173) {
tmp = (t * (1.0 + (230661.510616 * (y / t)))) / (i + (y * c));
} else if (y <= 3.9e-23) {
tmp = t / (i + (y * (c + (y * (b + (y * a))))));
} else if (y <= 1.7e+69) {
tmp = (230661.510616 + (y * (27464.7644705 + (y * ((x * y) + z))))) / c;
} else {
tmp = t_1;
}
return tmp;
}
real(8) function code(x, y, z, t, a, b, c, i)
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), intent (in) :: i
real(8) :: t_1
real(8) :: tmp
t_1 = x + ((z / y) - (a * (x / y)))
if (y <= (-7d+26)) then
tmp = t_1
else if (y <= (-1.2d-21)) then
tmp = t / (y * (c + (y * ((y * (y + a)) + b))))
else if (y <= 1.6d-173) then
tmp = (t * (1.0d0 + (230661.510616d0 * (y / t)))) / (i + (y * c))
else if (y <= 3.9d-23) then
tmp = t / (i + (y * (c + (y * (b + (y * a))))))
else if (y <= 1.7d+69) then
tmp = (230661.510616d0 + (y * (27464.7644705d0 + (y * ((x * y) + z))))) / c
else
tmp = t_1
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double t_1 = x + ((z / y) - (a * (x / y)));
double tmp;
if (y <= -7e+26) {
tmp = t_1;
} else if (y <= -1.2e-21) {
tmp = t / (y * (c + (y * ((y * (y + a)) + b))));
} else if (y <= 1.6e-173) {
tmp = (t * (1.0 + (230661.510616 * (y / t)))) / (i + (y * c));
} else if (y <= 3.9e-23) {
tmp = t / (i + (y * (c + (y * (b + (y * a))))));
} else if (y <= 1.7e+69) {
tmp = (230661.510616 + (y * (27464.7644705 + (y * ((x * y) + z))))) / c;
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b, c, i): t_1 = x + ((z / y) - (a * (x / y))) tmp = 0 if y <= -7e+26: tmp = t_1 elif y <= -1.2e-21: tmp = t / (y * (c + (y * ((y * (y + a)) + b)))) elif y <= 1.6e-173: tmp = (t * (1.0 + (230661.510616 * (y / t)))) / (i + (y * c)) elif y <= 3.9e-23: tmp = t / (i + (y * (c + (y * (b + (y * a)))))) elif y <= 1.7e+69: tmp = (230661.510616 + (y * (27464.7644705 + (y * ((x * y) + z))))) / c else: tmp = t_1 return tmp
function code(x, y, z, t, a, b, c, i) t_1 = Float64(x + Float64(Float64(z / y) - Float64(a * Float64(x / y)))) tmp = 0.0 if (y <= -7e+26) tmp = t_1; elseif (y <= -1.2e-21) tmp = Float64(t / Float64(y * Float64(c + Float64(y * Float64(Float64(y * Float64(y + a)) + b))))); elseif (y <= 1.6e-173) tmp = Float64(Float64(t * Float64(1.0 + Float64(230661.510616 * Float64(y / t)))) / Float64(i + Float64(y * c))); elseif (y <= 3.9e-23) tmp = Float64(t / Float64(i + Float64(y * Float64(c + Float64(y * Float64(b + Float64(y * a))))))); elseif (y <= 1.7e+69) tmp = Float64(Float64(230661.510616 + Float64(y * Float64(27464.7644705 + Float64(y * Float64(Float64(x * y) + z))))) / c); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c, i) t_1 = x + ((z / y) - (a * (x / y))); tmp = 0.0; if (y <= -7e+26) tmp = t_1; elseif (y <= -1.2e-21) tmp = t / (y * (c + (y * ((y * (y + a)) + b)))); elseif (y <= 1.6e-173) tmp = (t * (1.0 + (230661.510616 * (y / t)))) / (i + (y * c)); elseif (y <= 3.9e-23) tmp = t / (i + (y * (c + (y * (b + (y * a)))))); elseif (y <= 1.7e+69) tmp = (230661.510616 + (y * (27464.7644705 + (y * ((x * y) + z))))) / c; else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_, i_] := Block[{t$95$1 = N[(x + N[(N[(z / y), $MachinePrecision] - N[(a * N[(x / y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -7e+26], t$95$1, If[LessEqual[y, -1.2e-21], N[(t / N[(y * N[(c + N[(y * N[(N[(y * N[(y + a), $MachinePrecision]), $MachinePrecision] + b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 1.6e-173], N[(N[(t * N[(1.0 + N[(230661.510616 * N[(y / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(i + N[(y * c), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 3.9e-23], N[(t / N[(i + N[(y * N[(c + N[(y * N[(b + N[(y * a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 1.7e+69], N[(N[(230661.510616 + N[(y * N[(27464.7644705 + N[(y * N[(N[(x * y), $MachinePrecision] + z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / c), $MachinePrecision], t$95$1]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x + \left(\frac{z}{y} - a \cdot \frac{x}{y}\right)\\
\mathbf{if}\;y \leq -7 \cdot 10^{+26}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y \leq -1.2 \cdot 10^{-21}:\\
\;\;\;\;\frac{t}{y \cdot \left(c + y \cdot \left(y \cdot \left(y + a\right) + b\right)\right)}\\
\mathbf{elif}\;y \leq 1.6 \cdot 10^{-173}:\\
\;\;\;\;\frac{t \cdot \left(1 + 230661.510616 \cdot \frac{y}{t}\right)}{i + y \cdot c}\\
\mathbf{elif}\;y \leq 3.9 \cdot 10^{-23}:\\
\;\;\;\;\frac{t}{i + y \cdot \left(c + y \cdot \left(b + y \cdot a\right)\right)}\\
\mathbf{elif}\;y \leq 1.7 \cdot 10^{+69}:\\
\;\;\;\;\frac{230661.510616 + y \cdot \left(27464.7644705 + y \cdot \left(x \cdot y + z\right)\right)}{c}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if y < -6.9999999999999998e26 or 1.69999999999999993e69 < y Initial program 2.6%
Taylor expanded in y around inf 69.5%
associate--l+69.5%
associate-/l*77.6%
Simplified77.6%
if -6.9999999999999998e26 < y < -1.2e-21Initial program 90.1%
Taylor expanded in t around inf 47.9%
Taylor expanded in i around 0 47.9%
if -1.2e-21 < y < 1.6e-173Initial program 99.8%
Taylor expanded in y around 0 91.8%
*-commutative91.8%
Simplified91.8%
Taylor expanded in y around 0 87.6%
Taylor expanded in t around inf 87.6%
if 1.6e-173 < y < 3.9e-23Initial program 99.7%
Taylor expanded in t around inf 71.0%
Taylor expanded in y around 0 71.0%
*-commutative71.0%
Simplified71.0%
if 3.9e-23 < y < 1.69999999999999993e69Initial program 73.6%
Taylor expanded in t around 0 73.7%
Taylor expanded in i around 0 56.5%
Taylor expanded in c around inf 33.2%
Taylor expanded in t around 0 28.8%
Final simplification74.3%
(FPCore (x y z t a b c i)
:precision binary64
(if (or (<= y -6.2e+59) (not (<= y 1.8e-7)))
(+ x (- (/ z y) (* a (/ x y))))
(/
(+ t (* y (+ 230661.510616 (* y (+ 27464.7644705 (* y z))))))
(+ i (* y (+ c (* y (+ (* y (+ y a)) b))))))))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double tmp;
if ((y <= -6.2e+59) || !(y <= 1.8e-7)) {
tmp = x + ((z / y) - (a * (x / y)));
} else {
tmp = (t + (y * (230661.510616 + (y * (27464.7644705 + (y * z)))))) / (i + (y * (c + (y * ((y * (y + a)) + b)))));
}
return tmp;
}
real(8) function code(x, y, z, t, a, b, c, i)
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), intent (in) :: i
real(8) :: tmp
if ((y <= (-6.2d+59)) .or. (.not. (y <= 1.8d-7))) then
tmp = x + ((z / y) - (a * (x / y)))
else
tmp = (t + (y * (230661.510616d0 + (y * (27464.7644705d0 + (y * z)))))) / (i + (y * (c + (y * ((y * (y + 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 i) {
double tmp;
if ((y <= -6.2e+59) || !(y <= 1.8e-7)) {
tmp = x + ((z / y) - (a * (x / y)));
} else {
tmp = (t + (y * (230661.510616 + (y * (27464.7644705 + (y * z)))))) / (i + (y * (c + (y * ((y * (y + a)) + b)))));
}
return tmp;
}
def code(x, y, z, t, a, b, c, i): tmp = 0 if (y <= -6.2e+59) or not (y <= 1.8e-7): tmp = x + ((z / y) - (a * (x / y))) else: tmp = (t + (y * (230661.510616 + (y * (27464.7644705 + (y * z)))))) / (i + (y * (c + (y * ((y * (y + a)) + b))))) return tmp
function code(x, y, z, t, a, b, c, i) tmp = 0.0 if ((y <= -6.2e+59) || !(y <= 1.8e-7)) tmp = Float64(x + Float64(Float64(z / y) - Float64(a * Float64(x / y)))); else tmp = Float64(Float64(t + Float64(y * Float64(230661.510616 + Float64(y * Float64(27464.7644705 + Float64(y * z)))))) / Float64(i + Float64(y * Float64(c + Float64(y * Float64(Float64(y * Float64(y + a)) + b)))))); end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c, i) tmp = 0.0; if ((y <= -6.2e+59) || ~((y <= 1.8e-7))) tmp = x + ((z / y) - (a * (x / y))); else tmp = (t + (y * (230661.510616 + (y * (27464.7644705 + (y * z)))))) / (i + (y * (c + (y * ((y * (y + a)) + b))))); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_, i_] := If[Or[LessEqual[y, -6.2e+59], N[Not[LessEqual[y, 1.8e-7]], $MachinePrecision]], N[(x + N[(N[(z / y), $MachinePrecision] - N[(a * N[(x / y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(t + N[(y * N[(230661.510616 + N[(y * N[(27464.7644705 + N[(y * z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(i + N[(y * N[(c + N[(y * N[(N[(y * N[(y + a), $MachinePrecision]), $MachinePrecision] + b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -6.2 \cdot 10^{+59} \lor \neg \left(y \leq 1.8 \cdot 10^{-7}\right):\\
\;\;\;\;x + \left(\frac{z}{y} - a \cdot \frac{x}{y}\right)\\
\mathbf{else}:\\
\;\;\;\;\frac{t + y \cdot \left(230661.510616 + y \cdot \left(27464.7644705 + y \cdot z\right)\right)}{i + y \cdot \left(c + y \cdot \left(y \cdot \left(y + a\right) + b\right)\right)}\\
\end{array}
\end{array}
if y < -6.20000000000000029e59 or 1.79999999999999997e-7 < y Initial program 7.8%
Taylor expanded in y around inf 62.6%
associate--l+62.6%
associate-/l*69.8%
Simplified69.8%
if -6.20000000000000029e59 < y < 1.79999999999999997e-7Initial program 99.0%
Taylor expanded in x around 0 95.1%
Final simplification83.9%
(FPCore (x y z t a b c i)
:precision binary64
(let* ((t_1 (+ x (- (/ z y) (* a (/ x y))))))
(if (<= y -2.35e+27)
t_1
(if (<= y -0.25)
(/
(+ t (* y (+ 230661.510616 (* y 27464.7644705))))
(* y (+ c (* y (+ (* y (+ y a)) b)))))
(if (<= y 1700000.0)
(/
(+ t (* y (+ 230661.510616 (* y (+ 27464.7644705 (* y z))))))
(+ i (* y (+ c (* y b)))))
t_1)))))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double t_1 = x + ((z / y) - (a * (x / y)));
double tmp;
if (y <= -2.35e+27) {
tmp = t_1;
} else if (y <= -0.25) {
tmp = (t + (y * (230661.510616 + (y * 27464.7644705)))) / (y * (c + (y * ((y * (y + a)) + b))));
} else if (y <= 1700000.0) {
tmp = (t + (y * (230661.510616 + (y * (27464.7644705 + (y * z)))))) / (i + (y * (c + (y * b))));
} else {
tmp = t_1;
}
return tmp;
}
real(8) function code(x, y, z, t, a, b, c, i)
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), intent (in) :: i
real(8) :: t_1
real(8) :: tmp
t_1 = x + ((z / y) - (a * (x / y)))
if (y <= (-2.35d+27)) then
tmp = t_1
else if (y <= (-0.25d0)) then
tmp = (t + (y * (230661.510616d0 + (y * 27464.7644705d0)))) / (y * (c + (y * ((y * (y + a)) + b))))
else if (y <= 1700000.0d0) then
tmp = (t + (y * (230661.510616d0 + (y * (27464.7644705d0 + (y * z)))))) / (i + (y * (c + (y * b))))
else
tmp = t_1
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double t_1 = x + ((z / y) - (a * (x / y)));
double tmp;
if (y <= -2.35e+27) {
tmp = t_1;
} else if (y <= -0.25) {
tmp = (t + (y * (230661.510616 + (y * 27464.7644705)))) / (y * (c + (y * ((y * (y + a)) + b))));
} else if (y <= 1700000.0) {
tmp = (t + (y * (230661.510616 + (y * (27464.7644705 + (y * z)))))) / (i + (y * (c + (y * b))));
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b, c, i): t_1 = x + ((z / y) - (a * (x / y))) tmp = 0 if y <= -2.35e+27: tmp = t_1 elif y <= -0.25: tmp = (t + (y * (230661.510616 + (y * 27464.7644705)))) / (y * (c + (y * ((y * (y + a)) + b)))) elif y <= 1700000.0: tmp = (t + (y * (230661.510616 + (y * (27464.7644705 + (y * z)))))) / (i + (y * (c + (y * b)))) else: tmp = t_1 return tmp
function code(x, y, z, t, a, b, c, i) t_1 = Float64(x + Float64(Float64(z / y) - Float64(a * Float64(x / y)))) tmp = 0.0 if (y <= -2.35e+27) tmp = t_1; elseif (y <= -0.25) tmp = Float64(Float64(t + Float64(y * Float64(230661.510616 + Float64(y * 27464.7644705)))) / Float64(y * Float64(c + Float64(y * Float64(Float64(y * Float64(y + a)) + b))))); elseif (y <= 1700000.0) tmp = Float64(Float64(t + Float64(y * Float64(230661.510616 + Float64(y * Float64(27464.7644705 + Float64(y * z)))))) / Float64(i + Float64(y * Float64(c + Float64(y * b))))); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c, i) t_1 = x + ((z / y) - (a * (x / y))); tmp = 0.0; if (y <= -2.35e+27) tmp = t_1; elseif (y <= -0.25) tmp = (t + (y * (230661.510616 + (y * 27464.7644705)))) / (y * (c + (y * ((y * (y + a)) + b)))); elseif (y <= 1700000.0) tmp = (t + (y * (230661.510616 + (y * (27464.7644705 + (y * z)))))) / (i + (y * (c + (y * b)))); else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_, i_] := Block[{t$95$1 = N[(x + N[(N[(z / y), $MachinePrecision] - N[(a * N[(x / y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -2.35e+27], t$95$1, If[LessEqual[y, -0.25], N[(N[(t + N[(y * N[(230661.510616 + N[(y * 27464.7644705), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(y * N[(c + N[(y * N[(N[(y * N[(y + a), $MachinePrecision]), $MachinePrecision] + b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 1700000.0], N[(N[(t + N[(y * N[(230661.510616 + N[(y * N[(27464.7644705 + N[(y * z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(i + N[(y * N[(c + N[(y * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$1]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x + \left(\frac{z}{y} - a \cdot \frac{x}{y}\right)\\
\mathbf{if}\;y \leq -2.35 \cdot 10^{+27}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y \leq -0.25:\\
\;\;\;\;\frac{t + y \cdot \left(230661.510616 + y \cdot 27464.7644705\right)}{y \cdot \left(c + y \cdot \left(y \cdot \left(y + a\right) + b\right)\right)}\\
\mathbf{elif}\;y \leq 1700000:\\
\;\;\;\;\frac{t + y \cdot \left(230661.510616 + y \cdot \left(27464.7644705 + y \cdot z\right)\right)}{i + y \cdot \left(c + y \cdot b\right)}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if y < -2.34999999999999988e27 or 1.7e6 < y Initial program 7.0%
Taylor expanded in y around inf 64.9%
associate--l+64.9%
associate-/l*72.2%
Simplified72.2%
if -2.34999999999999988e27 < y < -0.25Initial program 83.8%
Taylor expanded in i around 0 84.0%
Taylor expanded in y around 0 75.7%
if -0.25 < y < 1.7e6Initial program 99.7%
Taylor expanded in x around 0 95.6%
Taylor expanded in y around 0 91.6%
Final simplification82.7%
(FPCore (x y z t a b c i)
:precision binary64
(let* ((t_1 (+ x (- (/ z y) (* a (/ x y))))))
(if (<= y -1.9e+26)
t_1
(if (<= y -4.5e-21)
(/ t (* y (+ c (* y (+ (* y (+ y a)) b)))))
(if (<= y 2.2e-173)
(/ (* t (+ 1.0 (* 230661.510616 (/ y t)))) (+ i (* y c)))
(if (<= y 3.05e-29)
(/ t (+ i (* y (+ c (* y (+ b (* y a)))))))
t_1))))))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double t_1 = x + ((z / y) - (a * (x / y)));
double tmp;
if (y <= -1.9e+26) {
tmp = t_1;
} else if (y <= -4.5e-21) {
tmp = t / (y * (c + (y * ((y * (y + a)) + b))));
} else if (y <= 2.2e-173) {
tmp = (t * (1.0 + (230661.510616 * (y / t)))) / (i + (y * c));
} else if (y <= 3.05e-29) {
tmp = t / (i + (y * (c + (y * (b + (y * a))))));
} else {
tmp = t_1;
}
return tmp;
}
real(8) function code(x, y, z, t, a, b, c, i)
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), intent (in) :: i
real(8) :: t_1
real(8) :: tmp
t_1 = x + ((z / y) - (a * (x / y)))
if (y <= (-1.9d+26)) then
tmp = t_1
else if (y <= (-4.5d-21)) then
tmp = t / (y * (c + (y * ((y * (y + a)) + b))))
else if (y <= 2.2d-173) then
tmp = (t * (1.0d0 + (230661.510616d0 * (y / t)))) / (i + (y * c))
else if (y <= 3.05d-29) then
tmp = t / (i + (y * (c + (y * (b + (y * a))))))
else
tmp = t_1
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double t_1 = x + ((z / y) - (a * (x / y)));
double tmp;
if (y <= -1.9e+26) {
tmp = t_1;
} else if (y <= -4.5e-21) {
tmp = t / (y * (c + (y * ((y * (y + a)) + b))));
} else if (y <= 2.2e-173) {
tmp = (t * (1.0 + (230661.510616 * (y / t)))) / (i + (y * c));
} else if (y <= 3.05e-29) {
tmp = t / (i + (y * (c + (y * (b + (y * a))))));
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b, c, i): t_1 = x + ((z / y) - (a * (x / y))) tmp = 0 if y <= -1.9e+26: tmp = t_1 elif y <= -4.5e-21: tmp = t / (y * (c + (y * ((y * (y + a)) + b)))) elif y <= 2.2e-173: tmp = (t * (1.0 + (230661.510616 * (y / t)))) / (i + (y * c)) elif y <= 3.05e-29: tmp = t / (i + (y * (c + (y * (b + (y * a)))))) else: tmp = t_1 return tmp
function code(x, y, z, t, a, b, c, i) t_1 = Float64(x + Float64(Float64(z / y) - Float64(a * Float64(x / y)))) tmp = 0.0 if (y <= -1.9e+26) tmp = t_1; elseif (y <= -4.5e-21) tmp = Float64(t / Float64(y * Float64(c + Float64(y * Float64(Float64(y * Float64(y + a)) + b))))); elseif (y <= 2.2e-173) tmp = Float64(Float64(t * Float64(1.0 + Float64(230661.510616 * Float64(y / t)))) / Float64(i + Float64(y * c))); elseif (y <= 3.05e-29) tmp = Float64(t / Float64(i + Float64(y * Float64(c + Float64(y * Float64(b + Float64(y * a))))))); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c, i) t_1 = x + ((z / y) - (a * (x / y))); tmp = 0.0; if (y <= -1.9e+26) tmp = t_1; elseif (y <= -4.5e-21) tmp = t / (y * (c + (y * ((y * (y + a)) + b)))); elseif (y <= 2.2e-173) tmp = (t * (1.0 + (230661.510616 * (y / t)))) / (i + (y * c)); elseif (y <= 3.05e-29) tmp = t / (i + (y * (c + (y * (b + (y * a)))))); else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_, i_] := Block[{t$95$1 = N[(x + N[(N[(z / y), $MachinePrecision] - N[(a * N[(x / y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -1.9e+26], t$95$1, If[LessEqual[y, -4.5e-21], N[(t / N[(y * N[(c + N[(y * N[(N[(y * N[(y + a), $MachinePrecision]), $MachinePrecision] + b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 2.2e-173], N[(N[(t * N[(1.0 + N[(230661.510616 * N[(y / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(i + N[(y * c), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 3.05e-29], N[(t / N[(i + N[(y * N[(c + N[(y * N[(b + N[(y * a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$1]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x + \left(\frac{z}{y} - a \cdot \frac{x}{y}\right)\\
\mathbf{if}\;y \leq -1.9 \cdot 10^{+26}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y \leq -4.5 \cdot 10^{-21}:\\
\;\;\;\;\frac{t}{y \cdot \left(c + y \cdot \left(y \cdot \left(y + a\right) + b\right)\right)}\\
\mathbf{elif}\;y \leq 2.2 \cdot 10^{-173}:\\
\;\;\;\;\frac{t \cdot \left(1 + 230661.510616 \cdot \frac{y}{t}\right)}{i + y \cdot c}\\
\mathbf{elif}\;y \leq 3.05 \cdot 10^{-29}:\\
\;\;\;\;\frac{t}{i + y \cdot \left(c + y \cdot \left(b + y \cdot a\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if y < -1.9000000000000001e26 or 3.05e-29 < y Initial program 16.5%
Taylor expanded in y around inf 58.5%
associate--l+58.5%
associate-/l*65.0%
Simplified65.0%
if -1.9000000000000001e26 < y < -4.49999999999999968e-21Initial program 90.1%
Taylor expanded in t around inf 47.9%
Taylor expanded in i around 0 47.9%
if -4.49999999999999968e-21 < y < 2.1999999999999999e-173Initial program 99.8%
Taylor expanded in y around 0 91.8%
*-commutative91.8%
Simplified91.8%
Taylor expanded in y around 0 87.6%
Taylor expanded in t around inf 87.6%
if 2.1999999999999999e-173 < y < 3.05e-29Initial program 99.6%
Taylor expanded in t around inf 72.0%
Taylor expanded in y around 0 72.0%
*-commutative72.0%
Simplified72.0%
Final simplification72.5%
(FPCore (x y z t a b c i)
:precision binary64
(let* ((t_1 (+ x (- (/ z y) (* a (/ x y))))))
(if (<= y -1.7e+27)
t_1
(if (<= y -2.05e-22)
(/ t (* y (+ c (* y (+ (* y (+ y a)) b)))))
(if (<= y 2e-173)
(/ (+ t (* y 230661.510616)) (+ i (* y c)))
(if (<= y 3.05e-29)
(/ t (+ i (* y (+ c (* y (+ b (* y a)))))))
t_1))))))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double t_1 = x + ((z / y) - (a * (x / y)));
double tmp;
if (y <= -1.7e+27) {
tmp = t_1;
} else if (y <= -2.05e-22) {
tmp = t / (y * (c + (y * ((y * (y + a)) + b))));
} else if (y <= 2e-173) {
tmp = (t + (y * 230661.510616)) / (i + (y * c));
} else if (y <= 3.05e-29) {
tmp = t / (i + (y * (c + (y * (b + (y * a))))));
} else {
tmp = t_1;
}
return tmp;
}
real(8) function code(x, y, z, t, a, b, c, i)
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), intent (in) :: i
real(8) :: t_1
real(8) :: tmp
t_1 = x + ((z / y) - (a * (x / y)))
if (y <= (-1.7d+27)) then
tmp = t_1
else if (y <= (-2.05d-22)) then
tmp = t / (y * (c + (y * ((y * (y + a)) + b))))
else if (y <= 2d-173) then
tmp = (t + (y * 230661.510616d0)) / (i + (y * c))
else if (y <= 3.05d-29) then
tmp = t / (i + (y * (c + (y * (b + (y * a))))))
else
tmp = t_1
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double t_1 = x + ((z / y) - (a * (x / y)));
double tmp;
if (y <= -1.7e+27) {
tmp = t_1;
} else if (y <= -2.05e-22) {
tmp = t / (y * (c + (y * ((y * (y + a)) + b))));
} else if (y <= 2e-173) {
tmp = (t + (y * 230661.510616)) / (i + (y * c));
} else if (y <= 3.05e-29) {
tmp = t / (i + (y * (c + (y * (b + (y * a))))));
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b, c, i): t_1 = x + ((z / y) - (a * (x / y))) tmp = 0 if y <= -1.7e+27: tmp = t_1 elif y <= -2.05e-22: tmp = t / (y * (c + (y * ((y * (y + a)) + b)))) elif y <= 2e-173: tmp = (t + (y * 230661.510616)) / (i + (y * c)) elif y <= 3.05e-29: tmp = t / (i + (y * (c + (y * (b + (y * a)))))) else: tmp = t_1 return tmp
function code(x, y, z, t, a, b, c, i) t_1 = Float64(x + Float64(Float64(z / y) - Float64(a * Float64(x / y)))) tmp = 0.0 if (y <= -1.7e+27) tmp = t_1; elseif (y <= -2.05e-22) tmp = Float64(t / Float64(y * Float64(c + Float64(y * Float64(Float64(y * Float64(y + a)) + b))))); elseif (y <= 2e-173) tmp = Float64(Float64(t + Float64(y * 230661.510616)) / Float64(i + Float64(y * c))); elseif (y <= 3.05e-29) tmp = Float64(t / Float64(i + Float64(y * Float64(c + Float64(y * Float64(b + Float64(y * a))))))); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c, i) t_1 = x + ((z / y) - (a * (x / y))); tmp = 0.0; if (y <= -1.7e+27) tmp = t_1; elseif (y <= -2.05e-22) tmp = t / (y * (c + (y * ((y * (y + a)) + b)))); elseif (y <= 2e-173) tmp = (t + (y * 230661.510616)) / (i + (y * c)); elseif (y <= 3.05e-29) tmp = t / (i + (y * (c + (y * (b + (y * a)))))); else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_, i_] := Block[{t$95$1 = N[(x + N[(N[(z / y), $MachinePrecision] - N[(a * N[(x / y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -1.7e+27], t$95$1, If[LessEqual[y, -2.05e-22], N[(t / N[(y * N[(c + N[(y * N[(N[(y * N[(y + a), $MachinePrecision]), $MachinePrecision] + b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 2e-173], N[(N[(t + N[(y * 230661.510616), $MachinePrecision]), $MachinePrecision] / N[(i + N[(y * c), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 3.05e-29], N[(t / N[(i + N[(y * N[(c + N[(y * N[(b + N[(y * a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$1]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x + \left(\frac{z}{y} - a \cdot \frac{x}{y}\right)\\
\mathbf{if}\;y \leq -1.7 \cdot 10^{+27}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y \leq -2.05 \cdot 10^{-22}:\\
\;\;\;\;\frac{t}{y \cdot \left(c + y \cdot \left(y \cdot \left(y + a\right) + b\right)\right)}\\
\mathbf{elif}\;y \leq 2 \cdot 10^{-173}:\\
\;\;\;\;\frac{t + y \cdot 230661.510616}{i + y \cdot c}\\
\mathbf{elif}\;y \leq 3.05 \cdot 10^{-29}:\\
\;\;\;\;\frac{t}{i + y \cdot \left(c + y \cdot \left(b + y \cdot a\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if y < -1.7e27 or 3.05e-29 < y Initial program 16.5%
Taylor expanded in y around inf 58.5%
associate--l+58.5%
associate-/l*65.0%
Simplified65.0%
if -1.7e27 < y < -2.05e-22Initial program 90.1%
Taylor expanded in t around inf 47.9%
Taylor expanded in i around 0 47.9%
if -2.05e-22 < y < 2.0000000000000001e-173Initial program 99.8%
Taylor expanded in y around 0 91.8%
*-commutative91.8%
Simplified91.8%
Taylor expanded in y around 0 87.6%
if 2.0000000000000001e-173 < y < 3.05e-29Initial program 99.6%
Taylor expanded in t around inf 72.0%
Taylor expanded in y around 0 72.0%
*-commutative72.0%
Simplified72.0%
Final simplification72.5%
(FPCore (x y z t a b c i)
:precision binary64
(let* ((t_1 (+ x (- (/ z y) (* a (/ x y))))))
(if (<= y -6.5e+26)
t_1
(if (<= y 2.6e-39)
(/
(+ t (* y 230661.510616))
(+ i (* y (+ c (* y (+ (* y (+ y a)) b))))))
(if (<= y 5e-10)
(/
(+ t (* y (+ 230661.510616 (* y (+ 27464.7644705 (* y z))))))
(+ i (* y c)))
t_1)))))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double t_1 = x + ((z / y) - (a * (x / y)));
double tmp;
if (y <= -6.5e+26) {
tmp = t_1;
} else if (y <= 2.6e-39) {
tmp = (t + (y * 230661.510616)) / (i + (y * (c + (y * ((y * (y + a)) + b)))));
} else if (y <= 5e-10) {
tmp = (t + (y * (230661.510616 + (y * (27464.7644705 + (y * z)))))) / (i + (y * c));
} else {
tmp = t_1;
}
return tmp;
}
real(8) function code(x, y, z, t, a, b, c, i)
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), intent (in) :: i
real(8) :: t_1
real(8) :: tmp
t_1 = x + ((z / y) - (a * (x / y)))
if (y <= (-6.5d+26)) then
tmp = t_1
else if (y <= 2.6d-39) then
tmp = (t + (y * 230661.510616d0)) / (i + (y * (c + (y * ((y * (y + a)) + b)))))
else if (y <= 5d-10) then
tmp = (t + (y * (230661.510616d0 + (y * (27464.7644705d0 + (y * z)))))) / (i + (y * c))
else
tmp = t_1
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double t_1 = x + ((z / y) - (a * (x / y)));
double tmp;
if (y <= -6.5e+26) {
tmp = t_1;
} else if (y <= 2.6e-39) {
tmp = (t + (y * 230661.510616)) / (i + (y * (c + (y * ((y * (y + a)) + b)))));
} else if (y <= 5e-10) {
tmp = (t + (y * (230661.510616 + (y * (27464.7644705 + (y * z)))))) / (i + (y * c));
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b, c, i): t_1 = x + ((z / y) - (a * (x / y))) tmp = 0 if y <= -6.5e+26: tmp = t_1 elif y <= 2.6e-39: tmp = (t + (y * 230661.510616)) / (i + (y * (c + (y * ((y * (y + a)) + b))))) elif y <= 5e-10: tmp = (t + (y * (230661.510616 + (y * (27464.7644705 + (y * z)))))) / (i + (y * c)) else: tmp = t_1 return tmp
function code(x, y, z, t, a, b, c, i) t_1 = Float64(x + Float64(Float64(z / y) - Float64(a * Float64(x / y)))) tmp = 0.0 if (y <= -6.5e+26) tmp = t_1; elseif (y <= 2.6e-39) tmp = Float64(Float64(t + Float64(y * 230661.510616)) / Float64(i + Float64(y * Float64(c + Float64(y * Float64(Float64(y * Float64(y + a)) + b)))))); elseif (y <= 5e-10) tmp = Float64(Float64(t + Float64(y * Float64(230661.510616 + Float64(y * Float64(27464.7644705 + Float64(y * z)))))) / Float64(i + Float64(y * c))); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c, i) t_1 = x + ((z / y) - (a * (x / y))); tmp = 0.0; if (y <= -6.5e+26) tmp = t_1; elseif (y <= 2.6e-39) tmp = (t + (y * 230661.510616)) / (i + (y * (c + (y * ((y * (y + a)) + b))))); elseif (y <= 5e-10) tmp = (t + (y * (230661.510616 + (y * (27464.7644705 + (y * z)))))) / (i + (y * c)); else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_, i_] := Block[{t$95$1 = N[(x + N[(N[(z / y), $MachinePrecision] - N[(a * N[(x / y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -6.5e+26], t$95$1, If[LessEqual[y, 2.6e-39], N[(N[(t + N[(y * 230661.510616), $MachinePrecision]), $MachinePrecision] / N[(i + N[(y * N[(c + N[(y * N[(N[(y * N[(y + a), $MachinePrecision]), $MachinePrecision] + b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 5e-10], N[(N[(t + N[(y * N[(230661.510616 + N[(y * N[(27464.7644705 + N[(y * z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(i + N[(y * c), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$1]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x + \left(\frac{z}{y} - a \cdot \frac{x}{y}\right)\\
\mathbf{if}\;y \leq -6.5 \cdot 10^{+26}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y \leq 2.6 \cdot 10^{-39}:\\
\;\;\;\;\frac{t + y \cdot 230661.510616}{i + y \cdot \left(c + y \cdot \left(y \cdot \left(y + a\right) + b\right)\right)}\\
\mathbf{elif}\;y \leq 5 \cdot 10^{-10}:\\
\;\;\;\;\frac{t + y \cdot \left(230661.510616 + y \cdot \left(27464.7644705 + y \cdot z\right)\right)}{i + y \cdot c}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if y < -6.50000000000000022e26 or 5.00000000000000031e-10 < y Initial program 10.1%
Taylor expanded in y around inf 62.7%
associate--l+62.7%
associate-/l*69.8%
Simplified69.8%
if -6.50000000000000022e26 < y < 2.6e-39Initial program 99.0%
Taylor expanded in y around 0 89.8%
*-commutative89.8%
Simplified89.8%
if 2.6e-39 < y < 5.00000000000000031e-10Initial program 99.6%
Taylor expanded in x around 0 84.9%
Taylor expanded in y around 0 80.6%
Taylor expanded in b around 0 66.0%
*-commutative66.0%
Simplified66.0%
Final simplification79.4%
(FPCore (x y z t a b c i)
:precision binary64
(let* ((t_1 (+ x (- (/ z y) (* a (/ x y))))))
(if (<= y -1.8e+27)
t_1
(if (<= y 5.9e-43)
(/ (+ t (* y 230661.510616)) (+ i (* y (+ c (* y (+ b (* y a)))))))
(if (<= y 5e-10)
(/
(+ t (* y (+ 230661.510616 (* y (+ 27464.7644705 (* y z))))))
(+ i (* y c)))
t_1)))))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double t_1 = x + ((z / y) - (a * (x / y)));
double tmp;
if (y <= -1.8e+27) {
tmp = t_1;
} else if (y <= 5.9e-43) {
tmp = (t + (y * 230661.510616)) / (i + (y * (c + (y * (b + (y * a))))));
} else if (y <= 5e-10) {
tmp = (t + (y * (230661.510616 + (y * (27464.7644705 + (y * z)))))) / (i + (y * c));
} else {
tmp = t_1;
}
return tmp;
}
real(8) function code(x, y, z, t, a, b, c, i)
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), intent (in) :: i
real(8) :: t_1
real(8) :: tmp
t_1 = x + ((z / y) - (a * (x / y)))
if (y <= (-1.8d+27)) then
tmp = t_1
else if (y <= 5.9d-43) then
tmp = (t + (y * 230661.510616d0)) / (i + (y * (c + (y * (b + (y * a))))))
else if (y <= 5d-10) then
tmp = (t + (y * (230661.510616d0 + (y * (27464.7644705d0 + (y * z)))))) / (i + (y * c))
else
tmp = t_1
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double t_1 = x + ((z / y) - (a * (x / y)));
double tmp;
if (y <= -1.8e+27) {
tmp = t_1;
} else if (y <= 5.9e-43) {
tmp = (t + (y * 230661.510616)) / (i + (y * (c + (y * (b + (y * a))))));
} else if (y <= 5e-10) {
tmp = (t + (y * (230661.510616 + (y * (27464.7644705 + (y * z)))))) / (i + (y * c));
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b, c, i): t_1 = x + ((z / y) - (a * (x / y))) tmp = 0 if y <= -1.8e+27: tmp = t_1 elif y <= 5.9e-43: tmp = (t + (y * 230661.510616)) / (i + (y * (c + (y * (b + (y * a)))))) elif y <= 5e-10: tmp = (t + (y * (230661.510616 + (y * (27464.7644705 + (y * z)))))) / (i + (y * c)) else: tmp = t_1 return tmp
function code(x, y, z, t, a, b, c, i) t_1 = Float64(x + Float64(Float64(z / y) - Float64(a * Float64(x / y)))) tmp = 0.0 if (y <= -1.8e+27) tmp = t_1; elseif (y <= 5.9e-43) tmp = Float64(Float64(t + Float64(y * 230661.510616)) / Float64(i + Float64(y * Float64(c + Float64(y * Float64(b + Float64(y * a))))))); elseif (y <= 5e-10) tmp = Float64(Float64(t + Float64(y * Float64(230661.510616 + Float64(y * Float64(27464.7644705 + Float64(y * z)))))) / Float64(i + Float64(y * c))); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c, i) t_1 = x + ((z / y) - (a * (x / y))); tmp = 0.0; if (y <= -1.8e+27) tmp = t_1; elseif (y <= 5.9e-43) tmp = (t + (y * 230661.510616)) / (i + (y * (c + (y * (b + (y * a)))))); elseif (y <= 5e-10) tmp = (t + (y * (230661.510616 + (y * (27464.7644705 + (y * z)))))) / (i + (y * c)); else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_, i_] := Block[{t$95$1 = N[(x + N[(N[(z / y), $MachinePrecision] - N[(a * N[(x / y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -1.8e+27], t$95$1, If[LessEqual[y, 5.9e-43], N[(N[(t + N[(y * 230661.510616), $MachinePrecision]), $MachinePrecision] / N[(i + N[(y * N[(c + N[(y * N[(b + N[(y * a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 5e-10], N[(N[(t + N[(y * N[(230661.510616 + N[(y * N[(27464.7644705 + N[(y * z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(i + N[(y * c), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$1]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x + \left(\frac{z}{y} - a \cdot \frac{x}{y}\right)\\
\mathbf{if}\;y \leq -1.8 \cdot 10^{+27}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y \leq 5.9 \cdot 10^{-43}:\\
\;\;\;\;\frac{t + y \cdot 230661.510616}{i + y \cdot \left(c + y \cdot \left(b + y \cdot a\right)\right)}\\
\mathbf{elif}\;y \leq 5 \cdot 10^{-10}:\\
\;\;\;\;\frac{t + y \cdot \left(230661.510616 + y \cdot \left(27464.7644705 + y \cdot z\right)\right)}{i + y \cdot c}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if y < -1.79999999999999991e27 or 5.00000000000000031e-10 < y Initial program 10.1%
Taylor expanded in y around inf 62.7%
associate--l+62.7%
associate-/l*69.8%
Simplified69.8%
if -1.79999999999999991e27 < y < 5.89999999999999976e-43Initial program 99.0%
Taylor expanded in y around 0 89.8%
*-commutative89.8%
Simplified89.8%
Taylor expanded in y around 0 89.0%
if 5.89999999999999976e-43 < y < 5.00000000000000031e-10Initial program 99.6%
Taylor expanded in x around 0 84.9%
Taylor expanded in y around 0 80.6%
Taylor expanded in b around 0 66.0%
*-commutative66.0%
Simplified66.0%
Final simplification79.1%
(FPCore (x y z t a b c i)
:precision binary64
(let* ((t_1 (+ x (- (/ z y) (* a (/ x y))))))
(if (<= y -1.45e+26)
t_1
(if (<= y 4.6e-23)
(/ (+ t (* y 230661.510616)) (+ i (* y (+ c (* y (+ b (* y a)))))))
(if (<= y 7.2e+23)
(/
(+
230661.510616
(+ (* y (+ 27464.7644705 (* y (+ (* x y) z)))) (/ t y)))
c)
t_1)))))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double t_1 = x + ((z / y) - (a * (x / y)));
double tmp;
if (y <= -1.45e+26) {
tmp = t_1;
} else if (y <= 4.6e-23) {
tmp = (t + (y * 230661.510616)) / (i + (y * (c + (y * (b + (y * a))))));
} else if (y <= 7.2e+23) {
tmp = (230661.510616 + ((y * (27464.7644705 + (y * ((x * y) + z)))) + (t / y))) / c;
} else {
tmp = t_1;
}
return tmp;
}
real(8) function code(x, y, z, t, a, b, c, i)
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), intent (in) :: i
real(8) :: t_1
real(8) :: tmp
t_1 = x + ((z / y) - (a * (x / y)))
if (y <= (-1.45d+26)) then
tmp = t_1
else if (y <= 4.6d-23) then
tmp = (t + (y * 230661.510616d0)) / (i + (y * (c + (y * (b + (y * a))))))
else if (y <= 7.2d+23) then
tmp = (230661.510616d0 + ((y * (27464.7644705d0 + (y * ((x * y) + z)))) + (t / y))) / c
else
tmp = t_1
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double t_1 = x + ((z / y) - (a * (x / y)));
double tmp;
if (y <= -1.45e+26) {
tmp = t_1;
} else if (y <= 4.6e-23) {
tmp = (t + (y * 230661.510616)) / (i + (y * (c + (y * (b + (y * a))))));
} else if (y <= 7.2e+23) {
tmp = (230661.510616 + ((y * (27464.7644705 + (y * ((x * y) + z)))) + (t / y))) / c;
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b, c, i): t_1 = x + ((z / y) - (a * (x / y))) tmp = 0 if y <= -1.45e+26: tmp = t_1 elif y <= 4.6e-23: tmp = (t + (y * 230661.510616)) / (i + (y * (c + (y * (b + (y * a)))))) elif y <= 7.2e+23: tmp = (230661.510616 + ((y * (27464.7644705 + (y * ((x * y) + z)))) + (t / y))) / c else: tmp = t_1 return tmp
function code(x, y, z, t, a, b, c, i) t_1 = Float64(x + Float64(Float64(z / y) - Float64(a * Float64(x / y)))) tmp = 0.0 if (y <= -1.45e+26) tmp = t_1; elseif (y <= 4.6e-23) tmp = Float64(Float64(t + Float64(y * 230661.510616)) / Float64(i + Float64(y * Float64(c + Float64(y * Float64(b + Float64(y * a))))))); elseif (y <= 7.2e+23) tmp = Float64(Float64(230661.510616 + Float64(Float64(y * Float64(27464.7644705 + Float64(y * Float64(Float64(x * y) + z)))) + Float64(t / y))) / c); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c, i) t_1 = x + ((z / y) - (a * (x / y))); tmp = 0.0; if (y <= -1.45e+26) tmp = t_1; elseif (y <= 4.6e-23) tmp = (t + (y * 230661.510616)) / (i + (y * (c + (y * (b + (y * a)))))); elseif (y <= 7.2e+23) tmp = (230661.510616 + ((y * (27464.7644705 + (y * ((x * y) + z)))) + (t / y))) / c; else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_, i_] := Block[{t$95$1 = N[(x + N[(N[(z / y), $MachinePrecision] - N[(a * N[(x / y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -1.45e+26], t$95$1, If[LessEqual[y, 4.6e-23], N[(N[(t + N[(y * 230661.510616), $MachinePrecision]), $MachinePrecision] / N[(i + N[(y * N[(c + N[(y * N[(b + N[(y * a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 7.2e+23], N[(N[(230661.510616 + N[(N[(y * N[(27464.7644705 + N[(y * N[(N[(x * y), $MachinePrecision] + z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(t / y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / c), $MachinePrecision], t$95$1]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x + \left(\frac{z}{y} - a \cdot \frac{x}{y}\right)\\
\mathbf{if}\;y \leq -1.45 \cdot 10^{+26}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y \leq 4.6 \cdot 10^{-23}:\\
\;\;\;\;\frac{t + y \cdot 230661.510616}{i + y \cdot \left(c + y \cdot \left(b + y \cdot a\right)\right)}\\
\mathbf{elif}\;y \leq 7.2 \cdot 10^{+23}:\\
\;\;\;\;\frac{230661.510616 + \left(y \cdot \left(27464.7644705 + y \cdot \left(x \cdot y + z\right)\right) + \frac{t}{y}\right)}{c}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if y < -1.45e26 or 7.1999999999999997e23 < y Initial program 6.1%
Taylor expanded in y around inf 65.4%
associate--l+65.4%
associate-/l*72.8%
Simplified72.8%
if -1.45e26 < y < 4.6000000000000002e-23Initial program 99.0%
Taylor expanded in y around 0 88.0%
*-commutative88.0%
Simplified88.0%
Taylor expanded in y around 0 87.3%
if 4.6000000000000002e-23 < y < 7.1999999999999997e23Initial program 99.5%
Taylor expanded in t around 0 99.6%
Taylor expanded in i around 0 76.0%
Taylor expanded in c around inf 51.3%
Final simplification79.3%
(FPCore (x y z t a b c i)
:precision binary64
(let* ((t_1 (+ x (- (/ z y) (* a (/ x y))))))
(if (<= y -5e+26)
t_1
(if (<= y 3.05e-29)
(/ (+ t (* y 230661.510616)) (+ i (* y (+ c (* y b)))))
(if (<= y 9e+19)
(/
(+
230661.510616
(+ (* y (+ 27464.7644705 (* y (+ (* x y) z)))) (/ t y)))
c)
t_1)))))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double t_1 = x + ((z / y) - (a * (x / y)));
double tmp;
if (y <= -5e+26) {
tmp = t_1;
} else if (y <= 3.05e-29) {
tmp = (t + (y * 230661.510616)) / (i + (y * (c + (y * b))));
} else if (y <= 9e+19) {
tmp = (230661.510616 + ((y * (27464.7644705 + (y * ((x * y) + z)))) + (t / y))) / c;
} else {
tmp = t_1;
}
return tmp;
}
real(8) function code(x, y, z, t, a, b, c, i)
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), intent (in) :: i
real(8) :: t_1
real(8) :: tmp
t_1 = x + ((z / y) - (a * (x / y)))
if (y <= (-5d+26)) then
tmp = t_1
else if (y <= 3.05d-29) then
tmp = (t + (y * 230661.510616d0)) / (i + (y * (c + (y * b))))
else if (y <= 9d+19) then
tmp = (230661.510616d0 + ((y * (27464.7644705d0 + (y * ((x * y) + z)))) + (t / y))) / c
else
tmp = t_1
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double t_1 = x + ((z / y) - (a * (x / y)));
double tmp;
if (y <= -5e+26) {
tmp = t_1;
} else if (y <= 3.05e-29) {
tmp = (t + (y * 230661.510616)) / (i + (y * (c + (y * b))));
} else if (y <= 9e+19) {
tmp = (230661.510616 + ((y * (27464.7644705 + (y * ((x * y) + z)))) + (t / y))) / c;
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b, c, i): t_1 = x + ((z / y) - (a * (x / y))) tmp = 0 if y <= -5e+26: tmp = t_1 elif y <= 3.05e-29: tmp = (t + (y * 230661.510616)) / (i + (y * (c + (y * b)))) elif y <= 9e+19: tmp = (230661.510616 + ((y * (27464.7644705 + (y * ((x * y) + z)))) + (t / y))) / c else: tmp = t_1 return tmp
function code(x, y, z, t, a, b, c, i) t_1 = Float64(x + Float64(Float64(z / y) - Float64(a * Float64(x / y)))) tmp = 0.0 if (y <= -5e+26) tmp = t_1; elseif (y <= 3.05e-29) tmp = Float64(Float64(t + Float64(y * 230661.510616)) / Float64(i + Float64(y * Float64(c + Float64(y * b))))); elseif (y <= 9e+19) tmp = Float64(Float64(230661.510616 + Float64(Float64(y * Float64(27464.7644705 + Float64(y * Float64(Float64(x * y) + z)))) + Float64(t / y))) / c); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c, i) t_1 = x + ((z / y) - (a * (x / y))); tmp = 0.0; if (y <= -5e+26) tmp = t_1; elseif (y <= 3.05e-29) tmp = (t + (y * 230661.510616)) / (i + (y * (c + (y * b)))); elseif (y <= 9e+19) tmp = (230661.510616 + ((y * (27464.7644705 + (y * ((x * y) + z)))) + (t / y))) / c; else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_, i_] := Block[{t$95$1 = N[(x + N[(N[(z / y), $MachinePrecision] - N[(a * N[(x / y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -5e+26], t$95$1, If[LessEqual[y, 3.05e-29], N[(N[(t + N[(y * 230661.510616), $MachinePrecision]), $MachinePrecision] / N[(i + N[(y * N[(c + N[(y * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 9e+19], N[(N[(230661.510616 + N[(N[(y * N[(27464.7644705 + N[(y * N[(N[(x * y), $MachinePrecision] + z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(t / y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / c), $MachinePrecision], t$95$1]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x + \left(\frac{z}{y} - a \cdot \frac{x}{y}\right)\\
\mathbf{if}\;y \leq -5 \cdot 10^{+26}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y \leq 3.05 \cdot 10^{-29}:\\
\;\;\;\;\frac{t + y \cdot 230661.510616}{i + y \cdot \left(c + y \cdot b\right)}\\
\mathbf{elif}\;y \leq 9 \cdot 10^{+19}:\\
\;\;\;\;\frac{230661.510616 + \left(y \cdot \left(27464.7644705 + y \cdot \left(x \cdot y + z\right)\right) + \frac{t}{y}\right)}{c}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if y < -5.0000000000000001e26 or 9e19 < y Initial program 6.1%
Taylor expanded in y around inf 65.4%
associate--l+65.4%
associate-/l*72.8%
Simplified72.8%
if -5.0000000000000001e26 < y < 3.05e-29Initial program 99.0%
Taylor expanded in y around 0 88.6%
*-commutative88.6%
Simplified88.6%
Taylor expanded in y around 0 84.2%
if 3.05e-29 < y < 9e19Initial program 99.6%
Taylor expanded in t around 0 99.7%
Taylor expanded in i around 0 72.3%
Taylor expanded in c around inf 44.7%
Final simplification77.0%
(FPCore (x y z t a b c i)
:precision binary64
(let* ((t_1 (+ x (/ (- z (* x a)) y))))
(if (<= y -3.2e+26)
t_1
(if (<= y -1.9e-34)
(/ t (* y (+ c (* y b))))
(if (<= y 2.5e-50) (/ t (+ i (* y c))) (if (<= y 1.85e+253) t_1 x))))))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double t_1 = x + ((z - (x * a)) / y);
double tmp;
if (y <= -3.2e+26) {
tmp = t_1;
} else if (y <= -1.9e-34) {
tmp = t / (y * (c + (y * b)));
} else if (y <= 2.5e-50) {
tmp = t / (i + (y * c));
} else if (y <= 1.85e+253) {
tmp = t_1;
} else {
tmp = x;
}
return tmp;
}
real(8) function code(x, y, z, t, a, b, c, i)
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), intent (in) :: i
real(8) :: t_1
real(8) :: tmp
t_1 = x + ((z - (x * a)) / y)
if (y <= (-3.2d+26)) then
tmp = t_1
else if (y <= (-1.9d-34)) then
tmp = t / (y * (c + (y * b)))
else if (y <= 2.5d-50) then
tmp = t / (i + (y * c))
else if (y <= 1.85d+253) then
tmp = t_1
else
tmp = 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 i) {
double t_1 = x + ((z - (x * a)) / y);
double tmp;
if (y <= -3.2e+26) {
tmp = t_1;
} else if (y <= -1.9e-34) {
tmp = t / (y * (c + (y * b)));
} else if (y <= 2.5e-50) {
tmp = t / (i + (y * c));
} else if (y <= 1.85e+253) {
tmp = t_1;
} else {
tmp = x;
}
return tmp;
}
def code(x, y, z, t, a, b, c, i): t_1 = x + ((z - (x * a)) / y) tmp = 0 if y <= -3.2e+26: tmp = t_1 elif y <= -1.9e-34: tmp = t / (y * (c + (y * b))) elif y <= 2.5e-50: tmp = t / (i + (y * c)) elif y <= 1.85e+253: tmp = t_1 else: tmp = x return tmp
function code(x, y, z, t, a, b, c, i) t_1 = Float64(x + Float64(Float64(z - Float64(x * a)) / y)) tmp = 0.0 if (y <= -3.2e+26) tmp = t_1; elseif (y <= -1.9e-34) tmp = Float64(t / Float64(y * Float64(c + Float64(y * b)))); elseif (y <= 2.5e-50) tmp = Float64(t / Float64(i + Float64(y * c))); elseif (y <= 1.85e+253) tmp = t_1; else tmp = x; end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c, i) t_1 = x + ((z - (x * a)) / y); tmp = 0.0; if (y <= -3.2e+26) tmp = t_1; elseif (y <= -1.9e-34) tmp = t / (y * (c + (y * b))); elseif (y <= 2.5e-50) tmp = t / (i + (y * c)); elseif (y <= 1.85e+253) tmp = t_1; else tmp = x; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_, i_] := Block[{t$95$1 = N[(x + N[(N[(z - N[(x * a), $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -3.2e+26], t$95$1, If[LessEqual[y, -1.9e-34], N[(t / N[(y * N[(c + N[(y * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 2.5e-50], N[(t / N[(i + N[(y * c), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 1.85e+253], t$95$1, x]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x + \frac{z - x \cdot a}{y}\\
\mathbf{if}\;y \leq -3.2 \cdot 10^{+26}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y \leq -1.9 \cdot 10^{-34}:\\
\;\;\;\;\frac{t}{y \cdot \left(c + y \cdot b\right)}\\
\mathbf{elif}\;y \leq 2.5 \cdot 10^{-50}:\\
\;\;\;\;\frac{t}{i + y \cdot c}\\
\mathbf{elif}\;y \leq 1.85 \cdot 10^{+253}:\\
\;\;\;\;t\_1\\
\mathbf{else}:\\
\;\;\;\;x\\
\end{array}
\end{array}
if y < -3.20000000000000029e26 or 2.49999999999999984e-50 < y < 1.85000000000000014e253Initial program 22.9%
Taylor expanded in i around 0 14.6%
Taylor expanded in y around inf 54.4%
associate--l+54.4%
div-sub54.4%
Simplified54.4%
if -3.20000000000000029e26 < y < -1.9000000000000001e-34Initial program 92.3%
Taylor expanded in t around inf 45.0%
Taylor expanded in y around 0 37.4%
*-commutative37.4%
Simplified37.4%
Taylor expanded in i around 0 36.7%
if -1.9000000000000001e-34 < y < 2.49999999999999984e-50Initial program 99.7%
Taylor expanded in t around inf 75.9%
Taylor expanded in y around 0 68.7%
*-commutative68.7%
Simplified68.7%
if 1.85000000000000014e253 < y Initial program 0.0%
Taylor expanded in y around inf 88.8%
Final simplification61.7%
(FPCore (x y z t a b c i)
:precision binary64
(let* ((t_1 (/ t (+ i (* y c)))))
(if (<= y -5.8e+31)
x
(if (<= y -1.15e-151)
t_1
(if (<= y 1e-107)
(/ (+ t (* y 230661.510616)) i)
(if (<= y 4.2e-75) t_1 x))))))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double t_1 = t / (i + (y * c));
double tmp;
if (y <= -5.8e+31) {
tmp = x;
} else if (y <= -1.15e-151) {
tmp = t_1;
} else if (y <= 1e-107) {
tmp = (t + (y * 230661.510616)) / i;
} else if (y <= 4.2e-75) {
tmp = t_1;
} else {
tmp = x;
}
return tmp;
}
real(8) function code(x, y, z, t, a, b, c, i)
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), intent (in) :: i
real(8) :: t_1
real(8) :: tmp
t_1 = t / (i + (y * c))
if (y <= (-5.8d+31)) then
tmp = x
else if (y <= (-1.15d-151)) then
tmp = t_1
else if (y <= 1d-107) then
tmp = (t + (y * 230661.510616d0)) / i
else if (y <= 4.2d-75) then
tmp = t_1
else
tmp = 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 i) {
double t_1 = t / (i + (y * c));
double tmp;
if (y <= -5.8e+31) {
tmp = x;
} else if (y <= -1.15e-151) {
tmp = t_1;
} else if (y <= 1e-107) {
tmp = (t + (y * 230661.510616)) / i;
} else if (y <= 4.2e-75) {
tmp = t_1;
} else {
tmp = x;
}
return tmp;
}
def code(x, y, z, t, a, b, c, i): t_1 = t / (i + (y * c)) tmp = 0 if y <= -5.8e+31: tmp = x elif y <= -1.15e-151: tmp = t_1 elif y <= 1e-107: tmp = (t + (y * 230661.510616)) / i elif y <= 4.2e-75: tmp = t_1 else: tmp = x return tmp
function code(x, y, z, t, a, b, c, i) t_1 = Float64(t / Float64(i + Float64(y * c))) tmp = 0.0 if (y <= -5.8e+31) tmp = x; elseif (y <= -1.15e-151) tmp = t_1; elseif (y <= 1e-107) tmp = Float64(Float64(t + Float64(y * 230661.510616)) / i); elseif (y <= 4.2e-75) tmp = t_1; else tmp = x; end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c, i) t_1 = t / (i + (y * c)); tmp = 0.0; if (y <= -5.8e+31) tmp = x; elseif (y <= -1.15e-151) tmp = t_1; elseif (y <= 1e-107) tmp = (t + (y * 230661.510616)) / i; elseif (y <= 4.2e-75) tmp = t_1; else tmp = x; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_, i_] := Block[{t$95$1 = N[(t / N[(i + N[(y * c), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -5.8e+31], x, If[LessEqual[y, -1.15e-151], t$95$1, If[LessEqual[y, 1e-107], N[(N[(t + N[(y * 230661.510616), $MachinePrecision]), $MachinePrecision] / i), $MachinePrecision], If[LessEqual[y, 4.2e-75], t$95$1, x]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{t}{i + y \cdot c}\\
\mathbf{if}\;y \leq -5.8 \cdot 10^{+31}:\\
\;\;\;\;x\\
\mathbf{elif}\;y \leq -1.15 \cdot 10^{-151}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y \leq 10^{-107}:\\
\;\;\;\;\frac{t + y \cdot 230661.510616}{i}\\
\mathbf{elif}\;y \leq 4.2 \cdot 10^{-75}:\\
\;\;\;\;t\_1\\
\mathbf{else}:\\
\;\;\;\;x\\
\end{array}
\end{array}
if y < -5.8000000000000001e31 or 4.2000000000000002e-75 < y Initial program 23.7%
Taylor expanded in y around inf 42.7%
if -5.8000000000000001e31 < y < -1.14999999999999998e-151 or 1e-107 < y < 4.2000000000000002e-75Initial program 97.4%
Taylor expanded in t around inf 59.3%
Taylor expanded in y around 0 45.4%
*-commutative45.4%
Simplified45.4%
if -1.14999999999999998e-151 < y < 1e-107Initial program 99.9%
Taylor expanded in y around 0 99.9%
*-commutative99.9%
Simplified99.9%
Taylor expanded in i around inf 79.0%
Final simplification53.8%
(FPCore (x y z t a b c i)
:precision binary64
(let* ((t_1 (+ x (- (/ z y) (* a (/ x y))))))
(if (<= y -2.7e+25)
t_1
(if (<= y -2e-23)
(/ t (* y (+ c (* y (+ (* y (+ y a)) b)))))
(if (<= y 2.5e-50) (/ (+ t (* y 230661.510616)) (+ i (* y c))) t_1)))))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double t_1 = x + ((z / y) - (a * (x / y)));
double tmp;
if (y <= -2.7e+25) {
tmp = t_1;
} else if (y <= -2e-23) {
tmp = t / (y * (c + (y * ((y * (y + a)) + b))));
} else if (y <= 2.5e-50) {
tmp = (t + (y * 230661.510616)) / (i + (y * c));
} else {
tmp = t_1;
}
return tmp;
}
real(8) function code(x, y, z, t, a, b, c, i)
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), intent (in) :: i
real(8) :: t_1
real(8) :: tmp
t_1 = x + ((z / y) - (a * (x / y)))
if (y <= (-2.7d+25)) then
tmp = t_1
else if (y <= (-2d-23)) then
tmp = t / (y * (c + (y * ((y * (y + a)) + b))))
else if (y <= 2.5d-50) then
tmp = (t + (y * 230661.510616d0)) / (i + (y * c))
else
tmp = t_1
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double t_1 = x + ((z / y) - (a * (x / y)));
double tmp;
if (y <= -2.7e+25) {
tmp = t_1;
} else if (y <= -2e-23) {
tmp = t / (y * (c + (y * ((y * (y + a)) + b))));
} else if (y <= 2.5e-50) {
tmp = (t + (y * 230661.510616)) / (i + (y * c));
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b, c, i): t_1 = x + ((z / y) - (a * (x / y))) tmp = 0 if y <= -2.7e+25: tmp = t_1 elif y <= -2e-23: tmp = t / (y * (c + (y * ((y * (y + a)) + b)))) elif y <= 2.5e-50: tmp = (t + (y * 230661.510616)) / (i + (y * c)) else: tmp = t_1 return tmp
function code(x, y, z, t, a, b, c, i) t_1 = Float64(x + Float64(Float64(z / y) - Float64(a * Float64(x / y)))) tmp = 0.0 if (y <= -2.7e+25) tmp = t_1; elseif (y <= -2e-23) tmp = Float64(t / Float64(y * Float64(c + Float64(y * Float64(Float64(y * Float64(y + a)) + b))))); elseif (y <= 2.5e-50) tmp = Float64(Float64(t + Float64(y * 230661.510616)) / Float64(i + Float64(y * c))); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c, i) t_1 = x + ((z / y) - (a * (x / y))); tmp = 0.0; if (y <= -2.7e+25) tmp = t_1; elseif (y <= -2e-23) tmp = t / (y * (c + (y * ((y * (y + a)) + b)))); elseif (y <= 2.5e-50) tmp = (t + (y * 230661.510616)) / (i + (y * c)); else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_, i_] := Block[{t$95$1 = N[(x + N[(N[(z / y), $MachinePrecision] - N[(a * N[(x / y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -2.7e+25], t$95$1, If[LessEqual[y, -2e-23], N[(t / N[(y * N[(c + N[(y * N[(N[(y * N[(y + a), $MachinePrecision]), $MachinePrecision] + b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 2.5e-50], N[(N[(t + N[(y * 230661.510616), $MachinePrecision]), $MachinePrecision] / N[(i + N[(y * c), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$1]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x + \left(\frac{z}{y} - a \cdot \frac{x}{y}\right)\\
\mathbf{if}\;y \leq -2.7 \cdot 10^{+25}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y \leq -2 \cdot 10^{-23}:\\
\;\;\;\;\frac{t}{y \cdot \left(c + y \cdot \left(y \cdot \left(y + a\right) + b\right)\right)}\\
\mathbf{elif}\;y \leq 2.5 \cdot 10^{-50}:\\
\;\;\;\;\frac{t + y \cdot 230661.510616}{i + y \cdot c}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if y < -2.7e25 or 2.49999999999999984e-50 < y Initial program 20.3%
Taylor expanded in y around inf 56.6%
associate--l+56.6%
associate-/l*62.9%
Simplified62.9%
if -2.7e25 < y < -1.99999999999999992e-23Initial program 90.1%
Taylor expanded in t around inf 47.9%
Taylor expanded in i around 0 47.9%
if -1.99999999999999992e-23 < y < 2.49999999999999984e-50Initial program 99.7%
Taylor expanded in y around 0 92.5%
*-commutative92.5%
Simplified92.5%
Taylor expanded in y around 0 82.0%
Final simplification70.8%
(FPCore (x y z t a b c i)
:precision binary64
(let* ((t_1 (+ t (* y 230661.510616))) (t_2 (+ x (- (/ z y) (* a (/ x y))))))
(if (<= y -6.8e+25)
t_2
(if (<= y -3e-30)
(/ t_1 (* y (+ c (* y b))))
(if (<= y 2.5e-50) (/ t_1 (+ i (* y c))) t_2)))))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double t_1 = t + (y * 230661.510616);
double t_2 = x + ((z / y) - (a * (x / y)));
double tmp;
if (y <= -6.8e+25) {
tmp = t_2;
} else if (y <= -3e-30) {
tmp = t_1 / (y * (c + (y * b)));
} else if (y <= 2.5e-50) {
tmp = t_1 / (i + (y * c));
} else {
tmp = t_2;
}
return tmp;
}
real(8) function code(x, y, z, t, a, b, c, i)
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), intent (in) :: i
real(8) :: t_1
real(8) :: t_2
real(8) :: tmp
t_1 = t + (y * 230661.510616d0)
t_2 = x + ((z / y) - (a * (x / y)))
if (y <= (-6.8d+25)) then
tmp = t_2
else if (y <= (-3d-30)) then
tmp = t_1 / (y * (c + (y * b)))
else if (y <= 2.5d-50) then
tmp = t_1 / (i + (y * c))
else
tmp = t_2
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 i) {
double t_1 = t + (y * 230661.510616);
double t_2 = x + ((z / y) - (a * (x / y)));
double tmp;
if (y <= -6.8e+25) {
tmp = t_2;
} else if (y <= -3e-30) {
tmp = t_1 / (y * (c + (y * b)));
} else if (y <= 2.5e-50) {
tmp = t_1 / (i + (y * c));
} else {
tmp = t_2;
}
return tmp;
}
def code(x, y, z, t, a, b, c, i): t_1 = t + (y * 230661.510616) t_2 = x + ((z / y) - (a * (x / y))) tmp = 0 if y <= -6.8e+25: tmp = t_2 elif y <= -3e-30: tmp = t_1 / (y * (c + (y * b))) elif y <= 2.5e-50: tmp = t_1 / (i + (y * c)) else: tmp = t_2 return tmp
function code(x, y, z, t, a, b, c, i) t_1 = Float64(t + Float64(y * 230661.510616)) t_2 = Float64(x + Float64(Float64(z / y) - Float64(a * Float64(x / y)))) tmp = 0.0 if (y <= -6.8e+25) tmp = t_2; elseif (y <= -3e-30) tmp = Float64(t_1 / Float64(y * Float64(c + Float64(y * b)))); elseif (y <= 2.5e-50) tmp = Float64(t_1 / Float64(i + Float64(y * c))); else tmp = t_2; end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c, i) t_1 = t + (y * 230661.510616); t_2 = x + ((z / y) - (a * (x / y))); tmp = 0.0; if (y <= -6.8e+25) tmp = t_2; elseif (y <= -3e-30) tmp = t_1 / (y * (c + (y * b))); elseif (y <= 2.5e-50) tmp = t_1 / (i + (y * c)); else tmp = t_2; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_, i_] := Block[{t$95$1 = N[(t + N[(y * 230661.510616), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(x + N[(N[(z / y), $MachinePrecision] - N[(a * N[(x / y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -6.8e+25], t$95$2, If[LessEqual[y, -3e-30], N[(t$95$1 / N[(y * N[(c + N[(y * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 2.5e-50], N[(t$95$1 / N[(i + N[(y * c), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$2]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := t + y \cdot 230661.510616\\
t_2 := x + \left(\frac{z}{y} - a \cdot \frac{x}{y}\right)\\
\mathbf{if}\;y \leq -6.8 \cdot 10^{+25}:\\
\;\;\;\;t\_2\\
\mathbf{elif}\;y \leq -3 \cdot 10^{-30}:\\
\;\;\;\;\frac{t\_1}{y \cdot \left(c + y \cdot b\right)}\\
\mathbf{elif}\;y \leq 2.5 \cdot 10^{-50}:\\
\;\;\;\;\frac{t\_1}{i + y \cdot c}\\
\mathbf{else}:\\
\;\;\;\;t\_2\\
\end{array}
\end{array}
if y < -6.79999999999999967e25 or 2.49999999999999984e-50 < y Initial program 20.3%
Taylor expanded in y around inf 56.6%
associate--l+56.6%
associate-/l*62.9%
Simplified62.9%
if -6.79999999999999967e25 < y < -2.9999999999999999e-30Initial program 91.0%
Taylor expanded in y around 0 69.0%
*-commutative69.0%
Simplified69.0%
Taylor expanded in y around 0 42.8%
Taylor expanded in i around 0 43.2%
if -2.9999999999999999e-30 < y < 2.49999999999999984e-50Initial program 99.7%
Taylor expanded in y around 0 92.4%
*-commutative92.4%
Simplified92.4%
Taylor expanded in y around 0 82.7%
Final simplification70.8%
(FPCore (x y z t a b c i) :precision binary64 (if (or (<= y -3.5e+26) (not (<= y 4.6e-23))) (+ x (- (/ z y) (* a (/ x y)))) (/ (+ t (* y 230661.510616)) (+ i (* y (+ c (* y b)))))))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double tmp;
if ((y <= -3.5e+26) || !(y <= 4.6e-23)) {
tmp = x + ((z / y) - (a * (x / y)));
} else {
tmp = (t + (y * 230661.510616)) / (i + (y * (c + (y * b))));
}
return tmp;
}
real(8) function code(x, y, z, t, a, b, c, i)
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), intent (in) :: i
real(8) :: tmp
if ((y <= (-3.5d+26)) .or. (.not. (y <= 4.6d-23))) then
tmp = x + ((z / y) - (a * (x / y)))
else
tmp = (t + (y * 230661.510616d0)) / (i + (y * (c + (y * 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 i) {
double tmp;
if ((y <= -3.5e+26) || !(y <= 4.6e-23)) {
tmp = x + ((z / y) - (a * (x / y)));
} else {
tmp = (t + (y * 230661.510616)) / (i + (y * (c + (y * b))));
}
return tmp;
}
def code(x, y, z, t, a, b, c, i): tmp = 0 if (y <= -3.5e+26) or not (y <= 4.6e-23): tmp = x + ((z / y) - (a * (x / y))) else: tmp = (t + (y * 230661.510616)) / (i + (y * (c + (y * b)))) return tmp
function code(x, y, z, t, a, b, c, i) tmp = 0.0 if ((y <= -3.5e+26) || !(y <= 4.6e-23)) tmp = Float64(x + Float64(Float64(z / y) - Float64(a * Float64(x / y)))); else tmp = Float64(Float64(t + Float64(y * 230661.510616)) / Float64(i + Float64(y * Float64(c + Float64(y * b))))); end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c, i) tmp = 0.0; if ((y <= -3.5e+26) || ~((y <= 4.6e-23))) tmp = x + ((z / y) - (a * (x / y))); else tmp = (t + (y * 230661.510616)) / (i + (y * (c + (y * b)))); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_, i_] := If[Or[LessEqual[y, -3.5e+26], N[Not[LessEqual[y, 4.6e-23]], $MachinePrecision]], N[(x + N[(N[(z / y), $MachinePrecision] - N[(a * N[(x / y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(t + N[(y * 230661.510616), $MachinePrecision]), $MachinePrecision] / N[(i + N[(y * N[(c + N[(y * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -3.5 \cdot 10^{+26} \lor \neg \left(y \leq 4.6 \cdot 10^{-23}\right):\\
\;\;\;\;x + \left(\frac{z}{y} - a \cdot \frac{x}{y}\right)\\
\mathbf{else}:\\
\;\;\;\;\frac{t + y \cdot 230661.510616}{i + y \cdot \left(c + y \cdot b\right)}\\
\end{array}
\end{array}
if y < -3.4999999999999999e26 or 4.6000000000000002e-23 < y Initial program 15.2%
Taylor expanded in y around inf 59.3%
associate--l+59.3%
associate-/l*65.9%
Simplified65.9%
if -3.4999999999999999e26 < y < 4.6000000000000002e-23Initial program 99.0%
Taylor expanded in y around 0 88.0%
*-commutative88.0%
Simplified88.0%
Taylor expanded in y around 0 83.2%
Final simplification74.9%
(FPCore (x y z t a b c i) :precision binary64 (if (or (<= y -1.6e+26) (not (<= y 2.5e-50))) (+ x (- (/ z y) (* a (/ x y)))) (/ (+ t (* y 230661.510616)) (+ i (* y c)))))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double tmp;
if ((y <= -1.6e+26) || !(y <= 2.5e-50)) {
tmp = x + ((z / y) - (a * (x / y)));
} else {
tmp = (t + (y * 230661.510616)) / (i + (y * c));
}
return tmp;
}
real(8) function code(x, y, z, t, a, b, c, i)
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), intent (in) :: i
real(8) :: tmp
if ((y <= (-1.6d+26)) .or. (.not. (y <= 2.5d-50))) then
tmp = x + ((z / y) - (a * (x / y)))
else
tmp = (t + (y * 230661.510616d0)) / (i + (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 i) {
double tmp;
if ((y <= -1.6e+26) || !(y <= 2.5e-50)) {
tmp = x + ((z / y) - (a * (x / y)));
} else {
tmp = (t + (y * 230661.510616)) / (i + (y * c));
}
return tmp;
}
def code(x, y, z, t, a, b, c, i): tmp = 0 if (y <= -1.6e+26) or not (y <= 2.5e-50): tmp = x + ((z / y) - (a * (x / y))) else: tmp = (t + (y * 230661.510616)) / (i + (y * c)) return tmp
function code(x, y, z, t, a, b, c, i) tmp = 0.0 if ((y <= -1.6e+26) || !(y <= 2.5e-50)) tmp = Float64(x + Float64(Float64(z / y) - Float64(a * Float64(x / y)))); else tmp = Float64(Float64(t + Float64(y * 230661.510616)) / Float64(i + Float64(y * c))); end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c, i) tmp = 0.0; if ((y <= -1.6e+26) || ~((y <= 2.5e-50))) tmp = x + ((z / y) - (a * (x / y))); else tmp = (t + (y * 230661.510616)) / (i + (y * c)); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_, i_] := If[Or[LessEqual[y, -1.6e+26], N[Not[LessEqual[y, 2.5e-50]], $MachinePrecision]], N[(x + N[(N[(z / y), $MachinePrecision] - N[(a * N[(x / y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(t + N[(y * 230661.510616), $MachinePrecision]), $MachinePrecision] / N[(i + N[(y * c), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -1.6 \cdot 10^{+26} \lor \neg \left(y \leq 2.5 \cdot 10^{-50}\right):\\
\;\;\;\;x + \left(\frac{z}{y} - a \cdot \frac{x}{y}\right)\\
\mathbf{else}:\\
\;\;\;\;\frac{t + y \cdot 230661.510616}{i + y \cdot c}\\
\end{array}
\end{array}
if y < -1.60000000000000014e26 or 2.49999999999999984e-50 < y Initial program 20.3%
Taylor expanded in y around inf 56.6%
associate--l+56.6%
associate-/l*62.9%
Simplified62.9%
if -1.60000000000000014e26 < y < 2.49999999999999984e-50Initial program 98.9%
Taylor expanded in y around 0 90.3%
*-commutative90.3%
Simplified90.3%
Taylor expanded in y around 0 75.7%
Final simplification69.1%
(FPCore (x y z t a b c i) :precision binary64 (if (or (<= y -1.15e+26) (not (<= y 2.5e-50))) (+ x (- (/ z y) (* a (/ x y)))) (/ t (+ i (* y (+ c (* y b)))))))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double tmp;
if ((y <= -1.15e+26) || !(y <= 2.5e-50)) {
tmp = x + ((z / y) - (a * (x / y)));
} else {
tmp = t / (i + (y * (c + (y * b))));
}
return tmp;
}
real(8) function code(x, y, z, t, a, b, c, i)
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), intent (in) :: i
real(8) :: tmp
if ((y <= (-1.15d+26)) .or. (.not. (y <= 2.5d-50))) then
tmp = x + ((z / y) - (a * (x / y)))
else
tmp = t / (i + (y * (c + (y * 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 i) {
double tmp;
if ((y <= -1.15e+26) || !(y <= 2.5e-50)) {
tmp = x + ((z / y) - (a * (x / y)));
} else {
tmp = t / (i + (y * (c + (y * b))));
}
return tmp;
}
def code(x, y, z, t, a, b, c, i): tmp = 0 if (y <= -1.15e+26) or not (y <= 2.5e-50): tmp = x + ((z / y) - (a * (x / y))) else: tmp = t / (i + (y * (c + (y * b)))) return tmp
function code(x, y, z, t, a, b, c, i) tmp = 0.0 if ((y <= -1.15e+26) || !(y <= 2.5e-50)) tmp = Float64(x + Float64(Float64(z / y) - Float64(a * Float64(x / y)))); else tmp = Float64(t / Float64(i + Float64(y * Float64(c + Float64(y * b))))); end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c, i) tmp = 0.0; if ((y <= -1.15e+26) || ~((y <= 2.5e-50))) tmp = x + ((z / y) - (a * (x / y))); else tmp = t / (i + (y * (c + (y * b)))); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_, i_] := If[Or[LessEqual[y, -1.15e+26], N[Not[LessEqual[y, 2.5e-50]], $MachinePrecision]], N[(x + N[(N[(z / y), $MachinePrecision] - N[(a * N[(x / y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(t / N[(i + N[(y * N[(c + N[(y * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -1.15 \cdot 10^{+26} \lor \neg \left(y \leq 2.5 \cdot 10^{-50}\right):\\
\;\;\;\;x + \left(\frac{z}{y} - a \cdot \frac{x}{y}\right)\\
\mathbf{else}:\\
\;\;\;\;\frac{t}{i + y \cdot \left(c + y \cdot b\right)}\\
\end{array}
\end{array}
if y < -1.15e26 or 2.49999999999999984e-50 < y Initial program 20.3%
Taylor expanded in y around inf 56.6%
associate--l+56.6%
associate-/l*62.9%
Simplified62.9%
if -1.15e26 < y < 2.49999999999999984e-50Initial program 98.9%
Taylor expanded in t around inf 72.6%
Taylor expanded in y around 0 70.3%
*-commutative70.3%
Simplified70.3%
Final simplification66.5%
(FPCore (x y z t a b c i) :precision binary64 (if (<= y -3.4e+63) x (if (<= y 3.5e-70) (/ t (+ i (* y c))) x)))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double tmp;
if (y <= -3.4e+63) {
tmp = x;
} else if (y <= 3.5e-70) {
tmp = t / (i + (y * c));
} else {
tmp = x;
}
return tmp;
}
real(8) function code(x, y, z, t, a, b, c, i)
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), intent (in) :: i
real(8) :: tmp
if (y <= (-3.4d+63)) then
tmp = x
else if (y <= 3.5d-70) then
tmp = t / (i + (y * c))
else
tmp = 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 i) {
double tmp;
if (y <= -3.4e+63) {
tmp = x;
} else if (y <= 3.5e-70) {
tmp = t / (i + (y * c));
} else {
tmp = x;
}
return tmp;
}
def code(x, y, z, t, a, b, c, i): tmp = 0 if y <= -3.4e+63: tmp = x elif y <= 3.5e-70: tmp = t / (i + (y * c)) else: tmp = x return tmp
function code(x, y, z, t, a, b, c, i) tmp = 0.0 if (y <= -3.4e+63) tmp = x; elseif (y <= 3.5e-70) tmp = Float64(t / Float64(i + Float64(y * c))); else tmp = x; end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c, i) tmp = 0.0; if (y <= -3.4e+63) tmp = x; elseif (y <= 3.5e-70) tmp = t / (i + (y * c)); else tmp = x; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_, i_] := If[LessEqual[y, -3.4e+63], x, If[LessEqual[y, 3.5e-70], N[(t / N[(i + N[(y * c), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], x]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -3.4 \cdot 10^{+63}:\\
\;\;\;\;x\\
\mathbf{elif}\;y \leq 3.5 \cdot 10^{-70}:\\
\;\;\;\;\frac{t}{i + y \cdot c}\\
\mathbf{else}:\\
\;\;\;\;x\\
\end{array}
\end{array}
if y < -3.3999999999999999e63 or 3.49999999999999974e-70 < y Initial program 21.6%
Taylor expanded in y around inf 43.5%
if -3.3999999999999999e63 < y < 3.49999999999999974e-70Initial program 98.2%
Taylor expanded in t around inf 71.7%
Taylor expanded in y around 0 62.4%
*-commutative62.4%
Simplified62.4%
(FPCore (x y z t a b c i) :precision binary64 (if (<= y -3.4e+63) x (if (<= y 4.2e-75) (/ t i) x)))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double tmp;
if (y <= -3.4e+63) {
tmp = x;
} else if (y <= 4.2e-75) {
tmp = t / i;
} else {
tmp = x;
}
return tmp;
}
real(8) function code(x, y, z, t, a, b, c, i)
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), intent (in) :: i
real(8) :: tmp
if (y <= (-3.4d+63)) then
tmp = x
else if (y <= 4.2d-75) then
tmp = t / i
else
tmp = 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 i) {
double tmp;
if (y <= -3.4e+63) {
tmp = x;
} else if (y <= 4.2e-75) {
tmp = t / i;
} else {
tmp = x;
}
return tmp;
}
def code(x, y, z, t, a, b, c, i): tmp = 0 if y <= -3.4e+63: tmp = x elif y <= 4.2e-75: tmp = t / i else: tmp = x return tmp
function code(x, y, z, t, a, b, c, i) tmp = 0.0 if (y <= -3.4e+63) tmp = x; elseif (y <= 4.2e-75) tmp = Float64(t / i); else tmp = x; end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c, i) tmp = 0.0; if (y <= -3.4e+63) tmp = x; elseif (y <= 4.2e-75) tmp = t / i; else tmp = x; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_, i_] := If[LessEqual[y, -3.4e+63], x, If[LessEqual[y, 4.2e-75], N[(t / i), $MachinePrecision], x]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -3.4 \cdot 10^{+63}:\\
\;\;\;\;x\\
\mathbf{elif}\;y \leq 4.2 \cdot 10^{-75}:\\
\;\;\;\;\frac{t}{i}\\
\mathbf{else}:\\
\;\;\;\;x\\
\end{array}
\end{array}
if y < -3.3999999999999999e63 or 4.2000000000000002e-75 < y Initial program 23.3%
Taylor expanded in y around inf 42.6%
if -3.3999999999999999e63 < y < 4.2000000000000002e-75Initial program 98.2%
Taylor expanded in y around 0 50.3%
(FPCore (x y z t a b c i) :precision binary64 x)
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
return x;
}
real(8) function code(x, y, z, t, a, b, c, i)
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), intent (in) :: i
code = x
end function
public static double code(double x, double y, double z, double t, double a, double b, double c, double i) {
return x;
}
def code(x, y, z, t, a, b, c, i): return x
function code(x, y, z, t, a, b, c, i) return x end
function tmp = code(x, y, z, t, a, b, c, i) tmp = x; end
code[x_, y_, z_, t_, a_, b_, c_, i_] := x
\begin{array}{l}
\\
x
\end{array}
Initial program 58.4%
Taylor expanded in y around inf 24.7%
herbie shell --seed 2024105
(FPCore (x y z t a b c i)
:name "Numeric.SpecFunctions:logGamma from math-functions-0.1.5.2"
:precision binary64
(/ (+ (* (+ (* (+ (* (+ (* x y) z) y) 27464.7644705) y) 230661.510616) y) t) (+ (* (+ (* (+ (* (+ y a) y) b) y) c) y) i)))