
(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 23 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 (* x1 (* x1 3.0)))
(t_1 (- (* 2.0 x2) (* x1 (- 1.0 (* x1 3.0)))))
(t_2 (+ (* x1 x1) 1.0))
(t_3 (/ (- (+ t_0 (* 2.0 x2)) x1) t_2)))
(if (<=
(+
x1
(+
(+
x1
(+
(+
(*
t_2
(+
(* (* (* x1 2.0) t_3) (- t_3 3.0))
(* (* x1 x1) (- (* t_3 4.0) 6.0))))
(* t_0 t_3))
(* x1 (* x1 x1))))
(* 3.0 (/ (- (- t_0 (* 2.0 x2)) x1) t_2))))
INFINITY)
(+
(+ (* x1 -3.0) (* x2 -6.0))
(+
(+
(* x1 2.0)
(*
t_2
(+
(* (* x1 x1) -6.0)
(/
(* x1 (+ (- -6.0 (/ t_1 (/ (- -1.0 (* x1 x1)) 2.0))) (* x1 4.0)))
(/ t_2 t_1)))))
(* x1 (* x1 (+ x1 9.0)))))
(*
(* (* x1 x1) (* x1 x1))
(+ 6.0 (/ (- (/ (- 15.0 (+ 18.0 (* x2 -8.0))) x1) 3.0) x1))))))
double code(double x1, double x2) {
double t_0 = x1 * (x1 * 3.0);
double t_1 = (2.0 * x2) - (x1 * (1.0 - (x1 * 3.0)));
double t_2 = (x1 * x1) + 1.0;
double t_3 = ((t_0 + (2.0 * x2)) - x1) / t_2;
double tmp;
if ((x1 + ((x1 + (((t_2 * ((((x1 * 2.0) * t_3) * (t_3 - 3.0)) + ((x1 * x1) * ((t_3 * 4.0) - 6.0)))) + (t_0 * t_3)) + (x1 * (x1 * x1)))) + (3.0 * (((t_0 - (2.0 * x2)) - x1) / t_2)))) <= ((double) INFINITY)) {
tmp = ((x1 * -3.0) + (x2 * -6.0)) + (((x1 * 2.0) + (t_2 * (((x1 * x1) * -6.0) + ((x1 * ((-6.0 - (t_1 / ((-1.0 - (x1 * x1)) / 2.0))) + (x1 * 4.0))) / (t_2 / t_1))))) + (x1 * (x1 * (x1 + 9.0))));
} else {
tmp = ((x1 * x1) * (x1 * x1)) * (6.0 + ((((15.0 - (18.0 + (x2 * -8.0))) / x1) - 3.0) / x1));
}
return tmp;
}
public static double code(double x1, double x2) {
double t_0 = x1 * (x1 * 3.0);
double t_1 = (2.0 * x2) - (x1 * (1.0 - (x1 * 3.0)));
double t_2 = (x1 * x1) + 1.0;
double t_3 = ((t_0 + (2.0 * x2)) - x1) / t_2;
double tmp;
if ((x1 + ((x1 + (((t_2 * ((((x1 * 2.0) * t_3) * (t_3 - 3.0)) + ((x1 * x1) * ((t_3 * 4.0) - 6.0)))) + (t_0 * t_3)) + (x1 * (x1 * x1)))) + (3.0 * (((t_0 - (2.0 * x2)) - x1) / t_2)))) <= Double.POSITIVE_INFINITY) {
tmp = ((x1 * -3.0) + (x2 * -6.0)) + (((x1 * 2.0) + (t_2 * (((x1 * x1) * -6.0) + ((x1 * ((-6.0 - (t_1 / ((-1.0 - (x1 * x1)) / 2.0))) + (x1 * 4.0))) / (t_2 / t_1))))) + (x1 * (x1 * (x1 + 9.0))));
} else {
tmp = ((x1 * x1) * (x1 * x1)) * (6.0 + ((((15.0 - (18.0 + (x2 * -8.0))) / x1) - 3.0) / x1));
}
return tmp;
}
def code(x1, x2): t_0 = x1 * (x1 * 3.0) t_1 = (2.0 * x2) - (x1 * (1.0 - (x1 * 3.0))) t_2 = (x1 * x1) + 1.0 t_3 = ((t_0 + (2.0 * x2)) - x1) / t_2 tmp = 0 if (x1 + ((x1 + (((t_2 * ((((x1 * 2.0) * t_3) * (t_3 - 3.0)) + ((x1 * x1) * ((t_3 * 4.0) - 6.0)))) + (t_0 * t_3)) + (x1 * (x1 * x1)))) + (3.0 * (((t_0 - (2.0 * x2)) - x1) / t_2)))) <= math.inf: tmp = ((x1 * -3.0) + (x2 * -6.0)) + (((x1 * 2.0) + (t_2 * (((x1 * x1) * -6.0) + ((x1 * ((-6.0 - (t_1 / ((-1.0 - (x1 * x1)) / 2.0))) + (x1 * 4.0))) / (t_2 / t_1))))) + (x1 * (x1 * (x1 + 9.0)))) else: tmp = ((x1 * x1) * (x1 * x1)) * (6.0 + ((((15.0 - (18.0 + (x2 * -8.0))) / x1) - 3.0) / x1)) return tmp
function code(x1, x2) t_0 = Float64(x1 * Float64(x1 * 3.0)) t_1 = Float64(Float64(2.0 * x2) - Float64(x1 * Float64(1.0 - Float64(x1 * 3.0)))) t_2 = Float64(Float64(x1 * x1) + 1.0) t_3 = Float64(Float64(Float64(t_0 + Float64(2.0 * x2)) - x1) / t_2) tmp = 0.0 if (Float64(x1 + Float64(Float64(x1 + Float64(Float64(Float64(t_2 * Float64(Float64(Float64(Float64(x1 * 2.0) * t_3) * Float64(t_3 - 3.0)) + Float64(Float64(x1 * x1) * Float64(Float64(t_3 * 4.0) - 6.0)))) + Float64(t_0 * t_3)) + Float64(x1 * Float64(x1 * x1)))) + Float64(3.0 * Float64(Float64(Float64(t_0 - Float64(2.0 * x2)) - x1) / t_2)))) <= Inf) tmp = Float64(Float64(Float64(x1 * -3.0) + Float64(x2 * -6.0)) + Float64(Float64(Float64(x1 * 2.0) + Float64(t_2 * Float64(Float64(Float64(x1 * x1) * -6.0) + Float64(Float64(x1 * Float64(Float64(-6.0 - Float64(t_1 / Float64(Float64(-1.0 - Float64(x1 * x1)) / 2.0))) + Float64(x1 * 4.0))) / Float64(t_2 / t_1))))) + Float64(x1 * Float64(x1 * Float64(x1 + 9.0))))); else tmp = Float64(Float64(Float64(x1 * x1) * Float64(x1 * x1)) * Float64(6.0 + Float64(Float64(Float64(Float64(15.0 - Float64(18.0 + Float64(x2 * -8.0))) / x1) - 3.0) / x1))); end return tmp end
function tmp_2 = code(x1, x2) t_0 = x1 * (x1 * 3.0); t_1 = (2.0 * x2) - (x1 * (1.0 - (x1 * 3.0))); t_2 = (x1 * x1) + 1.0; t_3 = ((t_0 + (2.0 * x2)) - x1) / t_2; tmp = 0.0; if ((x1 + ((x1 + (((t_2 * ((((x1 * 2.0) * t_3) * (t_3 - 3.0)) + ((x1 * x1) * ((t_3 * 4.0) - 6.0)))) + (t_0 * t_3)) + (x1 * (x1 * x1)))) + (3.0 * (((t_0 - (2.0 * x2)) - x1) / t_2)))) <= Inf) tmp = ((x1 * -3.0) + (x2 * -6.0)) + (((x1 * 2.0) + (t_2 * (((x1 * x1) * -6.0) + ((x1 * ((-6.0 - (t_1 / ((-1.0 - (x1 * x1)) / 2.0))) + (x1 * 4.0))) / (t_2 / t_1))))) + (x1 * (x1 * (x1 + 9.0)))); else tmp = ((x1 * x1) * (x1 * x1)) * (6.0 + ((((15.0 - (18.0 + (x2 * -8.0))) / x1) - 3.0) / x1)); end tmp_2 = tmp; end
code[x1_, x2_] := Block[{t$95$0 = N[(x1 * N[(x1 * 3.0), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(N[(2.0 * x2), $MachinePrecision] - N[(x1 * N[(1.0 - N[(x1 * 3.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(N[(x1 * x1), $MachinePrecision] + 1.0), $MachinePrecision]}, Block[{t$95$3 = N[(N[(N[(t$95$0 + N[(2.0 * x2), $MachinePrecision]), $MachinePrecision] - x1), $MachinePrecision] / t$95$2), $MachinePrecision]}, If[LessEqual[N[(x1 + N[(N[(x1 + N[(N[(N[(t$95$2 * N[(N[(N[(N[(x1 * 2.0), $MachinePrecision] * t$95$3), $MachinePrecision] * N[(t$95$3 - 3.0), $MachinePrecision]), $MachinePrecision] + N[(N[(x1 * x1), $MachinePrecision] * N[(N[(t$95$3 * 4.0), $MachinePrecision] - 6.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(t$95$0 * t$95$3), $MachinePrecision]), $MachinePrecision] + N[(x1 * N[(x1 * x1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(3.0 * N[(N[(N[(t$95$0 - N[(2.0 * x2), $MachinePrecision]), $MachinePrecision] - x1), $MachinePrecision] / t$95$2), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], Infinity], N[(N[(N[(x1 * -3.0), $MachinePrecision] + N[(x2 * -6.0), $MachinePrecision]), $MachinePrecision] + N[(N[(N[(x1 * 2.0), $MachinePrecision] + N[(t$95$2 * N[(N[(N[(x1 * x1), $MachinePrecision] * -6.0), $MachinePrecision] + N[(N[(x1 * N[(N[(-6.0 - N[(t$95$1 / N[(N[(-1.0 - N[(x1 * x1), $MachinePrecision]), $MachinePrecision] / 2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(x1 * 4.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(t$95$2 / t$95$1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(x1 * N[(x1 * N[(x1 + 9.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(N[(x1 * x1), $MachinePrecision] * N[(x1 * x1), $MachinePrecision]), $MachinePrecision] * N[(6.0 + N[(N[(N[(N[(15.0 - N[(18.0 + N[(x2 * -8.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / x1), $MachinePrecision] - 3.0), $MachinePrecision] / x1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := x1 \cdot \left(x1 \cdot 3\right)\\
t_1 := 2 \cdot x2 - x1 \cdot \left(1 - x1 \cdot 3\right)\\
t_2 := x1 \cdot x1 + 1\\
t_3 := \frac{\left(t\_0 + 2 \cdot x2\right) - x1}{t\_2}\\
\mathbf{if}\;x1 + \left(\left(x1 + \left(\left(t\_2 \cdot \left(\left(\left(x1 \cdot 2\right) \cdot t\_3\right) \cdot \left(t\_3 - 3\right) + \left(x1 \cdot x1\right) \cdot \left(t\_3 \cdot 4 - 6\right)\right) + t\_0 \cdot t\_3\right) + x1 \cdot \left(x1 \cdot x1\right)\right)\right) + 3 \cdot \frac{\left(t\_0 - 2 \cdot x2\right) - x1}{t\_2}\right) \leq \infty:\\
\;\;\;\;\left(x1 \cdot -3 + x2 \cdot -6\right) + \left(\left(x1 \cdot 2 + t\_2 \cdot \left(\left(x1 \cdot x1\right) \cdot -6 + \frac{x1 \cdot \left(\left(-6 - \frac{t\_1}{\frac{-1 - x1 \cdot x1}{2}}\right) + x1 \cdot 4\right)}{\frac{t\_2}{t\_1}}\right)\right) + x1 \cdot \left(x1 \cdot \left(x1 + 9\right)\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\left(\left(x1 \cdot x1\right) \cdot \left(x1 \cdot x1\right)\right) \cdot \left(6 + \frac{\frac{15 - \left(18 + x2 \cdot -8\right)}{x1} - 3}{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.5%
Simplified99.5%
Applied egg-rr99.5%
Taylor expanded in x1 around inf
Simplified97.5%
Taylor expanded in x1 around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f6499.5%
Simplified99.5%
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%
Simplified0.0%
Applied egg-rr0.0%
Taylor expanded in x1 around -inf
*-lowering-*.f64N/A
metadata-evalN/A
pow-sqrN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
/-lowering-/.f64N/A
Simplified100.0%
Final simplification99.7%
(FPCore (x1 x2)
:precision binary64
(let* ((t_0 (+ (* x1 x1) 1.0))
(t_1
(*
(* (* x1 x1) (* x1 x1))
(+ 6.0 (/ (- (/ (- 15.0 (+ 18.0 (* x2 -8.0))) x1) 3.0) x1))))
(t_2 (* x1 (- 1.0 (* x1 3.0)))))
(if (<= x1 -2.7e+19)
t_1
(if (<= x1 3000000000000.0)
(+
(/ (* 3.0 (- (* x2 -2.0) t_2)) t_0)
(+
(* t_0 (/ (* 8.0 (* x1 (* x2 x2))) (* t_0 t_0)))
(*
x1
(+
2.0
(* x1 (- x1 (* (- (* 2.0 x2) t_2) (/ 3.0 (- -1.0 (* x1 x1))))))))))
t_1))))
double code(double x1, double x2) {
double t_0 = (x1 * x1) + 1.0;
double t_1 = ((x1 * x1) * (x1 * x1)) * (6.0 + ((((15.0 - (18.0 + (x2 * -8.0))) / x1) - 3.0) / x1));
double t_2 = x1 * (1.0 - (x1 * 3.0));
double tmp;
if (x1 <= -2.7e+19) {
tmp = t_1;
} else if (x1 <= 3000000000000.0) {
tmp = ((3.0 * ((x2 * -2.0) - t_2)) / t_0) + ((t_0 * ((8.0 * (x1 * (x2 * x2))) / (t_0 * t_0))) + (x1 * (2.0 + (x1 * (x1 - (((2.0 * x2) - t_2) * (3.0 / (-1.0 - (x1 * x1)))))))));
} else {
tmp = t_1;
}
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 = (x1 * x1) + 1.0d0
t_1 = ((x1 * x1) * (x1 * x1)) * (6.0d0 + ((((15.0d0 - (18.0d0 + (x2 * (-8.0d0)))) / x1) - 3.0d0) / x1))
t_2 = x1 * (1.0d0 - (x1 * 3.0d0))
if (x1 <= (-2.7d+19)) then
tmp = t_1
else if (x1 <= 3000000000000.0d0) then
tmp = ((3.0d0 * ((x2 * (-2.0d0)) - t_2)) / t_0) + ((t_0 * ((8.0d0 * (x1 * (x2 * x2))) / (t_0 * t_0))) + (x1 * (2.0d0 + (x1 * (x1 - (((2.0d0 * x2) - t_2) * (3.0d0 / ((-1.0d0) - (x1 * x1)))))))))
else
tmp = t_1
end if
code = tmp
end function
public static double code(double x1, double x2) {
double t_0 = (x1 * x1) + 1.0;
double t_1 = ((x1 * x1) * (x1 * x1)) * (6.0 + ((((15.0 - (18.0 + (x2 * -8.0))) / x1) - 3.0) / x1));
double t_2 = x1 * (1.0 - (x1 * 3.0));
double tmp;
if (x1 <= -2.7e+19) {
tmp = t_1;
} else if (x1 <= 3000000000000.0) {
tmp = ((3.0 * ((x2 * -2.0) - t_2)) / t_0) + ((t_0 * ((8.0 * (x1 * (x2 * x2))) / (t_0 * t_0))) + (x1 * (2.0 + (x1 * (x1 - (((2.0 * x2) - t_2) * (3.0 / (-1.0 - (x1 * x1)))))))));
} else {
tmp = t_1;
}
return tmp;
}
def code(x1, x2): t_0 = (x1 * x1) + 1.0 t_1 = ((x1 * x1) * (x1 * x1)) * (6.0 + ((((15.0 - (18.0 + (x2 * -8.0))) / x1) - 3.0) / x1)) t_2 = x1 * (1.0 - (x1 * 3.0)) tmp = 0 if x1 <= -2.7e+19: tmp = t_1 elif x1 <= 3000000000000.0: tmp = ((3.0 * ((x2 * -2.0) - t_2)) / t_0) + ((t_0 * ((8.0 * (x1 * (x2 * x2))) / (t_0 * t_0))) + (x1 * (2.0 + (x1 * (x1 - (((2.0 * x2) - t_2) * (3.0 / (-1.0 - (x1 * x1))))))))) else: tmp = t_1 return tmp
function code(x1, x2) t_0 = Float64(Float64(x1 * x1) + 1.0) t_1 = Float64(Float64(Float64(x1 * x1) * Float64(x1 * x1)) * Float64(6.0 + Float64(Float64(Float64(Float64(15.0 - Float64(18.0 + Float64(x2 * -8.0))) / x1) - 3.0) / x1))) t_2 = Float64(x1 * Float64(1.0 - Float64(x1 * 3.0))) tmp = 0.0 if (x1 <= -2.7e+19) tmp = t_1; elseif (x1 <= 3000000000000.0) tmp = Float64(Float64(Float64(3.0 * Float64(Float64(x2 * -2.0) - t_2)) / t_0) + Float64(Float64(t_0 * Float64(Float64(8.0 * Float64(x1 * Float64(x2 * x2))) / Float64(t_0 * t_0))) + Float64(x1 * Float64(2.0 + Float64(x1 * Float64(x1 - Float64(Float64(Float64(2.0 * x2) - t_2) * Float64(3.0 / Float64(-1.0 - Float64(x1 * x1)))))))))); else tmp = t_1; end return tmp end
function tmp_2 = code(x1, x2) t_0 = (x1 * x1) + 1.0; t_1 = ((x1 * x1) * (x1 * x1)) * (6.0 + ((((15.0 - (18.0 + (x2 * -8.0))) / x1) - 3.0) / x1)); t_2 = x1 * (1.0 - (x1 * 3.0)); tmp = 0.0; if (x1 <= -2.7e+19) tmp = t_1; elseif (x1 <= 3000000000000.0) tmp = ((3.0 * ((x2 * -2.0) - t_2)) / t_0) + ((t_0 * ((8.0 * (x1 * (x2 * x2))) / (t_0 * t_0))) + (x1 * (2.0 + (x1 * (x1 - (((2.0 * x2) - t_2) * (3.0 / (-1.0 - (x1 * x1))))))))); else tmp = t_1; end tmp_2 = tmp; end
code[x1_, x2_] := Block[{t$95$0 = N[(N[(x1 * x1), $MachinePrecision] + 1.0), $MachinePrecision]}, Block[{t$95$1 = N[(N[(N[(x1 * x1), $MachinePrecision] * N[(x1 * x1), $MachinePrecision]), $MachinePrecision] * N[(6.0 + N[(N[(N[(N[(15.0 - N[(18.0 + N[(x2 * -8.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / x1), $MachinePrecision] - 3.0), $MachinePrecision] / x1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(x1 * N[(1.0 - N[(x1 * 3.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[x1, -2.7e+19], t$95$1, If[LessEqual[x1, 3000000000000.0], N[(N[(N[(3.0 * N[(N[(x2 * -2.0), $MachinePrecision] - t$95$2), $MachinePrecision]), $MachinePrecision] / t$95$0), $MachinePrecision] + N[(N[(t$95$0 * N[(N[(8.0 * N[(x1 * N[(x2 * x2), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(t$95$0 * t$95$0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(x1 * N[(2.0 + N[(x1 * N[(x1 - N[(N[(N[(2.0 * x2), $MachinePrecision] - t$95$2), $MachinePrecision] * N[(3.0 / N[(-1.0 - N[(x1 * x1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$1]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := x1 \cdot x1 + 1\\
t_1 := \left(\left(x1 \cdot x1\right) \cdot \left(x1 \cdot x1\right)\right) \cdot \left(6 + \frac{\frac{15 - \left(18 + x2 \cdot -8\right)}{x1} - 3}{x1}\right)\\
t_2 := x1 \cdot \left(1 - x1 \cdot 3\right)\\
\mathbf{if}\;x1 \leq -2.7 \cdot 10^{+19}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;x1 \leq 3000000000000:\\
\;\;\;\;\frac{3 \cdot \left(x2 \cdot -2 - t\_2\right)}{t\_0} + \left(t\_0 \cdot \frac{8 \cdot \left(x1 \cdot \left(x2 \cdot x2\right)\right)}{t\_0 \cdot t\_0} + x1 \cdot \left(2 + x1 \cdot \left(x1 - \left(2 \cdot x2 - t\_2\right) \cdot \frac{3}{-1 - x1 \cdot x1}\right)\right)\right)\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if x1 < -2.7e19 or 3e12 < x1 Initial program 35.6%
Simplified35.6%
Applied egg-rr35.6%
Taylor expanded in x1 around -inf
*-lowering-*.f64N/A
metadata-evalN/A
pow-sqrN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
/-lowering-/.f64N/A
Simplified97.0%
if -2.7e19 < x1 < 3e12Initial program 99.5%
Simplified99.5%
Taylor expanded in x2 around inf
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
unpow2N/A
*-lowering-*.f6488.4%
Simplified88.4%
Final simplification92.0%
(FPCore (x1 x2)
:precision binary64
(let* ((t_0
(*
(* (* x1 x1) (* x1 x1))
(+ 6.0 (/ (- (/ (- 15.0 (+ 18.0 (* x2 -8.0))) x1) 3.0) x1)))))
(if (<= x1 -1.5e+19)
t_0
(if (<= x1 580.0)
(+
(* x2 -6.0)
(*
x1
(+
(* (* 2.0 x2) (+ -6.0 (* x2 4.0)))
(+
-1.0
(*
x1
(+
-6.0
(+
(+ 6.0 (* x2 -4.0))
(+ (* 3.0 (+ 3.0 (* 2.0 x2))) (* x2 10.0)))))))))
t_0))))
double code(double x1, double x2) {
double t_0 = ((x1 * x1) * (x1 * x1)) * (6.0 + ((((15.0 - (18.0 + (x2 * -8.0))) / x1) - 3.0) / x1));
double tmp;
if (x1 <= -1.5e+19) {
tmp = t_0;
} else if (x1 <= 580.0) {
tmp = (x2 * -6.0) + (x1 * (((2.0 * x2) * (-6.0 + (x2 * 4.0))) + (-1.0 + (x1 * (-6.0 + ((6.0 + (x2 * -4.0)) + ((3.0 * (3.0 + (2.0 * x2))) + (x2 * 10.0))))))));
} else {
tmp = t_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) :: tmp
t_0 = ((x1 * x1) * (x1 * x1)) * (6.0d0 + ((((15.0d0 - (18.0d0 + (x2 * (-8.0d0)))) / x1) - 3.0d0) / x1))
if (x1 <= (-1.5d+19)) then
tmp = t_0
else if (x1 <= 580.0d0) then
tmp = (x2 * (-6.0d0)) + (x1 * (((2.0d0 * x2) * ((-6.0d0) + (x2 * 4.0d0))) + ((-1.0d0) + (x1 * ((-6.0d0) + ((6.0d0 + (x2 * (-4.0d0))) + ((3.0d0 * (3.0d0 + (2.0d0 * x2))) + (x2 * 10.0d0))))))))
else
tmp = t_0
end if
code = tmp
end function
public static double code(double x1, double x2) {
double t_0 = ((x1 * x1) * (x1 * x1)) * (6.0 + ((((15.0 - (18.0 + (x2 * -8.0))) / x1) - 3.0) / x1));
double tmp;
if (x1 <= -1.5e+19) {
tmp = t_0;
} else if (x1 <= 580.0) {
tmp = (x2 * -6.0) + (x1 * (((2.0 * x2) * (-6.0 + (x2 * 4.0))) + (-1.0 + (x1 * (-6.0 + ((6.0 + (x2 * -4.0)) + ((3.0 * (3.0 + (2.0 * x2))) + (x2 * 10.0))))))));
} else {
tmp = t_0;
}
return tmp;
}
def code(x1, x2): t_0 = ((x1 * x1) * (x1 * x1)) * (6.0 + ((((15.0 - (18.0 + (x2 * -8.0))) / x1) - 3.0) / x1)) tmp = 0 if x1 <= -1.5e+19: tmp = t_0 elif x1 <= 580.0: tmp = (x2 * -6.0) + (x1 * (((2.0 * x2) * (-6.0 + (x2 * 4.0))) + (-1.0 + (x1 * (-6.0 + ((6.0 + (x2 * -4.0)) + ((3.0 * (3.0 + (2.0 * x2))) + (x2 * 10.0)))))))) else: tmp = t_0 return tmp
function code(x1, x2) t_0 = Float64(Float64(Float64(x1 * x1) * Float64(x1 * x1)) * Float64(6.0 + Float64(Float64(Float64(Float64(15.0 - Float64(18.0 + Float64(x2 * -8.0))) / x1) - 3.0) / x1))) tmp = 0.0 if (x1 <= -1.5e+19) tmp = t_0; elseif (x1 <= 580.0) tmp = Float64(Float64(x2 * -6.0) + Float64(x1 * Float64(Float64(Float64(2.0 * x2) * Float64(-6.0 + Float64(x2 * 4.0))) + Float64(-1.0 + Float64(x1 * Float64(-6.0 + Float64(Float64(6.0 + Float64(x2 * -4.0)) + Float64(Float64(3.0 * Float64(3.0 + Float64(2.0 * x2))) + Float64(x2 * 10.0))))))))); else tmp = t_0; end return tmp end
function tmp_2 = code(x1, x2) t_0 = ((x1 * x1) * (x1 * x1)) * (6.0 + ((((15.0 - (18.0 + (x2 * -8.0))) / x1) - 3.0) / x1)); tmp = 0.0; if (x1 <= -1.5e+19) tmp = t_0; elseif (x1 <= 580.0) tmp = (x2 * -6.0) + (x1 * (((2.0 * x2) * (-6.0 + (x2 * 4.0))) + (-1.0 + (x1 * (-6.0 + ((6.0 + (x2 * -4.0)) + ((3.0 * (3.0 + (2.0 * x2))) + (x2 * 10.0)))))))); else tmp = t_0; end tmp_2 = tmp; end
code[x1_, x2_] := Block[{t$95$0 = N[(N[(N[(x1 * x1), $MachinePrecision] * N[(x1 * x1), $MachinePrecision]), $MachinePrecision] * N[(6.0 + N[(N[(N[(N[(15.0 - N[(18.0 + N[(x2 * -8.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / x1), $MachinePrecision] - 3.0), $MachinePrecision] / x1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[x1, -1.5e+19], t$95$0, If[LessEqual[x1, 580.0], N[(N[(x2 * -6.0), $MachinePrecision] + N[(x1 * N[(N[(N[(2.0 * x2), $MachinePrecision] * N[(-6.0 + N[(x2 * 4.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(-1.0 + N[(x1 * N[(-6.0 + N[(N[(6.0 + N[(x2 * -4.0), $MachinePrecision]), $MachinePrecision] + N[(N[(3.0 * N[(3.0 + N[(2.0 * x2), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(x2 * 10.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$0]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \left(\left(x1 \cdot x1\right) \cdot \left(x1 \cdot x1\right)\right) \cdot \left(6 + \frac{\frac{15 - \left(18 + x2 \cdot -8\right)}{x1} - 3}{x1}\right)\\
\mathbf{if}\;x1 \leq -1.5 \cdot 10^{+19}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;x1 \leq 580:\\
\;\;\;\;x2 \cdot -6 + x1 \cdot \left(\left(2 \cdot x2\right) \cdot \left(-6 + x2 \cdot 4\right) + \left(-1 + x1 \cdot \left(-6 + \left(\left(6 + x2 \cdot -4\right) + \left(3 \cdot \left(3 + 2 \cdot x2\right) + x2 \cdot 10\right)\right)\right)\right)\right)\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
if x1 < -1.5e19 or 580 < x1 Initial program 35.6%
Simplified35.6%
Applied egg-rr35.6%
Taylor expanded in x1 around -inf
*-lowering-*.f64N/A
metadata-evalN/A
pow-sqrN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
/-lowering-/.f64N/A
Simplified97.0%
if -1.5e19 < x1 < 580Initial program 99.5%
Simplified99.5%
Applied egg-rr99.5%
Taylor expanded in x1 around 0
Simplified86.9%
Final simplification91.2%
(FPCore (x1 x2)
:precision binary64
(let* ((t_0
(*
(* (* x1 x1) (* x1 x1))
(+ 6.0 (/ (- (/ (- 15.0 (+ 18.0 (* x2 -8.0))) x1) 3.0) x1)))))
(if (<= x1 -1.5e+19)
t_0
(if (<= x1 15000.0)
(+
(/ (* 3.0 (- (* x2 -2.0) (* x1 (- 1.0 (* x1 3.0))))) (+ (* x1 x1) 1.0))
(* x1 (+ 2.0 (* (* 2.0 x2) (+ -6.0 (* x2 4.0))))))
t_0))))
double code(double x1, double x2) {
double t_0 = ((x1 * x1) * (x1 * x1)) * (6.0 + ((((15.0 - (18.0 + (x2 * -8.0))) / x1) - 3.0) / x1));
double tmp;
if (x1 <= -1.5e+19) {
tmp = t_0;
} else if (x1 <= 15000.0) {
tmp = ((3.0 * ((x2 * -2.0) - (x1 * (1.0 - (x1 * 3.0))))) / ((x1 * x1) + 1.0)) + (x1 * (2.0 + ((2.0 * x2) * (-6.0 + (x2 * 4.0)))));
} else {
tmp = t_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) :: tmp
t_0 = ((x1 * x1) * (x1 * x1)) * (6.0d0 + ((((15.0d0 - (18.0d0 + (x2 * (-8.0d0)))) / x1) - 3.0d0) / x1))
if (x1 <= (-1.5d+19)) then
tmp = t_0
else if (x1 <= 15000.0d0) then
tmp = ((3.0d0 * ((x2 * (-2.0d0)) - (x1 * (1.0d0 - (x1 * 3.0d0))))) / ((x1 * x1) + 1.0d0)) + (x1 * (2.0d0 + ((2.0d0 * x2) * ((-6.0d0) + (x2 * 4.0d0)))))
else
tmp = t_0
end if
code = tmp
end function
public static double code(double x1, double x2) {
double t_0 = ((x1 * x1) * (x1 * x1)) * (6.0 + ((((15.0 - (18.0 + (x2 * -8.0))) / x1) - 3.0) / x1));
double tmp;
if (x1 <= -1.5e+19) {
tmp = t_0;
} else if (x1 <= 15000.0) {
tmp = ((3.0 * ((x2 * -2.0) - (x1 * (1.0 - (x1 * 3.0))))) / ((x1 * x1) + 1.0)) + (x1 * (2.0 + ((2.0 * x2) * (-6.0 + (x2 * 4.0)))));
} else {
tmp = t_0;
}
return tmp;
}
def code(x1, x2): t_0 = ((x1 * x1) * (x1 * x1)) * (6.0 + ((((15.0 - (18.0 + (x2 * -8.0))) / x1) - 3.0) / x1)) tmp = 0 if x1 <= -1.5e+19: tmp = t_0 elif x1 <= 15000.0: tmp = ((3.0 * ((x2 * -2.0) - (x1 * (1.0 - (x1 * 3.0))))) / ((x1 * x1) + 1.0)) + (x1 * (2.0 + ((2.0 * x2) * (-6.0 + (x2 * 4.0))))) else: tmp = t_0 return tmp
function code(x1, x2) t_0 = Float64(Float64(Float64(x1 * x1) * Float64(x1 * x1)) * Float64(6.0 + Float64(Float64(Float64(Float64(15.0 - Float64(18.0 + Float64(x2 * -8.0))) / x1) - 3.0) / x1))) tmp = 0.0 if (x1 <= -1.5e+19) tmp = t_0; elseif (x1 <= 15000.0) tmp = Float64(Float64(Float64(3.0 * Float64(Float64(x2 * -2.0) - Float64(x1 * Float64(1.0 - Float64(x1 * 3.0))))) / Float64(Float64(x1 * x1) + 1.0)) + Float64(x1 * Float64(2.0 + Float64(Float64(2.0 * x2) * Float64(-6.0 + Float64(x2 * 4.0)))))); else tmp = t_0; end return tmp end
function tmp_2 = code(x1, x2) t_0 = ((x1 * x1) * (x1 * x1)) * (6.0 + ((((15.0 - (18.0 + (x2 * -8.0))) / x1) - 3.0) / x1)); tmp = 0.0; if (x1 <= -1.5e+19) tmp = t_0; elseif (x1 <= 15000.0) tmp = ((3.0 * ((x2 * -2.0) - (x1 * (1.0 - (x1 * 3.0))))) / ((x1 * x1) + 1.0)) + (x1 * (2.0 + ((2.0 * x2) * (-6.0 + (x2 * 4.0))))); else tmp = t_0; end tmp_2 = tmp; end
code[x1_, x2_] := Block[{t$95$0 = N[(N[(N[(x1 * x1), $MachinePrecision] * N[(x1 * x1), $MachinePrecision]), $MachinePrecision] * N[(6.0 + N[(N[(N[(N[(15.0 - N[(18.0 + N[(x2 * -8.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / x1), $MachinePrecision] - 3.0), $MachinePrecision] / x1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[x1, -1.5e+19], t$95$0, If[LessEqual[x1, 15000.0], N[(N[(N[(3.0 * N[(N[(x2 * -2.0), $MachinePrecision] - N[(x1 * N[(1.0 - N[(x1 * 3.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(N[(x1 * x1), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision] + N[(x1 * N[(2.0 + N[(N[(2.0 * x2), $MachinePrecision] * N[(-6.0 + N[(x2 * 4.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$0]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \left(\left(x1 \cdot x1\right) \cdot \left(x1 \cdot x1\right)\right) \cdot \left(6 + \frac{\frac{15 - \left(18 + x2 \cdot -8\right)}{x1} - 3}{x1}\right)\\
\mathbf{if}\;x1 \leq -1.5 \cdot 10^{+19}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;x1 \leq 15000:\\
\;\;\;\;\frac{3 \cdot \left(x2 \cdot -2 - x1 \cdot \left(1 - x1 \cdot 3\right)\right)}{x1 \cdot x1 + 1} + x1 \cdot \left(2 + \left(2 \cdot x2\right) \cdot \left(-6 + x2 \cdot 4\right)\right)\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
if x1 < -1.5e19 or 15000 < x1 Initial program 35.6%
Simplified35.6%
Applied egg-rr35.6%
Taylor expanded in x1 around -inf
*-lowering-*.f64N/A
metadata-evalN/A
pow-sqrN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
/-lowering-/.f64N/A
Simplified97.0%
if -1.5e19 < x1 < 15000Initial program 99.5%
Simplified99.5%
Taylor expanded in x1 around 0
*-lowering-*.f64N/A
+-lowering-+.f64N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
*-lowering-*.f6486.7%
Simplified86.7%
Final simplification91.1%
(FPCore (x1 x2)
:precision binary64
(let* ((t_0
(*
(* (* x1 x1) (* x1 x1))
(+ 6.0 (/ (- (/ (- 15.0 (+ 18.0 (* x2 -8.0))) x1) 3.0) x1)))))
(if (<= x1 -1.55e+19)
t_0
(if (<= x1 95.0)
(+ (* x2 -6.0) (* x1 (+ -1.0 (* (* 2.0 x2) (+ -6.0 (* x2 4.0))))))
t_0))))
double code(double x1, double x2) {
double t_0 = ((x1 * x1) * (x1 * x1)) * (6.0 + ((((15.0 - (18.0 + (x2 * -8.0))) / x1) - 3.0) / x1));
double tmp;
if (x1 <= -1.55e+19) {
tmp = t_0;
} else if (x1 <= 95.0) {
tmp = (x2 * -6.0) + (x1 * (-1.0 + ((2.0 * x2) * (-6.0 + (x2 * 4.0)))));
} else {
tmp = t_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) :: tmp
t_0 = ((x1 * x1) * (x1 * x1)) * (6.0d0 + ((((15.0d0 - (18.0d0 + (x2 * (-8.0d0)))) / x1) - 3.0d0) / x1))
if (x1 <= (-1.55d+19)) then
tmp = t_0
else if (x1 <= 95.0d0) then
tmp = (x2 * (-6.0d0)) + (x1 * ((-1.0d0) + ((2.0d0 * x2) * ((-6.0d0) + (x2 * 4.0d0)))))
else
tmp = t_0
end if
code = tmp
end function
public static double code(double x1, double x2) {
double t_0 = ((x1 * x1) * (x1 * x1)) * (6.0 + ((((15.0 - (18.0 + (x2 * -8.0))) / x1) - 3.0) / x1));
double tmp;
if (x1 <= -1.55e+19) {
tmp = t_0;
} else if (x1 <= 95.0) {
tmp = (x2 * -6.0) + (x1 * (-1.0 + ((2.0 * x2) * (-6.0 + (x2 * 4.0)))));
} else {
tmp = t_0;
}
return tmp;
}
def code(x1, x2): t_0 = ((x1 * x1) * (x1 * x1)) * (6.0 + ((((15.0 - (18.0 + (x2 * -8.0))) / x1) - 3.0) / x1)) tmp = 0 if x1 <= -1.55e+19: tmp = t_0 elif x1 <= 95.0: tmp = (x2 * -6.0) + (x1 * (-1.0 + ((2.0 * x2) * (-6.0 + (x2 * 4.0))))) else: tmp = t_0 return tmp
function code(x1, x2) t_0 = Float64(Float64(Float64(x1 * x1) * Float64(x1 * x1)) * Float64(6.0 + Float64(Float64(Float64(Float64(15.0 - Float64(18.0 + Float64(x2 * -8.0))) / x1) - 3.0) / x1))) tmp = 0.0 if (x1 <= -1.55e+19) tmp = t_0; elseif (x1 <= 95.0) tmp = Float64(Float64(x2 * -6.0) + Float64(x1 * Float64(-1.0 + Float64(Float64(2.0 * x2) * Float64(-6.0 + Float64(x2 * 4.0)))))); else tmp = t_0; end return tmp end
function tmp_2 = code(x1, x2) t_0 = ((x1 * x1) * (x1 * x1)) * (6.0 + ((((15.0 - (18.0 + (x2 * -8.0))) / x1) - 3.0) / x1)); tmp = 0.0; if (x1 <= -1.55e+19) tmp = t_0; elseif (x1 <= 95.0) tmp = (x2 * -6.0) + (x1 * (-1.0 + ((2.0 * x2) * (-6.0 + (x2 * 4.0))))); else tmp = t_0; end tmp_2 = tmp; end
code[x1_, x2_] := Block[{t$95$0 = N[(N[(N[(x1 * x1), $MachinePrecision] * N[(x1 * x1), $MachinePrecision]), $MachinePrecision] * N[(6.0 + N[(N[(N[(N[(15.0 - N[(18.0 + N[(x2 * -8.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / x1), $MachinePrecision] - 3.0), $MachinePrecision] / x1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[x1, -1.55e+19], t$95$0, If[LessEqual[x1, 95.0], N[(N[(x2 * -6.0), $MachinePrecision] + N[(x1 * N[(-1.0 + N[(N[(2.0 * x2), $MachinePrecision] * N[(-6.0 + N[(x2 * 4.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$0]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \left(\left(x1 \cdot x1\right) \cdot \left(x1 \cdot x1\right)\right) \cdot \left(6 + \frac{\frac{15 - \left(18 + x2 \cdot -8\right)}{x1} - 3}{x1}\right)\\
\mathbf{if}\;x1 \leq -1.55 \cdot 10^{+19}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;x1 \leq 95:\\
\;\;\;\;x2 \cdot -6 + x1 \cdot \left(-1 + \left(2 \cdot x2\right) \cdot \left(-6 + x2 \cdot 4\right)\right)\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
if x1 < -1.55e19 or 95 < x1 Initial program 35.6%
Simplified35.6%
Applied egg-rr35.6%
Taylor expanded in x1 around -inf
*-lowering-*.f64N/A
metadata-evalN/A
pow-sqrN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
/-lowering-/.f64N/A
Simplified97.0%
if -1.55e19 < x1 < 95Initial program 99.5%
Simplified99.5%
Applied egg-rr99.5%
Taylor expanded in x1 around 0
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6485.9%
Simplified85.9%
Final simplification90.6%
(FPCore (x1 x2)
:precision binary64
(if (<= x1 -1.5e+19)
(+ x1 (* 6.0 (* x1 (* x1 (* x1 x1)))))
(if (<= x1 -7.6e-141)
(* (* x2 x2) (- (/ (* x1 8.0) (+ (* x1 x1) 1.0)) (/ 6.0 x2)))
(if (<= x1 2.1)
(+ (* x2 -6.0) (* x1 (+ -1.0 (* (* x1 3.0) (+ 3.0 (* 2.0 x2))))))
(* (* (* x1 x1) (* x1 x1)) (+ 6.0 (/ -3.0 x1)))))))
double code(double x1, double x2) {
double tmp;
if (x1 <= -1.5e+19) {
tmp = x1 + (6.0 * (x1 * (x1 * (x1 * x1))));
} else if (x1 <= -7.6e-141) {
tmp = (x2 * x2) * (((x1 * 8.0) / ((x1 * x1) + 1.0)) - (6.0 / x2));
} else if (x1 <= 2.1) {
tmp = (x2 * -6.0) + (x1 * (-1.0 + ((x1 * 3.0) * (3.0 + (2.0 * x2)))));
} else {
tmp = ((x1 * x1) * (x1 * x1)) * (6.0 + (-3.0 / x1));
}
return tmp;
}
real(8) function code(x1, x2)
real(8), intent (in) :: x1
real(8), intent (in) :: x2
real(8) :: tmp
if (x1 <= (-1.5d+19)) then
tmp = x1 + (6.0d0 * (x1 * (x1 * (x1 * x1))))
else if (x1 <= (-7.6d-141)) then
tmp = (x2 * x2) * (((x1 * 8.0d0) / ((x1 * x1) + 1.0d0)) - (6.0d0 / x2))
else if (x1 <= 2.1d0) then
tmp = (x2 * (-6.0d0)) + (x1 * ((-1.0d0) + ((x1 * 3.0d0) * (3.0d0 + (2.0d0 * x2)))))
else
tmp = ((x1 * x1) * (x1 * x1)) * (6.0d0 + ((-3.0d0) / x1))
end if
code = tmp
end function
public static double code(double x1, double x2) {
double tmp;
if (x1 <= -1.5e+19) {
tmp = x1 + (6.0 * (x1 * (x1 * (x1 * x1))));
} else if (x1 <= -7.6e-141) {
tmp = (x2 * x2) * (((x1 * 8.0) / ((x1 * x1) + 1.0)) - (6.0 / x2));
} else if (x1 <= 2.1) {
tmp = (x2 * -6.0) + (x1 * (-1.0 + ((x1 * 3.0) * (3.0 + (2.0 * x2)))));
} else {
tmp = ((x1 * x1) * (x1 * x1)) * (6.0 + (-3.0 / x1));
}
return tmp;
}
def code(x1, x2): tmp = 0 if x1 <= -1.5e+19: tmp = x1 + (6.0 * (x1 * (x1 * (x1 * x1)))) elif x1 <= -7.6e-141: tmp = (x2 * x2) * (((x1 * 8.0) / ((x1 * x1) + 1.0)) - (6.0 / x2)) elif x1 <= 2.1: tmp = (x2 * -6.0) + (x1 * (-1.0 + ((x1 * 3.0) * (3.0 + (2.0 * x2))))) else: tmp = ((x1 * x1) * (x1 * x1)) * (6.0 + (-3.0 / x1)) return tmp
function code(x1, x2) tmp = 0.0 if (x1 <= -1.5e+19) tmp = Float64(x1 + Float64(6.0 * Float64(x1 * Float64(x1 * Float64(x1 * x1))))); elseif (x1 <= -7.6e-141) tmp = Float64(Float64(x2 * x2) * Float64(Float64(Float64(x1 * 8.0) / Float64(Float64(x1 * x1) + 1.0)) - Float64(6.0 / x2))); elseif (x1 <= 2.1) tmp = Float64(Float64(x2 * -6.0) + Float64(x1 * Float64(-1.0 + Float64(Float64(x1 * 3.0) * Float64(3.0 + Float64(2.0 * x2)))))); else tmp = Float64(Float64(Float64(x1 * x1) * Float64(x1 * x1)) * Float64(6.0 + Float64(-3.0 / x1))); end return tmp end
function tmp_2 = code(x1, x2) tmp = 0.0; if (x1 <= -1.5e+19) tmp = x1 + (6.0 * (x1 * (x1 * (x1 * x1)))); elseif (x1 <= -7.6e-141) tmp = (x2 * x2) * (((x1 * 8.0) / ((x1 * x1) + 1.0)) - (6.0 / x2)); elseif (x1 <= 2.1) tmp = (x2 * -6.0) + (x1 * (-1.0 + ((x1 * 3.0) * (3.0 + (2.0 * x2))))); else tmp = ((x1 * x1) * (x1 * x1)) * (6.0 + (-3.0 / x1)); end tmp_2 = tmp; end
code[x1_, x2_] := If[LessEqual[x1, -1.5e+19], N[(x1 + N[(6.0 * N[(x1 * N[(x1 * N[(x1 * x1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[x1, -7.6e-141], N[(N[(x2 * x2), $MachinePrecision] * N[(N[(N[(x1 * 8.0), $MachinePrecision] / N[(N[(x1 * x1), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision] - N[(6.0 / x2), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[x1, 2.1], N[(N[(x2 * -6.0), $MachinePrecision] + N[(x1 * N[(-1.0 + N[(N[(x1 * 3.0), $MachinePrecision] * N[(3.0 + N[(2.0 * x2), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(N[(x1 * x1), $MachinePrecision] * N[(x1 * x1), $MachinePrecision]), $MachinePrecision] * N[(6.0 + N[(-3.0 / x1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x1 \leq -1.5 \cdot 10^{+19}:\\
\;\;\;\;x1 + 6 \cdot \left(x1 \cdot \left(x1 \cdot \left(x1 \cdot x1\right)\right)\right)\\
\mathbf{elif}\;x1 \leq -7.6 \cdot 10^{-141}:\\
\;\;\;\;\left(x2 \cdot x2\right) \cdot \left(\frac{x1 \cdot 8}{x1 \cdot x1 + 1} - \frac{6}{x2}\right)\\
\mathbf{elif}\;x1 \leq 2.1:\\
\;\;\;\;x2 \cdot -6 + x1 \cdot \left(-1 + \left(x1 \cdot 3\right) \cdot \left(3 + 2 \cdot x2\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\left(\left(x1 \cdot x1\right) \cdot \left(x1 \cdot x1\right)\right) \cdot \left(6 + \frac{-3}{x1}\right)\\
\end{array}
\end{array}
if x1 < -1.5e19Initial program 28.4%
Taylor expanded in x1 around inf
*-lowering-*.f64N/A
pow-lowering-pow.f6493.0%
Simplified93.0%
metadata-evalN/A
pow-prod-upN/A
pow2N/A
pow2N/A
*-commutativeN/A
*-lowering-*.f64N/A
associate-*l*N/A
cube-multN/A
*-lowering-*.f64N/A
cube-multN/A
*-lowering-*.f64N/A
*-lowering-*.f6493.0%
Applied egg-rr93.0%
if -1.5e19 < x1 < -7.59999999999999973e-141Initial program 99.4%
Simplified99.3%
Applied egg-rr99.3%
Taylor expanded in x2 around -inf
Simplified69.8%
Taylor expanded in x1 around 0
Simplified69.8%
if -7.59999999999999973e-141 < x1 < 2.10000000000000009Initial program 99.5%
Simplified99.6%
Taylor expanded in x1 around inf
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6478.9%
Simplified78.9%
Taylor expanded in x1 around 0
*-commutativeN/A
*-lowering-*.f6479.3%
Simplified79.3%
Taylor expanded in x1 around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
cancel-sign-sub-invN/A
metadata-evalN/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6479.5%
Simplified79.5%
if 2.10000000000000009 < x1 Initial program 44.3%
Simplified44.3%
Applied egg-rr44.2%
Taylor expanded in x1 around inf
*-lowering-*.f64N/A
metadata-evalN/A
pow-sqrN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
associate-*r/N/A
metadata-evalN/A
distribute-neg-fracN/A
metadata-evalN/A
/-lowering-/.f6490.9%
Simplified90.9%
Final simplification83.4%
(FPCore (x1 x2)
:precision binary64
(if (<= x1 -1.55e+22)
(+ x1 (* 6.0 (* x1 (* x1 (* x1 x1)))))
(if (<= x1 -2.8e-147)
(* (* x2 x2) (- (/ (* x1 8.0) (+ (* x1 x1) 1.0)) (/ 6.0 x2)))
(if (<= x1 6.5)
(- (* x2 -6.0) x1)
(* (* (* x1 x1) (* x1 x1)) (+ 6.0 (/ -3.0 x1)))))))
double code(double x1, double x2) {
double tmp;
if (x1 <= -1.55e+22) {
tmp = x1 + (6.0 * (x1 * (x1 * (x1 * x1))));
} else if (x1 <= -2.8e-147) {
tmp = (x2 * x2) * (((x1 * 8.0) / ((x1 * x1) + 1.0)) - (6.0 / x2));
} else if (x1 <= 6.5) {
tmp = (x2 * -6.0) - x1;
} else {
tmp = ((x1 * x1) * (x1 * x1)) * (6.0 + (-3.0 / x1));
}
return tmp;
}
real(8) function code(x1, x2)
real(8), intent (in) :: x1
real(8), intent (in) :: x2
real(8) :: tmp
if (x1 <= (-1.55d+22)) then
tmp = x1 + (6.0d0 * (x1 * (x1 * (x1 * x1))))
else if (x1 <= (-2.8d-147)) then
tmp = (x2 * x2) * (((x1 * 8.0d0) / ((x1 * x1) + 1.0d0)) - (6.0d0 / x2))
else if (x1 <= 6.5d0) then
tmp = (x2 * (-6.0d0)) - x1
else
tmp = ((x1 * x1) * (x1 * x1)) * (6.0d0 + ((-3.0d0) / x1))
end if
code = tmp
end function
public static double code(double x1, double x2) {
double tmp;
if (x1 <= -1.55e+22) {
tmp = x1 + (6.0 * (x1 * (x1 * (x1 * x1))));
} else if (x1 <= -2.8e-147) {
tmp = (x2 * x2) * (((x1 * 8.0) / ((x1 * x1) + 1.0)) - (6.0 / x2));
} else if (x1 <= 6.5) {
tmp = (x2 * -6.0) - x1;
} else {
tmp = ((x1 * x1) * (x1 * x1)) * (6.0 + (-3.0 / x1));
}
return tmp;
}
def code(x1, x2): tmp = 0 if x1 <= -1.55e+22: tmp = x1 + (6.0 * (x1 * (x1 * (x1 * x1)))) elif x1 <= -2.8e-147: tmp = (x2 * x2) * (((x1 * 8.0) / ((x1 * x1) + 1.0)) - (6.0 / x2)) elif x1 <= 6.5: tmp = (x2 * -6.0) - x1 else: tmp = ((x1 * x1) * (x1 * x1)) * (6.0 + (-3.0 / x1)) return tmp
function code(x1, x2) tmp = 0.0 if (x1 <= -1.55e+22) tmp = Float64(x1 + Float64(6.0 * Float64(x1 * Float64(x1 * Float64(x1 * x1))))); elseif (x1 <= -2.8e-147) tmp = Float64(Float64(x2 * x2) * Float64(Float64(Float64(x1 * 8.0) / Float64(Float64(x1 * x1) + 1.0)) - Float64(6.0 / x2))); elseif (x1 <= 6.5) tmp = Float64(Float64(x2 * -6.0) - x1); else tmp = Float64(Float64(Float64(x1 * x1) * Float64(x1 * x1)) * Float64(6.0 + Float64(-3.0 / x1))); end return tmp end
function tmp_2 = code(x1, x2) tmp = 0.0; if (x1 <= -1.55e+22) tmp = x1 + (6.0 * (x1 * (x1 * (x1 * x1)))); elseif (x1 <= -2.8e-147) tmp = (x2 * x2) * (((x1 * 8.0) / ((x1 * x1) + 1.0)) - (6.0 / x2)); elseif (x1 <= 6.5) tmp = (x2 * -6.0) - x1; else tmp = ((x1 * x1) * (x1 * x1)) * (6.0 + (-3.0 / x1)); end tmp_2 = tmp; end
code[x1_, x2_] := If[LessEqual[x1, -1.55e+22], N[(x1 + N[(6.0 * N[(x1 * N[(x1 * N[(x1 * x1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[x1, -2.8e-147], N[(N[(x2 * x2), $MachinePrecision] * N[(N[(N[(x1 * 8.0), $MachinePrecision] / N[(N[(x1 * x1), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision] - N[(6.0 / x2), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[x1, 6.5], N[(N[(x2 * -6.0), $MachinePrecision] - x1), $MachinePrecision], N[(N[(N[(x1 * x1), $MachinePrecision] * N[(x1 * x1), $MachinePrecision]), $MachinePrecision] * N[(6.0 + N[(-3.0 / x1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x1 \leq -1.55 \cdot 10^{+22}:\\
\;\;\;\;x1 + 6 \cdot \left(x1 \cdot \left(x1 \cdot \left(x1 \cdot x1\right)\right)\right)\\
\mathbf{elif}\;x1 \leq -2.8 \cdot 10^{-147}:\\
\;\;\;\;\left(x2 \cdot x2\right) \cdot \left(\frac{x1 \cdot 8}{x1 \cdot x1 + 1} - \frac{6}{x2}\right)\\
\mathbf{elif}\;x1 \leq 6.5:\\
\;\;\;\;x2 \cdot -6 - x1\\
\mathbf{else}:\\
\;\;\;\;\left(\left(x1 \cdot x1\right) \cdot \left(x1 \cdot x1\right)\right) \cdot \left(6 + \frac{-3}{x1}\right)\\
\end{array}
\end{array}
if x1 < -1.5500000000000001e22Initial program 28.4%
Taylor expanded in x1 around inf
*-lowering-*.f64N/A
pow-lowering-pow.f6493.0%
Simplified93.0%
metadata-evalN/A
pow-prod-upN/A
pow2N/A
pow2N/A
*-commutativeN/A
*-lowering-*.f64N/A
associate-*l*N/A
cube-multN/A
*-lowering-*.f64N/A
cube-multN/A
*-lowering-*.f64N/A
*-lowering-*.f6493.0%
Applied egg-rr93.0%
if -1.5500000000000001e22 < x1 < -2.8e-147Initial program 99.4%
Simplified99.3%
Applied egg-rr99.3%
Taylor expanded in x2 around -inf
Simplified69.8%
Taylor expanded in x1 around 0
Simplified69.8%
if -2.8e-147 < x1 < 6.5Initial program 99.5%
Simplified99.6%
Taylor expanded in x1 around inf
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6478.2%
Simplified78.2%
Taylor expanded in x1 around 0
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
*-commutativeN/A
*-lowering-*.f6478.4%
Simplified78.4%
if 6.5 < x1 Initial program 43.2%
Simplified43.2%
Applied egg-rr43.2%
Taylor expanded in x1 around inf
*-lowering-*.f64N/A
metadata-evalN/A
pow-sqrN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
associate-*r/N/A
metadata-evalN/A
distribute-neg-fracN/A
metadata-evalN/A
/-lowering-/.f6492.5%
Simplified92.5%
Final simplification83.2%
(FPCore (x1 x2)
:precision binary64
(if (<= x1 -1.5e+19)
(+ x1 (* 6.0 (* x1 (* x1 (* x1 x1)))))
(if (<= x1 -5e-145)
(+ (* x2 -6.0) (* x1 (+ -3.0 (* 8.0 (* x2 x2)))))
(if (<= x1 6.5)
(- (* x2 -6.0) x1)
(* (* (* x1 x1) (* x1 x1)) (+ 6.0 (/ -3.0 x1)))))))
double code(double x1, double x2) {
double tmp;
if (x1 <= -1.5e+19) {
tmp = x1 + (6.0 * (x1 * (x1 * (x1 * x1))));
} else if (x1 <= -5e-145) {
tmp = (x2 * -6.0) + (x1 * (-3.0 + (8.0 * (x2 * x2))));
} else if (x1 <= 6.5) {
tmp = (x2 * -6.0) - x1;
} else {
tmp = ((x1 * x1) * (x1 * x1)) * (6.0 + (-3.0 / x1));
}
return tmp;
}
real(8) function code(x1, x2)
real(8), intent (in) :: x1
real(8), intent (in) :: x2
real(8) :: tmp
if (x1 <= (-1.5d+19)) then
tmp = x1 + (6.0d0 * (x1 * (x1 * (x1 * x1))))
else if (x1 <= (-5d-145)) then
tmp = (x2 * (-6.0d0)) + (x1 * ((-3.0d0) + (8.0d0 * (x2 * x2))))
else if (x1 <= 6.5d0) then
tmp = (x2 * (-6.0d0)) - x1
else
tmp = ((x1 * x1) * (x1 * x1)) * (6.0d0 + ((-3.0d0) / x1))
end if
code = tmp
end function
public static double code(double x1, double x2) {
double tmp;
if (x1 <= -1.5e+19) {
tmp = x1 + (6.0 * (x1 * (x1 * (x1 * x1))));
} else if (x1 <= -5e-145) {
tmp = (x2 * -6.0) + (x1 * (-3.0 + (8.0 * (x2 * x2))));
} else if (x1 <= 6.5) {
tmp = (x2 * -6.0) - x1;
} else {
tmp = ((x1 * x1) * (x1 * x1)) * (6.0 + (-3.0 / x1));
}
return tmp;
}
def code(x1, x2): tmp = 0 if x1 <= -1.5e+19: tmp = x1 + (6.0 * (x1 * (x1 * (x1 * x1)))) elif x1 <= -5e-145: tmp = (x2 * -6.0) + (x1 * (-3.0 + (8.0 * (x2 * x2)))) elif x1 <= 6.5: tmp = (x2 * -6.0) - x1 else: tmp = ((x1 * x1) * (x1 * x1)) * (6.0 + (-3.0 / x1)) return tmp
function code(x1, x2) tmp = 0.0 if (x1 <= -1.5e+19) tmp = Float64(x1 + Float64(6.0 * Float64(x1 * Float64(x1 * Float64(x1 * x1))))); elseif (x1 <= -5e-145) tmp = Float64(Float64(x2 * -6.0) + Float64(x1 * Float64(-3.0 + Float64(8.0 * Float64(x2 * x2))))); elseif (x1 <= 6.5) tmp = Float64(Float64(x2 * -6.0) - x1); else tmp = Float64(Float64(Float64(x1 * x1) * Float64(x1 * x1)) * Float64(6.0 + Float64(-3.0 / x1))); end return tmp end
function tmp_2 = code(x1, x2) tmp = 0.0; if (x1 <= -1.5e+19) tmp = x1 + (6.0 * (x1 * (x1 * (x1 * x1)))); elseif (x1 <= -5e-145) tmp = (x2 * -6.0) + (x1 * (-3.0 + (8.0 * (x2 * x2)))); elseif (x1 <= 6.5) tmp = (x2 * -6.0) - x1; else tmp = ((x1 * x1) * (x1 * x1)) * (6.0 + (-3.0 / x1)); end tmp_2 = tmp; end
code[x1_, x2_] := If[LessEqual[x1, -1.5e+19], N[(x1 + N[(6.0 * N[(x1 * N[(x1 * N[(x1 * x1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[x1, -5e-145], N[(N[(x2 * -6.0), $MachinePrecision] + N[(x1 * N[(-3.0 + N[(8.0 * N[(x2 * x2), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[x1, 6.5], N[(N[(x2 * -6.0), $MachinePrecision] - x1), $MachinePrecision], N[(N[(N[(x1 * x1), $MachinePrecision] * N[(x1 * x1), $MachinePrecision]), $MachinePrecision] * N[(6.0 + N[(-3.0 / x1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x1 \leq -1.5 \cdot 10^{+19}:\\
\;\;\;\;x1 + 6 \cdot \left(x1 \cdot \left(x1 \cdot \left(x1 \cdot x1\right)\right)\right)\\
\mathbf{elif}\;x1 \leq -5 \cdot 10^{-145}:\\
\;\;\;\;x2 \cdot -6 + x1 \cdot \left(-3 + 8 \cdot \left(x2 \cdot x2\right)\right)\\
\mathbf{elif}\;x1 \leq 6.5:\\
\;\;\;\;x2 \cdot -6 - x1\\
\mathbf{else}:\\
\;\;\;\;\left(\left(x1 \cdot x1\right) \cdot \left(x1 \cdot x1\right)\right) \cdot \left(6 + \frac{-3}{x1}\right)\\
\end{array}
\end{array}
if x1 < -1.5e19Initial program 28.4%
Taylor expanded in x1 around inf
*-lowering-*.f64N/A
pow-lowering-pow.f6493.0%
Simplified93.0%
metadata-evalN/A
pow-prod-upN/A
pow2N/A
pow2N/A
*-commutativeN/A
*-lowering-*.f64N/A
associate-*l*N/A
cube-multN/A
*-lowering-*.f64N/A
cube-multN/A
*-lowering-*.f64N/A
*-lowering-*.f6493.0%
Applied egg-rr93.0%
if -1.5e19 < x1 < -4.9999999999999998e-145Initial program 99.4%
Simplified99.3%
Taylor expanded in x1 around 0
Simplified68.8%
Taylor expanded in x2 around inf
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6464.3%
Simplified64.3%
Taylor expanded in x1 around 0
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6469.0%
Simplified69.0%
if -4.9999999999999998e-145 < x1 < 6.5Initial program 99.5%
Simplified99.6%
Taylor expanded in x1 around inf
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6478.2%
Simplified78.2%
Taylor expanded in x1 around 0
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
*-commutativeN/A
*-lowering-*.f6478.4%
Simplified78.4%
if 6.5 < x1 Initial program 43.2%
Simplified43.2%
Applied egg-rr43.2%
Taylor expanded in x1 around inf
*-lowering-*.f64N/A
metadata-evalN/A
pow-sqrN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
associate-*r/N/A
metadata-evalN/A
distribute-neg-fracN/A
metadata-evalN/A
/-lowering-/.f6492.5%
Simplified92.5%
Final simplification83.1%
(FPCore (x1 x2)
:precision binary64
(if (<= x1 -1.5e+19)
(+ x1 (* 6.0 (* x1 (* x1 (* x1 x1)))))
(if (<= x1 1950.0)
(+ (* x2 -6.0) (* x1 (+ -1.0 (* (* 2.0 x2) (+ -6.0 (* x2 4.0))))))
(* (* (* x1 x1) (* x1 x1)) (+ 6.0 (/ -3.0 x1))))))
double code(double x1, double x2) {
double tmp;
if (x1 <= -1.5e+19) {
tmp = x1 + (6.0 * (x1 * (x1 * (x1 * x1))));
} else if (x1 <= 1950.0) {
tmp = (x2 * -6.0) + (x1 * (-1.0 + ((2.0 * x2) * (-6.0 + (x2 * 4.0)))));
} else {
tmp = ((x1 * x1) * (x1 * x1)) * (6.0 + (-3.0 / x1));
}
return tmp;
}
real(8) function code(x1, x2)
real(8), intent (in) :: x1
real(8), intent (in) :: x2
real(8) :: tmp
if (x1 <= (-1.5d+19)) then
tmp = x1 + (6.0d0 * (x1 * (x1 * (x1 * x1))))
else if (x1 <= 1950.0d0) then
tmp = (x2 * (-6.0d0)) + (x1 * ((-1.0d0) + ((2.0d0 * x2) * ((-6.0d0) + (x2 * 4.0d0)))))
else
tmp = ((x1 * x1) * (x1 * x1)) * (6.0d0 + ((-3.0d0) / x1))
end if
code = tmp
end function
public static double code(double x1, double x2) {
double tmp;
if (x1 <= -1.5e+19) {
tmp = x1 + (6.0 * (x1 * (x1 * (x1 * x1))));
} else if (x1 <= 1950.0) {
tmp = (x2 * -6.0) + (x1 * (-1.0 + ((2.0 * x2) * (-6.0 + (x2 * 4.0)))));
} else {
tmp = ((x1 * x1) * (x1 * x1)) * (6.0 + (-3.0 / x1));
}
return tmp;
}
def code(x1, x2): tmp = 0 if x1 <= -1.5e+19: tmp = x1 + (6.0 * (x1 * (x1 * (x1 * x1)))) elif x1 <= 1950.0: tmp = (x2 * -6.0) + (x1 * (-1.0 + ((2.0 * x2) * (-6.0 + (x2 * 4.0))))) else: tmp = ((x1 * x1) * (x1 * x1)) * (6.0 + (-3.0 / x1)) return tmp
function code(x1, x2) tmp = 0.0 if (x1 <= -1.5e+19) tmp = Float64(x1 + Float64(6.0 * Float64(x1 * Float64(x1 * Float64(x1 * x1))))); elseif (x1 <= 1950.0) tmp = Float64(Float64(x2 * -6.0) + Float64(x1 * Float64(-1.0 + Float64(Float64(2.0 * x2) * Float64(-6.0 + Float64(x2 * 4.0)))))); else tmp = Float64(Float64(Float64(x1 * x1) * Float64(x1 * x1)) * Float64(6.0 + Float64(-3.0 / x1))); end return tmp end
function tmp_2 = code(x1, x2) tmp = 0.0; if (x1 <= -1.5e+19) tmp = x1 + (6.0 * (x1 * (x1 * (x1 * x1)))); elseif (x1 <= 1950.0) tmp = (x2 * -6.0) + (x1 * (-1.0 + ((2.0 * x2) * (-6.0 + (x2 * 4.0))))); else tmp = ((x1 * x1) * (x1 * x1)) * (6.0 + (-3.0 / x1)); end tmp_2 = tmp; end
code[x1_, x2_] := If[LessEqual[x1, -1.5e+19], N[(x1 + N[(6.0 * N[(x1 * N[(x1 * N[(x1 * x1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[x1, 1950.0], N[(N[(x2 * -6.0), $MachinePrecision] + N[(x1 * N[(-1.0 + N[(N[(2.0 * x2), $MachinePrecision] * N[(-6.0 + N[(x2 * 4.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(N[(x1 * x1), $MachinePrecision] * N[(x1 * x1), $MachinePrecision]), $MachinePrecision] * N[(6.0 + N[(-3.0 / x1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x1 \leq -1.5 \cdot 10^{+19}:\\
\;\;\;\;x1 + 6 \cdot \left(x1 \cdot \left(x1 \cdot \left(x1 \cdot x1\right)\right)\right)\\
\mathbf{elif}\;x1 \leq 1950:\\
\;\;\;\;x2 \cdot -6 + x1 \cdot \left(-1 + \left(2 \cdot x2\right) \cdot \left(-6 + x2 \cdot 4\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\left(\left(x1 \cdot x1\right) \cdot \left(x1 \cdot x1\right)\right) \cdot \left(6 + \frac{-3}{x1}\right)\\
\end{array}
\end{array}
if x1 < -1.5e19Initial program 28.4%
Taylor expanded in x1 around inf
*-lowering-*.f64N/A
pow-lowering-pow.f6493.0%
Simplified93.0%
metadata-evalN/A
pow-prod-upN/A
pow2N/A
pow2N/A
*-commutativeN/A
*-lowering-*.f64N/A
associate-*l*N/A
cube-multN/A
*-lowering-*.f64N/A
cube-multN/A
*-lowering-*.f64N/A
*-lowering-*.f6493.0%
Applied egg-rr93.0%
if -1.5e19 < x1 < 1950Initial program 99.5%
Simplified99.5%
Applied egg-rr99.5%
Taylor expanded in x1 around 0
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6485.9%
Simplified85.9%
if 1950 < x1 Initial program 43.2%
Simplified43.2%
Applied egg-rr43.2%
Taylor expanded in x1 around inf
*-lowering-*.f64N/A
metadata-evalN/A
pow-sqrN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
associate-*r/N/A
metadata-evalN/A
distribute-neg-fracN/A
metadata-evalN/A
/-lowering-/.f6492.5%
Simplified92.5%
Final simplification88.8%
(FPCore (x1 x2)
:precision binary64
(if (<= x1 -1.5e+19)
(+ x1 (* 6.0 (* x1 (* x1 (* x1 x1)))))
(if (<= x1 6.5)
(- (* x2 -6.0) x1)
(* (* (* x1 x1) (* x1 x1)) (+ 6.0 (/ -3.0 x1))))))
double code(double x1, double x2) {
double tmp;
if (x1 <= -1.5e+19) {
tmp = x1 + (6.0 * (x1 * (x1 * (x1 * x1))));
} else if (x1 <= 6.5) {
tmp = (x2 * -6.0) - x1;
} else {
tmp = ((x1 * x1) * (x1 * x1)) * (6.0 + (-3.0 / x1));
}
return tmp;
}
real(8) function code(x1, x2)
real(8), intent (in) :: x1
real(8), intent (in) :: x2
real(8) :: tmp
if (x1 <= (-1.5d+19)) then
tmp = x1 + (6.0d0 * (x1 * (x1 * (x1 * x1))))
else if (x1 <= 6.5d0) then
tmp = (x2 * (-6.0d0)) - x1
else
tmp = ((x1 * x1) * (x1 * x1)) * (6.0d0 + ((-3.0d0) / x1))
end if
code = tmp
end function
public static double code(double x1, double x2) {
double tmp;
if (x1 <= -1.5e+19) {
tmp = x1 + (6.0 * (x1 * (x1 * (x1 * x1))));
} else if (x1 <= 6.5) {
tmp = (x2 * -6.0) - x1;
} else {
tmp = ((x1 * x1) * (x1 * x1)) * (6.0 + (-3.0 / x1));
}
return tmp;
}
def code(x1, x2): tmp = 0 if x1 <= -1.5e+19: tmp = x1 + (6.0 * (x1 * (x1 * (x1 * x1)))) elif x1 <= 6.5: tmp = (x2 * -6.0) - x1 else: tmp = ((x1 * x1) * (x1 * x1)) * (6.0 + (-3.0 / x1)) return tmp
function code(x1, x2) tmp = 0.0 if (x1 <= -1.5e+19) tmp = Float64(x1 + Float64(6.0 * Float64(x1 * Float64(x1 * Float64(x1 * x1))))); elseif (x1 <= 6.5) tmp = Float64(Float64(x2 * -6.0) - x1); else tmp = Float64(Float64(Float64(x1 * x1) * Float64(x1 * x1)) * Float64(6.0 + Float64(-3.0 / x1))); end return tmp end
function tmp_2 = code(x1, x2) tmp = 0.0; if (x1 <= -1.5e+19) tmp = x1 + (6.0 * (x1 * (x1 * (x1 * x1)))); elseif (x1 <= 6.5) tmp = (x2 * -6.0) - x1; else tmp = ((x1 * x1) * (x1 * x1)) * (6.0 + (-3.0 / x1)); end tmp_2 = tmp; end
code[x1_, x2_] := If[LessEqual[x1, -1.5e+19], N[(x1 + N[(6.0 * N[(x1 * N[(x1 * N[(x1 * x1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[x1, 6.5], N[(N[(x2 * -6.0), $MachinePrecision] - x1), $MachinePrecision], N[(N[(N[(x1 * x1), $MachinePrecision] * N[(x1 * x1), $MachinePrecision]), $MachinePrecision] * N[(6.0 + N[(-3.0 / x1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x1 \leq -1.5 \cdot 10^{+19}:\\
\;\;\;\;x1 + 6 \cdot \left(x1 \cdot \left(x1 \cdot \left(x1 \cdot x1\right)\right)\right)\\
\mathbf{elif}\;x1 \leq 6.5:\\
\;\;\;\;x2 \cdot -6 - x1\\
\mathbf{else}:\\
\;\;\;\;\left(\left(x1 \cdot x1\right) \cdot \left(x1 \cdot x1\right)\right) \cdot \left(6 + \frac{-3}{x1}\right)\\
\end{array}
\end{array}
if x1 < -1.5e19Initial program 28.4%
Taylor expanded in x1 around inf
*-lowering-*.f64N/A
pow-lowering-pow.f6493.0%
Simplified93.0%
metadata-evalN/A
pow-prod-upN/A
pow2N/A
pow2N/A
*-commutativeN/A
*-lowering-*.f64N/A
associate-*l*N/A
cube-multN/A
*-lowering-*.f64N/A
cube-multN/A
*-lowering-*.f64N/A
*-lowering-*.f6493.0%
Applied egg-rr93.0%
if -1.5e19 < x1 < 6.5Initial program 99.5%
Simplified99.5%
Taylor expanded in x1 around inf
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6470.7%
Simplified70.7%
Taylor expanded in x1 around 0
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
*-commutativeN/A
*-lowering-*.f6471.0%
Simplified71.0%
if 6.5 < x1 Initial program 43.2%
Simplified43.2%
Applied egg-rr43.2%
Taylor expanded in x1 around inf
*-lowering-*.f64N/A
metadata-evalN/A
pow-sqrN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
associate-*r/N/A
metadata-evalN/A
distribute-neg-fracN/A
metadata-evalN/A
/-lowering-/.f6492.5%
Simplified92.5%
Final simplification80.3%
(FPCore (x1 x2)
:precision binary64
(let* ((t_0 (* 8.0 (* x1 (* x2 x2)))))
(if (<= x1 -1e+102)
(* -3.0 (* x1 (* x1 x1)))
(if (<= x1 -0.00022) t_0 (if (<= x1 2.55e-14) (- (* x2 -6.0) x1) t_0)))))
double code(double x1, double x2) {
double t_0 = 8.0 * (x1 * (x2 * x2));
double tmp;
if (x1 <= -1e+102) {
tmp = -3.0 * (x1 * (x1 * x1));
} else if (x1 <= -0.00022) {
tmp = t_0;
} else if (x1 <= 2.55e-14) {
tmp = (x2 * -6.0) - x1;
} else {
tmp = t_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) :: tmp
t_0 = 8.0d0 * (x1 * (x2 * x2))
if (x1 <= (-1d+102)) then
tmp = (-3.0d0) * (x1 * (x1 * x1))
else if (x1 <= (-0.00022d0)) then
tmp = t_0
else if (x1 <= 2.55d-14) then
tmp = (x2 * (-6.0d0)) - x1
else
tmp = t_0
end if
code = tmp
end function
public static double code(double x1, double x2) {
double t_0 = 8.0 * (x1 * (x2 * x2));
double tmp;
if (x1 <= -1e+102) {
tmp = -3.0 * (x1 * (x1 * x1));
} else if (x1 <= -0.00022) {
tmp = t_0;
} else if (x1 <= 2.55e-14) {
tmp = (x2 * -6.0) - x1;
} else {
tmp = t_0;
}
return tmp;
}
def code(x1, x2): t_0 = 8.0 * (x1 * (x2 * x2)) tmp = 0 if x1 <= -1e+102: tmp = -3.0 * (x1 * (x1 * x1)) elif x1 <= -0.00022: tmp = t_0 elif x1 <= 2.55e-14: tmp = (x2 * -6.0) - x1 else: tmp = t_0 return tmp
function code(x1, x2) t_0 = Float64(8.0 * Float64(x1 * Float64(x2 * x2))) tmp = 0.0 if (x1 <= -1e+102) tmp = Float64(-3.0 * Float64(x1 * Float64(x1 * x1))); elseif (x1 <= -0.00022) tmp = t_0; elseif (x1 <= 2.55e-14) tmp = Float64(Float64(x2 * -6.0) - x1); else tmp = t_0; end return tmp end
function tmp_2 = code(x1, x2) t_0 = 8.0 * (x1 * (x2 * x2)); tmp = 0.0; if (x1 <= -1e+102) tmp = -3.0 * (x1 * (x1 * x1)); elseif (x1 <= -0.00022) tmp = t_0; elseif (x1 <= 2.55e-14) tmp = (x2 * -6.0) - x1; else tmp = t_0; end tmp_2 = tmp; end
code[x1_, x2_] := Block[{t$95$0 = N[(8.0 * N[(x1 * N[(x2 * x2), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[x1, -1e+102], N[(-3.0 * N[(x1 * N[(x1 * x1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[x1, -0.00022], t$95$0, If[LessEqual[x1, 2.55e-14], N[(N[(x2 * -6.0), $MachinePrecision] - x1), $MachinePrecision], t$95$0]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := 8 \cdot \left(x1 \cdot \left(x2 \cdot x2\right)\right)\\
\mathbf{if}\;x1 \leq -1 \cdot 10^{+102}:\\
\;\;\;\;-3 \cdot \left(x1 \cdot \left(x1 \cdot x1\right)\right)\\
\mathbf{elif}\;x1 \leq -0.00022:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;x1 \leq 2.55 \cdot 10^{-14}:\\
\;\;\;\;x2 \cdot -6 - x1\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
if x1 < -9.99999999999999977e101Initial program 0.0%
Simplified0.0%
Applied egg-rr0.0%
Taylor expanded in x1 around inf
*-lowering-*.f64N/A
metadata-evalN/A
pow-sqrN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
associate-*r/N/A
metadata-evalN/A
distribute-neg-fracN/A
metadata-evalN/A
/-lowering-/.f64100.0%
Simplified100.0%
Taylor expanded in x1 around 0
*-commutativeN/A
*-lowering-*.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64100.0%
Simplified100.0%
if -9.99999999999999977e101 < x1 < -2.20000000000000008e-4 or 2.5499999999999999e-14 < x1 Initial program 61.7%
Simplified61.7%
Taylor expanded in x1 around 0
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
*-lowering-*.f6415.1%
Simplified15.1%
Taylor expanded in x2 around inf
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6431.5%
Simplified31.5%
if -2.20000000000000008e-4 < x1 < 2.5499999999999999e-14Initial program 99.5%
Simplified99.5%
Taylor expanded in x1 around inf
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6475.4%
Simplified75.4%
Taylor expanded in x1 around 0
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
*-commutativeN/A
*-lowering-*.f6475.5%
Simplified75.5%
Final simplification65.8%
(FPCore (x1 x2) :precision binary64 (if (<= x1 -1.5e+19) (+ x1 (* 6.0 (* x1 (* x1 (* x1 x1))))) (if (<= x1 6.5) (- (* x2 -6.0) x1) (+ x1 (* x1 (* x1 (* (* x1 x1) 6.0)))))))
double code(double x1, double x2) {
double tmp;
if (x1 <= -1.5e+19) {
tmp = x1 + (6.0 * (x1 * (x1 * (x1 * x1))));
} else if (x1 <= 6.5) {
tmp = (x2 * -6.0) - x1;
} else {
tmp = x1 + (x1 * (x1 * ((x1 * x1) * 6.0)));
}
return tmp;
}
real(8) function code(x1, x2)
real(8), intent (in) :: x1
real(8), intent (in) :: x2
real(8) :: tmp
if (x1 <= (-1.5d+19)) then
tmp = x1 + (6.0d0 * (x1 * (x1 * (x1 * x1))))
else if (x1 <= 6.5d0) then
tmp = (x2 * (-6.0d0)) - x1
else
tmp = x1 + (x1 * (x1 * ((x1 * x1) * 6.0d0)))
end if
code = tmp
end function
public static double code(double x1, double x2) {
double tmp;
if (x1 <= -1.5e+19) {
tmp = x1 + (6.0 * (x1 * (x1 * (x1 * x1))));
} else if (x1 <= 6.5) {
tmp = (x2 * -6.0) - x1;
} else {
tmp = x1 + (x1 * (x1 * ((x1 * x1) * 6.0)));
}
return tmp;
}
def code(x1, x2): tmp = 0 if x1 <= -1.5e+19: tmp = x1 + (6.0 * (x1 * (x1 * (x1 * x1)))) elif x1 <= 6.5: tmp = (x2 * -6.0) - x1 else: tmp = x1 + (x1 * (x1 * ((x1 * x1) * 6.0))) return tmp
function code(x1, x2) tmp = 0.0 if (x1 <= -1.5e+19) tmp = Float64(x1 + Float64(6.0 * Float64(x1 * Float64(x1 * Float64(x1 * x1))))); elseif (x1 <= 6.5) tmp = Float64(Float64(x2 * -6.0) - x1); else tmp = Float64(x1 + Float64(x1 * Float64(x1 * Float64(Float64(x1 * x1) * 6.0)))); end return tmp end
function tmp_2 = code(x1, x2) tmp = 0.0; if (x1 <= -1.5e+19) tmp = x1 + (6.0 * (x1 * (x1 * (x1 * x1)))); elseif (x1 <= 6.5) tmp = (x2 * -6.0) - x1; else tmp = x1 + (x1 * (x1 * ((x1 * x1) * 6.0))); end tmp_2 = tmp; end
code[x1_, x2_] := If[LessEqual[x1, -1.5e+19], N[(x1 + N[(6.0 * N[(x1 * N[(x1 * N[(x1 * x1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[x1, 6.5], N[(N[(x2 * -6.0), $MachinePrecision] - x1), $MachinePrecision], N[(x1 + N[(x1 * N[(x1 * N[(N[(x1 * x1), $MachinePrecision] * 6.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x1 \leq -1.5 \cdot 10^{+19}:\\
\;\;\;\;x1 + 6 \cdot \left(x1 \cdot \left(x1 \cdot \left(x1 \cdot x1\right)\right)\right)\\
\mathbf{elif}\;x1 \leq 6.5:\\
\;\;\;\;x2 \cdot -6 - x1\\
\mathbf{else}:\\
\;\;\;\;x1 + x1 \cdot \left(x1 \cdot \left(\left(x1 \cdot x1\right) \cdot 6\right)\right)\\
\end{array}
\end{array}
if x1 < -1.5e19Initial program 28.4%
Taylor expanded in x1 around inf
*-lowering-*.f64N/A
pow-lowering-pow.f6493.0%
Simplified93.0%
metadata-evalN/A
pow-prod-upN/A
pow2N/A
pow2N/A
*-commutativeN/A
*-lowering-*.f64N/A
associate-*l*N/A
cube-multN/A
*-lowering-*.f64N/A
cube-multN/A
*-lowering-*.f64N/A
*-lowering-*.f6493.0%
Applied egg-rr93.0%
if -1.5e19 < x1 < 6.5Initial program 99.5%
Simplified99.5%
Taylor expanded in x1 around inf
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6470.7%
Simplified70.7%
Taylor expanded in x1 around 0
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
*-commutativeN/A
*-lowering-*.f6471.0%
Simplified71.0%
if 6.5 < x1 Initial program 43.2%
Taylor expanded in x1 around inf
*-lowering-*.f64N/A
pow-lowering-pow.f6492.5%
Simplified92.5%
+-commutativeN/A
+-lowering-+.f64N/A
metadata-evalN/A
pow-prod-upN/A
pow2N/A
pow2N/A
*-commutativeN/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f6492.4%
Applied egg-rr92.4%
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
metadata-evalN/A
distribute-lft-neg-inN/A
*-commutativeN/A
*-lowering-*.f64N/A
distribute-rgt-neg-inN/A
metadata-evalN/A
*-lowering-*.f64N/A
*-lowering-*.f6492.4%
Applied egg-rr92.4%
Final simplification80.3%
(FPCore (x1 x2)
:precision binary64
(let* ((t_0 (* x1 (* x1 x1))))
(if (<= x1 -1.5e+19)
(+ x1 (* 6.0 (* x1 t_0)))
(if (<= x1 6.5) (- (* x2 -6.0) x1) (* t_0 (+ -3.0 (* x1 6.0)))))))
double code(double x1, double x2) {
double t_0 = x1 * (x1 * x1);
double tmp;
if (x1 <= -1.5e+19) {
tmp = x1 + (6.0 * (x1 * t_0));
} else if (x1 <= 6.5) {
tmp = (x2 * -6.0) - x1;
} else {
tmp = t_0 * (-3.0 + (x1 * 6.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) :: tmp
t_0 = x1 * (x1 * x1)
if (x1 <= (-1.5d+19)) then
tmp = x1 + (6.0d0 * (x1 * t_0))
else if (x1 <= 6.5d0) then
tmp = (x2 * (-6.0d0)) - x1
else
tmp = t_0 * ((-3.0d0) + (x1 * 6.0d0))
end if
code = tmp
end function
public static double code(double x1, double x2) {
double t_0 = x1 * (x1 * x1);
double tmp;
if (x1 <= -1.5e+19) {
tmp = x1 + (6.0 * (x1 * t_0));
} else if (x1 <= 6.5) {
tmp = (x2 * -6.0) - x1;
} else {
tmp = t_0 * (-3.0 + (x1 * 6.0));
}
return tmp;
}
def code(x1, x2): t_0 = x1 * (x1 * x1) tmp = 0 if x1 <= -1.5e+19: tmp = x1 + (6.0 * (x1 * t_0)) elif x1 <= 6.5: tmp = (x2 * -6.0) - x1 else: tmp = t_0 * (-3.0 + (x1 * 6.0)) return tmp
function code(x1, x2) t_0 = Float64(x1 * Float64(x1 * x1)) tmp = 0.0 if (x1 <= -1.5e+19) tmp = Float64(x1 + Float64(6.0 * Float64(x1 * t_0))); elseif (x1 <= 6.5) tmp = Float64(Float64(x2 * -6.0) - x1); else tmp = Float64(t_0 * Float64(-3.0 + Float64(x1 * 6.0))); end return tmp end
function tmp_2 = code(x1, x2) t_0 = x1 * (x1 * x1); tmp = 0.0; if (x1 <= -1.5e+19) tmp = x1 + (6.0 * (x1 * t_0)); elseif (x1 <= 6.5) tmp = (x2 * -6.0) - x1; else tmp = t_0 * (-3.0 + (x1 * 6.0)); end tmp_2 = tmp; end
code[x1_, x2_] := Block[{t$95$0 = N[(x1 * N[(x1 * x1), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[x1, -1.5e+19], N[(x1 + N[(6.0 * N[(x1 * t$95$0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[x1, 6.5], N[(N[(x2 * -6.0), $MachinePrecision] - x1), $MachinePrecision], N[(t$95$0 * N[(-3.0 + N[(x1 * 6.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := x1 \cdot \left(x1 \cdot x1\right)\\
\mathbf{if}\;x1 \leq -1.5 \cdot 10^{+19}:\\
\;\;\;\;x1 + 6 \cdot \left(x1 \cdot t\_0\right)\\
\mathbf{elif}\;x1 \leq 6.5:\\
\;\;\;\;x2 \cdot -6 - x1\\
\mathbf{else}:\\
\;\;\;\;t\_0 \cdot \left(-3 + x1 \cdot 6\right)\\
\end{array}
\end{array}
if x1 < -1.5e19Initial program 28.4%
Taylor expanded in x1 around inf
*-lowering-*.f64N/A
pow-lowering-pow.f6493.0%
Simplified93.0%
metadata-evalN/A
pow-prod-upN/A
pow2N/A
pow2N/A
*-commutativeN/A
*-lowering-*.f64N/A
associate-*l*N/A
cube-multN/A
*-lowering-*.f64N/A
cube-multN/A
*-lowering-*.f64N/A
*-lowering-*.f6493.0%
Applied egg-rr93.0%
if -1.5e19 < x1 < 6.5Initial program 99.5%
Simplified99.5%
Taylor expanded in x1 around inf
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6470.7%
Simplified70.7%
Taylor expanded in x1 around 0
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
*-commutativeN/A
*-lowering-*.f6471.0%
Simplified71.0%
if 6.5 < x1 Initial program 43.2%
Simplified43.2%
Applied egg-rr43.2%
Taylor expanded in x1 around inf
*-lowering-*.f64N/A
metadata-evalN/A
pow-sqrN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
associate-*r/N/A
metadata-evalN/A
distribute-neg-fracN/A
metadata-evalN/A
/-lowering-/.f6492.5%
Simplified92.5%
Taylor expanded in x1 around 0
*-lowering-*.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6492.4%
Simplified92.4%
Final simplification80.3%
(FPCore (x1 x2) :precision binary64 (let* ((t_0 (* (* x1 (* x1 x1)) (+ -3.0 (* x1 6.0))))) (if (<= x1 -1.5e+19) t_0 (if (<= x1 6.5) (- (* x2 -6.0) x1) t_0))))
double code(double x1, double x2) {
double t_0 = (x1 * (x1 * x1)) * (-3.0 + (x1 * 6.0));
double tmp;
if (x1 <= -1.5e+19) {
tmp = t_0;
} else if (x1 <= 6.5) {
tmp = (x2 * -6.0) - x1;
} else {
tmp = t_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) :: tmp
t_0 = (x1 * (x1 * x1)) * ((-3.0d0) + (x1 * 6.0d0))
if (x1 <= (-1.5d+19)) then
tmp = t_0
else if (x1 <= 6.5d0) then
tmp = (x2 * (-6.0d0)) - x1
else
tmp = t_0
end if
code = tmp
end function
public static double code(double x1, double x2) {
double t_0 = (x1 * (x1 * x1)) * (-3.0 + (x1 * 6.0));
double tmp;
if (x1 <= -1.5e+19) {
tmp = t_0;
} else if (x1 <= 6.5) {
tmp = (x2 * -6.0) - x1;
} else {
tmp = t_0;
}
return tmp;
}
def code(x1, x2): t_0 = (x1 * (x1 * x1)) * (-3.0 + (x1 * 6.0)) tmp = 0 if x1 <= -1.5e+19: tmp = t_0 elif x1 <= 6.5: tmp = (x2 * -6.0) - x1 else: tmp = t_0 return tmp
function code(x1, x2) t_0 = Float64(Float64(x1 * Float64(x1 * x1)) * Float64(-3.0 + Float64(x1 * 6.0))) tmp = 0.0 if (x1 <= -1.5e+19) tmp = t_0; elseif (x1 <= 6.5) tmp = Float64(Float64(x2 * -6.0) - x1); else tmp = t_0; end return tmp end
function tmp_2 = code(x1, x2) t_0 = (x1 * (x1 * x1)) * (-3.0 + (x1 * 6.0)); tmp = 0.0; if (x1 <= -1.5e+19) tmp = t_0; elseif (x1 <= 6.5) tmp = (x2 * -6.0) - x1; else tmp = t_0; end tmp_2 = tmp; end
code[x1_, x2_] := Block[{t$95$0 = N[(N[(x1 * N[(x1 * x1), $MachinePrecision]), $MachinePrecision] * N[(-3.0 + N[(x1 * 6.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[x1, -1.5e+19], t$95$0, If[LessEqual[x1, 6.5], N[(N[(x2 * -6.0), $MachinePrecision] - x1), $MachinePrecision], t$95$0]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \left(x1 \cdot \left(x1 \cdot x1\right)\right) \cdot \left(-3 + x1 \cdot 6\right)\\
\mathbf{if}\;x1 \leq -1.5 \cdot 10^{+19}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;x1 \leq 6.5:\\
\;\;\;\;x2 \cdot -6 - x1\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
if x1 < -1.5e19 or 6.5 < x1 Initial program 35.6%
Simplified35.6%
Applied egg-rr35.6%
Taylor expanded in x1 around inf
*-lowering-*.f64N/A
metadata-evalN/A
pow-sqrN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
sub-negN/A
+-lowering-+.f64N/A
associate-*r/N/A
metadata-evalN/A
distribute-neg-fracN/A
metadata-evalN/A
/-lowering-/.f6492.7%
Simplified92.7%
Taylor expanded in x1 around 0
*-lowering-*.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6492.7%
Simplified92.7%
if -1.5e19 < x1 < 6.5Initial program 99.5%
Simplified99.5%
Taylor expanded in x1 around inf
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6470.7%
Simplified70.7%
Taylor expanded in x1 around 0
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
*-commutativeN/A
*-lowering-*.f6471.0%
Simplified71.0%
Final simplification80.2%
(FPCore (x1 x2) :precision binary64 (let* ((t_0 (* x1 (+ (* 6.0 (* x1 (* x1 x1))) 1.0)))) (if (<= x1 -1.5e+19) t_0 (if (<= x1 6.5) (- (* x2 -6.0) x1) t_0))))
double code(double x1, double x2) {
double t_0 = x1 * ((6.0 * (x1 * (x1 * x1))) + 1.0);
double tmp;
if (x1 <= -1.5e+19) {
tmp = t_0;
} else if (x1 <= 6.5) {
tmp = (x2 * -6.0) - x1;
} else {
tmp = t_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) :: tmp
t_0 = x1 * ((6.0d0 * (x1 * (x1 * x1))) + 1.0d0)
if (x1 <= (-1.5d+19)) then
tmp = t_0
else if (x1 <= 6.5d0) then
tmp = (x2 * (-6.0d0)) - x1
else
tmp = t_0
end if
code = tmp
end function
public static double code(double x1, double x2) {
double t_0 = x1 * ((6.0 * (x1 * (x1 * x1))) + 1.0);
double tmp;
if (x1 <= -1.5e+19) {
tmp = t_0;
} else if (x1 <= 6.5) {
tmp = (x2 * -6.0) - x1;
} else {
tmp = t_0;
}
return tmp;
}
def code(x1, x2): t_0 = x1 * ((6.0 * (x1 * (x1 * x1))) + 1.0) tmp = 0 if x1 <= -1.5e+19: tmp = t_0 elif x1 <= 6.5: tmp = (x2 * -6.0) - x1 else: tmp = t_0 return tmp
function code(x1, x2) t_0 = Float64(x1 * Float64(Float64(6.0 * Float64(x1 * Float64(x1 * x1))) + 1.0)) tmp = 0.0 if (x1 <= -1.5e+19) tmp = t_0; elseif (x1 <= 6.5) tmp = Float64(Float64(x2 * -6.0) - x1); else tmp = t_0; end return tmp end
function tmp_2 = code(x1, x2) t_0 = x1 * ((6.0 * (x1 * (x1 * x1))) + 1.0); tmp = 0.0; if (x1 <= -1.5e+19) tmp = t_0; elseif (x1 <= 6.5) tmp = (x2 * -6.0) - x1; else tmp = t_0; end tmp_2 = tmp; end
code[x1_, x2_] := Block[{t$95$0 = N[(x1 * N[(N[(6.0 * N[(x1 * N[(x1 * x1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[x1, -1.5e+19], t$95$0, If[LessEqual[x1, 6.5], N[(N[(x2 * -6.0), $MachinePrecision] - x1), $MachinePrecision], t$95$0]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := x1 \cdot \left(6 \cdot \left(x1 \cdot \left(x1 \cdot x1\right)\right) + 1\right)\\
\mathbf{if}\;x1 \leq -1.5 \cdot 10^{+19}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;x1 \leq 6.5:\\
\;\;\;\;x2 \cdot -6 - x1\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
if x1 < -1.5e19 or 6.5 < x1 Initial program 35.6%
Taylor expanded in x1 around inf
*-lowering-*.f64N/A
pow-lowering-pow.f6492.7%
Simplified92.7%
Taylor expanded in x1 around 0
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6492.7%
Simplified92.7%
if -1.5e19 < x1 < 6.5Initial program 99.5%
Simplified99.5%
Taylor expanded in x1 around inf
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6470.7%
Simplified70.7%
Taylor expanded in x1 around 0
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
*-commutativeN/A
*-lowering-*.f6471.0%
Simplified71.0%
Final simplification80.2%
(FPCore (x1 x2) :precision binary64 (let* ((t_0 (* 6.0 (* (* x1 x1) (* x1 x1))))) (if (<= x1 -1.5e+19) t_0 (if (<= x1 6.5) (- (* x2 -6.0) x1) t_0))))
double code(double x1, double x2) {
double t_0 = 6.0 * ((x1 * x1) * (x1 * x1));
double tmp;
if (x1 <= -1.5e+19) {
tmp = t_0;
} else if (x1 <= 6.5) {
tmp = (x2 * -6.0) - x1;
} else {
tmp = t_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) :: tmp
t_0 = 6.0d0 * ((x1 * x1) * (x1 * x1))
if (x1 <= (-1.5d+19)) then
tmp = t_0
else if (x1 <= 6.5d0) then
tmp = (x2 * (-6.0d0)) - x1
else
tmp = t_0
end if
code = tmp
end function
public static double code(double x1, double x2) {
double t_0 = 6.0 * ((x1 * x1) * (x1 * x1));
double tmp;
if (x1 <= -1.5e+19) {
tmp = t_0;
} else if (x1 <= 6.5) {
tmp = (x2 * -6.0) - x1;
} else {
tmp = t_0;
}
return tmp;
}
def code(x1, x2): t_0 = 6.0 * ((x1 * x1) * (x1 * x1)) tmp = 0 if x1 <= -1.5e+19: tmp = t_0 elif x1 <= 6.5: tmp = (x2 * -6.0) - x1 else: tmp = t_0 return tmp
function code(x1, x2) t_0 = Float64(6.0 * Float64(Float64(x1 * x1) * Float64(x1 * x1))) tmp = 0.0 if (x1 <= -1.5e+19) tmp = t_0; elseif (x1 <= 6.5) tmp = Float64(Float64(x2 * -6.0) - x1); else tmp = t_0; end return tmp end
function tmp_2 = code(x1, x2) t_0 = 6.0 * ((x1 * x1) * (x1 * x1)); tmp = 0.0; if (x1 <= -1.5e+19) tmp = t_0; elseif (x1 <= 6.5) tmp = (x2 * -6.0) - x1; else tmp = t_0; end tmp_2 = tmp; end
code[x1_, x2_] := Block[{t$95$0 = N[(6.0 * N[(N[(x1 * x1), $MachinePrecision] * N[(x1 * x1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[x1, -1.5e+19], t$95$0, If[LessEqual[x1, 6.5], N[(N[(x2 * -6.0), $MachinePrecision] - x1), $MachinePrecision], t$95$0]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := 6 \cdot \left(\left(x1 \cdot x1\right) \cdot \left(x1 \cdot x1\right)\right)\\
\mathbf{if}\;x1 \leq -1.5 \cdot 10^{+19}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;x1 \leq 6.5:\\
\;\;\;\;x2 \cdot -6 - x1\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
if x1 < -1.5e19 or 6.5 < x1 Initial program 35.6%
Simplified35.6%
Applied egg-rr35.6%
Taylor expanded in x1 around inf
*-lowering-*.f64N/A
metadata-evalN/A
pow-sqrN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6492.6%
Simplified92.6%
if -1.5e19 < x1 < 6.5Initial program 99.5%
Simplified99.5%
Taylor expanded in x1 around inf
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6470.7%
Simplified70.7%
Taylor expanded in x1 around 0
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
*-commutativeN/A
*-lowering-*.f6471.0%
Simplified71.0%
(FPCore (x1 x2) :precision binary64 (if (<= x1 -3.05e+35) (* 8.0 (* x2 (* x1 x1))) (if (<= x1 5.2e-14) (- (* x2 -6.0) x1) (* 8.0 (* x1 (* x2 x2))))))
double code(double x1, double x2) {
double tmp;
if (x1 <= -3.05e+35) {
tmp = 8.0 * (x2 * (x1 * x1));
} else if (x1 <= 5.2e-14) {
tmp = (x2 * -6.0) - x1;
} else {
tmp = 8.0 * (x1 * (x2 * x2));
}
return tmp;
}
real(8) function code(x1, x2)
real(8), intent (in) :: x1
real(8), intent (in) :: x2
real(8) :: tmp
if (x1 <= (-3.05d+35)) then
tmp = 8.0d0 * (x2 * (x1 * x1))
else if (x1 <= 5.2d-14) then
tmp = (x2 * (-6.0d0)) - x1
else
tmp = 8.0d0 * (x1 * (x2 * x2))
end if
code = tmp
end function
public static double code(double x1, double x2) {
double tmp;
if (x1 <= -3.05e+35) {
tmp = 8.0 * (x2 * (x1 * x1));
} else if (x1 <= 5.2e-14) {
tmp = (x2 * -6.0) - x1;
} else {
tmp = 8.0 * (x1 * (x2 * x2));
}
return tmp;
}
def code(x1, x2): tmp = 0 if x1 <= -3.05e+35: tmp = 8.0 * (x2 * (x1 * x1)) elif x1 <= 5.2e-14: tmp = (x2 * -6.0) - x1 else: tmp = 8.0 * (x1 * (x2 * x2)) return tmp
function code(x1, x2) tmp = 0.0 if (x1 <= -3.05e+35) tmp = Float64(8.0 * Float64(x2 * Float64(x1 * x1))); elseif (x1 <= 5.2e-14) tmp = Float64(Float64(x2 * -6.0) - x1); else tmp = Float64(8.0 * Float64(x1 * Float64(x2 * x2))); end return tmp end
function tmp_2 = code(x1, x2) tmp = 0.0; if (x1 <= -3.05e+35) tmp = 8.0 * (x2 * (x1 * x1)); elseif (x1 <= 5.2e-14) tmp = (x2 * -6.0) - x1; else tmp = 8.0 * (x1 * (x2 * x2)); end tmp_2 = tmp; end
code[x1_, x2_] := If[LessEqual[x1, -3.05e+35], N[(8.0 * N[(x2 * N[(x1 * x1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[x1, 5.2e-14], N[(N[(x2 * -6.0), $MachinePrecision] - x1), $MachinePrecision], N[(8.0 * N[(x1 * N[(x2 * x2), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x1 \leq -3.05 \cdot 10^{+35}:\\
\;\;\;\;8 \cdot \left(x2 \cdot \left(x1 \cdot x1\right)\right)\\
\mathbf{elif}\;x1 \leq 5.2 \cdot 10^{-14}:\\
\;\;\;\;x2 \cdot -6 - x1\\
\mathbf{else}:\\
\;\;\;\;8 \cdot \left(x1 \cdot \left(x2 \cdot x2\right)\right)\\
\end{array}
\end{array}
if x1 < -3.04999999999999989e35Initial program 25.8%
Simplified25.8%
Applied egg-rr25.7%
Taylor expanded in x2 around -inf
Simplified13.2%
Taylor expanded in x1 around inf
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6428.1%
Simplified28.1%
if -3.04999999999999989e35 < x1 < 5.19999999999999993e-14Initial program 99.5%
Simplified99.5%
Taylor expanded in x1 around inf
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6473.6%
Simplified73.6%
Taylor expanded in x1 around 0
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
*-commutativeN/A
*-lowering-*.f6472.6%
Simplified72.6%
if 5.19999999999999993e-14 < x1 Initial program 48.9%
Simplified48.9%
Taylor expanded in x1 around 0
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
*-lowering-*.f649.2%
Simplified9.2%
Taylor expanded in x2 around inf
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6432.7%
Simplified32.7%
(FPCore (x1 x2) :precision binary64 (let* ((t_0 (* 8.0 (* x1 (* x2 x2))))) (if (<= x2 -8.5e+79) t_0 (if (<= x2 1.55e+90) (- (* x2 -6.0) x1) t_0))))
double code(double x1, double x2) {
double t_0 = 8.0 * (x1 * (x2 * x2));
double tmp;
if (x2 <= -8.5e+79) {
tmp = t_0;
} else if (x2 <= 1.55e+90) {
tmp = (x2 * -6.0) - x1;
} else {
tmp = t_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) :: tmp
t_0 = 8.0d0 * (x1 * (x2 * x2))
if (x2 <= (-8.5d+79)) then
tmp = t_0
else if (x2 <= 1.55d+90) then
tmp = (x2 * (-6.0d0)) - x1
else
tmp = t_0
end if
code = tmp
end function
public static double code(double x1, double x2) {
double t_0 = 8.0 * (x1 * (x2 * x2));
double tmp;
if (x2 <= -8.5e+79) {
tmp = t_0;
} else if (x2 <= 1.55e+90) {
tmp = (x2 * -6.0) - x1;
} else {
tmp = t_0;
}
return tmp;
}
def code(x1, x2): t_0 = 8.0 * (x1 * (x2 * x2)) tmp = 0 if x2 <= -8.5e+79: tmp = t_0 elif x2 <= 1.55e+90: tmp = (x2 * -6.0) - x1 else: tmp = t_0 return tmp
function code(x1, x2) t_0 = Float64(8.0 * Float64(x1 * Float64(x2 * x2))) tmp = 0.0 if (x2 <= -8.5e+79) tmp = t_0; elseif (x2 <= 1.55e+90) tmp = Float64(Float64(x2 * -6.0) - x1); else tmp = t_0; end return tmp end
function tmp_2 = code(x1, x2) t_0 = 8.0 * (x1 * (x2 * x2)); tmp = 0.0; if (x2 <= -8.5e+79) tmp = t_0; elseif (x2 <= 1.55e+90) tmp = (x2 * -6.0) - x1; else tmp = t_0; end tmp_2 = tmp; end
code[x1_, x2_] := Block[{t$95$0 = N[(8.0 * N[(x1 * N[(x2 * x2), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[x2, -8.5e+79], t$95$0, If[LessEqual[x2, 1.55e+90], N[(N[(x2 * -6.0), $MachinePrecision] - x1), $MachinePrecision], t$95$0]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := 8 \cdot \left(x1 \cdot \left(x2 \cdot x2\right)\right)\\
\mathbf{if}\;x2 \leq -8.5 \cdot 10^{+79}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;x2 \leq 1.55 \cdot 10^{+90}:\\
\;\;\;\;x2 \cdot -6 - x1\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
if x2 < -8.4999999999999998e79 or 1.54999999999999994e90 < x2 Initial program 74.0%
Simplified74.0%
Taylor expanded in x1 around 0
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-lowering-+.f64N/A
*-lowering-*.f6464.2%
Simplified64.2%
Taylor expanded in x2 around inf
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6451.9%
Simplified51.9%
if -8.4999999999999998e79 < x2 < 1.54999999999999994e90Initial program 71.4%
Simplified71.4%
Taylor expanded in x1 around inf
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6467.7%
Simplified67.7%
Taylor expanded in x1 around 0
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
*-commutativeN/A
*-lowering-*.f6453.1%
Simplified53.1%
(FPCore (x1 x2) :precision binary64 (- (* x2 -6.0) x1))
double code(double x1, double x2) {
return (x2 * -6.0) - x1;
}
real(8) function code(x1, x2)
real(8), intent (in) :: x1
real(8), intent (in) :: x2
code = (x2 * (-6.0d0)) - x1
end function
public static double code(double x1, double x2) {
return (x2 * -6.0) - x1;
}
def code(x1, x2): return (x2 * -6.0) - x1
function code(x1, x2) return Float64(Float64(x2 * -6.0) - x1) end
function tmp = code(x1, x2) tmp = (x2 * -6.0) - x1; end
code[x1_, x2_] := N[(N[(x2 * -6.0), $MachinePrecision] - x1), $MachinePrecision]
\begin{array}{l}
\\
x2 \cdot -6 - x1
\end{array}
Initial program 72.3%
Simplified72.3%
Taylor expanded in x1 around inf
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6453.1%
Simplified53.1%
Taylor expanded in x1 around 0
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
*-commutativeN/A
*-lowering-*.f6442.2%
Simplified42.2%
(FPCore (x1 x2) :precision binary64 (+ x1 (* x2 -6.0)))
double code(double x1, double x2) {
return x1 + (x2 * -6.0);
}
real(8) function code(x1, x2)
real(8), intent (in) :: x1
real(8), intent (in) :: x2
code = x1 + (x2 * (-6.0d0))
end function
public static double code(double x1, double x2) {
return x1 + (x2 * -6.0);
}
def code(x1, x2): return x1 + (x2 * -6.0)
function code(x1, x2) return Float64(x1 + Float64(x2 * -6.0)) end
function tmp = code(x1, x2) tmp = x1 + (x2 * -6.0); end
code[x1_, x2_] := N[(x1 + N[(x2 * -6.0), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
x1 + x2 \cdot -6
\end{array}
Initial program 72.3%
Taylor expanded in x1 around 0
*-commutativeN/A
*-lowering-*.f6430.9%
Simplified30.9%
(FPCore (x1 x2) :precision binary64 (* x2 -6.0))
double code(double x1, double x2) {
return x2 * -6.0;
}
real(8) function code(x1, x2)
real(8), intent (in) :: x1
real(8), intent (in) :: x2
code = x2 * (-6.0d0)
end function
public static double code(double x1, double x2) {
return x2 * -6.0;
}
def code(x1, x2): return x2 * -6.0
function code(x1, x2) return Float64(x2 * -6.0) end
function tmp = code(x1, x2) tmp = x2 * -6.0; end
code[x1_, x2_] := N[(x2 * -6.0), $MachinePrecision]
\begin{array}{l}
\\
x2 \cdot -6
\end{array}
Initial program 72.3%
Simplified72.3%
Applied egg-rr72.3%
Taylor expanded in x1 around 0
*-commutativeN/A
*-lowering-*.f6430.9%
Simplified30.9%
(FPCore (x1 x2) :precision binary64 (* x1 2.0))
double code(double x1, double x2) {
return x1 * 2.0;
}
real(8) function code(x1, x2)
real(8), intent (in) :: x1
real(8), intent (in) :: x2
code = x1 * 2.0d0
end function
public static double code(double x1, double x2) {
return x1 * 2.0;
}
def code(x1, x2): return x1 * 2.0
function code(x1, x2) return Float64(x1 * 2.0) end
function tmp = code(x1, x2) tmp = x1 * 2.0; end
code[x1_, x2_] := N[(x1 * 2.0), $MachinePrecision]
\begin{array}{l}
\\
x1 \cdot 2
\end{array}
Initial program 72.3%
Simplified72.3%
Taylor expanded in x1 around inf
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6453.1%
Simplified53.1%
Taylor expanded in x1 around 0
*-commutativeN/A
*-lowering-*.f6441.7%
Simplified41.7%
Taylor expanded in x1 around inf
*-commutativeN/A
*-lowering-*.f643.2%
Simplified3.2%
(FPCore (x1 x2) :precision binary64 x1)
double code(double x1, double x2) {
return x1;
}
real(8) function code(x1, x2)
real(8), intent (in) :: x1
real(8), intent (in) :: x2
code = x1
end function
public static double code(double x1, double x2) {
return x1;
}
def code(x1, x2): return x1
function code(x1, x2) return x1 end
function tmp = code(x1, x2) tmp = x1; end
code[x1_, x2_] := x1
\begin{array}{l}
\\
x1
\end{array}
Initial program 72.3%
Taylor expanded in x1 around 0
*-commutativeN/A
*-lowering-*.f6430.9%
Simplified30.9%
Taylor expanded in x1 around inf
Simplified3.1%
herbie shell --seed 2024161
(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))))))