
(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 19 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
(if (<=
(/
(*
(+
b
(* (+ a (* (+ t (* (+ 11.1667541262 (* 3.13060547623 z)) z)) z)) z))
y)
(+
0.607771387771
(*
(+ 11.9400905721 (* (+ 31.4690115749 (* (+ 15.234687407 z) z)) z))
z)))
INFINITY)
(+
(/
y
(/
(fma
(fma (fma (+ 15.234687407 z) z 31.4690115749) z 11.9400905721)
z
0.607771387771)
(fma (fma (fma (fma 3.13060547623 z 11.1667541262) z t) z a) z b)))
x)
(+
(-
(fma (/ t z) (/ y z) (fma (/ y z) 11.1667541262 (* 3.13060547623 y)))
(fma
(/ (* -36.52704169880642 y) z)
(/ 15.234687407 z)
(fma (/ y (* z z)) 98.5170599679272 (* 47.69379582500642 (/ y z)))))
x)))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if ((((b + ((a + ((t + ((11.1667541262 + (3.13060547623 * z)) * z)) * z)) * z)) * y) / (0.607771387771 + ((11.9400905721 + ((31.4690115749 + ((15.234687407 + z) * z)) * z)) * z))) <= ((double) INFINITY)) {
tmp = (y / (fma(fma(fma((15.234687407 + z), z, 31.4690115749), z, 11.9400905721), z, 0.607771387771) / fma(fma(fma(fma(3.13060547623, z, 11.1667541262), z, t), z, a), z, b))) + x;
} else {
tmp = (fma((t / z), (y / z), fma((y / z), 11.1667541262, (3.13060547623 * y))) - fma(((-36.52704169880642 * y) / z), (15.234687407 / z), fma((y / (z * z)), 98.5170599679272, (47.69379582500642 * (y / z))))) + x;
}
return tmp;
}
function code(x, y, z, t, a, b) tmp = 0.0 if (Float64(Float64(Float64(b + Float64(Float64(a + Float64(Float64(t + Float64(Float64(11.1667541262 + Float64(3.13060547623 * z)) * z)) * z)) * z)) * y) / Float64(0.607771387771 + Float64(Float64(11.9400905721 + Float64(Float64(31.4690115749 + Float64(Float64(15.234687407 + z) * z)) * z)) * z))) <= Inf) tmp = Float64(Float64(y / Float64(fma(fma(fma(Float64(15.234687407 + z), z, 31.4690115749), z, 11.9400905721), z, 0.607771387771) / fma(fma(fma(fma(3.13060547623, z, 11.1667541262), z, t), z, a), z, b))) + x); else tmp = Float64(Float64(fma(Float64(t / z), Float64(y / z), fma(Float64(y / z), 11.1667541262, Float64(3.13060547623 * y))) - fma(Float64(Float64(-36.52704169880642 * y) / z), Float64(15.234687407 / z), fma(Float64(y / Float64(z * z)), 98.5170599679272, Float64(47.69379582500642 * Float64(y / z))))) + x); end return tmp end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[N[(N[(N[(b + N[(N[(a + N[(N[(t + N[(N[(11.1667541262 + N[(3.13060547623 * z), $MachinePrecision]), $MachinePrecision] * z), $MachinePrecision]), $MachinePrecision] * z), $MachinePrecision]), $MachinePrecision] * z), $MachinePrecision]), $MachinePrecision] * y), $MachinePrecision] / N[(0.607771387771 + N[(N[(11.9400905721 + N[(N[(31.4690115749 + N[(N[(15.234687407 + z), $MachinePrecision] * z), $MachinePrecision]), $MachinePrecision] * z), $MachinePrecision]), $MachinePrecision] * z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], Infinity], N[(N[(y / N[(N[(N[(N[(N[(15.234687407 + z), $MachinePrecision] * z + 31.4690115749), $MachinePrecision] * z + 11.9400905721), $MachinePrecision] * z + 0.607771387771), $MachinePrecision] / N[(N[(N[(N[(3.13060547623 * z + 11.1667541262), $MachinePrecision] * z + t), $MachinePrecision] * z + a), $MachinePrecision] * z + b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + x), $MachinePrecision], N[(N[(N[(N[(t / z), $MachinePrecision] * N[(y / z), $MachinePrecision] + N[(N[(y / z), $MachinePrecision] * 11.1667541262 + N[(3.13060547623 * y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - N[(N[(N[(-36.52704169880642 * y), $MachinePrecision] / z), $MachinePrecision] * N[(15.234687407 / z), $MachinePrecision] + N[(N[(y / N[(z * z), $MachinePrecision]), $MachinePrecision] * 98.5170599679272 + N[(47.69379582500642 * N[(y / z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + x), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\frac{\left(b + \left(a + \left(t + \left(11.1667541262 + 3.13060547623 \cdot z\right) \cdot z\right) \cdot z\right) \cdot z\right) \cdot y}{0.607771387771 + \left(11.9400905721 + \left(31.4690115749 + \left(15.234687407 + z\right) \cdot z\right) \cdot z\right) \cdot z} \leq \infty:\\
\;\;\;\;\frac{y}{\frac{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(15.234687407 + z, z, 31.4690115749\right), z, 11.9400905721\right), z, 0.607771387771\right)}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(3.13060547623, z, 11.1667541262\right), z, t\right), z, a\right), z, b\right)}} + x\\
\mathbf{else}:\\
\;\;\;\;\left(\mathsf{fma}\left(\frac{t}{z}, \frac{y}{z}, \mathsf{fma}\left(\frac{y}{z}, 11.1667541262, 3.13060547623 \cdot y\right)\right) - \mathsf{fma}\left(\frac{-36.52704169880642 \cdot y}{z}, \frac{15.234687407}{z}, \mathsf{fma}\left(\frac{y}{z \cdot z}, 98.5170599679272, 47.69379582500642 \cdot \frac{y}{z}\right)\right)\right) + x\\
\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 88.3%
lift-/.f64N/A
lift-*.f64N/A
associate-/l*N/A
clear-numN/A
un-div-invN/A
lower-/.f64N/A
lower-/.f6495.4
Applied rewrites95.4%
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 t around inf
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
unpow2N/A
lower-*.f6410.5
Applied rewrites10.5%
Taylor expanded in z around 0
Applied rewrites8.6%
Taylor expanded in z around inf
Applied rewrites99.8%
Final simplification97.1%
(FPCore (x y z t a b)
:precision binary64
(if (<=
(/
(*
(+
b
(* (+ a (* (+ t (* (+ 11.1667541262 (* 3.13060547623 z)) z)) z)) z))
y)
(+
0.607771387771
(*
(+ 11.9400905721 (* (+ 31.4690115749 (* (+ 15.234687407 z) z)) z))
z)))
INFINITY)
(+
(/
y
(/
(fma
(fma (fma (+ 15.234687407 z) z 31.4690115749) z 11.9400905721)
z
0.607771387771)
(fma (fma (fma (fma 3.13060547623 z 11.1667541262) z t) z a) z b)))
x)
(fma -36.52704169880642 (/ y z) (fma 3.13060547623 y x))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if ((((b + ((a + ((t + ((11.1667541262 + (3.13060547623 * z)) * z)) * z)) * z)) * y) / (0.607771387771 + ((11.9400905721 + ((31.4690115749 + ((15.234687407 + z) * z)) * z)) * z))) <= ((double) INFINITY)) {
tmp = (y / (fma(fma(fma((15.234687407 + z), z, 31.4690115749), z, 11.9400905721), z, 0.607771387771) / fma(fma(fma(fma(3.13060547623, z, 11.1667541262), z, t), z, a), z, b))) + x;
} else {
tmp = fma(-36.52704169880642, (y / z), fma(3.13060547623, y, x));
}
return tmp;
}
function code(x, y, z, t, a, b) tmp = 0.0 if (Float64(Float64(Float64(b + Float64(Float64(a + Float64(Float64(t + Float64(Float64(11.1667541262 + Float64(3.13060547623 * z)) * z)) * z)) * z)) * y) / Float64(0.607771387771 + Float64(Float64(11.9400905721 + Float64(Float64(31.4690115749 + Float64(Float64(15.234687407 + z) * z)) * z)) * z))) <= Inf) tmp = Float64(Float64(y / Float64(fma(fma(fma(Float64(15.234687407 + z), z, 31.4690115749), z, 11.9400905721), z, 0.607771387771) / fma(fma(fma(fma(3.13060547623, z, 11.1667541262), z, t), z, a), z, b))) + x); else tmp = fma(-36.52704169880642, Float64(y / z), fma(3.13060547623, y, x)); end return tmp end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[N[(N[(N[(b + N[(N[(a + N[(N[(t + N[(N[(11.1667541262 + N[(3.13060547623 * z), $MachinePrecision]), $MachinePrecision] * z), $MachinePrecision]), $MachinePrecision] * z), $MachinePrecision]), $MachinePrecision] * z), $MachinePrecision]), $MachinePrecision] * y), $MachinePrecision] / N[(0.607771387771 + N[(N[(11.9400905721 + N[(N[(31.4690115749 + N[(N[(15.234687407 + z), $MachinePrecision] * z), $MachinePrecision]), $MachinePrecision] * z), $MachinePrecision]), $MachinePrecision] * z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], Infinity], N[(N[(y / N[(N[(N[(N[(N[(15.234687407 + z), $MachinePrecision] * z + 31.4690115749), $MachinePrecision] * z + 11.9400905721), $MachinePrecision] * z + 0.607771387771), $MachinePrecision] / N[(N[(N[(N[(3.13060547623 * z + 11.1667541262), $MachinePrecision] * z + t), $MachinePrecision] * z + a), $MachinePrecision] * z + b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + x), $MachinePrecision], N[(-36.52704169880642 * N[(y / z), $MachinePrecision] + N[(3.13060547623 * y + x), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\frac{\left(b + \left(a + \left(t + \left(11.1667541262 + 3.13060547623 \cdot z\right) \cdot z\right) \cdot z\right) \cdot z\right) \cdot y}{0.607771387771 + \left(11.9400905721 + \left(31.4690115749 + \left(15.234687407 + z\right) \cdot z\right) \cdot z\right) \cdot z} \leq \infty:\\
\;\;\;\;\frac{y}{\frac{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(15.234687407 + z, z, 31.4690115749\right), z, 11.9400905721\right), z, 0.607771387771\right)}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(3.13060547623, z, 11.1667541262\right), z, t\right), z, a\right), z, b\right)}} + x\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(-36.52704169880642, \frac{y}{z}, \mathsf{fma}\left(3.13060547623, y, x\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 88.3%
lift-/.f64N/A
lift-*.f64N/A
associate-/l*N/A
clear-numN/A
un-div-invN/A
lower-/.f64N/A
lower-/.f6495.4
Applied rewrites95.4%
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
associate--l+N/A
+-commutativeN/A
associate--l+N/A
+-commutativeN/A
distribute-rgt-out--N/A
metadata-evalN/A
metadata-evalN/A
metadata-evalN/A
times-fracN/A
distribute-rgt-out--N/A
*-commutativeN/A
mul-1-negN/A
distribute-neg-frac2N/A
mul-1-negN/A
associate-+l+N/A
Applied rewrites95.1%
Final simplification95.3%
(FPCore (x y z t a b)
:precision binary64
(if (<=
(/
(*
(+
b
(* (+ a (* (+ t (* (+ 11.1667541262 (* 3.13060547623 z)) z)) z)) z))
y)
(+
0.607771387771
(*
(+ 11.9400905721 (* (+ 31.4690115749 (* (+ 15.234687407 z) z)) z))
z)))
INFINITY)
(fma
(/
(fma (fma (fma (fma 3.13060547623 z 11.1667541262) z t) z a) z b)
(fma
(fma (fma (+ 15.234687407 z) z 31.4690115749) z 11.9400905721)
z
0.607771387771))
y
x)
(fma -36.52704169880642 (/ y z) (fma 3.13060547623 y x))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if ((((b + ((a + ((t + ((11.1667541262 + (3.13060547623 * z)) * z)) * z)) * z)) * y) / (0.607771387771 + ((11.9400905721 + ((31.4690115749 + ((15.234687407 + z) * z)) * z)) * z))) <= ((double) INFINITY)) {
tmp = fma((fma(fma(fma(fma(3.13060547623, z, 11.1667541262), z, t), z, a), z, b) / fma(fma(fma((15.234687407 + z), z, 31.4690115749), z, 11.9400905721), z, 0.607771387771)), y, x);
} else {
tmp = fma(-36.52704169880642, (y / z), fma(3.13060547623, y, x));
}
return tmp;
}
function code(x, y, z, t, a, b) tmp = 0.0 if (Float64(Float64(Float64(b + Float64(Float64(a + Float64(Float64(t + Float64(Float64(11.1667541262 + Float64(3.13060547623 * z)) * z)) * z)) * z)) * y) / Float64(0.607771387771 + Float64(Float64(11.9400905721 + Float64(Float64(31.4690115749 + Float64(Float64(15.234687407 + z) * z)) * z)) * z))) <= Inf) tmp = fma(Float64(fma(fma(fma(fma(3.13060547623, z, 11.1667541262), z, t), z, a), z, b) / fma(fma(fma(Float64(15.234687407 + z), z, 31.4690115749), z, 11.9400905721), z, 0.607771387771)), y, x); else tmp = fma(-36.52704169880642, Float64(y / z), fma(3.13060547623, y, x)); end return tmp end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[N[(N[(N[(b + N[(N[(a + N[(N[(t + N[(N[(11.1667541262 + N[(3.13060547623 * z), $MachinePrecision]), $MachinePrecision] * z), $MachinePrecision]), $MachinePrecision] * z), $MachinePrecision]), $MachinePrecision] * z), $MachinePrecision]), $MachinePrecision] * y), $MachinePrecision] / N[(0.607771387771 + N[(N[(11.9400905721 + N[(N[(31.4690115749 + N[(N[(15.234687407 + z), $MachinePrecision] * z), $MachinePrecision]), $MachinePrecision] * z), $MachinePrecision]), $MachinePrecision] * z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], Infinity], N[(N[(N[(N[(N[(N[(3.13060547623 * z + 11.1667541262), $MachinePrecision] * z + t), $MachinePrecision] * z + a), $MachinePrecision] * z + b), $MachinePrecision] / N[(N[(N[(N[(15.234687407 + z), $MachinePrecision] * z + 31.4690115749), $MachinePrecision] * z + 11.9400905721), $MachinePrecision] * z + 0.607771387771), $MachinePrecision]), $MachinePrecision] * y + x), $MachinePrecision], N[(-36.52704169880642 * N[(y / z), $MachinePrecision] + N[(3.13060547623 * y + x), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\frac{\left(b + \left(a + \left(t + \left(11.1667541262 + 3.13060547623 \cdot z\right) \cdot z\right) \cdot z\right) \cdot z\right) \cdot y}{0.607771387771 + \left(11.9400905721 + \left(31.4690115749 + \left(15.234687407 + z\right) \cdot z\right) \cdot z\right) \cdot z} \leq \infty:\\
\;\;\;\;\mathsf{fma}\left(\frac{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(3.13060547623, z, 11.1667541262\right), z, t\right), z, a\right), z, b\right)}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(15.234687407 + z, z, 31.4690115749\right), z, 11.9400905721\right), z, 0.607771387771\right)}, y, x\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(-36.52704169880642, \frac{y}{z}, \mathsf{fma}\left(3.13060547623, y, x\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 88.3%
lift-+.f64N/A
+-commutativeN/A
lift-/.f64N/A
lift-*.f64N/A
associate-/l*N/A
*-commutativeN/A
lower-fma.f64N/A
Applied rewrites95.4%
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
associate--l+N/A
+-commutativeN/A
associate--l+N/A
+-commutativeN/A
distribute-rgt-out--N/A
metadata-evalN/A
metadata-evalN/A
metadata-evalN/A
times-fracN/A
distribute-rgt-out--N/A
*-commutativeN/A
mul-1-negN/A
distribute-neg-frac2N/A
mul-1-negN/A
associate-+l+N/A
Applied rewrites95.1%
Final simplification95.3%
(FPCore (x y z t a b)
:precision binary64
(if (<=
(/
(*
(+
b
(* (+ a (* (+ t (* (+ 11.1667541262 (* 3.13060547623 z)) z)) z)) z))
y)
(+
0.607771387771
(*
(+ 11.9400905721 (* (+ 31.4690115749 (* (+ 15.234687407 z) z)) z))
z)))
-5e+67)
(* 1.6453555072203998 (* b y))
(fma 3.13060547623 y x)))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if ((((b + ((a + ((t + ((11.1667541262 + (3.13060547623 * z)) * z)) * z)) * z)) * y) / (0.607771387771 + ((11.9400905721 + ((31.4690115749 + ((15.234687407 + z) * z)) * z)) * z))) <= -5e+67) {
tmp = 1.6453555072203998 * (b * y);
} else {
tmp = fma(3.13060547623, y, x);
}
return tmp;
}
function code(x, y, z, t, a, b) tmp = 0.0 if (Float64(Float64(Float64(b + Float64(Float64(a + Float64(Float64(t + Float64(Float64(11.1667541262 + Float64(3.13060547623 * z)) * z)) * z)) * z)) * y) / Float64(0.607771387771 + Float64(Float64(11.9400905721 + Float64(Float64(31.4690115749 + Float64(Float64(15.234687407 + z) * z)) * z)) * z))) <= -5e+67) tmp = Float64(1.6453555072203998 * Float64(b * y)); else tmp = fma(3.13060547623, y, x); end return tmp end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[N[(N[(N[(b + N[(N[(a + N[(N[(t + N[(N[(11.1667541262 + N[(3.13060547623 * z), $MachinePrecision]), $MachinePrecision] * z), $MachinePrecision]), $MachinePrecision] * z), $MachinePrecision]), $MachinePrecision] * z), $MachinePrecision]), $MachinePrecision] * y), $MachinePrecision] / N[(0.607771387771 + N[(N[(11.9400905721 + N[(N[(31.4690115749 + N[(N[(15.234687407 + z), $MachinePrecision] * z), $MachinePrecision]), $MachinePrecision] * z), $MachinePrecision]), $MachinePrecision] * z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], -5e+67], N[(1.6453555072203998 * N[(b * y), $MachinePrecision]), $MachinePrecision], N[(3.13060547623 * y + x), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\frac{\left(b + \left(a + \left(t + \left(11.1667541262 + 3.13060547623 \cdot z\right) \cdot z\right) \cdot z\right) \cdot z\right) \cdot y}{0.607771387771 + \left(11.9400905721 + \left(31.4690115749 + \left(15.234687407 + z\right) \cdot z\right) \cdot z\right) \cdot z} \leq -5 \cdot 10^{+67}:\\
\;\;\;\;1.6453555072203998 \cdot \left(b \cdot y\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(3.13060547623, y, x\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))) < -4.99999999999999976e67Initial program 75.1%
Taylor expanded in z around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower-*.f6442.6
Applied rewrites42.6%
Taylor expanded in x around 0
Applied rewrites34.8%
if -4.99999999999999976e67 < (/.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 50.5%
Taylor expanded in z around inf
+-commutativeN/A
lower-fma.f6471.7
Applied rewrites71.7%
Final simplification65.6%
(FPCore (x y z t a b)
:precision binary64
(if (<= z -6.2e+69)
(fma 3.13060547623 y x)
(if (<= z 9e+42)
(+
(/
(* (fma (fma t z a) z b) y)
(+
0.607771387771
(*
(+ 11.9400905721 (* (+ 31.4690115749 (* (+ 15.234687407 z) z)) z))
z)))
x)
(fma -36.52704169880642 (/ y z) (fma 3.13060547623 y x)))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (z <= -6.2e+69) {
tmp = fma(3.13060547623, y, x);
} else if (z <= 9e+42) {
tmp = ((fma(fma(t, z, a), z, b) * y) / (0.607771387771 + ((11.9400905721 + ((31.4690115749 + ((15.234687407 + z) * z)) * z)) * z))) + x;
} else {
tmp = fma(-36.52704169880642, (y / z), fma(3.13060547623, y, x));
}
return tmp;
}
function code(x, y, z, t, a, b) tmp = 0.0 if (z <= -6.2e+69) tmp = fma(3.13060547623, y, x); elseif (z <= 9e+42) tmp = Float64(Float64(Float64(fma(fma(t, z, a), z, b) * y) / Float64(0.607771387771 + Float64(Float64(11.9400905721 + Float64(Float64(31.4690115749 + Float64(Float64(15.234687407 + z) * z)) * z)) * z))) + x); else tmp = fma(-36.52704169880642, Float64(y / z), fma(3.13060547623, y, x)); end return tmp end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[z, -6.2e+69], N[(3.13060547623 * y + x), $MachinePrecision], If[LessEqual[z, 9e+42], N[(N[(N[(N[(N[(t * z + a), $MachinePrecision] * z + b), $MachinePrecision] * y), $MachinePrecision] / N[(0.607771387771 + N[(N[(11.9400905721 + N[(N[(31.4690115749 + N[(N[(15.234687407 + z), $MachinePrecision] * z), $MachinePrecision]), $MachinePrecision] * z), $MachinePrecision]), $MachinePrecision] * z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + x), $MachinePrecision], N[(-36.52704169880642 * N[(y / z), $MachinePrecision] + N[(3.13060547623 * y + x), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;z \leq -6.2 \cdot 10^{+69}:\\
\;\;\;\;\mathsf{fma}\left(3.13060547623, y, x\right)\\
\mathbf{elif}\;z \leq 9 \cdot 10^{+42}:\\
\;\;\;\;\frac{\mathsf{fma}\left(\mathsf{fma}\left(t, z, a\right), z, b\right) \cdot y}{0.607771387771 + \left(11.9400905721 + \left(31.4690115749 + \left(15.234687407 + z\right) \cdot z\right) \cdot z\right) \cdot z} + x\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(-36.52704169880642, \frac{y}{z}, \mathsf{fma}\left(3.13060547623, y, x\right)\right)\\
\end{array}
\end{array}
if z < -6.1999999999999997e69Initial program 0.3%
Taylor expanded in z around inf
+-commutativeN/A
lower-fma.f6493.7
Applied rewrites93.7%
if -6.1999999999999997e69 < z < 9.00000000000000025e42Initial program 92.1%
Taylor expanded in z around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f6491.4
Applied rewrites91.4%
if 9.00000000000000025e42 < z Initial program 3.9%
Taylor expanded in z around inf
associate--l+N/A
+-commutativeN/A
associate--l+N/A
+-commutativeN/A
distribute-rgt-out--N/A
metadata-evalN/A
metadata-evalN/A
metadata-evalN/A
times-fracN/A
distribute-rgt-out--N/A
*-commutativeN/A
mul-1-negN/A
distribute-neg-frac2N/A
mul-1-negN/A
associate-+l+N/A
Applied rewrites94.9%
Final simplification92.6%
(FPCore (x y z t a b)
:precision binary64
(if (<= z -9.5e+31)
(fma 3.13060547623 y x)
(if (<= z 1.1e+20)
(+
(/
y
(/
0.607771387771
(fma (fma (fma (fma 3.13060547623 z 11.1667541262) z t) z a) z b)))
x)
(fma -36.52704169880642 (/ y z) (fma 3.13060547623 y x)))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (z <= -9.5e+31) {
tmp = fma(3.13060547623, y, x);
} else if (z <= 1.1e+20) {
tmp = (y / (0.607771387771 / fma(fma(fma(fma(3.13060547623, z, 11.1667541262), z, t), z, a), z, b))) + x;
} else {
tmp = fma(-36.52704169880642, (y / z), fma(3.13060547623, y, x));
}
return tmp;
}
function code(x, y, z, t, a, b) tmp = 0.0 if (z <= -9.5e+31) tmp = fma(3.13060547623, y, x); elseif (z <= 1.1e+20) tmp = Float64(Float64(y / Float64(0.607771387771 / fma(fma(fma(fma(3.13060547623, z, 11.1667541262), z, t), z, a), z, b))) + x); else tmp = fma(-36.52704169880642, Float64(y / z), fma(3.13060547623, y, x)); end return tmp end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[z, -9.5e+31], N[(3.13060547623 * y + x), $MachinePrecision], If[LessEqual[z, 1.1e+20], N[(N[(y / N[(0.607771387771 / N[(N[(N[(N[(3.13060547623 * z + 11.1667541262), $MachinePrecision] * z + t), $MachinePrecision] * z + a), $MachinePrecision] * z + b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + x), $MachinePrecision], N[(-36.52704169880642 * N[(y / z), $MachinePrecision] + N[(3.13060547623 * y + x), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;z \leq -9.5 \cdot 10^{+31}:\\
\;\;\;\;\mathsf{fma}\left(3.13060547623, y, x\right)\\
\mathbf{elif}\;z \leq 1.1 \cdot 10^{+20}:\\
\;\;\;\;\frac{y}{\frac{0.607771387771}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(3.13060547623, z, 11.1667541262\right), z, t\right), z, a\right), z, b\right)}} + x\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(-36.52704169880642, \frac{y}{z}, \mathsf{fma}\left(3.13060547623, y, x\right)\right)\\
\end{array}
\end{array}
if z < -9.5000000000000008e31Initial program 7.2%
Taylor expanded in z around inf
+-commutativeN/A
lower-fma.f6481.3
Applied rewrites81.3%
if -9.5000000000000008e31 < z < 1.1e20Initial program 99.0%
lift-/.f64N/A
lift-*.f64N/A
associate-/l*N/A
clear-numN/A
un-div-invN/A
lower-/.f64N/A
lower-/.f6499.7
Applied rewrites99.8%
Taylor expanded in z around 0
Applied rewrites95.8%
if 1.1e20 < z Initial program 8.6%
Taylor expanded in z around inf
associate--l+N/A
+-commutativeN/A
associate--l+N/A
+-commutativeN/A
distribute-rgt-out--N/A
metadata-evalN/A
metadata-evalN/A
metadata-evalN/A
times-fracN/A
distribute-rgt-out--N/A
*-commutativeN/A
mul-1-negN/A
distribute-neg-frac2N/A
mul-1-negN/A
associate-+l+N/A
Applied rewrites88.1%
Final simplification90.3%
(FPCore (x y z t a b)
:precision binary64
(if (<= z -4.05e+43)
(fma 3.13060547623 y x)
(if (<= z 1e-92)
(+ (/ (fma (fma (* z y) t (* a y)) z (* b y)) 0.607771387771) x)
(if (<= z 1.1e+20)
(+
(/
(* (fma (* z z) (fma (fma 3.13060547623 z 11.1667541262) z t) b) y)
0.607771387771)
x)
(fma -36.52704169880642 (/ y z) (fma 3.13060547623 y x))))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (z <= -4.05e+43) {
tmp = fma(3.13060547623, y, x);
} else if (z <= 1e-92) {
tmp = (fma(fma((z * y), t, (a * y)), z, (b * y)) / 0.607771387771) + x;
} else if (z <= 1.1e+20) {
tmp = ((fma((z * z), fma(fma(3.13060547623, z, 11.1667541262), z, t), b) * y) / 0.607771387771) + x;
} else {
tmp = fma(-36.52704169880642, (y / z), fma(3.13060547623, y, x));
}
return tmp;
}
function code(x, y, z, t, a, b) tmp = 0.0 if (z <= -4.05e+43) tmp = fma(3.13060547623, y, x); elseif (z <= 1e-92) tmp = Float64(Float64(fma(fma(Float64(z * y), t, Float64(a * y)), z, Float64(b * y)) / 0.607771387771) + x); elseif (z <= 1.1e+20) tmp = Float64(Float64(Float64(fma(Float64(z * z), fma(fma(3.13060547623, z, 11.1667541262), z, t), b) * y) / 0.607771387771) + x); else tmp = fma(-36.52704169880642, Float64(y / z), fma(3.13060547623, y, x)); end return tmp end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[z, -4.05e+43], N[(3.13060547623 * y + x), $MachinePrecision], If[LessEqual[z, 1e-92], N[(N[(N[(N[(N[(z * y), $MachinePrecision] * t + N[(a * y), $MachinePrecision]), $MachinePrecision] * z + N[(b * y), $MachinePrecision]), $MachinePrecision] / 0.607771387771), $MachinePrecision] + x), $MachinePrecision], If[LessEqual[z, 1.1e+20], N[(N[(N[(N[(N[(z * z), $MachinePrecision] * N[(N[(3.13060547623 * z + 11.1667541262), $MachinePrecision] * z + t), $MachinePrecision] + b), $MachinePrecision] * y), $MachinePrecision] / 0.607771387771), $MachinePrecision] + x), $MachinePrecision], N[(-36.52704169880642 * N[(y / z), $MachinePrecision] + N[(3.13060547623 * y + x), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;z \leq -4.05 \cdot 10^{+43}:\\
\;\;\;\;\mathsf{fma}\left(3.13060547623, y, x\right)\\
\mathbf{elif}\;z \leq 10^{-92}:\\
\;\;\;\;\frac{\mathsf{fma}\left(\mathsf{fma}\left(z \cdot y, t, a \cdot y\right), z, b \cdot y\right)}{0.607771387771} + x\\
\mathbf{elif}\;z \leq 1.1 \cdot 10^{+20}:\\
\;\;\;\;\frac{\mathsf{fma}\left(z \cdot z, \mathsf{fma}\left(\mathsf{fma}\left(3.13060547623, z, 11.1667541262\right), z, t\right), b\right) \cdot y}{0.607771387771} + x\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(-36.52704169880642, \frac{y}{z}, \mathsf{fma}\left(3.13060547623, y, x\right)\right)\\
\end{array}
\end{array}
if z < -4.0499999999999998e43Initial program 2.6%
Taylor expanded in z around inf
+-commutativeN/A
lower-fma.f6484.6
Applied rewrites84.6%
if -4.0499999999999998e43 < z < 9.99999999999999988e-93Initial program 97.3%
Taylor expanded in z around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-*.f64N/A
lower-*.f6488.6
Applied rewrites88.6%
Taylor expanded in z around 0
Applied rewrites85.6%
if 9.99999999999999988e-93 < z < 1.1e20Initial program 99.6%
Taylor expanded in t around inf
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
unpow2N/A
lower-*.f6471.4
Applied rewrites71.4%
Taylor expanded in z around 0
Applied rewrites66.7%
Taylor expanded in a around 0
*-commutativeN/A
lower-*.f64N/A
+-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f6483.3
Applied rewrites83.3%
if 1.1e20 < z Initial program 8.6%
Taylor expanded in z around inf
associate--l+N/A
+-commutativeN/A
associate--l+N/A
+-commutativeN/A
distribute-rgt-out--N/A
metadata-evalN/A
metadata-evalN/A
metadata-evalN/A
times-fracN/A
distribute-rgt-out--N/A
*-commutativeN/A
mul-1-negN/A
distribute-neg-frac2N/A
mul-1-negN/A
associate-+l+N/A
Applied rewrites88.1%
Final simplification85.8%
(FPCore (x y z t a b)
:precision binary64
(if (<= z -1.6e+25)
(fma 3.13060547623 y x)
(if (<= z 2.25e+19)
(+
(/
(* (fma (fma (* z z) (fma 3.13060547623 z 11.1667541262) a) z b) y)
0.607771387771)
x)
(fma -36.52704169880642 (/ y z) (fma 3.13060547623 y x)))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (z <= -1.6e+25) {
tmp = fma(3.13060547623, y, x);
} else if (z <= 2.25e+19) {
tmp = ((fma(fma((z * z), fma(3.13060547623, z, 11.1667541262), a), z, b) * y) / 0.607771387771) + x;
} else {
tmp = fma(-36.52704169880642, (y / z), fma(3.13060547623, y, x));
}
return tmp;
}
function code(x, y, z, t, a, b) tmp = 0.0 if (z <= -1.6e+25) tmp = fma(3.13060547623, y, x); elseif (z <= 2.25e+19) tmp = Float64(Float64(Float64(fma(fma(Float64(z * z), fma(3.13060547623, z, 11.1667541262), a), z, b) * y) / 0.607771387771) + x); else tmp = fma(-36.52704169880642, Float64(y / z), fma(3.13060547623, y, x)); end return tmp end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[z, -1.6e+25], N[(3.13060547623 * y + x), $MachinePrecision], If[LessEqual[z, 2.25e+19], N[(N[(N[(N[(N[(N[(z * z), $MachinePrecision] * N[(3.13060547623 * z + 11.1667541262), $MachinePrecision] + a), $MachinePrecision] * z + b), $MachinePrecision] * y), $MachinePrecision] / 0.607771387771), $MachinePrecision] + x), $MachinePrecision], N[(-36.52704169880642 * N[(y / z), $MachinePrecision] + N[(3.13060547623 * y + x), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;z \leq -1.6 \cdot 10^{+25}:\\
\;\;\;\;\mathsf{fma}\left(3.13060547623, y, x\right)\\
\mathbf{elif}\;z \leq 2.25 \cdot 10^{+19}:\\
\;\;\;\;\frac{\mathsf{fma}\left(\mathsf{fma}\left(z \cdot z, \mathsf{fma}\left(3.13060547623, z, 11.1667541262\right), a\right), z, b\right) \cdot y}{0.607771387771} + x\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(-36.52704169880642, \frac{y}{z}, \mathsf{fma}\left(3.13060547623, y, x\right)\right)\\
\end{array}
\end{array}
if z < -1.6e25Initial program 7.2%
Taylor expanded in z around inf
+-commutativeN/A
lower-fma.f6481.3
Applied rewrites81.3%
if -1.6e25 < z < 2.25e19Initial program 99.0%
Taylor expanded in t around inf
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
unpow2N/A
lower-*.f6456.1
Applied rewrites56.1%
Taylor expanded in z around 0
Applied rewrites55.0%
Taylor expanded in t around 0
*-commutativeN/A
lower-*.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
unpow2N/A
lower-*.f64N/A
+-commutativeN/A
lower-fma.f6489.4
Applied rewrites89.4%
if 2.25e19 < z Initial program 8.6%
Taylor expanded in z around inf
associate--l+N/A
+-commutativeN/A
associate--l+N/A
+-commutativeN/A
distribute-rgt-out--N/A
metadata-evalN/A
metadata-evalN/A
metadata-evalN/A
times-fracN/A
distribute-rgt-out--N/A
*-commutativeN/A
mul-1-negN/A
distribute-neg-frac2N/A
mul-1-negN/A
associate-+l+N/A
Applied rewrites88.1%
Final simplification87.1%
(FPCore (x y z t a b)
:precision binary64
(if (<= z -4.05e+43)
(fma 3.13060547623 y x)
(if (<= z 1.85e+24)
(+ (/ (fma (fma (* z y) t (* a y)) z (* b y)) 0.607771387771) x)
(fma -36.52704169880642 (/ y z) (fma 3.13060547623 y x)))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (z <= -4.05e+43) {
tmp = fma(3.13060547623, y, x);
} else if (z <= 1.85e+24) {
tmp = (fma(fma((z * y), t, (a * y)), z, (b * y)) / 0.607771387771) + x;
} else {
tmp = fma(-36.52704169880642, (y / z), fma(3.13060547623, y, x));
}
return tmp;
}
function code(x, y, z, t, a, b) tmp = 0.0 if (z <= -4.05e+43) tmp = fma(3.13060547623, y, x); elseif (z <= 1.85e+24) tmp = Float64(Float64(fma(fma(Float64(z * y), t, Float64(a * y)), z, Float64(b * y)) / 0.607771387771) + x); else tmp = fma(-36.52704169880642, Float64(y / z), fma(3.13060547623, y, x)); end return tmp end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[z, -4.05e+43], N[(3.13060547623 * y + x), $MachinePrecision], If[LessEqual[z, 1.85e+24], N[(N[(N[(N[(N[(z * y), $MachinePrecision] * t + N[(a * y), $MachinePrecision]), $MachinePrecision] * z + N[(b * y), $MachinePrecision]), $MachinePrecision] / 0.607771387771), $MachinePrecision] + x), $MachinePrecision], N[(-36.52704169880642 * N[(y / z), $MachinePrecision] + N[(3.13060547623 * y + x), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;z \leq -4.05 \cdot 10^{+43}:\\
\;\;\;\;\mathsf{fma}\left(3.13060547623, y, x\right)\\
\mathbf{elif}\;z \leq 1.85 \cdot 10^{+24}:\\
\;\;\;\;\frac{\mathsf{fma}\left(\mathsf{fma}\left(z \cdot y, t, a \cdot y\right), z, b \cdot y\right)}{0.607771387771} + x\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(-36.52704169880642, \frac{y}{z}, \mathsf{fma}\left(3.13060547623, y, x\right)\right)\\
\end{array}
\end{array}
if z < -4.0499999999999998e43Initial program 2.6%
Taylor expanded in z around inf
+-commutativeN/A
lower-fma.f6484.6
Applied rewrites84.6%
if -4.0499999999999998e43 < z < 1.85e24Initial program 96.9%
Taylor expanded in z around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-*.f64N/A
lower-*.f6486.6
Applied rewrites86.6%
Taylor expanded in z around 0
Applied rewrites81.4%
if 1.85e24 < z Initial program 8.7%
Taylor expanded in z around inf
associate--l+N/A
+-commutativeN/A
associate--l+N/A
+-commutativeN/A
distribute-rgt-out--N/A
metadata-evalN/A
metadata-evalN/A
metadata-evalN/A
times-fracN/A
distribute-rgt-out--N/A
*-commutativeN/A
mul-1-negN/A
distribute-neg-frac2N/A
mul-1-negN/A
associate-+l+N/A
Applied rewrites89.5%
Final simplification84.1%
(FPCore (x y z t a b)
:precision binary64
(if (<= z -8.8e+64)
(fma 3.13060547623 y x)
(if (<= z 3.3e+22)
(+
(fma
(* 1.6453555072203998 y)
b
(* (* (fma 1.6453555072203998 a (* -32.324150453290734 b)) y) z))
x)
(fma -36.52704169880642 (/ y z) (fma 3.13060547623 y x)))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (z <= -8.8e+64) {
tmp = fma(3.13060547623, y, x);
} else if (z <= 3.3e+22) {
tmp = fma((1.6453555072203998 * y), b, ((fma(1.6453555072203998, a, (-32.324150453290734 * b)) * y) * z)) + x;
} else {
tmp = fma(-36.52704169880642, (y / z), fma(3.13060547623, y, x));
}
return tmp;
}
function code(x, y, z, t, a, b) tmp = 0.0 if (z <= -8.8e+64) tmp = fma(3.13060547623, y, x); elseif (z <= 3.3e+22) tmp = Float64(fma(Float64(1.6453555072203998 * y), b, Float64(Float64(fma(1.6453555072203998, a, Float64(-32.324150453290734 * b)) * y) * z)) + x); else tmp = fma(-36.52704169880642, Float64(y / z), fma(3.13060547623, y, x)); end return tmp end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[z, -8.8e+64], N[(3.13060547623 * y + x), $MachinePrecision], If[LessEqual[z, 3.3e+22], N[(N[(N[(1.6453555072203998 * y), $MachinePrecision] * b + N[(N[(N[(1.6453555072203998 * a + N[(-32.324150453290734 * b), $MachinePrecision]), $MachinePrecision] * y), $MachinePrecision] * z), $MachinePrecision]), $MachinePrecision] + x), $MachinePrecision], N[(-36.52704169880642 * N[(y / z), $MachinePrecision] + N[(3.13060547623 * y + x), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;z \leq -8.8 \cdot 10^{+64}:\\
\;\;\;\;\mathsf{fma}\left(3.13060547623, y, x\right)\\
\mathbf{elif}\;z \leq 3.3 \cdot 10^{+22}:\\
\;\;\;\;\mathsf{fma}\left(1.6453555072203998 \cdot y, b, \left(\mathsf{fma}\left(1.6453555072203998, a, -32.324150453290734 \cdot b\right) \cdot y\right) \cdot z\right) + x\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(-36.52704169880642, \frac{y}{z}, \mathsf{fma}\left(3.13060547623, y, x\right)\right)\\
\end{array}
\end{array}
if z < -8.80000000000000007e64Initial program 0.3%
Taylor expanded in z around inf
+-commutativeN/A
lower-fma.f6493.8
Applied rewrites93.8%
if -8.80000000000000007e64 < z < 3.2999999999999998e22Initial program 93.2%
Taylor expanded in z around 0
*-commutativeN/A
associate-*r*N/A
lower-fma.f64N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
associate-*r*N/A
associate-*r*N/A
distribute-rgt-out--N/A
lower-*.f64N/A
sub-negN/A
lower-fma.f64N/A
distribute-lft-neg-inN/A
lower-*.f64N/A
metadata-eval75.3
Applied rewrites75.3%
if 3.2999999999999998e22 < z Initial program 8.7%
Taylor expanded in z around inf
associate--l+N/A
+-commutativeN/A
associate--l+N/A
+-commutativeN/A
distribute-rgt-out--N/A
metadata-evalN/A
metadata-evalN/A
metadata-evalN/A
times-fracN/A
distribute-rgt-out--N/A
*-commutativeN/A
mul-1-negN/A
distribute-neg-frac2N/A
mul-1-negN/A
associate-+l+N/A
Applied rewrites89.5%
Final simplification82.4%
(FPCore (x y z t a b)
:precision binary64
(if (<= z -4.05e+43)
(fma 3.13060547623 y x)
(if (<= z 2.1e+22)
(fma
(* (fma 1.6453555072203998 a (* -32.324150453290734 b)) y)
z
(fma (* b y) 1.6453555072203998 x))
(fma -36.52704169880642 (/ y z) (fma 3.13060547623 y x)))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (z <= -4.05e+43) {
tmp = fma(3.13060547623, y, x);
} else if (z <= 2.1e+22) {
tmp = fma((fma(1.6453555072203998, a, (-32.324150453290734 * b)) * y), z, fma((b * y), 1.6453555072203998, x));
} else {
tmp = fma(-36.52704169880642, (y / z), fma(3.13060547623, y, x));
}
return tmp;
}
function code(x, y, z, t, a, b) tmp = 0.0 if (z <= -4.05e+43) tmp = fma(3.13060547623, y, x); elseif (z <= 2.1e+22) tmp = fma(Float64(fma(1.6453555072203998, a, Float64(-32.324150453290734 * b)) * y), z, fma(Float64(b * y), 1.6453555072203998, x)); else tmp = fma(-36.52704169880642, Float64(y / z), fma(3.13060547623, y, x)); end return tmp end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[z, -4.05e+43], N[(3.13060547623 * y + x), $MachinePrecision], If[LessEqual[z, 2.1e+22], N[(N[(N[(1.6453555072203998 * a + N[(-32.324150453290734 * b), $MachinePrecision]), $MachinePrecision] * y), $MachinePrecision] * z + N[(N[(b * y), $MachinePrecision] * 1.6453555072203998 + x), $MachinePrecision]), $MachinePrecision], N[(-36.52704169880642 * N[(y / z), $MachinePrecision] + N[(3.13060547623 * y + x), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;z \leq -4.05 \cdot 10^{+43}:\\
\;\;\;\;\mathsf{fma}\left(3.13060547623, y, x\right)\\
\mathbf{elif}\;z \leq 2.1 \cdot 10^{+22}:\\
\;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(1.6453555072203998, a, -32.324150453290734 \cdot b\right) \cdot y, z, \mathsf{fma}\left(b \cdot y, 1.6453555072203998, x\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(-36.52704169880642, \frac{y}{z}, \mathsf{fma}\left(3.13060547623, y, x\right)\right)\\
\end{array}
\end{array}
if z < -4.0499999999999998e43Initial program 2.6%
Taylor expanded in z around inf
+-commutativeN/A
lower-fma.f6484.6
Applied rewrites84.6%
if -4.0499999999999998e43 < z < 2.0999999999999998e22Initial program 96.9%
Taylor expanded in z around 0
associate-+r+N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
associate-*r*N/A
associate-*r*N/A
distribute-rgt-out--N/A
lower-*.f64N/A
sub-negN/A
lower-fma.f64N/A
distribute-lft-neg-inN/A
lower-*.f64N/A
metadata-evalN/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower-*.f6478.1
Applied rewrites78.1%
if 2.0999999999999998e22 < z Initial program 8.7%
Taylor expanded in z around inf
associate--l+N/A
+-commutativeN/A
associate--l+N/A
+-commutativeN/A
distribute-rgt-out--N/A
metadata-evalN/A
metadata-evalN/A
metadata-evalN/A
times-fracN/A
distribute-rgt-out--N/A
*-commutativeN/A
mul-1-negN/A
distribute-neg-frac2N/A
mul-1-negN/A
associate-+l+N/A
Applied rewrites89.5%
Final simplification82.3%
(FPCore (x y z t a b)
:precision binary64
(if (<= z -2.3e+53)
(fma 3.13060547623 y x)
(if (<= z -1.4e-111)
(+ (/ (* (* a z) y) 0.607771387771) x)
(if (<= z 7.6e+21)
(+ (* (fma (* b z) -32.324150453290734 (* 1.6453555072203998 b)) y) x)
(fma -36.52704169880642 (/ y z) (fma 3.13060547623 y x))))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (z <= -2.3e+53) {
tmp = fma(3.13060547623, y, x);
} else if (z <= -1.4e-111) {
tmp = (((a * z) * y) / 0.607771387771) + x;
} else if (z <= 7.6e+21) {
tmp = (fma((b * z), -32.324150453290734, (1.6453555072203998 * b)) * y) + x;
} else {
tmp = fma(-36.52704169880642, (y / z), fma(3.13060547623, y, x));
}
return tmp;
}
function code(x, y, z, t, a, b) tmp = 0.0 if (z <= -2.3e+53) tmp = fma(3.13060547623, y, x); elseif (z <= -1.4e-111) tmp = Float64(Float64(Float64(Float64(a * z) * y) / 0.607771387771) + x); elseif (z <= 7.6e+21) tmp = Float64(Float64(fma(Float64(b * z), -32.324150453290734, Float64(1.6453555072203998 * b)) * y) + x); else tmp = fma(-36.52704169880642, Float64(y / z), fma(3.13060547623, y, x)); end return tmp end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[z, -2.3e+53], N[(3.13060547623 * y + x), $MachinePrecision], If[LessEqual[z, -1.4e-111], N[(N[(N[(N[(a * z), $MachinePrecision] * y), $MachinePrecision] / 0.607771387771), $MachinePrecision] + x), $MachinePrecision], If[LessEqual[z, 7.6e+21], N[(N[(N[(N[(b * z), $MachinePrecision] * -32.324150453290734 + N[(1.6453555072203998 * b), $MachinePrecision]), $MachinePrecision] * y), $MachinePrecision] + x), $MachinePrecision], N[(-36.52704169880642 * N[(y / z), $MachinePrecision] + N[(3.13060547623 * y + x), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;z \leq -2.3 \cdot 10^{+53}:\\
\;\;\;\;\mathsf{fma}\left(3.13060547623, y, x\right)\\
\mathbf{elif}\;z \leq -1.4 \cdot 10^{-111}:\\
\;\;\;\;\frac{\left(a \cdot z\right) \cdot y}{0.607771387771} + x\\
\mathbf{elif}\;z \leq 7.6 \cdot 10^{+21}:\\
\;\;\;\;\mathsf{fma}\left(b \cdot z, -32.324150453290734, 1.6453555072203998 \cdot b\right) \cdot y + x\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(-36.52704169880642, \frac{y}{z}, \mathsf{fma}\left(3.13060547623, y, x\right)\right)\\
\end{array}
\end{array}
if z < -2.3000000000000002e53Initial program 0.4%
Taylor expanded in z around inf
+-commutativeN/A
lower-fma.f6492.1
Applied rewrites92.1%
if -2.3000000000000002e53 < z < -1.39999999999999998e-111Initial program 83.8%
Taylor expanded in a around inf
lower-*.f6460.4
Applied rewrites60.4%
Taylor expanded in z around 0
Applied rewrites60.2%
if -1.39999999999999998e-111 < z < 7.6e21Initial program 98.7%
Taylor expanded in b around inf
associate-*l/N/A
lower-*.f64N/A
lower-/.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower-+.f6477.5
Applied rewrites77.5%
Taylor expanded in z around 0
Applied rewrites76.6%
if 7.6e21 < z Initial program 8.7%
Taylor expanded in z around inf
associate--l+N/A
+-commutativeN/A
associate--l+N/A
+-commutativeN/A
distribute-rgt-out--N/A
metadata-evalN/A
metadata-evalN/A
metadata-evalN/A
times-fracN/A
distribute-rgt-out--N/A
*-commutativeN/A
mul-1-negN/A
distribute-neg-frac2N/A
mul-1-negN/A
associate-+l+N/A
Applied rewrites89.5%
Final simplification79.8%
(FPCore (x y z t a b)
:precision binary64
(if (<= z -2.3e+53)
(fma 3.13060547623 y x)
(if (<= z -3.5e-110)
(+ (/ (* (* z y) a) 0.607771387771) x)
(if (<= z 7.6e+21)
(+ (* (fma (* b z) -32.324150453290734 (* 1.6453555072203998 b)) y) x)
(fma -36.52704169880642 (/ y z) (fma 3.13060547623 y x))))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (z <= -2.3e+53) {
tmp = fma(3.13060547623, y, x);
} else if (z <= -3.5e-110) {
tmp = (((z * y) * a) / 0.607771387771) + x;
} else if (z <= 7.6e+21) {
tmp = (fma((b * z), -32.324150453290734, (1.6453555072203998 * b)) * y) + x;
} else {
tmp = fma(-36.52704169880642, (y / z), fma(3.13060547623, y, x));
}
return tmp;
}
function code(x, y, z, t, a, b) tmp = 0.0 if (z <= -2.3e+53) tmp = fma(3.13060547623, y, x); elseif (z <= -3.5e-110) tmp = Float64(Float64(Float64(Float64(z * y) * a) / 0.607771387771) + x); elseif (z <= 7.6e+21) tmp = Float64(Float64(fma(Float64(b * z), -32.324150453290734, Float64(1.6453555072203998 * b)) * y) + x); else tmp = fma(-36.52704169880642, Float64(y / z), fma(3.13060547623, y, x)); end return tmp end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[z, -2.3e+53], N[(3.13060547623 * y + x), $MachinePrecision], If[LessEqual[z, -3.5e-110], N[(N[(N[(N[(z * y), $MachinePrecision] * a), $MachinePrecision] / 0.607771387771), $MachinePrecision] + x), $MachinePrecision], If[LessEqual[z, 7.6e+21], N[(N[(N[(N[(b * z), $MachinePrecision] * -32.324150453290734 + N[(1.6453555072203998 * b), $MachinePrecision]), $MachinePrecision] * y), $MachinePrecision] + x), $MachinePrecision], N[(-36.52704169880642 * N[(y / z), $MachinePrecision] + N[(3.13060547623 * y + x), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;z \leq -2.3 \cdot 10^{+53}:\\
\;\;\;\;\mathsf{fma}\left(3.13060547623, y, x\right)\\
\mathbf{elif}\;z \leq -3.5 \cdot 10^{-110}:\\
\;\;\;\;\frac{\left(z \cdot y\right) \cdot a}{0.607771387771} + x\\
\mathbf{elif}\;z \leq 7.6 \cdot 10^{+21}:\\
\;\;\;\;\mathsf{fma}\left(b \cdot z, -32.324150453290734, 1.6453555072203998 \cdot b\right) \cdot y + x\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(-36.52704169880642, \frac{y}{z}, \mathsf{fma}\left(3.13060547623, y, x\right)\right)\\
\end{array}
\end{array}
if z < -2.3000000000000002e53Initial program 0.4%
Taylor expanded in z around inf
+-commutativeN/A
lower-fma.f6492.1
Applied rewrites92.1%
if -2.3000000000000002e53 < z < -3.49999999999999974e-110Initial program 83.8%
Taylor expanded in t around inf
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
unpow2N/A
lower-*.f6444.7
Applied rewrites44.7%
Taylor expanded in z around 0
Applied rewrites39.5%
Taylor expanded in a around inf
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f6458.3
Applied rewrites58.3%
if -3.49999999999999974e-110 < z < 7.6e21Initial program 98.7%
Taylor expanded in b around inf
associate-*l/N/A
lower-*.f64N/A
lower-/.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower-+.f6477.5
Applied rewrites77.5%
Taylor expanded in z around 0
Applied rewrites76.6%
if 7.6e21 < z Initial program 8.7%
Taylor expanded in z around inf
associate--l+N/A
+-commutativeN/A
associate--l+N/A
+-commutativeN/A
distribute-rgt-out--N/A
metadata-evalN/A
metadata-evalN/A
metadata-evalN/A
times-fracN/A
distribute-rgt-out--N/A
*-commutativeN/A
mul-1-negN/A
distribute-neg-frac2N/A
mul-1-negN/A
associate-+l+N/A
Applied rewrites89.5%
Final simplification79.5%
(FPCore (x y z t a b)
:precision binary64
(if (<= z -6.2e+69)
(fma 3.13060547623 y x)
(if (<= z -3.5e-110)
(+ (/ (* (* a y) z) 0.607771387771) x)
(if (<= z 7.6e+21)
(+ (* (fma (* b z) -32.324150453290734 (* 1.6453555072203998 b)) y) x)
(fma -36.52704169880642 (/ y z) (fma 3.13060547623 y x))))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (z <= -6.2e+69) {
tmp = fma(3.13060547623, y, x);
} else if (z <= -3.5e-110) {
tmp = (((a * y) * z) / 0.607771387771) + x;
} else if (z <= 7.6e+21) {
tmp = (fma((b * z), -32.324150453290734, (1.6453555072203998 * b)) * y) + x;
} else {
tmp = fma(-36.52704169880642, (y / z), fma(3.13060547623, y, x));
}
return tmp;
}
function code(x, y, z, t, a, b) tmp = 0.0 if (z <= -6.2e+69) tmp = fma(3.13060547623, y, x); elseif (z <= -3.5e-110) tmp = Float64(Float64(Float64(Float64(a * y) * z) / 0.607771387771) + x); elseif (z <= 7.6e+21) tmp = Float64(Float64(fma(Float64(b * z), -32.324150453290734, Float64(1.6453555072203998 * b)) * y) + x); else tmp = fma(-36.52704169880642, Float64(y / z), fma(3.13060547623, y, x)); end return tmp end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[z, -6.2e+69], N[(3.13060547623 * y + x), $MachinePrecision], If[LessEqual[z, -3.5e-110], N[(N[(N[(N[(a * y), $MachinePrecision] * z), $MachinePrecision] / 0.607771387771), $MachinePrecision] + x), $MachinePrecision], If[LessEqual[z, 7.6e+21], N[(N[(N[(N[(b * z), $MachinePrecision] * -32.324150453290734 + N[(1.6453555072203998 * b), $MachinePrecision]), $MachinePrecision] * y), $MachinePrecision] + x), $MachinePrecision], N[(-36.52704169880642 * N[(y / z), $MachinePrecision] + N[(3.13060547623 * y + x), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;z \leq -6.2 \cdot 10^{+69}:\\
\;\;\;\;\mathsf{fma}\left(3.13060547623, y, x\right)\\
\mathbf{elif}\;z \leq -3.5 \cdot 10^{-110}:\\
\;\;\;\;\frac{\left(a \cdot y\right) \cdot z}{0.607771387771} + x\\
\mathbf{elif}\;z \leq 7.6 \cdot 10^{+21}:\\
\;\;\;\;\mathsf{fma}\left(b \cdot z, -32.324150453290734, 1.6453555072203998 \cdot b\right) \cdot y + x\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(-36.52704169880642, \frac{y}{z}, \mathsf{fma}\left(3.13060547623, y, x\right)\right)\\
\end{array}
\end{array}
if z < -6.1999999999999997e69Initial program 0.3%
Taylor expanded in z around inf
+-commutativeN/A
lower-fma.f6493.7
Applied rewrites93.7%
if -6.1999999999999997e69 < z < -3.49999999999999974e-110Initial program 80.5%
Taylor expanded in t around inf
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
unpow2N/A
lower-*.f6445.0
Applied rewrites45.0%
Taylor expanded in z around 0
Applied rewrites40.0%
Applied rewrites36.2%
Taylor expanded in a around inf
associate-*r*N/A
lower-*.f64N/A
lower-*.f6454.3
Applied rewrites54.3%
if -3.49999999999999974e-110 < z < 7.6e21Initial program 98.7%
Taylor expanded in b around inf
associate-*l/N/A
lower-*.f64N/A
lower-/.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower-+.f6477.5
Applied rewrites77.5%
Taylor expanded in z around 0
Applied rewrites76.6%
if 7.6e21 < z Initial program 8.7%
Taylor expanded in z around inf
associate--l+N/A
+-commutativeN/A
associate--l+N/A
+-commutativeN/A
distribute-rgt-out--N/A
metadata-evalN/A
metadata-evalN/A
metadata-evalN/A
times-fracN/A
distribute-rgt-out--N/A
*-commutativeN/A
mul-1-negN/A
distribute-neg-frac2N/A
mul-1-negN/A
associate-+l+N/A
Applied rewrites89.5%
Final simplification78.8%
(FPCore (x y z t a b)
:precision binary64
(if (<= z -4.05e+43)
(fma 3.13060547623 y x)
(if (<= z 7.6e+21)
(+ (* (fma (* b z) -32.324150453290734 (* 1.6453555072203998 b)) y) x)
(fma -36.52704169880642 (/ y z) (fma 3.13060547623 y x)))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (z <= -4.05e+43) {
tmp = fma(3.13060547623, y, x);
} else if (z <= 7.6e+21) {
tmp = (fma((b * z), -32.324150453290734, (1.6453555072203998 * b)) * y) + x;
} else {
tmp = fma(-36.52704169880642, (y / z), fma(3.13060547623, y, x));
}
return tmp;
}
function code(x, y, z, t, a, b) tmp = 0.0 if (z <= -4.05e+43) tmp = fma(3.13060547623, y, x); elseif (z <= 7.6e+21) tmp = Float64(Float64(fma(Float64(b * z), -32.324150453290734, Float64(1.6453555072203998 * b)) * y) + x); else tmp = fma(-36.52704169880642, Float64(y / z), fma(3.13060547623, y, x)); end return tmp end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[z, -4.05e+43], N[(3.13060547623 * y + x), $MachinePrecision], If[LessEqual[z, 7.6e+21], N[(N[(N[(N[(b * z), $MachinePrecision] * -32.324150453290734 + N[(1.6453555072203998 * b), $MachinePrecision]), $MachinePrecision] * y), $MachinePrecision] + x), $MachinePrecision], N[(-36.52704169880642 * N[(y / z), $MachinePrecision] + N[(3.13060547623 * y + x), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;z \leq -4.05 \cdot 10^{+43}:\\
\;\;\;\;\mathsf{fma}\left(3.13060547623, y, x\right)\\
\mathbf{elif}\;z \leq 7.6 \cdot 10^{+21}:\\
\;\;\;\;\mathsf{fma}\left(b \cdot z, -32.324150453290734, 1.6453555072203998 \cdot b\right) \cdot y + x\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(-36.52704169880642, \frac{y}{z}, \mathsf{fma}\left(3.13060547623, y, x\right)\right)\\
\end{array}
\end{array}
if z < -4.0499999999999998e43Initial program 2.6%
Taylor expanded in z around inf
+-commutativeN/A
lower-fma.f6484.6
Applied rewrites84.6%
if -4.0499999999999998e43 < z < 7.6e21Initial program 96.9%
Taylor expanded in b around inf
associate-*l/N/A
lower-*.f64N/A
lower-/.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower-+.f6468.0
Applied rewrites68.0%
Taylor expanded in z around 0
Applied rewrites67.4%
if 7.6e21 < z Initial program 8.7%
Taylor expanded in z around inf
associate--l+N/A
+-commutativeN/A
associate--l+N/A
+-commutativeN/A
distribute-rgt-out--N/A
metadata-evalN/A
metadata-evalN/A
metadata-evalN/A
times-fracN/A
distribute-rgt-out--N/A
*-commutativeN/A
mul-1-negN/A
distribute-neg-frac2N/A
mul-1-negN/A
associate-+l+N/A
Applied rewrites89.5%
Final simplification76.5%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (fma -36.52704169880642 (/ y z) (fma 3.13060547623 y x))))
(if (<= z -70000000000000.0)
t_1
(if (<= z 1.06e+14) (fma (* 1.6453555072203998 b) y x) t_1))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = fma(-36.52704169880642, (y / z), fma(3.13060547623, y, x));
double tmp;
if (z <= -70000000000000.0) {
tmp = t_1;
} else if (z <= 1.06e+14) {
tmp = fma((1.6453555072203998 * b), y, x);
} else {
tmp = t_1;
}
return tmp;
}
function code(x, y, z, t, a, b) t_1 = fma(-36.52704169880642, Float64(y / z), fma(3.13060547623, y, x)) tmp = 0.0 if (z <= -70000000000000.0) tmp = t_1; elseif (z <= 1.06e+14) tmp = fma(Float64(1.6453555072203998 * b), y, x); else tmp = t_1; end return tmp end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(-36.52704169880642 * N[(y / z), $MachinePrecision] + N[(3.13060547623 * y + x), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[z, -70000000000000.0], t$95$1, If[LessEqual[z, 1.06e+14], N[(N[(1.6453555072203998 * b), $MachinePrecision] * y + x), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \mathsf{fma}\left(-36.52704169880642, \frac{y}{z}, \mathsf{fma}\left(3.13060547623, y, x\right)\right)\\
\mathbf{if}\;z \leq -70000000000000:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;z \leq 1.06 \cdot 10^{+14}:\\
\;\;\;\;\mathsf{fma}\left(1.6453555072203998 \cdot b, y, x\right)\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if z < -7e13 or 1.06e14 < z Initial program 11.4%
Taylor expanded in z around inf
associate--l+N/A
+-commutativeN/A
associate--l+N/A
+-commutativeN/A
distribute-rgt-out--N/A
metadata-evalN/A
metadata-evalN/A
metadata-evalN/A
times-fracN/A
distribute-rgt-out--N/A
*-commutativeN/A
mul-1-negN/A
distribute-neg-frac2N/A
mul-1-negN/A
associate-+l+N/A
Applied rewrites82.4%
if -7e13 < z < 1.06e14Initial program 99.7%
Taylor expanded in z around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower-*.f6470.3
Applied rewrites70.3%
Applied rewrites70.4%
(FPCore (x y z t a b)
:precision binary64
(if (<= z -70000000000000.0)
(fma 3.13060547623 y x)
(if (<= z 1.06e+14)
(fma (* 1.6453555072203998 b) y x)
(fma 3.13060547623 y x))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (z <= -70000000000000.0) {
tmp = fma(3.13060547623, y, x);
} else if (z <= 1.06e+14) {
tmp = fma((1.6453555072203998 * b), y, x);
} else {
tmp = fma(3.13060547623, y, x);
}
return tmp;
}
function code(x, y, z, t, a, b) tmp = 0.0 if (z <= -70000000000000.0) tmp = fma(3.13060547623, y, x); elseif (z <= 1.06e+14) tmp = fma(Float64(1.6453555072203998 * b), y, x); else tmp = fma(3.13060547623, y, x); end return tmp end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[z, -70000000000000.0], N[(3.13060547623 * y + x), $MachinePrecision], If[LessEqual[z, 1.06e+14], N[(N[(1.6453555072203998 * b), $MachinePrecision] * y + x), $MachinePrecision], N[(3.13060547623 * y + x), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;z \leq -70000000000000:\\
\;\;\;\;\mathsf{fma}\left(3.13060547623, y, x\right)\\
\mathbf{elif}\;z \leq 1.06 \cdot 10^{+14}:\\
\;\;\;\;\mathsf{fma}\left(1.6453555072203998 \cdot b, y, x\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(3.13060547623, y, x\right)\\
\end{array}
\end{array}
if z < -7e13 or 1.06e14 < z Initial program 11.4%
Taylor expanded in z around inf
+-commutativeN/A
lower-fma.f6482.3
Applied rewrites82.3%
if -7e13 < z < 1.06e14Initial program 99.7%
Taylor expanded in z around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower-*.f6470.3
Applied rewrites70.3%
Applied rewrites70.4%
(FPCore (x y z t a b) :precision binary64 (fma 3.13060547623 y x))
double code(double x, double y, double z, double t, double a, double b) {
return fma(3.13060547623, y, x);
}
function code(x, y, z, t, a, b) return fma(3.13060547623, y, x) end
code[x_, y_, z_, t_, a_, b_] := N[(3.13060547623 * y + x), $MachinePrecision]
\begin{array}{l}
\\
\mathsf{fma}\left(3.13060547623, y, x\right)
\end{array}
Initial program 54.5%
Taylor expanded in z around inf
+-commutativeN/A
lower-fma.f6462.2
Applied rewrites62.2%
(FPCore (x y z t a b) :precision binary64 (* 3.13060547623 y))
double code(double x, double y, double z, double t, double a, double b) {
return 3.13060547623 * y;
}
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 = 3.13060547623d0 * y
end function
public static double code(double x, double y, double z, double t, double a, double b) {
return 3.13060547623 * y;
}
def code(x, y, z, t, a, b): return 3.13060547623 * y
function code(x, y, z, t, a, b) return Float64(3.13060547623 * y) end
function tmp = code(x, y, z, t, a, b) tmp = 3.13060547623 * y; end
code[x_, y_, z_, t_, a_, b_] := N[(3.13060547623 * y), $MachinePrecision]
\begin{array}{l}
\\
3.13060547623 \cdot y
\end{array}
Initial program 54.5%
Taylor expanded in z around inf
+-commutativeN/A
lower-fma.f6462.2
Applied rewrites62.2%
Taylor expanded in x around 0
Applied rewrites25.0%
(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 2024276
(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))))