
(FPCore (x y z t a b)
:precision binary64
(+
x
(/
(*
y
(+ (* (+ (* (+ (* (+ (* z 3.13060547623) 11.1667541262) z) t) z) a) z) b))
(+
(* (+ (* (+ (* (+ z 15.234687407) z) 31.4690115749) z) 11.9400905721) z)
0.607771387771))))
double code(double x, double y, double z, double t, double a, double b) {
return x + ((y * ((((((((z * 3.13060547623) + 11.1667541262) * z) + t) * z) + a) * z) + b)) / (((((((z + 15.234687407) * z) + 31.4690115749) * z) + 11.9400905721) * z) + 0.607771387771));
}
real(8) function code(x, y, z, t, a, b)
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
code = x + ((y * ((((((((z * 3.13060547623d0) + 11.1667541262d0) * z) + t) * z) + a) * z) + b)) / (((((((z + 15.234687407d0) * z) + 31.4690115749d0) * z) + 11.9400905721d0) * z) + 0.607771387771d0))
end function
public static double code(double x, double y, double z, double t, double a, double b) {
return x + ((y * ((((((((z * 3.13060547623) + 11.1667541262) * z) + t) * z) + a) * z) + b)) / (((((((z + 15.234687407) * z) + 31.4690115749) * z) + 11.9400905721) * z) + 0.607771387771));
}
def code(x, y, z, t, a, b): return x + ((y * ((((((((z * 3.13060547623) + 11.1667541262) * z) + t) * z) + a) * z) + b)) / (((((((z + 15.234687407) * z) + 31.4690115749) * z) + 11.9400905721) * z) + 0.607771387771))
function code(x, y, z, t, a, b) return Float64(x + Float64(Float64(y * Float64(Float64(Float64(Float64(Float64(Float64(Float64(Float64(z * 3.13060547623) + 11.1667541262) * z) + t) * z) + a) * z) + b)) / Float64(Float64(Float64(Float64(Float64(Float64(Float64(z + 15.234687407) * z) + 31.4690115749) * z) + 11.9400905721) * z) + 0.607771387771))) end
function tmp = code(x, y, z, t, a, b) tmp = x + ((y * ((((((((z * 3.13060547623) + 11.1667541262) * z) + t) * z) + a) * z) + b)) / (((((((z + 15.234687407) * z) + 31.4690115749) * z) + 11.9400905721) * z) + 0.607771387771)); end
code[x_, y_, z_, t_, a_, b_] := N[(x + N[(N[(y * N[(N[(N[(N[(N[(N[(N[(N[(z * 3.13060547623), $MachinePrecision] + 11.1667541262), $MachinePrecision] * z), $MachinePrecision] + t), $MachinePrecision] * z), $MachinePrecision] + a), $MachinePrecision] * z), $MachinePrecision] + b), $MachinePrecision]), $MachinePrecision] / N[(N[(N[(N[(N[(N[(N[(z + 15.234687407), $MachinePrecision] * z), $MachinePrecision] + 31.4690115749), $MachinePrecision] * z), $MachinePrecision] + 11.9400905721), $MachinePrecision] * z), $MachinePrecision] + 0.607771387771), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
x + \frac{y \cdot \left(\left(\left(\left(z \cdot 3.13060547623 + 11.1667541262\right) \cdot z + t\right) \cdot z + a\right) \cdot z + b\right)}{\left(\left(\left(z + 15.234687407\right) \cdot z + 31.4690115749\right) \cdot z + 11.9400905721\right) \cdot z + 0.607771387771}
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 17 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (x y z t a b)
:precision binary64
(+
x
(/
(*
y
(+ (* (+ (* (+ (* (+ (* z 3.13060547623) 11.1667541262) z) t) z) a) z) b))
(+
(* (+ (* (+ (* (+ z 15.234687407) z) 31.4690115749) z) 11.9400905721) z)
0.607771387771))))
double code(double x, double y, double z, double t, double a, double b) {
return x + ((y * ((((((((z * 3.13060547623) + 11.1667541262) * z) + t) * z) + a) * z) + b)) / (((((((z + 15.234687407) * z) + 31.4690115749) * z) + 11.9400905721) * z) + 0.607771387771));
}
real(8) function code(x, y, z, t, a, b)
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
code = x + ((y * ((((((((z * 3.13060547623d0) + 11.1667541262d0) * z) + t) * z) + a) * z) + b)) / (((((((z + 15.234687407d0) * z) + 31.4690115749d0) * z) + 11.9400905721d0) * z) + 0.607771387771d0))
end function
public static double code(double x, double y, double z, double t, double a, double b) {
return x + ((y * ((((((((z * 3.13060547623) + 11.1667541262) * z) + t) * z) + a) * z) + b)) / (((((((z + 15.234687407) * z) + 31.4690115749) * z) + 11.9400905721) * z) + 0.607771387771));
}
def code(x, y, z, t, a, b): return x + ((y * ((((((((z * 3.13060547623) + 11.1667541262) * z) + t) * z) + a) * z) + b)) / (((((((z + 15.234687407) * z) + 31.4690115749) * z) + 11.9400905721) * z) + 0.607771387771))
function code(x, y, z, t, a, b) return Float64(x + Float64(Float64(y * Float64(Float64(Float64(Float64(Float64(Float64(Float64(Float64(z * 3.13060547623) + 11.1667541262) * z) + t) * z) + a) * z) + b)) / Float64(Float64(Float64(Float64(Float64(Float64(Float64(z + 15.234687407) * z) + 31.4690115749) * z) + 11.9400905721) * z) + 0.607771387771))) end
function tmp = code(x, y, z, t, a, b) tmp = x + ((y * ((((((((z * 3.13060547623) + 11.1667541262) * z) + t) * z) + a) * z) + b)) / (((((((z + 15.234687407) * z) + 31.4690115749) * z) + 11.9400905721) * z) + 0.607771387771)); end
code[x_, y_, z_, t_, a_, b_] := N[(x + N[(N[(y * N[(N[(N[(N[(N[(N[(N[(N[(z * 3.13060547623), $MachinePrecision] + 11.1667541262), $MachinePrecision] * z), $MachinePrecision] + t), $MachinePrecision] * z), $MachinePrecision] + a), $MachinePrecision] * z), $MachinePrecision] + b), $MachinePrecision]), $MachinePrecision] / N[(N[(N[(N[(N[(N[(N[(z + 15.234687407), $MachinePrecision] * z), $MachinePrecision] + 31.4690115749), $MachinePrecision] * z), $MachinePrecision] + 11.9400905721), $MachinePrecision] * z), $MachinePrecision] + 0.607771387771), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
x + \frac{y \cdot \left(\left(\left(\left(z \cdot 3.13060547623 + 11.1667541262\right) \cdot z + t\right) \cdot z + a\right) \cdot z + b\right)}{\left(\left(\left(z + 15.234687407\right) \cdot z + 31.4690115749\right) \cdot z + 11.9400905721\right) \cdot z + 0.607771387771}
\end{array}
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1
(+
0.607771387771
(*
z
(+
11.9400905721
(* z (+ 31.4690115749 (* z (+ z 15.234687407))))))))
(t_2
(+
b
(*
z
(+ a (* z (+ (* z (+ (* z 3.13060547623) 11.1667541262)) t)))))))
(if (<= (/ (* y t_2) t_1) INFINITY)
(+ x (/ y (/ t_1 t_2)))
(+
(* -36.52704169880642 (/ y z))
(+ x (* y (+ 3.13060547623 (/ (+ t 457.9610022158428) (* z z)))))))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = 0.607771387771 + (z * (11.9400905721 + (z * (31.4690115749 + (z * (z + 15.234687407))))));
double t_2 = b + (z * (a + (z * ((z * ((z * 3.13060547623) + 11.1667541262)) + t))));
double tmp;
if (((y * t_2) / t_1) <= ((double) INFINITY)) {
tmp = x + (y / (t_1 / t_2));
} else {
tmp = (-36.52704169880642 * (y / z)) + (x + (y * (3.13060547623 + ((t + 457.9610022158428) / (z * z)))));
}
return tmp;
}
public static double code(double x, double y, double z, double t, double a, double b) {
double t_1 = 0.607771387771 + (z * (11.9400905721 + (z * (31.4690115749 + (z * (z + 15.234687407))))));
double t_2 = b + (z * (a + (z * ((z * ((z * 3.13060547623) + 11.1667541262)) + t))));
double tmp;
if (((y * t_2) / t_1) <= Double.POSITIVE_INFINITY) {
tmp = x + (y / (t_1 / t_2));
} else {
tmp = (-36.52704169880642 * (y / z)) + (x + (y * (3.13060547623 + ((t + 457.9610022158428) / (z * z)))));
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = 0.607771387771 + (z * (11.9400905721 + (z * (31.4690115749 + (z * (z + 15.234687407)))))) t_2 = b + (z * (a + (z * ((z * ((z * 3.13060547623) + 11.1667541262)) + t)))) tmp = 0 if ((y * t_2) / t_1) <= math.inf: tmp = x + (y / (t_1 / t_2)) else: tmp = (-36.52704169880642 * (y / z)) + (x + (y * (3.13060547623 + ((t + 457.9610022158428) / (z * z))))) return tmp
function code(x, y, z, t, a, b) t_1 = Float64(0.607771387771 + Float64(z * Float64(11.9400905721 + Float64(z * Float64(31.4690115749 + Float64(z * Float64(z + 15.234687407))))))) t_2 = Float64(b + Float64(z * Float64(a + Float64(z * Float64(Float64(z * Float64(Float64(z * 3.13060547623) + 11.1667541262)) + t))))) tmp = 0.0 if (Float64(Float64(y * t_2) / t_1) <= Inf) tmp = Float64(x + Float64(y / Float64(t_1 / t_2))); else tmp = Float64(Float64(-36.52704169880642 * Float64(y / z)) + Float64(x + Float64(y * Float64(3.13060547623 + Float64(Float64(t + 457.9610022158428) / Float64(z * z)))))); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) t_1 = 0.607771387771 + (z * (11.9400905721 + (z * (31.4690115749 + (z * (z + 15.234687407)))))); t_2 = b + (z * (a + (z * ((z * ((z * 3.13060547623) + 11.1667541262)) + t)))); tmp = 0.0; if (((y * t_2) / t_1) <= Inf) tmp = x + (y / (t_1 / t_2)); else tmp = (-36.52704169880642 * (y / z)) + (x + (y * (3.13060547623 + ((t + 457.9610022158428) / (z * z))))); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(0.607771387771 + N[(z * N[(11.9400905721 + N[(z * N[(31.4690115749 + N[(z * N[(z + 15.234687407), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(b + N[(z * N[(a + N[(z * N[(N[(z * N[(N[(z * 3.13060547623), $MachinePrecision] + 11.1667541262), $MachinePrecision]), $MachinePrecision] + t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[N[(N[(y * t$95$2), $MachinePrecision] / t$95$1), $MachinePrecision], Infinity], N[(x + N[(y / N[(t$95$1 / t$95$2), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(-36.52704169880642 * N[(y / z), $MachinePrecision]), $MachinePrecision] + N[(x + N[(y * N[(3.13060547623 + N[(N[(t + 457.9610022158428), $MachinePrecision] / N[(z * z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := 0.607771387771 + z \cdot \left(11.9400905721 + z \cdot \left(31.4690115749 + z \cdot \left(z + 15.234687407\right)\right)\right)\\
t_2 := b + z \cdot \left(a + z \cdot \left(z \cdot \left(z \cdot 3.13060547623 + 11.1667541262\right) + t\right)\right)\\
\mathbf{if}\;\frac{y \cdot t\_2}{t\_1} \leq \infty:\\
\;\;\;\;x + \frac{y}{\frac{t\_1}{t\_2}}\\
\mathbf{else}:\\
\;\;\;\;-36.52704169880642 \cdot \frac{y}{z} + \left(x + y \cdot \left(3.13060547623 + \frac{t + 457.9610022158428}{z \cdot z}\right)\right)\\
\end{array}
\end{array}
if (/.f64 (*.f64 y (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 z #s(literal 313060547623/100000000000 binary64)) #s(literal 55833770631/5000000000 binary64)) z) t) z) a) z) b)) (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 z #s(literal 15234687407/1000000000 binary64)) z) #s(literal 314690115749/10000000000 binary64)) z) #s(literal 119400905721/10000000000 binary64)) z) #s(literal 607771387771/1000000000000 binary64))) < +inf.0Initial program 97.2%
associate-/l*N/A
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
Applied egg-rr99.0%
if +inf.0 < (/.f64 (*.f64 y (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 z #s(literal 313060547623/100000000000 binary64)) #s(literal 55833770631/5000000000 binary64)) z) t) z) a) z) b)) (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 z #s(literal 15234687407/1000000000 binary64)) z) #s(literal 314690115749/10000000000 binary64)) z) #s(literal 119400905721/10000000000 binary64)) z) #s(literal 607771387771/1000000000000 binary64))) Initial program 0.0%
Taylor expanded in z around inf
Simplified98.9%
Taylor expanded in z around 0
/-lowering-/.f64N/A
associate-+r+N/A
associate-+r+N/A
+-lowering-+.f64N/A
associate-+r+N/A
distribute-rgt-outN/A
*-commutativeN/A
distribute-lft-outN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
metadata-evalN/A
*-lowering-*.f64N/A
distribute-rgt-outN/A
*-lowering-*.f64N/A
metadata-evalN/A
unpow2N/A
*-lowering-*.f6463.7%
Simplified63.7%
Taylor expanded in x around 0
+-commutativeN/A
associate-+l+N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
associate-/l*N/A
distribute-lft-outN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
unpow2N/A
*-lowering-*.f6499.9%
Simplified99.9%
Final simplification99.4%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1
(+
(* -36.52704169880642 (/ y z))
(+ x (* y (+ 3.13060547623 (/ (+ t 457.9610022158428) (* z z))))))))
(if (<= z -5.5e+62)
t_1
(if (<= z 3.35e+43)
(+
x
(*
(+
b
(* z (+ a (* z (+ (* z (+ (* z 3.13060547623) 11.1667541262)) t)))))
(/
y
(+
0.607771387771
(*
z
(+
11.9400905721
(* z (+ 31.4690115749 (* z (+ z 15.234687407))))))))))
t_1))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = (-36.52704169880642 * (y / z)) + (x + (y * (3.13060547623 + ((t + 457.9610022158428) / (z * z)))));
double tmp;
if (z <= -5.5e+62) {
tmp = t_1;
} else if (z <= 3.35e+43) {
tmp = x + ((b + (z * (a + (z * ((z * ((z * 3.13060547623) + 11.1667541262)) + t))))) * (y / (0.607771387771 + (z * (11.9400905721 + (z * (31.4690115749 + (z * (z + 15.234687407)))))))));
} else {
tmp = t_1;
}
return tmp;
}
real(8) function code(x, y, z, t, a, b)
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) :: t_1
real(8) :: tmp
t_1 = ((-36.52704169880642d0) * (y / z)) + (x + (y * (3.13060547623d0 + ((t + 457.9610022158428d0) / (z * z)))))
if (z <= (-5.5d+62)) then
tmp = t_1
else if (z <= 3.35d+43) then
tmp = x + ((b + (z * (a + (z * ((z * ((z * 3.13060547623d0) + 11.1667541262d0)) + t))))) * (y / (0.607771387771d0 + (z * (11.9400905721d0 + (z * (31.4690115749d0 + (z * (z + 15.234687407d0)))))))))
else
tmp = t_1
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b) {
double t_1 = (-36.52704169880642 * (y / z)) + (x + (y * (3.13060547623 + ((t + 457.9610022158428) / (z * z)))));
double tmp;
if (z <= -5.5e+62) {
tmp = t_1;
} else if (z <= 3.35e+43) {
tmp = x + ((b + (z * (a + (z * ((z * ((z * 3.13060547623) + 11.1667541262)) + t))))) * (y / (0.607771387771 + (z * (11.9400905721 + (z * (31.4690115749 + (z * (z + 15.234687407)))))))));
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = (-36.52704169880642 * (y / z)) + (x + (y * (3.13060547623 + ((t + 457.9610022158428) / (z * z))))) tmp = 0 if z <= -5.5e+62: tmp = t_1 elif z <= 3.35e+43: tmp = x + ((b + (z * (a + (z * ((z * ((z * 3.13060547623) + 11.1667541262)) + t))))) * (y / (0.607771387771 + (z * (11.9400905721 + (z * (31.4690115749 + (z * (z + 15.234687407))))))))) else: tmp = t_1 return tmp
function code(x, y, z, t, a, b) t_1 = Float64(Float64(-36.52704169880642 * Float64(y / z)) + Float64(x + Float64(y * Float64(3.13060547623 + Float64(Float64(t + 457.9610022158428) / Float64(z * z)))))) tmp = 0.0 if (z <= -5.5e+62) tmp = t_1; elseif (z <= 3.35e+43) tmp = Float64(x + Float64(Float64(b + Float64(z * Float64(a + Float64(z * Float64(Float64(z * Float64(Float64(z * 3.13060547623) + 11.1667541262)) + t))))) * Float64(y / Float64(0.607771387771 + Float64(z * Float64(11.9400905721 + Float64(z * Float64(31.4690115749 + Float64(z * Float64(z + 15.234687407)))))))))); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a, b) t_1 = (-36.52704169880642 * (y / z)) + (x + (y * (3.13060547623 + ((t + 457.9610022158428) / (z * z))))); tmp = 0.0; if (z <= -5.5e+62) tmp = t_1; elseif (z <= 3.35e+43) tmp = x + ((b + (z * (a + (z * ((z * ((z * 3.13060547623) + 11.1667541262)) + t))))) * (y / (0.607771387771 + (z * (11.9400905721 + (z * (31.4690115749 + (z * (z + 15.234687407))))))))); else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(N[(-36.52704169880642 * N[(y / z), $MachinePrecision]), $MachinePrecision] + N[(x + N[(y * N[(3.13060547623 + N[(N[(t + 457.9610022158428), $MachinePrecision] / N[(z * z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[z, -5.5e+62], t$95$1, If[LessEqual[z, 3.35e+43], N[(x + N[(N[(b + N[(z * N[(a + N[(z * N[(N[(z * N[(N[(z * 3.13060547623), $MachinePrecision] + 11.1667541262), $MachinePrecision]), $MachinePrecision] + t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(y / N[(0.607771387771 + N[(z * N[(11.9400905721 + N[(z * N[(31.4690115749 + N[(z * N[(z + 15.234687407), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := -36.52704169880642 \cdot \frac{y}{z} + \left(x + y \cdot \left(3.13060547623 + \frac{t + 457.9610022158428}{z \cdot z}\right)\right)\\
\mathbf{if}\;z \leq -5.5 \cdot 10^{+62}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;z \leq 3.35 \cdot 10^{+43}:\\
\;\;\;\;x + \left(b + z \cdot \left(a + z \cdot \left(z \cdot \left(z \cdot 3.13060547623 + 11.1667541262\right) + t\right)\right)\right) \cdot \frac{y}{0.607771387771 + z \cdot \left(11.9400905721 + z \cdot \left(31.4690115749 + z \cdot \left(z + 15.234687407\right)\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if z < -5.4999999999999997e62 or 3.34999999999999988e43 < z Initial program 4.6%
Taylor expanded in z around inf
Simplified98.2%
Taylor expanded in z around 0
/-lowering-/.f64N/A
associate-+r+N/A
associate-+r+N/A
+-lowering-+.f64N/A
associate-+r+N/A
distribute-rgt-outN/A
*-commutativeN/A
distribute-lft-outN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
metadata-evalN/A
*-lowering-*.f64N/A
distribute-rgt-outN/A
*-lowering-*.f64N/A
metadata-evalN/A
unpow2N/A
*-lowering-*.f6465.1%
Simplified65.1%
Taylor expanded in x around 0
+-commutativeN/A
associate-+l+N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
associate-/l*N/A
distribute-lft-outN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
unpow2N/A
*-lowering-*.f6499.1%
Simplified99.1%
if -5.4999999999999997e62 < z < 3.34999999999999988e43Initial program 98.4%
*-commutativeN/A
associate-/l*N/A
*-lowering-*.f64N/A
Applied egg-rr99.1%
Final simplification99.1%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1
(+
(* -36.52704169880642 (/ y z))
(+ x (* y (+ 3.13060547623 (/ (+ t 457.9610022158428) (* z z))))))))
(if (<= z -6.5e+22)
t_1
(if (<= z 10000.0)
(+
x
(/
(*
y
(+
b
(*
z
(+ a (* z (+ (* z (+ (* z 3.13060547623) 11.1667541262)) t))))))
(+ 0.607771387771 (* z (+ 11.9400905721 (* z 31.4690115749))))))
t_1))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = (-36.52704169880642 * (y / z)) + (x + (y * (3.13060547623 + ((t + 457.9610022158428) / (z * z)))));
double tmp;
if (z <= -6.5e+22) {
tmp = t_1;
} else if (z <= 10000.0) {
tmp = x + ((y * (b + (z * (a + (z * ((z * ((z * 3.13060547623) + 11.1667541262)) + t)))))) / (0.607771387771 + (z * (11.9400905721 + (z * 31.4690115749)))));
} else {
tmp = t_1;
}
return tmp;
}
real(8) function code(x, y, z, t, a, b)
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) :: t_1
real(8) :: tmp
t_1 = ((-36.52704169880642d0) * (y / z)) + (x + (y * (3.13060547623d0 + ((t + 457.9610022158428d0) / (z * z)))))
if (z <= (-6.5d+22)) then
tmp = t_1
else if (z <= 10000.0d0) then
tmp = x + ((y * (b + (z * (a + (z * ((z * ((z * 3.13060547623d0) + 11.1667541262d0)) + t)))))) / (0.607771387771d0 + (z * (11.9400905721d0 + (z * 31.4690115749d0)))))
else
tmp = t_1
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b) {
double t_1 = (-36.52704169880642 * (y / z)) + (x + (y * (3.13060547623 + ((t + 457.9610022158428) / (z * z)))));
double tmp;
if (z <= -6.5e+22) {
tmp = t_1;
} else if (z <= 10000.0) {
tmp = x + ((y * (b + (z * (a + (z * ((z * ((z * 3.13060547623) + 11.1667541262)) + t)))))) / (0.607771387771 + (z * (11.9400905721 + (z * 31.4690115749)))));
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = (-36.52704169880642 * (y / z)) + (x + (y * (3.13060547623 + ((t + 457.9610022158428) / (z * z))))) tmp = 0 if z <= -6.5e+22: tmp = t_1 elif z <= 10000.0: tmp = x + ((y * (b + (z * (a + (z * ((z * ((z * 3.13060547623) + 11.1667541262)) + t)))))) / (0.607771387771 + (z * (11.9400905721 + (z * 31.4690115749))))) else: tmp = t_1 return tmp
function code(x, y, z, t, a, b) t_1 = Float64(Float64(-36.52704169880642 * Float64(y / z)) + Float64(x + Float64(y * Float64(3.13060547623 + Float64(Float64(t + 457.9610022158428) / Float64(z * z)))))) tmp = 0.0 if (z <= -6.5e+22) tmp = t_1; elseif (z <= 10000.0) tmp = Float64(x + Float64(Float64(y * Float64(b + Float64(z * Float64(a + Float64(z * Float64(Float64(z * Float64(Float64(z * 3.13060547623) + 11.1667541262)) + t)))))) / Float64(0.607771387771 + Float64(z * Float64(11.9400905721 + Float64(z * 31.4690115749)))))); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a, b) t_1 = (-36.52704169880642 * (y / z)) + (x + (y * (3.13060547623 + ((t + 457.9610022158428) / (z * z))))); tmp = 0.0; if (z <= -6.5e+22) tmp = t_1; elseif (z <= 10000.0) tmp = x + ((y * (b + (z * (a + (z * ((z * ((z * 3.13060547623) + 11.1667541262)) + t)))))) / (0.607771387771 + (z * (11.9400905721 + (z * 31.4690115749))))); else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(N[(-36.52704169880642 * N[(y / z), $MachinePrecision]), $MachinePrecision] + N[(x + N[(y * N[(3.13060547623 + N[(N[(t + 457.9610022158428), $MachinePrecision] / N[(z * z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[z, -6.5e+22], t$95$1, If[LessEqual[z, 10000.0], N[(x + N[(N[(y * N[(b + N[(z * N[(a + N[(z * N[(N[(z * N[(N[(z * 3.13060547623), $MachinePrecision] + 11.1667541262), $MachinePrecision]), $MachinePrecision] + t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(0.607771387771 + N[(z * N[(11.9400905721 + N[(z * 31.4690115749), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := -36.52704169880642 \cdot \frac{y}{z} + \left(x + y \cdot \left(3.13060547623 + \frac{t + 457.9610022158428}{z \cdot z}\right)\right)\\
\mathbf{if}\;z \leq -6.5 \cdot 10^{+22}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;z \leq 10000:\\
\;\;\;\;x + \frac{y \cdot \left(b + z \cdot \left(a + z \cdot \left(z \cdot \left(z \cdot 3.13060547623 + 11.1667541262\right) + t\right)\right)\right)}{0.607771387771 + z \cdot \left(11.9400905721 + z \cdot 31.4690115749\right)}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if z < -6.49999999999999979e22 or 1e4 < z Initial program 15.8%
Taylor expanded in z around inf
Simplified95.3%
Taylor expanded in z around 0
/-lowering-/.f64N/A
associate-+r+N/A
associate-+r+N/A
+-lowering-+.f64N/A
associate-+r+N/A
distribute-rgt-outN/A
*-commutativeN/A
distribute-lft-outN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
metadata-evalN/A
*-lowering-*.f64N/A
distribute-rgt-outN/A
*-lowering-*.f64N/A
metadata-evalN/A
unpow2N/A
*-lowering-*.f6465.8%
Simplified65.8%
Taylor expanded in x around 0
+-commutativeN/A
associate-+l+N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
associate-/l*N/A
distribute-lft-outN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
unpow2N/A
*-lowering-*.f6496.1%
Simplified96.1%
if -6.49999999999999979e22 < z < 1e4Initial program 99.7%
Taylor expanded in z around 0
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6496.3%
Simplified96.3%
Final simplification96.2%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1
(+
(* -36.52704169880642 (/ y z))
(+ x (* y (+ 3.13060547623 (/ (+ t 457.9610022158428) (* z z))))))))
(if (<= z -0.05)
t_1
(if (<= z 80000.0)
(+
x
(*
(+
b
(* z (+ a (* z (+ (* z (+ (* z 3.13060547623) 11.1667541262)) t)))))
(* y (/ 1.0 (+ 0.607771387771 (* z 11.9400905721))))))
t_1))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = (-36.52704169880642 * (y / z)) + (x + (y * (3.13060547623 + ((t + 457.9610022158428) / (z * z)))));
double tmp;
if (z <= -0.05) {
tmp = t_1;
} else if (z <= 80000.0) {
tmp = x + ((b + (z * (a + (z * ((z * ((z * 3.13060547623) + 11.1667541262)) + t))))) * (y * (1.0 / (0.607771387771 + (z * 11.9400905721)))));
} else {
tmp = t_1;
}
return tmp;
}
real(8) function code(x, y, z, t, a, b)
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) :: t_1
real(8) :: tmp
t_1 = ((-36.52704169880642d0) * (y / z)) + (x + (y * (3.13060547623d0 + ((t + 457.9610022158428d0) / (z * z)))))
if (z <= (-0.05d0)) then
tmp = t_1
else if (z <= 80000.0d0) then
tmp = x + ((b + (z * (a + (z * ((z * ((z * 3.13060547623d0) + 11.1667541262d0)) + t))))) * (y * (1.0d0 / (0.607771387771d0 + (z * 11.9400905721d0)))))
else
tmp = t_1
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b) {
double t_1 = (-36.52704169880642 * (y / z)) + (x + (y * (3.13060547623 + ((t + 457.9610022158428) / (z * z)))));
double tmp;
if (z <= -0.05) {
tmp = t_1;
} else if (z <= 80000.0) {
tmp = x + ((b + (z * (a + (z * ((z * ((z * 3.13060547623) + 11.1667541262)) + t))))) * (y * (1.0 / (0.607771387771 + (z * 11.9400905721)))));
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = (-36.52704169880642 * (y / z)) + (x + (y * (3.13060547623 + ((t + 457.9610022158428) / (z * z))))) tmp = 0 if z <= -0.05: tmp = t_1 elif z <= 80000.0: tmp = x + ((b + (z * (a + (z * ((z * ((z * 3.13060547623) + 11.1667541262)) + t))))) * (y * (1.0 / (0.607771387771 + (z * 11.9400905721))))) else: tmp = t_1 return tmp
function code(x, y, z, t, a, b) t_1 = Float64(Float64(-36.52704169880642 * Float64(y / z)) + Float64(x + Float64(y * Float64(3.13060547623 + Float64(Float64(t + 457.9610022158428) / Float64(z * z)))))) tmp = 0.0 if (z <= -0.05) tmp = t_1; elseif (z <= 80000.0) tmp = Float64(x + Float64(Float64(b + Float64(z * Float64(a + Float64(z * Float64(Float64(z * Float64(Float64(z * 3.13060547623) + 11.1667541262)) + t))))) * Float64(y * Float64(1.0 / Float64(0.607771387771 + Float64(z * 11.9400905721)))))); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a, b) t_1 = (-36.52704169880642 * (y / z)) + (x + (y * (3.13060547623 + ((t + 457.9610022158428) / (z * z))))); tmp = 0.0; if (z <= -0.05) tmp = t_1; elseif (z <= 80000.0) tmp = x + ((b + (z * (a + (z * ((z * ((z * 3.13060547623) + 11.1667541262)) + t))))) * (y * (1.0 / (0.607771387771 + (z * 11.9400905721))))); else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(N[(-36.52704169880642 * N[(y / z), $MachinePrecision]), $MachinePrecision] + N[(x + N[(y * N[(3.13060547623 + N[(N[(t + 457.9610022158428), $MachinePrecision] / N[(z * z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[z, -0.05], t$95$1, If[LessEqual[z, 80000.0], N[(x + N[(N[(b + N[(z * N[(a + N[(z * N[(N[(z * N[(N[(z * 3.13060547623), $MachinePrecision] + 11.1667541262), $MachinePrecision]), $MachinePrecision] + t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(y * N[(1.0 / N[(0.607771387771 + N[(z * 11.9400905721), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := -36.52704169880642 \cdot \frac{y}{z} + \left(x + y \cdot \left(3.13060547623 + \frac{t + 457.9610022158428}{z \cdot z}\right)\right)\\
\mathbf{if}\;z \leq -0.05:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;z \leq 80000:\\
\;\;\;\;x + \left(b + z \cdot \left(a + z \cdot \left(z \cdot \left(z \cdot 3.13060547623 + 11.1667541262\right) + t\right)\right)\right) \cdot \left(y \cdot \frac{1}{0.607771387771 + z \cdot 11.9400905721}\right)\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if z < -0.050000000000000003 or 8e4 < z Initial program 19.5%
Taylor expanded in z around inf
Simplified92.7%
Taylor expanded in z around 0
/-lowering-/.f64N/A
associate-+r+N/A
associate-+r+N/A
+-lowering-+.f64N/A
associate-+r+N/A
distribute-rgt-outN/A
*-commutativeN/A
distribute-lft-outN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
metadata-evalN/A
*-lowering-*.f64N/A
distribute-rgt-outN/A
*-lowering-*.f64N/A
metadata-evalN/A
unpow2N/A
*-lowering-*.f6464.6%
Simplified64.6%
Taylor expanded in x around 0
+-commutativeN/A
associate-+l+N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
associate-/l*N/A
distribute-lft-outN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
unpow2N/A
*-lowering-*.f6493.5%
Simplified93.5%
if -0.050000000000000003 < z < 8e4Initial program 99.7%
Taylor expanded in z around 0
*-commutativeN/A
*-lowering-*.f6498.8%
Simplified98.8%
div-invN/A
*-commutativeN/A
associate-*l*N/A
*-lowering-*.f64N/A
Applied egg-rr98.9%
Final simplification96.0%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1
(+
(* -36.52704169880642 (/ y z))
(+ x (* y (+ 3.13060547623 (/ (+ t 457.9610022158428) (* z z))))))))
(if (<= z -0.05)
t_1
(if (<= z 66000.0)
(+
x
(*
(+
b
(* z (+ a (* z (+ (* z (+ (* z 3.13060547623) 11.1667541262)) t)))))
(/ y (+ 0.607771387771 (* z 11.9400905721)))))
t_1))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = (-36.52704169880642 * (y / z)) + (x + (y * (3.13060547623 + ((t + 457.9610022158428) / (z * z)))));
double tmp;
if (z <= -0.05) {
tmp = t_1;
} else if (z <= 66000.0) {
tmp = x + ((b + (z * (a + (z * ((z * ((z * 3.13060547623) + 11.1667541262)) + t))))) * (y / (0.607771387771 + (z * 11.9400905721))));
} else {
tmp = t_1;
}
return tmp;
}
real(8) function code(x, y, z, t, a, b)
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) :: t_1
real(8) :: tmp
t_1 = ((-36.52704169880642d0) * (y / z)) + (x + (y * (3.13060547623d0 + ((t + 457.9610022158428d0) / (z * z)))))
if (z <= (-0.05d0)) then
tmp = t_1
else if (z <= 66000.0d0) then
tmp = x + ((b + (z * (a + (z * ((z * ((z * 3.13060547623d0) + 11.1667541262d0)) + t))))) * (y / (0.607771387771d0 + (z * 11.9400905721d0))))
else
tmp = t_1
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b) {
double t_1 = (-36.52704169880642 * (y / z)) + (x + (y * (3.13060547623 + ((t + 457.9610022158428) / (z * z)))));
double tmp;
if (z <= -0.05) {
tmp = t_1;
} else if (z <= 66000.0) {
tmp = x + ((b + (z * (a + (z * ((z * ((z * 3.13060547623) + 11.1667541262)) + t))))) * (y / (0.607771387771 + (z * 11.9400905721))));
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = (-36.52704169880642 * (y / z)) + (x + (y * (3.13060547623 + ((t + 457.9610022158428) / (z * z))))) tmp = 0 if z <= -0.05: tmp = t_1 elif z <= 66000.0: tmp = x + ((b + (z * (a + (z * ((z * ((z * 3.13060547623) + 11.1667541262)) + t))))) * (y / (0.607771387771 + (z * 11.9400905721)))) else: tmp = t_1 return tmp
function code(x, y, z, t, a, b) t_1 = Float64(Float64(-36.52704169880642 * Float64(y / z)) + Float64(x + Float64(y * Float64(3.13060547623 + Float64(Float64(t + 457.9610022158428) / Float64(z * z)))))) tmp = 0.0 if (z <= -0.05) tmp = t_1; elseif (z <= 66000.0) tmp = Float64(x + Float64(Float64(b + Float64(z * Float64(a + Float64(z * Float64(Float64(z * Float64(Float64(z * 3.13060547623) + 11.1667541262)) + t))))) * Float64(y / Float64(0.607771387771 + Float64(z * 11.9400905721))))); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a, b) t_1 = (-36.52704169880642 * (y / z)) + (x + (y * (3.13060547623 + ((t + 457.9610022158428) / (z * z))))); tmp = 0.0; if (z <= -0.05) tmp = t_1; elseif (z <= 66000.0) tmp = x + ((b + (z * (a + (z * ((z * ((z * 3.13060547623) + 11.1667541262)) + t))))) * (y / (0.607771387771 + (z * 11.9400905721)))); else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(N[(-36.52704169880642 * N[(y / z), $MachinePrecision]), $MachinePrecision] + N[(x + N[(y * N[(3.13060547623 + N[(N[(t + 457.9610022158428), $MachinePrecision] / N[(z * z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[z, -0.05], t$95$1, If[LessEqual[z, 66000.0], N[(x + N[(N[(b + N[(z * N[(a + N[(z * N[(N[(z * N[(N[(z * 3.13060547623), $MachinePrecision] + 11.1667541262), $MachinePrecision]), $MachinePrecision] + t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(y / N[(0.607771387771 + N[(z * 11.9400905721), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := -36.52704169880642 \cdot \frac{y}{z} + \left(x + y \cdot \left(3.13060547623 + \frac{t + 457.9610022158428}{z \cdot z}\right)\right)\\
\mathbf{if}\;z \leq -0.05:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;z \leq 66000:\\
\;\;\;\;x + \left(b + z \cdot \left(a + z \cdot \left(z \cdot \left(z \cdot 3.13060547623 + 11.1667541262\right) + t\right)\right)\right) \cdot \frac{y}{0.607771387771 + z \cdot 11.9400905721}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if z < -0.050000000000000003 or 66000 < z Initial program 19.5%
Taylor expanded in z around inf
Simplified92.7%
Taylor expanded in z around 0
/-lowering-/.f64N/A
associate-+r+N/A
associate-+r+N/A
+-lowering-+.f64N/A
associate-+r+N/A
distribute-rgt-outN/A
*-commutativeN/A
distribute-lft-outN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
metadata-evalN/A
*-lowering-*.f64N/A
distribute-rgt-outN/A
*-lowering-*.f64N/A
metadata-evalN/A
unpow2N/A
*-lowering-*.f6464.6%
Simplified64.6%
Taylor expanded in x around 0
+-commutativeN/A
associate-+l+N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
associate-/l*N/A
distribute-lft-outN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
unpow2N/A
*-lowering-*.f6493.5%
Simplified93.5%
if -0.050000000000000003 < z < 66000Initial program 99.7%
Taylor expanded in z around 0
*-commutativeN/A
*-lowering-*.f6498.8%
Simplified98.8%
*-commutativeN/A
associate-/l*N/A
*-lowering-*.f64N/A
Applied egg-rr98.8%
Final simplification96.0%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1
(+
(* -36.52704169880642 (/ y z))
(+ x (* y (+ 3.13060547623 (/ (+ t 457.9610022158428) (* z z))))))))
(if (<= z -0.05)
t_1
(if (<= z 72000.0)
(+
x
(/
(* y (+ b (* z (+ a (* z (+ t (* z 11.1667541262)))))))
(+ 0.607771387771 (* z 11.9400905721))))
t_1))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = (-36.52704169880642 * (y / z)) + (x + (y * (3.13060547623 + ((t + 457.9610022158428) / (z * z)))));
double tmp;
if (z <= -0.05) {
tmp = t_1;
} else if (z <= 72000.0) {
tmp = x + ((y * (b + (z * (a + (z * (t + (z * 11.1667541262))))))) / (0.607771387771 + (z * 11.9400905721)));
} else {
tmp = t_1;
}
return tmp;
}
real(8) function code(x, y, z, t, a, b)
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) :: t_1
real(8) :: tmp
t_1 = ((-36.52704169880642d0) * (y / z)) + (x + (y * (3.13060547623d0 + ((t + 457.9610022158428d0) / (z * z)))))
if (z <= (-0.05d0)) then
tmp = t_1
else if (z <= 72000.0d0) then
tmp = x + ((y * (b + (z * (a + (z * (t + (z * 11.1667541262d0))))))) / (0.607771387771d0 + (z * 11.9400905721d0)))
else
tmp = t_1
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b) {
double t_1 = (-36.52704169880642 * (y / z)) + (x + (y * (3.13060547623 + ((t + 457.9610022158428) / (z * z)))));
double tmp;
if (z <= -0.05) {
tmp = t_1;
} else if (z <= 72000.0) {
tmp = x + ((y * (b + (z * (a + (z * (t + (z * 11.1667541262))))))) / (0.607771387771 + (z * 11.9400905721)));
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = (-36.52704169880642 * (y / z)) + (x + (y * (3.13060547623 + ((t + 457.9610022158428) / (z * z))))) tmp = 0 if z <= -0.05: tmp = t_1 elif z <= 72000.0: tmp = x + ((y * (b + (z * (a + (z * (t + (z * 11.1667541262))))))) / (0.607771387771 + (z * 11.9400905721))) else: tmp = t_1 return tmp
function code(x, y, z, t, a, b) t_1 = Float64(Float64(-36.52704169880642 * Float64(y / z)) + Float64(x + Float64(y * Float64(3.13060547623 + Float64(Float64(t + 457.9610022158428) / Float64(z * z)))))) tmp = 0.0 if (z <= -0.05) tmp = t_1; elseif (z <= 72000.0) tmp = Float64(x + Float64(Float64(y * Float64(b + Float64(z * Float64(a + Float64(z * Float64(t + Float64(z * 11.1667541262))))))) / Float64(0.607771387771 + Float64(z * 11.9400905721)))); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a, b) t_1 = (-36.52704169880642 * (y / z)) + (x + (y * (3.13060547623 + ((t + 457.9610022158428) / (z * z))))); tmp = 0.0; if (z <= -0.05) tmp = t_1; elseif (z <= 72000.0) tmp = x + ((y * (b + (z * (a + (z * (t + (z * 11.1667541262))))))) / (0.607771387771 + (z * 11.9400905721))); else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(N[(-36.52704169880642 * N[(y / z), $MachinePrecision]), $MachinePrecision] + N[(x + N[(y * N[(3.13060547623 + N[(N[(t + 457.9610022158428), $MachinePrecision] / N[(z * z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[z, -0.05], t$95$1, If[LessEqual[z, 72000.0], N[(x + N[(N[(y * N[(b + N[(z * N[(a + N[(z * N[(t + N[(z * 11.1667541262), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(0.607771387771 + N[(z * 11.9400905721), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := -36.52704169880642 \cdot \frac{y}{z} + \left(x + y \cdot \left(3.13060547623 + \frac{t + 457.9610022158428}{z \cdot z}\right)\right)\\
\mathbf{if}\;z \leq -0.05:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;z \leq 72000:\\
\;\;\;\;x + \frac{y \cdot \left(b + z \cdot \left(a + z \cdot \left(t + z \cdot 11.1667541262\right)\right)\right)}{0.607771387771 + z \cdot 11.9400905721}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if z < -0.050000000000000003 or 72000 < z Initial program 19.5%
Taylor expanded in z around inf
Simplified92.7%
Taylor expanded in z around 0
/-lowering-/.f64N/A
associate-+r+N/A
associate-+r+N/A
+-lowering-+.f64N/A
associate-+r+N/A
distribute-rgt-outN/A
*-commutativeN/A
distribute-lft-outN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
metadata-evalN/A
*-lowering-*.f64N/A
distribute-rgt-outN/A
*-lowering-*.f64N/A
metadata-evalN/A
unpow2N/A
*-lowering-*.f6464.6%
Simplified64.6%
Taylor expanded in x around 0
+-commutativeN/A
associate-+l+N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
associate-/l*N/A
distribute-lft-outN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
unpow2N/A
*-lowering-*.f6493.5%
Simplified93.5%
if -0.050000000000000003 < z < 72000Initial program 99.7%
Taylor expanded in z around 0
*-commutativeN/A
*-lowering-*.f6498.8%
Simplified98.8%
Taylor expanded in z around 0
*-commutativeN/A
*-lowering-*.f6498.5%
Simplified98.5%
Final simplification95.9%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1
(+
(* -36.52704169880642 (/ y z))
(+ x (* y (+ 3.13060547623 (/ (+ t 457.9610022158428) (* z z))))))))
(if (<= z -0.05)
t_1
(if (<= z 55000.0)
(+
x
(/
(* y (+ b (* z (+ a (* z t)))))
(+ 0.607771387771 (* z 11.9400905721))))
t_1))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = (-36.52704169880642 * (y / z)) + (x + (y * (3.13060547623 + ((t + 457.9610022158428) / (z * z)))));
double tmp;
if (z <= -0.05) {
tmp = t_1;
} else if (z <= 55000.0) {
tmp = x + ((y * (b + (z * (a + (z * t))))) / (0.607771387771 + (z * 11.9400905721)));
} else {
tmp = t_1;
}
return tmp;
}
real(8) function code(x, y, z, t, a, b)
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) :: t_1
real(8) :: tmp
t_1 = ((-36.52704169880642d0) * (y / z)) + (x + (y * (3.13060547623d0 + ((t + 457.9610022158428d0) / (z * z)))))
if (z <= (-0.05d0)) then
tmp = t_1
else if (z <= 55000.0d0) then
tmp = x + ((y * (b + (z * (a + (z * t))))) / (0.607771387771d0 + (z * 11.9400905721d0)))
else
tmp = t_1
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b) {
double t_1 = (-36.52704169880642 * (y / z)) + (x + (y * (3.13060547623 + ((t + 457.9610022158428) / (z * z)))));
double tmp;
if (z <= -0.05) {
tmp = t_1;
} else if (z <= 55000.0) {
tmp = x + ((y * (b + (z * (a + (z * t))))) / (0.607771387771 + (z * 11.9400905721)));
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = (-36.52704169880642 * (y / z)) + (x + (y * (3.13060547623 + ((t + 457.9610022158428) / (z * z))))) tmp = 0 if z <= -0.05: tmp = t_1 elif z <= 55000.0: tmp = x + ((y * (b + (z * (a + (z * t))))) / (0.607771387771 + (z * 11.9400905721))) else: tmp = t_1 return tmp
function code(x, y, z, t, a, b) t_1 = Float64(Float64(-36.52704169880642 * Float64(y / z)) + Float64(x + Float64(y * Float64(3.13060547623 + Float64(Float64(t + 457.9610022158428) / Float64(z * z)))))) tmp = 0.0 if (z <= -0.05) tmp = t_1; elseif (z <= 55000.0) tmp = Float64(x + Float64(Float64(y * Float64(b + Float64(z * Float64(a + Float64(z * t))))) / Float64(0.607771387771 + Float64(z * 11.9400905721)))); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a, b) t_1 = (-36.52704169880642 * (y / z)) + (x + (y * (3.13060547623 + ((t + 457.9610022158428) / (z * z))))); tmp = 0.0; if (z <= -0.05) tmp = t_1; elseif (z <= 55000.0) tmp = x + ((y * (b + (z * (a + (z * t))))) / (0.607771387771 + (z * 11.9400905721))); else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(N[(-36.52704169880642 * N[(y / z), $MachinePrecision]), $MachinePrecision] + N[(x + N[(y * N[(3.13060547623 + N[(N[(t + 457.9610022158428), $MachinePrecision] / N[(z * z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[z, -0.05], t$95$1, If[LessEqual[z, 55000.0], N[(x + N[(N[(y * N[(b + N[(z * N[(a + N[(z * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(0.607771387771 + N[(z * 11.9400905721), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := -36.52704169880642 \cdot \frac{y}{z} + \left(x + y \cdot \left(3.13060547623 + \frac{t + 457.9610022158428}{z \cdot z}\right)\right)\\
\mathbf{if}\;z \leq -0.05:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;z \leq 55000:\\
\;\;\;\;x + \frac{y \cdot \left(b + z \cdot \left(a + z \cdot t\right)\right)}{0.607771387771 + z \cdot 11.9400905721}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if z < -0.050000000000000003 or 55000 < z Initial program 19.5%
Taylor expanded in z around inf
Simplified92.7%
Taylor expanded in z around 0
/-lowering-/.f64N/A
associate-+r+N/A
associate-+r+N/A
+-lowering-+.f64N/A
associate-+r+N/A
distribute-rgt-outN/A
*-commutativeN/A
distribute-lft-outN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
metadata-evalN/A
*-lowering-*.f64N/A
distribute-rgt-outN/A
*-lowering-*.f64N/A
metadata-evalN/A
unpow2N/A
*-lowering-*.f6464.6%
Simplified64.6%
Taylor expanded in x around 0
+-commutativeN/A
associate-+l+N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
associate-/l*N/A
distribute-lft-outN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
unpow2N/A
*-lowering-*.f6493.5%
Simplified93.5%
if -0.050000000000000003 < z < 55000Initial program 99.7%
Taylor expanded in z around 0
*-commutativeN/A
*-lowering-*.f6498.8%
Simplified98.8%
Taylor expanded in z around 0
+-commutativeN/A
distribute-lft-inN/A
*-commutativeN/A
associate-*r*N/A
unpow2N/A
*-commutativeN/A
*-commutativeN/A
+-lowering-+.f64N/A
*-commutativeN/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
*-commutativeN/A
distribute-lft-inN/A
+-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6496.6%
Simplified96.6%
Final simplification95.0%
(FPCore (x y z t a b)
:precision binary64
(if (<= z -0.05)
(+ x (+ (* y 3.13060547623) (/ (* y -36.52704169880642) z)))
(if (<= z 105000.0)
(+
x
(/
(* y (+ b (* z (+ a (* z t)))))
(+ 0.607771387771 (* z 11.9400905721))))
(+ x (* y 3.13060547623)))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (z <= -0.05) {
tmp = x + ((y * 3.13060547623) + ((y * -36.52704169880642) / z));
} else if (z <= 105000.0) {
tmp = x + ((y * (b + (z * (a + (z * t))))) / (0.607771387771 + (z * 11.9400905721)));
} else {
tmp = x + (y * 3.13060547623);
}
return tmp;
}
real(8) function code(x, y, z, t, a, b)
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) :: tmp
if (z <= (-0.05d0)) then
tmp = x + ((y * 3.13060547623d0) + ((y * (-36.52704169880642d0)) / z))
else if (z <= 105000.0d0) then
tmp = x + ((y * (b + (z * (a + (z * t))))) / (0.607771387771d0 + (z * 11.9400905721d0)))
else
tmp = x + (y * 3.13060547623d0)
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (z <= -0.05) {
tmp = x + ((y * 3.13060547623) + ((y * -36.52704169880642) / z));
} else if (z <= 105000.0) {
tmp = x + ((y * (b + (z * (a + (z * t))))) / (0.607771387771 + (z * 11.9400905721)));
} else {
tmp = x + (y * 3.13060547623);
}
return tmp;
}
def code(x, y, z, t, a, b): tmp = 0 if z <= -0.05: tmp = x + ((y * 3.13060547623) + ((y * -36.52704169880642) / z)) elif z <= 105000.0: tmp = x + ((y * (b + (z * (a + (z * t))))) / (0.607771387771 + (z * 11.9400905721))) else: tmp = x + (y * 3.13060547623) return tmp
function code(x, y, z, t, a, b) tmp = 0.0 if (z <= -0.05) tmp = Float64(x + Float64(Float64(y * 3.13060547623) + Float64(Float64(y * -36.52704169880642) / z))); elseif (z <= 105000.0) tmp = Float64(x + Float64(Float64(y * Float64(b + Float64(z * Float64(a + Float64(z * t))))) / Float64(0.607771387771 + Float64(z * 11.9400905721)))); else tmp = Float64(x + Float64(y * 3.13060547623)); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) tmp = 0.0; if (z <= -0.05) tmp = x + ((y * 3.13060547623) + ((y * -36.52704169880642) / z)); elseif (z <= 105000.0) tmp = x + ((y * (b + (z * (a + (z * t))))) / (0.607771387771 + (z * 11.9400905721))); else tmp = x + (y * 3.13060547623); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[z, -0.05], N[(x + N[(N[(y * 3.13060547623), $MachinePrecision] + N[(N[(y * -36.52704169880642), $MachinePrecision] / z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[z, 105000.0], N[(x + N[(N[(y * N[(b + N[(z * N[(a + N[(z * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(0.607771387771 + N[(z * 11.9400905721), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x + N[(y * 3.13060547623), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;z \leq -0.05:\\
\;\;\;\;x + \left(y \cdot 3.13060547623 + \frac{y \cdot -36.52704169880642}{z}\right)\\
\mathbf{elif}\;z \leq 105000:\\
\;\;\;\;x + \frac{y \cdot \left(b + z \cdot \left(a + z \cdot t\right)\right)}{0.607771387771 + z \cdot 11.9400905721}\\
\mathbf{else}:\\
\;\;\;\;x + y \cdot 3.13060547623\\
\end{array}
\end{array}
if z < -0.050000000000000003Initial program 25.0%
Taylor expanded in z around inf
associate--l+N/A
associate--l+N/A
+-commutativeN/A
distribute-rgt-out--N/A
metadata-evalN/A
metadata-evalN/A
metadata-evalN/A
times-fracN/A
*-commutativeN/A
distribute-rgt-out--N/A
mul-1-negN/A
distribute-neg-frac2N/A
mul-1-negN/A
+-lowering-+.f64N/A
Simplified85.8%
if -0.050000000000000003 < z < 105000Initial program 99.7%
Taylor expanded in z around 0
*-commutativeN/A
*-lowering-*.f6498.8%
Simplified98.8%
Taylor expanded in z around 0
+-commutativeN/A
distribute-lft-inN/A
*-commutativeN/A
associate-*r*N/A
unpow2N/A
*-commutativeN/A
*-commutativeN/A
+-lowering-+.f64N/A
*-commutativeN/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
*-commutativeN/A
distribute-lft-inN/A
+-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6496.6%
Simplified96.6%
if 105000 < z Initial program 14.4%
Taylor expanded in z around inf
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6486.6%
Simplified86.6%
Final simplification91.2%
(FPCore (x y z t a b)
:precision binary64
(if (<= z -6500.0)
(+ x (+ (* y 3.13060547623) (/ (* y -36.52704169880642) z)))
(if (<= z 1.5e-94)
(+ x (* y (* b 1.6453555072203998)))
(if (<= z 58000.0)
(+ x (/ (* y (* z a)) (+ 0.607771387771 (* z 11.9400905721))))
(+ x (* y 3.13060547623))))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (z <= -6500.0) {
tmp = x + ((y * 3.13060547623) + ((y * -36.52704169880642) / z));
} else if (z <= 1.5e-94) {
tmp = x + (y * (b * 1.6453555072203998));
} else if (z <= 58000.0) {
tmp = x + ((y * (z * a)) / (0.607771387771 + (z * 11.9400905721)));
} else {
tmp = x + (y * 3.13060547623);
}
return tmp;
}
real(8) function code(x, y, z, t, a, b)
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) :: tmp
if (z <= (-6500.0d0)) then
tmp = x + ((y * 3.13060547623d0) + ((y * (-36.52704169880642d0)) / z))
else if (z <= 1.5d-94) then
tmp = x + (y * (b * 1.6453555072203998d0))
else if (z <= 58000.0d0) then
tmp = x + ((y * (z * a)) / (0.607771387771d0 + (z * 11.9400905721d0)))
else
tmp = x + (y * 3.13060547623d0)
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (z <= -6500.0) {
tmp = x + ((y * 3.13060547623) + ((y * -36.52704169880642) / z));
} else if (z <= 1.5e-94) {
tmp = x + (y * (b * 1.6453555072203998));
} else if (z <= 58000.0) {
tmp = x + ((y * (z * a)) / (0.607771387771 + (z * 11.9400905721)));
} else {
tmp = x + (y * 3.13060547623);
}
return tmp;
}
def code(x, y, z, t, a, b): tmp = 0 if z <= -6500.0: tmp = x + ((y * 3.13060547623) + ((y * -36.52704169880642) / z)) elif z <= 1.5e-94: tmp = x + (y * (b * 1.6453555072203998)) elif z <= 58000.0: tmp = x + ((y * (z * a)) / (0.607771387771 + (z * 11.9400905721))) else: tmp = x + (y * 3.13060547623) return tmp
function code(x, y, z, t, a, b) tmp = 0.0 if (z <= -6500.0) tmp = Float64(x + Float64(Float64(y * 3.13060547623) + Float64(Float64(y * -36.52704169880642) / z))); elseif (z <= 1.5e-94) tmp = Float64(x + Float64(y * Float64(b * 1.6453555072203998))); elseif (z <= 58000.0) tmp = Float64(x + Float64(Float64(y * Float64(z * a)) / Float64(0.607771387771 + Float64(z * 11.9400905721)))); else tmp = Float64(x + Float64(y * 3.13060547623)); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) tmp = 0.0; if (z <= -6500.0) tmp = x + ((y * 3.13060547623) + ((y * -36.52704169880642) / z)); elseif (z <= 1.5e-94) tmp = x + (y * (b * 1.6453555072203998)); elseif (z <= 58000.0) tmp = x + ((y * (z * a)) / (0.607771387771 + (z * 11.9400905721))); else tmp = x + (y * 3.13060547623); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[z, -6500.0], N[(x + N[(N[(y * 3.13060547623), $MachinePrecision] + N[(N[(y * -36.52704169880642), $MachinePrecision] / z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[z, 1.5e-94], N[(x + N[(y * N[(b * 1.6453555072203998), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[z, 58000.0], N[(x + N[(N[(y * N[(z * a), $MachinePrecision]), $MachinePrecision] / N[(0.607771387771 + N[(z * 11.9400905721), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x + N[(y * 3.13060547623), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;z \leq -6500:\\
\;\;\;\;x + \left(y \cdot 3.13060547623 + \frac{y \cdot -36.52704169880642}{z}\right)\\
\mathbf{elif}\;z \leq 1.5 \cdot 10^{-94}:\\
\;\;\;\;x + y \cdot \left(b \cdot 1.6453555072203998\right)\\
\mathbf{elif}\;z \leq 58000:\\
\;\;\;\;x + \frac{y \cdot \left(z \cdot a\right)}{0.607771387771 + z \cdot 11.9400905721}\\
\mathbf{else}:\\
\;\;\;\;x + y \cdot 3.13060547623\\
\end{array}
\end{array}
if z < -6500Initial program 23.9%
Taylor expanded in z around inf
associate--l+N/A
associate--l+N/A
+-commutativeN/A
distribute-rgt-out--N/A
metadata-evalN/A
metadata-evalN/A
metadata-evalN/A
times-fracN/A
*-commutativeN/A
distribute-rgt-out--N/A
mul-1-negN/A
distribute-neg-frac2N/A
mul-1-negN/A
+-lowering-+.f64N/A
Simplified87.1%
if -6500 < z < 1.5000000000000001e-94Initial program 99.7%
Taylor expanded in z around 0
+-lowering-+.f64N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f6483.3%
Simplified83.3%
if 1.5000000000000001e-94 < z < 58000Initial program 99.8%
Taylor expanded in z around 0
*-commutativeN/A
*-lowering-*.f6495.9%
Simplified95.9%
Taylor expanded in a around inf
associate-*r*N/A
*-commutativeN/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f6482.0%
Simplified82.0%
if 58000 < z Initial program 14.4%
Taylor expanded in z around inf
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6486.6%
Simplified86.6%
Final simplification85.0%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (+ x (* y 3.13060547623))))
(if (<= z -1.55e+18)
t_1
(if (<= z 7200.0)
(+
x
(*
y
(+
(* b 1.6453555072203998)
(* z (+ (* a 1.6453555072203998) (* b -32.324150453290734))))))
t_1))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = x + (y * 3.13060547623);
double tmp;
if (z <= -1.55e+18) {
tmp = t_1;
} else if (z <= 7200.0) {
tmp = x + (y * ((b * 1.6453555072203998) + (z * ((a * 1.6453555072203998) + (b * -32.324150453290734)))));
} else {
tmp = t_1;
}
return tmp;
}
real(8) function code(x, y, z, t, a, b)
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) :: t_1
real(8) :: tmp
t_1 = x + (y * 3.13060547623d0)
if (z <= (-1.55d+18)) then
tmp = t_1
else if (z <= 7200.0d0) then
tmp = x + (y * ((b * 1.6453555072203998d0) + (z * ((a * 1.6453555072203998d0) + (b * (-32.324150453290734d0))))))
else
tmp = t_1
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b) {
double t_1 = x + (y * 3.13060547623);
double tmp;
if (z <= -1.55e+18) {
tmp = t_1;
} else if (z <= 7200.0) {
tmp = x + (y * ((b * 1.6453555072203998) + (z * ((a * 1.6453555072203998) + (b * -32.324150453290734)))));
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = x + (y * 3.13060547623) tmp = 0 if z <= -1.55e+18: tmp = t_1 elif z <= 7200.0: tmp = x + (y * ((b * 1.6453555072203998) + (z * ((a * 1.6453555072203998) + (b * -32.324150453290734))))) else: tmp = t_1 return tmp
function code(x, y, z, t, a, b) t_1 = Float64(x + Float64(y * 3.13060547623)) tmp = 0.0 if (z <= -1.55e+18) tmp = t_1; elseif (z <= 7200.0) tmp = Float64(x + Float64(y * Float64(Float64(b * 1.6453555072203998) + Float64(z * Float64(Float64(a * 1.6453555072203998) + Float64(b * -32.324150453290734)))))); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a, b) t_1 = x + (y * 3.13060547623); tmp = 0.0; if (z <= -1.55e+18) tmp = t_1; elseif (z <= 7200.0) tmp = x + (y * ((b * 1.6453555072203998) + (z * ((a * 1.6453555072203998) + (b * -32.324150453290734))))); else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(x + N[(y * 3.13060547623), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[z, -1.55e+18], t$95$1, If[LessEqual[z, 7200.0], N[(x + N[(y * N[(N[(b * 1.6453555072203998), $MachinePrecision] + N[(z * N[(N[(a * 1.6453555072203998), $MachinePrecision] + N[(b * -32.324150453290734), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x + y \cdot 3.13060547623\\
\mathbf{if}\;z \leq -1.55 \cdot 10^{+18}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;z \leq 7200:\\
\;\;\;\;x + y \cdot \left(b \cdot 1.6453555072203998 + z \cdot \left(a \cdot 1.6453555072203998 + b \cdot -32.324150453290734\right)\right)\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if z < -1.55e18 or 7200 < z Initial program 17.1%
Taylor expanded in z around inf
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6488.4%
Simplified88.4%
if -1.55e18 < z < 7200Initial program 99.7%
Taylor expanded in z around 0
+-commutativeN/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
associate-*r*N/A
*-commutativeN/A
associate-*l*N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f6491.2%
Simplified91.2%
div-invN/A
distribute-lft-outN/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
Applied egg-rr91.3%
Taylor expanded in z around 0
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
distribute-rgt-neg-inN/A
*-lowering-*.f64N/A
metadata-eval90.5%
Simplified90.5%
(FPCore (x y z t a b)
:precision binary64
(if (<= z -0.05)
(+ x (+ (* y 3.13060547623) (/ (* y -36.52704169880642) z)))
(if (<= z 105000.0)
(+ x (/ (* y (+ b (* z a))) (+ 0.607771387771 (* z 11.9400905721))))
(+ x (* y 3.13060547623)))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (z <= -0.05) {
tmp = x + ((y * 3.13060547623) + ((y * -36.52704169880642) / z));
} else if (z <= 105000.0) {
tmp = x + ((y * (b + (z * a))) / (0.607771387771 + (z * 11.9400905721)));
} else {
tmp = x + (y * 3.13060547623);
}
return tmp;
}
real(8) function code(x, y, z, t, a, b)
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) :: tmp
if (z <= (-0.05d0)) then
tmp = x + ((y * 3.13060547623d0) + ((y * (-36.52704169880642d0)) / z))
else if (z <= 105000.0d0) then
tmp = x + ((y * (b + (z * a))) / (0.607771387771d0 + (z * 11.9400905721d0)))
else
tmp = x + (y * 3.13060547623d0)
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (z <= -0.05) {
tmp = x + ((y * 3.13060547623) + ((y * -36.52704169880642) / z));
} else if (z <= 105000.0) {
tmp = x + ((y * (b + (z * a))) / (0.607771387771 + (z * 11.9400905721)));
} else {
tmp = x + (y * 3.13060547623);
}
return tmp;
}
def code(x, y, z, t, a, b): tmp = 0 if z <= -0.05: tmp = x + ((y * 3.13060547623) + ((y * -36.52704169880642) / z)) elif z <= 105000.0: tmp = x + ((y * (b + (z * a))) / (0.607771387771 + (z * 11.9400905721))) else: tmp = x + (y * 3.13060547623) return tmp
function code(x, y, z, t, a, b) tmp = 0.0 if (z <= -0.05) tmp = Float64(x + Float64(Float64(y * 3.13060547623) + Float64(Float64(y * -36.52704169880642) / z))); elseif (z <= 105000.0) tmp = Float64(x + Float64(Float64(y * Float64(b + Float64(z * a))) / Float64(0.607771387771 + Float64(z * 11.9400905721)))); else tmp = Float64(x + Float64(y * 3.13060547623)); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) tmp = 0.0; if (z <= -0.05) tmp = x + ((y * 3.13060547623) + ((y * -36.52704169880642) / z)); elseif (z <= 105000.0) tmp = x + ((y * (b + (z * a))) / (0.607771387771 + (z * 11.9400905721))); else tmp = x + (y * 3.13060547623); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[z, -0.05], N[(x + N[(N[(y * 3.13060547623), $MachinePrecision] + N[(N[(y * -36.52704169880642), $MachinePrecision] / z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[z, 105000.0], N[(x + N[(N[(y * N[(b + N[(z * a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(0.607771387771 + N[(z * 11.9400905721), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x + N[(y * 3.13060547623), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;z \leq -0.05:\\
\;\;\;\;x + \left(y \cdot 3.13060547623 + \frac{y \cdot -36.52704169880642}{z}\right)\\
\mathbf{elif}\;z \leq 105000:\\
\;\;\;\;x + \frac{y \cdot \left(b + z \cdot a\right)}{0.607771387771 + z \cdot 11.9400905721}\\
\mathbf{else}:\\
\;\;\;\;x + y \cdot 3.13060547623\\
\end{array}
\end{array}
if z < -0.050000000000000003Initial program 25.0%
Taylor expanded in z around inf
associate--l+N/A
associate--l+N/A
+-commutativeN/A
distribute-rgt-out--N/A
metadata-evalN/A
metadata-evalN/A
metadata-evalN/A
times-fracN/A
*-commutativeN/A
distribute-rgt-out--N/A
mul-1-negN/A
distribute-neg-frac2N/A
mul-1-negN/A
+-lowering-+.f64N/A
Simplified85.8%
if -0.050000000000000003 < z < 105000Initial program 99.7%
Taylor expanded in z around 0
*-commutativeN/A
*-lowering-*.f6498.8%
Simplified98.8%
Taylor expanded in z around 0
*-commutativeN/A
associate-*r*N/A
distribute-rgt-inN/A
+-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f6492.6%
Simplified92.6%
if 105000 < z Initial program 14.4%
Taylor expanded in z around inf
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6486.6%
Simplified86.6%
Final simplification89.3%
(FPCore (x y z t a b)
:precision binary64
(if (<= z -720.0)
(+ x (+ (* y 3.13060547623) (/ (* y -36.52704169880642) z)))
(if (<= z 1.3e-14)
(+ x (* y (* b 1.6453555072203998)))
(+ x (* y 3.13060547623)))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (z <= -720.0) {
tmp = x + ((y * 3.13060547623) + ((y * -36.52704169880642) / z));
} else if (z <= 1.3e-14) {
tmp = x + (y * (b * 1.6453555072203998));
} else {
tmp = x + (y * 3.13060547623);
}
return tmp;
}
real(8) function code(x, y, z, t, a, b)
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) :: tmp
if (z <= (-720.0d0)) then
tmp = x + ((y * 3.13060547623d0) + ((y * (-36.52704169880642d0)) / z))
else if (z <= 1.3d-14) then
tmp = x + (y * (b * 1.6453555072203998d0))
else
tmp = x + (y * 3.13060547623d0)
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (z <= -720.0) {
tmp = x + ((y * 3.13060547623) + ((y * -36.52704169880642) / z));
} else if (z <= 1.3e-14) {
tmp = x + (y * (b * 1.6453555072203998));
} else {
tmp = x + (y * 3.13060547623);
}
return tmp;
}
def code(x, y, z, t, a, b): tmp = 0 if z <= -720.0: tmp = x + ((y * 3.13060547623) + ((y * -36.52704169880642) / z)) elif z <= 1.3e-14: tmp = x + (y * (b * 1.6453555072203998)) else: tmp = x + (y * 3.13060547623) return tmp
function code(x, y, z, t, a, b) tmp = 0.0 if (z <= -720.0) tmp = Float64(x + Float64(Float64(y * 3.13060547623) + Float64(Float64(y * -36.52704169880642) / z))); elseif (z <= 1.3e-14) tmp = Float64(x + Float64(y * Float64(b * 1.6453555072203998))); else tmp = Float64(x + Float64(y * 3.13060547623)); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) tmp = 0.0; if (z <= -720.0) tmp = x + ((y * 3.13060547623) + ((y * -36.52704169880642) / z)); elseif (z <= 1.3e-14) tmp = x + (y * (b * 1.6453555072203998)); else tmp = x + (y * 3.13060547623); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[z, -720.0], N[(x + N[(N[(y * 3.13060547623), $MachinePrecision] + N[(N[(y * -36.52704169880642), $MachinePrecision] / z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[z, 1.3e-14], N[(x + N[(y * N[(b * 1.6453555072203998), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x + N[(y * 3.13060547623), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;z \leq -720:\\
\;\;\;\;x + \left(y \cdot 3.13060547623 + \frac{y \cdot -36.52704169880642}{z}\right)\\
\mathbf{elif}\;z \leq 1.3 \cdot 10^{-14}:\\
\;\;\;\;x + y \cdot \left(b \cdot 1.6453555072203998\right)\\
\mathbf{else}:\\
\;\;\;\;x + y \cdot 3.13060547623\\
\end{array}
\end{array}
if z < -720Initial program 23.9%
Taylor expanded in z around inf
associate--l+N/A
associate--l+N/A
+-commutativeN/A
distribute-rgt-out--N/A
metadata-evalN/A
metadata-evalN/A
metadata-evalN/A
times-fracN/A
*-commutativeN/A
distribute-rgt-out--N/A
mul-1-negN/A
distribute-neg-frac2N/A
mul-1-negN/A
+-lowering-+.f64N/A
Simplified87.1%
if -720 < z < 1.29999999999999998e-14Initial program 99.7%
Taylor expanded in z around 0
+-lowering-+.f64N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f6478.1%
Simplified78.1%
if 1.29999999999999998e-14 < z Initial program 20.2%
Taylor expanded in z around inf
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6484.4%
Simplified84.4%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (+ x (* y 3.13060547623))))
(if (<= z -400.0)
t_1
(if (<= z 1.3e-14) (+ x (* y (* b 1.6453555072203998))) t_1))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = x + (y * 3.13060547623);
double tmp;
if (z <= -400.0) {
tmp = t_1;
} else if (z <= 1.3e-14) {
tmp = x + (y * (b * 1.6453555072203998));
} else {
tmp = t_1;
}
return tmp;
}
real(8) function code(x, y, z, t, a, b)
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) :: t_1
real(8) :: tmp
t_1 = x + (y * 3.13060547623d0)
if (z <= (-400.0d0)) then
tmp = t_1
else if (z <= 1.3d-14) then
tmp = x + (y * (b * 1.6453555072203998d0))
else
tmp = t_1
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b) {
double t_1 = x + (y * 3.13060547623);
double tmp;
if (z <= -400.0) {
tmp = t_1;
} else if (z <= 1.3e-14) {
tmp = x + (y * (b * 1.6453555072203998));
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = x + (y * 3.13060547623) tmp = 0 if z <= -400.0: tmp = t_1 elif z <= 1.3e-14: tmp = x + (y * (b * 1.6453555072203998)) else: tmp = t_1 return tmp
function code(x, y, z, t, a, b) t_1 = Float64(x + Float64(y * 3.13060547623)) tmp = 0.0 if (z <= -400.0) tmp = t_1; elseif (z <= 1.3e-14) tmp = Float64(x + Float64(y * Float64(b * 1.6453555072203998))); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a, b) t_1 = x + (y * 3.13060547623); tmp = 0.0; if (z <= -400.0) tmp = t_1; elseif (z <= 1.3e-14) tmp = x + (y * (b * 1.6453555072203998)); else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(x + N[(y * 3.13060547623), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[z, -400.0], t$95$1, If[LessEqual[z, 1.3e-14], N[(x + N[(y * N[(b * 1.6453555072203998), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x + y \cdot 3.13060547623\\
\mathbf{if}\;z \leq -400:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;z \leq 1.3 \cdot 10^{-14}:\\
\;\;\;\;x + y \cdot \left(b \cdot 1.6453555072203998\right)\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if z < -400 or 1.29999999999999998e-14 < z Initial program 21.8%
Taylor expanded in z around inf
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6485.6%
Simplified85.6%
if -400 < z < 1.29999999999999998e-14Initial program 99.7%
Taylor expanded in z around 0
+-lowering-+.f64N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f6478.1%
Simplified78.1%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (+ x (* y 3.13060547623))))
(if (<= z -60.0)
t_1
(if (<= z 1.3e-14) (+ x (* 1.6453555072203998 (* y b))) t_1))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = x + (y * 3.13060547623);
double tmp;
if (z <= -60.0) {
tmp = t_1;
} else if (z <= 1.3e-14) {
tmp = x + (1.6453555072203998 * (y * b));
} else {
tmp = t_1;
}
return tmp;
}
real(8) function code(x, y, z, t, a, b)
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) :: t_1
real(8) :: tmp
t_1 = x + (y * 3.13060547623d0)
if (z <= (-60.0d0)) then
tmp = t_1
else if (z <= 1.3d-14) then
tmp = x + (1.6453555072203998d0 * (y * b))
else
tmp = t_1
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b) {
double t_1 = x + (y * 3.13060547623);
double tmp;
if (z <= -60.0) {
tmp = t_1;
} else if (z <= 1.3e-14) {
tmp = x + (1.6453555072203998 * (y * b));
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = x + (y * 3.13060547623) tmp = 0 if z <= -60.0: tmp = t_1 elif z <= 1.3e-14: tmp = x + (1.6453555072203998 * (y * b)) else: tmp = t_1 return tmp
function code(x, y, z, t, a, b) t_1 = Float64(x + Float64(y * 3.13060547623)) tmp = 0.0 if (z <= -60.0) tmp = t_1; elseif (z <= 1.3e-14) tmp = Float64(x + Float64(1.6453555072203998 * Float64(y * b))); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a, b) t_1 = x + (y * 3.13060547623); tmp = 0.0; if (z <= -60.0) tmp = t_1; elseif (z <= 1.3e-14) tmp = x + (1.6453555072203998 * (y * b)); else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(x + N[(y * 3.13060547623), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[z, -60.0], t$95$1, If[LessEqual[z, 1.3e-14], N[(x + N[(1.6453555072203998 * N[(y * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x + y \cdot 3.13060547623\\
\mathbf{if}\;z \leq -60:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;z \leq 1.3 \cdot 10^{-14}:\\
\;\;\;\;x + 1.6453555072203998 \cdot \left(y \cdot b\right)\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if z < -60 or 1.29999999999999998e-14 < z Initial program 21.8%
Taylor expanded in z around inf
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6485.6%
Simplified85.6%
if -60 < z < 1.29999999999999998e-14Initial program 99.7%
Taylor expanded in z around 0
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
associate-*r*N/A
*-commutativeN/A
associate-*l*N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f6493.0%
Simplified93.0%
Taylor expanded in z around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f6478.0%
Simplified78.0%
(FPCore (x y z t a b) :precision binary64 (let* ((t_1 (+ x (* y 3.13060547623)))) (if (<= z -7.5e-36) t_1 (if (<= z 3.5e-119) x t_1))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = x + (y * 3.13060547623);
double tmp;
if (z <= -7.5e-36) {
tmp = t_1;
} else if (z <= 3.5e-119) {
tmp = x;
} else {
tmp = t_1;
}
return tmp;
}
real(8) function code(x, y, z, t, a, b)
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) :: t_1
real(8) :: tmp
t_1 = x + (y * 3.13060547623d0)
if (z <= (-7.5d-36)) then
tmp = t_1
else if (z <= 3.5d-119) then
tmp = x
else
tmp = t_1
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b) {
double t_1 = x + (y * 3.13060547623);
double tmp;
if (z <= -7.5e-36) {
tmp = t_1;
} else if (z <= 3.5e-119) {
tmp = x;
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = x + (y * 3.13060547623) tmp = 0 if z <= -7.5e-36: tmp = t_1 elif z <= 3.5e-119: tmp = x else: tmp = t_1 return tmp
function code(x, y, z, t, a, b) t_1 = Float64(x + Float64(y * 3.13060547623)) tmp = 0.0 if (z <= -7.5e-36) tmp = t_1; elseif (z <= 3.5e-119) tmp = x; else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a, b) t_1 = x + (y * 3.13060547623); tmp = 0.0; if (z <= -7.5e-36) tmp = t_1; elseif (z <= 3.5e-119) tmp = x; else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(x + N[(y * 3.13060547623), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[z, -7.5e-36], t$95$1, If[LessEqual[z, 3.5e-119], x, t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x + y \cdot 3.13060547623\\
\mathbf{if}\;z \leq -7.5 \cdot 10^{-36}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;z \leq 3.5 \cdot 10^{-119}:\\
\;\;\;\;x\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if z < -7.49999999999999972e-36 or 3.5e-119 < z Initial program 37.2%
Taylor expanded in z around inf
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6475.2%
Simplified75.2%
if -7.49999999999999972e-36 < z < 3.5e-119Initial program 99.7%
Taylor expanded in x around inf
Simplified50.1%
(FPCore (x y z t a b) :precision binary64 (if (<= y -5.5e+131) (* y 3.13060547623) (if (<= y 1.3e+104) x (* y 3.13060547623))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (y <= -5.5e+131) {
tmp = y * 3.13060547623;
} else if (y <= 1.3e+104) {
tmp = x;
} else {
tmp = y * 3.13060547623;
}
return tmp;
}
real(8) function code(x, y, z, t, a, b)
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) :: tmp
if (y <= (-5.5d+131)) then
tmp = y * 3.13060547623d0
else if (y <= 1.3d+104) then
tmp = x
else
tmp = y * 3.13060547623d0
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (y <= -5.5e+131) {
tmp = y * 3.13060547623;
} else if (y <= 1.3e+104) {
tmp = x;
} else {
tmp = y * 3.13060547623;
}
return tmp;
}
def code(x, y, z, t, a, b): tmp = 0 if y <= -5.5e+131: tmp = y * 3.13060547623 elif y <= 1.3e+104: tmp = x else: tmp = y * 3.13060547623 return tmp
function code(x, y, z, t, a, b) tmp = 0.0 if (y <= -5.5e+131) tmp = Float64(y * 3.13060547623); elseif (y <= 1.3e+104) tmp = x; else tmp = Float64(y * 3.13060547623); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) tmp = 0.0; if (y <= -5.5e+131) tmp = y * 3.13060547623; elseif (y <= 1.3e+104) tmp = x; else tmp = y * 3.13060547623; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[y, -5.5e+131], N[(y * 3.13060547623), $MachinePrecision], If[LessEqual[y, 1.3e+104], x, N[(y * 3.13060547623), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -5.5 \cdot 10^{+131}:\\
\;\;\;\;y \cdot 3.13060547623\\
\mathbf{elif}\;y \leq 1.3 \cdot 10^{+104}:\\
\;\;\;\;x\\
\mathbf{else}:\\
\;\;\;\;y \cdot 3.13060547623\\
\end{array}
\end{array}
if y < -5.49999999999999971e131 or 1.3e104 < y Initial program 49.5%
Taylor expanded in z around inf
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6446.2%
Simplified46.2%
Taylor expanded in x around 0
*-lowering-*.f6440.2%
Simplified40.2%
if -5.49999999999999971e131 < y < 1.3e104Initial program 61.9%
Taylor expanded in x around inf
Simplified62.6%
Final simplification55.0%
(FPCore (x y z t a b) :precision binary64 x)
double code(double x, double y, double z, double t, double a, double b) {
return x;
}
real(8) function code(x, y, z, t, a, b)
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
code = x
end function
public static double code(double x, double y, double z, double t, double a, double b) {
return x;
}
def code(x, y, z, t, a, b): return x
function code(x, y, z, t, a, b) return x end
function tmp = code(x, y, z, t, a, b) tmp = x; end
code[x_, y_, z_, t_, a_, b_] := x
\begin{array}{l}
\\
x
\end{array}
Initial program 57.7%
Taylor expanded in x around inf
Simplified45.2%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1
(+
x
(*
(+ (- 3.13060547623 (/ 36.527041698806414 z)) (/ t (* z z)))
(/ y 1.0)))))
(if (< z -6.499344996252632e+53)
t_1
(if (< z 7.066965436914287e+59)
(+
x
(/
y
(/
(+
(*
(+ (* (+ (* (+ z 15.234687407) z) 31.4690115749) z) 11.9400905721)
z)
0.607771387771)
(+
(* (+ (* (+ (* (+ (* z 3.13060547623) 11.1667541262) z) t) z) a) z)
b))))
t_1))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = x + (((3.13060547623 - (36.527041698806414 / z)) + (t / (z * z))) * (y / 1.0));
double tmp;
if (z < -6.499344996252632e+53) {
tmp = t_1;
} else if (z < 7.066965436914287e+59) {
tmp = x + (y / ((((((((z + 15.234687407) * z) + 31.4690115749) * z) + 11.9400905721) * z) + 0.607771387771) / ((((((((z * 3.13060547623) + 11.1667541262) * z) + t) * z) + a) * z) + b)));
} else {
tmp = t_1;
}
return tmp;
}
real(8) function code(x, y, z, t, a, b)
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) :: t_1
real(8) :: tmp
t_1 = x + (((3.13060547623d0 - (36.527041698806414d0 / z)) + (t / (z * z))) * (y / 1.0d0))
if (z < (-6.499344996252632d+53)) then
tmp = t_1
else if (z < 7.066965436914287d+59) then
tmp = x + (y / ((((((((z + 15.234687407d0) * z) + 31.4690115749d0) * z) + 11.9400905721d0) * z) + 0.607771387771d0) / ((((((((z * 3.13060547623d0) + 11.1667541262d0) * z) + t) * z) + a) * z) + b)))
else
tmp = t_1
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b) {
double t_1 = x + (((3.13060547623 - (36.527041698806414 / z)) + (t / (z * z))) * (y / 1.0));
double tmp;
if (z < -6.499344996252632e+53) {
tmp = t_1;
} else if (z < 7.066965436914287e+59) {
tmp = x + (y / ((((((((z + 15.234687407) * z) + 31.4690115749) * z) + 11.9400905721) * z) + 0.607771387771) / ((((((((z * 3.13060547623) + 11.1667541262) * z) + t) * z) + a) * z) + b)));
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = x + (((3.13060547623 - (36.527041698806414 / z)) + (t / (z * z))) * (y / 1.0)) tmp = 0 if z < -6.499344996252632e+53: tmp = t_1 elif z < 7.066965436914287e+59: tmp = x + (y / ((((((((z + 15.234687407) * z) + 31.4690115749) * z) + 11.9400905721) * z) + 0.607771387771) / ((((((((z * 3.13060547623) + 11.1667541262) * z) + t) * z) + a) * z) + b))) else: tmp = t_1 return tmp
function code(x, y, z, t, a, b) t_1 = Float64(x + Float64(Float64(Float64(3.13060547623 - Float64(36.527041698806414 / z)) + Float64(t / Float64(z * z))) * Float64(y / 1.0))) tmp = 0.0 if (z < -6.499344996252632e+53) tmp = t_1; elseif (z < 7.066965436914287e+59) tmp = Float64(x + Float64(y / Float64(Float64(Float64(Float64(Float64(Float64(Float64(Float64(z + 15.234687407) * z) + 31.4690115749) * z) + 11.9400905721) * z) + 0.607771387771) / Float64(Float64(Float64(Float64(Float64(Float64(Float64(Float64(z * 3.13060547623) + 11.1667541262) * z) + t) * z) + a) * z) + b)))); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a, b) t_1 = x + (((3.13060547623 - (36.527041698806414 / z)) + (t / (z * z))) * (y / 1.0)); tmp = 0.0; if (z < -6.499344996252632e+53) tmp = t_1; elseif (z < 7.066965436914287e+59) tmp = x + (y / ((((((((z + 15.234687407) * z) + 31.4690115749) * z) + 11.9400905721) * z) + 0.607771387771) / ((((((((z * 3.13060547623) + 11.1667541262) * z) + t) * z) + a) * z) + b))); else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(x + N[(N[(N[(3.13060547623 - N[(36.527041698806414 / z), $MachinePrecision]), $MachinePrecision] + N[(t / N[(z * z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(y / 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[Less[z, -6.499344996252632e+53], t$95$1, If[Less[z, 7.066965436914287e+59], N[(x + N[(y / N[(N[(N[(N[(N[(N[(N[(N[(z + 15.234687407), $MachinePrecision] * z), $MachinePrecision] + 31.4690115749), $MachinePrecision] * z), $MachinePrecision] + 11.9400905721), $MachinePrecision] * z), $MachinePrecision] + 0.607771387771), $MachinePrecision] / N[(N[(N[(N[(N[(N[(N[(N[(z * 3.13060547623), $MachinePrecision] + 11.1667541262), $MachinePrecision] * z), $MachinePrecision] + t), $MachinePrecision] * z), $MachinePrecision] + a), $MachinePrecision] * z), $MachinePrecision] + b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x + \left(\left(3.13060547623 - \frac{36.527041698806414}{z}\right) + \frac{t}{z \cdot z}\right) \cdot \frac{y}{1}\\
\mathbf{if}\;z < -6.499344996252632 \cdot 10^{+53}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;z < 7.066965436914287 \cdot 10^{+59}:\\
\;\;\;\;x + \frac{y}{\frac{\left(\left(\left(z + 15.234687407\right) \cdot z + 31.4690115749\right) \cdot z + 11.9400905721\right) \cdot z + 0.607771387771}{\left(\left(\left(z \cdot 3.13060547623 + 11.1667541262\right) \cdot z + t\right) \cdot z + a\right) \cdot z + b}}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
herbie shell --seed 2024191
(FPCore (x y z t a b)
:name "Numeric.SpecFunctions:logGamma from math-functions-0.1.5.2, D"
:precision binary64
:alt
(! :herbie-platform default (if (< z -649934499625263200000000000000000000000000000000000000) (+ x (* (+ (- 313060547623/100000000000 (/ 18263520849403207/500000000000000 z)) (/ t (* z z))) (/ y 1))) (if (< z 706696543691428700000000000000000000000000000000000000000000) (+ x (/ y (/ (+ (* (+ (* (+ (* (+ z 15234687407/1000000000) z) 314690115749/10000000000) z) 119400905721/10000000000) z) 607771387771/1000000000000) (+ (* (+ (* (+ (* (+ (* z 313060547623/100000000000) 55833770631/5000000000) z) t) z) a) z) b)))) (+ x (* (+ (- 313060547623/100000000000 (/ 18263520849403207/500000000000000 z)) (/ t (* z z))) (/ y 1))))))
(+ x (/ (* y (+ (* (+ (* (+ (* (+ (* z 3.13060547623) 11.1667541262) z) t) z) a) z) b)) (+ (* (+ (* (+ (* (+ z 15.234687407) z) 31.4690115749) z) 11.9400905721) z) 0.607771387771))))