
(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 18 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
(+
(/ 3.7269864963038164 z)
(+ 0.31942702700572795 (/ -3.241970391368047 (* z z))))))
(if (<= z -1.45e+33)
(+
x
(fma
0.10203362558171805
(/ (* (/ y z) (/ t z)) (pow t_1 2.0))
(/ y t_1)))
(if (<= z 2.1e+23)
(+
x
(/
y
(/
(fma
(fma (fma (+ z 15.234687407) z 31.4690115749) z 11.9400905721)
z
0.607771387771)
(fma (fma (fma (fma z 3.13060547623 11.1667541262) z t) z a) z b))))
(+
x
(-
(-
(fma y 3.13060547623 (/ y (/ (* z z) t)))
(/ (* y 36.52704169880642) z))
(fma
98.5170599679272
(/ y (* z z))
(/ (* (* y 36.52704169880642) -15.234687407) (* z z)))))))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = (3.7269864963038164 / z) + (0.31942702700572795 + (-3.241970391368047 / (z * z)));
double tmp;
if (z <= -1.45e+33) {
tmp = x + fma(0.10203362558171805, (((y / z) * (t / z)) / pow(t_1, 2.0)), (y / t_1));
} else if (z <= 2.1e+23) {
tmp = x + (y / (fma(fma(fma((z + 15.234687407), z, 31.4690115749), z, 11.9400905721), z, 0.607771387771) / fma(fma(fma(fma(z, 3.13060547623, 11.1667541262), z, t), z, a), z, b)));
} else {
tmp = x + ((fma(y, 3.13060547623, (y / ((z * z) / t))) - ((y * 36.52704169880642) / z)) - fma(98.5170599679272, (y / (z * z)), (((y * 36.52704169880642) * -15.234687407) / (z * z))));
}
return tmp;
}
function code(x, y, z, t, a, b) t_1 = Float64(Float64(3.7269864963038164 / z) + Float64(0.31942702700572795 + Float64(-3.241970391368047 / Float64(z * z)))) tmp = 0.0 if (z <= -1.45e+33) tmp = Float64(x + fma(0.10203362558171805, Float64(Float64(Float64(y / z) * Float64(t / z)) / (t_1 ^ 2.0)), Float64(y / t_1))); elseif (z <= 2.1e+23) tmp = Float64(x + Float64(y / Float64(fma(fma(fma(Float64(z + 15.234687407), z, 31.4690115749), z, 11.9400905721), z, 0.607771387771) / fma(fma(fma(fma(z, 3.13060547623, 11.1667541262), z, t), z, a), z, b)))); else tmp = Float64(x + Float64(Float64(fma(y, 3.13060547623, Float64(y / Float64(Float64(z * z) / t))) - Float64(Float64(y * 36.52704169880642) / z)) - fma(98.5170599679272, Float64(y / Float64(z * z)), Float64(Float64(Float64(y * 36.52704169880642) * -15.234687407) / Float64(z * z))))); end return tmp end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(N[(3.7269864963038164 / z), $MachinePrecision] + N[(0.31942702700572795 + N[(-3.241970391368047 / N[(z * z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[z, -1.45e+33], N[(x + N[(0.10203362558171805 * N[(N[(N[(y / z), $MachinePrecision] * N[(t / z), $MachinePrecision]), $MachinePrecision] / N[Power[t$95$1, 2.0], $MachinePrecision]), $MachinePrecision] + N[(y / t$95$1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[z, 2.1e+23], N[(x + N[(y / N[(N[(N[(N[(N[(z + 15.234687407), $MachinePrecision] * z + 31.4690115749), $MachinePrecision] * z + 11.9400905721), $MachinePrecision] * z + 0.607771387771), $MachinePrecision] / N[(N[(N[(N[(z * 3.13060547623 + 11.1667541262), $MachinePrecision] * z + t), $MachinePrecision] * z + a), $MachinePrecision] * z + b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x + N[(N[(N[(y * 3.13060547623 + N[(y / N[(N[(z * z), $MachinePrecision] / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - N[(N[(y * 36.52704169880642), $MachinePrecision] / z), $MachinePrecision]), $MachinePrecision] - N[(98.5170599679272 * N[(y / N[(z * z), $MachinePrecision]), $MachinePrecision] + N[(N[(N[(y * 36.52704169880642), $MachinePrecision] * -15.234687407), $MachinePrecision] / N[(z * z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{3.7269864963038164}{z} + \left(0.31942702700572795 + \frac{-3.241970391368047}{z \cdot z}\right)\\
\mathbf{if}\;z \leq -1.45 \cdot 10^{+33}:\\
\;\;\;\;x + \mathsf{fma}\left(0.10203362558171805, \frac{\frac{y}{z} \cdot \frac{t}{z}}{{t_1}^{2}}, \frac{y}{t_1}\right)\\
\mathbf{elif}\;z \leq 2.1 \cdot 10^{+23}:\\
\;\;\;\;x + \frac{y}{\frac{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(z + 15.234687407, z, 31.4690115749\right), z, 11.9400905721\right), z, 0.607771387771\right)}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(z, 3.13060547623, 11.1667541262\right), z, t\right), z, a\right), z, b\right)}}\\
\mathbf{else}:\\
\;\;\;\;x + \left(\left(\mathsf{fma}\left(y, 3.13060547623, \frac{y}{\frac{z \cdot z}{t}}\right) - \frac{y \cdot 36.52704169880642}{z}\right) - \mathsf{fma}\left(98.5170599679272, \frac{y}{z \cdot z}, \frac{\left(y \cdot 36.52704169880642\right) \cdot -15.234687407}{z \cdot z}\right)\right)\\
\end{array}
\end{array}
if z < -1.45000000000000012e33Initial program 2.2%
associate-/l*8.1%
fma-def8.1%
fma-def8.1%
fma-def8.1%
fma-def8.1%
fma-def8.1%
fma-def8.1%
fma-def8.1%
Simplified8.1%
Taylor expanded in z around inf 86.3%
associate-*r/86.3%
metadata-eval86.3%
mul-1-neg86.3%
*-commutative86.3%
unpow286.3%
Simplified86.3%
Taylor expanded in t around 0 90.9%
+-commutative90.9%
fma-def90.9%
Simplified98.5%
if -1.45000000000000012e33 < z < 2.1000000000000001e23Initial program 99.0%
associate-/l*99.7%
fma-def99.7%
fma-def99.7%
fma-def99.7%
fma-def99.7%
fma-def99.7%
fma-def99.7%
fma-def99.7%
Simplified99.7%
if 2.1000000000000001e23 < z Initial program 6.9%
associate-*l/12.7%
*-commutative12.7%
fma-def12.7%
*-commutative12.7%
fma-def12.7%
*-commutative12.7%
fma-def12.7%
*-commutative12.7%
fma-def12.7%
Simplified12.7%
Taylor expanded in z around -inf 89.1%
+-commutative89.1%
mul-1-neg89.1%
unsub-neg89.1%
+-commutative89.1%
*-commutative89.1%
fma-def89.1%
associate-/l*98.3%
unpow298.3%
distribute-rgt-out--98.3%
metadata-eval98.3%
+-commutative98.3%
fma-def98.3%
Simplified98.3%
Final simplification99.1%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1
(+
(/ 3.7269864963038164 z)
(+ 0.31942702700572795 (/ -3.241970391368047 (* z z))))))
(if (<= z -1.35e+33)
(+
x
(fma
0.10203362558171805
(/ (* (/ y z) (/ t z)) (pow t_1 2.0))
(/ y t_1)))
(if (<= z 2.1e+23)
(+
x
(*
(/
y
(fma
z
(fma z (fma z (+ z 15.234687407) 31.4690115749) 11.9400905721)
0.607771387771))
(fma z (fma z (fma z (fma z 3.13060547623 11.1667541262) t) a) b)))
(+
x
(-
(-
(fma y 3.13060547623 (/ y (/ (* z z) t)))
(/ (* y 36.52704169880642) z))
(fma
98.5170599679272
(/ y (* z z))
(/ (* (* y 36.52704169880642) -15.234687407) (* z z)))))))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = (3.7269864963038164 / z) + (0.31942702700572795 + (-3.241970391368047 / (z * z)));
double tmp;
if (z <= -1.35e+33) {
tmp = x + fma(0.10203362558171805, (((y / z) * (t / z)) / pow(t_1, 2.0)), (y / t_1));
} else if (z <= 2.1e+23) {
tmp = x + ((y / fma(z, fma(z, fma(z, (z + 15.234687407), 31.4690115749), 11.9400905721), 0.607771387771)) * fma(z, fma(z, fma(z, fma(z, 3.13060547623, 11.1667541262), t), a), b));
} else {
tmp = x + ((fma(y, 3.13060547623, (y / ((z * z) / t))) - ((y * 36.52704169880642) / z)) - fma(98.5170599679272, (y / (z * z)), (((y * 36.52704169880642) * -15.234687407) / (z * z))));
}
return tmp;
}
function code(x, y, z, t, a, b) t_1 = Float64(Float64(3.7269864963038164 / z) + Float64(0.31942702700572795 + Float64(-3.241970391368047 / Float64(z * z)))) tmp = 0.0 if (z <= -1.35e+33) tmp = Float64(x + fma(0.10203362558171805, Float64(Float64(Float64(y / z) * Float64(t / z)) / (t_1 ^ 2.0)), Float64(y / t_1))); elseif (z <= 2.1e+23) tmp = Float64(x + Float64(Float64(y / fma(z, fma(z, fma(z, Float64(z + 15.234687407), 31.4690115749), 11.9400905721), 0.607771387771)) * fma(z, fma(z, fma(z, fma(z, 3.13060547623, 11.1667541262), t), a), b))); else tmp = Float64(x + Float64(Float64(fma(y, 3.13060547623, Float64(y / Float64(Float64(z * z) / t))) - Float64(Float64(y * 36.52704169880642) / z)) - fma(98.5170599679272, Float64(y / Float64(z * z)), Float64(Float64(Float64(y * 36.52704169880642) * -15.234687407) / Float64(z * z))))); end return tmp end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(N[(3.7269864963038164 / z), $MachinePrecision] + N[(0.31942702700572795 + N[(-3.241970391368047 / N[(z * z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[z, -1.35e+33], N[(x + N[(0.10203362558171805 * N[(N[(N[(y / z), $MachinePrecision] * N[(t / z), $MachinePrecision]), $MachinePrecision] / N[Power[t$95$1, 2.0], $MachinePrecision]), $MachinePrecision] + N[(y / t$95$1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[z, 2.1e+23], N[(x + N[(N[(y / N[(z * N[(z * N[(z * N[(z + 15.234687407), $MachinePrecision] + 31.4690115749), $MachinePrecision] + 11.9400905721), $MachinePrecision] + 0.607771387771), $MachinePrecision]), $MachinePrecision] * N[(z * N[(z * N[(z * N[(z * 3.13060547623 + 11.1667541262), $MachinePrecision] + t), $MachinePrecision] + a), $MachinePrecision] + b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x + N[(N[(N[(y * 3.13060547623 + N[(y / N[(N[(z * z), $MachinePrecision] / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - N[(N[(y * 36.52704169880642), $MachinePrecision] / z), $MachinePrecision]), $MachinePrecision] - N[(98.5170599679272 * N[(y / N[(z * z), $MachinePrecision]), $MachinePrecision] + N[(N[(N[(y * 36.52704169880642), $MachinePrecision] * -15.234687407), $MachinePrecision] / N[(z * z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{3.7269864963038164}{z} + \left(0.31942702700572795 + \frac{-3.241970391368047}{z \cdot z}\right)\\
\mathbf{if}\;z \leq -1.35 \cdot 10^{+33}:\\
\;\;\;\;x + \mathsf{fma}\left(0.10203362558171805, \frac{\frac{y}{z} \cdot \frac{t}{z}}{{t_1}^{2}}, \frac{y}{t_1}\right)\\
\mathbf{elif}\;z \leq 2.1 \cdot 10^{+23}:\\
\;\;\;\;x + \frac{y}{\mathsf{fma}\left(z, \mathsf{fma}\left(z, \mathsf{fma}\left(z, z + 15.234687407, 31.4690115749\right), 11.9400905721\right), 0.607771387771\right)} \cdot \mathsf{fma}\left(z, \mathsf{fma}\left(z, \mathsf{fma}\left(z, \mathsf{fma}\left(z, 3.13060547623, 11.1667541262\right), t\right), a\right), b\right)\\
\mathbf{else}:\\
\;\;\;\;x + \left(\left(\mathsf{fma}\left(y, 3.13060547623, \frac{y}{\frac{z \cdot z}{t}}\right) - \frac{y \cdot 36.52704169880642}{z}\right) - \mathsf{fma}\left(98.5170599679272, \frac{y}{z \cdot z}, \frac{\left(y \cdot 36.52704169880642\right) \cdot -15.234687407}{z \cdot z}\right)\right)\\
\end{array}
\end{array}
if z < -1.34999999999999996e33Initial program 2.2%
associate-/l*8.1%
fma-def8.1%
fma-def8.1%
fma-def8.1%
fma-def8.1%
fma-def8.1%
fma-def8.1%
fma-def8.1%
Simplified8.1%
Taylor expanded in z around inf 86.3%
associate-*r/86.3%
metadata-eval86.3%
mul-1-neg86.3%
*-commutative86.3%
unpow286.3%
Simplified86.3%
Taylor expanded in t around 0 90.9%
+-commutative90.9%
fma-def90.9%
Simplified98.5%
if -1.34999999999999996e33 < z < 2.1000000000000001e23Initial program 99.0%
associate-*l/99.6%
*-commutative99.6%
fma-def99.6%
*-commutative99.6%
fma-def99.6%
*-commutative99.6%
fma-def99.6%
*-commutative99.6%
fma-def99.6%
Simplified99.6%
if 2.1000000000000001e23 < z Initial program 6.9%
associate-*l/12.7%
*-commutative12.7%
fma-def12.7%
*-commutative12.7%
fma-def12.7%
*-commutative12.7%
fma-def12.7%
*-commutative12.7%
fma-def12.7%
Simplified12.7%
Taylor expanded in z around -inf 89.1%
+-commutative89.1%
mul-1-neg89.1%
unsub-neg89.1%
+-commutative89.1%
*-commutative89.1%
fma-def89.1%
associate-/l*98.3%
unpow298.3%
distribute-rgt-out--98.3%
metadata-eval98.3%
+-commutative98.3%
fma-def98.3%
Simplified98.3%
Final simplification99.0%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1
(+
(/ 3.7269864963038164 z)
(+ 0.31942702700572795 (/ -3.241970391368047 (* z z))))))
(if (<= z -1.15e+33)
(+
x
(fma
0.10203362558171805
(/ (* (/ y z) (/ t z)) (pow t_1 2.0))
(/ y t_1)))
(if (<= z 1.45e+23)
(+
x
(/
(*
y
(+
b
(*
z
(+ a (* z (+ t (* z (+ 11.1667541262 (* z 3.13060547623)))))))))
(+
0.607771387771
(*
z
(+
11.9400905721
(* z (+ 31.4690115749 (* z (+ z 15.234687407)))))))))
(+
x
(-
(-
(fma y 3.13060547623 (/ y (/ (* z z) t)))
(/ (* y 36.52704169880642) z))
(fma
98.5170599679272
(/ y (* z z))
(/ (* (* y 36.52704169880642) -15.234687407) (* z z)))))))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = (3.7269864963038164 / z) + (0.31942702700572795 + (-3.241970391368047 / (z * z)));
double tmp;
if (z <= -1.15e+33) {
tmp = x + fma(0.10203362558171805, (((y / z) * (t / z)) / pow(t_1, 2.0)), (y / t_1));
} else if (z <= 1.45e+23) {
tmp = x + ((y * (b + (z * (a + (z * (t + (z * (11.1667541262 + (z * 3.13060547623))))))))) / (0.607771387771 + (z * (11.9400905721 + (z * (31.4690115749 + (z * (z + 15.234687407))))))));
} else {
tmp = x + ((fma(y, 3.13060547623, (y / ((z * z) / t))) - ((y * 36.52704169880642) / z)) - fma(98.5170599679272, (y / (z * z)), (((y * 36.52704169880642) * -15.234687407) / (z * z))));
}
return tmp;
}
function code(x, y, z, t, a, b) t_1 = Float64(Float64(3.7269864963038164 / z) + Float64(0.31942702700572795 + Float64(-3.241970391368047 / Float64(z * z)))) tmp = 0.0 if (z <= -1.15e+33) tmp = Float64(x + fma(0.10203362558171805, Float64(Float64(Float64(y / z) * Float64(t / z)) / (t_1 ^ 2.0)), Float64(y / t_1))); elseif (z <= 1.45e+23) tmp = Float64(x + Float64(Float64(y * Float64(b + Float64(z * Float64(a + Float64(z * Float64(t + Float64(z * Float64(11.1667541262 + Float64(z * 3.13060547623))))))))) / Float64(0.607771387771 + Float64(z * Float64(11.9400905721 + Float64(z * Float64(31.4690115749 + Float64(z * Float64(z + 15.234687407))))))))); else tmp = Float64(x + Float64(Float64(fma(y, 3.13060547623, Float64(y / Float64(Float64(z * z) / t))) - Float64(Float64(y * 36.52704169880642) / z)) - fma(98.5170599679272, Float64(y / Float64(z * z)), Float64(Float64(Float64(y * 36.52704169880642) * -15.234687407) / Float64(z * z))))); end return tmp end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(N[(3.7269864963038164 / z), $MachinePrecision] + N[(0.31942702700572795 + N[(-3.241970391368047 / N[(z * z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[z, -1.15e+33], N[(x + N[(0.10203362558171805 * N[(N[(N[(y / z), $MachinePrecision] * N[(t / z), $MachinePrecision]), $MachinePrecision] / N[Power[t$95$1, 2.0], $MachinePrecision]), $MachinePrecision] + N[(y / t$95$1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[z, 1.45e+23], N[(x + N[(N[(y * N[(b + N[(z * N[(a + N[(z * N[(t + N[(z * N[(11.1667541262 + N[(z * 3.13060547623), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / 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], N[(x + N[(N[(N[(y * 3.13060547623 + N[(y / N[(N[(z * z), $MachinePrecision] / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - N[(N[(y * 36.52704169880642), $MachinePrecision] / z), $MachinePrecision]), $MachinePrecision] - N[(98.5170599679272 * N[(y / N[(z * z), $MachinePrecision]), $MachinePrecision] + N[(N[(N[(y * 36.52704169880642), $MachinePrecision] * -15.234687407), $MachinePrecision] / N[(z * z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{3.7269864963038164}{z} + \left(0.31942702700572795 + \frac{-3.241970391368047}{z \cdot z}\right)\\
\mathbf{if}\;z \leq -1.15 \cdot 10^{+33}:\\
\;\;\;\;x + \mathsf{fma}\left(0.10203362558171805, \frac{\frac{y}{z} \cdot \frac{t}{z}}{{t_1}^{2}}, \frac{y}{t_1}\right)\\
\mathbf{elif}\;z \leq 1.45 \cdot 10^{+23}:\\
\;\;\;\;x + \frac{y \cdot \left(b + z \cdot \left(a + z \cdot \left(t + z \cdot \left(11.1667541262 + z \cdot 3.13060547623\right)\right)\right)\right)}{0.607771387771 + z \cdot \left(11.9400905721 + z \cdot \left(31.4690115749 + z \cdot \left(z + 15.234687407\right)\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;x + \left(\left(\mathsf{fma}\left(y, 3.13060547623, \frac{y}{\frac{z \cdot z}{t}}\right) - \frac{y \cdot 36.52704169880642}{z}\right) - \mathsf{fma}\left(98.5170599679272, \frac{y}{z \cdot z}, \frac{\left(y \cdot 36.52704169880642\right) \cdot -15.234687407}{z \cdot z}\right)\right)\\
\end{array}
\end{array}
if z < -1.15000000000000005e33Initial program 2.2%
associate-/l*8.1%
fma-def8.1%
fma-def8.1%
fma-def8.1%
fma-def8.1%
fma-def8.1%
fma-def8.1%
fma-def8.1%
Simplified8.1%
Taylor expanded in z around inf 86.3%
associate-*r/86.3%
metadata-eval86.3%
mul-1-neg86.3%
*-commutative86.3%
unpow286.3%
Simplified86.3%
Taylor expanded in t around 0 90.9%
+-commutative90.9%
fma-def90.9%
Simplified98.5%
if -1.15000000000000005e33 < z < 1.45000000000000006e23Initial program 99.0%
if 1.45000000000000006e23 < z Initial program 6.9%
associate-*l/12.7%
*-commutative12.7%
fma-def12.7%
*-commutative12.7%
fma-def12.7%
*-commutative12.7%
fma-def12.7%
*-commutative12.7%
fma-def12.7%
Simplified12.7%
Taylor expanded in z around -inf 89.1%
+-commutative89.1%
mul-1-neg89.1%
unsub-neg89.1%
+-commutative89.1%
*-commutative89.1%
fma-def89.1%
associate-/l*98.3%
unpow298.3%
distribute-rgt-out--98.3%
metadata-eval98.3%
+-commutative98.3%
fma-def98.3%
Simplified98.3%
Final simplification98.7%
(FPCore (x y z t a b)
:precision binary64
(if (or (<= z -1.45e+33) (not (<= z 4.6e+22)))
(+
x
(-
(-
(fma y 3.13060547623 (/ y (/ (* z z) t)))
(/ (* y 36.52704169880642) z))
(fma
98.5170599679272
(/ y (* z z))
(/ (* (* y 36.52704169880642) -15.234687407) (* z z)))))
(+
x
(/
(*
y
(+
b
(* z (+ a (* z (+ t (* z (+ 11.1667541262 (* z 3.13060547623)))))))))
(+
0.607771387771
(*
z
(+ 11.9400905721 (* z (+ 31.4690115749 (* z (+ z 15.234687407)))))))))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if ((z <= -1.45e+33) || !(z <= 4.6e+22)) {
tmp = x + ((fma(y, 3.13060547623, (y / ((z * z) / t))) - ((y * 36.52704169880642) / z)) - fma(98.5170599679272, (y / (z * z)), (((y * 36.52704169880642) * -15.234687407) / (z * z))));
} else {
tmp = x + ((y * (b + (z * (a + (z * (t + (z * (11.1667541262 + (z * 3.13060547623))))))))) / (0.607771387771 + (z * (11.9400905721 + (z * (31.4690115749 + (z * (z + 15.234687407))))))));
}
return tmp;
}
function code(x, y, z, t, a, b) tmp = 0.0 if ((z <= -1.45e+33) || !(z <= 4.6e+22)) tmp = Float64(x + Float64(Float64(fma(y, 3.13060547623, Float64(y / Float64(Float64(z * z) / t))) - Float64(Float64(y * 36.52704169880642) / z)) - fma(98.5170599679272, Float64(y / Float64(z * z)), Float64(Float64(Float64(y * 36.52704169880642) * -15.234687407) / Float64(z * z))))); else tmp = Float64(x + Float64(Float64(y * Float64(b + Float64(z * Float64(a + Float64(z * Float64(t + Float64(z * Float64(11.1667541262 + Float64(z * 3.13060547623))))))))) / Float64(0.607771387771 + Float64(z * Float64(11.9400905721 + Float64(z * Float64(31.4690115749 + Float64(z * Float64(z + 15.234687407))))))))); end return tmp end
code[x_, y_, z_, t_, a_, b_] := If[Or[LessEqual[z, -1.45e+33], N[Not[LessEqual[z, 4.6e+22]], $MachinePrecision]], N[(x + N[(N[(N[(y * 3.13060547623 + N[(y / N[(N[(z * z), $MachinePrecision] / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - N[(N[(y * 36.52704169880642), $MachinePrecision] / z), $MachinePrecision]), $MachinePrecision] - N[(98.5170599679272 * N[(y / N[(z * z), $MachinePrecision]), $MachinePrecision] + N[(N[(N[(y * 36.52704169880642), $MachinePrecision] * -15.234687407), $MachinePrecision] / N[(z * z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x + N[(N[(y * N[(b + N[(z * N[(a + N[(z * N[(t + N[(z * N[(11.1667541262 + N[(z * 3.13060547623), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / 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]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;z \leq -1.45 \cdot 10^{+33} \lor \neg \left(z \leq 4.6 \cdot 10^{+22}\right):\\
\;\;\;\;x + \left(\left(\mathsf{fma}\left(y, 3.13060547623, \frac{y}{\frac{z \cdot z}{t}}\right) - \frac{y \cdot 36.52704169880642}{z}\right) - \mathsf{fma}\left(98.5170599679272, \frac{y}{z \cdot z}, \frac{\left(y \cdot 36.52704169880642\right) \cdot -15.234687407}{z \cdot z}\right)\right)\\
\mathbf{else}:\\
\;\;\;\;x + \frac{y \cdot \left(b + z \cdot \left(a + z \cdot \left(t + z \cdot \left(11.1667541262 + z \cdot 3.13060547623\right)\right)\right)\right)}{0.607771387771 + z \cdot \left(11.9400905721 + z \cdot \left(31.4690115749 + z \cdot \left(z + 15.234687407\right)\right)\right)}\\
\end{array}
\end{array}
if z < -1.45000000000000012e33 or 4.6000000000000004e22 < z Initial program 4.5%
associate-*l/10.4%
*-commutative10.4%
fma-def10.4%
*-commutative10.4%
fma-def10.4%
*-commutative10.4%
fma-def10.4%
*-commutative10.4%
fma-def10.4%
Simplified10.4%
Taylor expanded in z around -inf 90.1%
+-commutative90.1%
mul-1-neg90.1%
unsub-neg90.1%
+-commutative90.1%
*-commutative90.1%
fma-def90.1%
associate-/l*98.4%
unpow298.4%
distribute-rgt-out--98.4%
metadata-eval98.4%
+-commutative98.4%
fma-def98.4%
Simplified98.4%
if -1.45000000000000012e33 < z < 4.6000000000000004e22Initial program 99.0%
Final simplification98.7%
(FPCore (x y z t a b)
:precision binary64
(if (or (<= z -1.5e+33) (not (<= z 1.36e+23)))
(+
x
(*
y
(-
(- (/ (+ t 457.9610022158428) (* z z)) (/ 36.52704169880642 z))
-3.13060547623)))
(+
x
(/
(*
y
(+
b
(* z (+ a (* z (+ t (* z (+ 11.1667541262 (* z 3.13060547623)))))))))
(+
0.607771387771
(*
z
(+ 11.9400905721 (* z (+ 31.4690115749 (* z (+ z 15.234687407)))))))))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if ((z <= -1.5e+33) || !(z <= 1.36e+23)) {
tmp = x + (y * ((((t + 457.9610022158428) / (z * z)) - (36.52704169880642 / z)) - -3.13060547623));
} else {
tmp = x + ((y * (b + (z * (a + (z * (t + (z * (11.1667541262 + (z * 3.13060547623))))))))) / (0.607771387771 + (z * (11.9400905721 + (z * (31.4690115749 + (z * (z + 15.234687407))))))));
}
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 <= (-1.5d+33)) .or. (.not. (z <= 1.36d+23))) then
tmp = x + (y * ((((t + 457.9610022158428d0) / (z * z)) - (36.52704169880642d0 / z)) - (-3.13060547623d0)))
else
tmp = x + ((y * (b + (z * (a + (z * (t + (z * (11.1667541262d0 + (z * 3.13060547623d0))))))))) / (0.607771387771d0 + (z * (11.9400905721d0 + (z * (31.4690115749d0 + (z * (z + 15.234687407d0))))))))
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 <= -1.5e+33) || !(z <= 1.36e+23)) {
tmp = x + (y * ((((t + 457.9610022158428) / (z * z)) - (36.52704169880642 / z)) - -3.13060547623));
} else {
tmp = x + ((y * (b + (z * (a + (z * (t + (z * (11.1667541262 + (z * 3.13060547623))))))))) / (0.607771387771 + (z * (11.9400905721 + (z * (31.4690115749 + (z * (z + 15.234687407))))))));
}
return tmp;
}
def code(x, y, z, t, a, b): tmp = 0 if (z <= -1.5e+33) or not (z <= 1.36e+23): tmp = x + (y * ((((t + 457.9610022158428) / (z * z)) - (36.52704169880642 / z)) - -3.13060547623)) else: tmp = x + ((y * (b + (z * (a + (z * (t + (z * (11.1667541262 + (z * 3.13060547623))))))))) / (0.607771387771 + (z * (11.9400905721 + (z * (31.4690115749 + (z * (z + 15.234687407)))))))) return tmp
function code(x, y, z, t, a, b) tmp = 0.0 if ((z <= -1.5e+33) || !(z <= 1.36e+23)) tmp = Float64(x + Float64(y * Float64(Float64(Float64(Float64(t + 457.9610022158428) / Float64(z * z)) - Float64(36.52704169880642 / z)) - -3.13060547623))); else tmp = Float64(x + Float64(Float64(y * Float64(b + Float64(z * Float64(a + Float64(z * Float64(t + Float64(z * Float64(11.1667541262 + Float64(z * 3.13060547623))))))))) / Float64(0.607771387771 + Float64(z * Float64(11.9400905721 + Float64(z * Float64(31.4690115749 + Float64(z * Float64(z + 15.234687407))))))))); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) tmp = 0.0; if ((z <= -1.5e+33) || ~((z <= 1.36e+23))) tmp = x + (y * ((((t + 457.9610022158428) / (z * z)) - (36.52704169880642 / z)) - -3.13060547623)); else tmp = x + ((y * (b + (z * (a + (z * (t + (z * (11.1667541262 + (z * 3.13060547623))))))))) / (0.607771387771 + (z * (11.9400905721 + (z * (31.4690115749 + (z * (z + 15.234687407)))))))); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := If[Or[LessEqual[z, -1.5e+33], N[Not[LessEqual[z, 1.36e+23]], $MachinePrecision]], N[(x + N[(y * N[(N[(N[(N[(t + 457.9610022158428), $MachinePrecision] / N[(z * z), $MachinePrecision]), $MachinePrecision] - N[(36.52704169880642 / z), $MachinePrecision]), $MachinePrecision] - -3.13060547623), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x + N[(N[(y * N[(b + N[(z * N[(a + N[(z * N[(t + N[(z * N[(11.1667541262 + N[(z * 3.13060547623), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / 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]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;z \leq -1.5 \cdot 10^{+33} \lor \neg \left(z \leq 1.36 \cdot 10^{+23}\right):\\
\;\;\;\;x + y \cdot \left(\left(\frac{t + 457.9610022158428}{z \cdot z} - \frac{36.52704169880642}{z}\right) - -3.13060547623\right)\\
\mathbf{else}:\\
\;\;\;\;x + \frac{y \cdot \left(b + z \cdot \left(a + z \cdot \left(t + z \cdot \left(11.1667541262 + z \cdot 3.13060547623\right)\right)\right)\right)}{0.607771387771 + z \cdot \left(11.9400905721 + z \cdot \left(31.4690115749 + z \cdot \left(z + 15.234687407\right)\right)\right)}\\
\end{array}
\end{array}
if z < -1.49999999999999992e33 or 1.36e23 < z Initial program 4.5%
associate-/l*10.4%
fma-def10.4%
fma-def10.4%
fma-def10.4%
fma-def10.4%
fma-def10.4%
fma-def10.4%
fma-def10.4%
Simplified10.4%
Taylor expanded in z around inf 87.6%
associate-*r/87.6%
metadata-eval87.6%
mul-1-neg87.6%
*-commutative87.6%
unpow287.6%
Simplified87.6%
Taylor expanded in z around inf 90.0%
Taylor expanded in y around -inf 98.4%
Simplified98.4%
if -1.49999999999999992e33 < z < 1.36e23Initial program 99.0%
Final simplification98.7%
(FPCore (x y z t a b)
:precision binary64
(if (or (<= z -1.22e+33) (not (<= z 1.12e+23)))
(+
x
(*
y
(-
(- (/ (+ t 457.9610022158428) (* z z)) (/ 36.52704169880642 z))
-3.13060547623)))
(+
x
(/
(* y (+ b (* z (+ a (* z (+ t (* z 11.1667541262)))))))
(+
0.607771387771
(*
z
(+ 11.9400905721 (* z (+ 31.4690115749 (* z (+ z 15.234687407)))))))))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if ((z <= -1.22e+33) || !(z <= 1.12e+23)) {
tmp = x + (y * ((((t + 457.9610022158428) / (z * z)) - (36.52704169880642 / z)) - -3.13060547623));
} else {
tmp = x + ((y * (b + (z * (a + (z * (t + (z * 11.1667541262))))))) / (0.607771387771 + (z * (11.9400905721 + (z * (31.4690115749 + (z * (z + 15.234687407))))))));
}
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 <= (-1.22d+33)) .or. (.not. (z <= 1.12d+23))) then
tmp = x + (y * ((((t + 457.9610022158428d0) / (z * z)) - (36.52704169880642d0 / z)) - (-3.13060547623d0)))
else
tmp = x + ((y * (b + (z * (a + (z * (t + (z * 11.1667541262d0))))))) / (0.607771387771d0 + (z * (11.9400905721d0 + (z * (31.4690115749d0 + (z * (z + 15.234687407d0))))))))
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 <= -1.22e+33) || !(z <= 1.12e+23)) {
tmp = x + (y * ((((t + 457.9610022158428) / (z * z)) - (36.52704169880642 / z)) - -3.13060547623));
} else {
tmp = x + ((y * (b + (z * (a + (z * (t + (z * 11.1667541262))))))) / (0.607771387771 + (z * (11.9400905721 + (z * (31.4690115749 + (z * (z + 15.234687407))))))));
}
return tmp;
}
def code(x, y, z, t, a, b): tmp = 0 if (z <= -1.22e+33) or not (z <= 1.12e+23): tmp = x + (y * ((((t + 457.9610022158428) / (z * z)) - (36.52704169880642 / z)) - -3.13060547623)) else: tmp = x + ((y * (b + (z * (a + (z * (t + (z * 11.1667541262))))))) / (0.607771387771 + (z * (11.9400905721 + (z * (31.4690115749 + (z * (z + 15.234687407)))))))) return tmp
function code(x, y, z, t, a, b) tmp = 0.0 if ((z <= -1.22e+33) || !(z <= 1.12e+23)) tmp = Float64(x + Float64(y * Float64(Float64(Float64(Float64(t + 457.9610022158428) / Float64(z * z)) - Float64(36.52704169880642 / z)) - -3.13060547623))); else 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 * Float64(11.9400905721 + Float64(z * Float64(31.4690115749 + Float64(z * Float64(z + 15.234687407))))))))); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) tmp = 0.0; if ((z <= -1.22e+33) || ~((z <= 1.12e+23))) tmp = x + (y * ((((t + 457.9610022158428) / (z * z)) - (36.52704169880642 / z)) - -3.13060547623)); else tmp = x + ((y * (b + (z * (a + (z * (t + (z * 11.1667541262))))))) / (0.607771387771 + (z * (11.9400905721 + (z * (31.4690115749 + (z * (z + 15.234687407)))))))); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := If[Or[LessEqual[z, -1.22e+33], N[Not[LessEqual[z, 1.12e+23]], $MachinePrecision]], N[(x + N[(y * N[(N[(N[(N[(t + 457.9610022158428), $MachinePrecision] / N[(z * z), $MachinePrecision]), $MachinePrecision] - N[(36.52704169880642 / z), $MachinePrecision]), $MachinePrecision] - -3.13060547623), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 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 * N[(11.9400905721 + N[(z * N[(31.4690115749 + N[(z * N[(z + 15.234687407), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;z \leq -1.22 \cdot 10^{+33} \lor \neg \left(z \leq 1.12 \cdot 10^{+23}\right):\\
\;\;\;\;x + y \cdot \left(\left(\frac{t + 457.9610022158428}{z \cdot z} - \frac{36.52704169880642}{z}\right) - -3.13060547623\right)\\
\mathbf{else}:\\
\;\;\;\;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 \left(11.9400905721 + z \cdot \left(31.4690115749 + z \cdot \left(z + 15.234687407\right)\right)\right)}\\
\end{array}
\end{array}
if z < -1.22000000000000005e33 or 1.12e23 < z Initial program 4.5%
associate-/l*10.4%
fma-def10.4%
fma-def10.4%
fma-def10.4%
fma-def10.4%
fma-def10.4%
fma-def10.4%
fma-def10.4%
Simplified10.4%
Taylor expanded in z around inf 87.6%
associate-*r/87.6%
metadata-eval87.6%
mul-1-neg87.6%
*-commutative87.6%
unpow287.6%
Simplified87.6%
Taylor expanded in z around inf 90.0%
Taylor expanded in y around -inf 98.4%
Simplified98.4%
if -1.22000000000000005e33 < z < 1.12e23Initial program 99.0%
Taylor expanded in z around 0 98.4%
*-commutative95.8%
Simplified98.4%
Final simplification98.4%
(FPCore (x y z t a b)
:precision binary64
(if (or (<= z -1.18e+33) (not (<= z 2.9e+22)))
(+
x
(*
y
(-
(- (/ (+ t 457.9610022158428) (* z z)) (/ 36.52704169880642 z))
-3.13060547623)))
(+
x
(/
(* y (+ b (* z (+ a (* z (+ t (* z 11.1667541262)))))))
(+ 0.607771387771 (* z (+ 11.9400905721 (* z 31.4690115749))))))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if ((z <= -1.18e+33) || !(z <= 2.9e+22)) {
tmp = x + (y * ((((t + 457.9610022158428) / (z * z)) - (36.52704169880642 / z)) - -3.13060547623));
} else {
tmp = x + ((y * (b + (z * (a + (z * (t + (z * 11.1667541262))))))) / (0.607771387771 + (z * (11.9400905721 + (z * 31.4690115749)))));
}
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 <= (-1.18d+33)) .or. (.not. (z <= 2.9d+22))) then
tmp = x + (y * ((((t + 457.9610022158428d0) / (z * z)) - (36.52704169880642d0 / z)) - (-3.13060547623d0)))
else
tmp = x + ((y * (b + (z * (a + (z * (t + (z * 11.1667541262d0))))))) / (0.607771387771d0 + (z * (11.9400905721d0 + (z * 31.4690115749d0)))))
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 <= -1.18e+33) || !(z <= 2.9e+22)) {
tmp = x + (y * ((((t + 457.9610022158428) / (z * z)) - (36.52704169880642 / z)) - -3.13060547623));
} else {
tmp = x + ((y * (b + (z * (a + (z * (t + (z * 11.1667541262))))))) / (0.607771387771 + (z * (11.9400905721 + (z * 31.4690115749)))));
}
return tmp;
}
def code(x, y, z, t, a, b): tmp = 0 if (z <= -1.18e+33) or not (z <= 2.9e+22): tmp = x + (y * ((((t + 457.9610022158428) / (z * z)) - (36.52704169880642 / z)) - -3.13060547623)) else: tmp = x + ((y * (b + (z * (a + (z * (t + (z * 11.1667541262))))))) / (0.607771387771 + (z * (11.9400905721 + (z * 31.4690115749))))) return tmp
function code(x, y, z, t, a, b) tmp = 0.0 if ((z <= -1.18e+33) || !(z <= 2.9e+22)) tmp = Float64(x + Float64(y * Float64(Float64(Float64(Float64(t + 457.9610022158428) / Float64(z * z)) - Float64(36.52704169880642 / z)) - -3.13060547623))); else 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 * Float64(11.9400905721 + Float64(z * 31.4690115749)))))); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) tmp = 0.0; if ((z <= -1.18e+33) || ~((z <= 2.9e+22))) tmp = x + (y * ((((t + 457.9610022158428) / (z * z)) - (36.52704169880642 / z)) - -3.13060547623)); else tmp = x + ((y * (b + (z * (a + (z * (t + (z * 11.1667541262))))))) / (0.607771387771 + (z * (11.9400905721 + (z * 31.4690115749))))); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := If[Or[LessEqual[z, -1.18e+33], N[Not[LessEqual[z, 2.9e+22]], $MachinePrecision]], N[(x + N[(y * N[(N[(N[(N[(t + 457.9610022158428), $MachinePrecision] / N[(z * z), $MachinePrecision]), $MachinePrecision] - N[(36.52704169880642 / z), $MachinePrecision]), $MachinePrecision] - -3.13060547623), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 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 * N[(11.9400905721 + N[(z * 31.4690115749), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;z \leq -1.18 \cdot 10^{+33} \lor \neg \left(z \leq 2.9 \cdot 10^{+22}\right):\\
\;\;\;\;x + y \cdot \left(\left(\frac{t + 457.9610022158428}{z \cdot z} - \frac{36.52704169880642}{z}\right) - -3.13060547623\right)\\
\mathbf{else}:\\
\;\;\;\;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 \left(11.9400905721 + z \cdot 31.4690115749\right)}\\
\end{array}
\end{array}
if z < -1.17999999999999993e33 or 2.9e22 < z Initial program 4.5%
associate-/l*10.4%
fma-def10.4%
fma-def10.4%
fma-def10.4%
fma-def10.4%
fma-def10.4%
fma-def10.4%
fma-def10.4%
Simplified10.4%
Taylor expanded in z around inf 87.6%
associate-*r/87.6%
metadata-eval87.6%
mul-1-neg87.6%
*-commutative87.6%
unpow287.6%
Simplified87.6%
Taylor expanded in z around inf 90.0%
Taylor expanded in y around -inf 98.4%
Simplified98.4%
if -1.17999999999999993e33 < z < 2.9e22Initial program 99.0%
Taylor expanded in z around 0 95.8%
*-commutative95.8%
Simplified95.8%
Taylor expanded in z around 0 95.8%
*-commutative95.8%
Simplified95.8%
Final simplification97.1%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (+ x (/ y 0.31942702700572795))))
(if (<= z -3.5e+125)
t_1
(if (<= z -1380.0)
(+ x (* (/ y z) (/ t z)))
(if (<= z 2.2e+23)
(+
x
(*
y
(+
(* z (- (* a 1.6453555072203998) (* b 32.324150453290734)))
(* b 1.6453555072203998))))
t_1)))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = x + (y / 0.31942702700572795);
double tmp;
if (z <= -3.5e+125) {
tmp = t_1;
} else if (z <= -1380.0) {
tmp = x + ((y / z) * (t / z));
} else if (z <= 2.2e+23) {
tmp = x + (y * ((z * ((a * 1.6453555072203998) - (b * 32.324150453290734))) + (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 / 0.31942702700572795d0)
if (z <= (-3.5d+125)) then
tmp = t_1
else if (z <= (-1380.0d0)) then
tmp = x + ((y / z) * (t / z))
else if (z <= 2.2d+23) then
tmp = x + (y * ((z * ((a * 1.6453555072203998d0) - (b * 32.324150453290734d0))) + (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 / 0.31942702700572795);
double tmp;
if (z <= -3.5e+125) {
tmp = t_1;
} else if (z <= -1380.0) {
tmp = x + ((y / z) * (t / z));
} else if (z <= 2.2e+23) {
tmp = x + (y * ((z * ((a * 1.6453555072203998) - (b * 32.324150453290734))) + (b * 1.6453555072203998)));
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = x + (y / 0.31942702700572795) tmp = 0 if z <= -3.5e+125: tmp = t_1 elif z <= -1380.0: tmp = x + ((y / z) * (t / z)) elif z <= 2.2e+23: tmp = x + (y * ((z * ((a * 1.6453555072203998) - (b * 32.324150453290734))) + (b * 1.6453555072203998))) else: tmp = t_1 return tmp
function code(x, y, z, t, a, b) t_1 = Float64(x + Float64(y / 0.31942702700572795)) tmp = 0.0 if (z <= -3.5e+125) tmp = t_1; elseif (z <= -1380.0) tmp = Float64(x + Float64(Float64(y / z) * Float64(t / z))); elseif (z <= 2.2e+23) tmp = Float64(x + Float64(y * Float64(Float64(z * Float64(Float64(a * 1.6453555072203998) - Float64(b * 32.324150453290734))) + 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 / 0.31942702700572795); tmp = 0.0; if (z <= -3.5e+125) tmp = t_1; elseif (z <= -1380.0) tmp = x + ((y / z) * (t / z)); elseif (z <= 2.2e+23) tmp = x + (y * ((z * ((a * 1.6453555072203998) - (b * 32.324150453290734))) + (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 / 0.31942702700572795), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[z, -3.5e+125], t$95$1, If[LessEqual[z, -1380.0], N[(x + N[(N[(y / z), $MachinePrecision] * N[(t / z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[z, 2.2e+23], N[(x + N[(y * N[(N[(z * N[(N[(a * 1.6453555072203998), $MachinePrecision] - N[(b * 32.324150453290734), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(b * 1.6453555072203998), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$1]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x + \frac{y}{0.31942702700572795}\\
\mathbf{if}\;z \leq -3.5 \cdot 10^{+125}:\\
\;\;\;\;t_1\\
\mathbf{elif}\;z \leq -1380:\\
\;\;\;\;x + \frac{y}{z} \cdot \frac{t}{z}\\
\mathbf{elif}\;z \leq 2.2 \cdot 10^{+23}:\\
\;\;\;\;x + y \cdot \left(z \cdot \left(a \cdot 1.6453555072203998 - b \cdot 32.324150453290734\right) + b \cdot 1.6453555072203998\right)\\
\mathbf{else}:\\
\;\;\;\;t_1\\
\end{array}
\end{array}
if z < -3.50000000000000011e125 or 2.20000000000000008e23 < z Initial program 3.2%
associate-/l*6.7%
fma-def6.7%
fma-def6.7%
fma-def6.7%
fma-def6.7%
fma-def6.7%
fma-def6.7%
fma-def6.7%
Simplified6.7%
Taylor expanded in z around inf 94.0%
if -3.50000000000000011e125 < z < -1380Initial program 21.6%
associate-/l*36.7%
fma-def36.7%
fma-def36.7%
fma-def36.7%
fma-def36.7%
fma-def36.7%
fma-def36.7%
fma-def36.7%
Simplified36.7%
Taylor expanded in z around inf 61.1%
associate-*r/61.1%
metadata-eval61.1%
mul-1-neg61.1%
*-commutative61.1%
unpow261.1%
Simplified61.1%
Taylor expanded in z around inf 88.7%
Taylor expanded in t around inf 80.1%
unpow280.1%
times-frac87.5%
Simplified87.5%
if -1380 < z < 2.20000000000000008e23Initial program 99.0%
associate-*l/99.6%
*-commutative99.6%
fma-def99.6%
*-commutative99.6%
fma-def99.6%
*-commutative99.6%
fma-def99.6%
*-commutative99.6%
fma-def99.6%
Simplified99.6%
Taylor expanded in z around 0 76.8%
Taylor expanded in y around 0 89.2%
Final simplification91.0%
(FPCore (x y z t a b)
:precision binary64
(if (<= z -3.5e+125)
(+ x (/ y 0.31942702700572795))
(if (<= z -65.0)
(+ x (* (/ y z) (/ t z)))
(if (<= z 2.6e+23)
(+
x
(*
y
(+
(* z (- (* a 1.6453555072203998) (* b 32.324150453290734)))
(* b 1.6453555072203998))))
(+
x
(/
y
(+
(/ 3.7269864963038164 z)
(+ 0.31942702700572795 (/ (* t -0.10203362558171805) (* z z))))))))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (z <= -3.5e+125) {
tmp = x + (y / 0.31942702700572795);
} else if (z <= -65.0) {
tmp = x + ((y / z) * (t / z));
} else if (z <= 2.6e+23) {
tmp = x + (y * ((z * ((a * 1.6453555072203998) - (b * 32.324150453290734))) + (b * 1.6453555072203998)));
} else {
tmp = x + (y / ((3.7269864963038164 / z) + (0.31942702700572795 + ((t * -0.10203362558171805) / (z * z)))));
}
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 <= (-3.5d+125)) then
tmp = x + (y / 0.31942702700572795d0)
else if (z <= (-65.0d0)) then
tmp = x + ((y / z) * (t / z))
else if (z <= 2.6d+23) then
tmp = x + (y * ((z * ((a * 1.6453555072203998d0) - (b * 32.324150453290734d0))) + (b * 1.6453555072203998d0)))
else
tmp = x + (y / ((3.7269864963038164d0 / z) + (0.31942702700572795d0 + ((t * (-0.10203362558171805d0)) / (z * z)))))
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 <= -3.5e+125) {
tmp = x + (y / 0.31942702700572795);
} else if (z <= -65.0) {
tmp = x + ((y / z) * (t / z));
} else if (z <= 2.6e+23) {
tmp = x + (y * ((z * ((a * 1.6453555072203998) - (b * 32.324150453290734))) + (b * 1.6453555072203998)));
} else {
tmp = x + (y / ((3.7269864963038164 / z) + (0.31942702700572795 + ((t * -0.10203362558171805) / (z * z)))));
}
return tmp;
}
def code(x, y, z, t, a, b): tmp = 0 if z <= -3.5e+125: tmp = x + (y / 0.31942702700572795) elif z <= -65.0: tmp = x + ((y / z) * (t / z)) elif z <= 2.6e+23: tmp = x + (y * ((z * ((a * 1.6453555072203998) - (b * 32.324150453290734))) + (b * 1.6453555072203998))) else: tmp = x + (y / ((3.7269864963038164 / z) + (0.31942702700572795 + ((t * -0.10203362558171805) / (z * z))))) return tmp
function code(x, y, z, t, a, b) tmp = 0.0 if (z <= -3.5e+125) tmp = Float64(x + Float64(y / 0.31942702700572795)); elseif (z <= -65.0) tmp = Float64(x + Float64(Float64(y / z) * Float64(t / z))); elseif (z <= 2.6e+23) tmp = Float64(x + Float64(y * Float64(Float64(z * Float64(Float64(a * 1.6453555072203998) - Float64(b * 32.324150453290734))) + Float64(b * 1.6453555072203998)))); else tmp = Float64(x + Float64(y / Float64(Float64(3.7269864963038164 / z) + Float64(0.31942702700572795 + Float64(Float64(t * -0.10203362558171805) / Float64(z * z)))))); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) tmp = 0.0; if (z <= -3.5e+125) tmp = x + (y / 0.31942702700572795); elseif (z <= -65.0) tmp = x + ((y / z) * (t / z)); elseif (z <= 2.6e+23) tmp = x + (y * ((z * ((a * 1.6453555072203998) - (b * 32.324150453290734))) + (b * 1.6453555072203998))); else tmp = x + (y / ((3.7269864963038164 / z) + (0.31942702700572795 + ((t * -0.10203362558171805) / (z * z))))); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[z, -3.5e+125], N[(x + N[(y / 0.31942702700572795), $MachinePrecision]), $MachinePrecision], If[LessEqual[z, -65.0], N[(x + N[(N[(y / z), $MachinePrecision] * N[(t / z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[z, 2.6e+23], N[(x + N[(y * N[(N[(z * N[(N[(a * 1.6453555072203998), $MachinePrecision] - N[(b * 32.324150453290734), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(b * 1.6453555072203998), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x + N[(y / N[(N[(3.7269864963038164 / z), $MachinePrecision] + N[(0.31942702700572795 + N[(N[(t * -0.10203362558171805), $MachinePrecision] / N[(z * z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;z \leq -3.5 \cdot 10^{+125}:\\
\;\;\;\;x + \frac{y}{0.31942702700572795}\\
\mathbf{elif}\;z \leq -65:\\
\;\;\;\;x + \frac{y}{z} \cdot \frac{t}{z}\\
\mathbf{elif}\;z \leq 2.6 \cdot 10^{+23}:\\
\;\;\;\;x + y \cdot \left(z \cdot \left(a \cdot 1.6453555072203998 - b \cdot 32.324150453290734\right) + b \cdot 1.6453555072203998\right)\\
\mathbf{else}:\\
\;\;\;\;x + \frac{y}{\frac{3.7269864963038164}{z} + \left(0.31942702700572795 + \frac{t \cdot -0.10203362558171805}{z \cdot z}\right)}\\
\end{array}
\end{array}
if z < -3.50000000000000011e125Initial program 0.0%
associate-/l*0.0%
fma-def0.0%
fma-def0.0%
fma-def0.0%
fma-def0.0%
fma-def0.0%
fma-def0.0%
fma-def0.0%
Simplified0.0%
Taylor expanded in z around inf 100.0%
if -3.50000000000000011e125 < z < -65Initial program 21.6%
associate-/l*36.7%
fma-def36.7%
fma-def36.7%
fma-def36.7%
fma-def36.7%
fma-def36.7%
fma-def36.7%
fma-def36.7%
Simplified36.7%
Taylor expanded in z around inf 61.1%
associate-*r/61.1%
metadata-eval61.1%
mul-1-neg61.1%
*-commutative61.1%
unpow261.1%
Simplified61.1%
Taylor expanded in z around inf 88.7%
Taylor expanded in t around inf 80.1%
unpow280.1%
times-frac87.5%
Simplified87.5%
if -65 < z < 2.59999999999999992e23Initial program 99.0%
associate-*l/99.6%
*-commutative99.6%
fma-def99.6%
*-commutative99.6%
fma-def99.6%
*-commutative99.6%
fma-def99.6%
*-commutative99.6%
fma-def99.6%
Simplified99.6%
Taylor expanded in z around 0 76.8%
Taylor expanded in y around 0 89.2%
if 2.59999999999999992e23 < z Initial program 5.4%
associate-/l*11.3%
fma-def11.3%
fma-def11.3%
fma-def11.3%
fma-def11.3%
fma-def11.3%
fma-def11.3%
fma-def11.3%
Simplified11.3%
Taylor expanded in z around inf 90.2%
associate-*r/90.2%
metadata-eval90.2%
mul-1-neg90.2%
*-commutative90.2%
unpow290.2%
Simplified90.2%
Taylor expanded in t around inf 90.2%
associate-*r/90.2%
unpow290.2%
Simplified90.2%
Taylor expanded in y around 0 90.2%
associate--l+90.2%
associate-*r/90.2%
metadata-eval90.2%
cancel-sign-sub-inv90.2%
associate-*r/90.2%
metadata-eval90.2%
unpow290.2%
Simplified90.2%
Final simplification91.1%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (+ x (/ y 0.31942702700572795))))
(if (<= z -3.5e+125)
t_1
(if (<= z -30000.0)
(+ x (* (/ y z) (/ t z)))
(if (<= z 4.2e+15)
(+
x
(+
(* 1.6453555072203998 (* y (* z a)))
(* 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 / 0.31942702700572795);
double tmp;
if (z <= -3.5e+125) {
tmp = t_1;
} else if (z <= -30000.0) {
tmp = x + ((y / z) * (t / z));
} else if (z <= 4.2e+15) {
tmp = x + ((1.6453555072203998 * (y * (z * a))) + (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 / 0.31942702700572795d0)
if (z <= (-3.5d+125)) then
tmp = t_1
else if (z <= (-30000.0d0)) then
tmp = x + ((y / z) * (t / z))
else if (z <= 4.2d+15) then
tmp = x + ((1.6453555072203998d0 * (y * (z * a))) + (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 / 0.31942702700572795);
double tmp;
if (z <= -3.5e+125) {
tmp = t_1;
} else if (z <= -30000.0) {
tmp = x + ((y / z) * (t / z));
} else if (z <= 4.2e+15) {
tmp = x + ((1.6453555072203998 * (y * (z * a))) + (1.6453555072203998 * (y * b)));
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = x + (y / 0.31942702700572795) tmp = 0 if z <= -3.5e+125: tmp = t_1 elif z <= -30000.0: tmp = x + ((y / z) * (t / z)) elif z <= 4.2e+15: tmp = x + ((1.6453555072203998 * (y * (z * a))) + (1.6453555072203998 * (y * b))) else: tmp = t_1 return tmp
function code(x, y, z, t, a, b) t_1 = Float64(x + Float64(y / 0.31942702700572795)) tmp = 0.0 if (z <= -3.5e+125) tmp = t_1; elseif (z <= -30000.0) tmp = Float64(x + Float64(Float64(y / z) * Float64(t / z))); elseif (z <= 4.2e+15) tmp = Float64(x + Float64(Float64(1.6453555072203998 * Float64(y * Float64(z * a))) + 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 / 0.31942702700572795); tmp = 0.0; if (z <= -3.5e+125) tmp = t_1; elseif (z <= -30000.0) tmp = x + ((y / z) * (t / z)); elseif (z <= 4.2e+15) tmp = x + ((1.6453555072203998 * (y * (z * a))) + (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 / 0.31942702700572795), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[z, -3.5e+125], t$95$1, If[LessEqual[z, -30000.0], N[(x + N[(N[(y / z), $MachinePrecision] * N[(t / z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[z, 4.2e+15], N[(x + N[(N[(1.6453555072203998 * N[(y * N[(z * a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(1.6453555072203998 * N[(y * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$1]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x + \frac{y}{0.31942702700572795}\\
\mathbf{if}\;z \leq -3.5 \cdot 10^{+125}:\\
\;\;\;\;t_1\\
\mathbf{elif}\;z \leq -30000:\\
\;\;\;\;x + \frac{y}{z} \cdot \frac{t}{z}\\
\mathbf{elif}\;z \leq 4.2 \cdot 10^{+15}:\\
\;\;\;\;x + \left(1.6453555072203998 \cdot \left(y \cdot \left(z \cdot a\right)\right) + 1.6453555072203998 \cdot \left(y \cdot b\right)\right)\\
\mathbf{else}:\\
\;\;\;\;t_1\\
\end{array}
\end{array}
if z < -3.50000000000000011e125 or 4.2e15 < z Initial program 6.8%
associate-/l*10.9%
fma-def10.9%
fma-def10.9%
fma-def10.9%
fma-def10.9%
fma-def10.9%
fma-def10.9%
fma-def10.9%
Simplified10.9%
Taylor expanded in z around inf 91.0%
if -3.50000000000000011e125 < z < -3e4Initial program 21.6%
associate-/l*36.7%
fma-def36.7%
fma-def36.7%
fma-def36.7%
fma-def36.7%
fma-def36.7%
fma-def36.7%
fma-def36.7%
Simplified36.7%
Taylor expanded in z around inf 61.1%
associate-*r/61.1%
metadata-eval61.1%
mul-1-neg61.1%
*-commutative61.1%
unpow261.1%
Simplified61.1%
Taylor expanded in z around inf 88.7%
Taylor expanded in t around inf 80.1%
unpow280.1%
times-frac87.5%
Simplified87.5%
if -3e4 < z < 4.2e15Initial program 99.7%
associate-*l/99.6%
*-commutative99.6%
fma-def99.6%
*-commutative99.6%
fma-def99.7%
*-commutative99.7%
fma-def99.6%
*-commutative99.6%
fma-def99.6%
Simplified99.6%
Taylor expanded in z around 0 79.1%
Taylor expanded in a around inf 90.2%
*-commutative90.2%
associate-*r*90.2%
Simplified90.2%
Final simplification90.3%
(FPCore (x y z t a b)
:precision binary64
(if (or (<= z -13200.0) (not (<= z 0.6)))
(+
x
(*
y
(-
(- (/ (+ t 457.9610022158428) (* z z)) (/ 36.52704169880642 z))
-3.13060547623)))
(+
x
(*
y
(+
(* z (- (* a 1.6453555072203998) (* b 32.324150453290734)))
(* b 1.6453555072203998))))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if ((z <= -13200.0) || !(z <= 0.6)) {
tmp = x + (y * ((((t + 457.9610022158428) / (z * z)) - (36.52704169880642 / z)) - -3.13060547623));
} else {
tmp = x + (y * ((z * ((a * 1.6453555072203998) - (b * 32.324150453290734))) + (b * 1.6453555072203998)));
}
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 <= (-13200.0d0)) .or. (.not. (z <= 0.6d0))) then
tmp = x + (y * ((((t + 457.9610022158428d0) / (z * z)) - (36.52704169880642d0 / z)) - (-3.13060547623d0)))
else
tmp = x + (y * ((z * ((a * 1.6453555072203998d0) - (b * 32.324150453290734d0))) + (b * 1.6453555072203998d0)))
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 <= -13200.0) || !(z <= 0.6)) {
tmp = x + (y * ((((t + 457.9610022158428) / (z * z)) - (36.52704169880642 / z)) - -3.13060547623));
} else {
tmp = x + (y * ((z * ((a * 1.6453555072203998) - (b * 32.324150453290734))) + (b * 1.6453555072203998)));
}
return tmp;
}
def code(x, y, z, t, a, b): tmp = 0 if (z <= -13200.0) or not (z <= 0.6): tmp = x + (y * ((((t + 457.9610022158428) / (z * z)) - (36.52704169880642 / z)) - -3.13060547623)) else: tmp = x + (y * ((z * ((a * 1.6453555072203998) - (b * 32.324150453290734))) + (b * 1.6453555072203998))) return tmp
function code(x, y, z, t, a, b) tmp = 0.0 if ((z <= -13200.0) || !(z <= 0.6)) tmp = Float64(x + Float64(y * Float64(Float64(Float64(Float64(t + 457.9610022158428) / Float64(z * z)) - Float64(36.52704169880642 / z)) - -3.13060547623))); else tmp = Float64(x + Float64(y * Float64(Float64(z * Float64(Float64(a * 1.6453555072203998) - Float64(b * 32.324150453290734))) + Float64(b * 1.6453555072203998)))); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) tmp = 0.0; if ((z <= -13200.0) || ~((z <= 0.6))) tmp = x + (y * ((((t + 457.9610022158428) / (z * z)) - (36.52704169880642 / z)) - -3.13060547623)); else tmp = x + (y * ((z * ((a * 1.6453555072203998) - (b * 32.324150453290734))) + (b * 1.6453555072203998))); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := If[Or[LessEqual[z, -13200.0], N[Not[LessEqual[z, 0.6]], $MachinePrecision]], N[(x + N[(y * N[(N[(N[(N[(t + 457.9610022158428), $MachinePrecision] / N[(z * z), $MachinePrecision]), $MachinePrecision] - N[(36.52704169880642 / z), $MachinePrecision]), $MachinePrecision] - -3.13060547623), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x + N[(y * N[(N[(z * N[(N[(a * 1.6453555072203998), $MachinePrecision] - N[(b * 32.324150453290734), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(b * 1.6453555072203998), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;z \leq -13200 \lor \neg \left(z \leq 0.6\right):\\
\;\;\;\;x + y \cdot \left(\left(\frac{t + 457.9610022158428}{z \cdot z} - \frac{36.52704169880642}{z}\right) - -3.13060547623\right)\\
\mathbf{else}:\\
\;\;\;\;x + y \cdot \left(z \cdot \left(a \cdot 1.6453555072203998 - b \cdot 32.324150453290734\right) + b \cdot 1.6453555072203998\right)\\
\end{array}
\end{array}
if z < -13200 or 0.599999999999999978 < z Initial program 10.2%
associate-/l*16.3%
fma-def16.3%
fma-def16.3%
fma-def16.3%
fma-def16.3%
fma-def16.3%
fma-def16.3%
fma-def16.3%
Simplified16.3%
Taylor expanded in z around inf 85.1%
associate-*r/85.1%
metadata-eval85.1%
mul-1-neg85.1%
*-commutative85.1%
unpow285.1%
Simplified85.1%
Taylor expanded in z around inf 88.1%
Taylor expanded in y around -inf 95.9%
Simplified95.9%
if -13200 < z < 0.599999999999999978Initial program 99.7%
associate-*l/99.7%
*-commutative99.7%
fma-def99.7%
*-commutative99.7%
fma-def99.7%
*-commutative99.7%
fma-def99.6%
*-commutative99.6%
fma-def99.6%
Simplified99.6%
Taylor expanded in z around 0 79.8%
Taylor expanded in y around 0 91.9%
Final simplification94.0%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (+ x (/ y 0.31942702700572795))))
(if (<= z -3.8e+125)
t_1
(if (<= z -170000.0)
(+ x (* (/ y z) (/ t z)))
(if (<= z 3.4e+22)
(+
x
(* b (+ (* -32.324150453290734 (* z y)) (* y 1.6453555072203998))))
t_1)))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = x + (y / 0.31942702700572795);
double tmp;
if (z <= -3.8e+125) {
tmp = t_1;
} else if (z <= -170000.0) {
tmp = x + ((y / z) * (t / z));
} else if (z <= 3.4e+22) {
tmp = x + (b * ((-32.324150453290734 * (z * y)) + (y * 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 / 0.31942702700572795d0)
if (z <= (-3.8d+125)) then
tmp = t_1
else if (z <= (-170000.0d0)) then
tmp = x + ((y / z) * (t / z))
else if (z <= 3.4d+22) then
tmp = x + (b * (((-32.324150453290734d0) * (z * y)) + (y * 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 / 0.31942702700572795);
double tmp;
if (z <= -3.8e+125) {
tmp = t_1;
} else if (z <= -170000.0) {
tmp = x + ((y / z) * (t / z));
} else if (z <= 3.4e+22) {
tmp = x + (b * ((-32.324150453290734 * (z * y)) + (y * 1.6453555072203998)));
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = x + (y / 0.31942702700572795) tmp = 0 if z <= -3.8e+125: tmp = t_1 elif z <= -170000.0: tmp = x + ((y / z) * (t / z)) elif z <= 3.4e+22: tmp = x + (b * ((-32.324150453290734 * (z * y)) + (y * 1.6453555072203998))) else: tmp = t_1 return tmp
function code(x, y, z, t, a, b) t_1 = Float64(x + Float64(y / 0.31942702700572795)) tmp = 0.0 if (z <= -3.8e+125) tmp = t_1; elseif (z <= -170000.0) tmp = Float64(x + Float64(Float64(y / z) * Float64(t / z))); elseif (z <= 3.4e+22) tmp = Float64(x + Float64(b * Float64(Float64(-32.324150453290734 * Float64(z * y)) + Float64(y * 1.6453555072203998)))); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a, b) t_1 = x + (y / 0.31942702700572795); tmp = 0.0; if (z <= -3.8e+125) tmp = t_1; elseif (z <= -170000.0) tmp = x + ((y / z) * (t / z)); elseif (z <= 3.4e+22) tmp = x + (b * ((-32.324150453290734 * (z * y)) + (y * 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 / 0.31942702700572795), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[z, -3.8e+125], t$95$1, If[LessEqual[z, -170000.0], N[(x + N[(N[(y / z), $MachinePrecision] * N[(t / z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[z, 3.4e+22], N[(x + N[(b * N[(N[(-32.324150453290734 * N[(z * y), $MachinePrecision]), $MachinePrecision] + N[(y * 1.6453555072203998), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$1]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x + \frac{y}{0.31942702700572795}\\
\mathbf{if}\;z \leq -3.8 \cdot 10^{+125}:\\
\;\;\;\;t_1\\
\mathbf{elif}\;z \leq -170000:\\
\;\;\;\;x + \frac{y}{z} \cdot \frac{t}{z}\\
\mathbf{elif}\;z \leq 3.4 \cdot 10^{+22}:\\
\;\;\;\;x + b \cdot \left(-32.324150453290734 \cdot \left(z \cdot y\right) + y \cdot 1.6453555072203998\right)\\
\mathbf{else}:\\
\;\;\;\;t_1\\
\end{array}
\end{array}
if z < -3.80000000000000002e125 or 3.4e22 < z Initial program 4.1%
associate-/l*7.6%
fma-def7.6%
fma-def7.6%
fma-def7.6%
fma-def7.6%
fma-def7.6%
fma-def7.6%
fma-def7.6%
Simplified7.6%
Taylor expanded in z around inf 93.2%
if -3.80000000000000002e125 < z < -1.7e5Initial program 21.6%
associate-/l*36.7%
fma-def36.7%
fma-def36.7%
fma-def36.7%
fma-def36.7%
fma-def36.7%
fma-def36.7%
fma-def36.7%
Simplified36.7%
Taylor expanded in z around inf 61.1%
associate-*r/61.1%
metadata-eval61.1%
mul-1-neg61.1%
*-commutative61.1%
unpow261.1%
Simplified61.1%
Taylor expanded in z around inf 88.7%
Taylor expanded in t around inf 80.1%
unpow280.1%
times-frac87.5%
Simplified87.5%
if -1.7e5 < z < 3.4e22Initial program 99.0%
associate-*l/99.6%
*-commutative99.6%
fma-def99.6%
*-commutative99.6%
fma-def99.6%
*-commutative99.6%
fma-def99.6%
*-commutative99.6%
fma-def99.6%
Simplified99.6%
Taylor expanded in z around 0 76.6%
Taylor expanded in b around inf 78.5%
Final simplification85.4%
(FPCore (x y z t a b)
:precision binary64
(if (or (<= y -2e+22)
(and (not (<= y 3.8e+29))
(or (<= y 1.65e+125) (not (<= y 5.8e+189)))))
(* 1.6453555072203998 (* y b))
x))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if ((y <= -2e+22) || (!(y <= 3.8e+29) && ((y <= 1.65e+125) || !(y <= 5.8e+189)))) {
tmp = 1.6453555072203998 * (y * b);
} else {
tmp = x;
}
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 <= (-2d+22)) .or. (.not. (y <= 3.8d+29)) .and. (y <= 1.65d+125) .or. (.not. (y <= 5.8d+189))) then
tmp = 1.6453555072203998d0 * (y * b)
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 tmp;
if ((y <= -2e+22) || (!(y <= 3.8e+29) && ((y <= 1.65e+125) || !(y <= 5.8e+189)))) {
tmp = 1.6453555072203998 * (y * b);
} else {
tmp = x;
}
return tmp;
}
def code(x, y, z, t, a, b): tmp = 0 if (y <= -2e+22) or (not (y <= 3.8e+29) and ((y <= 1.65e+125) or not (y <= 5.8e+189))): tmp = 1.6453555072203998 * (y * b) else: tmp = x return tmp
function code(x, y, z, t, a, b) tmp = 0.0 if ((y <= -2e+22) || (!(y <= 3.8e+29) && ((y <= 1.65e+125) || !(y <= 5.8e+189)))) tmp = Float64(1.6453555072203998 * Float64(y * b)); else tmp = x; end return tmp end
function tmp_2 = code(x, y, z, t, a, b) tmp = 0.0; if ((y <= -2e+22) || (~((y <= 3.8e+29)) && ((y <= 1.65e+125) || ~((y <= 5.8e+189))))) tmp = 1.6453555072203998 * (y * b); else tmp = x; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := If[Or[LessEqual[y, -2e+22], And[N[Not[LessEqual[y, 3.8e+29]], $MachinePrecision], Or[LessEqual[y, 1.65e+125], N[Not[LessEqual[y, 5.8e+189]], $MachinePrecision]]]], N[(1.6453555072203998 * N[(y * b), $MachinePrecision]), $MachinePrecision], x]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -2 \cdot 10^{+22} \lor \neg \left(y \leq 3.8 \cdot 10^{+29}\right) \land \left(y \leq 1.65 \cdot 10^{+125} \lor \neg \left(y \leq 5.8 \cdot 10^{+189}\right)\right):\\
\;\;\;\;1.6453555072203998 \cdot \left(y \cdot b\right)\\
\mathbf{else}:\\
\;\;\;\;x\\
\end{array}
\end{array}
if y < -2e22 or 3.79999999999999971e29 < y < 1.65000000000000003e125 or 5.80000000000000038e189 < y Initial program 57.2%
associate-*l/63.8%
*-commutative63.8%
fma-def63.8%
*-commutative63.8%
fma-def63.8%
*-commutative63.8%
fma-def63.8%
*-commutative63.8%
fma-def63.8%
Simplified63.8%
Taylor expanded in z around 0 45.2%
Taylor expanded in x around 0 35.6%
if -2e22 < y < 3.79999999999999971e29 or 1.65000000000000003e125 < y < 5.80000000000000038e189Initial program 48.3%
associate-/l*49.0%
fma-def49.0%
fma-def49.0%
fma-def49.0%
fma-def49.0%
fma-def49.0%
fma-def49.0%
fma-def49.0%
Simplified49.0%
Taylor expanded in z around inf 77.2%
Taylor expanded in x around inf 67.5%
Final simplification53.7%
(FPCore (x y z t a b)
:precision binary64
(if (or (<= z -7.6e-152)
(and (not (<= z -1.05e-223))
(or (<= z -1.6e-273) (not (<= z 4.4e-172)))))
(+ x (/ y 0.31942702700572795))
(* 1.6453555072203998 (* y b))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if ((z <= -7.6e-152) || (!(z <= -1.05e-223) && ((z <= -1.6e-273) || !(z <= 4.4e-172)))) {
tmp = x + (y / 0.31942702700572795);
} else {
tmp = 1.6453555072203998 * (y * b);
}
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 <= (-7.6d-152)) .or. (.not. (z <= (-1.05d-223))) .and. (z <= (-1.6d-273)) .or. (.not. (z <= 4.4d-172))) then
tmp = x + (y / 0.31942702700572795d0)
else
tmp = 1.6453555072203998d0 * (y * b)
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if ((z <= -7.6e-152) || (!(z <= -1.05e-223) && ((z <= -1.6e-273) || !(z <= 4.4e-172)))) {
tmp = x + (y / 0.31942702700572795);
} else {
tmp = 1.6453555072203998 * (y * b);
}
return tmp;
}
def code(x, y, z, t, a, b): tmp = 0 if (z <= -7.6e-152) or (not (z <= -1.05e-223) and ((z <= -1.6e-273) or not (z <= 4.4e-172))): tmp = x + (y / 0.31942702700572795) else: tmp = 1.6453555072203998 * (y * b) return tmp
function code(x, y, z, t, a, b) tmp = 0.0 if ((z <= -7.6e-152) || (!(z <= -1.05e-223) && ((z <= -1.6e-273) || !(z <= 4.4e-172)))) tmp = Float64(x + Float64(y / 0.31942702700572795)); else tmp = Float64(1.6453555072203998 * Float64(y * b)); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) tmp = 0.0; if ((z <= -7.6e-152) || (~((z <= -1.05e-223)) && ((z <= -1.6e-273) || ~((z <= 4.4e-172))))) tmp = x + (y / 0.31942702700572795); else tmp = 1.6453555072203998 * (y * b); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := If[Or[LessEqual[z, -7.6e-152], And[N[Not[LessEqual[z, -1.05e-223]], $MachinePrecision], Or[LessEqual[z, -1.6e-273], N[Not[LessEqual[z, 4.4e-172]], $MachinePrecision]]]], N[(x + N[(y / 0.31942702700572795), $MachinePrecision]), $MachinePrecision], N[(1.6453555072203998 * N[(y * b), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;z \leq -7.6 \cdot 10^{-152} \lor \neg \left(z \leq -1.05 \cdot 10^{-223}\right) \land \left(z \leq -1.6 \cdot 10^{-273} \lor \neg \left(z \leq 4.4 \cdot 10^{-172}\right)\right):\\
\;\;\;\;x + \frac{y}{0.31942702700572795}\\
\mathbf{else}:\\
\;\;\;\;1.6453555072203998 \cdot \left(y \cdot b\right)\\
\end{array}
\end{array}
if z < -7.60000000000000024e-152 or -1.04999999999999991e-223 < z < -1.59999999999999995e-273 or 4.40000000000000018e-172 < z Initial program 40.6%
associate-/l*44.7%
fma-def44.6%
fma-def44.6%
fma-def44.6%
fma-def44.6%
fma-def44.6%
fma-def44.6%
fma-def44.6%
Simplified44.6%
Taylor expanded in z around inf 71.1%
if -7.60000000000000024e-152 < z < -1.04999999999999991e-223 or -1.59999999999999995e-273 < z < 4.40000000000000018e-172Initial program 99.8%
associate-*l/99.7%
*-commutative99.7%
fma-def99.7%
*-commutative99.7%
fma-def99.7%
*-commutative99.7%
fma-def99.7%
*-commutative99.7%
fma-def99.7%
Simplified99.7%
Taylor expanded in z around 0 97.9%
Taylor expanded in x around 0 66.8%
Final simplification70.2%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (* 1.6453555072203998 (* y b))))
(if (<= y -2.4e+22)
t_1
(if (<= y 8.6e+29)
x
(if (<= y 3.8e+124)
(* y (* b 1.6453555072203998))
(if (<= y 8.5e+188) x t_1))))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = 1.6453555072203998 * (y * b);
double tmp;
if (y <= -2.4e+22) {
tmp = t_1;
} else if (y <= 8.6e+29) {
tmp = x;
} else if (y <= 3.8e+124) {
tmp = y * (b * 1.6453555072203998);
} else if (y <= 8.5e+188) {
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 = 1.6453555072203998d0 * (y * b)
if (y <= (-2.4d+22)) then
tmp = t_1
else if (y <= 8.6d+29) then
tmp = x
else if (y <= 3.8d+124) then
tmp = y * (b * 1.6453555072203998d0)
else if (y <= 8.5d+188) 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 = 1.6453555072203998 * (y * b);
double tmp;
if (y <= -2.4e+22) {
tmp = t_1;
} else if (y <= 8.6e+29) {
tmp = x;
} else if (y <= 3.8e+124) {
tmp = y * (b * 1.6453555072203998);
} else if (y <= 8.5e+188) {
tmp = x;
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = 1.6453555072203998 * (y * b) tmp = 0 if y <= -2.4e+22: tmp = t_1 elif y <= 8.6e+29: tmp = x elif y <= 3.8e+124: tmp = y * (b * 1.6453555072203998) elif y <= 8.5e+188: tmp = x else: tmp = t_1 return tmp
function code(x, y, z, t, a, b) t_1 = Float64(1.6453555072203998 * Float64(y * b)) tmp = 0.0 if (y <= -2.4e+22) tmp = t_1; elseif (y <= 8.6e+29) tmp = x; elseif (y <= 3.8e+124) tmp = Float64(y * Float64(b * 1.6453555072203998)); elseif (y <= 8.5e+188) tmp = x; else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a, b) t_1 = 1.6453555072203998 * (y * b); tmp = 0.0; if (y <= -2.4e+22) tmp = t_1; elseif (y <= 8.6e+29) tmp = x; elseif (y <= 3.8e+124) tmp = y * (b * 1.6453555072203998); elseif (y <= 8.5e+188) tmp = x; else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(1.6453555072203998 * N[(y * b), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -2.4e+22], t$95$1, If[LessEqual[y, 8.6e+29], x, If[LessEqual[y, 3.8e+124], N[(y * N[(b * 1.6453555072203998), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 8.5e+188], x, t$95$1]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := 1.6453555072203998 \cdot \left(y \cdot b\right)\\
\mathbf{if}\;y \leq -2.4 \cdot 10^{+22}:\\
\;\;\;\;t_1\\
\mathbf{elif}\;y \leq 8.6 \cdot 10^{+29}:\\
\;\;\;\;x\\
\mathbf{elif}\;y \leq 3.8 \cdot 10^{+124}:\\
\;\;\;\;y \cdot \left(b \cdot 1.6453555072203998\right)\\
\mathbf{elif}\;y \leq 8.5 \cdot 10^{+188}:\\
\;\;\;\;x\\
\mathbf{else}:\\
\;\;\;\;t_1\\
\end{array}
\end{array}
if y < -2.4e22 or 8.49999999999999958e188 < y Initial program 58.4%
associate-*l/64.6%
*-commutative64.6%
fma-def64.6%
*-commutative64.6%
fma-def64.6%
*-commutative64.6%
fma-def64.6%
*-commutative64.6%
fma-def64.6%
Simplified64.6%
Taylor expanded in z around 0 43.2%
Taylor expanded in x around 0 34.4%
if -2.4e22 < y < 8.6000000000000006e29 or 3.7999999999999998e124 < y < 8.49999999999999958e188Initial program 48.3%
associate-/l*49.0%
fma-def49.0%
fma-def49.0%
fma-def49.0%
fma-def49.0%
fma-def49.0%
fma-def49.0%
fma-def49.0%
Simplified49.0%
Taylor expanded in z around inf 77.2%
Taylor expanded in x around inf 67.5%
if 8.6000000000000006e29 < y < 3.7999999999999998e124Initial program 52.7%
associate-*l/60.8%
*-commutative60.8%
fma-def60.8%
*-commutative60.8%
fma-def60.8%
*-commutative60.8%
fma-def60.8%
*-commutative60.8%
fma-def60.8%
Simplified60.8%
Taylor expanded in z around 0 52.8%
Taylor expanded in x around 0 40.0%
*-commutative40.0%
associate-*l*40.0%
Simplified40.0%
Final simplification53.7%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (+ x (/ y 0.31942702700572795))))
(if (<= z -3.5e+125)
t_1
(if (<= z -86.0)
(+ x (* (/ y z) (/ t z)))
(if (<= z 1.5e-9) (+ 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 / 0.31942702700572795);
double tmp;
if (z <= -3.5e+125) {
tmp = t_1;
} else if (z <= -86.0) {
tmp = x + ((y / z) * (t / z));
} else if (z <= 1.5e-9) {
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 / 0.31942702700572795d0)
if (z <= (-3.5d+125)) then
tmp = t_1
else if (z <= (-86.0d0)) then
tmp = x + ((y / z) * (t / z))
else if (z <= 1.5d-9) 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 / 0.31942702700572795);
double tmp;
if (z <= -3.5e+125) {
tmp = t_1;
} else if (z <= -86.0) {
tmp = x + ((y / z) * (t / z));
} else if (z <= 1.5e-9) {
tmp = x + (1.6453555072203998 * (y * b));
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = x + (y / 0.31942702700572795) tmp = 0 if z <= -3.5e+125: tmp = t_1 elif z <= -86.0: tmp = x + ((y / z) * (t / z)) elif z <= 1.5e-9: 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 / 0.31942702700572795)) tmp = 0.0 if (z <= -3.5e+125) tmp = t_1; elseif (z <= -86.0) tmp = Float64(x + Float64(Float64(y / z) * Float64(t / z))); elseif (z <= 1.5e-9) 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 / 0.31942702700572795); tmp = 0.0; if (z <= -3.5e+125) tmp = t_1; elseif (z <= -86.0) tmp = x + ((y / z) * (t / z)); elseif (z <= 1.5e-9) 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 / 0.31942702700572795), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[z, -3.5e+125], t$95$1, If[LessEqual[z, -86.0], N[(x + N[(N[(y / z), $MachinePrecision] * N[(t / z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[z, 1.5e-9], N[(x + N[(1.6453555072203998 * N[(y * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$1]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x + \frac{y}{0.31942702700572795}\\
\mathbf{if}\;z \leq -3.5 \cdot 10^{+125}:\\
\;\;\;\;t_1\\
\mathbf{elif}\;z \leq -86:\\
\;\;\;\;x + \frac{y}{z} \cdot \frac{t}{z}\\
\mathbf{elif}\;z \leq 1.5 \cdot 10^{-9}:\\
\;\;\;\;x + 1.6453555072203998 \cdot \left(y \cdot b\right)\\
\mathbf{else}:\\
\;\;\;\;t_1\\
\end{array}
\end{array}
if z < -3.50000000000000011e125 or 1.49999999999999999e-9 < z Initial program 10.0%
associate-/l*14.0%
fma-def14.0%
fma-def14.0%
fma-def14.0%
fma-def14.0%
fma-def14.0%
fma-def14.0%
fma-def14.0%
Simplified14.0%
Taylor expanded in z around inf 88.9%
if -3.50000000000000011e125 < z < -86Initial program 21.6%
associate-/l*36.7%
fma-def36.7%
fma-def36.7%
fma-def36.7%
fma-def36.7%
fma-def36.7%
fma-def36.7%
fma-def36.7%
Simplified36.7%
Taylor expanded in z around inf 61.1%
associate-*r/61.1%
metadata-eval61.1%
mul-1-neg61.1%
*-commutative61.1%
unpow261.1%
Simplified61.1%
Taylor expanded in z around inf 88.7%
Taylor expanded in t around inf 80.1%
unpow280.1%
times-frac87.5%
Simplified87.5%
if -86 < z < 1.49999999999999999e-9Initial program 99.7%
associate-*l/99.7%
*-commutative99.7%
fma-def99.7%
*-commutative99.7%
fma-def99.7%
*-commutative99.7%
fma-def99.6%
*-commutative99.6%
fma-def99.7%
Simplified99.7%
Taylor expanded in z around 0 81.0%
Final simplification85.1%
(FPCore (x y z t a b) :precision binary64 (if (or (<= z -580000.0) (not (<= z 0.17))) (+ x (/ y 0.31942702700572795)) (+ x (* 1.6453555072203998 (* y b)))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if ((z <= -580000.0) || !(z <= 0.17)) {
tmp = x + (y / 0.31942702700572795);
} else {
tmp = x + (1.6453555072203998 * (y * b));
}
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 <= (-580000.0d0)) .or. (.not. (z <= 0.17d0))) then
tmp = x + (y / 0.31942702700572795d0)
else
tmp = x + (1.6453555072203998d0 * (y * b))
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if ((z <= -580000.0) || !(z <= 0.17)) {
tmp = x + (y / 0.31942702700572795);
} else {
tmp = x + (1.6453555072203998 * (y * b));
}
return tmp;
}
def code(x, y, z, t, a, b): tmp = 0 if (z <= -580000.0) or not (z <= 0.17): tmp = x + (y / 0.31942702700572795) else: tmp = x + (1.6453555072203998 * (y * b)) return tmp
function code(x, y, z, t, a, b) tmp = 0.0 if ((z <= -580000.0) || !(z <= 0.17)) tmp = Float64(x + Float64(y / 0.31942702700572795)); else tmp = Float64(x + Float64(1.6453555072203998 * Float64(y * b))); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) tmp = 0.0; if ((z <= -580000.0) || ~((z <= 0.17))) tmp = x + (y / 0.31942702700572795); else tmp = x + (1.6453555072203998 * (y * b)); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := If[Or[LessEqual[z, -580000.0], N[Not[LessEqual[z, 0.17]], $MachinePrecision]], N[(x + N[(y / 0.31942702700572795), $MachinePrecision]), $MachinePrecision], N[(x + N[(1.6453555072203998 * N[(y * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;z \leq -580000 \lor \neg \left(z \leq 0.17\right):\\
\;\;\;\;x + \frac{y}{0.31942702700572795}\\
\mathbf{else}:\\
\;\;\;\;x + 1.6453555072203998 \cdot \left(y \cdot b\right)\\
\end{array}
\end{array}
if z < -5.8e5 or 0.170000000000000012 < z Initial program 10.8%
associate-/l*16.9%
fma-def16.9%
fma-def16.9%
fma-def16.9%
fma-def16.9%
fma-def16.9%
fma-def16.9%
fma-def16.9%
Simplified16.9%
Taylor expanded in z around inf 84.5%
if -5.8e5 < z < 0.170000000000000012Initial program 99.7%
associate-*l/99.7%
*-commutative99.7%
fma-def99.7%
*-commutative99.7%
fma-def99.7%
*-commutative99.7%
fma-def99.6%
*-commutative99.6%
fma-def99.6%
Simplified99.6%
Taylor expanded in z around 0 80.5%
Final simplification82.6%
(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 52.1%
associate-/l*55.4%
fma-def55.4%
fma-def55.4%
fma-def55.4%
fma-def55.4%
fma-def55.4%
fma-def55.4%
fma-def55.4%
Simplified55.4%
Taylor expanded in z around inf 62.2%
Taylor expanded in x around inf 44.6%
Final simplification44.6%
(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 2023257
(FPCore (x y z t a b)
:name "Numeric.SpecFunctions:logGamma from math-functions-0.1.5.2, D"
:precision binary64
:herbie-target
(if (< z -6.499344996252632e+53) (+ x (* (+ (- 3.13060547623 (/ 36.527041698806414 z)) (/ t (* z z))) (/ y 1.0))) (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)))) (+ x (* (+ (- 3.13060547623 (/ 36.527041698806414 z)) (/ t (* z z))) (/ y 1.0)))))
(+ 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))))