
(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 21 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 (+ 1.0 (/ a y)))
(t_2 (+ (* y (+ (* y (+ (* x y) z)) 27464.7644705)) 230661.510616))
(t_3 (* y t_2))
(t_4 (+ t t_3))
(t_5 (+ (* y (+ y a)) b))
(t_6 (+ i (* y (+ c (* y t_5)))))
(t_7 (/ t_4 t_6))
(t_8 (* x t_6)))
(if (<= t_7 (- INFINITY))
(*
x
(+
(/ t t_8)
(+
(/ (* y (+ 230661.510616 (* y (+ 27464.7644705 (* y z))))) t_8)
(/ (pow y 4.0) t_6))))
(if (<= t_7 -2e-320)
t_7
(if (<= t_7 0.0)
(/
1.0
(* c (+ (/ 1.0 t_2) (+ (/ i (* c t_3)) (* (/ y c) (/ t_5 t_2))))))
(if (<= t_7 INFINITY)
(/ t_4 (+ i (+ (* y c) (* t_5 (pow y 2.0)))))
(* x (+ (/ 1.0 t_1) (/ z (* x (* y (pow t_1 2.0))))))))))))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double t_1 = 1.0 + (a / y);
double t_2 = (y * ((y * ((x * y) + z)) + 27464.7644705)) + 230661.510616;
double t_3 = y * t_2;
double t_4 = t + t_3;
double t_5 = (y * (y + a)) + b;
double t_6 = i + (y * (c + (y * t_5)));
double t_7 = t_4 / t_6;
double t_8 = x * t_6;
double tmp;
if (t_7 <= -((double) INFINITY)) {
tmp = x * ((t / t_8) + (((y * (230661.510616 + (y * (27464.7644705 + (y * z))))) / t_8) + (pow(y, 4.0) / t_6)));
} else if (t_7 <= -2e-320) {
tmp = t_7;
} else if (t_7 <= 0.0) {
tmp = 1.0 / (c * ((1.0 / t_2) + ((i / (c * t_3)) + ((y / c) * (t_5 / t_2)))));
} else if (t_7 <= ((double) INFINITY)) {
tmp = t_4 / (i + ((y * c) + (t_5 * pow(y, 2.0))));
} else {
tmp = x * ((1.0 / t_1) + (z / (x * (y * pow(t_1, 2.0)))));
}
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 = 1.0 + (a / y);
double t_2 = (y * ((y * ((x * y) + z)) + 27464.7644705)) + 230661.510616;
double t_3 = y * t_2;
double t_4 = t + t_3;
double t_5 = (y * (y + a)) + b;
double t_6 = i + (y * (c + (y * t_5)));
double t_7 = t_4 / t_6;
double t_8 = x * t_6;
double tmp;
if (t_7 <= -Double.POSITIVE_INFINITY) {
tmp = x * ((t / t_8) + (((y * (230661.510616 + (y * (27464.7644705 + (y * z))))) / t_8) + (Math.pow(y, 4.0) / t_6)));
} else if (t_7 <= -2e-320) {
tmp = t_7;
} else if (t_7 <= 0.0) {
tmp = 1.0 / (c * ((1.0 / t_2) + ((i / (c * t_3)) + ((y / c) * (t_5 / t_2)))));
} else if (t_7 <= Double.POSITIVE_INFINITY) {
tmp = t_4 / (i + ((y * c) + (t_5 * Math.pow(y, 2.0))));
} else {
tmp = x * ((1.0 / t_1) + (z / (x * (y * Math.pow(t_1, 2.0)))));
}
return tmp;
}
def code(x, y, z, t, a, b, c, i): t_1 = 1.0 + (a / y) t_2 = (y * ((y * ((x * y) + z)) + 27464.7644705)) + 230661.510616 t_3 = y * t_2 t_4 = t + t_3 t_5 = (y * (y + a)) + b t_6 = i + (y * (c + (y * t_5))) t_7 = t_4 / t_6 t_8 = x * t_6 tmp = 0 if t_7 <= -math.inf: tmp = x * ((t / t_8) + (((y * (230661.510616 + (y * (27464.7644705 + (y * z))))) / t_8) + (math.pow(y, 4.0) / t_6))) elif t_7 <= -2e-320: tmp = t_7 elif t_7 <= 0.0: tmp = 1.0 / (c * ((1.0 / t_2) + ((i / (c * t_3)) + ((y / c) * (t_5 / t_2))))) elif t_7 <= math.inf: tmp = t_4 / (i + ((y * c) + (t_5 * math.pow(y, 2.0)))) else: tmp = x * ((1.0 / t_1) + (z / (x * (y * math.pow(t_1, 2.0))))) return tmp
function code(x, y, z, t, a, b, c, i) t_1 = Float64(1.0 + Float64(a / y)) t_2 = Float64(Float64(y * Float64(Float64(y * Float64(Float64(x * y) + z)) + 27464.7644705)) + 230661.510616) t_3 = Float64(y * t_2) t_4 = Float64(t + t_3) t_5 = Float64(Float64(y * Float64(y + a)) + b) t_6 = Float64(i + Float64(y * Float64(c + Float64(y * t_5)))) t_7 = Float64(t_4 / t_6) t_8 = Float64(x * t_6) tmp = 0.0 if (t_7 <= Float64(-Inf)) tmp = Float64(x * Float64(Float64(t / t_8) + Float64(Float64(Float64(y * Float64(230661.510616 + Float64(y * Float64(27464.7644705 + Float64(y * z))))) / t_8) + Float64((y ^ 4.0) / t_6)))); elseif (t_7 <= -2e-320) tmp = t_7; elseif (t_7 <= 0.0) tmp = Float64(1.0 / Float64(c * Float64(Float64(1.0 / t_2) + Float64(Float64(i / Float64(c * t_3)) + Float64(Float64(y / c) * Float64(t_5 / t_2)))))); elseif (t_7 <= Inf) tmp = Float64(t_4 / Float64(i + Float64(Float64(y * c) + Float64(t_5 * (y ^ 2.0))))); else tmp = Float64(x * Float64(Float64(1.0 / t_1) + Float64(z / Float64(x * Float64(y * (t_1 ^ 2.0)))))); end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c, i) t_1 = 1.0 + (a / y); t_2 = (y * ((y * ((x * y) + z)) + 27464.7644705)) + 230661.510616; t_3 = y * t_2; t_4 = t + t_3; t_5 = (y * (y + a)) + b; t_6 = i + (y * (c + (y * t_5))); t_7 = t_4 / t_6; t_8 = x * t_6; tmp = 0.0; if (t_7 <= -Inf) tmp = x * ((t / t_8) + (((y * (230661.510616 + (y * (27464.7644705 + (y * z))))) / t_8) + ((y ^ 4.0) / t_6))); elseif (t_7 <= -2e-320) tmp = t_7; elseif (t_7 <= 0.0) tmp = 1.0 / (c * ((1.0 / t_2) + ((i / (c * t_3)) + ((y / c) * (t_5 / t_2))))); elseif (t_7 <= Inf) tmp = t_4 / (i + ((y * c) + (t_5 * (y ^ 2.0)))); else tmp = x * ((1.0 / t_1) + (z / (x * (y * (t_1 ^ 2.0))))); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_, i_] := Block[{t$95$1 = N[(1.0 + N[(a / y), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(N[(y * N[(N[(y * N[(N[(x * y), $MachinePrecision] + z), $MachinePrecision]), $MachinePrecision] + 27464.7644705), $MachinePrecision]), $MachinePrecision] + 230661.510616), $MachinePrecision]}, Block[{t$95$3 = N[(y * t$95$2), $MachinePrecision]}, Block[{t$95$4 = N[(t + t$95$3), $MachinePrecision]}, Block[{t$95$5 = N[(N[(y * N[(y + a), $MachinePrecision]), $MachinePrecision] + b), $MachinePrecision]}, Block[{t$95$6 = N[(i + N[(y * N[(c + N[(y * t$95$5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$7 = N[(t$95$4 / t$95$6), $MachinePrecision]}, Block[{t$95$8 = N[(x * t$95$6), $MachinePrecision]}, If[LessEqual[t$95$7, (-Infinity)], N[(x * N[(N[(t / t$95$8), $MachinePrecision] + N[(N[(N[(y * N[(230661.510616 + N[(y * N[(27464.7644705 + N[(y * z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / t$95$8), $MachinePrecision] + N[(N[Power[y, 4.0], $MachinePrecision] / t$95$6), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$7, -2e-320], t$95$7, If[LessEqual[t$95$7, 0.0], N[(1.0 / N[(c * N[(N[(1.0 / t$95$2), $MachinePrecision] + N[(N[(i / N[(c * t$95$3), $MachinePrecision]), $MachinePrecision] + N[(N[(y / c), $MachinePrecision] * N[(t$95$5 / t$95$2), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$7, Infinity], N[(t$95$4 / N[(i + N[(N[(y * c), $MachinePrecision] + N[(t$95$5 * N[Power[y, 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x * N[(N[(1.0 / t$95$1), $MachinePrecision] + N[(z / N[(x * N[(y * N[Power[t$95$1, 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]]]]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := 1 + \frac{a}{y}\\
t_2 := y \cdot \left(y \cdot \left(x \cdot y + z\right) + 27464.7644705\right) + 230661.510616\\
t_3 := y \cdot t\_2\\
t_4 := t + t\_3\\
t_5 := y \cdot \left(y + a\right) + b\\
t_6 := i + y \cdot \left(c + y \cdot t\_5\right)\\
t_7 := \frac{t\_4}{t\_6}\\
t_8 := x \cdot t\_6\\
\mathbf{if}\;t\_7 \leq -\infty:\\
\;\;\;\;x \cdot \left(\frac{t}{t\_8} + \left(\frac{y \cdot \left(230661.510616 + y \cdot \left(27464.7644705 + y \cdot z\right)\right)}{t\_8} + \frac{{y}^{4}}{t\_6}\right)\right)\\
\mathbf{elif}\;t\_7 \leq -2 \cdot 10^{-320}:\\
\;\;\;\;t\_7\\
\mathbf{elif}\;t\_7 \leq 0:\\
\;\;\;\;\frac{1}{c \cdot \left(\frac{1}{t\_2} + \left(\frac{i}{c \cdot t\_3} + \frac{y}{c} \cdot \frac{t\_5}{t\_2}\right)\right)}\\
\mathbf{elif}\;t\_7 \leq \infty:\\
\;\;\;\;\frac{t\_4}{i + \left(y \cdot c + t\_5 \cdot {y}^{2}\right)}\\
\mathbf{else}:\\
\;\;\;\;x \cdot \left(\frac{1}{t\_1} + \frac{z}{x \cdot \left(y \cdot {t\_1}^{2}\right)}\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 26.2%
fma-define26.2%
fma-define26.2%
fma-define26.2%
fma-define26.2%
fma-define26.2%
fma-define26.2%
fma-define26.2%
Simplified26.2%
Taylor expanded in x around inf 89.2%
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)) < -1.99998e-320Initial program 99.6%
if -1.99998e-320 < (/.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)) < 0.0Initial program 48.0%
fma-define48.0%
fma-define48.0%
fma-define48.0%
fma-define48.0%
fma-define48.0%
fma-define48.0%
fma-define48.0%
Simplified48.0%
Applied egg-rr48.0%
unpow-148.0%
+-commutative48.0%
Simplified48.0%
Taylor expanded in t around 0 48.0%
associate-/r*48.0%
Simplified48.0%
Taylor expanded in c around inf 76.5%
times-frac94.8%
Simplified94.8%
if 0.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)) < +inf.0Initial program 97.0%
Taylor expanded in c around 0 97.1%
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)) Initial program 0.0%
fma-define0.0%
fma-define0.0%
fma-define0.0%
fma-define0.0%
fma-define0.0%
fma-define0.0%
fma-define0.0%
Simplified0.0%
Applied egg-rr0.0%
unpow-10.0%
+-commutative0.0%
Simplified0.0%
Taylor expanded in y around inf 56.3%
Taylor expanded in x around inf 86.6%
Final simplification92.9%
(FPCore (x y z t a b c i)
:precision binary64
(let* ((t_1 (+ 1.0 (/ a y)))
(t_2 (+ (* y (+ (* y (+ (* x y) z)) 27464.7644705)) 230661.510616))
(t_3 (* y t_2))
(t_4 (+ t t_3))
(t_5 (+ (* y (+ y a)) b))
(t_6 (+ i (* y (+ c (* y t_5)))))
(t_7 (/ t_4 t_6)))
(if (<= t_7 -2e-320)
(* t (+ (/ 1.0 t_6) (* (/ y t) (/ t_2 t_6))))
(if (<= t_7 0.0)
(/
1.0
(* c (+ (/ 1.0 t_2) (+ (/ i (* c t_3)) (* (/ y c) (/ t_5 t_2))))))
(if (<= t_7 INFINITY)
(/ t_4 (+ i (+ (* y c) (* t_5 (pow y 2.0)))))
(* x (+ (/ 1.0 t_1) (/ z (* x (* y (pow t_1 2.0)))))))))))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double t_1 = 1.0 + (a / y);
double t_2 = (y * ((y * ((x * y) + z)) + 27464.7644705)) + 230661.510616;
double t_3 = y * t_2;
double t_4 = t + t_3;
double t_5 = (y * (y + a)) + b;
double t_6 = i + (y * (c + (y * t_5)));
double t_7 = t_4 / t_6;
double tmp;
if (t_7 <= -2e-320) {
tmp = t * ((1.0 / t_6) + ((y / t) * (t_2 / t_6)));
} else if (t_7 <= 0.0) {
tmp = 1.0 / (c * ((1.0 / t_2) + ((i / (c * t_3)) + ((y / c) * (t_5 / t_2)))));
} else if (t_7 <= ((double) INFINITY)) {
tmp = t_4 / (i + ((y * c) + (t_5 * pow(y, 2.0))));
} else {
tmp = x * ((1.0 / t_1) + (z / (x * (y * pow(t_1, 2.0)))));
}
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 = 1.0 + (a / y);
double t_2 = (y * ((y * ((x * y) + z)) + 27464.7644705)) + 230661.510616;
double t_3 = y * t_2;
double t_4 = t + t_3;
double t_5 = (y * (y + a)) + b;
double t_6 = i + (y * (c + (y * t_5)));
double t_7 = t_4 / t_6;
double tmp;
if (t_7 <= -2e-320) {
tmp = t * ((1.0 / t_6) + ((y / t) * (t_2 / t_6)));
} else if (t_7 <= 0.0) {
tmp = 1.0 / (c * ((1.0 / t_2) + ((i / (c * t_3)) + ((y / c) * (t_5 / t_2)))));
} else if (t_7 <= Double.POSITIVE_INFINITY) {
tmp = t_4 / (i + ((y * c) + (t_5 * Math.pow(y, 2.0))));
} else {
tmp = x * ((1.0 / t_1) + (z / (x * (y * Math.pow(t_1, 2.0)))));
}
return tmp;
}
def code(x, y, z, t, a, b, c, i): t_1 = 1.0 + (a / y) t_2 = (y * ((y * ((x * y) + z)) + 27464.7644705)) + 230661.510616 t_3 = y * t_2 t_4 = t + t_3 t_5 = (y * (y + a)) + b t_6 = i + (y * (c + (y * t_5))) t_7 = t_4 / t_6 tmp = 0 if t_7 <= -2e-320: tmp = t * ((1.0 / t_6) + ((y / t) * (t_2 / t_6))) elif t_7 <= 0.0: tmp = 1.0 / (c * ((1.0 / t_2) + ((i / (c * t_3)) + ((y / c) * (t_5 / t_2))))) elif t_7 <= math.inf: tmp = t_4 / (i + ((y * c) + (t_5 * math.pow(y, 2.0)))) else: tmp = x * ((1.0 / t_1) + (z / (x * (y * math.pow(t_1, 2.0))))) return tmp
function code(x, y, z, t, a, b, c, i) t_1 = Float64(1.0 + Float64(a / y)) t_2 = Float64(Float64(y * Float64(Float64(y * Float64(Float64(x * y) + z)) + 27464.7644705)) + 230661.510616) t_3 = Float64(y * t_2) t_4 = Float64(t + t_3) t_5 = Float64(Float64(y * Float64(y + a)) + b) t_6 = Float64(i + Float64(y * Float64(c + Float64(y * t_5)))) t_7 = Float64(t_4 / t_6) tmp = 0.0 if (t_7 <= -2e-320) tmp = Float64(t * Float64(Float64(1.0 / t_6) + Float64(Float64(y / t) * Float64(t_2 / t_6)))); elseif (t_7 <= 0.0) tmp = Float64(1.0 / Float64(c * Float64(Float64(1.0 / t_2) + Float64(Float64(i / Float64(c * t_3)) + Float64(Float64(y / c) * Float64(t_5 / t_2)))))); elseif (t_7 <= Inf) tmp = Float64(t_4 / Float64(i + Float64(Float64(y * c) + Float64(t_5 * (y ^ 2.0))))); else tmp = Float64(x * Float64(Float64(1.0 / t_1) + Float64(z / Float64(x * Float64(y * (t_1 ^ 2.0)))))); end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c, i) t_1 = 1.0 + (a / y); t_2 = (y * ((y * ((x * y) + z)) + 27464.7644705)) + 230661.510616; t_3 = y * t_2; t_4 = t + t_3; t_5 = (y * (y + a)) + b; t_6 = i + (y * (c + (y * t_5))); t_7 = t_4 / t_6; tmp = 0.0; if (t_7 <= -2e-320) tmp = t * ((1.0 / t_6) + ((y / t) * (t_2 / t_6))); elseif (t_7 <= 0.0) tmp = 1.0 / (c * ((1.0 / t_2) + ((i / (c * t_3)) + ((y / c) * (t_5 / t_2))))); elseif (t_7 <= Inf) tmp = t_4 / (i + ((y * c) + (t_5 * (y ^ 2.0)))); else tmp = x * ((1.0 / t_1) + (z / (x * (y * (t_1 ^ 2.0))))); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_, i_] := Block[{t$95$1 = N[(1.0 + N[(a / y), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(N[(y * N[(N[(y * N[(N[(x * y), $MachinePrecision] + z), $MachinePrecision]), $MachinePrecision] + 27464.7644705), $MachinePrecision]), $MachinePrecision] + 230661.510616), $MachinePrecision]}, Block[{t$95$3 = N[(y * t$95$2), $MachinePrecision]}, Block[{t$95$4 = N[(t + t$95$3), $MachinePrecision]}, Block[{t$95$5 = N[(N[(y * N[(y + a), $MachinePrecision]), $MachinePrecision] + b), $MachinePrecision]}, Block[{t$95$6 = N[(i + N[(y * N[(c + N[(y * t$95$5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$7 = N[(t$95$4 / t$95$6), $MachinePrecision]}, If[LessEqual[t$95$7, -2e-320], N[(t * N[(N[(1.0 / t$95$6), $MachinePrecision] + N[(N[(y / t), $MachinePrecision] * N[(t$95$2 / t$95$6), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$7, 0.0], N[(1.0 / N[(c * N[(N[(1.0 / t$95$2), $MachinePrecision] + N[(N[(i / N[(c * t$95$3), $MachinePrecision]), $MachinePrecision] + N[(N[(y / c), $MachinePrecision] * N[(t$95$5 / t$95$2), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$7, Infinity], N[(t$95$4 / N[(i + N[(N[(y * c), $MachinePrecision] + N[(t$95$5 * N[Power[y, 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x * N[(N[(1.0 / t$95$1), $MachinePrecision] + N[(z / N[(x * N[(y * N[Power[t$95$1, 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := 1 + \frac{a}{y}\\
t_2 := y \cdot \left(y \cdot \left(x \cdot y + z\right) + 27464.7644705\right) + 230661.510616\\
t_3 := y \cdot t\_2\\
t_4 := t + t\_3\\
t_5 := y \cdot \left(y + a\right) + b\\
t_6 := i + y \cdot \left(c + y \cdot t\_5\right)\\
t_7 := \frac{t\_4}{t\_6}\\
\mathbf{if}\;t\_7 \leq -2 \cdot 10^{-320}:\\
\;\;\;\;t \cdot \left(\frac{1}{t\_6} + \frac{y}{t} \cdot \frac{t\_2}{t\_6}\right)\\
\mathbf{elif}\;t\_7 \leq 0:\\
\;\;\;\;\frac{1}{c \cdot \left(\frac{1}{t\_2} + \left(\frac{i}{c \cdot t\_3} + \frac{y}{c} \cdot \frac{t\_5}{t\_2}\right)\right)}\\
\mathbf{elif}\;t\_7 \leq \infty:\\
\;\;\;\;\frac{t\_4}{i + \left(y \cdot c + t\_5 \cdot {y}^{2}\right)}\\
\mathbf{else}:\\
\;\;\;\;x \cdot \left(\frac{1}{t\_1} + \frac{z}{x \cdot \left(y \cdot {t\_1}^{2}\right)}\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)) < -1.99998e-320Initial program 87.8%
fma-define87.8%
fma-define87.8%
fma-define87.8%
fma-define87.8%
fma-define87.8%
fma-define87.8%
fma-define87.8%
Simplified87.8%
Applied egg-rr87.3%
unpow-187.3%
+-commutative87.3%
Simplified87.3%
Taylor expanded in t around inf 80.7%
times-frac89.3%
Simplified89.3%
if -1.99998e-320 < (/.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)) < 0.0Initial program 48.0%
fma-define48.0%
fma-define48.0%
fma-define48.0%
fma-define48.0%
fma-define48.0%
fma-define48.0%
fma-define48.0%
Simplified48.0%
Applied egg-rr48.0%
unpow-148.0%
+-commutative48.0%
Simplified48.0%
Taylor expanded in t around 0 48.0%
associate-/r*48.0%
Simplified48.0%
Taylor expanded in c around inf 76.5%
times-frac94.8%
Simplified94.8%
if 0.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)) < +inf.0Initial program 97.0%
Taylor expanded in c around 0 97.1%
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)) Initial program 0.0%
fma-define0.0%
fma-define0.0%
fma-define0.0%
fma-define0.0%
fma-define0.0%
fma-define0.0%
fma-define0.0%
Simplified0.0%
Applied egg-rr0.0%
unpow-10.0%
+-commutative0.0%
Simplified0.0%
Taylor expanded in y around inf 56.3%
Taylor expanded in x around inf 86.6%
Final simplification91.0%
(FPCore (x y z t a b c i)
:precision binary64
(let* ((t_1 (+ (* y (+ (* y (+ (* x y) z)) 27464.7644705)) 230661.510616))
(t_2 (* y t_1))
(t_3 (+ (* y (+ y a)) b))
(t_4 (+ i (* y (+ c (* y t_3)))))
(t_5 (/ (+ t t_2) t_4))
(t_6 (+ 1.0 (/ a y))))
(if (<= t_5 -2e-320)
(* t (+ (/ 1.0 t_4) (* (/ y t) (/ t_1 t_4))))
(if (<= t_5 0.0)
(/
1.0
(* c (+ (/ 1.0 t_1) (+ (/ i (* c t_2)) (* (/ y c) (/ t_3 t_1))))))
(if (<= t_5 INFINITY)
t_5
(* x (+ (/ 1.0 t_6) (/ z (* x (* y (pow t_6 2.0)))))))))))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double t_1 = (y * ((y * ((x * y) + z)) + 27464.7644705)) + 230661.510616;
double t_2 = y * t_1;
double t_3 = (y * (y + a)) + b;
double t_4 = i + (y * (c + (y * t_3)));
double t_5 = (t + t_2) / t_4;
double t_6 = 1.0 + (a / y);
double tmp;
if (t_5 <= -2e-320) {
tmp = t * ((1.0 / t_4) + ((y / t) * (t_1 / t_4)));
} else if (t_5 <= 0.0) {
tmp = 1.0 / (c * ((1.0 / t_1) + ((i / (c * t_2)) + ((y / c) * (t_3 / t_1)))));
} else if (t_5 <= ((double) INFINITY)) {
tmp = t_5;
} else {
tmp = x * ((1.0 / t_6) + (z / (x * (y * pow(t_6, 2.0)))));
}
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 * ((y * ((x * y) + z)) + 27464.7644705)) + 230661.510616;
double t_2 = y * t_1;
double t_3 = (y * (y + a)) + b;
double t_4 = i + (y * (c + (y * t_3)));
double t_5 = (t + t_2) / t_4;
double t_6 = 1.0 + (a / y);
double tmp;
if (t_5 <= -2e-320) {
tmp = t * ((1.0 / t_4) + ((y / t) * (t_1 / t_4)));
} else if (t_5 <= 0.0) {
tmp = 1.0 / (c * ((1.0 / t_1) + ((i / (c * t_2)) + ((y / c) * (t_3 / t_1)))));
} else if (t_5 <= Double.POSITIVE_INFINITY) {
tmp = t_5;
} else {
tmp = x * ((1.0 / t_6) + (z / (x * (y * Math.pow(t_6, 2.0)))));
}
return tmp;
}
def code(x, y, z, t, a, b, c, i): t_1 = (y * ((y * ((x * y) + z)) + 27464.7644705)) + 230661.510616 t_2 = y * t_1 t_3 = (y * (y + a)) + b t_4 = i + (y * (c + (y * t_3))) t_5 = (t + t_2) / t_4 t_6 = 1.0 + (a / y) tmp = 0 if t_5 <= -2e-320: tmp = t * ((1.0 / t_4) + ((y / t) * (t_1 / t_4))) elif t_5 <= 0.0: tmp = 1.0 / (c * ((1.0 / t_1) + ((i / (c * t_2)) + ((y / c) * (t_3 / t_1))))) elif t_5 <= math.inf: tmp = t_5 else: tmp = x * ((1.0 / t_6) + (z / (x * (y * math.pow(t_6, 2.0))))) return tmp
function code(x, y, z, t, a, b, c, i) t_1 = Float64(Float64(y * Float64(Float64(y * Float64(Float64(x * y) + z)) + 27464.7644705)) + 230661.510616) t_2 = Float64(y * t_1) t_3 = Float64(Float64(y * Float64(y + a)) + b) t_4 = Float64(i + Float64(y * Float64(c + Float64(y * t_3)))) t_5 = Float64(Float64(t + t_2) / t_4) t_6 = Float64(1.0 + Float64(a / y)) tmp = 0.0 if (t_5 <= -2e-320) tmp = Float64(t * Float64(Float64(1.0 / t_4) + Float64(Float64(y / t) * Float64(t_1 / t_4)))); elseif (t_5 <= 0.0) tmp = Float64(1.0 / Float64(c * Float64(Float64(1.0 / t_1) + Float64(Float64(i / Float64(c * t_2)) + Float64(Float64(y / c) * Float64(t_3 / t_1)))))); elseif (t_5 <= Inf) tmp = t_5; else tmp = Float64(x * Float64(Float64(1.0 / t_6) + Float64(z / Float64(x * Float64(y * (t_6 ^ 2.0)))))); end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c, i) t_1 = (y * ((y * ((x * y) + z)) + 27464.7644705)) + 230661.510616; t_2 = y * t_1; t_3 = (y * (y + a)) + b; t_4 = i + (y * (c + (y * t_3))); t_5 = (t + t_2) / t_4; t_6 = 1.0 + (a / y); tmp = 0.0; if (t_5 <= -2e-320) tmp = t * ((1.0 / t_4) + ((y / t) * (t_1 / t_4))); elseif (t_5 <= 0.0) tmp = 1.0 / (c * ((1.0 / t_1) + ((i / (c * t_2)) + ((y / c) * (t_3 / t_1))))); elseif (t_5 <= Inf) tmp = t_5; else tmp = x * ((1.0 / t_6) + (z / (x * (y * (t_6 ^ 2.0))))); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_, i_] := Block[{t$95$1 = N[(N[(y * N[(N[(y * N[(N[(x * y), $MachinePrecision] + z), $MachinePrecision]), $MachinePrecision] + 27464.7644705), $MachinePrecision]), $MachinePrecision] + 230661.510616), $MachinePrecision]}, Block[{t$95$2 = N[(y * t$95$1), $MachinePrecision]}, Block[{t$95$3 = N[(N[(y * N[(y + a), $MachinePrecision]), $MachinePrecision] + b), $MachinePrecision]}, Block[{t$95$4 = N[(i + N[(y * N[(c + N[(y * t$95$3), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$5 = N[(N[(t + t$95$2), $MachinePrecision] / t$95$4), $MachinePrecision]}, Block[{t$95$6 = N[(1.0 + N[(a / y), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$5, -2e-320], N[(t * N[(N[(1.0 / t$95$4), $MachinePrecision] + N[(N[(y / t), $MachinePrecision] * N[(t$95$1 / t$95$4), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$5, 0.0], N[(1.0 / N[(c * N[(N[(1.0 / t$95$1), $MachinePrecision] + N[(N[(i / N[(c * t$95$2), $MachinePrecision]), $MachinePrecision] + N[(N[(y / c), $MachinePrecision] * N[(t$95$3 / t$95$1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$5, Infinity], t$95$5, N[(x * N[(N[(1.0 / t$95$6), $MachinePrecision] + N[(z / N[(x * N[(y * N[Power[t$95$6, 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := y \cdot \left(y \cdot \left(x \cdot y + z\right) + 27464.7644705\right) + 230661.510616\\
t_2 := y \cdot t\_1\\
t_3 := y \cdot \left(y + a\right) + b\\
t_4 := i + y \cdot \left(c + y \cdot t\_3\right)\\
t_5 := \frac{t + t\_2}{t\_4}\\
t_6 := 1 + \frac{a}{y}\\
\mathbf{if}\;t\_5 \leq -2 \cdot 10^{-320}:\\
\;\;\;\;t \cdot \left(\frac{1}{t\_4} + \frac{y}{t} \cdot \frac{t\_1}{t\_4}\right)\\
\mathbf{elif}\;t\_5 \leq 0:\\
\;\;\;\;\frac{1}{c \cdot \left(\frac{1}{t\_1} + \left(\frac{i}{c \cdot t\_2} + \frac{y}{c} \cdot \frac{t\_3}{t\_1}\right)\right)}\\
\mathbf{elif}\;t\_5 \leq \infty:\\
\;\;\;\;t\_5\\
\mathbf{else}:\\
\;\;\;\;x \cdot \left(\frac{1}{t\_6} + \frac{z}{x \cdot \left(y \cdot {t\_6}^{2}\right)}\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)) < -1.99998e-320Initial program 87.8%
fma-define87.8%
fma-define87.8%
fma-define87.8%
fma-define87.8%
fma-define87.8%
fma-define87.8%
fma-define87.8%
Simplified87.8%
Applied egg-rr87.3%
unpow-187.3%
+-commutative87.3%
Simplified87.3%
Taylor expanded in t around inf 80.7%
times-frac89.3%
Simplified89.3%
if -1.99998e-320 < (/.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)) < 0.0Initial program 48.0%
fma-define48.0%
fma-define48.0%
fma-define48.0%
fma-define48.0%
fma-define48.0%
fma-define48.0%
fma-define48.0%
Simplified48.0%
Applied egg-rr48.0%
unpow-148.0%
+-commutative48.0%
Simplified48.0%
Taylor expanded in t around 0 48.0%
associate-/r*48.0%
Simplified48.0%
Taylor expanded in c around inf 76.5%
times-frac94.8%
Simplified94.8%
if 0.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)) < +inf.0Initial program 97.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)) Initial program 0.0%
fma-define0.0%
fma-define0.0%
fma-define0.0%
fma-define0.0%
fma-define0.0%
fma-define0.0%
fma-define0.0%
Simplified0.0%
Applied egg-rr0.0%
unpow-10.0%
+-commutative0.0%
Simplified0.0%
Taylor expanded in y around inf 56.3%
Taylor expanded in x around inf 86.6%
Final simplification91.0%
(FPCore (x y z t a b c i)
:precision binary64
(let* ((t_1 (+ (* y (+ (* y (+ (* x y) z)) 27464.7644705)) 230661.510616))
(t_2 (* y t_1))
(t_3 (+ (* y (+ y a)) b))
(t_4 (+ i (* y (+ c (* y t_3)))))
(t_5 (/ (+ t t_2) t_4)))
(if (<= t_5 -2e-320)
(* t (+ (/ 1.0 t_4) (* (/ y t) (/ t_1 t_4))))
(if (<= t_5 0.0)
(/
1.0
(* c (+ (/ 1.0 t_1) (+ (/ i (* c t_2)) (* (/ y c) (/ t_3 t_1))))))
(if (<= t_5 INFINITY) t_5 (+ 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 * ((y * ((x * y) + z)) + 27464.7644705)) + 230661.510616;
double t_2 = y * t_1;
double t_3 = (y * (y + a)) + b;
double t_4 = i + (y * (c + (y * t_3)));
double t_5 = (t + t_2) / t_4;
double tmp;
if (t_5 <= -2e-320) {
tmp = t * ((1.0 / t_4) + ((y / t) * (t_1 / t_4)));
} else if (t_5 <= 0.0) {
tmp = 1.0 / (c * ((1.0 / t_1) + ((i / (c * t_2)) + ((y / c) * (t_3 / t_1)))));
} else if (t_5 <= ((double) INFINITY)) {
tmp = t_5;
} 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 * ((y * ((x * y) + z)) + 27464.7644705)) + 230661.510616;
double t_2 = y * t_1;
double t_3 = (y * (y + a)) + b;
double t_4 = i + (y * (c + (y * t_3)));
double t_5 = (t + t_2) / t_4;
double tmp;
if (t_5 <= -2e-320) {
tmp = t * ((1.0 / t_4) + ((y / t) * (t_1 / t_4)));
} else if (t_5 <= 0.0) {
tmp = 1.0 / (c * ((1.0 / t_1) + ((i / (c * t_2)) + ((y / c) * (t_3 / t_1)))));
} else if (t_5 <= Double.POSITIVE_INFINITY) {
tmp = t_5;
} else {
tmp = x + ((z / y) - (a * (x / y)));
}
return tmp;
}
def code(x, y, z, t, a, b, c, i): t_1 = (y * ((y * ((x * y) + z)) + 27464.7644705)) + 230661.510616 t_2 = y * t_1 t_3 = (y * (y + a)) + b t_4 = i + (y * (c + (y * t_3))) t_5 = (t + t_2) / t_4 tmp = 0 if t_5 <= -2e-320: tmp = t * ((1.0 / t_4) + ((y / t) * (t_1 / t_4))) elif t_5 <= 0.0: tmp = 1.0 / (c * ((1.0 / t_1) + ((i / (c * t_2)) + ((y / c) * (t_3 / t_1))))) elif t_5 <= math.inf: tmp = t_5 else: tmp = x + ((z / y) - (a * (x / y))) return tmp
function code(x, y, z, t, a, b, c, i) t_1 = Float64(Float64(y * Float64(Float64(y * Float64(Float64(x * y) + z)) + 27464.7644705)) + 230661.510616) t_2 = Float64(y * t_1) t_3 = Float64(Float64(y * Float64(y + a)) + b) t_4 = Float64(i + Float64(y * Float64(c + Float64(y * t_3)))) t_5 = Float64(Float64(t + t_2) / t_4) tmp = 0.0 if (t_5 <= -2e-320) tmp = Float64(t * Float64(Float64(1.0 / t_4) + Float64(Float64(y / t) * Float64(t_1 / t_4)))); elseif (t_5 <= 0.0) tmp = Float64(1.0 / Float64(c * Float64(Float64(1.0 / t_1) + Float64(Float64(i / Float64(c * t_2)) + Float64(Float64(y / c) * Float64(t_3 / t_1)))))); elseif (t_5 <= Inf) tmp = t_5; 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 * ((y * ((x * y) + z)) + 27464.7644705)) + 230661.510616; t_2 = y * t_1; t_3 = (y * (y + a)) + b; t_4 = i + (y * (c + (y * t_3))); t_5 = (t + t_2) / t_4; tmp = 0.0; if (t_5 <= -2e-320) tmp = t * ((1.0 / t_4) + ((y / t) * (t_1 / t_4))); elseif (t_5 <= 0.0) tmp = 1.0 / (c * ((1.0 / t_1) + ((i / (c * t_2)) + ((y / c) * (t_3 / t_1))))); elseif (t_5 <= Inf) tmp = t_5; 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[(N[(y * N[(N[(y * N[(N[(x * y), $MachinePrecision] + z), $MachinePrecision]), $MachinePrecision] + 27464.7644705), $MachinePrecision]), $MachinePrecision] + 230661.510616), $MachinePrecision]}, Block[{t$95$2 = N[(y * t$95$1), $MachinePrecision]}, Block[{t$95$3 = N[(N[(y * N[(y + a), $MachinePrecision]), $MachinePrecision] + b), $MachinePrecision]}, Block[{t$95$4 = N[(i + N[(y * N[(c + N[(y * t$95$3), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$5 = N[(N[(t + t$95$2), $MachinePrecision] / t$95$4), $MachinePrecision]}, If[LessEqual[t$95$5, -2e-320], N[(t * N[(N[(1.0 / t$95$4), $MachinePrecision] + N[(N[(y / t), $MachinePrecision] * N[(t$95$1 / t$95$4), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$5, 0.0], N[(1.0 / N[(c * N[(N[(1.0 / t$95$1), $MachinePrecision] + N[(N[(i / N[(c * t$95$2), $MachinePrecision]), $MachinePrecision] + N[(N[(y / c), $MachinePrecision] * N[(t$95$3 / t$95$1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$5, Infinity], t$95$5, 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(y \cdot \left(x \cdot y + z\right) + 27464.7644705\right) + 230661.510616\\
t_2 := y \cdot t\_1\\
t_3 := y \cdot \left(y + a\right) + b\\
t_4 := i + y \cdot \left(c + y \cdot t\_3\right)\\
t_5 := \frac{t + t\_2}{t\_4}\\
\mathbf{if}\;t\_5 \leq -2 \cdot 10^{-320}:\\
\;\;\;\;t \cdot \left(\frac{1}{t\_4} + \frac{y}{t} \cdot \frac{t\_1}{t\_4}\right)\\
\mathbf{elif}\;t\_5 \leq 0:\\
\;\;\;\;\frac{1}{c \cdot \left(\frac{1}{t\_1} + \left(\frac{i}{c \cdot t\_2} + \frac{y}{c} \cdot \frac{t\_3}{t\_1}\right)\right)}\\
\mathbf{elif}\;t\_5 \leq \infty:\\
\;\;\;\;t\_5\\
\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)) < -1.99998e-320Initial program 87.8%
fma-define87.8%
fma-define87.8%
fma-define87.8%
fma-define87.8%
fma-define87.8%
fma-define87.8%
fma-define87.8%
Simplified87.8%
Applied egg-rr87.3%
unpow-187.3%
+-commutative87.3%
Simplified87.3%
Taylor expanded in t around inf 80.7%
times-frac89.3%
Simplified89.3%
if -1.99998e-320 < (/.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)) < 0.0Initial program 48.0%
fma-define48.0%
fma-define48.0%
fma-define48.0%
fma-define48.0%
fma-define48.0%
fma-define48.0%
fma-define48.0%
Simplified48.0%
Applied egg-rr48.0%
unpow-148.0%
+-commutative48.0%
Simplified48.0%
Taylor expanded in t around 0 48.0%
associate-/r*48.0%
Simplified48.0%
Taylor expanded in c around inf 76.5%
times-frac94.8%
Simplified94.8%
if 0.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)) < +inf.0Initial program 97.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)) Initial program 0.0%
fma-define0.0%
fma-define0.0%
fma-define0.0%
fma-define0.0%
fma-define0.0%
fma-define0.0%
fma-define0.0%
Simplified0.0%
Taylor expanded in y around inf 73.4%
associate--l+73.4%
associate-/l*76.0%
Simplified76.0%
Final simplification86.7%
(FPCore (x y z t a b c i)
:precision binary64
(let* ((t_1 (+ (* y (+ (* y (+ (* x y) z)) 27464.7644705)) 230661.510616)))
(if (<= y -3.3e+40)
(/ x (+ 1.0 (/ a y)))
(if (<= y 2.4e-10)
(/ (+ t (* y t_1)) (+ i (* y (+ c (* y (+ b (* y a)))))))
(if (<= y 7.6e+61)
(/ 1.0 (/ (+ c (+ (* y b) (* y (* y (+ y a))))) 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 = (y * ((y * ((x * y) + z)) + 27464.7644705)) + 230661.510616;
double tmp;
if (y <= -3.3e+40) {
tmp = x / (1.0 + (a / y));
} else if (y <= 2.4e-10) {
tmp = (t + (y * t_1)) / (i + (y * (c + (y * (b + (y * a))))));
} else if (y <= 7.6e+61) {
tmp = 1.0 / ((c + ((y * b) + (y * (y * (y + a))))) / t_1);
} else {
tmp = x + ((z / y) - (a * (x / y)));
}
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 = (y * ((y * ((x * y) + z)) + 27464.7644705d0)) + 230661.510616d0
if (y <= (-3.3d+40)) then
tmp = x / (1.0d0 + (a / y))
else if (y <= 2.4d-10) then
tmp = (t + (y * t_1)) / (i + (y * (c + (y * (b + (y * a))))))
else if (y <= 7.6d+61) then
tmp = 1.0d0 / ((c + ((y * b) + (y * (y * (y + a))))) / t_1)
else
tmp = x + ((z / y) - (a * (x / y)))
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 * ((y * ((x * y) + z)) + 27464.7644705)) + 230661.510616;
double tmp;
if (y <= -3.3e+40) {
tmp = x / (1.0 + (a / y));
} else if (y <= 2.4e-10) {
tmp = (t + (y * t_1)) / (i + (y * (c + (y * (b + (y * a))))));
} else if (y <= 7.6e+61) {
tmp = 1.0 / ((c + ((y * b) + (y * (y * (y + a))))) / t_1);
} else {
tmp = x + ((z / y) - (a * (x / y)));
}
return tmp;
}
def code(x, y, z, t, a, b, c, i): t_1 = (y * ((y * ((x * y) + z)) + 27464.7644705)) + 230661.510616 tmp = 0 if y <= -3.3e+40: tmp = x / (1.0 + (a / y)) elif y <= 2.4e-10: tmp = (t + (y * t_1)) / (i + (y * (c + (y * (b + (y * a)))))) elif y <= 7.6e+61: tmp = 1.0 / ((c + ((y * b) + (y * (y * (y + a))))) / 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(Float64(y * Float64(Float64(y * Float64(Float64(x * y) + z)) + 27464.7644705)) + 230661.510616) tmp = 0.0 if (y <= -3.3e+40) tmp = Float64(x / Float64(1.0 + Float64(a / y))); elseif (y <= 2.4e-10) tmp = Float64(Float64(t + Float64(y * t_1)) / Float64(i + Float64(y * Float64(c + Float64(y * Float64(b + Float64(y * a))))))); elseif (y <= 7.6e+61) tmp = Float64(1.0 / Float64(Float64(c + Float64(Float64(y * b) + Float64(y * Float64(y * Float64(y + a))))) / 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 = (y * ((y * ((x * y) + z)) + 27464.7644705)) + 230661.510616; tmp = 0.0; if (y <= -3.3e+40) tmp = x / (1.0 + (a / y)); elseif (y <= 2.4e-10) tmp = (t + (y * t_1)) / (i + (y * (c + (y * (b + (y * a)))))); elseif (y <= 7.6e+61) tmp = 1.0 / ((c + ((y * b) + (y * (y * (y + a))))) / 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[(N[(y * N[(N[(y * N[(N[(x * y), $MachinePrecision] + z), $MachinePrecision]), $MachinePrecision] + 27464.7644705), $MachinePrecision]), $MachinePrecision] + 230661.510616), $MachinePrecision]}, If[LessEqual[y, -3.3e+40], N[(x / N[(1.0 + N[(a / y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 2.4e-10], N[(N[(t + N[(y * t$95$1), $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.6e+61], N[(1.0 / N[(N[(c + N[(N[(y * b), $MachinePrecision] + N[(y * N[(y * N[(y + a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / 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 := y \cdot \left(y \cdot \left(x \cdot y + z\right) + 27464.7644705\right) + 230661.510616\\
\mathbf{if}\;y \leq -3.3 \cdot 10^{+40}:\\
\;\;\;\;\frac{x}{1 + \frac{a}{y}}\\
\mathbf{elif}\;y \leq 2.4 \cdot 10^{-10}:\\
\;\;\;\;\frac{t + y \cdot t\_1}{i + y \cdot \left(c + y \cdot \left(b + y \cdot a\right)\right)}\\
\mathbf{elif}\;y \leq 7.6 \cdot 10^{+61}:\\
\;\;\;\;\frac{1}{\frac{c + \left(y \cdot b + y \cdot \left(y \cdot \left(y + a\right)\right)\right)}{t\_1}}\\
\mathbf{else}:\\
\;\;\;\;x + \left(\frac{z}{y} - a \cdot \frac{x}{y}\right)\\
\end{array}
\end{array}
if y < -3.2999999999999998e40Initial program 4.2%
fma-define4.2%
fma-define4.2%
fma-define4.2%
fma-define4.2%
fma-define4.2%
fma-define4.2%
fma-define4.2%
Simplified4.2%
Applied egg-rr4.2%
unpow-14.2%
+-commutative4.2%
Simplified4.2%
Taylor expanded in y around inf 58.7%
Taylor expanded in x around inf 71.6%
if -3.2999999999999998e40 < y < 2.4e-10Initial program 98.9%
Taylor expanded in y around 0 98.8%
if 2.4e-10 < y < 7.5999999999999999e61Initial program 51.7%
fma-define51.7%
fma-define51.7%
fma-define51.7%
fma-define51.7%
fma-define51.7%
fma-define51.7%
fma-define51.7%
Simplified51.7%
Applied egg-rr51.8%
unpow-151.8%
+-commutative51.8%
Simplified51.8%
Taylor expanded in t around 0 45.1%
associate-/r*51.9%
Simplified51.9%
Taylor expanded in i around 0 65.7%
distribute-lft-in65.8%
+-commutative65.8%
Applied egg-rr65.8%
if 7.5999999999999999e61 < y Initial program 0.4%
fma-define0.4%
fma-define0.4%
fma-define0.4%
fma-define0.4%
fma-define0.4%
fma-define0.4%
fma-define0.4%
Simplified0.4%
Taylor expanded in y around inf 72.3%
associate--l+72.3%
associate-/l*73.8%
Simplified73.8%
Final simplification85.1%
(FPCore (x y z t a b c i)
:precision binary64
(let* ((t_1 (* y (+ y a))))
(if (<= y -3.5e+37)
(/ x (+ 1.0 (/ a y)))
(if (<= y 2.4e-10)
(/
(+ t (* y (+ 230661.510616 (* y (+ 27464.7644705 (* y z))))))
(+ i (* y (+ c (* y (+ t_1 b))))))
(if (<= y 8.2e+63)
(/
1.0
(/
(+ c (+ (* y b) (* y t_1)))
(+ (* y (+ (* y (+ (* x y) z)) 27464.7644705)) 230661.510616)))
(+ 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 * (y + a);
double tmp;
if (y <= -3.5e+37) {
tmp = x / (1.0 + (a / y));
} else if (y <= 2.4e-10) {
tmp = (t + (y * (230661.510616 + (y * (27464.7644705 + (y * z)))))) / (i + (y * (c + (y * (t_1 + b)))));
} else if (y <= 8.2e+63) {
tmp = 1.0 / ((c + ((y * b) + (y * t_1))) / ((y * ((y * ((x * y) + z)) + 27464.7644705)) + 230661.510616));
} else {
tmp = x + ((z / y) - (a * (x / y)));
}
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 = y * (y + a)
if (y <= (-3.5d+37)) then
tmp = x / (1.0d0 + (a / y))
else if (y <= 2.4d-10) then
tmp = (t + (y * (230661.510616d0 + (y * (27464.7644705d0 + (y * z)))))) / (i + (y * (c + (y * (t_1 + b)))))
else if (y <= 8.2d+63) then
tmp = 1.0d0 / ((c + ((y * b) + (y * t_1))) / ((y * ((y * ((x * y) + z)) + 27464.7644705d0)) + 230661.510616d0))
else
tmp = x + ((z / y) - (a * (x / y)))
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 * (y + a);
double tmp;
if (y <= -3.5e+37) {
tmp = x / (1.0 + (a / y));
} else if (y <= 2.4e-10) {
tmp = (t + (y * (230661.510616 + (y * (27464.7644705 + (y * z)))))) / (i + (y * (c + (y * (t_1 + b)))));
} else if (y <= 8.2e+63) {
tmp = 1.0 / ((c + ((y * b) + (y * t_1))) / ((y * ((y * ((x * y) + z)) + 27464.7644705)) + 230661.510616));
} else {
tmp = x + ((z / y) - (a * (x / y)));
}
return tmp;
}
def code(x, y, z, t, a, b, c, i): t_1 = y * (y + a) tmp = 0 if y <= -3.5e+37: tmp = x / (1.0 + (a / y)) elif y <= 2.4e-10: tmp = (t + (y * (230661.510616 + (y * (27464.7644705 + (y * z)))))) / (i + (y * (c + (y * (t_1 + b))))) elif y <= 8.2e+63: tmp = 1.0 / ((c + ((y * b) + (y * t_1))) / ((y * ((y * ((x * y) + z)) + 27464.7644705)) + 230661.510616)) 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(y + a)) tmp = 0.0 if (y <= -3.5e+37) tmp = Float64(x / Float64(1.0 + Float64(a / y))); elseif (y <= 2.4e-10) 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(t_1 + b)))))); elseif (y <= 8.2e+63) tmp = Float64(1.0 / Float64(Float64(c + Float64(Float64(y * b) + Float64(y * t_1))) / Float64(Float64(y * Float64(Float64(y * Float64(Float64(x * y) + z)) + 27464.7644705)) + 230661.510616))); 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 * (y + a); tmp = 0.0; if (y <= -3.5e+37) tmp = x / (1.0 + (a / y)); elseif (y <= 2.4e-10) tmp = (t + (y * (230661.510616 + (y * (27464.7644705 + (y * z)))))) / (i + (y * (c + (y * (t_1 + b))))); elseif (y <= 8.2e+63) tmp = 1.0 / ((c + ((y * b) + (y * t_1))) / ((y * ((y * ((x * y) + z)) + 27464.7644705)) + 230661.510616)); 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[(y + a), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -3.5e+37], N[(x / N[(1.0 + N[(a / y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 2.4e-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 * N[(c + N[(y * N[(t$95$1 + b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 8.2e+63], N[(1.0 / N[(N[(c + N[(N[(y * b), $MachinePrecision] + N[(y * t$95$1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(N[(y * N[(N[(y * N[(N[(x * y), $MachinePrecision] + z), $MachinePrecision]), $MachinePrecision] + 27464.7644705), $MachinePrecision]), $MachinePrecision] + 230661.510616), $MachinePrecision]), $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(y + a\right)\\
\mathbf{if}\;y \leq -3.5 \cdot 10^{+37}:\\
\;\;\;\;\frac{x}{1 + \frac{a}{y}}\\
\mathbf{elif}\;y \leq 2.4 \cdot 10^{-10}:\\
\;\;\;\;\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(t\_1 + b\right)\right)}\\
\mathbf{elif}\;y \leq 8.2 \cdot 10^{+63}:\\
\;\;\;\;\frac{1}{\frac{c + \left(y \cdot b + y \cdot t\_1\right)}{y \cdot \left(y \cdot \left(x \cdot y + z\right) + 27464.7644705\right) + 230661.510616}}\\
\mathbf{else}:\\
\;\;\;\;x + \left(\frac{z}{y} - a \cdot \frac{x}{y}\right)\\
\end{array}
\end{array}
if y < -3.5e37Initial program 4.2%
fma-define4.2%
fma-define4.2%
fma-define4.2%
fma-define4.2%
fma-define4.2%
fma-define4.2%
fma-define4.2%
Simplified4.2%
Applied egg-rr4.2%
unpow-14.2%
+-commutative4.2%
Simplified4.2%
Taylor expanded in y around inf 58.7%
Taylor expanded in x around inf 71.6%
if -3.5e37 < y < 2.4e-10Initial program 98.9%
Taylor expanded in x around 0 95.8%
if 2.4e-10 < y < 8.19999999999999985e63Initial program 51.7%
fma-define51.7%
fma-define51.7%
fma-define51.7%
fma-define51.7%
fma-define51.7%
fma-define51.7%
fma-define51.7%
Simplified51.7%
Applied egg-rr51.8%
unpow-151.8%
+-commutative51.8%
Simplified51.8%
Taylor expanded in t around 0 45.1%
associate-/r*51.9%
Simplified51.9%
Taylor expanded in i around 0 65.7%
distribute-lft-in65.8%
+-commutative65.8%
Applied egg-rr65.8%
if 8.19999999999999985e63 < y Initial program 0.4%
fma-define0.4%
fma-define0.4%
fma-define0.4%
fma-define0.4%
fma-define0.4%
fma-define0.4%
fma-define0.4%
Simplified0.4%
Taylor expanded in y around inf 72.3%
associate--l+72.3%
associate-/l*73.8%
Simplified73.8%
Final simplification83.6%
(FPCore (x y z t a b c i)
:precision binary64
(let* ((t_1 (+ (* y (+ (* y (+ (* x y) z)) 27464.7644705)) 230661.510616)))
(if (<= y -2.15e+38)
(/ x (+ 1.0 (/ a y)))
(if (<= y 2.4e-10)
(/ (+ t (* y t_1)) (+ i (* y (+ c (* y b)))))
(if (<= y 6.3e+63)
(/ 1.0 (/ (+ c (+ (* y b) (* y (* y (+ y a))))) 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 = (y * ((y * ((x * y) + z)) + 27464.7644705)) + 230661.510616;
double tmp;
if (y <= -2.15e+38) {
tmp = x / (1.0 + (a / y));
} else if (y <= 2.4e-10) {
tmp = (t + (y * t_1)) / (i + (y * (c + (y * b))));
} else if (y <= 6.3e+63) {
tmp = 1.0 / ((c + ((y * b) + (y * (y * (y + a))))) / t_1);
} else {
tmp = x + ((z / y) - (a * (x / y)));
}
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 = (y * ((y * ((x * y) + z)) + 27464.7644705d0)) + 230661.510616d0
if (y <= (-2.15d+38)) then
tmp = x / (1.0d0 + (a / y))
else if (y <= 2.4d-10) then
tmp = (t + (y * t_1)) / (i + (y * (c + (y * b))))
else if (y <= 6.3d+63) then
tmp = 1.0d0 / ((c + ((y * b) + (y * (y * (y + a))))) / t_1)
else
tmp = x + ((z / y) - (a * (x / y)))
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 * ((y * ((x * y) + z)) + 27464.7644705)) + 230661.510616;
double tmp;
if (y <= -2.15e+38) {
tmp = x / (1.0 + (a / y));
} else if (y <= 2.4e-10) {
tmp = (t + (y * t_1)) / (i + (y * (c + (y * b))));
} else if (y <= 6.3e+63) {
tmp = 1.0 / ((c + ((y * b) + (y * (y * (y + a))))) / t_1);
} else {
tmp = x + ((z / y) - (a * (x / y)));
}
return tmp;
}
def code(x, y, z, t, a, b, c, i): t_1 = (y * ((y * ((x * y) + z)) + 27464.7644705)) + 230661.510616 tmp = 0 if y <= -2.15e+38: tmp = x / (1.0 + (a / y)) elif y <= 2.4e-10: tmp = (t + (y * t_1)) / (i + (y * (c + (y * b)))) elif y <= 6.3e+63: tmp = 1.0 / ((c + ((y * b) + (y * (y * (y + a))))) / 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(Float64(y * Float64(Float64(y * Float64(Float64(x * y) + z)) + 27464.7644705)) + 230661.510616) tmp = 0.0 if (y <= -2.15e+38) tmp = Float64(x / Float64(1.0 + Float64(a / y))); elseif (y <= 2.4e-10) tmp = Float64(Float64(t + Float64(y * t_1)) / Float64(i + Float64(y * Float64(c + Float64(y * b))))); elseif (y <= 6.3e+63) tmp = Float64(1.0 / Float64(Float64(c + Float64(Float64(y * b) + Float64(y * Float64(y * Float64(y + a))))) / 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 = (y * ((y * ((x * y) + z)) + 27464.7644705)) + 230661.510616; tmp = 0.0; if (y <= -2.15e+38) tmp = x / (1.0 + (a / y)); elseif (y <= 2.4e-10) tmp = (t + (y * t_1)) / (i + (y * (c + (y * b)))); elseif (y <= 6.3e+63) tmp = 1.0 / ((c + ((y * b) + (y * (y * (y + a))))) / 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[(N[(y * N[(N[(y * N[(N[(x * y), $MachinePrecision] + z), $MachinePrecision]), $MachinePrecision] + 27464.7644705), $MachinePrecision]), $MachinePrecision] + 230661.510616), $MachinePrecision]}, If[LessEqual[y, -2.15e+38], N[(x / N[(1.0 + N[(a / y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 2.4e-10], N[(N[(t + N[(y * t$95$1), $MachinePrecision]), $MachinePrecision] / N[(i + N[(y * N[(c + N[(y * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 6.3e+63], N[(1.0 / N[(N[(c + N[(N[(y * b), $MachinePrecision] + N[(y * N[(y * N[(y + a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / 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 := y \cdot \left(y \cdot \left(x \cdot y + z\right) + 27464.7644705\right) + 230661.510616\\
\mathbf{if}\;y \leq -2.15 \cdot 10^{+38}:\\
\;\;\;\;\frac{x}{1 + \frac{a}{y}}\\
\mathbf{elif}\;y \leq 2.4 \cdot 10^{-10}:\\
\;\;\;\;\frac{t + y \cdot t\_1}{i + y \cdot \left(c + y \cdot b\right)}\\
\mathbf{elif}\;y \leq 6.3 \cdot 10^{+63}:\\
\;\;\;\;\frac{1}{\frac{c + \left(y \cdot b + y \cdot \left(y \cdot \left(y + a\right)\right)\right)}{t\_1}}\\
\mathbf{else}:\\
\;\;\;\;x + \left(\frac{z}{y} - a \cdot \frac{x}{y}\right)\\
\end{array}
\end{array}
if y < -2.1499999999999998e38Initial program 4.2%
fma-define4.2%
fma-define4.2%
fma-define4.2%
fma-define4.2%
fma-define4.2%
fma-define4.2%
fma-define4.2%
Simplified4.2%
Applied egg-rr4.2%
unpow-14.2%
+-commutative4.2%
Simplified4.2%
Taylor expanded in y around inf 58.7%
Taylor expanded in x around inf 71.6%
if -2.1499999999999998e38 < y < 2.4e-10Initial program 98.9%
Taylor expanded in y around 0 94.9%
*-commutative94.9%
Simplified94.9%
if 2.4e-10 < y < 6.2999999999999998e63Initial program 51.7%
fma-define51.7%
fma-define51.7%
fma-define51.7%
fma-define51.7%
fma-define51.7%
fma-define51.7%
fma-define51.7%
Simplified51.7%
Applied egg-rr51.8%
unpow-151.8%
+-commutative51.8%
Simplified51.8%
Taylor expanded in t around 0 45.1%
associate-/r*51.9%
Simplified51.9%
Taylor expanded in i around 0 65.7%
distribute-lft-in65.8%
+-commutative65.8%
Applied egg-rr65.8%
if 6.2999999999999998e63 < y Initial program 0.4%
fma-define0.4%
fma-define0.4%
fma-define0.4%
fma-define0.4%
fma-define0.4%
fma-define0.4%
fma-define0.4%
Simplified0.4%
Taylor expanded in y around inf 72.3%
associate--l+72.3%
associate-/l*73.8%
Simplified73.8%
Final simplification83.2%
(FPCore (x y z t a b c i)
:precision binary64
(if (<= y -4.4e+40)
(/ x (+ 1.0 (/ a y)))
(if (<= y 1.32e+68)
(/
(+ t (* y (+ (* y (+ (* y (+ (* x y) z)) 27464.7644705)) 230661.510616)))
(+ i (* y (+ c (* y (+ (* y (+ y a)) b))))))
(+ 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 tmp;
if (y <= -4.4e+40) {
tmp = x / (1.0 + (a / y));
} else if (y <= 1.32e+68) {
tmp = (t + (y * ((y * ((y * ((x * y) + z)) + 27464.7644705)) + 230661.510616))) / (i + (y * (c + (y * ((y * (y + a)) + b)))));
} else {
tmp = x + ((z / y) - (a * (x / y)));
}
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 <= (-4.4d+40)) then
tmp = x / (1.0d0 + (a / y))
else if (y <= 1.32d+68) then
tmp = (t + (y * ((y * ((y * ((x * y) + z)) + 27464.7644705d0)) + 230661.510616d0))) / (i + (y * (c + (y * ((y * (y + a)) + b)))))
else
tmp = x + ((z / y) - (a * (x / y)))
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 <= -4.4e+40) {
tmp = x / (1.0 + (a / y));
} else if (y <= 1.32e+68) {
tmp = (t + (y * ((y * ((y * ((x * y) + z)) + 27464.7644705)) + 230661.510616))) / (i + (y * (c + (y * ((y * (y + a)) + b)))));
} else {
tmp = x + ((z / y) - (a * (x / y)));
}
return tmp;
}
def code(x, y, z, t, a, b, c, i): tmp = 0 if y <= -4.4e+40: tmp = x / (1.0 + (a / y)) elif y <= 1.32e+68: tmp = (t + (y * ((y * ((y * ((x * y) + z)) + 27464.7644705)) + 230661.510616))) / (i + (y * (c + (y * ((y * (y + a)) + b))))) else: tmp = x + ((z / y) - (a * (x / y))) return tmp
function code(x, y, z, t, a, b, c, i) tmp = 0.0 if (y <= -4.4e+40) tmp = Float64(x / Float64(1.0 + Float64(a / y))); elseif (y <= 1.32e+68) tmp = Float64(Float64(t + Float64(y * Float64(Float64(y * Float64(Float64(y * Float64(Float64(x * y) + z)) + 27464.7644705)) + 230661.510616))) / Float64(i + Float64(y * Float64(c + Float64(y * Float64(Float64(y * Float64(y + a)) + b)))))); 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) tmp = 0.0; if (y <= -4.4e+40) tmp = x / (1.0 + (a / y)); elseif (y <= 1.32e+68) tmp = (t + (y * ((y * ((y * ((x * y) + z)) + 27464.7644705)) + 230661.510616))) / (i + (y * (c + (y * ((y * (y + a)) + b))))); else tmp = x + ((z / y) - (a * (x / y))); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_, i_] := If[LessEqual[y, -4.4e+40], N[(x / N[(1.0 + N[(a / y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 1.32e+68], N[(N[(t + N[(y * N[(N[(y * N[(N[(y * N[(N[(x * y), $MachinePrecision] + z), $MachinePrecision]), $MachinePrecision] + 27464.7644705), $MachinePrecision]), $MachinePrecision] + 230661.510616), $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], N[(x + N[(N[(z / y), $MachinePrecision] - N[(a * N[(x / y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -4.4 \cdot 10^{+40}:\\
\;\;\;\;\frac{x}{1 + \frac{a}{y}}\\
\mathbf{elif}\;y \leq 1.32 \cdot 10^{+68}:\\
\;\;\;\;\frac{t + y \cdot \left(y \cdot \left(y \cdot \left(x \cdot y + z\right) + 27464.7644705\right) + 230661.510616\right)}{i + y \cdot \left(c + y \cdot \left(y \cdot \left(y + a\right) + b\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;x + \left(\frac{z}{y} - a \cdot \frac{x}{y}\right)\\
\end{array}
\end{array}
if y < -4.3999999999999998e40Initial program 4.2%
fma-define4.2%
fma-define4.2%
fma-define4.2%
fma-define4.2%
fma-define4.2%
fma-define4.2%
fma-define4.2%
Simplified4.2%
Applied egg-rr4.2%
unpow-14.2%
+-commutative4.2%
Simplified4.2%
Taylor expanded in y around inf 58.7%
Taylor expanded in x around inf 71.6%
if -4.3999999999999998e40 < y < 1.3200000000000001e68Initial program 93.5%
if 1.3200000000000001e68 < y Initial program 0.3%
fma-define0.3%
fma-define0.3%
fma-define0.3%
fma-define0.3%
fma-define0.3%
fma-define0.3%
fma-define0.3%
Simplified0.3%
Taylor expanded in y around inf 73.5%
associate--l+73.5%
associate-/l*75.1%
Simplified75.1%
Final simplification84.3%
(FPCore (x y z t a b c i)
:precision binary64
(let* ((t_1 (+ (* y (+ (* y (+ (* x y) z)) 27464.7644705)) 230661.510616)))
(if (<= y -4.1e+29)
(/ x (+ 1.0 (/ a y)))
(if (<= y 2.4e-10)
(/ (+ t (* y t_1)) (+ i (* y (+ c (* y b)))))
(if (<= y 1.75e+62)
(/ t_1 (+ c (* y (+ (* y (+ y a)) b))))
(+ 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 * ((y * ((x * y) + z)) + 27464.7644705)) + 230661.510616;
double tmp;
if (y <= -4.1e+29) {
tmp = x / (1.0 + (a / y));
} else if (y <= 2.4e-10) {
tmp = (t + (y * t_1)) / (i + (y * (c + (y * b))));
} else if (y <= 1.75e+62) {
tmp = t_1 / (c + (y * ((y * (y + a)) + b)));
} else {
tmp = x + ((z / y) - (a * (x / y)));
}
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 = (y * ((y * ((x * y) + z)) + 27464.7644705d0)) + 230661.510616d0
if (y <= (-4.1d+29)) then
tmp = x / (1.0d0 + (a / y))
else if (y <= 2.4d-10) then
tmp = (t + (y * t_1)) / (i + (y * (c + (y * b))))
else if (y <= 1.75d+62) then
tmp = t_1 / (c + (y * ((y * (y + a)) + b)))
else
tmp = x + ((z / y) - (a * (x / y)))
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 * ((y * ((x * y) + z)) + 27464.7644705)) + 230661.510616;
double tmp;
if (y <= -4.1e+29) {
tmp = x / (1.0 + (a / y));
} else if (y <= 2.4e-10) {
tmp = (t + (y * t_1)) / (i + (y * (c + (y * b))));
} else if (y <= 1.75e+62) {
tmp = t_1 / (c + (y * ((y * (y + a)) + b)));
} else {
tmp = x + ((z / y) - (a * (x / y)));
}
return tmp;
}
def code(x, y, z, t, a, b, c, i): t_1 = (y * ((y * ((x * y) + z)) + 27464.7644705)) + 230661.510616 tmp = 0 if y <= -4.1e+29: tmp = x / (1.0 + (a / y)) elif y <= 2.4e-10: tmp = (t + (y * t_1)) / (i + (y * (c + (y * b)))) elif y <= 1.75e+62: tmp = t_1 / (c + (y * ((y * (y + a)) + b))) else: tmp = x + ((z / y) - (a * (x / y))) return tmp
function code(x, y, z, t, a, b, c, i) t_1 = Float64(Float64(y * Float64(Float64(y * Float64(Float64(x * y) + z)) + 27464.7644705)) + 230661.510616) tmp = 0.0 if (y <= -4.1e+29) tmp = Float64(x / Float64(1.0 + Float64(a / y))); elseif (y <= 2.4e-10) tmp = Float64(Float64(t + Float64(y * t_1)) / Float64(i + Float64(y * Float64(c + Float64(y * b))))); elseif (y <= 1.75e+62) tmp = Float64(t_1 / Float64(c + Float64(y * Float64(Float64(y * Float64(y + a)) + b)))); 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 * ((y * ((x * y) + z)) + 27464.7644705)) + 230661.510616; tmp = 0.0; if (y <= -4.1e+29) tmp = x / (1.0 + (a / y)); elseif (y <= 2.4e-10) tmp = (t + (y * t_1)) / (i + (y * (c + (y * b)))); elseif (y <= 1.75e+62) tmp = t_1 / (c + (y * ((y * (y + a)) + b))); 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[(N[(y * N[(N[(y * N[(N[(x * y), $MachinePrecision] + z), $MachinePrecision]), $MachinePrecision] + 27464.7644705), $MachinePrecision]), $MachinePrecision] + 230661.510616), $MachinePrecision]}, If[LessEqual[y, -4.1e+29], N[(x / N[(1.0 + N[(a / y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 2.4e-10], N[(N[(t + N[(y * t$95$1), $MachinePrecision]), $MachinePrecision] / N[(i + N[(y * N[(c + N[(y * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 1.75e+62], N[(t$95$1 / N[(c + N[(y * N[(N[(y * N[(y + a), $MachinePrecision]), $MachinePrecision] + b), $MachinePrecision]), $MachinePrecision]), $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(y \cdot \left(x \cdot y + z\right) + 27464.7644705\right) + 230661.510616\\
\mathbf{if}\;y \leq -4.1 \cdot 10^{+29}:\\
\;\;\;\;\frac{x}{1 + \frac{a}{y}}\\
\mathbf{elif}\;y \leq 2.4 \cdot 10^{-10}:\\
\;\;\;\;\frac{t + y \cdot t\_1}{i + y \cdot \left(c + y \cdot b\right)}\\
\mathbf{elif}\;y \leq 1.75 \cdot 10^{+62}:\\
\;\;\;\;\frac{t\_1}{c + y \cdot \left(y \cdot \left(y + a\right) + b\right)}\\
\mathbf{else}:\\
\;\;\;\;x + \left(\frac{z}{y} - a \cdot \frac{x}{y}\right)\\
\end{array}
\end{array}
if y < -4.1000000000000003e29Initial program 4.2%
fma-define4.2%
fma-define4.2%
fma-define4.2%
fma-define4.2%
fma-define4.2%
fma-define4.2%
fma-define4.2%
Simplified4.2%
Applied egg-rr4.2%
unpow-14.2%
+-commutative4.2%
Simplified4.2%
Taylor expanded in y around inf 58.7%
Taylor expanded in x around inf 71.6%
if -4.1000000000000003e29 < y < 2.4e-10Initial program 98.9%
Taylor expanded in y around 0 94.9%
*-commutative94.9%
Simplified94.9%
if 2.4e-10 < y < 1.74999999999999992e62Initial program 51.7%
fma-define51.7%
fma-define51.7%
fma-define51.7%
fma-define51.7%
fma-define51.7%
fma-define51.7%
fma-define51.7%
Simplified51.7%
Applied egg-rr51.8%
unpow-151.8%
+-commutative51.8%
Simplified51.8%
Taylor expanded in t around 0 45.1%
associate-/r*51.9%
Simplified51.9%
Taylor expanded in i around 0 65.7%
if 1.74999999999999992e62 < y Initial program 0.4%
fma-define0.4%
fma-define0.4%
fma-define0.4%
fma-define0.4%
fma-define0.4%
fma-define0.4%
fma-define0.4%
Simplified0.4%
Taylor expanded in y around inf 72.3%
associate--l+72.3%
associate-/l*73.8%
Simplified73.8%
Final simplification83.2%
(FPCore (x y z t a b c i)
:precision binary64
(if (<= y -1.45e+33)
(/ x (+ 1.0 (/ a y)))
(if (<= y 2.3e-10)
(/
(+ t (* y (+ 230661.510616 (* y (+ 27464.7644705 (* y z))))))
(+ i (* y (+ c (* y b)))))
(if (<= y 1.9e+62)
(/
(+ (* y (+ (* y (+ (* x y) z)) 27464.7644705)) 230661.510616)
(+ c (* y (+ (* y (+ y a)) b))))
(+ 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 tmp;
if (y <= -1.45e+33) {
tmp = x / (1.0 + (a / y));
} else if (y <= 2.3e-10) {
tmp = (t + (y * (230661.510616 + (y * (27464.7644705 + (y * z)))))) / (i + (y * (c + (y * b))));
} else if (y <= 1.9e+62) {
tmp = ((y * ((y * ((x * y) + z)) + 27464.7644705)) + 230661.510616) / (c + (y * ((y * (y + a)) + b)));
} else {
tmp = x + ((z / y) - (a * (x / y)));
}
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.45d+33)) then
tmp = x / (1.0d0 + (a / y))
else if (y <= 2.3d-10) then
tmp = (t + (y * (230661.510616d0 + (y * (27464.7644705d0 + (y * z)))))) / (i + (y * (c + (y * b))))
else if (y <= 1.9d+62) then
tmp = ((y * ((y * ((x * y) + z)) + 27464.7644705d0)) + 230661.510616d0) / (c + (y * ((y * (y + a)) + b)))
else
tmp = x + ((z / y) - (a * (x / y)))
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.45e+33) {
tmp = x / (1.0 + (a / y));
} else if (y <= 2.3e-10) {
tmp = (t + (y * (230661.510616 + (y * (27464.7644705 + (y * z)))))) / (i + (y * (c + (y * b))));
} else if (y <= 1.9e+62) {
tmp = ((y * ((y * ((x * y) + z)) + 27464.7644705)) + 230661.510616) / (c + (y * ((y * (y + a)) + b)));
} else {
tmp = x + ((z / y) - (a * (x / y)));
}
return tmp;
}
def code(x, y, z, t, a, b, c, i): tmp = 0 if y <= -1.45e+33: tmp = x / (1.0 + (a / y)) elif y <= 2.3e-10: tmp = (t + (y * (230661.510616 + (y * (27464.7644705 + (y * z)))))) / (i + (y * (c + (y * b)))) elif y <= 1.9e+62: tmp = ((y * ((y * ((x * y) + z)) + 27464.7644705)) + 230661.510616) / (c + (y * ((y * (y + a)) + b))) else: tmp = x + ((z / y) - (a * (x / y))) return tmp
function code(x, y, z, t, a, b, c, i) tmp = 0.0 if (y <= -1.45e+33) tmp = Float64(x / Float64(1.0 + Float64(a / y))); elseif (y <= 2.3e-10) 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))))); elseif (y <= 1.9e+62) tmp = Float64(Float64(Float64(y * Float64(Float64(y * Float64(Float64(x * y) + z)) + 27464.7644705)) + 230661.510616) / Float64(c + Float64(y * Float64(Float64(y * Float64(y + a)) + b)))); 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) tmp = 0.0; if (y <= -1.45e+33) tmp = x / (1.0 + (a / y)); elseif (y <= 2.3e-10) tmp = (t + (y * (230661.510616 + (y * (27464.7644705 + (y * z)))))) / (i + (y * (c + (y * b)))); elseif (y <= 1.9e+62) tmp = ((y * ((y * ((x * y) + z)) + 27464.7644705)) + 230661.510616) / (c + (y * ((y * (y + a)) + b))); else tmp = x + ((z / y) - (a * (x / y))); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_, i_] := If[LessEqual[y, -1.45e+33], N[(x / N[(1.0 + N[(a / y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 2.3e-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 * N[(c + N[(y * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 1.9e+62], N[(N[(N[(y * N[(N[(y * N[(N[(x * y), $MachinePrecision] + z), $MachinePrecision]), $MachinePrecision] + 27464.7644705), $MachinePrecision]), $MachinePrecision] + 230661.510616), $MachinePrecision] / N[(c + N[(y * N[(N[(y * N[(y + a), $MachinePrecision]), $MachinePrecision] + b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x + N[(N[(z / y), $MachinePrecision] - N[(a * N[(x / y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -1.45 \cdot 10^{+33}:\\
\;\;\;\;\frac{x}{1 + \frac{a}{y}}\\
\mathbf{elif}\;y \leq 2.3 \cdot 10^{-10}:\\
\;\;\;\;\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{elif}\;y \leq 1.9 \cdot 10^{+62}:\\
\;\;\;\;\frac{y \cdot \left(y \cdot \left(x \cdot y + z\right) + 27464.7644705\right) + 230661.510616}{c + y \cdot \left(y \cdot \left(y + a\right) + b\right)}\\
\mathbf{else}:\\
\;\;\;\;x + \left(\frac{z}{y} - a \cdot \frac{x}{y}\right)\\
\end{array}
\end{array}
if y < -1.45000000000000012e33Initial program 4.2%
fma-define4.2%
fma-define4.2%
fma-define4.2%
fma-define4.2%
fma-define4.2%
fma-define4.2%
fma-define4.2%
Simplified4.2%
Applied egg-rr4.2%
unpow-14.2%
+-commutative4.2%
Simplified4.2%
Taylor expanded in y around inf 58.7%
Taylor expanded in x around inf 71.6%
if -1.45000000000000012e33 < y < 2.30000000000000007e-10Initial program 98.9%
Taylor expanded in x around 0 95.8%
Taylor expanded in y around 0 91.9%
*-commutative94.9%
Simplified91.9%
if 2.30000000000000007e-10 < y < 1.89999999999999992e62Initial program 51.7%
fma-define51.7%
fma-define51.7%
fma-define51.7%
fma-define51.7%
fma-define51.7%
fma-define51.7%
fma-define51.7%
Simplified51.7%
Applied egg-rr51.8%
unpow-151.8%
+-commutative51.8%
Simplified51.8%
Taylor expanded in t around 0 45.1%
associate-/r*51.9%
Simplified51.9%
Taylor expanded in i around 0 65.7%
if 1.89999999999999992e62 < y Initial program 0.4%
fma-define0.4%
fma-define0.4%
fma-define0.4%
fma-define0.4%
fma-define0.4%
fma-define0.4%
fma-define0.4%
Simplified0.4%
Taylor expanded in y around inf 72.3%
associate--l+72.3%
associate-/l*73.8%
Simplified73.8%
Final simplification81.7%
(FPCore (x y z t a b c i)
:precision binary64
(if (<= y -1.26e+35)
(/ x (+ 1.0 (/ a y)))
(if (<= y 2.4e-10)
(/
(+ t (* y (+ 230661.510616 (* y (+ 27464.7644705 (* y z))))))
(+ i (* y (+ c (* y b)))))
(if (<= y 5e+58)
(/
1.0
(*
y
(/
(+ (* y (+ y a)) b)
(+ (* y (+ (* y (+ (* x y) z)) 27464.7644705)) 230661.510616))))
(+ 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 tmp;
if (y <= -1.26e+35) {
tmp = x / (1.0 + (a / y));
} else if (y <= 2.4e-10) {
tmp = (t + (y * (230661.510616 + (y * (27464.7644705 + (y * z)))))) / (i + (y * (c + (y * b))));
} else if (y <= 5e+58) {
tmp = 1.0 / (y * (((y * (y + a)) + b) / ((y * ((y * ((x * y) + z)) + 27464.7644705)) + 230661.510616)));
} else {
tmp = x + ((z / y) - (a * (x / y)));
}
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.26d+35)) then
tmp = x / (1.0d0 + (a / y))
else if (y <= 2.4d-10) then
tmp = (t + (y * (230661.510616d0 + (y * (27464.7644705d0 + (y * z)))))) / (i + (y * (c + (y * b))))
else if (y <= 5d+58) then
tmp = 1.0d0 / (y * (((y * (y + a)) + b) / ((y * ((y * ((x * y) + z)) + 27464.7644705d0)) + 230661.510616d0)))
else
tmp = x + ((z / y) - (a * (x / y)))
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.26e+35) {
tmp = x / (1.0 + (a / y));
} else if (y <= 2.4e-10) {
tmp = (t + (y * (230661.510616 + (y * (27464.7644705 + (y * z)))))) / (i + (y * (c + (y * b))));
} else if (y <= 5e+58) {
tmp = 1.0 / (y * (((y * (y + a)) + b) / ((y * ((y * ((x * y) + z)) + 27464.7644705)) + 230661.510616)));
} else {
tmp = x + ((z / y) - (a * (x / y)));
}
return tmp;
}
def code(x, y, z, t, a, b, c, i): tmp = 0 if y <= -1.26e+35: tmp = x / (1.0 + (a / y)) elif y <= 2.4e-10: tmp = (t + (y * (230661.510616 + (y * (27464.7644705 + (y * z)))))) / (i + (y * (c + (y * b)))) elif y <= 5e+58: tmp = 1.0 / (y * (((y * (y + a)) + b) / ((y * ((y * ((x * y) + z)) + 27464.7644705)) + 230661.510616))) else: tmp = x + ((z / y) - (a * (x / y))) return tmp
function code(x, y, z, t, a, b, c, i) tmp = 0.0 if (y <= -1.26e+35) tmp = Float64(x / Float64(1.0 + Float64(a / y))); elseif (y <= 2.4e-10) 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))))); elseif (y <= 5e+58) tmp = Float64(1.0 / Float64(y * Float64(Float64(Float64(y * Float64(y + a)) + b) / Float64(Float64(y * Float64(Float64(y * Float64(Float64(x * y) + z)) + 27464.7644705)) + 230661.510616)))); 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) tmp = 0.0; if (y <= -1.26e+35) tmp = x / (1.0 + (a / y)); elseif (y <= 2.4e-10) tmp = (t + (y * (230661.510616 + (y * (27464.7644705 + (y * z)))))) / (i + (y * (c + (y * b)))); elseif (y <= 5e+58) tmp = 1.0 / (y * (((y * (y + a)) + b) / ((y * ((y * ((x * y) + z)) + 27464.7644705)) + 230661.510616))); else tmp = x + ((z / y) - (a * (x / y))); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_, i_] := If[LessEqual[y, -1.26e+35], N[(x / N[(1.0 + N[(a / y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 2.4e-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 * N[(c + N[(y * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 5e+58], N[(1.0 / N[(y * N[(N[(N[(y * N[(y + a), $MachinePrecision]), $MachinePrecision] + b), $MachinePrecision] / N[(N[(y * N[(N[(y * N[(N[(x * y), $MachinePrecision] + z), $MachinePrecision]), $MachinePrecision] + 27464.7644705), $MachinePrecision]), $MachinePrecision] + 230661.510616), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x + N[(N[(z / y), $MachinePrecision] - N[(a * N[(x / y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -1.26 \cdot 10^{+35}:\\
\;\;\;\;\frac{x}{1 + \frac{a}{y}}\\
\mathbf{elif}\;y \leq 2.4 \cdot 10^{-10}:\\
\;\;\;\;\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{elif}\;y \leq 5 \cdot 10^{+58}:\\
\;\;\;\;\frac{1}{y \cdot \frac{y \cdot \left(y + a\right) + b}{y \cdot \left(y \cdot \left(x \cdot y + z\right) + 27464.7644705\right) + 230661.510616}}\\
\mathbf{else}:\\
\;\;\;\;x + \left(\frac{z}{y} - a \cdot \frac{x}{y}\right)\\
\end{array}
\end{array}
if y < -1.26e35Initial program 4.2%
fma-define4.2%
fma-define4.2%
fma-define4.2%
fma-define4.2%
fma-define4.2%
fma-define4.2%
fma-define4.2%
Simplified4.2%
Applied egg-rr4.2%
unpow-14.2%
+-commutative4.2%
Simplified4.2%
Taylor expanded in y around inf 58.7%
Taylor expanded in x around inf 71.6%
if -1.26e35 < y < 2.4e-10Initial program 98.9%
Taylor expanded in x around 0 95.8%
Taylor expanded in y around 0 91.9%
*-commutative94.9%
Simplified91.9%
if 2.4e-10 < y < 4.99999999999999986e58Initial program 55.2%
fma-define55.2%
fma-define55.2%
fma-define55.2%
fma-define55.2%
fma-define55.2%
fma-define55.2%
fma-define55.2%
Simplified55.2%
Applied egg-rr55.4%
unpow-155.4%
+-commutative55.4%
Simplified55.4%
Taylor expanded in t around 0 48.2%
associate-/r*55.4%
Simplified55.4%
Taylor expanded in i around 0 70.3%
Taylor expanded in c around 0 62.9%
associate-/l*70.3%
Simplified70.3%
if 4.99999999999999986e58 < y Initial program 0.5%
fma-define0.5%
fma-define0.5%
fma-define0.5%
fma-define0.5%
fma-define0.5%
fma-define0.5%
fma-define0.5%
Simplified0.5%
Taylor expanded in y around inf 71.1%
associate--l+71.1%
associate-/l*72.6%
Simplified72.6%
Final simplification81.7%
(FPCore (x y z t a b c i)
:precision binary64
(if (<= y -1.58e+32)
(/ x (+ 1.0 (/ a y)))
(if (<= y 2.4e-10)
(/
(+ t (* y (+ 230661.510616 (* y (+ 27464.7644705 (* y z))))))
(+ i (* y (+ c (* y b)))))
(+ 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 tmp;
if (y <= -1.58e+32) {
tmp = x / (1.0 + (a / y));
} else if (y <= 2.4e-10) {
tmp = (t + (y * (230661.510616 + (y * (27464.7644705 + (y * z)))))) / (i + (y * (c + (y * b))));
} else {
tmp = x + ((z / y) - (a * (x / y)));
}
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.58d+32)) then
tmp = x / (1.0d0 + (a / y))
else if (y <= 2.4d-10) then
tmp = (t + (y * (230661.510616d0 + (y * (27464.7644705d0 + (y * z)))))) / (i + (y * (c + (y * b))))
else
tmp = x + ((z / y) - (a * (x / y)))
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.58e+32) {
tmp = x / (1.0 + (a / y));
} else if (y <= 2.4e-10) {
tmp = (t + (y * (230661.510616 + (y * (27464.7644705 + (y * z)))))) / (i + (y * (c + (y * b))));
} else {
tmp = x + ((z / y) - (a * (x / y)));
}
return tmp;
}
def code(x, y, z, t, a, b, c, i): tmp = 0 if y <= -1.58e+32: tmp = x / (1.0 + (a / y)) elif y <= 2.4e-10: tmp = (t + (y * (230661.510616 + (y * (27464.7644705 + (y * z)))))) / (i + (y * (c + (y * b)))) else: tmp = x + ((z / y) - (a * (x / y))) return tmp
function code(x, y, z, t, a, b, c, i) tmp = 0.0 if (y <= -1.58e+32) tmp = Float64(x / Float64(1.0 + Float64(a / y))); elseif (y <= 2.4e-10) 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 = 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) tmp = 0.0; if (y <= -1.58e+32) tmp = x / (1.0 + (a / y)); elseif (y <= 2.4e-10) tmp = (t + (y * (230661.510616 + (y * (27464.7644705 + (y * z)))))) / (i + (y * (c + (y * b)))); else tmp = x + ((z / y) - (a * (x / y))); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_, i_] := If[LessEqual[y, -1.58e+32], N[(x / N[(1.0 + N[(a / y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 2.4e-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 * N[(c + N[(y * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x + N[(N[(z / y), $MachinePrecision] - N[(a * N[(x / y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -1.58 \cdot 10^{+32}:\\
\;\;\;\;\frac{x}{1 + \frac{a}{y}}\\
\mathbf{elif}\;y \leq 2.4 \cdot 10^{-10}:\\
\;\;\;\;\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}:\\
\;\;\;\;x + \left(\frac{z}{y} - a \cdot \frac{x}{y}\right)\\
\end{array}
\end{array}
if y < -1.58000000000000006e32Initial program 4.2%
fma-define4.2%
fma-define4.2%
fma-define4.2%
fma-define4.2%
fma-define4.2%
fma-define4.2%
fma-define4.2%
Simplified4.2%
Applied egg-rr4.2%
unpow-14.2%
+-commutative4.2%
Simplified4.2%
Taylor expanded in y around inf 58.7%
Taylor expanded in x around inf 71.6%
if -1.58000000000000006e32 < y < 2.4e-10Initial program 98.9%
Taylor expanded in x around 0 95.8%
Taylor expanded in y around 0 91.9%
*-commutative94.9%
Simplified91.9%
if 2.4e-10 < y Initial program 10.2%
fma-define10.2%
fma-define10.2%
fma-define10.2%
fma-define10.2%
fma-define10.2%
fma-define10.2%
fma-define10.2%
Simplified10.2%
Taylor expanded in y around inf 62.9%
associate--l+62.9%
associate-/l*64.1%
Simplified64.1%
Final simplification79.4%
(FPCore (x y z t a b c i)
:precision binary64
(if (<= y -6.8e+33)
(/ x (+ 1.0 (/ a y)))
(if (<= y 2.4e-10)
(/ (+ t (* y 230661.510616)) (+ i (* y (+ c (* y (+ (* y (+ y a)) b))))))
(+ 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 tmp;
if (y <= -6.8e+33) {
tmp = x / (1.0 + (a / y));
} else if (y <= 2.4e-10) {
tmp = (t + (y * 230661.510616)) / (i + (y * (c + (y * ((y * (y + a)) + b)))));
} else {
tmp = x + ((z / y) - (a * (x / y)));
}
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.8d+33)) then
tmp = x / (1.0d0 + (a / y))
else if (y <= 2.4d-10) then
tmp = (t + (y * 230661.510616d0)) / (i + (y * (c + (y * ((y * (y + a)) + b)))))
else
tmp = x + ((z / y) - (a * (x / y)))
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.8e+33) {
tmp = x / (1.0 + (a / y));
} else if (y <= 2.4e-10) {
tmp = (t + (y * 230661.510616)) / (i + (y * (c + (y * ((y * (y + a)) + b)))));
} else {
tmp = x + ((z / y) - (a * (x / y)));
}
return tmp;
}
def code(x, y, z, t, a, b, c, i): tmp = 0 if y <= -6.8e+33: tmp = x / (1.0 + (a / y)) elif y <= 2.4e-10: tmp = (t + (y * 230661.510616)) / (i + (y * (c + (y * ((y * (y + a)) + b))))) else: tmp = x + ((z / y) - (a * (x / y))) return tmp
function code(x, y, z, t, a, b, c, i) tmp = 0.0 if (y <= -6.8e+33) tmp = Float64(x / Float64(1.0 + Float64(a / y))); elseif (y <= 2.4e-10) tmp = Float64(Float64(t + Float64(y * 230661.510616)) / Float64(i + Float64(y * Float64(c + Float64(y * Float64(Float64(y * Float64(y + a)) + b)))))); 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) tmp = 0.0; if (y <= -6.8e+33) tmp = x / (1.0 + (a / y)); elseif (y <= 2.4e-10) tmp = (t + (y * 230661.510616)) / (i + (y * (c + (y * ((y * (y + a)) + b))))); else tmp = x + ((z / y) - (a * (x / y))); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_, i_] := If[LessEqual[y, -6.8e+33], N[(x / N[(1.0 + N[(a / y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 2.4e-10], 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], N[(x + N[(N[(z / y), $MachinePrecision] - N[(a * N[(x / y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -6.8 \cdot 10^{+33}:\\
\;\;\;\;\frac{x}{1 + \frac{a}{y}}\\
\mathbf{elif}\;y \leq 2.4 \cdot 10^{-10}:\\
\;\;\;\;\frac{t + y \cdot 230661.510616}{i + y \cdot \left(c + y \cdot \left(y \cdot \left(y + a\right) + b\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;x + \left(\frac{z}{y} - a \cdot \frac{x}{y}\right)\\
\end{array}
\end{array}
if y < -6.7999999999999999e33Initial program 4.2%
fma-define4.2%
fma-define4.2%
fma-define4.2%
fma-define4.2%
fma-define4.2%
fma-define4.2%
fma-define4.2%
Simplified4.2%
Applied egg-rr4.2%
unpow-14.2%
+-commutative4.2%
Simplified4.2%
Taylor expanded in y around inf 58.7%
Taylor expanded in x around inf 71.6%
if -6.7999999999999999e33 < y < 2.4e-10Initial program 98.9%
Taylor expanded in y around 0 90.5%
*-commutative90.5%
Simplified90.5%
if 2.4e-10 < y Initial program 10.2%
fma-define10.2%
fma-define10.2%
fma-define10.2%
fma-define10.2%
fma-define10.2%
fma-define10.2%
fma-define10.2%
Simplified10.2%
Taylor expanded in y around inf 62.9%
associate--l+62.9%
associate-/l*64.1%
Simplified64.1%
Final simplification78.7%
(FPCore (x y z t a b c i)
:precision binary64
(if (<= y -1.78e+31)
(/ x (+ 1.0 (/ a y)))
(if (<= y 2.4e-10)
(/ t (+ i (* y (+ c (* y (+ (* y (+ y a)) b))))))
(+ 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 tmp;
if (y <= -1.78e+31) {
tmp = x / (1.0 + (a / y));
} else if (y <= 2.4e-10) {
tmp = t / (i + (y * (c + (y * ((y * (y + a)) + b)))));
} else {
tmp = x + ((z / y) - (a * (x / y)));
}
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.78d+31)) then
tmp = x / (1.0d0 + (a / y))
else if (y <= 2.4d-10) then
tmp = t / (i + (y * (c + (y * ((y * (y + a)) + b)))))
else
tmp = x + ((z / y) - (a * (x / y)))
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.78e+31) {
tmp = x / (1.0 + (a / y));
} else if (y <= 2.4e-10) {
tmp = t / (i + (y * (c + (y * ((y * (y + a)) + b)))));
} else {
tmp = x + ((z / y) - (a * (x / y)));
}
return tmp;
}
def code(x, y, z, t, a, b, c, i): tmp = 0 if y <= -1.78e+31: tmp = x / (1.0 + (a / y)) elif y <= 2.4e-10: tmp = t / (i + (y * (c + (y * ((y * (y + a)) + b))))) else: tmp = x + ((z / y) - (a * (x / y))) return tmp
function code(x, y, z, t, a, b, c, i) tmp = 0.0 if (y <= -1.78e+31) tmp = Float64(x / Float64(1.0 + Float64(a / y))); elseif (y <= 2.4e-10) tmp = Float64(t / Float64(i + Float64(y * Float64(c + Float64(y * Float64(Float64(y * Float64(y + a)) + b)))))); 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) tmp = 0.0; if (y <= -1.78e+31) tmp = x / (1.0 + (a / y)); elseif (y <= 2.4e-10) tmp = t / (i + (y * (c + (y * ((y * (y + a)) + b))))); else tmp = x + ((z / y) - (a * (x / y))); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_, i_] := If[LessEqual[y, -1.78e+31], N[(x / N[(1.0 + N[(a / y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 2.4e-10], N[(t / N[(i + N[(y * N[(c + N[(y * N[(N[(y * N[(y + a), $MachinePrecision]), $MachinePrecision] + b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x + N[(N[(z / y), $MachinePrecision] - N[(a * N[(x / y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -1.78 \cdot 10^{+31}:\\
\;\;\;\;\frac{x}{1 + \frac{a}{y}}\\
\mathbf{elif}\;y \leq 2.4 \cdot 10^{-10}:\\
\;\;\;\;\frac{t}{i + y \cdot \left(c + y \cdot \left(y \cdot \left(y + a\right) + b\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;x + \left(\frac{z}{y} - a \cdot \frac{x}{y}\right)\\
\end{array}
\end{array}
if y < -1.7800000000000001e31Initial program 4.2%
fma-define4.2%
fma-define4.2%
fma-define4.2%
fma-define4.2%
fma-define4.2%
fma-define4.2%
fma-define4.2%
Simplified4.2%
Applied egg-rr4.2%
unpow-14.2%
+-commutative4.2%
Simplified4.2%
Taylor expanded in y around inf 58.7%
Taylor expanded in x around inf 71.6%
if -1.7800000000000001e31 < y < 2.4e-10Initial program 98.9%
fma-define98.9%
fma-define98.9%
fma-define98.9%
fma-define98.9%
fma-define98.9%
fma-define98.9%
fma-define98.9%
Simplified98.9%
Taylor expanded in t around inf 75.7%
if 2.4e-10 < y Initial program 10.2%
fma-define10.2%
fma-define10.2%
fma-define10.2%
fma-define10.2%
fma-define10.2%
fma-define10.2%
fma-define10.2%
Simplified10.2%
Taylor expanded in y around inf 62.9%
associate--l+62.9%
associate-/l*64.1%
Simplified64.1%
Final simplification71.5%
(FPCore (x y z t a b c i)
:precision binary64
(if (<= y -1.9e+17)
(/ x (+ 1.0 (/ a y)))
(if (<= y 2.4e-10)
(/ (+ t (* y (+ 230661.510616 (* y (+ 27464.7644705 (* y z)))))) i)
(+ 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 tmp;
if (y <= -1.9e+17) {
tmp = x / (1.0 + (a / y));
} else if (y <= 2.4e-10) {
tmp = (t + (y * (230661.510616 + (y * (27464.7644705 + (y * z)))))) / i;
} else {
tmp = x + ((z / y) - (a * (x / y)));
}
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.9d+17)) then
tmp = x / (1.0d0 + (a / y))
else if (y <= 2.4d-10) then
tmp = (t + (y * (230661.510616d0 + (y * (27464.7644705d0 + (y * z)))))) / i
else
tmp = x + ((z / y) - (a * (x / y)))
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.9e+17) {
tmp = x / (1.0 + (a / y));
} else if (y <= 2.4e-10) {
tmp = (t + (y * (230661.510616 + (y * (27464.7644705 + (y * z)))))) / i;
} else {
tmp = x + ((z / y) - (a * (x / y)));
}
return tmp;
}
def code(x, y, z, t, a, b, c, i): tmp = 0 if y <= -1.9e+17: tmp = x / (1.0 + (a / y)) elif y <= 2.4e-10: tmp = (t + (y * (230661.510616 + (y * (27464.7644705 + (y * z)))))) / i else: tmp = x + ((z / y) - (a * (x / y))) return tmp
function code(x, y, z, t, a, b, c, i) tmp = 0.0 if (y <= -1.9e+17) tmp = Float64(x / Float64(1.0 + Float64(a / y))); elseif (y <= 2.4e-10) tmp = Float64(Float64(t + Float64(y * Float64(230661.510616 + Float64(y * Float64(27464.7644705 + Float64(y * z)))))) / i); 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) tmp = 0.0; if (y <= -1.9e+17) tmp = x / (1.0 + (a / y)); elseif (y <= 2.4e-10) tmp = (t + (y * (230661.510616 + (y * (27464.7644705 + (y * z)))))) / i; else tmp = x + ((z / y) - (a * (x / y))); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_, i_] := If[LessEqual[y, -1.9e+17], N[(x / N[(1.0 + N[(a / y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 2.4e-10], N[(N[(t + N[(y * N[(230661.510616 + N[(y * N[(27464.7644705 + N[(y * z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / i), $MachinePrecision], N[(x + N[(N[(z / y), $MachinePrecision] - N[(a * N[(x / y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -1.9 \cdot 10^{+17}:\\
\;\;\;\;\frac{x}{1 + \frac{a}{y}}\\
\mathbf{elif}\;y \leq 2.4 \cdot 10^{-10}:\\
\;\;\;\;\frac{t + y \cdot \left(230661.510616 + y \cdot \left(27464.7644705 + y \cdot z\right)\right)}{i}\\
\mathbf{else}:\\
\;\;\;\;x + \left(\frac{z}{y} - a \cdot \frac{x}{y}\right)\\
\end{array}
\end{array}
if y < -1.9e17Initial program 7.4%
fma-define7.4%
fma-define7.4%
fma-define7.4%
fma-define7.4%
fma-define7.4%
fma-define7.4%
fma-define7.4%
Simplified7.4%
Applied egg-rr7.4%
unpow-17.4%
+-commutative7.4%
Simplified7.4%
Taylor expanded in y around inf 56.9%
Taylor expanded in x around inf 69.4%
if -1.9e17 < y < 2.4e-10Initial program 98.9%
Taylor expanded in x around 0 95.8%
Taylor expanded in i around inf 61.9%
if 2.4e-10 < y Initial program 10.2%
fma-define10.2%
fma-define10.2%
fma-define10.2%
fma-define10.2%
fma-define10.2%
fma-define10.2%
fma-define10.2%
Simplified10.2%
Taylor expanded in y around inf 62.9%
associate--l+62.9%
associate-/l*64.1%
Simplified64.1%
(FPCore (x y z t a b c i)
:precision binary64
(if (<= y -0.19)
(/ x (+ 1.0 (/ a y)))
(if (<= y 1.5e-10)
(/ (+ t (* y 230661.510616)) i)
(+ 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 tmp;
if (y <= -0.19) {
tmp = x / (1.0 + (a / y));
} else if (y <= 1.5e-10) {
tmp = (t + (y * 230661.510616)) / i;
} else {
tmp = x + ((z / y) - (a * (x / y)));
}
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 <= (-0.19d0)) then
tmp = x / (1.0d0 + (a / y))
else if (y <= 1.5d-10) then
tmp = (t + (y * 230661.510616d0)) / i
else
tmp = x + ((z / y) - (a * (x / y)))
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 <= -0.19) {
tmp = x / (1.0 + (a / y));
} else if (y <= 1.5e-10) {
tmp = (t + (y * 230661.510616)) / i;
} else {
tmp = x + ((z / y) - (a * (x / y)));
}
return tmp;
}
def code(x, y, z, t, a, b, c, i): tmp = 0 if y <= -0.19: tmp = x / (1.0 + (a / y)) elif y <= 1.5e-10: tmp = (t + (y * 230661.510616)) / i else: tmp = x + ((z / y) - (a * (x / y))) return tmp
function code(x, y, z, t, a, b, c, i) tmp = 0.0 if (y <= -0.19) tmp = Float64(x / Float64(1.0 + Float64(a / y))); elseif (y <= 1.5e-10) tmp = Float64(Float64(t + Float64(y * 230661.510616)) / i); 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) tmp = 0.0; if (y <= -0.19) tmp = x / (1.0 + (a / y)); elseif (y <= 1.5e-10) tmp = (t + (y * 230661.510616)) / i; else tmp = x + ((z / y) - (a * (x / y))); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_, i_] := If[LessEqual[y, -0.19], N[(x / N[(1.0 + N[(a / y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 1.5e-10], N[(N[(t + N[(y * 230661.510616), $MachinePrecision]), $MachinePrecision] / i), $MachinePrecision], N[(x + N[(N[(z / y), $MachinePrecision] - N[(a * N[(x / y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -0.19:\\
\;\;\;\;\frac{x}{1 + \frac{a}{y}}\\
\mathbf{elif}\;y \leq 1.5 \cdot 10^{-10}:\\
\;\;\;\;\frac{t + y \cdot 230661.510616}{i}\\
\mathbf{else}:\\
\;\;\;\;x + \left(\frac{z}{y} - a \cdot \frac{x}{y}\right)\\
\end{array}
\end{array}
if y < -0.19Initial program 8.9%
fma-define8.9%
fma-define8.9%
fma-define8.9%
fma-define8.9%
fma-define8.9%
fma-define8.9%
fma-define8.9%
Simplified8.9%
Applied egg-rr8.9%
unpow-18.9%
+-commutative8.9%
Simplified8.9%
Taylor expanded in y around inf 55.1%
Taylor expanded in x around inf 67.2%
if -0.19 < y < 1.5e-10Initial program 99.6%
Taylor expanded in x around 0 96.5%
Taylor expanded in i around inf 62.9%
Taylor expanded in y around 0 59.9%
if 1.5e-10 < y Initial program 10.2%
fma-define10.2%
fma-define10.2%
fma-define10.2%
fma-define10.2%
fma-define10.2%
fma-define10.2%
fma-define10.2%
Simplified10.2%
Taylor expanded in y around inf 62.9%
associate--l+62.9%
associate-/l*64.1%
Simplified64.1%
Final simplification62.9%
(FPCore (x y z t a b c i) :precision binary64 (if (or (<= y -0.85) (not (<= y 8e-14))) (/ x (+ 1.0 (/ a y))) (/ (+ t (* y 230661.510616)) i)))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double tmp;
if ((y <= -0.85) || !(y <= 8e-14)) {
tmp = x / (1.0 + (a / y));
} else {
tmp = (t + (y * 230661.510616)) / i;
}
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 <= (-0.85d0)) .or. (.not. (y <= 8d-14))) then
tmp = x / (1.0d0 + (a / y))
else
tmp = (t + (y * 230661.510616d0)) / i
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 <= -0.85) || !(y <= 8e-14)) {
tmp = x / (1.0 + (a / y));
} else {
tmp = (t + (y * 230661.510616)) / i;
}
return tmp;
}
def code(x, y, z, t, a, b, c, i): tmp = 0 if (y <= -0.85) or not (y <= 8e-14): tmp = x / (1.0 + (a / y)) else: tmp = (t + (y * 230661.510616)) / i return tmp
function code(x, y, z, t, a, b, c, i) tmp = 0.0 if ((y <= -0.85) || !(y <= 8e-14)) tmp = Float64(x / Float64(1.0 + Float64(a / y))); else tmp = Float64(Float64(t + Float64(y * 230661.510616)) / i); end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c, i) tmp = 0.0; if ((y <= -0.85) || ~((y <= 8e-14))) tmp = x / (1.0 + (a / y)); else tmp = (t + (y * 230661.510616)) / i; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_, i_] := If[Or[LessEqual[y, -0.85], N[Not[LessEqual[y, 8e-14]], $MachinePrecision]], N[(x / N[(1.0 + N[(a / y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(t + N[(y * 230661.510616), $MachinePrecision]), $MachinePrecision] / i), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -0.85 \lor \neg \left(y \leq 8 \cdot 10^{-14}\right):\\
\;\;\;\;\frac{x}{1 + \frac{a}{y}}\\
\mathbf{else}:\\
\;\;\;\;\frac{t + y \cdot 230661.510616}{i}\\
\end{array}
\end{array}
if y < -0.849999999999999978 or 7.99999999999999999e-14 < y Initial program 9.6%
fma-define9.6%
fma-define9.6%
fma-define9.6%
fma-define9.6%
fma-define9.6%
fma-define9.6%
fma-define9.6%
Simplified9.6%
Applied egg-rr9.6%
unpow-19.6%
+-commutative9.6%
Simplified9.6%
Taylor expanded in y around inf 48.8%
Taylor expanded in x around inf 61.9%
if -0.849999999999999978 < y < 7.99999999999999999e-14Initial program 99.6%
Taylor expanded in x around 0 96.5%
Taylor expanded in i around inf 62.9%
Taylor expanded in y around 0 59.9%
Final simplification61.0%
(FPCore (x y z t a b c i) :precision binary64 (if (or (<= y -5.7e-36) (not (<= y 1.9e-10))) (/ x (+ 1.0 (/ a y))) (/ t i)))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double tmp;
if ((y <= -5.7e-36) || !(y <= 1.9e-10)) {
tmp = x / (1.0 + (a / y));
} else {
tmp = t / i;
}
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 <= (-5.7d-36)) .or. (.not. (y <= 1.9d-10))) then
tmp = x / (1.0d0 + (a / y))
else
tmp = t / i
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 <= -5.7e-36) || !(y <= 1.9e-10)) {
tmp = x / (1.0 + (a / y));
} else {
tmp = t / i;
}
return tmp;
}
def code(x, y, z, t, a, b, c, i): tmp = 0 if (y <= -5.7e-36) or not (y <= 1.9e-10): tmp = x / (1.0 + (a / y)) else: tmp = t / i return tmp
function code(x, y, z, t, a, b, c, i) tmp = 0.0 if ((y <= -5.7e-36) || !(y <= 1.9e-10)) tmp = Float64(x / Float64(1.0 + Float64(a / y))); else tmp = Float64(t / i); end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c, i) tmp = 0.0; if ((y <= -5.7e-36) || ~((y <= 1.9e-10))) tmp = x / (1.0 + (a / y)); else tmp = t / i; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_, i_] := If[Or[LessEqual[y, -5.7e-36], N[Not[LessEqual[y, 1.9e-10]], $MachinePrecision]], N[(x / N[(1.0 + N[(a / y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(t / i), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -5.7 \cdot 10^{-36} \lor \neg \left(y \leq 1.9 \cdot 10^{-10}\right):\\
\;\;\;\;\frac{x}{1 + \frac{a}{y}}\\
\mathbf{else}:\\
\;\;\;\;\frac{t}{i}\\
\end{array}
\end{array}
if y < -5.6999999999999999e-36 or 1.8999999999999999e-10 < y Initial program 15.8%
fma-define15.8%
fma-define15.8%
fma-define15.8%
fma-define15.8%
fma-define15.8%
fma-define15.8%
fma-define15.8%
Simplified15.8%
Applied egg-rr15.8%
unpow-115.8%
+-commutative15.8%
Simplified15.8%
Taylor expanded in y around inf 45.7%
Taylor expanded in x around inf 57.9%
if -5.6999999999999999e-36 < y < 1.8999999999999999e-10Initial program 99.6%
fma-define99.6%
fma-define99.6%
fma-define99.6%
fma-define99.6%
fma-define99.6%
fma-define99.6%
fma-define99.6%
Simplified99.6%
Taylor expanded in y around 0 56.8%
Final simplification57.4%
(FPCore (x y z t a b c i) :precision binary64 (if (or (<= y -2.6e-12) (not (<= y 2.4e-10))) (* x (- 1.0 (/ a y))) (/ t i)))
double code(double x, double y, double z, double t, double a, double b, double c, double i) {
double tmp;
if ((y <= -2.6e-12) || !(y <= 2.4e-10)) {
tmp = x * (1.0 - (a / y));
} else {
tmp = t / i;
}
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 <= (-2.6d-12)) .or. (.not. (y <= 2.4d-10))) then
tmp = x * (1.0d0 - (a / y))
else
tmp = t / i
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 <= -2.6e-12) || !(y <= 2.4e-10)) {
tmp = x * (1.0 - (a / y));
} else {
tmp = t / i;
}
return tmp;
}
def code(x, y, z, t, a, b, c, i): tmp = 0 if (y <= -2.6e-12) or not (y <= 2.4e-10): tmp = x * (1.0 - (a / y)) else: tmp = t / i return tmp
function code(x, y, z, t, a, b, c, i) tmp = 0.0 if ((y <= -2.6e-12) || !(y <= 2.4e-10)) tmp = Float64(x * Float64(1.0 - Float64(a / y))); else tmp = Float64(t / i); end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c, i) tmp = 0.0; if ((y <= -2.6e-12) || ~((y <= 2.4e-10))) tmp = x * (1.0 - (a / y)); else tmp = t / i; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_, i_] := If[Or[LessEqual[y, -2.6e-12], N[Not[LessEqual[y, 2.4e-10]], $MachinePrecision]], N[(x * N[(1.0 - N[(a / y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(t / i), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -2.6 \cdot 10^{-12} \lor \neg \left(y \leq 2.4 \cdot 10^{-10}\right):\\
\;\;\;\;x \cdot \left(1 - \frac{a}{y}\right)\\
\mathbf{else}:\\
\;\;\;\;\frac{t}{i}\\
\end{array}
\end{array}
if y < -2.59999999999999983e-12 or 2.4e-10 < y Initial program 11.6%
fma-define11.6%
fma-define11.6%
fma-define11.6%
fma-define11.6%
fma-define11.6%
fma-define11.6%
fma-define11.6%
Simplified11.6%
Taylor expanded in x around inf 3.6%
associate-/l*8.0%
+-commutative8.0%
fma-undefine8.0%
+-commutative8.0%
fma-undefine8.0%
+-commutative8.0%
+-commutative8.0%
fma-undefine8.0%
+-commutative8.0%
Simplified8.0%
Taylor expanded in y around inf 49.3%
mul-1-neg49.3%
sub-neg49.3%
Simplified49.3%
if -2.59999999999999983e-12 < y < 2.4e-10Initial program 99.6%
fma-define99.6%
fma-define99.6%
fma-define99.6%
fma-define99.6%
fma-define99.6%
fma-define99.6%
fma-define99.6%
Simplified99.6%
Taylor expanded in y around 0 53.6%
Final simplification51.3%
(FPCore (x y z t a b c i) :precision binary64 (if (<= y -1.95e-13) x (if (<= y 3.4e-11) (/ 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 <= -1.95e-13) {
tmp = x;
} else if (y <= 3.4e-11) {
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 <= (-1.95d-13)) then
tmp = x
else if (y <= 3.4d-11) 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 <= -1.95e-13) {
tmp = x;
} else if (y <= 3.4e-11) {
tmp = t / i;
} else {
tmp = x;
}
return tmp;
}
def code(x, y, z, t, a, b, c, i): tmp = 0 if y <= -1.95e-13: tmp = x elif y <= 3.4e-11: tmp = t / i else: tmp = x return tmp
function code(x, y, z, t, a, b, c, i) tmp = 0.0 if (y <= -1.95e-13) tmp = x; elseif (y <= 3.4e-11) 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 <= -1.95e-13) tmp = x; elseif (y <= 3.4e-11) tmp = t / i; else tmp = x; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_, i_] := If[LessEqual[y, -1.95e-13], x, If[LessEqual[y, 3.4e-11], N[(t / i), $MachinePrecision], x]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -1.95 \cdot 10^{-13}:\\
\;\;\;\;x\\
\mathbf{elif}\;y \leq 3.4 \cdot 10^{-11}:\\
\;\;\;\;\frac{t}{i}\\
\mathbf{else}:\\
\;\;\;\;x\\
\end{array}
\end{array}
if y < -1.95000000000000002e-13 or 3.3999999999999999e-11 < y Initial program 11.6%
fma-define11.6%
fma-define11.6%
fma-define11.6%
fma-define11.6%
fma-define11.6%
fma-define11.6%
fma-define11.6%
Simplified11.6%
Taylor expanded in y around inf 49.1%
if -1.95000000000000002e-13 < y < 3.3999999999999999e-11Initial program 99.6%
fma-define99.6%
fma-define99.6%
fma-define99.6%
fma-define99.6%
fma-define99.6%
fma-define99.6%
fma-define99.6%
Simplified99.6%
Taylor expanded in y around 0 53.6%
(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 52.1%
fma-define52.1%
fma-define52.1%
fma-define52.1%
fma-define52.1%
fma-define52.1%
fma-define52.1%
fma-define52.1%
Simplified52.1%
Taylor expanded in y around inf 28.4%
herbie shell --seed 2024151
(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)))