
(FPCore (x1 x2)
:precision binary64
(let* ((t_0 (* (* 3.0 x1) x1))
(t_1 (+ (* x1 x1) 1.0))
(t_2 (/ (- (+ t_0 (* 2.0 x2)) x1) t_1)))
(+
x1
(+
(+
(+
(+
(*
(+
(* (* (* 2.0 x1) t_2) (- t_2 3.0))
(* (* x1 x1) (- (* 4.0 t_2) 6.0)))
t_1)
(* t_0 t_2))
(* (* x1 x1) x1))
x1)
(* 3.0 (/ (- (- t_0 (* 2.0 x2)) x1) t_1))))))
double code(double x1, double x2) {
double t_0 = (3.0 * x1) * x1;
double t_1 = (x1 * x1) + 1.0;
double t_2 = ((t_0 + (2.0 * x2)) - x1) / t_1;
return x1 + (((((((((2.0 * x1) * t_2) * (t_2 - 3.0)) + ((x1 * x1) * ((4.0 * t_2) - 6.0))) * t_1) + (t_0 * t_2)) + ((x1 * x1) * x1)) + x1) + (3.0 * (((t_0 - (2.0 * x2)) - x1) / t_1)));
}
real(8) function code(x1, x2)
real(8), intent (in) :: x1
real(8), intent (in) :: x2
real(8) :: t_0
real(8) :: t_1
real(8) :: t_2
t_0 = (3.0d0 * x1) * x1
t_1 = (x1 * x1) + 1.0d0
t_2 = ((t_0 + (2.0d0 * x2)) - x1) / t_1
code = x1 + (((((((((2.0d0 * x1) * t_2) * (t_2 - 3.0d0)) + ((x1 * x1) * ((4.0d0 * t_2) - 6.0d0))) * t_1) + (t_0 * t_2)) + ((x1 * x1) * x1)) + x1) + (3.0d0 * (((t_0 - (2.0d0 * x2)) - x1) / t_1)))
end function
public static double code(double x1, double x2) {
double t_0 = (3.0 * x1) * x1;
double t_1 = (x1 * x1) + 1.0;
double t_2 = ((t_0 + (2.0 * x2)) - x1) / t_1;
return x1 + (((((((((2.0 * x1) * t_2) * (t_2 - 3.0)) + ((x1 * x1) * ((4.0 * t_2) - 6.0))) * t_1) + (t_0 * t_2)) + ((x1 * x1) * x1)) + x1) + (3.0 * (((t_0 - (2.0 * x2)) - x1) / t_1)));
}
def code(x1, x2): t_0 = (3.0 * x1) * x1 t_1 = (x1 * x1) + 1.0 t_2 = ((t_0 + (2.0 * x2)) - x1) / t_1 return x1 + (((((((((2.0 * x1) * t_2) * (t_2 - 3.0)) + ((x1 * x1) * ((4.0 * t_2) - 6.0))) * t_1) + (t_0 * t_2)) + ((x1 * x1) * x1)) + x1) + (3.0 * (((t_0 - (2.0 * x2)) - x1) / t_1)))
function code(x1, x2) t_0 = Float64(Float64(3.0 * x1) * x1) t_1 = Float64(Float64(x1 * x1) + 1.0) t_2 = Float64(Float64(Float64(t_0 + Float64(2.0 * x2)) - x1) / t_1) return Float64(x1 + Float64(Float64(Float64(Float64(Float64(Float64(Float64(Float64(Float64(2.0 * x1) * t_2) * Float64(t_2 - 3.0)) + Float64(Float64(x1 * x1) * Float64(Float64(4.0 * t_2) - 6.0))) * t_1) + Float64(t_0 * t_2)) + Float64(Float64(x1 * x1) * x1)) + x1) + Float64(3.0 * Float64(Float64(Float64(t_0 - Float64(2.0 * x2)) - x1) / t_1)))) end
function tmp = code(x1, x2) t_0 = (3.0 * x1) * x1; t_1 = (x1 * x1) + 1.0; t_2 = ((t_0 + (2.0 * x2)) - x1) / t_1; tmp = x1 + (((((((((2.0 * x1) * t_2) * (t_2 - 3.0)) + ((x1 * x1) * ((4.0 * t_2) - 6.0))) * t_1) + (t_0 * t_2)) + ((x1 * x1) * x1)) + x1) + (3.0 * (((t_0 - (2.0 * x2)) - x1) / t_1))); end
code[x1_, x2_] := Block[{t$95$0 = N[(N[(3.0 * x1), $MachinePrecision] * x1), $MachinePrecision]}, Block[{t$95$1 = N[(N[(x1 * x1), $MachinePrecision] + 1.0), $MachinePrecision]}, Block[{t$95$2 = N[(N[(N[(t$95$0 + N[(2.0 * x2), $MachinePrecision]), $MachinePrecision] - x1), $MachinePrecision] / t$95$1), $MachinePrecision]}, N[(x1 + N[(N[(N[(N[(N[(N[(N[(N[(N[(2.0 * x1), $MachinePrecision] * t$95$2), $MachinePrecision] * N[(t$95$2 - 3.0), $MachinePrecision]), $MachinePrecision] + N[(N[(x1 * x1), $MachinePrecision] * N[(N[(4.0 * t$95$2), $MachinePrecision] - 6.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * t$95$1), $MachinePrecision] + N[(t$95$0 * t$95$2), $MachinePrecision]), $MachinePrecision] + N[(N[(x1 * x1), $MachinePrecision] * x1), $MachinePrecision]), $MachinePrecision] + x1), $MachinePrecision] + N[(3.0 * N[(N[(N[(t$95$0 - N[(2.0 * x2), $MachinePrecision]), $MachinePrecision] - x1), $MachinePrecision] / t$95$1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \left(3 \cdot x1\right) \cdot x1\\
t_1 := x1 \cdot x1 + 1\\
t_2 := \frac{\left(t\_0 + 2 \cdot x2\right) - x1}{t\_1}\\
x1 + \left(\left(\left(\left(\left(\left(\left(2 \cdot x1\right) \cdot t\_2\right) \cdot \left(t\_2 - 3\right) + \left(x1 \cdot x1\right) \cdot \left(4 \cdot t\_2 - 6\right)\right) \cdot t\_1 + t\_0 \cdot t\_2\right) + \left(x1 \cdot x1\right) \cdot x1\right) + x1\right) + 3 \cdot \frac{\left(t\_0 - 2 \cdot x2\right) - x1}{t\_1}\right)
\end{array}
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 16 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (x1 x2)
:precision binary64
(let* ((t_0 (* (* 3.0 x1) x1))
(t_1 (+ (* x1 x1) 1.0))
(t_2 (/ (- (+ t_0 (* 2.0 x2)) x1) t_1)))
(+
x1
(+
(+
(+
(+
(*
(+
(* (* (* 2.0 x1) t_2) (- t_2 3.0))
(* (* x1 x1) (- (* 4.0 t_2) 6.0)))
t_1)
(* t_0 t_2))
(* (* x1 x1) x1))
x1)
(* 3.0 (/ (- (- t_0 (* 2.0 x2)) x1) t_1))))))
double code(double x1, double x2) {
double t_0 = (3.0 * x1) * x1;
double t_1 = (x1 * x1) + 1.0;
double t_2 = ((t_0 + (2.0 * x2)) - x1) / t_1;
return x1 + (((((((((2.0 * x1) * t_2) * (t_2 - 3.0)) + ((x1 * x1) * ((4.0 * t_2) - 6.0))) * t_1) + (t_0 * t_2)) + ((x1 * x1) * x1)) + x1) + (3.0 * (((t_0 - (2.0 * x2)) - x1) / t_1)));
}
real(8) function code(x1, x2)
real(8), intent (in) :: x1
real(8), intent (in) :: x2
real(8) :: t_0
real(8) :: t_1
real(8) :: t_2
t_0 = (3.0d0 * x1) * x1
t_1 = (x1 * x1) + 1.0d0
t_2 = ((t_0 + (2.0d0 * x2)) - x1) / t_1
code = x1 + (((((((((2.0d0 * x1) * t_2) * (t_2 - 3.0d0)) + ((x1 * x1) * ((4.0d0 * t_2) - 6.0d0))) * t_1) + (t_0 * t_2)) + ((x1 * x1) * x1)) + x1) + (3.0d0 * (((t_0 - (2.0d0 * x2)) - x1) / t_1)))
end function
public static double code(double x1, double x2) {
double t_0 = (3.0 * x1) * x1;
double t_1 = (x1 * x1) + 1.0;
double t_2 = ((t_0 + (2.0 * x2)) - x1) / t_1;
return x1 + (((((((((2.0 * x1) * t_2) * (t_2 - 3.0)) + ((x1 * x1) * ((4.0 * t_2) - 6.0))) * t_1) + (t_0 * t_2)) + ((x1 * x1) * x1)) + x1) + (3.0 * (((t_0 - (2.0 * x2)) - x1) / t_1)));
}
def code(x1, x2): t_0 = (3.0 * x1) * x1 t_1 = (x1 * x1) + 1.0 t_2 = ((t_0 + (2.0 * x2)) - x1) / t_1 return x1 + (((((((((2.0 * x1) * t_2) * (t_2 - 3.0)) + ((x1 * x1) * ((4.0 * t_2) - 6.0))) * t_1) + (t_0 * t_2)) + ((x1 * x1) * x1)) + x1) + (3.0 * (((t_0 - (2.0 * x2)) - x1) / t_1)))
function code(x1, x2) t_0 = Float64(Float64(3.0 * x1) * x1) t_1 = Float64(Float64(x1 * x1) + 1.0) t_2 = Float64(Float64(Float64(t_0 + Float64(2.0 * x2)) - x1) / t_1) return Float64(x1 + Float64(Float64(Float64(Float64(Float64(Float64(Float64(Float64(Float64(2.0 * x1) * t_2) * Float64(t_2 - 3.0)) + Float64(Float64(x1 * x1) * Float64(Float64(4.0 * t_2) - 6.0))) * t_1) + Float64(t_0 * t_2)) + Float64(Float64(x1 * x1) * x1)) + x1) + Float64(3.0 * Float64(Float64(Float64(t_0 - Float64(2.0 * x2)) - x1) / t_1)))) end
function tmp = code(x1, x2) t_0 = (3.0 * x1) * x1; t_1 = (x1 * x1) + 1.0; t_2 = ((t_0 + (2.0 * x2)) - x1) / t_1; tmp = x1 + (((((((((2.0 * x1) * t_2) * (t_2 - 3.0)) + ((x1 * x1) * ((4.0 * t_2) - 6.0))) * t_1) + (t_0 * t_2)) + ((x1 * x1) * x1)) + x1) + (3.0 * (((t_0 - (2.0 * x2)) - x1) / t_1))); end
code[x1_, x2_] := Block[{t$95$0 = N[(N[(3.0 * x1), $MachinePrecision] * x1), $MachinePrecision]}, Block[{t$95$1 = N[(N[(x1 * x1), $MachinePrecision] + 1.0), $MachinePrecision]}, Block[{t$95$2 = N[(N[(N[(t$95$0 + N[(2.0 * x2), $MachinePrecision]), $MachinePrecision] - x1), $MachinePrecision] / t$95$1), $MachinePrecision]}, N[(x1 + N[(N[(N[(N[(N[(N[(N[(N[(N[(2.0 * x1), $MachinePrecision] * t$95$2), $MachinePrecision] * N[(t$95$2 - 3.0), $MachinePrecision]), $MachinePrecision] + N[(N[(x1 * x1), $MachinePrecision] * N[(N[(4.0 * t$95$2), $MachinePrecision] - 6.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * t$95$1), $MachinePrecision] + N[(t$95$0 * t$95$2), $MachinePrecision]), $MachinePrecision] + N[(N[(x1 * x1), $MachinePrecision] * x1), $MachinePrecision]), $MachinePrecision] + x1), $MachinePrecision] + N[(3.0 * N[(N[(N[(t$95$0 - N[(2.0 * x2), $MachinePrecision]), $MachinePrecision] - x1), $MachinePrecision] / t$95$1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \left(3 \cdot x1\right) \cdot x1\\
t_1 := x1 \cdot x1 + 1\\
t_2 := \frac{\left(t\_0 + 2 \cdot x2\right) - x1}{t\_1}\\
x1 + \left(\left(\left(\left(\left(\left(\left(2 \cdot x1\right) \cdot t\_2\right) \cdot \left(t\_2 - 3\right) + \left(x1 \cdot x1\right) \cdot \left(4 \cdot t\_2 - 6\right)\right) \cdot t\_1 + t\_0 \cdot t\_2\right) + \left(x1 \cdot x1\right) \cdot x1\right) + x1\right) + 3 \cdot \frac{\left(t\_0 - 2 \cdot x2\right) - x1}{t\_1}\right)
\end{array}
\end{array}
(FPCore (x1 x2)
:precision binary64
(let* ((t_0 (* (* 3.0 x1) x1))
(t_1 (+ (* x1 x1) 1.0))
(t_2 (/ (- (+ t_0 (* 2.0 x2)) x1) t_1))
(t_3 (/ (- (fma x2 2.0 t_0) x1) (fma x1 x1 1.0))))
(if (<=
(+
x1
(+
(+
(+
(+
(*
(+
(* (* (* 2.0 x1) t_2) (- t_2 3.0))
(* (* x1 x1) (- (* 4.0 t_2) 6.0)))
t_1)
(* t_0 t_2))
(* (* x1 x1) x1))
x1)
(* 3.0 (/ (- (- t_0 (* 2.0 x2)) x1) t_1))))
INFINITY)
(+
x1
(fma
(* x1 x1)
x1
(+
(fma
(fma (fma 4.0 t_3 -6.0) (* x1 x1) (* (- t_3 3.0) (* t_3 (* x1 2.0))))
(fma x1 x1 1.0)
(* 3.0 t_0))
(fma (fma -2.0 x2 (- x1)) 3.0 x1))))
(+ x1 (* (fma (fma 6.0 x1 -3.0) x1 -3.0) (* x1 x1))))))
double code(double x1, double x2) {
double t_0 = (3.0 * x1) * x1;
double t_1 = (x1 * x1) + 1.0;
double t_2 = ((t_0 + (2.0 * x2)) - x1) / t_1;
double t_3 = (fma(x2, 2.0, t_0) - x1) / fma(x1, x1, 1.0);
double tmp;
if ((x1 + (((((((((2.0 * x1) * t_2) * (t_2 - 3.0)) + ((x1 * x1) * ((4.0 * t_2) - 6.0))) * t_1) + (t_0 * t_2)) + ((x1 * x1) * x1)) + x1) + (3.0 * (((t_0 - (2.0 * x2)) - x1) / t_1)))) <= ((double) INFINITY)) {
tmp = x1 + fma((x1 * x1), x1, (fma(fma(fma(4.0, t_3, -6.0), (x1 * x1), ((t_3 - 3.0) * (t_3 * (x1 * 2.0)))), fma(x1, x1, 1.0), (3.0 * t_0)) + fma(fma(-2.0, x2, -x1), 3.0, x1)));
} else {
tmp = x1 + (fma(fma(6.0, x1, -3.0), x1, -3.0) * (x1 * x1));
}
return tmp;
}
function code(x1, x2) t_0 = Float64(Float64(3.0 * x1) * x1) t_1 = Float64(Float64(x1 * x1) + 1.0) t_2 = Float64(Float64(Float64(t_0 + Float64(2.0 * x2)) - x1) / t_1) t_3 = Float64(Float64(fma(x2, 2.0, t_0) - x1) / fma(x1, x1, 1.0)) tmp = 0.0 if (Float64(x1 + Float64(Float64(Float64(Float64(Float64(Float64(Float64(Float64(Float64(2.0 * x1) * t_2) * Float64(t_2 - 3.0)) + Float64(Float64(x1 * x1) * Float64(Float64(4.0 * t_2) - 6.0))) * t_1) + Float64(t_0 * t_2)) + Float64(Float64(x1 * x1) * x1)) + x1) + Float64(3.0 * Float64(Float64(Float64(t_0 - Float64(2.0 * x2)) - x1) / t_1)))) <= Inf) tmp = Float64(x1 + fma(Float64(x1 * x1), x1, Float64(fma(fma(fma(4.0, t_3, -6.0), Float64(x1 * x1), Float64(Float64(t_3 - 3.0) * Float64(t_3 * Float64(x1 * 2.0)))), fma(x1, x1, 1.0), Float64(3.0 * t_0)) + fma(fma(-2.0, x2, Float64(-x1)), 3.0, x1)))); else tmp = Float64(x1 + Float64(fma(fma(6.0, x1, -3.0), x1, -3.0) * Float64(x1 * x1))); end return tmp end
code[x1_, x2_] := Block[{t$95$0 = N[(N[(3.0 * x1), $MachinePrecision] * x1), $MachinePrecision]}, Block[{t$95$1 = N[(N[(x1 * x1), $MachinePrecision] + 1.0), $MachinePrecision]}, Block[{t$95$2 = N[(N[(N[(t$95$0 + N[(2.0 * x2), $MachinePrecision]), $MachinePrecision] - x1), $MachinePrecision] / t$95$1), $MachinePrecision]}, Block[{t$95$3 = N[(N[(N[(x2 * 2.0 + t$95$0), $MachinePrecision] - x1), $MachinePrecision] / N[(x1 * x1 + 1.0), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[N[(x1 + N[(N[(N[(N[(N[(N[(N[(N[(N[(2.0 * x1), $MachinePrecision] * t$95$2), $MachinePrecision] * N[(t$95$2 - 3.0), $MachinePrecision]), $MachinePrecision] + N[(N[(x1 * x1), $MachinePrecision] * N[(N[(4.0 * t$95$2), $MachinePrecision] - 6.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * t$95$1), $MachinePrecision] + N[(t$95$0 * t$95$2), $MachinePrecision]), $MachinePrecision] + N[(N[(x1 * x1), $MachinePrecision] * x1), $MachinePrecision]), $MachinePrecision] + x1), $MachinePrecision] + N[(3.0 * N[(N[(N[(t$95$0 - N[(2.0 * x2), $MachinePrecision]), $MachinePrecision] - x1), $MachinePrecision] / t$95$1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], Infinity], N[(x1 + N[(N[(x1 * x1), $MachinePrecision] * x1 + N[(N[(N[(N[(4.0 * t$95$3 + -6.0), $MachinePrecision] * N[(x1 * x1), $MachinePrecision] + N[(N[(t$95$3 - 3.0), $MachinePrecision] * N[(t$95$3 * N[(x1 * 2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(x1 * x1 + 1.0), $MachinePrecision] + N[(3.0 * t$95$0), $MachinePrecision]), $MachinePrecision] + N[(N[(-2.0 * x2 + (-x1)), $MachinePrecision] * 3.0 + x1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x1 + N[(N[(N[(6.0 * x1 + -3.0), $MachinePrecision] * x1 + -3.0), $MachinePrecision] * N[(x1 * x1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \left(3 \cdot x1\right) \cdot x1\\
t_1 := x1 \cdot x1 + 1\\
t_2 := \frac{\left(t\_0 + 2 \cdot x2\right) - x1}{t\_1}\\
t_3 := \frac{\mathsf{fma}\left(x2, 2, t\_0\right) - x1}{\mathsf{fma}\left(x1, x1, 1\right)}\\
\mathbf{if}\;x1 + \left(\left(\left(\left(\left(\left(\left(2 \cdot x1\right) \cdot t\_2\right) \cdot \left(t\_2 - 3\right) + \left(x1 \cdot x1\right) \cdot \left(4 \cdot t\_2 - 6\right)\right) \cdot t\_1 + t\_0 \cdot t\_2\right) + \left(x1 \cdot x1\right) \cdot x1\right) + x1\right) + 3 \cdot \frac{\left(t\_0 - 2 \cdot x2\right) - x1}{t\_1}\right) \leq \infty:\\
\;\;\;\;x1 + \mathsf{fma}\left(x1 \cdot x1, x1, \mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(4, t\_3, -6\right), x1 \cdot x1, \left(t\_3 - 3\right) \cdot \left(t\_3 \cdot \left(x1 \cdot 2\right)\right)\right), \mathsf{fma}\left(x1, x1, 1\right), 3 \cdot t\_0\right) + \mathsf{fma}\left(\mathsf{fma}\left(-2, x2, -x1\right), 3, x1\right)\right)\\
\mathbf{else}:\\
\;\;\;\;x1 + \mathsf{fma}\left(\mathsf{fma}\left(6, x1, -3\right), x1, -3\right) \cdot \left(x1 \cdot x1\right)\\
\end{array}
\end{array}
if (+.f64 x1 (+.f64 (+.f64 (+.f64 (+.f64 (*.f64 (+.f64 (*.f64 (*.f64 (*.f64 #s(literal 2 binary64) x1) (/.f64 (-.f64 (+.f64 (*.f64 (*.f64 #s(literal 3 binary64) x1) x1) (*.f64 #s(literal 2 binary64) x2)) x1) (+.f64 (*.f64 x1 x1) #s(literal 1 binary64)))) (-.f64 (/.f64 (-.f64 (+.f64 (*.f64 (*.f64 #s(literal 3 binary64) x1) x1) (*.f64 #s(literal 2 binary64) x2)) x1) (+.f64 (*.f64 x1 x1) #s(literal 1 binary64))) #s(literal 3 binary64))) (*.f64 (*.f64 x1 x1) (-.f64 (*.f64 #s(literal 4 binary64) (/.f64 (-.f64 (+.f64 (*.f64 (*.f64 #s(literal 3 binary64) x1) x1) (*.f64 #s(literal 2 binary64) x2)) x1) (+.f64 (*.f64 x1 x1) #s(literal 1 binary64)))) #s(literal 6 binary64)))) (+.f64 (*.f64 x1 x1) #s(literal 1 binary64))) (*.f64 (*.f64 (*.f64 #s(literal 3 binary64) x1) x1) (/.f64 (-.f64 (+.f64 (*.f64 (*.f64 #s(literal 3 binary64) x1) x1) (*.f64 #s(literal 2 binary64) x2)) x1) (+.f64 (*.f64 x1 x1) #s(literal 1 binary64))))) (*.f64 (*.f64 x1 x1) x1)) x1) (*.f64 #s(literal 3 binary64) (/.f64 (-.f64 (-.f64 (*.f64 (*.f64 #s(literal 3 binary64) x1) x1) (*.f64 #s(literal 2 binary64) x2)) x1) (+.f64 (*.f64 x1 x1) #s(literal 1 binary64)))))) < +inf.0Initial program 99.6%
Applied rewrites99.8%
Taylor expanded in x1 around inf
Applied rewrites99.1%
Taylor expanded in x1 around 0
lower-fma.f64N/A
mul-1-negN/A
lower-neg.f6499.8
Applied rewrites99.8%
if +inf.0 < (+.f64 x1 (+.f64 (+.f64 (+.f64 (+.f64 (*.f64 (+.f64 (*.f64 (*.f64 (*.f64 #s(literal 2 binary64) x1) (/.f64 (-.f64 (+.f64 (*.f64 (*.f64 #s(literal 3 binary64) x1) x1) (*.f64 #s(literal 2 binary64) x2)) x1) (+.f64 (*.f64 x1 x1) #s(literal 1 binary64)))) (-.f64 (/.f64 (-.f64 (+.f64 (*.f64 (*.f64 #s(literal 3 binary64) x1) x1) (*.f64 #s(literal 2 binary64) x2)) x1) (+.f64 (*.f64 x1 x1) #s(literal 1 binary64))) #s(literal 3 binary64))) (*.f64 (*.f64 x1 x1) (-.f64 (*.f64 #s(literal 4 binary64) (/.f64 (-.f64 (+.f64 (*.f64 (*.f64 #s(literal 3 binary64) x1) x1) (*.f64 #s(literal 2 binary64) x2)) x1) (+.f64 (*.f64 x1 x1) #s(literal 1 binary64)))) #s(literal 6 binary64)))) (+.f64 (*.f64 x1 x1) #s(literal 1 binary64))) (*.f64 (*.f64 (*.f64 #s(literal 3 binary64) x1) x1) (/.f64 (-.f64 (+.f64 (*.f64 (*.f64 #s(literal 3 binary64) x1) x1) (*.f64 #s(literal 2 binary64) x2)) x1) (+.f64 (*.f64 x1 x1) #s(literal 1 binary64))))) (*.f64 (*.f64 x1 x1) x1)) x1) (*.f64 #s(literal 3 binary64) (/.f64 (-.f64 (-.f64 (*.f64 (*.f64 #s(literal 3 binary64) x1) x1) (*.f64 #s(literal 2 binary64) x2)) x1) (+.f64 (*.f64 x1 x1) #s(literal 1 binary64)))))) Initial program 0.0%
Applied rewrites13.5%
Taylor expanded in x1 around -inf
*-commutativeN/A
lower-*.f64N/A
Applied rewrites100.0%
Taylor expanded in x1 around 0
Applied rewrites100.0%
Taylor expanded in x2 around 0
Applied rewrites100.0%
(FPCore (x1 x2)
:precision binary64
(let* ((t_0 (* (* 3.0 x1) x1))
(t_1 (+ (* x1 x1) 1.0))
(t_2 (/ (- (+ t_0 (* 2.0 x2)) x1) t_1))
(t_3
(+
x1
(+
(+
(+
(+
(*
(+
(* (* (* 2.0 x1) t_2) (- t_2 3.0))
(* (* x1 x1) (- (* 4.0 t_2) 6.0)))
t_1)
(* t_0 t_2))
(* (* x1 x1) x1))
x1)
(* 3.0 (/ (- (- t_0 (* 2.0 x2)) x1) t_1)))))
(t_4 (+ x1 (* (* (* x2 x2) x1) 8.0))))
(if (<= t_3 -1e+229)
t_4
(if (<= t_3 4e+177)
(+ x1 (* -6.0 x2))
(if (<= t_3 INFINITY) t_4 (+ x1 (* (* (* x1 x1) x2) 8.0)))))))
double code(double x1, double x2) {
double t_0 = (3.0 * x1) * x1;
double t_1 = (x1 * x1) + 1.0;
double t_2 = ((t_0 + (2.0 * x2)) - x1) / t_1;
double t_3 = x1 + (((((((((2.0 * x1) * t_2) * (t_2 - 3.0)) + ((x1 * x1) * ((4.0 * t_2) - 6.0))) * t_1) + (t_0 * t_2)) + ((x1 * x1) * x1)) + x1) + (3.0 * (((t_0 - (2.0 * x2)) - x1) / t_1)));
double t_4 = x1 + (((x2 * x2) * x1) * 8.0);
double tmp;
if (t_3 <= -1e+229) {
tmp = t_4;
} else if (t_3 <= 4e+177) {
tmp = x1 + (-6.0 * x2);
} else if (t_3 <= ((double) INFINITY)) {
tmp = t_4;
} else {
tmp = x1 + (((x1 * x1) * x2) * 8.0);
}
return tmp;
}
public static double code(double x1, double x2) {
double t_0 = (3.0 * x1) * x1;
double t_1 = (x1 * x1) + 1.0;
double t_2 = ((t_0 + (2.0 * x2)) - x1) / t_1;
double t_3 = x1 + (((((((((2.0 * x1) * t_2) * (t_2 - 3.0)) + ((x1 * x1) * ((4.0 * t_2) - 6.0))) * t_1) + (t_0 * t_2)) + ((x1 * x1) * x1)) + x1) + (3.0 * (((t_0 - (2.0 * x2)) - x1) / t_1)));
double t_4 = x1 + (((x2 * x2) * x1) * 8.0);
double tmp;
if (t_3 <= -1e+229) {
tmp = t_4;
} else if (t_3 <= 4e+177) {
tmp = x1 + (-6.0 * x2);
} else if (t_3 <= Double.POSITIVE_INFINITY) {
tmp = t_4;
} else {
tmp = x1 + (((x1 * x1) * x2) * 8.0);
}
return tmp;
}
def code(x1, x2): t_0 = (3.0 * x1) * x1 t_1 = (x1 * x1) + 1.0 t_2 = ((t_0 + (2.0 * x2)) - x1) / t_1 t_3 = x1 + (((((((((2.0 * x1) * t_2) * (t_2 - 3.0)) + ((x1 * x1) * ((4.0 * t_2) - 6.0))) * t_1) + (t_0 * t_2)) + ((x1 * x1) * x1)) + x1) + (3.0 * (((t_0 - (2.0 * x2)) - x1) / t_1))) t_4 = x1 + (((x2 * x2) * x1) * 8.0) tmp = 0 if t_3 <= -1e+229: tmp = t_4 elif t_3 <= 4e+177: tmp = x1 + (-6.0 * x2) elif t_3 <= math.inf: tmp = t_4 else: tmp = x1 + (((x1 * x1) * x2) * 8.0) return tmp
function code(x1, x2) t_0 = Float64(Float64(3.0 * x1) * x1) t_1 = Float64(Float64(x1 * x1) + 1.0) t_2 = Float64(Float64(Float64(t_0 + Float64(2.0 * x2)) - x1) / t_1) t_3 = Float64(x1 + Float64(Float64(Float64(Float64(Float64(Float64(Float64(Float64(Float64(2.0 * x1) * t_2) * Float64(t_2 - 3.0)) + Float64(Float64(x1 * x1) * Float64(Float64(4.0 * t_2) - 6.0))) * t_1) + Float64(t_0 * t_2)) + Float64(Float64(x1 * x1) * x1)) + x1) + Float64(3.0 * Float64(Float64(Float64(t_0 - Float64(2.0 * x2)) - x1) / t_1)))) t_4 = Float64(x1 + Float64(Float64(Float64(x2 * x2) * x1) * 8.0)) tmp = 0.0 if (t_3 <= -1e+229) tmp = t_4; elseif (t_3 <= 4e+177) tmp = Float64(x1 + Float64(-6.0 * x2)); elseif (t_3 <= Inf) tmp = t_4; else tmp = Float64(x1 + Float64(Float64(Float64(x1 * x1) * x2) * 8.0)); end return tmp end
function tmp_2 = code(x1, x2) t_0 = (3.0 * x1) * x1; t_1 = (x1 * x1) + 1.0; t_2 = ((t_0 + (2.0 * x2)) - x1) / t_1; t_3 = x1 + (((((((((2.0 * x1) * t_2) * (t_2 - 3.0)) + ((x1 * x1) * ((4.0 * t_2) - 6.0))) * t_1) + (t_0 * t_2)) + ((x1 * x1) * x1)) + x1) + (3.0 * (((t_0 - (2.0 * x2)) - x1) / t_1))); t_4 = x1 + (((x2 * x2) * x1) * 8.0); tmp = 0.0; if (t_3 <= -1e+229) tmp = t_4; elseif (t_3 <= 4e+177) tmp = x1 + (-6.0 * x2); elseif (t_3 <= Inf) tmp = t_4; else tmp = x1 + (((x1 * x1) * x2) * 8.0); end tmp_2 = tmp; end
code[x1_, x2_] := Block[{t$95$0 = N[(N[(3.0 * x1), $MachinePrecision] * x1), $MachinePrecision]}, Block[{t$95$1 = N[(N[(x1 * x1), $MachinePrecision] + 1.0), $MachinePrecision]}, Block[{t$95$2 = N[(N[(N[(t$95$0 + N[(2.0 * x2), $MachinePrecision]), $MachinePrecision] - x1), $MachinePrecision] / t$95$1), $MachinePrecision]}, Block[{t$95$3 = N[(x1 + N[(N[(N[(N[(N[(N[(N[(N[(N[(2.0 * x1), $MachinePrecision] * t$95$2), $MachinePrecision] * N[(t$95$2 - 3.0), $MachinePrecision]), $MachinePrecision] + N[(N[(x1 * x1), $MachinePrecision] * N[(N[(4.0 * t$95$2), $MachinePrecision] - 6.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * t$95$1), $MachinePrecision] + N[(t$95$0 * t$95$2), $MachinePrecision]), $MachinePrecision] + N[(N[(x1 * x1), $MachinePrecision] * x1), $MachinePrecision]), $MachinePrecision] + x1), $MachinePrecision] + N[(3.0 * N[(N[(N[(t$95$0 - N[(2.0 * x2), $MachinePrecision]), $MachinePrecision] - x1), $MachinePrecision] / t$95$1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$4 = N[(x1 + N[(N[(N[(x2 * x2), $MachinePrecision] * x1), $MachinePrecision] * 8.0), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$3, -1e+229], t$95$4, If[LessEqual[t$95$3, 4e+177], N[(x1 + N[(-6.0 * x2), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$3, Infinity], t$95$4, N[(x1 + N[(N[(N[(x1 * x1), $MachinePrecision] * x2), $MachinePrecision] * 8.0), $MachinePrecision]), $MachinePrecision]]]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \left(3 \cdot x1\right) \cdot x1\\
t_1 := x1 \cdot x1 + 1\\
t_2 := \frac{\left(t\_0 + 2 \cdot x2\right) - x1}{t\_1}\\
t_3 := x1 + \left(\left(\left(\left(\left(\left(\left(2 \cdot x1\right) \cdot t\_2\right) \cdot \left(t\_2 - 3\right) + \left(x1 \cdot x1\right) \cdot \left(4 \cdot t\_2 - 6\right)\right) \cdot t\_1 + t\_0 \cdot t\_2\right) + \left(x1 \cdot x1\right) \cdot x1\right) + x1\right) + 3 \cdot \frac{\left(t\_0 - 2 \cdot x2\right) - x1}{t\_1}\right)\\
t_4 := x1 + \left(\left(x2 \cdot x2\right) \cdot x1\right) \cdot 8\\
\mathbf{if}\;t\_3 \leq -1 \cdot 10^{+229}:\\
\;\;\;\;t\_4\\
\mathbf{elif}\;t\_3 \leq 4 \cdot 10^{+177}:\\
\;\;\;\;x1 + -6 \cdot x2\\
\mathbf{elif}\;t\_3 \leq \infty:\\
\;\;\;\;t\_4\\
\mathbf{else}:\\
\;\;\;\;x1 + \left(\left(x1 \cdot x1\right) \cdot x2\right) \cdot 8\\
\end{array}
\end{array}
if (+.f64 x1 (+.f64 (+.f64 (+.f64 (+.f64 (*.f64 (+.f64 (*.f64 (*.f64 (*.f64 #s(literal 2 binary64) x1) (/.f64 (-.f64 (+.f64 (*.f64 (*.f64 #s(literal 3 binary64) x1) x1) (*.f64 #s(literal 2 binary64) x2)) x1) (+.f64 (*.f64 x1 x1) #s(literal 1 binary64)))) (-.f64 (/.f64 (-.f64 (+.f64 (*.f64 (*.f64 #s(literal 3 binary64) x1) x1) (*.f64 #s(literal 2 binary64) x2)) x1) (+.f64 (*.f64 x1 x1) #s(literal 1 binary64))) #s(literal 3 binary64))) (*.f64 (*.f64 x1 x1) (-.f64 (*.f64 #s(literal 4 binary64) (/.f64 (-.f64 (+.f64 (*.f64 (*.f64 #s(literal 3 binary64) x1) x1) (*.f64 #s(literal 2 binary64) x2)) x1) (+.f64 (*.f64 x1 x1) #s(literal 1 binary64)))) #s(literal 6 binary64)))) (+.f64 (*.f64 x1 x1) #s(literal 1 binary64))) (*.f64 (*.f64 (*.f64 #s(literal 3 binary64) x1) x1) (/.f64 (-.f64 (+.f64 (*.f64 (*.f64 #s(literal 3 binary64) x1) x1) (*.f64 #s(literal 2 binary64) x2)) x1) (+.f64 (*.f64 x1 x1) #s(literal 1 binary64))))) (*.f64 (*.f64 x1 x1) x1)) x1) (*.f64 #s(literal 3 binary64) (/.f64 (-.f64 (-.f64 (*.f64 (*.f64 #s(literal 3 binary64) x1) x1) (*.f64 #s(literal 2 binary64) x2)) x1) (+.f64 (*.f64 x1 x1) #s(literal 1 binary64)))))) < -9.9999999999999999e228 or 4e177 < (+.f64 x1 (+.f64 (+.f64 (+.f64 (+.f64 (*.f64 (+.f64 (*.f64 (*.f64 (*.f64 #s(literal 2 binary64) x1) (/.f64 (-.f64 (+.f64 (*.f64 (*.f64 #s(literal 3 binary64) x1) x1) (*.f64 #s(literal 2 binary64) x2)) x1) (+.f64 (*.f64 x1 x1) #s(literal 1 binary64)))) (-.f64 (/.f64 (-.f64 (+.f64 (*.f64 (*.f64 #s(literal 3 binary64) x1) x1) (*.f64 #s(literal 2 binary64) x2)) x1) (+.f64 (*.f64 x1 x1) #s(literal 1 binary64))) #s(literal 3 binary64))) (*.f64 (*.f64 x1 x1) (-.f64 (*.f64 #s(literal 4 binary64) (/.f64 (-.f64 (+.f64 (*.f64 (*.f64 #s(literal 3 binary64) x1) x1) (*.f64 #s(literal 2 binary64) x2)) x1) (+.f64 (*.f64 x1 x1) #s(literal 1 binary64)))) #s(literal 6 binary64)))) (+.f64 (*.f64 x1 x1) #s(literal 1 binary64))) (*.f64 (*.f64 (*.f64 #s(literal 3 binary64) x1) x1) (/.f64 (-.f64 (+.f64 (*.f64 (*.f64 #s(literal 3 binary64) x1) x1) (*.f64 #s(literal 2 binary64) x2)) x1) (+.f64 (*.f64 x1 x1) #s(literal 1 binary64))))) (*.f64 (*.f64 x1 x1) x1)) x1) (*.f64 #s(literal 3 binary64) (/.f64 (-.f64 (-.f64 (*.f64 (*.f64 #s(literal 3 binary64) x1) x1) (*.f64 #s(literal 2 binary64) x2)) x1) (+.f64 (*.f64 x1 x1) #s(literal 1 binary64)))))) < +inf.0Initial program 99.9%
Taylor expanded in x1 around 0
Applied rewrites7.6%
Taylor expanded in x2 around 0
Applied rewrites50.9%
Taylor expanded in x2 around inf
Applied rewrites36.8%
Taylor expanded in x1 around 0
Applied rewrites57.1%
if -9.9999999999999999e228 < (+.f64 x1 (+.f64 (+.f64 (+.f64 (+.f64 (*.f64 (+.f64 (*.f64 (*.f64 (*.f64 #s(literal 2 binary64) x1) (/.f64 (-.f64 (+.f64 (*.f64 (*.f64 #s(literal 3 binary64) x1) x1) (*.f64 #s(literal 2 binary64) x2)) x1) (+.f64 (*.f64 x1 x1) #s(literal 1 binary64)))) (-.f64 (/.f64 (-.f64 (+.f64 (*.f64 (*.f64 #s(literal 3 binary64) x1) x1) (*.f64 #s(literal 2 binary64) x2)) x1) (+.f64 (*.f64 x1 x1) #s(literal 1 binary64))) #s(literal 3 binary64))) (*.f64 (*.f64 x1 x1) (-.f64 (*.f64 #s(literal 4 binary64) (/.f64 (-.f64 (+.f64 (*.f64 (*.f64 #s(literal 3 binary64) x1) x1) (*.f64 #s(literal 2 binary64) x2)) x1) (+.f64 (*.f64 x1 x1) #s(literal 1 binary64)))) #s(literal 6 binary64)))) (+.f64 (*.f64 x1 x1) #s(literal 1 binary64))) (*.f64 (*.f64 (*.f64 #s(literal 3 binary64) x1) x1) (/.f64 (-.f64 (+.f64 (*.f64 (*.f64 #s(literal 3 binary64) x1) x1) (*.f64 #s(literal 2 binary64) x2)) x1) (+.f64 (*.f64 x1 x1) #s(literal 1 binary64))))) (*.f64 (*.f64 x1 x1) x1)) x1) (*.f64 #s(literal 3 binary64) (/.f64 (-.f64 (-.f64 (*.f64 (*.f64 #s(literal 3 binary64) x1) x1) (*.f64 #s(literal 2 binary64) x2)) x1) (+.f64 (*.f64 x1 x1) #s(literal 1 binary64)))))) < 4e177Initial program 99.4%
Taylor expanded in x1 around 0
lower-*.f6458.2
Applied rewrites58.2%
if +inf.0 < (+.f64 x1 (+.f64 (+.f64 (+.f64 (+.f64 (*.f64 (+.f64 (*.f64 (*.f64 (*.f64 #s(literal 2 binary64) x1) (/.f64 (-.f64 (+.f64 (*.f64 (*.f64 #s(literal 3 binary64) x1) x1) (*.f64 #s(literal 2 binary64) x2)) x1) (+.f64 (*.f64 x1 x1) #s(literal 1 binary64)))) (-.f64 (/.f64 (-.f64 (+.f64 (*.f64 (*.f64 #s(literal 3 binary64) x1) x1) (*.f64 #s(literal 2 binary64) x2)) x1) (+.f64 (*.f64 x1 x1) #s(literal 1 binary64))) #s(literal 3 binary64))) (*.f64 (*.f64 x1 x1) (-.f64 (*.f64 #s(literal 4 binary64) (/.f64 (-.f64 (+.f64 (*.f64 (*.f64 #s(literal 3 binary64) x1) x1) (*.f64 #s(literal 2 binary64) x2)) x1) (+.f64 (*.f64 x1 x1) #s(literal 1 binary64)))) #s(literal 6 binary64)))) (+.f64 (*.f64 x1 x1) #s(literal 1 binary64))) (*.f64 (*.f64 (*.f64 #s(literal 3 binary64) x1) x1) (/.f64 (-.f64 (+.f64 (*.f64 (*.f64 #s(literal 3 binary64) x1) x1) (*.f64 #s(literal 2 binary64) x2)) x1) (+.f64 (*.f64 x1 x1) #s(literal 1 binary64))))) (*.f64 (*.f64 x1 x1) x1)) x1) (*.f64 #s(literal 3 binary64) (/.f64 (-.f64 (-.f64 (*.f64 (*.f64 #s(literal 3 binary64) x1) x1) (*.f64 #s(literal 2 binary64) x2)) x1) (+.f64 (*.f64 x1 x1) #s(literal 1 binary64)))))) Initial program 0.0%
Applied rewrites13.5%
Taylor expanded in x1 around -inf
*-commutativeN/A
lower-*.f64N/A
Applied rewrites100.0%
Taylor expanded in x2 around inf
Applied rewrites46.0%
(FPCore (x1 x2)
:precision binary64
(let* ((t_0 (+ (* x1 x1) 1.0))
(t_1 (* (* 3.0 x1) x1))
(t_2 (/ (- (+ t_1 (* 2.0 x2)) x1) t_0))
(t_3 (/ (- (fma x2 2.0 t_1) x1) (fma x1 x1 1.0)))
(t_4 (* 3.0 (/ (- (- t_1 (* 2.0 x2)) x1) t_0)))
(t_5
(+
x1
(+
(+
(+
(+
(*
(+
(* (* (* 2.0 x1) t_2) (- t_2 3.0))
(* (* x1 x1) (- (* 4.0 t_2) 6.0)))
t_0)
(* t_1 t_2))
(* (* x1 x1) x1))
x1)
t_4))))
(if (<= t_5 4e-7)
(+ x1 (+ (+ (* (* (* (/ x1 (fma x1 x1 1.0)) 8.0) x2) x2) x1) t_4))
(if (<= t_5 INFINITY)
(+
x1
(fma
(* x1 x1)
x1
(+
(fma
(fma
(fma 4.0 t_3 -6.0)
(* x1 x1)
(* (- t_3 3.0) (* t_3 (* x1 2.0))))
(fma x1 x1 1.0)
(* 3.0 t_1))
(fma (* -2.0 x2) 3.0 x1))))
(+ x1 (* (fma (fma 6.0 x1 -3.0) x1 -3.0) (* x1 x1)))))))
double code(double x1, double x2) {
double t_0 = (x1 * x1) + 1.0;
double t_1 = (3.0 * x1) * x1;
double t_2 = ((t_1 + (2.0 * x2)) - x1) / t_0;
double t_3 = (fma(x2, 2.0, t_1) - x1) / fma(x1, x1, 1.0);
double t_4 = 3.0 * (((t_1 - (2.0 * x2)) - x1) / t_0);
double t_5 = x1 + (((((((((2.0 * x1) * t_2) * (t_2 - 3.0)) + ((x1 * x1) * ((4.0 * t_2) - 6.0))) * t_0) + (t_1 * t_2)) + ((x1 * x1) * x1)) + x1) + t_4);
double tmp;
if (t_5 <= 4e-7) {
tmp = x1 + ((((((x1 / fma(x1, x1, 1.0)) * 8.0) * x2) * x2) + x1) + t_4);
} else if (t_5 <= ((double) INFINITY)) {
tmp = x1 + fma((x1 * x1), x1, (fma(fma(fma(4.0, t_3, -6.0), (x1 * x1), ((t_3 - 3.0) * (t_3 * (x1 * 2.0)))), fma(x1, x1, 1.0), (3.0 * t_1)) + fma((-2.0 * x2), 3.0, x1)));
} else {
tmp = x1 + (fma(fma(6.0, x1, -3.0), x1, -3.0) * (x1 * x1));
}
return tmp;
}
function code(x1, x2) t_0 = Float64(Float64(x1 * x1) + 1.0) t_1 = Float64(Float64(3.0 * x1) * x1) t_2 = Float64(Float64(Float64(t_1 + Float64(2.0 * x2)) - x1) / t_0) t_3 = Float64(Float64(fma(x2, 2.0, t_1) - x1) / fma(x1, x1, 1.0)) t_4 = Float64(3.0 * Float64(Float64(Float64(t_1 - Float64(2.0 * x2)) - x1) / t_0)) t_5 = Float64(x1 + Float64(Float64(Float64(Float64(Float64(Float64(Float64(Float64(Float64(2.0 * x1) * t_2) * Float64(t_2 - 3.0)) + Float64(Float64(x1 * x1) * Float64(Float64(4.0 * t_2) - 6.0))) * t_0) + Float64(t_1 * t_2)) + Float64(Float64(x1 * x1) * x1)) + x1) + t_4)) tmp = 0.0 if (t_5 <= 4e-7) tmp = Float64(x1 + Float64(Float64(Float64(Float64(Float64(Float64(x1 / fma(x1, x1, 1.0)) * 8.0) * x2) * x2) + x1) + t_4)); elseif (t_5 <= Inf) tmp = Float64(x1 + fma(Float64(x1 * x1), x1, Float64(fma(fma(fma(4.0, t_3, -6.0), Float64(x1 * x1), Float64(Float64(t_3 - 3.0) * Float64(t_3 * Float64(x1 * 2.0)))), fma(x1, x1, 1.0), Float64(3.0 * t_1)) + fma(Float64(-2.0 * x2), 3.0, x1)))); else tmp = Float64(x1 + Float64(fma(fma(6.0, x1, -3.0), x1, -3.0) * Float64(x1 * x1))); end return tmp end
code[x1_, x2_] := Block[{t$95$0 = N[(N[(x1 * x1), $MachinePrecision] + 1.0), $MachinePrecision]}, Block[{t$95$1 = N[(N[(3.0 * x1), $MachinePrecision] * x1), $MachinePrecision]}, Block[{t$95$2 = N[(N[(N[(t$95$1 + N[(2.0 * x2), $MachinePrecision]), $MachinePrecision] - x1), $MachinePrecision] / t$95$0), $MachinePrecision]}, Block[{t$95$3 = N[(N[(N[(x2 * 2.0 + t$95$1), $MachinePrecision] - x1), $MachinePrecision] / N[(x1 * x1 + 1.0), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$4 = N[(3.0 * N[(N[(N[(t$95$1 - N[(2.0 * x2), $MachinePrecision]), $MachinePrecision] - x1), $MachinePrecision] / t$95$0), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$5 = N[(x1 + N[(N[(N[(N[(N[(N[(N[(N[(N[(2.0 * x1), $MachinePrecision] * t$95$2), $MachinePrecision] * N[(t$95$2 - 3.0), $MachinePrecision]), $MachinePrecision] + N[(N[(x1 * x1), $MachinePrecision] * N[(N[(4.0 * t$95$2), $MachinePrecision] - 6.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * t$95$0), $MachinePrecision] + N[(t$95$1 * t$95$2), $MachinePrecision]), $MachinePrecision] + N[(N[(x1 * x1), $MachinePrecision] * x1), $MachinePrecision]), $MachinePrecision] + x1), $MachinePrecision] + t$95$4), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$5, 4e-7], N[(x1 + N[(N[(N[(N[(N[(N[(x1 / N[(x1 * x1 + 1.0), $MachinePrecision]), $MachinePrecision] * 8.0), $MachinePrecision] * x2), $MachinePrecision] * x2), $MachinePrecision] + x1), $MachinePrecision] + t$95$4), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$5, Infinity], N[(x1 + N[(N[(x1 * x1), $MachinePrecision] * x1 + N[(N[(N[(N[(4.0 * t$95$3 + -6.0), $MachinePrecision] * N[(x1 * x1), $MachinePrecision] + N[(N[(t$95$3 - 3.0), $MachinePrecision] * N[(t$95$3 * N[(x1 * 2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(x1 * x1 + 1.0), $MachinePrecision] + N[(3.0 * t$95$1), $MachinePrecision]), $MachinePrecision] + N[(N[(-2.0 * x2), $MachinePrecision] * 3.0 + x1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x1 + N[(N[(N[(6.0 * x1 + -3.0), $MachinePrecision] * x1 + -3.0), $MachinePrecision] * N[(x1 * x1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := x1 \cdot x1 + 1\\
t_1 := \left(3 \cdot x1\right) \cdot x1\\
t_2 := \frac{\left(t\_1 + 2 \cdot x2\right) - x1}{t\_0}\\
t_3 := \frac{\mathsf{fma}\left(x2, 2, t\_1\right) - x1}{\mathsf{fma}\left(x1, x1, 1\right)}\\
t_4 := 3 \cdot \frac{\left(t\_1 - 2 \cdot x2\right) - x1}{t\_0}\\
t_5 := x1 + \left(\left(\left(\left(\left(\left(\left(2 \cdot x1\right) \cdot t\_2\right) \cdot \left(t\_2 - 3\right) + \left(x1 \cdot x1\right) \cdot \left(4 \cdot t\_2 - 6\right)\right) \cdot t\_0 + t\_1 \cdot t\_2\right) + \left(x1 \cdot x1\right) \cdot x1\right) + x1\right) + t\_4\right)\\
\mathbf{if}\;t\_5 \leq 4 \cdot 10^{-7}:\\
\;\;\;\;x1 + \left(\left(\left(\left(\frac{x1}{\mathsf{fma}\left(x1, x1, 1\right)} \cdot 8\right) \cdot x2\right) \cdot x2 + x1\right) + t\_4\right)\\
\mathbf{elif}\;t\_5 \leq \infty:\\
\;\;\;\;x1 + \mathsf{fma}\left(x1 \cdot x1, x1, \mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(4, t\_3, -6\right), x1 \cdot x1, \left(t\_3 - 3\right) \cdot \left(t\_3 \cdot \left(x1 \cdot 2\right)\right)\right), \mathsf{fma}\left(x1, x1, 1\right), 3 \cdot t\_1\right) + \mathsf{fma}\left(-2 \cdot x2, 3, x1\right)\right)\\
\mathbf{else}:\\
\;\;\;\;x1 + \mathsf{fma}\left(\mathsf{fma}\left(6, x1, -3\right), x1, -3\right) \cdot \left(x1 \cdot x1\right)\\
\end{array}
\end{array}
if (+.f64 x1 (+.f64 (+.f64 (+.f64 (+.f64 (*.f64 (+.f64 (*.f64 (*.f64 (*.f64 #s(literal 2 binary64) x1) (/.f64 (-.f64 (+.f64 (*.f64 (*.f64 #s(literal 3 binary64) x1) x1) (*.f64 #s(literal 2 binary64) x2)) x1) (+.f64 (*.f64 x1 x1) #s(literal 1 binary64)))) (-.f64 (/.f64 (-.f64 (+.f64 (*.f64 (*.f64 #s(literal 3 binary64) x1) x1) (*.f64 #s(literal 2 binary64) x2)) x1) (+.f64 (*.f64 x1 x1) #s(literal 1 binary64))) #s(literal 3 binary64))) (*.f64 (*.f64 x1 x1) (-.f64 (*.f64 #s(literal 4 binary64) (/.f64 (-.f64 (+.f64 (*.f64 (*.f64 #s(literal 3 binary64) x1) x1) (*.f64 #s(literal 2 binary64) x2)) x1) (+.f64 (*.f64 x1 x1) #s(literal 1 binary64)))) #s(literal 6 binary64)))) (+.f64 (*.f64 x1 x1) #s(literal 1 binary64))) (*.f64 (*.f64 (*.f64 #s(literal 3 binary64) x1) x1) (/.f64 (-.f64 (+.f64 (*.f64 (*.f64 #s(literal 3 binary64) x1) x1) (*.f64 #s(literal 2 binary64) x2)) x1) (+.f64 (*.f64 x1 x1) #s(literal 1 binary64))))) (*.f64 (*.f64 x1 x1) x1)) x1) (*.f64 #s(literal 3 binary64) (/.f64 (-.f64 (-.f64 (*.f64 (*.f64 #s(literal 3 binary64) x1) x1) (*.f64 #s(literal 2 binary64) x2)) x1) (+.f64 (*.f64 x1 x1) #s(literal 1 binary64)))))) < 3.9999999999999998e-7Initial program 99.5%
Taylor expanded in x2 around inf
associate-*r/N/A
associate-*r*N/A
associate-*l/N/A
associate-*r/N/A
unpow2N/A
associate-*r*N/A
lower-*.f64N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-/.f64N/A
+-commutativeN/A
unpow2N/A
lower-fma.f6499.3
Applied rewrites99.3%
if 3.9999999999999998e-7 < (+.f64 x1 (+.f64 (+.f64 (+.f64 (+.f64 (*.f64 (+.f64 (*.f64 (*.f64 (*.f64 #s(literal 2 binary64) x1) (/.f64 (-.f64 (+.f64 (*.f64 (*.f64 #s(literal 3 binary64) x1) x1) (*.f64 #s(literal 2 binary64) x2)) x1) (+.f64 (*.f64 x1 x1) #s(literal 1 binary64)))) (-.f64 (/.f64 (-.f64 (+.f64 (*.f64 (*.f64 #s(literal 3 binary64) x1) x1) (*.f64 #s(literal 2 binary64) x2)) x1) (+.f64 (*.f64 x1 x1) #s(literal 1 binary64))) #s(literal 3 binary64))) (*.f64 (*.f64 x1 x1) (-.f64 (*.f64 #s(literal 4 binary64) (/.f64 (-.f64 (+.f64 (*.f64 (*.f64 #s(literal 3 binary64) x1) x1) (*.f64 #s(literal 2 binary64) x2)) x1) (+.f64 (*.f64 x1 x1) #s(literal 1 binary64)))) #s(literal 6 binary64)))) (+.f64 (*.f64 x1 x1) #s(literal 1 binary64))) (*.f64 (*.f64 (*.f64 #s(literal 3 binary64) x1) x1) (/.f64 (-.f64 (+.f64 (*.f64 (*.f64 #s(literal 3 binary64) x1) x1) (*.f64 #s(literal 2 binary64) x2)) x1) (+.f64 (*.f64 x1 x1) #s(literal 1 binary64))))) (*.f64 (*.f64 x1 x1) x1)) x1) (*.f64 #s(literal 3 binary64) (/.f64 (-.f64 (-.f64 (*.f64 (*.f64 #s(literal 3 binary64) x1) x1) (*.f64 #s(literal 2 binary64) x2)) x1) (+.f64 (*.f64 x1 x1) #s(literal 1 binary64)))))) < +inf.0Initial program 99.7%
Applied rewrites99.8%
Taylor expanded in x1 around inf
Applied rewrites99.1%
Taylor expanded in x1 around 0
lower-*.f6499.1
Applied rewrites99.1%
if +inf.0 < (+.f64 x1 (+.f64 (+.f64 (+.f64 (+.f64 (*.f64 (+.f64 (*.f64 (*.f64 (*.f64 #s(literal 2 binary64) x1) (/.f64 (-.f64 (+.f64 (*.f64 (*.f64 #s(literal 3 binary64) x1) x1) (*.f64 #s(literal 2 binary64) x2)) x1) (+.f64 (*.f64 x1 x1) #s(literal 1 binary64)))) (-.f64 (/.f64 (-.f64 (+.f64 (*.f64 (*.f64 #s(literal 3 binary64) x1) x1) (*.f64 #s(literal 2 binary64) x2)) x1) (+.f64 (*.f64 x1 x1) #s(literal 1 binary64))) #s(literal 3 binary64))) (*.f64 (*.f64 x1 x1) (-.f64 (*.f64 #s(literal 4 binary64) (/.f64 (-.f64 (+.f64 (*.f64 (*.f64 #s(literal 3 binary64) x1) x1) (*.f64 #s(literal 2 binary64) x2)) x1) (+.f64 (*.f64 x1 x1) #s(literal 1 binary64)))) #s(literal 6 binary64)))) (+.f64 (*.f64 x1 x1) #s(literal 1 binary64))) (*.f64 (*.f64 (*.f64 #s(literal 3 binary64) x1) x1) (/.f64 (-.f64 (+.f64 (*.f64 (*.f64 #s(literal 3 binary64) x1) x1) (*.f64 #s(literal 2 binary64) x2)) x1) (+.f64 (*.f64 x1 x1) #s(literal 1 binary64))))) (*.f64 (*.f64 x1 x1) x1)) x1) (*.f64 #s(literal 3 binary64) (/.f64 (-.f64 (-.f64 (*.f64 (*.f64 #s(literal 3 binary64) x1) x1) (*.f64 #s(literal 2 binary64) x2)) x1) (+.f64 (*.f64 x1 x1) #s(literal 1 binary64)))))) Initial program 0.0%
Applied rewrites13.5%
Taylor expanded in x1 around -inf
*-commutativeN/A
lower-*.f64N/A
Applied rewrites100.0%
Taylor expanded in x1 around 0
Applied rewrites100.0%
Taylor expanded in x2 around 0
Applied rewrites100.0%
(FPCore (x1 x2)
:precision binary64
(let* ((t_0 (* (* 3.0 x1) x1))
(t_1 (+ (* x1 x1) 1.0))
(t_2 (/ (- (+ t_0 (* 2.0 x2)) x1) t_1))
(t_3
(+
x1
(+
(+
(+
(+
(*
(+
(* (* (* 2.0 x1) t_2) (- t_2 3.0))
(* (* x1 x1) (- (* 4.0 t_2) 6.0)))
t_1)
(* t_0 t_2))
(* (* x1 x1) x1))
x1)
(* 3.0 (/ (- (- t_0 (* 2.0 x2)) x1) t_1))))))
(if (<= t_3 -1e+229)
(+ x1 (* (* (* x2 x2) x1) 8.0))
(if (<= t_3 4e+177)
(+ x1 (fma (fma (fma -19.0 x1 9.0) x1 -2.0) x1 (* -6.0 x2)))
(+ x1 (* (fma (fma 6.0 x1 -3.0) x1 -3.0) (* x1 x1)))))))
double code(double x1, double x2) {
double t_0 = (3.0 * x1) * x1;
double t_1 = (x1 * x1) + 1.0;
double t_2 = ((t_0 + (2.0 * x2)) - x1) / t_1;
double t_3 = x1 + (((((((((2.0 * x1) * t_2) * (t_2 - 3.0)) + ((x1 * x1) * ((4.0 * t_2) - 6.0))) * t_1) + (t_0 * t_2)) + ((x1 * x1) * x1)) + x1) + (3.0 * (((t_0 - (2.0 * x2)) - x1) / t_1)));
double tmp;
if (t_3 <= -1e+229) {
tmp = x1 + (((x2 * x2) * x1) * 8.0);
} else if (t_3 <= 4e+177) {
tmp = x1 + fma(fma(fma(-19.0, x1, 9.0), x1, -2.0), x1, (-6.0 * x2));
} else {
tmp = x1 + (fma(fma(6.0, x1, -3.0), x1, -3.0) * (x1 * x1));
}
return tmp;
}
function code(x1, x2) t_0 = Float64(Float64(3.0 * x1) * x1) t_1 = Float64(Float64(x1 * x1) + 1.0) t_2 = Float64(Float64(Float64(t_0 + Float64(2.0 * x2)) - x1) / t_1) t_3 = Float64(x1 + Float64(Float64(Float64(Float64(Float64(Float64(Float64(Float64(Float64(2.0 * x1) * t_2) * Float64(t_2 - 3.0)) + Float64(Float64(x1 * x1) * Float64(Float64(4.0 * t_2) - 6.0))) * t_1) + Float64(t_0 * t_2)) + Float64(Float64(x1 * x1) * x1)) + x1) + Float64(3.0 * Float64(Float64(Float64(t_0 - Float64(2.0 * x2)) - x1) / t_1)))) tmp = 0.0 if (t_3 <= -1e+229) tmp = Float64(x1 + Float64(Float64(Float64(x2 * x2) * x1) * 8.0)); elseif (t_3 <= 4e+177) tmp = Float64(x1 + fma(fma(fma(-19.0, x1, 9.0), x1, -2.0), x1, Float64(-6.0 * x2))); else tmp = Float64(x1 + Float64(fma(fma(6.0, x1, -3.0), x1, -3.0) * Float64(x1 * x1))); end return tmp end
code[x1_, x2_] := Block[{t$95$0 = N[(N[(3.0 * x1), $MachinePrecision] * x1), $MachinePrecision]}, Block[{t$95$1 = N[(N[(x1 * x1), $MachinePrecision] + 1.0), $MachinePrecision]}, Block[{t$95$2 = N[(N[(N[(t$95$0 + N[(2.0 * x2), $MachinePrecision]), $MachinePrecision] - x1), $MachinePrecision] / t$95$1), $MachinePrecision]}, Block[{t$95$3 = N[(x1 + N[(N[(N[(N[(N[(N[(N[(N[(N[(2.0 * x1), $MachinePrecision] * t$95$2), $MachinePrecision] * N[(t$95$2 - 3.0), $MachinePrecision]), $MachinePrecision] + N[(N[(x1 * x1), $MachinePrecision] * N[(N[(4.0 * t$95$2), $MachinePrecision] - 6.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * t$95$1), $MachinePrecision] + N[(t$95$0 * t$95$2), $MachinePrecision]), $MachinePrecision] + N[(N[(x1 * x1), $MachinePrecision] * x1), $MachinePrecision]), $MachinePrecision] + x1), $MachinePrecision] + N[(3.0 * N[(N[(N[(t$95$0 - N[(2.0 * x2), $MachinePrecision]), $MachinePrecision] - x1), $MachinePrecision] / t$95$1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$3, -1e+229], N[(x1 + N[(N[(N[(x2 * x2), $MachinePrecision] * x1), $MachinePrecision] * 8.0), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$3, 4e+177], N[(x1 + N[(N[(N[(-19.0 * x1 + 9.0), $MachinePrecision] * x1 + -2.0), $MachinePrecision] * x1 + N[(-6.0 * x2), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x1 + N[(N[(N[(6.0 * x1 + -3.0), $MachinePrecision] * x1 + -3.0), $MachinePrecision] * N[(x1 * x1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \left(3 \cdot x1\right) \cdot x1\\
t_1 := x1 \cdot x1 + 1\\
t_2 := \frac{\left(t\_0 + 2 \cdot x2\right) - x1}{t\_1}\\
t_3 := x1 + \left(\left(\left(\left(\left(\left(\left(2 \cdot x1\right) \cdot t\_2\right) \cdot \left(t\_2 - 3\right) + \left(x1 \cdot x1\right) \cdot \left(4 \cdot t\_2 - 6\right)\right) \cdot t\_1 + t\_0 \cdot t\_2\right) + \left(x1 \cdot x1\right) \cdot x1\right) + x1\right) + 3 \cdot \frac{\left(t\_0 - 2 \cdot x2\right) - x1}{t\_1}\right)\\
\mathbf{if}\;t\_3 \leq -1 \cdot 10^{+229}:\\
\;\;\;\;x1 + \left(\left(x2 \cdot x2\right) \cdot x1\right) \cdot 8\\
\mathbf{elif}\;t\_3 \leq 4 \cdot 10^{+177}:\\
\;\;\;\;x1 + \mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(-19, x1, 9\right), x1, -2\right), x1, -6 \cdot x2\right)\\
\mathbf{else}:\\
\;\;\;\;x1 + \mathsf{fma}\left(\mathsf{fma}\left(6, x1, -3\right), x1, -3\right) \cdot \left(x1 \cdot x1\right)\\
\end{array}
\end{array}
if (+.f64 x1 (+.f64 (+.f64 (+.f64 (+.f64 (*.f64 (+.f64 (*.f64 (*.f64 (*.f64 #s(literal 2 binary64) x1) (/.f64 (-.f64 (+.f64 (*.f64 (*.f64 #s(literal 3 binary64) x1) x1) (*.f64 #s(literal 2 binary64) x2)) x1) (+.f64 (*.f64 x1 x1) #s(literal 1 binary64)))) (-.f64 (/.f64 (-.f64 (+.f64 (*.f64 (*.f64 #s(literal 3 binary64) x1) x1) (*.f64 #s(literal 2 binary64) x2)) x1) (+.f64 (*.f64 x1 x1) #s(literal 1 binary64))) #s(literal 3 binary64))) (*.f64 (*.f64 x1 x1) (-.f64 (*.f64 #s(literal 4 binary64) (/.f64 (-.f64 (+.f64 (*.f64 (*.f64 #s(literal 3 binary64) x1) x1) (*.f64 #s(literal 2 binary64) x2)) x1) (+.f64 (*.f64 x1 x1) #s(literal 1 binary64)))) #s(literal 6 binary64)))) (+.f64 (*.f64 x1 x1) #s(literal 1 binary64))) (*.f64 (*.f64 (*.f64 #s(literal 3 binary64) x1) x1) (/.f64 (-.f64 (+.f64 (*.f64 (*.f64 #s(literal 3 binary64) x1) x1) (*.f64 #s(literal 2 binary64) x2)) x1) (+.f64 (*.f64 x1 x1) #s(literal 1 binary64))))) (*.f64 (*.f64 x1 x1) x1)) x1) (*.f64 #s(literal 3 binary64) (/.f64 (-.f64 (-.f64 (*.f64 (*.f64 #s(literal 3 binary64) x1) x1) (*.f64 #s(literal 2 binary64) x2)) x1) (+.f64 (*.f64 x1 x1) #s(literal 1 binary64)))))) < -9.9999999999999999e228Initial program 99.9%
Taylor expanded in x1 around 0
Applied rewrites4.2%
Taylor expanded in x2 around 0
Applied rewrites73.4%
Taylor expanded in x2 around inf
Applied rewrites55.9%
Taylor expanded in x1 around 0
Applied rewrites77.5%
if -9.9999999999999999e228 < (+.f64 x1 (+.f64 (+.f64 (+.f64 (+.f64 (*.f64 (+.f64 (*.f64 (*.f64 (*.f64 #s(literal 2 binary64) x1) (/.f64 (-.f64 (+.f64 (*.f64 (*.f64 #s(literal 3 binary64) x1) x1) (*.f64 #s(literal 2 binary64) x2)) x1) (+.f64 (*.f64 x1 x1) #s(literal 1 binary64)))) (-.f64 (/.f64 (-.f64 (+.f64 (*.f64 (*.f64 #s(literal 3 binary64) x1) x1) (*.f64 #s(literal 2 binary64) x2)) x1) (+.f64 (*.f64 x1 x1) #s(literal 1 binary64))) #s(literal 3 binary64))) (*.f64 (*.f64 x1 x1) (-.f64 (*.f64 #s(literal 4 binary64) (/.f64 (-.f64 (+.f64 (*.f64 (*.f64 #s(literal 3 binary64) x1) x1) (*.f64 #s(literal 2 binary64) x2)) x1) (+.f64 (*.f64 x1 x1) #s(literal 1 binary64)))) #s(literal 6 binary64)))) (+.f64 (*.f64 x1 x1) #s(literal 1 binary64))) (*.f64 (*.f64 (*.f64 #s(literal 3 binary64) x1) x1) (/.f64 (-.f64 (+.f64 (*.f64 (*.f64 #s(literal 3 binary64) x1) x1) (*.f64 #s(literal 2 binary64) x2)) x1) (+.f64 (*.f64 x1 x1) #s(literal 1 binary64))))) (*.f64 (*.f64 x1 x1) x1)) x1) (*.f64 #s(literal 3 binary64) (/.f64 (-.f64 (-.f64 (*.f64 (*.f64 #s(literal 3 binary64) x1) x1) (*.f64 #s(literal 2 binary64) x2)) x1) (+.f64 (*.f64 x1 x1) #s(literal 1 binary64)))))) < 4e177Initial program 99.4%
Taylor expanded in x1 around 0
Applied rewrites78.5%
Taylor expanded in x2 around 0
Applied rewrites79.4%
if 4e177 < (+.f64 x1 (+.f64 (+.f64 (+.f64 (+.f64 (*.f64 (+.f64 (*.f64 (*.f64 (*.f64 #s(literal 2 binary64) x1) (/.f64 (-.f64 (+.f64 (*.f64 (*.f64 #s(literal 3 binary64) x1) x1) (*.f64 #s(literal 2 binary64) x2)) x1) (+.f64 (*.f64 x1 x1) #s(literal 1 binary64)))) (-.f64 (/.f64 (-.f64 (+.f64 (*.f64 (*.f64 #s(literal 3 binary64) x1) x1) (*.f64 #s(literal 2 binary64) x2)) x1) (+.f64 (*.f64 x1 x1) #s(literal 1 binary64))) #s(literal 3 binary64))) (*.f64 (*.f64 x1 x1) (-.f64 (*.f64 #s(literal 4 binary64) (/.f64 (-.f64 (+.f64 (*.f64 (*.f64 #s(literal 3 binary64) x1) x1) (*.f64 #s(literal 2 binary64) x2)) x1) (+.f64 (*.f64 x1 x1) #s(literal 1 binary64)))) #s(literal 6 binary64)))) (+.f64 (*.f64 x1 x1) #s(literal 1 binary64))) (*.f64 (*.f64 (*.f64 #s(literal 3 binary64) x1) x1) (/.f64 (-.f64 (+.f64 (*.f64 (*.f64 #s(literal 3 binary64) x1) x1) (*.f64 #s(literal 2 binary64) x2)) x1) (+.f64 (*.f64 x1 x1) #s(literal 1 binary64))))) (*.f64 (*.f64 x1 x1) x1)) x1) (*.f64 #s(literal 3 binary64) (/.f64 (-.f64 (-.f64 (*.f64 (*.f64 #s(literal 3 binary64) x1) x1) (*.f64 #s(literal 2 binary64) x2)) x1) (+.f64 (*.f64 x1 x1) #s(literal 1 binary64)))))) Initial program 38.8%
Applied rewrites47.0%
Taylor expanded in x1 around -inf
*-commutativeN/A
lower-*.f64N/A
Applied rewrites87.6%
Taylor expanded in x1 around 0
Applied rewrites87.7%
Taylor expanded in x2 around 0
Applied rewrites87.7%
(FPCore (x1 x2)
:precision binary64
(let* ((t_0 (* (* 3.0 x1) x1))
(t_1 (+ (* x1 x1) 1.0))
(t_2 (/ (- (+ t_0 (* 2.0 x2)) x1) t_1)))
(if (<=
(+
x1
(+
(+
(+
(+
(*
(+
(* (* (* 2.0 x1) t_2) (- t_2 3.0))
(* (* x1 x1) (- (* 4.0 t_2) 6.0)))
t_1)
(* t_0 t_2))
(* (* x1 x1) x1))
x1)
(* 3.0 (/ (- (- t_0 (* 2.0 x2)) x1) t_1))))
5e+217)
(+ x1 (* -6.0 x2))
(+ x1 (* (* (* x1 x1) x2) 8.0)))))
double code(double x1, double x2) {
double t_0 = (3.0 * x1) * x1;
double t_1 = (x1 * x1) + 1.0;
double t_2 = ((t_0 + (2.0 * x2)) - x1) / t_1;
double tmp;
if ((x1 + (((((((((2.0 * x1) * t_2) * (t_2 - 3.0)) + ((x1 * x1) * ((4.0 * t_2) - 6.0))) * t_1) + (t_0 * t_2)) + ((x1 * x1) * x1)) + x1) + (3.0 * (((t_0 - (2.0 * x2)) - x1) / t_1)))) <= 5e+217) {
tmp = x1 + (-6.0 * x2);
} else {
tmp = x1 + (((x1 * x1) * x2) * 8.0);
}
return tmp;
}
real(8) function code(x1, x2)
real(8), intent (in) :: x1
real(8), intent (in) :: x2
real(8) :: t_0
real(8) :: t_1
real(8) :: t_2
real(8) :: tmp
t_0 = (3.0d0 * x1) * x1
t_1 = (x1 * x1) + 1.0d0
t_2 = ((t_0 + (2.0d0 * x2)) - x1) / t_1
if ((x1 + (((((((((2.0d0 * x1) * t_2) * (t_2 - 3.0d0)) + ((x1 * x1) * ((4.0d0 * t_2) - 6.0d0))) * t_1) + (t_0 * t_2)) + ((x1 * x1) * x1)) + x1) + (3.0d0 * (((t_0 - (2.0d0 * x2)) - x1) / t_1)))) <= 5d+217) then
tmp = x1 + ((-6.0d0) * x2)
else
tmp = x1 + (((x1 * x1) * x2) * 8.0d0)
end if
code = tmp
end function
public static double code(double x1, double x2) {
double t_0 = (3.0 * x1) * x1;
double t_1 = (x1 * x1) + 1.0;
double t_2 = ((t_0 + (2.0 * x2)) - x1) / t_1;
double tmp;
if ((x1 + (((((((((2.0 * x1) * t_2) * (t_2 - 3.0)) + ((x1 * x1) * ((4.0 * t_2) - 6.0))) * t_1) + (t_0 * t_2)) + ((x1 * x1) * x1)) + x1) + (3.0 * (((t_0 - (2.0 * x2)) - x1) / t_1)))) <= 5e+217) {
tmp = x1 + (-6.0 * x2);
} else {
tmp = x1 + (((x1 * x1) * x2) * 8.0);
}
return tmp;
}
def code(x1, x2): t_0 = (3.0 * x1) * x1 t_1 = (x1 * x1) + 1.0 t_2 = ((t_0 + (2.0 * x2)) - x1) / t_1 tmp = 0 if (x1 + (((((((((2.0 * x1) * t_2) * (t_2 - 3.0)) + ((x1 * x1) * ((4.0 * t_2) - 6.0))) * t_1) + (t_0 * t_2)) + ((x1 * x1) * x1)) + x1) + (3.0 * (((t_0 - (2.0 * x2)) - x1) / t_1)))) <= 5e+217: tmp = x1 + (-6.0 * x2) else: tmp = x1 + (((x1 * x1) * x2) * 8.0) return tmp
function code(x1, x2) t_0 = Float64(Float64(3.0 * x1) * x1) t_1 = Float64(Float64(x1 * x1) + 1.0) t_2 = Float64(Float64(Float64(t_0 + Float64(2.0 * x2)) - x1) / t_1) tmp = 0.0 if (Float64(x1 + Float64(Float64(Float64(Float64(Float64(Float64(Float64(Float64(Float64(2.0 * x1) * t_2) * Float64(t_2 - 3.0)) + Float64(Float64(x1 * x1) * Float64(Float64(4.0 * t_2) - 6.0))) * t_1) + Float64(t_0 * t_2)) + Float64(Float64(x1 * x1) * x1)) + x1) + Float64(3.0 * Float64(Float64(Float64(t_0 - Float64(2.0 * x2)) - x1) / t_1)))) <= 5e+217) tmp = Float64(x1 + Float64(-6.0 * x2)); else tmp = Float64(x1 + Float64(Float64(Float64(x1 * x1) * x2) * 8.0)); end return tmp end
function tmp_2 = code(x1, x2) t_0 = (3.0 * x1) * x1; t_1 = (x1 * x1) + 1.0; t_2 = ((t_0 + (2.0 * x2)) - x1) / t_1; tmp = 0.0; if ((x1 + (((((((((2.0 * x1) * t_2) * (t_2 - 3.0)) + ((x1 * x1) * ((4.0 * t_2) - 6.0))) * t_1) + (t_0 * t_2)) + ((x1 * x1) * x1)) + x1) + (3.0 * (((t_0 - (2.0 * x2)) - x1) / t_1)))) <= 5e+217) tmp = x1 + (-6.0 * x2); else tmp = x1 + (((x1 * x1) * x2) * 8.0); end tmp_2 = tmp; end
code[x1_, x2_] := Block[{t$95$0 = N[(N[(3.0 * x1), $MachinePrecision] * x1), $MachinePrecision]}, Block[{t$95$1 = N[(N[(x1 * x1), $MachinePrecision] + 1.0), $MachinePrecision]}, Block[{t$95$2 = N[(N[(N[(t$95$0 + N[(2.0 * x2), $MachinePrecision]), $MachinePrecision] - x1), $MachinePrecision] / t$95$1), $MachinePrecision]}, If[LessEqual[N[(x1 + N[(N[(N[(N[(N[(N[(N[(N[(N[(2.0 * x1), $MachinePrecision] * t$95$2), $MachinePrecision] * N[(t$95$2 - 3.0), $MachinePrecision]), $MachinePrecision] + N[(N[(x1 * x1), $MachinePrecision] * N[(N[(4.0 * t$95$2), $MachinePrecision] - 6.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * t$95$1), $MachinePrecision] + N[(t$95$0 * t$95$2), $MachinePrecision]), $MachinePrecision] + N[(N[(x1 * x1), $MachinePrecision] * x1), $MachinePrecision]), $MachinePrecision] + x1), $MachinePrecision] + N[(3.0 * N[(N[(N[(t$95$0 - N[(2.0 * x2), $MachinePrecision]), $MachinePrecision] - x1), $MachinePrecision] / t$95$1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 5e+217], N[(x1 + N[(-6.0 * x2), $MachinePrecision]), $MachinePrecision], N[(x1 + N[(N[(N[(x1 * x1), $MachinePrecision] * x2), $MachinePrecision] * 8.0), $MachinePrecision]), $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \left(3 \cdot x1\right) \cdot x1\\
t_1 := x1 \cdot x1 + 1\\
t_2 := \frac{\left(t\_0 + 2 \cdot x2\right) - x1}{t\_1}\\
\mathbf{if}\;x1 + \left(\left(\left(\left(\left(\left(\left(2 \cdot x1\right) \cdot t\_2\right) \cdot \left(t\_2 - 3\right) + \left(x1 \cdot x1\right) \cdot \left(4 \cdot t\_2 - 6\right)\right) \cdot t\_1 + t\_0 \cdot t\_2\right) + \left(x1 \cdot x1\right) \cdot x1\right) + x1\right) + 3 \cdot \frac{\left(t\_0 - 2 \cdot x2\right) - x1}{t\_1}\right) \leq 5 \cdot 10^{+217}:\\
\;\;\;\;x1 + -6 \cdot x2\\
\mathbf{else}:\\
\;\;\;\;x1 + \left(\left(x1 \cdot x1\right) \cdot x2\right) \cdot 8\\
\end{array}
\end{array}
if (+.f64 x1 (+.f64 (+.f64 (+.f64 (+.f64 (*.f64 (+.f64 (*.f64 (*.f64 (*.f64 #s(literal 2 binary64) x1) (/.f64 (-.f64 (+.f64 (*.f64 (*.f64 #s(literal 3 binary64) x1) x1) (*.f64 #s(literal 2 binary64) x2)) x1) (+.f64 (*.f64 x1 x1) #s(literal 1 binary64)))) (-.f64 (/.f64 (-.f64 (+.f64 (*.f64 (*.f64 #s(literal 3 binary64) x1) x1) (*.f64 #s(literal 2 binary64) x2)) x1) (+.f64 (*.f64 x1 x1) #s(literal 1 binary64))) #s(literal 3 binary64))) (*.f64 (*.f64 x1 x1) (-.f64 (*.f64 #s(literal 4 binary64) (/.f64 (-.f64 (+.f64 (*.f64 (*.f64 #s(literal 3 binary64) x1) x1) (*.f64 #s(literal 2 binary64) x2)) x1) (+.f64 (*.f64 x1 x1) #s(literal 1 binary64)))) #s(literal 6 binary64)))) (+.f64 (*.f64 x1 x1) #s(literal 1 binary64))) (*.f64 (*.f64 (*.f64 #s(literal 3 binary64) x1) x1) (/.f64 (-.f64 (+.f64 (*.f64 (*.f64 #s(literal 3 binary64) x1) x1) (*.f64 #s(literal 2 binary64) x2)) x1) (+.f64 (*.f64 x1 x1) #s(literal 1 binary64))))) (*.f64 (*.f64 x1 x1) x1)) x1) (*.f64 #s(literal 3 binary64) (/.f64 (-.f64 (-.f64 (*.f64 (*.f64 #s(literal 3 binary64) x1) x1) (*.f64 #s(literal 2 binary64) x2)) x1) (+.f64 (*.f64 x1 x1) #s(literal 1 binary64)))))) < 5.00000000000000041e217Initial program 99.5%
Taylor expanded in x1 around 0
lower-*.f6447.2
Applied rewrites47.2%
if 5.00000000000000041e217 < (+.f64 x1 (+.f64 (+.f64 (+.f64 (+.f64 (*.f64 (+.f64 (*.f64 (*.f64 (*.f64 #s(literal 2 binary64) x1) (/.f64 (-.f64 (+.f64 (*.f64 (*.f64 #s(literal 3 binary64) x1) x1) (*.f64 #s(literal 2 binary64) x2)) x1) (+.f64 (*.f64 x1 x1) #s(literal 1 binary64)))) (-.f64 (/.f64 (-.f64 (+.f64 (*.f64 (*.f64 #s(literal 3 binary64) x1) x1) (*.f64 #s(literal 2 binary64) x2)) x1) (+.f64 (*.f64 x1 x1) #s(literal 1 binary64))) #s(literal 3 binary64))) (*.f64 (*.f64 x1 x1) (-.f64 (*.f64 #s(literal 4 binary64) (/.f64 (-.f64 (+.f64 (*.f64 (*.f64 #s(literal 3 binary64) x1) x1) (*.f64 #s(literal 2 binary64) x2)) x1) (+.f64 (*.f64 x1 x1) #s(literal 1 binary64)))) #s(literal 6 binary64)))) (+.f64 (*.f64 x1 x1) #s(literal 1 binary64))) (*.f64 (*.f64 (*.f64 #s(literal 3 binary64) x1) x1) (/.f64 (-.f64 (+.f64 (*.f64 (*.f64 #s(literal 3 binary64) x1) x1) (*.f64 #s(literal 2 binary64) x2)) x1) (+.f64 (*.f64 x1 x1) #s(literal 1 binary64))))) (*.f64 (*.f64 x1 x1) x1)) x1) (*.f64 #s(literal 3 binary64) (/.f64 (-.f64 (-.f64 (*.f64 (*.f64 #s(literal 3 binary64) x1) x1) (*.f64 #s(literal 2 binary64) x2)) x1) (+.f64 (*.f64 x1 x1) #s(literal 1 binary64)))))) Initial program 35.6%
Applied rewrites44.3%
Taylor expanded in x1 around -inf
*-commutativeN/A
lower-*.f64N/A
Applied rewrites88.7%
Taylor expanded in x2 around inf
Applied rewrites33.9%
(FPCore (x1 x2)
:precision binary64
(if (<= x1 -3e+34)
(+ x1 (* (fma (fma 6.0 x1 -3.0) x1 (* 8.0 x2)) (* x1 x1)))
(if (<= x1 4100000.0)
(+
x1
(+
(+ (* (* (* (/ x1 (fma x1 x1 1.0)) 8.0) x2) x2) x1)
(* 3.0 (/ (- (- (* (* 3.0 x1) x1) (* 2.0 x2)) x1) (+ (* x1 x1) 1.0)))))
(+
(* (* (fma (fma 6.0 x1 -3.0) x1 (fma (fma 2.0 x2 -3.0) 4.0 9.0)) x1) x1)
x1))))
double code(double x1, double x2) {
double tmp;
if (x1 <= -3e+34) {
tmp = x1 + (fma(fma(6.0, x1, -3.0), x1, (8.0 * x2)) * (x1 * x1));
} else if (x1 <= 4100000.0) {
tmp = x1 + ((((((x1 / fma(x1, x1, 1.0)) * 8.0) * x2) * x2) + x1) + (3.0 * (((((3.0 * x1) * x1) - (2.0 * x2)) - x1) / ((x1 * x1) + 1.0))));
} else {
tmp = ((fma(fma(6.0, x1, -3.0), x1, fma(fma(2.0, x2, -3.0), 4.0, 9.0)) * x1) * x1) + x1;
}
return tmp;
}
function code(x1, x2) tmp = 0.0 if (x1 <= -3e+34) tmp = Float64(x1 + Float64(fma(fma(6.0, x1, -3.0), x1, Float64(8.0 * x2)) * Float64(x1 * x1))); elseif (x1 <= 4100000.0) tmp = Float64(x1 + Float64(Float64(Float64(Float64(Float64(Float64(x1 / fma(x1, x1, 1.0)) * 8.0) * x2) * x2) + x1) + Float64(3.0 * Float64(Float64(Float64(Float64(Float64(3.0 * x1) * x1) - Float64(2.0 * x2)) - x1) / Float64(Float64(x1 * x1) + 1.0))))); else tmp = Float64(Float64(Float64(fma(fma(6.0, x1, -3.0), x1, fma(fma(2.0, x2, -3.0), 4.0, 9.0)) * x1) * x1) + x1); end return tmp end
code[x1_, x2_] := If[LessEqual[x1, -3e+34], N[(x1 + N[(N[(N[(6.0 * x1 + -3.0), $MachinePrecision] * x1 + N[(8.0 * x2), $MachinePrecision]), $MachinePrecision] * N[(x1 * x1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[x1, 4100000.0], N[(x1 + N[(N[(N[(N[(N[(N[(x1 / N[(x1 * x1 + 1.0), $MachinePrecision]), $MachinePrecision] * 8.0), $MachinePrecision] * x2), $MachinePrecision] * x2), $MachinePrecision] + x1), $MachinePrecision] + N[(3.0 * N[(N[(N[(N[(N[(3.0 * x1), $MachinePrecision] * x1), $MachinePrecision] - N[(2.0 * x2), $MachinePrecision]), $MachinePrecision] - x1), $MachinePrecision] / N[(N[(x1 * x1), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(N[(N[(N[(6.0 * x1 + -3.0), $MachinePrecision] * x1 + N[(N[(2.0 * x2 + -3.0), $MachinePrecision] * 4.0 + 9.0), $MachinePrecision]), $MachinePrecision] * x1), $MachinePrecision] * x1), $MachinePrecision] + x1), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x1 \leq -3 \cdot 10^{+34}:\\
\;\;\;\;x1 + \mathsf{fma}\left(\mathsf{fma}\left(6, x1, -3\right), x1, 8 \cdot x2\right) \cdot \left(x1 \cdot x1\right)\\
\mathbf{elif}\;x1 \leq 4100000:\\
\;\;\;\;x1 + \left(\left(\left(\left(\frac{x1}{\mathsf{fma}\left(x1, x1, 1\right)} \cdot 8\right) \cdot x2\right) \cdot x2 + x1\right) + 3 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 - 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right)\\
\mathbf{else}:\\
\;\;\;\;\left(\mathsf{fma}\left(\mathsf{fma}\left(6, x1, -3\right), x1, \mathsf{fma}\left(\mathsf{fma}\left(2, x2, -3\right), 4, 9\right)\right) \cdot x1\right) \cdot x1 + x1\\
\end{array}
\end{array}
if x1 < -3.00000000000000018e34Initial program 19.6%
Applied rewrites37.4%
Taylor expanded in x1 around -inf
*-commutativeN/A
lower-*.f64N/A
Applied rewrites96.7%
Taylor expanded in x1 around 0
Applied rewrites96.7%
Taylor expanded in x2 around inf
Applied rewrites96.7%
if -3.00000000000000018e34 < x1 < 4.1e6Initial program 99.5%
Taylor expanded in x2 around inf
associate-*r/N/A
associate-*r*N/A
associate-*l/N/A
associate-*r/N/A
unpow2N/A
associate-*r*N/A
lower-*.f64N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-/.f64N/A
+-commutativeN/A
unpow2N/A
lower-fma.f6497.3
Applied rewrites97.3%
if 4.1e6 < x1 Initial program 53.9%
Applied rewrites53.9%
Taylor expanded in x1 around -inf
*-commutativeN/A
lower-*.f64N/A
Applied rewrites98.2%
Taylor expanded in x1 around 0
Applied rewrites98.3%
lift-+.f64N/A
Applied rewrites98.3%
(FPCore (x1 x2)
:precision binary64
(if (or (<= x1 -1.0) (not (<= x1 1.36)))
(+
(* (* (fma (fma 6.0 x1 -3.0) x1 (fma (fma 2.0 x2 -3.0) 4.0 9.0)) x1) x1)
x1)
(+
x1
(fma
(fma
x1
(fma (fma -8.0 (* x1 x1) 8.0) x2 (fma (fma 24.0 x1 12.0) x1 -12.0))
-6.0)
x2
(* (fma (fma -19.0 x1 9.0) x1 -2.0) x1)))))
double code(double x1, double x2) {
double tmp;
if ((x1 <= -1.0) || !(x1 <= 1.36)) {
tmp = ((fma(fma(6.0, x1, -3.0), x1, fma(fma(2.0, x2, -3.0), 4.0, 9.0)) * x1) * x1) + x1;
} else {
tmp = x1 + fma(fma(x1, fma(fma(-8.0, (x1 * x1), 8.0), x2, fma(fma(24.0, x1, 12.0), x1, -12.0)), -6.0), x2, (fma(fma(-19.0, x1, 9.0), x1, -2.0) * x1));
}
return tmp;
}
function code(x1, x2) tmp = 0.0 if ((x1 <= -1.0) || !(x1 <= 1.36)) tmp = Float64(Float64(Float64(fma(fma(6.0, x1, -3.0), x1, fma(fma(2.0, x2, -3.0), 4.0, 9.0)) * x1) * x1) + x1); else tmp = Float64(x1 + fma(fma(x1, fma(fma(-8.0, Float64(x1 * x1), 8.0), x2, fma(fma(24.0, x1, 12.0), x1, -12.0)), -6.0), x2, Float64(fma(fma(-19.0, x1, 9.0), x1, -2.0) * x1))); end return tmp end
code[x1_, x2_] := If[Or[LessEqual[x1, -1.0], N[Not[LessEqual[x1, 1.36]], $MachinePrecision]], N[(N[(N[(N[(N[(6.0 * x1 + -3.0), $MachinePrecision] * x1 + N[(N[(2.0 * x2 + -3.0), $MachinePrecision] * 4.0 + 9.0), $MachinePrecision]), $MachinePrecision] * x1), $MachinePrecision] * x1), $MachinePrecision] + x1), $MachinePrecision], N[(x1 + N[(N[(x1 * N[(N[(-8.0 * N[(x1 * x1), $MachinePrecision] + 8.0), $MachinePrecision] * x2 + N[(N[(24.0 * x1 + 12.0), $MachinePrecision] * x1 + -12.0), $MachinePrecision]), $MachinePrecision] + -6.0), $MachinePrecision] * x2 + N[(N[(N[(-19.0 * x1 + 9.0), $MachinePrecision] * x1 + -2.0), $MachinePrecision] * x1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x1 \leq -1 \lor \neg \left(x1 \leq 1.36\right):\\
\;\;\;\;\left(\mathsf{fma}\left(\mathsf{fma}\left(6, x1, -3\right), x1, \mathsf{fma}\left(\mathsf{fma}\left(2, x2, -3\right), 4, 9\right)\right) \cdot x1\right) \cdot x1 + x1\\
\mathbf{else}:\\
\;\;\;\;x1 + \mathsf{fma}\left(\mathsf{fma}\left(x1, \mathsf{fma}\left(\mathsf{fma}\left(-8, x1 \cdot x1, 8\right), x2, \mathsf{fma}\left(\mathsf{fma}\left(24, x1, 12\right), x1, -12\right)\right), -6\right), x2, \mathsf{fma}\left(\mathsf{fma}\left(-19, x1, 9\right), x1, -2\right) \cdot x1\right)\\
\end{array}
\end{array}
if x1 < -1 or 1.3600000000000001 < x1 Initial program 42.9%
Applied rewrites50.6%
Taylor expanded in x1 around -inf
*-commutativeN/A
lower-*.f64N/A
Applied rewrites91.2%
Taylor expanded in x1 around 0
Applied rewrites91.3%
lift-+.f64N/A
Applied rewrites91.3%
if -1 < x1 < 1.3600000000000001Initial program 99.5%
Taylor expanded in x1 around 0
Applied rewrites71.1%
Taylor expanded in x2 around 0
Applied rewrites86.6%
Taylor expanded in x2 around 0
Applied rewrites99.6%
Final simplification95.4%
(FPCore (x1 x2)
:precision binary64
(if (or (<= x1 -4300000000.0) (not (<= x1 3700000.0)))
(+
(* (* (fma (fma 6.0 x1 -3.0) x1 (fma (fma 2.0 x2 -3.0) 4.0 9.0)) x1) x1)
x1)
(+
x1
(fma
x2
-6.0
(*
(fma
(fma (fma 24.0 x1 12.0) x2 (fma -19.0 x1 9.0))
x1
(fma (* (fma 2.0 x2 -3.0) x2) 4.0 -2.0))
x1)))))
double code(double x1, double x2) {
double tmp;
if ((x1 <= -4300000000.0) || !(x1 <= 3700000.0)) {
tmp = ((fma(fma(6.0, x1, -3.0), x1, fma(fma(2.0, x2, -3.0), 4.0, 9.0)) * x1) * x1) + x1;
} else {
tmp = x1 + fma(x2, -6.0, (fma(fma(fma(24.0, x1, 12.0), x2, fma(-19.0, x1, 9.0)), x1, fma((fma(2.0, x2, -3.0) * x2), 4.0, -2.0)) * x1));
}
return tmp;
}
function code(x1, x2) tmp = 0.0 if ((x1 <= -4300000000.0) || !(x1 <= 3700000.0)) tmp = Float64(Float64(Float64(fma(fma(6.0, x1, -3.0), x1, fma(fma(2.0, x2, -3.0), 4.0, 9.0)) * x1) * x1) + x1); else tmp = Float64(x1 + fma(x2, -6.0, Float64(fma(fma(fma(24.0, x1, 12.0), x2, fma(-19.0, x1, 9.0)), x1, fma(Float64(fma(2.0, x2, -3.0) * x2), 4.0, -2.0)) * x1))); end return tmp end
code[x1_, x2_] := If[Or[LessEqual[x1, -4300000000.0], N[Not[LessEqual[x1, 3700000.0]], $MachinePrecision]], N[(N[(N[(N[(N[(6.0 * x1 + -3.0), $MachinePrecision] * x1 + N[(N[(2.0 * x2 + -3.0), $MachinePrecision] * 4.0 + 9.0), $MachinePrecision]), $MachinePrecision] * x1), $MachinePrecision] * x1), $MachinePrecision] + x1), $MachinePrecision], N[(x1 + N[(x2 * -6.0 + N[(N[(N[(N[(24.0 * x1 + 12.0), $MachinePrecision] * x2 + N[(-19.0 * x1 + 9.0), $MachinePrecision]), $MachinePrecision] * x1 + N[(N[(N[(2.0 * x2 + -3.0), $MachinePrecision] * x2), $MachinePrecision] * 4.0 + -2.0), $MachinePrecision]), $MachinePrecision] * x1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x1 \leq -4300000000 \lor \neg \left(x1 \leq 3700000\right):\\
\;\;\;\;\left(\mathsf{fma}\left(\mathsf{fma}\left(6, x1, -3\right), x1, \mathsf{fma}\left(\mathsf{fma}\left(2, x2, -3\right), 4, 9\right)\right) \cdot x1\right) \cdot x1 + x1\\
\mathbf{else}:\\
\;\;\;\;x1 + \mathsf{fma}\left(x2, -6, \mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(24, x1, 12\right), x2, \mathsf{fma}\left(-19, x1, 9\right)\right), x1, \mathsf{fma}\left(\mathsf{fma}\left(2, x2, -3\right) \cdot x2, 4, -2\right)\right) \cdot x1\right)\\
\end{array}
\end{array}
if x1 < -4.3e9 or 3.7e6 < x1 Initial program 40.2%
Applied rewrites48.3%
Taylor expanded in x1 around -inf
*-commutativeN/A
lower-*.f64N/A
Applied rewrites95.2%
Taylor expanded in x1 around 0
Applied rewrites95.2%
lift-+.f64N/A
Applied rewrites95.3%
if -4.3e9 < x1 < 3.7e6Initial program 99.5%
Taylor expanded in x1 around 0
Applied rewrites68.0%
Taylor expanded in x2 around 0
Applied rewrites85.2%
Applied rewrites85.3%
Final simplification90.1%
(FPCore (x1 x2)
:precision binary64
(let* ((t_0 (+ x1 (* (* (* x2 x2) x1) 8.0)))
(t_1 (+ x1 (* (fma (fma 6.0 x1 -3.0) x1 -3.0) (* x1 x1)))))
(if (<= x1 -4300000000.0)
t_1
(if (<= x1 -4.2e-88)
t_0
(if (<= x1 4.2e-93)
(+ x1 (* -6.0 x2))
(if (<= x1 3800000.0) t_0 t_1))))))
double code(double x1, double x2) {
double t_0 = x1 + (((x2 * x2) * x1) * 8.0);
double t_1 = x1 + (fma(fma(6.0, x1, -3.0), x1, -3.0) * (x1 * x1));
double tmp;
if (x1 <= -4300000000.0) {
tmp = t_1;
} else if (x1 <= -4.2e-88) {
tmp = t_0;
} else if (x1 <= 4.2e-93) {
tmp = x1 + (-6.0 * x2);
} else if (x1 <= 3800000.0) {
tmp = t_0;
} else {
tmp = t_1;
}
return tmp;
}
function code(x1, x2) t_0 = Float64(x1 + Float64(Float64(Float64(x2 * x2) * x1) * 8.0)) t_1 = Float64(x1 + Float64(fma(fma(6.0, x1, -3.0), x1, -3.0) * Float64(x1 * x1))) tmp = 0.0 if (x1 <= -4300000000.0) tmp = t_1; elseif (x1 <= -4.2e-88) tmp = t_0; elseif (x1 <= 4.2e-93) tmp = Float64(x1 + Float64(-6.0 * x2)); elseif (x1 <= 3800000.0) tmp = t_0; else tmp = t_1; end return tmp end
code[x1_, x2_] := Block[{t$95$0 = N[(x1 + N[(N[(N[(x2 * x2), $MachinePrecision] * x1), $MachinePrecision] * 8.0), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(x1 + N[(N[(N[(6.0 * x1 + -3.0), $MachinePrecision] * x1 + -3.0), $MachinePrecision] * N[(x1 * x1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[x1, -4300000000.0], t$95$1, If[LessEqual[x1, -4.2e-88], t$95$0, If[LessEqual[x1, 4.2e-93], N[(x1 + N[(-6.0 * x2), $MachinePrecision]), $MachinePrecision], If[LessEqual[x1, 3800000.0], t$95$0, t$95$1]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := x1 + \left(\left(x2 \cdot x2\right) \cdot x1\right) \cdot 8\\
t_1 := x1 + \mathsf{fma}\left(\mathsf{fma}\left(6, x1, -3\right), x1, -3\right) \cdot \left(x1 \cdot x1\right)\\
\mathbf{if}\;x1 \leq -4300000000:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;x1 \leq -4.2 \cdot 10^{-88}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;x1 \leq 4.2 \cdot 10^{-93}:\\
\;\;\;\;x1 + -6 \cdot x2\\
\mathbf{elif}\;x1 \leq 3800000:\\
\;\;\;\;t\_0\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if x1 < -4.3e9 or 3.8e6 < x1 Initial program 40.2%
Applied rewrites48.3%
Taylor expanded in x1 around -inf
*-commutativeN/A
lower-*.f64N/A
Applied rewrites95.2%
Taylor expanded in x1 around 0
Applied rewrites95.2%
Taylor expanded in x2 around 0
Applied rewrites93.9%
if -4.3e9 < x1 < -4.1999999999999999e-88 or 4.2000000000000002e-93 < x1 < 3.8e6Initial program 99.4%
Taylor expanded in x1 around 0
Applied rewrites53.5%
Taylor expanded in x2 around 0
Applied rewrites84.6%
Taylor expanded in x2 around inf
Applied rewrites46.0%
Taylor expanded in x1 around 0
Applied rewrites53.4%
if -4.1999999999999999e-88 < x1 < 4.2000000000000002e-93Initial program 99.6%
Taylor expanded in x1 around 0
lower-*.f6466.9
Applied rewrites66.9%
(FPCore (x1 x2)
:precision binary64
(if (or (<= x1 -4300000000.0) (not (<= x1 3700000.0)))
(+
(* (* (fma (fma 6.0 x1 -3.0) x1 (fma (fma 2.0 x2 -3.0) 4.0 9.0)) x1) x1)
x1)
(fma (fma (* (fma 2.0 x2 -3.0) x2) 4.0 -1.0) x1 (* -6.0 x2))))
double code(double x1, double x2) {
double tmp;
if ((x1 <= -4300000000.0) || !(x1 <= 3700000.0)) {
tmp = ((fma(fma(6.0, x1, -3.0), x1, fma(fma(2.0, x2, -3.0), 4.0, 9.0)) * x1) * x1) + x1;
} else {
tmp = fma(fma((fma(2.0, x2, -3.0) * x2), 4.0, -1.0), x1, (-6.0 * x2));
}
return tmp;
}
function code(x1, x2) tmp = 0.0 if ((x1 <= -4300000000.0) || !(x1 <= 3700000.0)) tmp = Float64(Float64(Float64(fma(fma(6.0, x1, -3.0), x1, fma(fma(2.0, x2, -3.0), 4.0, 9.0)) * x1) * x1) + x1); else tmp = fma(fma(Float64(fma(2.0, x2, -3.0) * x2), 4.0, -1.0), x1, Float64(-6.0 * x2)); end return tmp end
code[x1_, x2_] := If[Or[LessEqual[x1, -4300000000.0], N[Not[LessEqual[x1, 3700000.0]], $MachinePrecision]], N[(N[(N[(N[(N[(6.0 * x1 + -3.0), $MachinePrecision] * x1 + N[(N[(2.0 * x2 + -3.0), $MachinePrecision] * 4.0 + 9.0), $MachinePrecision]), $MachinePrecision] * x1), $MachinePrecision] * x1), $MachinePrecision] + x1), $MachinePrecision], N[(N[(N[(N[(2.0 * x2 + -3.0), $MachinePrecision] * x2), $MachinePrecision] * 4.0 + -1.0), $MachinePrecision] * x1 + N[(-6.0 * x2), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x1 \leq -4300000000 \lor \neg \left(x1 \leq 3700000\right):\\
\;\;\;\;\left(\mathsf{fma}\left(\mathsf{fma}\left(6, x1, -3\right), x1, \mathsf{fma}\left(\mathsf{fma}\left(2, x2, -3\right), 4, 9\right)\right) \cdot x1\right) \cdot x1 + x1\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(2, x2, -3\right) \cdot x2, 4, -1\right), x1, -6 \cdot x2\right)\\
\end{array}
\end{array}
if x1 < -4.3e9 or 3.7e6 < x1 Initial program 40.2%
Applied rewrites48.3%
Taylor expanded in x1 around -inf
*-commutativeN/A
lower-*.f64N/A
Applied rewrites95.2%
Taylor expanded in x1 around 0
Applied rewrites95.2%
lift-+.f64N/A
Applied rewrites95.3%
if -4.3e9 < x1 < 3.7e6Initial program 99.5%
Applied rewrites99.8%
Taylor expanded in x1 around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
sub-negN/A
*-commutativeN/A
metadata-evalN/A
lower-fma.f64N/A
*-commutativeN/A
lower-*.f64N/A
sub-negN/A
metadata-evalN/A
lower-fma.f64N/A
lower-*.f6484.6
Applied rewrites84.6%
Final simplification89.8%
(FPCore (x1 x2)
:precision binary64
(if (<= x1 -4300000000.0)
(+ x1 (* (fma (fma 6.0 x1 -3.0) x1 (* 8.0 x2)) (* x1 x1)))
(if (<= x1 3700000.0)
(fma (fma (* (fma 2.0 x2 -3.0) x2) 4.0 -1.0) x1 (* -6.0 x2))
(+
x1
(*
(fma (fma 6.0 x1 -3.0) x1 (fma (fma 2.0 x2 -3.0) 4.0 9.0))
(* x1 x1))))))
double code(double x1, double x2) {
double tmp;
if (x1 <= -4300000000.0) {
tmp = x1 + (fma(fma(6.0, x1, -3.0), x1, (8.0 * x2)) * (x1 * x1));
} else if (x1 <= 3700000.0) {
tmp = fma(fma((fma(2.0, x2, -3.0) * x2), 4.0, -1.0), x1, (-6.0 * x2));
} else {
tmp = x1 + (fma(fma(6.0, x1, -3.0), x1, fma(fma(2.0, x2, -3.0), 4.0, 9.0)) * (x1 * x1));
}
return tmp;
}
function code(x1, x2) tmp = 0.0 if (x1 <= -4300000000.0) tmp = Float64(x1 + Float64(fma(fma(6.0, x1, -3.0), x1, Float64(8.0 * x2)) * Float64(x1 * x1))); elseif (x1 <= 3700000.0) tmp = fma(fma(Float64(fma(2.0, x2, -3.0) * x2), 4.0, -1.0), x1, Float64(-6.0 * x2)); else tmp = Float64(x1 + Float64(fma(fma(6.0, x1, -3.0), x1, fma(fma(2.0, x2, -3.0), 4.0, 9.0)) * Float64(x1 * x1))); end return tmp end
code[x1_, x2_] := If[LessEqual[x1, -4300000000.0], N[(x1 + N[(N[(N[(6.0 * x1 + -3.0), $MachinePrecision] * x1 + N[(8.0 * x2), $MachinePrecision]), $MachinePrecision] * N[(x1 * x1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[x1, 3700000.0], N[(N[(N[(N[(2.0 * x2 + -3.0), $MachinePrecision] * x2), $MachinePrecision] * 4.0 + -1.0), $MachinePrecision] * x1 + N[(-6.0 * x2), $MachinePrecision]), $MachinePrecision], N[(x1 + N[(N[(N[(6.0 * x1 + -3.0), $MachinePrecision] * x1 + N[(N[(2.0 * x2 + -3.0), $MachinePrecision] * 4.0 + 9.0), $MachinePrecision]), $MachinePrecision] * N[(x1 * x1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x1 \leq -4300000000:\\
\;\;\;\;x1 + \mathsf{fma}\left(\mathsf{fma}\left(6, x1, -3\right), x1, 8 \cdot x2\right) \cdot \left(x1 \cdot x1\right)\\
\mathbf{elif}\;x1 \leq 3700000:\\
\;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(2, x2, -3\right) \cdot x2, 4, -1\right), x1, -6 \cdot x2\right)\\
\mathbf{else}:\\
\;\;\;\;x1 + \mathsf{fma}\left(\mathsf{fma}\left(6, x1, -3\right), x1, \mathsf{fma}\left(\mathsf{fma}\left(2, x2, -3\right), 4, 9\right)\right) \cdot \left(x1 \cdot x1\right)\\
\end{array}
\end{array}
if x1 < -4.3e9Initial program 26.1%
Applied rewrites42.5%
Taylor expanded in x1 around -inf
*-commutativeN/A
lower-*.f64N/A
Applied rewrites92.0%
Taylor expanded in x1 around 0
Applied rewrites92.1%
Taylor expanded in x2 around inf
Applied rewrites92.1%
if -4.3e9 < x1 < 3.7e6Initial program 99.5%
Applied rewrites99.8%
Taylor expanded in x1 around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
sub-negN/A
*-commutativeN/A
metadata-evalN/A
lower-fma.f64N/A
*-commutativeN/A
lower-*.f64N/A
sub-negN/A
metadata-evalN/A
lower-fma.f64N/A
lower-*.f6484.6
Applied rewrites84.6%
if 3.7e6 < x1 Initial program 53.9%
Applied rewrites53.9%
Taylor expanded in x1 around -inf
*-commutativeN/A
lower-*.f64N/A
Applied rewrites98.2%
Taylor expanded in x1 around 0
Applied rewrites98.3%
(FPCore (x1 x2)
:precision binary64
(let* ((t_0 (+ x1 (* (* (* x2 x2) x1) 8.0))))
(if (<= x1 -2.5e+98)
(+ x1 (* (fma (fma -19.0 x1 9.0) x1 -2.0) x1))
(if (<= x1 -4.2e-88)
t_0
(if (<= x1 4.2e-93)
(+ x1 (* -6.0 x2))
(if (<= x1 1.4e+100) t_0 (+ x1 (fma (* x1 x1) x1 (* -6.0 x2)))))))))
double code(double x1, double x2) {
double t_0 = x1 + (((x2 * x2) * x1) * 8.0);
double tmp;
if (x1 <= -2.5e+98) {
tmp = x1 + (fma(fma(-19.0, x1, 9.0), x1, -2.0) * x1);
} else if (x1 <= -4.2e-88) {
tmp = t_0;
} else if (x1 <= 4.2e-93) {
tmp = x1 + (-6.0 * x2);
} else if (x1 <= 1.4e+100) {
tmp = t_0;
} else {
tmp = x1 + fma((x1 * x1), x1, (-6.0 * x2));
}
return tmp;
}
function code(x1, x2) t_0 = Float64(x1 + Float64(Float64(Float64(x2 * x2) * x1) * 8.0)) tmp = 0.0 if (x1 <= -2.5e+98) tmp = Float64(x1 + Float64(fma(fma(-19.0, x1, 9.0), x1, -2.0) * x1)); elseif (x1 <= -4.2e-88) tmp = t_0; elseif (x1 <= 4.2e-93) tmp = Float64(x1 + Float64(-6.0 * x2)); elseif (x1 <= 1.4e+100) tmp = t_0; else tmp = Float64(x1 + fma(Float64(x1 * x1), x1, Float64(-6.0 * x2))); end return tmp end
code[x1_, x2_] := Block[{t$95$0 = N[(x1 + N[(N[(N[(x2 * x2), $MachinePrecision] * x1), $MachinePrecision] * 8.0), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[x1, -2.5e+98], N[(x1 + N[(N[(N[(-19.0 * x1 + 9.0), $MachinePrecision] * x1 + -2.0), $MachinePrecision] * x1), $MachinePrecision]), $MachinePrecision], If[LessEqual[x1, -4.2e-88], t$95$0, If[LessEqual[x1, 4.2e-93], N[(x1 + N[(-6.0 * x2), $MachinePrecision]), $MachinePrecision], If[LessEqual[x1, 1.4e+100], t$95$0, N[(x1 + N[(N[(x1 * x1), $MachinePrecision] * x1 + N[(-6.0 * x2), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := x1 + \left(\left(x2 \cdot x2\right) \cdot x1\right) \cdot 8\\
\mathbf{if}\;x1 \leq -2.5 \cdot 10^{+98}:\\
\;\;\;\;x1 + \mathsf{fma}\left(\mathsf{fma}\left(-19, x1, 9\right), x1, -2\right) \cdot x1\\
\mathbf{elif}\;x1 \leq -4.2 \cdot 10^{-88}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;x1 \leq 4.2 \cdot 10^{-93}:\\
\;\;\;\;x1 + -6 \cdot x2\\
\mathbf{elif}\;x1 \leq 1.4 \cdot 10^{+100}:\\
\;\;\;\;t\_0\\
\mathbf{else}:\\
\;\;\;\;x1 + \mathsf{fma}\left(x1 \cdot x1, x1, -6 \cdot x2\right)\\
\end{array}
\end{array}
if x1 < -2.4999999999999999e98Initial program 0.0%
Taylor expanded in x1 around 0
Applied rewrites73.3%
Taylor expanded in x2 around 0
Applied rewrites62.2%
Taylor expanded in x2 around 0
Applied rewrites100.0%
if -2.4999999999999999e98 < x1 < -4.1999999999999999e-88 or 4.2000000000000002e-93 < x1 < 1.3999999999999999e100Initial program 99.5%
Taylor expanded in x1 around 0
Applied rewrites33.8%
Taylor expanded in x2 around 0
Applied rewrites55.8%
Taylor expanded in x2 around inf
Applied rewrites29.6%
Taylor expanded in x1 around 0
Applied rewrites41.4%
if -4.1999999999999999e-88 < x1 < 4.2000000000000002e-93Initial program 99.6%
Taylor expanded in x1 around 0
lower-*.f6466.9
Applied rewrites66.9%
if 1.3999999999999999e100 < x1 Initial program 39.6%
Applied rewrites39.6%
Taylor expanded in x1 around 0
lower-*.f6496.4
Applied rewrites96.4%
(FPCore (x1 x2)
:precision binary64
(let* ((t_0 (+ x1 (* (* (* x2 x2) x1) 8.0))))
(if (<= x1 -2.8e+58)
(+ x1 (* (* (* x1 x1) x2) 8.0))
(if (<= x1 -4.2e-88)
t_0
(if (<= x1 4.2e-93)
(+ x1 (* -6.0 x2))
(if (<= x1 1.4e+100) t_0 (+ x1 (fma (* x1 x1) x1 (* -6.0 x2)))))))))
double code(double x1, double x2) {
double t_0 = x1 + (((x2 * x2) * x1) * 8.0);
double tmp;
if (x1 <= -2.8e+58) {
tmp = x1 + (((x1 * x1) * x2) * 8.0);
} else if (x1 <= -4.2e-88) {
tmp = t_0;
} else if (x1 <= 4.2e-93) {
tmp = x1 + (-6.0 * x2);
} else if (x1 <= 1.4e+100) {
tmp = t_0;
} else {
tmp = x1 + fma((x1 * x1), x1, (-6.0 * x2));
}
return tmp;
}
function code(x1, x2) t_0 = Float64(x1 + Float64(Float64(Float64(x2 * x2) * x1) * 8.0)) tmp = 0.0 if (x1 <= -2.8e+58) tmp = Float64(x1 + Float64(Float64(Float64(x1 * x1) * x2) * 8.0)); elseif (x1 <= -4.2e-88) tmp = t_0; elseif (x1 <= 4.2e-93) tmp = Float64(x1 + Float64(-6.0 * x2)); elseif (x1 <= 1.4e+100) tmp = t_0; else tmp = Float64(x1 + fma(Float64(x1 * x1), x1, Float64(-6.0 * x2))); end return tmp end
code[x1_, x2_] := Block[{t$95$0 = N[(x1 + N[(N[(N[(x2 * x2), $MachinePrecision] * x1), $MachinePrecision] * 8.0), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[x1, -2.8e+58], N[(x1 + N[(N[(N[(x1 * x1), $MachinePrecision] * x2), $MachinePrecision] * 8.0), $MachinePrecision]), $MachinePrecision], If[LessEqual[x1, -4.2e-88], t$95$0, If[LessEqual[x1, 4.2e-93], N[(x1 + N[(-6.0 * x2), $MachinePrecision]), $MachinePrecision], If[LessEqual[x1, 1.4e+100], t$95$0, N[(x1 + N[(N[(x1 * x1), $MachinePrecision] * x1 + N[(-6.0 * x2), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := x1 + \left(\left(x2 \cdot x2\right) \cdot x1\right) \cdot 8\\
\mathbf{if}\;x1 \leq -2.8 \cdot 10^{+58}:\\
\;\;\;\;x1 + \left(\left(x1 \cdot x1\right) \cdot x2\right) \cdot 8\\
\mathbf{elif}\;x1 \leq -4.2 \cdot 10^{-88}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;x1 \leq 4.2 \cdot 10^{-93}:\\
\;\;\;\;x1 + -6 \cdot x2\\
\mathbf{elif}\;x1 \leq 1.4 \cdot 10^{+100}:\\
\;\;\;\;t\_0\\
\mathbf{else}:\\
\;\;\;\;x1 + \mathsf{fma}\left(x1 \cdot x1, x1, -6 \cdot x2\right)\\
\end{array}
\end{array}
if x1 < -2.7999999999999998e58Initial program 10.0%
Applied rewrites30.0%
Taylor expanded in x1 around -inf
*-commutativeN/A
lower-*.f64N/A
Applied rewrites100.0%
Taylor expanded in x2 around inf
Applied rewrites38.4%
if -2.7999999999999998e58 < x1 < -4.1999999999999999e-88 or 4.2000000000000002e-93 < x1 < 1.3999999999999999e100Initial program 99.5%
Taylor expanded in x1 around 0
Applied rewrites33.1%
Taylor expanded in x2 around 0
Applied rewrites58.4%
Taylor expanded in x2 around inf
Applied rewrites28.6%
Taylor expanded in x1 around 0
Applied rewrites43.0%
if -4.1999999999999999e-88 < x1 < 4.2000000000000002e-93Initial program 99.6%
Taylor expanded in x1 around 0
lower-*.f6466.9
Applied rewrites66.9%
if 1.3999999999999999e100 < x1 Initial program 39.6%
Applied rewrites39.6%
Taylor expanded in x1 around 0
lower-*.f6496.4
Applied rewrites96.4%
(FPCore (x1 x2) :precision binary64 (if (or (<= x1 -4300000000.0) (not (<= x1 3800000.0))) (+ x1 (* (fma (fma 6.0 x1 -3.0) x1 (* 8.0 x2)) (* x1 x1))) (fma (fma (* (fma 2.0 x2 -3.0) x2) 4.0 -1.0) x1 (* -6.0 x2))))
double code(double x1, double x2) {
double tmp;
if ((x1 <= -4300000000.0) || !(x1 <= 3800000.0)) {
tmp = x1 + (fma(fma(6.0, x1, -3.0), x1, (8.0 * x2)) * (x1 * x1));
} else {
tmp = fma(fma((fma(2.0, x2, -3.0) * x2), 4.0, -1.0), x1, (-6.0 * x2));
}
return tmp;
}
function code(x1, x2) tmp = 0.0 if ((x1 <= -4300000000.0) || !(x1 <= 3800000.0)) tmp = Float64(x1 + Float64(fma(fma(6.0, x1, -3.0), x1, Float64(8.0 * x2)) * Float64(x1 * x1))); else tmp = fma(fma(Float64(fma(2.0, x2, -3.0) * x2), 4.0, -1.0), x1, Float64(-6.0 * x2)); end return tmp end
code[x1_, x2_] := If[Or[LessEqual[x1, -4300000000.0], N[Not[LessEqual[x1, 3800000.0]], $MachinePrecision]], N[(x1 + N[(N[(N[(6.0 * x1 + -3.0), $MachinePrecision] * x1 + N[(8.0 * x2), $MachinePrecision]), $MachinePrecision] * N[(x1 * x1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(N[(N[(2.0 * x2 + -3.0), $MachinePrecision] * x2), $MachinePrecision] * 4.0 + -1.0), $MachinePrecision] * x1 + N[(-6.0 * x2), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x1 \leq -4300000000 \lor \neg \left(x1 \leq 3800000\right):\\
\;\;\;\;x1 + \mathsf{fma}\left(\mathsf{fma}\left(6, x1, -3\right), x1, 8 \cdot x2\right) \cdot \left(x1 \cdot x1\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(2, x2, -3\right) \cdot x2, 4, -1\right), x1, -6 \cdot x2\right)\\
\end{array}
\end{array}
if x1 < -4.3e9 or 3.8e6 < x1 Initial program 40.2%
Applied rewrites48.3%
Taylor expanded in x1 around -inf
*-commutativeN/A
lower-*.f64N/A
Applied rewrites95.2%
Taylor expanded in x1 around 0
Applied rewrites95.2%
Taylor expanded in x2 around inf
Applied rewrites95.1%
if -4.3e9 < x1 < 3.8e6Initial program 99.5%
Applied rewrites99.8%
Taylor expanded in x1 around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
sub-negN/A
*-commutativeN/A
metadata-evalN/A
lower-fma.f64N/A
*-commutativeN/A
lower-*.f64N/A
sub-negN/A
metadata-evalN/A
lower-fma.f64N/A
lower-*.f6484.6
Applied rewrites84.6%
Final simplification89.7%
(FPCore (x1 x2) :precision binary64 (if (or (<= x1 -4300000000.0) (not (<= x1 3800000.0))) (+ x1 (* (fma (fma 6.0 x1 -3.0) x1 -3.0) (* x1 x1))) (fma (fma (* (fma 2.0 x2 -3.0) x2) 4.0 -1.0) x1 (* -6.0 x2))))
double code(double x1, double x2) {
double tmp;
if ((x1 <= -4300000000.0) || !(x1 <= 3800000.0)) {
tmp = x1 + (fma(fma(6.0, x1, -3.0), x1, -3.0) * (x1 * x1));
} else {
tmp = fma(fma((fma(2.0, x2, -3.0) * x2), 4.0, -1.0), x1, (-6.0 * x2));
}
return tmp;
}
function code(x1, x2) tmp = 0.0 if ((x1 <= -4300000000.0) || !(x1 <= 3800000.0)) tmp = Float64(x1 + Float64(fma(fma(6.0, x1, -3.0), x1, -3.0) * Float64(x1 * x1))); else tmp = fma(fma(Float64(fma(2.0, x2, -3.0) * x2), 4.0, -1.0), x1, Float64(-6.0 * x2)); end return tmp end
code[x1_, x2_] := If[Or[LessEqual[x1, -4300000000.0], N[Not[LessEqual[x1, 3800000.0]], $MachinePrecision]], N[(x1 + N[(N[(N[(6.0 * x1 + -3.0), $MachinePrecision] * x1 + -3.0), $MachinePrecision] * N[(x1 * x1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(N[(N[(2.0 * x2 + -3.0), $MachinePrecision] * x2), $MachinePrecision] * 4.0 + -1.0), $MachinePrecision] * x1 + N[(-6.0 * x2), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x1 \leq -4300000000 \lor \neg \left(x1 \leq 3800000\right):\\
\;\;\;\;x1 + \mathsf{fma}\left(\mathsf{fma}\left(6, x1, -3\right), x1, -3\right) \cdot \left(x1 \cdot x1\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(2, x2, -3\right) \cdot x2, 4, -1\right), x1, -6 \cdot x2\right)\\
\end{array}
\end{array}
if x1 < -4.3e9 or 3.8e6 < x1 Initial program 40.2%
Applied rewrites48.3%
Taylor expanded in x1 around -inf
*-commutativeN/A
lower-*.f64N/A
Applied rewrites95.2%
Taylor expanded in x1 around 0
Applied rewrites95.2%
Taylor expanded in x2 around 0
Applied rewrites93.9%
if -4.3e9 < x1 < 3.8e6Initial program 99.5%
Applied rewrites99.8%
Taylor expanded in x1 around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
sub-negN/A
*-commutativeN/A
metadata-evalN/A
lower-fma.f64N/A
*-commutativeN/A
lower-*.f64N/A
sub-negN/A
metadata-evalN/A
lower-fma.f64N/A
lower-*.f6484.6
Applied rewrites84.6%
Final simplification89.1%
(FPCore (x1 x2) :precision binary64 (+ x1 (* -6.0 x2)))
double code(double x1, double x2) {
return x1 + (-6.0 * x2);
}
real(8) function code(x1, x2)
real(8), intent (in) :: x1
real(8), intent (in) :: x2
code = x1 + ((-6.0d0) * x2)
end function
public static double code(double x1, double x2) {
return x1 + (-6.0 * x2);
}
def code(x1, x2): return x1 + (-6.0 * x2)
function code(x1, x2) return Float64(x1 + Float64(-6.0 * x2)) end
function tmp = code(x1, x2) tmp = x1 + (-6.0 * x2); end
code[x1_, x2_] := N[(x1 + N[(-6.0 * x2), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
x1 + -6 \cdot x2
\end{array}
Initial program 70.8%
Taylor expanded in x1 around 0
lower-*.f6427.6
Applied rewrites27.6%
herbie shell --seed 2024317
(FPCore (x1 x2)
:name "Rosa's FloatVsDoubleBenchmark"
:precision binary64
(+ x1 (+ (+ (+ (+ (* (+ (* (* (* 2.0 x1) (/ (- (+ (* (* 3.0 x1) x1) (* 2.0 x2)) x1) (+ (* x1 x1) 1.0))) (- (/ (- (+ (* (* 3.0 x1) x1) (* 2.0 x2)) x1) (+ (* x1 x1) 1.0)) 3.0)) (* (* x1 x1) (- (* 4.0 (/ (- (+ (* (* 3.0 x1) x1) (* 2.0 x2)) x1) (+ (* x1 x1) 1.0))) 6.0))) (+ (* x1 x1) 1.0)) (* (* (* 3.0 x1) x1) (/ (- (+ (* (* 3.0 x1) x1) (* 2.0 x2)) x1) (+ (* x1 x1) 1.0)))) (* (* x1 x1) x1)) x1) (* 3.0 (/ (- (- (* (* 3.0 x1) x1) (* 2.0 x2)) x1) (+ (* x1 x1) 1.0))))))