Rosa's FloatVsDoubleBenchmark

Percentage Accurate: 70.6% → 99.5%

Time: 17.4s

Alternatives: 13

Speedup: 6.0×

• 46.3% of points produce a very large (infinite) output. You may want to add a precondition.

Specification

Math

\[\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

(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))))))

C

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)));
}

Fortran

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

Java

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)));
}

Python

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)))

Julia

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

MATLAB

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

Wolfram

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]]]]

Initial Program: 70.6% accurate, 1.0× speedup

Math

\[\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

(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))))))

C

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)));
}

Fortran

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

Java

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)));
}

Python

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)))

Julia

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

MATLAB

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

Wolfram

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]]]]

Alternative 1: 99.5% accurate, 0.5× speedup

Math

\[\begin{array}{l} \\ \begin{array}{l} t_0 := \left(3 \cdot x1\right) \cdot x1\\ t_1 := x1 \cdot x1 + 1\\ t_2 := \frac{\left(t\_0 + 2 \cdot x2\right) - x1}{t\_1}\\ t_3 := \frac{\mathsf{fma}\left(x1 \cdot x1, 3, 2 \cdot x2\right) - x1}{\mathsf{fma}\left(x1, x1, 1\right)}\\ \mathbf{if}\;x1 + \left(\left(\left(\left(\left(\left(\left(2 \cdot x1\right) \cdot t\_2\right) \cdot \left(t\_2 - 3\right) + \left(x1 \cdot x1\right) \cdot \left(4 \cdot t\_2 - 6\right)\right) \cdot t\_1 + t\_0 \cdot t\_2\right) + \left(x1 \cdot x1\right) \cdot x1\right) + x1\right) + 3 \cdot \frac{\left(t\_0 - 2 \cdot x2\right) - x1}{t\_1}\right) \leq \infty:\\ \;\;\;\;x1 + \mathsf{fma}\left(x1 \cdot x1, x1, \mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(t\_3, 4, -6\right) \cdot x1, x1, \left(t\_3 - 3\right) \cdot \left(\left(x1 \cdot 2\right) \cdot t\_3\right)\right), \mathsf{fma}\left(x1, x1, 1\right), \frac{\mathsf{fma}\left(x2, 2, t\_0\right) - x1}{\mathsf{fma}\left(x1, x1, 1\right)} \cdot t\_0\right) + \mathsf{fma}\left(\frac{\mathsf{fma}\left(-2, x2, t\_0\right) - x1}{\mathsf{fma}\left(x1, x1, 1\right)}, 3, x1\right)\right)\\ \mathbf{else}:\\ \;\;\;\;x1 + \left(\left(6 - \frac{3}{x1}\right) \cdot \left(x1 \cdot x1\right)\right) \cdot \left(x1 \cdot x1\right)\\ \end{array} \end{array} \]

FPCore

(FPCore (x1 x2)
 :precision binary64
 (let* ((t_0 (* (* 3.0 x1) x1))
        (t_1 (+ (* x1 x1) 1.0))
        (t_2 (/ (- (+ t_0 (* 2.0 x2)) x1) t_1))
        (t_3 (/ (- (fma (* x1 x1) 3.0 (* 2.0 x2)) x1) (fma x1 x1 1.0))))
   (if (<=
        (+
         x1
         (+
          (+
           (+
            (+
             (*
              (+
               (* (* (* 2.0 x1) t_2) (- t_2 3.0))
               (* (* x1 x1) (- (* 4.0 t_2) 6.0)))
              t_1)
             (* t_0 t_2))
            (* (* x1 x1) x1))
           x1)
          (* 3.0 (/ (- (- t_0 (* 2.0 x2)) x1) t_1))))
        INFINITY)
     (+
      x1
      (fma
       (* x1 x1)
       x1
       (+
        (fma
         (fma (* (fma t_3 4.0 -6.0) x1) x1 (* (- t_3 3.0) (* (* x1 2.0) t_3)))
         (fma x1 x1 1.0)
         (* (/ (- (fma x2 2.0 t_0) x1) (fma x1 x1 1.0)) t_0))
        (fma (/ (- (fma -2.0 x2 t_0) x1) (fma x1 x1 1.0)) 3.0 x1))))
     (+ x1 (* (* (- 6.0 (/ 3.0 x1)) (* x1 x1)) (* x1 x1))))))

C

double code(double x1, double x2) {
	double t_0 = (3.0 * x1) * x1;
	double t_1 = (x1 * x1) + 1.0;
	double t_2 = ((t_0 + (2.0 * x2)) - x1) / t_1;
	double t_3 = (fma((x1 * x1), 3.0, (2.0 * x2)) - x1) / fma(x1, x1, 1.0);
	double tmp;
	if ((x1 + (((((((((2.0 * x1) * t_2) * (t_2 - 3.0)) + ((x1 * x1) * ((4.0 * t_2) - 6.0))) * t_1) + (t_0 * t_2)) + ((x1 * x1) * x1)) + x1) + (3.0 * (((t_0 - (2.0 * x2)) - x1) / t_1)))) <= ((double) INFINITY)) {
		tmp = x1 + fma((x1 * x1), x1, (fma(fma((fma(t_3, 4.0, -6.0) * x1), x1, ((t_3 - 3.0) * ((x1 * 2.0) * t_3))), fma(x1, x1, 1.0), (((fma(x2, 2.0, t_0) - x1) / fma(x1, x1, 1.0)) * t_0)) + fma(((fma(-2.0, x2, t_0) - x1) / fma(x1, x1, 1.0)), 3.0, x1)));
	} else {
		tmp = x1 + (((6.0 - (3.0 / x1)) * (x1 * x1)) * (x1 * x1));
	}
	return tmp;
}

Julia

function code(x1, x2)
	t_0 = Float64(Float64(3.0 * x1) * x1)
	t_1 = Float64(Float64(x1 * x1) + 1.0)
	t_2 = Float64(Float64(Float64(t_0 + Float64(2.0 * x2)) - x1) / t_1)
	t_3 = Float64(Float64(fma(Float64(x1 * x1), 3.0, Float64(2.0 * x2)) - x1) / fma(x1, x1, 1.0))
	tmp = 0.0
	if (Float64(x1 + Float64(Float64(Float64(Float64(Float64(Float64(Float64(Float64(Float64(2.0 * x1) * t_2) * Float64(t_2 - 3.0)) + Float64(Float64(x1 * x1) * Float64(Float64(4.0 * t_2) - 6.0))) * t_1) + Float64(t_0 * t_2)) + Float64(Float64(x1 * x1) * x1)) + x1) + Float64(3.0 * Float64(Float64(Float64(t_0 - Float64(2.0 * x2)) - x1) / t_1)))) <= Inf)
		tmp = Float64(x1 + fma(Float64(x1 * x1), x1, Float64(fma(fma(Float64(fma(t_3, 4.0, -6.0) * x1), x1, Float64(Float64(t_3 - 3.0) * Float64(Float64(x1 * 2.0) * t_3))), fma(x1, x1, 1.0), Float64(Float64(Float64(fma(x2, 2.0, t_0) - x1) / fma(x1, x1, 1.0)) * t_0)) + fma(Float64(Float64(fma(-2.0, x2, t_0) - x1) / fma(x1, x1, 1.0)), 3.0, x1))));
	else
		tmp = Float64(x1 + Float64(Float64(Float64(6.0 - Float64(3.0 / x1)) * Float64(x1 * x1)) * Float64(x1 * x1)));
	end
	return tmp
end

Wolfram

code[x1_, x2_] := Block[{t$95$0 = N[(N[(3.0 * x1), $MachinePrecision] * x1), $MachinePrecision]}, Block[{t$95$1 = N[(N[(x1 * x1), $MachinePrecision] + 1.0), $MachinePrecision]}, Block[{t$95$2 = N[(N[(N[(t$95$0 + N[(2.0 * x2), $MachinePrecision]), $MachinePrecision] - x1), $MachinePrecision] / t$95$1), $MachinePrecision]}, Block[{t$95$3 = N[(N[(N[(N[(x1 * x1), $MachinePrecision] * 3.0 + N[(2.0 * x2), $MachinePrecision]), $MachinePrecision] - x1), $MachinePrecision] / N[(x1 * x1 + 1.0), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[N[(x1 + N[(N[(N[(N[(N[(N[(N[(N[(N[(2.0 * x1), $MachinePrecision] * t$95$2), $MachinePrecision] * N[(t$95$2 - 3.0), $MachinePrecision]), $MachinePrecision] + N[(N[(x1 * x1), $MachinePrecision] * N[(N[(4.0 * t$95$2), $MachinePrecision] - 6.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * t$95$1), $MachinePrecision] + N[(t$95$0 * t$95$2), $MachinePrecision]), $MachinePrecision] + N[(N[(x1 * x1), $MachinePrecision] * x1), $MachinePrecision]), $MachinePrecision] + x1), $MachinePrecision] + N[(3.0 * N[(N[(N[(t$95$0 - N[(2.0 * x2), $MachinePrecision]), $MachinePrecision] - x1), $MachinePrecision] / t$95$1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], Infinity], N[(x1 + N[(N[(x1 * x1), $MachinePrecision] * x1 + N[(N[(N[(N[(N[(t$95$3 * 4.0 + -6.0), $MachinePrecision] * x1), $MachinePrecision] * x1 + N[(N[(t$95$3 - 3.0), $MachinePrecision] * N[(N[(x1 * 2.0), $MachinePrecision] * t$95$3), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(x1 * x1 + 1.0), $MachinePrecision] + N[(N[(N[(N[(x2 * 2.0 + t$95$0), $MachinePrecision] - x1), $MachinePrecision] / N[(x1 * x1 + 1.0), $MachinePrecision]), $MachinePrecision] * t$95$0), $MachinePrecision]), $MachinePrecision] + N[(N[(N[(N[(-2.0 * x2 + t$95$0), $MachinePrecision] - x1), $MachinePrecision] / N[(x1 * x1 + 1.0), $MachinePrecision]), $MachinePrecision] * 3.0 + x1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x1 + N[(N[(N[(6.0 - N[(3.0 / x1), $MachinePrecision]), $MachinePrecision] * N[(x1 * x1), $MachinePrecision]), $MachinePrecision] * N[(x1 * x1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]]]

Derivation

1. Split input into 2 regimes

2. 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.0

   1. Initial program 99.4%

      \[x1 + \left(\left(\left(\left(\left(\left(\left(2 \cdot x1\right) \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 + 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \cdot \left(\frac{\left(\left(3 \cdot x1\right) \cdot x1 + 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1} - 3\right) + \left(x1 \cdot x1\right) \cdot \left(4 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 + 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1} - 6\right)\right) \cdot \left(x1 \cdot x1 + 1\right) + \left(\left(3 \cdot x1\right) \cdot x1\right) \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 + 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) + \left(x1 \cdot x1\right) \cdot x1\right) + x1\right) + 3 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 - 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \]
   2. Add Preprocessing

   4. Applied rewrites 99.7%

      \[\leadsto x1 + \color{blue}{\mathsf{fma}\left(x1 \cdot x1, x1, \mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(4, \frac{\mathsf{fma}\left(x2, 2, \left(3 \cdot x1\right) \cdot x1\right) - x1}{\mathsf{fma}\left(x1, x1, 1\right)}, -6\right), x1 \cdot x1, \left(\frac{\mathsf{fma}\left(x2, 2, \left(3 \cdot x1\right) \cdot x1\right) - x1}{\mathsf{fma}\left(x1, x1, 1\right)} - 3\right) \cdot \left(\frac{\mathsf{fma}\left(x2, 2, \left(3 \cdot x1\right) \cdot x1\right) - x1}{\mathsf{fma}\left(x1, x1, 1\right)} \cdot \left(x1 \cdot 2\right)\right)\right), \mathsf{fma}\left(x1, x1, 1\right), \frac{\mathsf{fma}\left(x2, 2, \left(3 \cdot x1\right) \cdot x1\right) - x1}{\mathsf{fma}\left(x1, x1, 1\right)} \cdot \left(\left(3 \cdot x1\right) \cdot x1\right)\right) + \mathsf{fma}\left(\frac{\mathsf{fma}\left(-2, x2, \left(3 \cdot x1\right) \cdot x1\right) - x1}{\mathsf{fma}\left(x1, x1, 1\right)}, 3, x1\right)\right)} \]

   6. Applied rewrites 99.7%

      \[\leadsto x1 + \mathsf{fma}\left(x1 \cdot x1, x1, \mathsf{fma}\left(\color{blue}{\mathsf{fma}\left(\mathsf{fma}\left(\frac{\mathsf{fma}\left(x1 \cdot x1, 3, 2 \cdot x2\right) - x1}{\mathsf{fma}\left(x1, x1, 1\right)}, 4, -6\right) \cdot x1, x1, \left(\frac{\mathsf{fma}\left(x1 \cdot x1, 3, 2 \cdot x2\right) - x1}{\mathsf{fma}\left(x1, x1, 1\right)} - 3\right) \cdot \left(\left(x1 \cdot 2\right) \cdot \frac{\mathsf{fma}\left(x1 \cdot x1, 3, 2 \cdot x2\right) - x1}{\mathsf{fma}\left(x1, x1, 1\right)}\right)\right)}, \mathsf{fma}\left(x1, x1, 1\right), \frac{\mathsf{fma}\left(x2, 2, \left(3 \cdot x1\right) \cdot x1\right) - x1}{\mathsf{fma}\left(x1, x1, 1\right)} \cdot \left(\left(3 \cdot x1\right) \cdot x1\right)\right) + \mathsf{fma}\left(\frac{\mathsf{fma}\left(-2, x2, \left(3 \cdot x1\right) \cdot x1\right) - x1}{\mathsf{fma}\left(x1, x1, 1\right)}, 3, x1\right)\right) \]

   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))))))

   1. Initial program 0.0%

      \[x1 + \left(\left(\left(\left(\left(\left(\left(2 \cdot x1\right) \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 + 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \cdot \left(\frac{\left(\left(3 \cdot x1\right) \cdot x1 + 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1} - 3\right) + \left(x1 \cdot x1\right) \cdot \left(4 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 + 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1} - 6\right)\right) \cdot \left(x1 \cdot x1 + 1\right) + \left(\left(3 \cdot x1\right) \cdot x1\right) \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 + 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) + \left(x1 \cdot x1\right) \cdot x1\right) + x1\right) + 3 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 - 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \]
   2. Add Preprocessing

   4. Applied rewrites 11.9%

      \[\leadsto x1 + \color{blue}{\mathsf{fma}\left(x1 \cdot x1, x1, \mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(4, \frac{\mathsf{fma}\left(x2, 2, \left(3 \cdot x1\right) \cdot x1\right) - x1}{\mathsf{fma}\left(x1, x1, 1\right)}, -6\right), x1 \cdot x1, \left(\frac{\mathsf{fma}\left(x2, 2, \left(3 \cdot x1\right) \cdot x1\right) - x1}{\mathsf{fma}\left(x1, x1, 1\right)} - 3\right) \cdot \left(\frac{\mathsf{fma}\left(x2, 2, \left(3 \cdot x1\right) \cdot x1\right) - x1}{\mathsf{fma}\left(x1, x1, 1\right)} \cdot \left(x1 \cdot 2\right)\right)\right), \mathsf{fma}\left(x1, x1, 1\right), \frac{\mathsf{fma}\left(x2, 2, \left(3 \cdot x1\right) \cdot x1\right) - x1}{\mathsf{fma}\left(x1, x1, 1\right)} \cdot \left(\left(3 \cdot x1\right) \cdot x1\right)\right) + \mathsf{fma}\left(\frac{\mathsf{fma}\left(-2, x2, \left(3 \cdot x1\right) \cdot x1\right) - x1}{\mathsf{fma}\left(x1, x1, 1\right)}, 3, x1\right)\right)} \]

   5. Taylor expanded in x1 around inf

      \[\leadsto x1 + \color{blue}{{x1}^{4} \cdot \left(6 - 3 \cdot \frac{1}{x1}\right)} \]
   6. Step-by-step derivation

      1. *-commutative N/A

         \[\leadsto x1 + \color{blue}{\left(6 - 3 \cdot \frac{1}{x1}\right) \cdot {x1}^{4}} \]

      2. lower-*.f64 N/A

         \[\leadsto x1 + \color{blue}{\left(6 - 3 \cdot \frac{1}{x1}\right) \cdot {x1}^{4}} \]

      3. lower--.f64 N/A

         \[\leadsto x1 + \color{blue}{\left(6 - 3 \cdot \frac{1}{x1}\right)} \cdot {x1}^{4} \]

      4. associate-*r/ N/A

         \[\leadsto x1 + \left(6 - \color{blue}{\frac{3 \cdot 1}{x1}}\right) \cdot {x1}^{4} \]

      5. metadata-eval N/A

         \[\leadsto x1 + \left(6 - \frac{\color{blue}{3}}{x1}\right) \cdot {x1}^{4} \]

      6. lower-/.f64 N/A

         \[\leadsto x1 + \left(6 - \color{blue}{\frac{3}{x1}}\right) \cdot {x1}^{4} \]

      7. lower-pow.f64 100.0

         \[\leadsto x1 + \left(6 - \frac{3}{x1}\right) \cdot \color{blue}{{x1}^{4}} \]

   7. Applied rewrites 100.0%

      \[\leadsto x1 + \color{blue}{\left(6 - \frac{3}{x1}\right) \cdot {x1}^{4}} \]
   8. Step-by-step derivation

   9. Applied rewrites 100.0%

      \[\leadsto x1 + \left(\left(6 - \frac{3}{x1}\right) \cdot \left(x1 \cdot x1\right)\right) \cdot \color{blue}{\left(x1 \cdot x1\right)} \]
3. Recombined 2 regimes into one program.
4. Add Preprocessing

Alternative 2: 94.5% accurate, 1.8× speedup

Math

\[\begin{array}{l} \\ \begin{array}{l} t_0 := \left(6 - \frac{3 - \frac{\mathsf{fma}\left(\mathsf{fma}\left(2, x2, -3\right), 4, 9\right)}{x1}}{x1}\right) \cdot {x1}^{4}\\ \mathbf{if}\;x1 \leq -2.8 \cdot 10^{+32}:\\ \;\;\;\;t\_0\\ \mathbf{elif}\;x1 \leq 1.65 \cdot 10^{+40}:\\ \;\;\;\;x1 + \left(\left(\left(\left(\frac{x1}{\mathsf{fma}\left(x1, x1, 1\right)} \cdot 8\right) \cdot x2\right) \cdot x2 + x1\right) + 3 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 - 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right)\\ \mathbf{else}:\\ \;\;\;\;x1 + t\_0\\ \end{array} \end{array} \]

FPCore

(FPCore (x1 x2)
 :precision binary64
 (let* ((t_0
         (*
          (- 6.0 (/ (- 3.0 (/ (fma (fma 2.0 x2 -3.0) 4.0 9.0) x1)) x1))
          (pow x1 4.0))))
   (if (<= x1 -2.8e+32)
     t_0
     (if (<= x1 1.65e+40)
       (+
        x1
        (+
         (+ (* (* (* (/ x1 (fma x1 x1 1.0)) 8.0) x2) x2) x1)
         (*
          3.0
          (/ (- (- (* (* 3.0 x1) x1) (* 2.0 x2)) x1) (+ (* x1 x1) 1.0)))))
       (+ x1 t_0)))))

C

double code(double x1, double x2) {
	double t_0 = (6.0 - ((3.0 - (fma(fma(2.0, x2, -3.0), 4.0, 9.0) / x1)) / x1)) * pow(x1, 4.0);
	double tmp;
	if (x1 <= -2.8e+32) {
		tmp = t_0;
	} else if (x1 <= 1.65e+40) {
		tmp = x1 + ((((((x1 / fma(x1, x1, 1.0)) * 8.0) * x2) * x2) + x1) + (3.0 * (((((3.0 * x1) * x1) - (2.0 * x2)) - x1) / ((x1 * x1) + 1.0))));
	} else {
		tmp = x1 + t_0;
	}
	return tmp;
}

Julia

function code(x1, x2)
	t_0 = Float64(Float64(6.0 - Float64(Float64(3.0 - Float64(fma(fma(2.0, x2, -3.0), 4.0, 9.0) / x1)) / x1)) * (x1 ^ 4.0))
	tmp = 0.0
	if (x1 <= -2.8e+32)
		tmp = t_0;
	elseif (x1 <= 1.65e+40)
		tmp = Float64(x1 + Float64(Float64(Float64(Float64(Float64(Float64(x1 / fma(x1, x1, 1.0)) * 8.0) * x2) * x2) + x1) + Float64(3.0 * Float64(Float64(Float64(Float64(Float64(3.0 * x1) * x1) - Float64(2.0 * x2)) - x1) / Float64(Float64(x1 * x1) + 1.0)))));
	else
		tmp = Float64(x1 + t_0);
	end
	return tmp
end

Wolfram

code[x1_, x2_] := Block[{t$95$0 = N[(N[(6.0 - N[(N[(3.0 - N[(N[(N[(2.0 * x2 + -3.0), $MachinePrecision] * 4.0 + 9.0), $MachinePrecision] / x1), $MachinePrecision]), $MachinePrecision] / x1), $MachinePrecision]), $MachinePrecision] * N[Power[x1, 4.0], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[x1, -2.8e+32], t$95$0, If[LessEqual[x1, 1.65e+40], N[(x1 + N[(N[(N[(N[(N[(N[(x1 / N[(x1 * x1 + 1.0), $MachinePrecision]), $MachinePrecision] * 8.0), $MachinePrecision] * x2), $MachinePrecision] * x2), $MachinePrecision] + x1), $MachinePrecision] + N[(3.0 * N[(N[(N[(N[(N[(3.0 * x1), $MachinePrecision] * x1), $MachinePrecision] - N[(2.0 * x2), $MachinePrecision]), $MachinePrecision] - x1), $MachinePrecision] / N[(N[(x1 * x1), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x1 + t$95$0), $MachinePrecision]]]]

Derivation

1. Split input into 3 regimes

2. if x1 < -2.8e32

   1. Initial program 21.5%

      \[x1 + \left(\left(\left(\left(\left(\left(\left(2 \cdot x1\right) \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 + 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \cdot \left(\frac{\left(\left(3 \cdot x1\right) \cdot x1 + 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1} - 3\right) + \left(x1 \cdot x1\right) \cdot \left(4 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 + 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1} - 6\right)\right) \cdot \left(x1 \cdot x1 + 1\right) + \left(\left(3 \cdot x1\right) \cdot x1\right) \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 + 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) + \left(x1 \cdot x1\right) \cdot x1\right) + x1\right) + 3 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 - 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \]
   2. Add Preprocessing

   4. Applied rewrites 36.9%

      \[\leadsto x1 + \color{blue}{\mathsf{fma}\left(x1 \cdot x1, x1, \mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(4, \frac{\mathsf{fma}\left(x2, 2, \left(3 \cdot x1\right) \cdot x1\right) - x1}{\mathsf{fma}\left(x1, x1, 1\right)}, -6\right), x1 \cdot x1, \left(\frac{\mathsf{fma}\left(x2, 2, \left(3 \cdot x1\right) \cdot x1\right) - x1}{\mathsf{fma}\left(x1, x1, 1\right)} - 3\right) \cdot \left(\frac{\mathsf{fma}\left(x2, 2, \left(3 \cdot x1\right) \cdot x1\right) - x1}{\mathsf{fma}\left(x1, x1, 1\right)} \cdot \left(x1 \cdot 2\right)\right)\right), \mathsf{fma}\left(x1, x1, 1\right), \frac{\mathsf{fma}\left(x2, 2, \left(3 \cdot x1\right) \cdot x1\right) - x1}{\mathsf{fma}\left(x1, x1, 1\right)} \cdot \left(\left(3 \cdot x1\right) \cdot x1\right)\right) + \mathsf{fma}\left(\frac{\mathsf{fma}\left(-2, x2, \left(3 \cdot x1\right) \cdot x1\right) - x1}{\mathsf{fma}\left(x1, x1, 1\right)}, 3, x1\right)\right)} \]

   5. Taylor expanded in x1 around 0

      \[\leadsto x1 + \mathsf{fma}\left(x1 \cdot x1, x1, \color{blue}{-6 \cdot x2}\right) \]
   6. Step-by-step derivation

      1. lower-*.f64 1.7

         \[\leadsto x1 + \mathsf{fma}\left(x1 \cdot x1, x1, \color{blue}{-6 \cdot x2}\right) \]

   7. Applied rewrites 1.7%

      \[\leadsto x1 + \mathsf{fma}\left(x1 \cdot x1, x1, \color{blue}{-6 \cdot x2}\right) \]
   8. Step-by-step derivation

      1. lift-+.f64 N/A

         \[\leadsto \color{blue}{x1 + \mathsf{fma}\left(x1 \cdot x1, x1, -6 \cdot x2\right)} \]

   9. Applied rewrites 1.7%

      \[\leadsto \color{blue}{\mathsf{fma}\left(\mathsf{fma}\left(x1, x1, 1\right), x1, -6 \cdot x2\right)} \]

   10. Taylor expanded in x1 around -inf

       \[\leadsto \color{blue}{{x1}^{4} \cdot \left(6 + -1 \cdot \frac{3 + -1 \cdot \frac{9 + 4 \cdot \left(2 \cdot x2 - 3\right)}{x1}}{x1}\right)} \]
   11. Step-by-step derivation

       1. *-commutative N/A

          \[\leadsto \color{blue}{\left(6 + -1 \cdot \frac{3 + -1 \cdot \frac{9 + 4 \cdot \left(2 \cdot x2 - 3\right)}{x1}}{x1}\right) \cdot {x1}^{4}} \]

       2. lower-*.f64 N/A

          \[\leadsto \color{blue}{\left(6 + -1 \cdot \frac{3 + -1 \cdot \frac{9 + 4 \cdot \left(2 \cdot x2 - 3\right)}{x1}}{x1}\right) \cdot {x1}^{4}} \]

   12. Applied rewrites 97.0%

       \[\leadsto \color{blue}{\left(6 - \frac{3 - \frac{\mathsf{fma}\left(\mathsf{fma}\left(2, x2, -3\right), 4, 9\right)}{x1}}{x1}\right) \cdot {x1}^{4}} \]

   if -2.8e32 < x1 < 1.6499999999999999e40

   1. Initial program 99.3%

      \[x1 + \left(\left(\left(\left(\left(\left(\left(2 \cdot x1\right) \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 + 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \cdot \left(\frac{\left(\left(3 \cdot x1\right) \cdot x1 + 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1} - 3\right) + \left(x1 \cdot x1\right) \cdot \left(4 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 + 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1} - 6\right)\right) \cdot \left(x1 \cdot x1 + 1\right) + \left(\left(3 \cdot x1\right) \cdot x1\right) \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 + 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) + \left(x1 \cdot x1\right) \cdot x1\right) + x1\right) + 3 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 - 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \]
   2. Add Preprocessing

   3. Taylor expanded in x2 around inf

      \[\leadsto x1 + \left(\left(\color{blue}{8 \cdot \frac{x1 \cdot {x2}^{2}}{1 + {x1}^{2}}} + x1\right) + 3 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 - 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \]
   4. Step-by-step derivation

      1. associate-*r/ N/A

         \[\leadsto x1 + \left(\left(\color{blue}{\frac{8 \cdot \left(x1 \cdot {x2}^{2}\right)}{1 + {x1}^{2}}} + x1\right) + 3 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 - 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \]

      2. associate-*r* N/A

         \[\leadsto x1 + \left(\left(\frac{\color{blue}{\left(8 \cdot x1\right) \cdot {x2}^{2}}}{1 + {x1}^{2}} + x1\right) + 3 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 - 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \]

      3. associate-*l/ N/A

         \[\leadsto x1 + \left(\left(\color{blue}{\frac{8 \cdot x1}{1 + {x1}^{2}} \cdot {x2}^{2}} + x1\right) + 3 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 - 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \]

      4. associate-*r/ N/A

         \[\leadsto x1 + \left(\left(\color{blue}{\left(8 \cdot \frac{x1}{1 + {x1}^{2}}\right)} \cdot {x2}^{2} + x1\right) + 3 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 - 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \]

      5. unpow2 N/A

         \[\leadsto x1 + \left(\left(\left(8 \cdot \frac{x1}{1 + {x1}^{2}}\right) \cdot \color{blue}{\left(x2 \cdot x2\right)} + x1\right) + 3 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 - 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \]

      6. associate-*r* N/A

         \[\leadsto x1 + \left(\left(\color{blue}{\left(\left(8 \cdot \frac{x1}{1 + {x1}^{2}}\right) \cdot x2\right) \cdot x2} + x1\right) + 3 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 - 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \]

      7. lower-*.f64 N/A

         \[\leadsto x1 + \left(\left(\color{blue}{\left(\left(8 \cdot \frac{x1}{1 + {x1}^{2}}\right) \cdot x2\right) \cdot x2} + x1\right) + 3 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 - 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \]

      8. lower-*.f64 N/A

         \[\leadsto x1 + \left(\left(\color{blue}{\left(\left(8 \cdot \frac{x1}{1 + {x1}^{2}}\right) \cdot x2\right)} \cdot x2 + x1\right) + 3 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 - 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \]

      9. *-commutative N/A

         \[\leadsto x1 + \left(\left(\left(\color{blue}{\left(\frac{x1}{1 + {x1}^{2}} \cdot 8\right)} \cdot x2\right) \cdot x2 + x1\right) + 3 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 - 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \]

      10. lower-*.f64 N/A

          \[\leadsto x1 + \left(\left(\left(\color{blue}{\left(\frac{x1}{1 + {x1}^{2}} \cdot 8\right)} \cdot x2\right) \cdot x2 + x1\right) + 3 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 - 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \]

      11. lower-/.f64 N/A

          \[\leadsto x1 + \left(\left(\left(\left(\color{blue}{\frac{x1}{1 + {x1}^{2}}} \cdot 8\right) \cdot x2\right) \cdot x2 + x1\right) + 3 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 - 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \]

      12. +-commutative N/A

          \[\leadsto x1 + \left(\left(\left(\left(\frac{x1}{\color{blue}{{x1}^{2} + 1}} \cdot 8\right) \cdot x2\right) \cdot x2 + x1\right) + 3 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 - 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \]

      13. unpow2 N/A

          \[\leadsto x1 + \left(\left(\left(\left(\frac{x1}{\color{blue}{x1 \cdot x1} + 1} \cdot 8\right) \cdot x2\right) \cdot x2 + x1\right) + 3 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 - 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \]

      14. lower-fma.f64 97.6

          \[\leadsto x1 + \left(\left(\left(\left(\frac{x1}{\color{blue}{\mathsf{fma}\left(x1, x1, 1\right)}} \cdot 8\right) \cdot x2\right) \cdot x2 + x1\right) + 3 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 - 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \]

   5. Applied rewrites 97.6%

      \[\leadsto x1 + \left(\left(\color{blue}{\left(\left(\frac{x1}{\mathsf{fma}\left(x1, x1, 1\right)} \cdot 8\right) \cdot x2\right) \cdot x2} + x1\right) + 3 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 - 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \]

   if 1.6499999999999999e40 < x1

   1. Initial program 39.9%

      \[x1 + \left(\left(\left(\left(\left(\left(\left(2 \cdot x1\right) \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 + 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \cdot \left(\frac{\left(\left(3 \cdot x1\right) \cdot x1 + 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1} - 3\right) + \left(x1 \cdot x1\right) \cdot \left(4 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 + 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1} - 6\right)\right) \cdot \left(x1 \cdot x1 + 1\right) + \left(\left(3 \cdot x1\right) \cdot x1\right) \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 + 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) + \left(x1 \cdot x1\right) \cdot x1\right) + x1\right) + 3 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 - 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \]
   2. Add Preprocessing

   4. Applied rewrites 39.9%

      \[\leadsto x1 + \color{blue}{\mathsf{fma}\left(x1 \cdot x1, x1, \mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(4, \frac{\mathsf{fma}\left(x2, 2, \left(3 \cdot x1\right) \cdot x1\right) - x1}{\mathsf{fma}\left(x1, x1, 1\right)}, -6\right), x1 \cdot x1, \left(\frac{\mathsf{fma}\left(x2, 2, \left(3 \cdot x1\right) \cdot x1\right) - x1}{\mathsf{fma}\left(x1, x1, 1\right)} - 3\right) \cdot \left(\frac{\mathsf{fma}\left(x2, 2, \left(3 \cdot x1\right) \cdot x1\right) - x1}{\mathsf{fma}\left(x1, x1, 1\right)} \cdot \left(x1 \cdot 2\right)\right)\right), \mathsf{fma}\left(x1, x1, 1\right), \frac{\mathsf{fma}\left(x2, 2, \left(3 \cdot x1\right) \cdot x1\right) - x1}{\mathsf{fma}\left(x1, x1, 1\right)} \cdot \left(\left(3 \cdot x1\right) \cdot x1\right)\right) + \mathsf{fma}\left(\frac{\mathsf{fma}\left(-2, x2, \left(3 \cdot x1\right) \cdot x1\right) - x1}{\mathsf{fma}\left(x1, x1, 1\right)}, 3, x1\right)\right)} \]

   5. Taylor expanded in x1 around -inf

      \[\leadsto x1 + \color{blue}{{x1}^{4} \cdot \left(6 + -1 \cdot \frac{3 + -1 \cdot \frac{9 + 4 \cdot \left(2 \cdot x2 - 3\right)}{x1}}{x1}\right)} \]
   6. Step-by-step derivation

      1. *-commutative N/A

         \[\leadsto x1 + \color{blue}{\left(6 + -1 \cdot \frac{3 + -1 \cdot \frac{9 + 4 \cdot \left(2 \cdot x2 - 3\right)}{x1}}{x1}\right) \cdot {x1}^{4}} \]

      2. lower-*.f64 N/A

         \[\leadsto x1 + \color{blue}{\left(6 + -1 \cdot \frac{3 + -1 \cdot \frac{9 + 4 \cdot \left(2 \cdot x2 - 3\right)}{x1}}{x1}\right) \cdot {x1}^{4}} \]

   7. Applied rewrites 98.1%

      \[\leadsto x1 + \color{blue}{\left(6 - \frac{3 - \frac{\mathsf{fma}\left(\mathsf{fma}\left(2, x2, -3\right), 4, 9\right)}{x1}}{x1}\right) \cdot {x1}^{4}} \]
3. Recombined 3 regimes into one program.

4. Final simplification 97.6%

   \[\leadsto \begin{array}{l} \mathbf{if}\;x1 \leq -2.8 \cdot 10^{+32}:\\ \;\;\;\;\left(6 - \frac{3 - \frac{\mathsf{fma}\left(\mathsf{fma}\left(2, x2, -3\right), 4, 9\right)}{x1}}{x1}\right) \cdot {x1}^{4}\\ \mathbf{elif}\;x1 \leq 1.65 \cdot 10^{+40}:\\ \;\;\;\;x1 + \left(\left(\left(\left(\frac{x1}{\mathsf{fma}\left(x1, x1, 1\right)} \cdot 8\right) \cdot x2\right) \cdot x2 + x1\right) + 3 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 - 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right)\\ \mathbf{else}:\\ \;\;\;\;x1 + \left(6 - \frac{3 - \frac{\mathsf{fma}\left(\mathsf{fma}\left(2, x2, -3\right), 4, 9\right)}{x1}}{x1}\right) \cdot {x1}^{4}\\ \end{array} \]
5. Add Preprocessing

Alternative 3: 94.5% accurate, 1.9× speedup

Math

\[\begin{array}{l} \\ \begin{array}{l} \mathbf{if}\;x1 \leq -2.8 \cdot 10^{+32} \lor \neg \left(x1 \leq 1.65 \cdot 10^{+40}\right):\\ \;\;\;\;\left(6 - \frac{3 - \frac{\mathsf{fma}\left(\mathsf{fma}\left(2, x2, -3\right), 4, 9\right)}{x1}}{x1}\right) \cdot {x1}^{4}\\ \mathbf{else}:\\ \;\;\;\;x1 + \left(\left(\left(\left(\frac{x1}{\mathsf{fma}\left(x1, x1, 1\right)} \cdot 8\right) \cdot x2\right) \cdot x2 + x1\right) + 3 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 - 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right)\\ \end{array} \end{array} \]

FPCore

(FPCore (x1 x2)
 :precision binary64
 (if (or (<= x1 -2.8e+32) (not (<= x1 1.65e+40)))
   (*
    (- 6.0 (/ (- 3.0 (/ (fma (fma 2.0 x2 -3.0) 4.0 9.0) x1)) x1))
    (pow x1 4.0))
   (+
    x1
    (+
     (+ (* (* (* (/ x1 (fma x1 x1 1.0)) 8.0) x2) x2) x1)
     (* 3.0 (/ (- (- (* (* 3.0 x1) x1) (* 2.0 x2)) x1) (+ (* x1 x1) 1.0)))))))

C

double code(double x1, double x2) {
	double tmp;
	if ((x1 <= -2.8e+32) || !(x1 <= 1.65e+40)) {
		tmp = (6.0 - ((3.0 - (fma(fma(2.0, x2, -3.0), 4.0, 9.0) / x1)) / x1)) * pow(x1, 4.0);
	} else {
		tmp = x1 + ((((((x1 / fma(x1, x1, 1.0)) * 8.0) * x2) * x2) + x1) + (3.0 * (((((3.0 * x1) * x1) - (2.0 * x2)) - x1) / ((x1 * x1) + 1.0))));
	}
	return tmp;
}

Julia

function code(x1, x2)
	tmp = 0.0
	if ((x1 <= -2.8e+32) || !(x1 <= 1.65e+40))
		tmp = Float64(Float64(6.0 - Float64(Float64(3.0 - Float64(fma(fma(2.0, x2, -3.0), 4.0, 9.0) / x1)) / x1)) * (x1 ^ 4.0));
	else
		tmp = Float64(x1 + Float64(Float64(Float64(Float64(Float64(Float64(x1 / fma(x1, x1, 1.0)) * 8.0) * x2) * x2) + x1) + Float64(3.0 * Float64(Float64(Float64(Float64(Float64(3.0 * x1) * x1) - Float64(2.0 * x2)) - x1) / Float64(Float64(x1 * x1) + 1.0)))));
	end
	return tmp
end

Wolfram

code[x1_, x2_] := If[Or[LessEqual[x1, -2.8e+32], N[Not[LessEqual[x1, 1.65e+40]], $MachinePrecision]], N[(N[(6.0 - N[(N[(3.0 - N[(N[(N[(2.0 * x2 + -3.0), $MachinePrecision] * 4.0 + 9.0), $MachinePrecision] / x1), $MachinePrecision]), $MachinePrecision] / x1), $MachinePrecision]), $MachinePrecision] * N[Power[x1, 4.0], $MachinePrecision]), $MachinePrecision], N[(x1 + N[(N[(N[(N[(N[(N[(x1 / N[(x1 * x1 + 1.0), $MachinePrecision]), $MachinePrecision] * 8.0), $MachinePrecision] * x2), $MachinePrecision] * x2), $MachinePrecision] + x1), $MachinePrecision] + N[(3.0 * N[(N[(N[(N[(N[(3.0 * x1), $MachinePrecision] * x1), $MachinePrecision] - N[(2.0 * x2), $MachinePrecision]), $MachinePrecision] - x1), $MachinePrecision] / N[(N[(x1 * x1), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]

Derivation

1. Split input into 2 regimes

2. if x1 < -2.8e32 or 1.6499999999999999e40 < x1

   1. Initial program 29.9%

      \[x1 + \left(\left(\left(\left(\left(\left(\left(2 \cdot x1\right) \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 + 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \cdot \left(\frac{\left(\left(3 \cdot x1\right) \cdot x1 + 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1} - 3\right) + \left(x1 \cdot x1\right) \cdot \left(4 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 + 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1} - 6\right)\right) \cdot \left(x1 \cdot x1 + 1\right) + \left(\left(3 \cdot x1\right) \cdot x1\right) \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 + 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) + \left(x1 \cdot x1\right) \cdot x1\right) + x1\right) + 3 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 - 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \]
   2. Add Preprocessing

   4. Applied rewrites 38.2%

      \[\leadsto x1 + \color{blue}{\mathsf{fma}\left(x1 \cdot x1, x1, \mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(4, \frac{\mathsf{fma}\left(x2, 2, \left(3 \cdot x1\right) \cdot x1\right) - x1}{\mathsf{fma}\left(x1, x1, 1\right)}, -6\right), x1 \cdot x1, \left(\frac{\mathsf{fma}\left(x2, 2, \left(3 \cdot x1\right) \cdot x1\right) - x1}{\mathsf{fma}\left(x1, x1, 1\right)} - 3\right) \cdot \left(\frac{\mathsf{fma}\left(x2, 2, \left(3 \cdot x1\right) \cdot x1\right) - x1}{\mathsf{fma}\left(x1, x1, 1\right)} \cdot \left(x1 \cdot 2\right)\right)\right), \mathsf{fma}\left(x1, x1, 1\right), \frac{\mathsf{fma}\left(x2, 2, \left(3 \cdot x1\right) \cdot x1\right) - x1}{\mathsf{fma}\left(x1, x1, 1\right)} \cdot \left(\left(3 \cdot x1\right) \cdot x1\right)\right) + \mathsf{fma}\left(\frac{\mathsf{fma}\left(-2, x2, \left(3 \cdot x1\right) \cdot x1\right) - x1}{\mathsf{fma}\left(x1, x1, 1\right)}, 3, x1\right)\right)} \]

   5. Taylor expanded in x1 around 0

      \[\leadsto x1 + \mathsf{fma}\left(x1 \cdot x1, x1, \color{blue}{-6 \cdot x2}\right) \]
   6. Step-by-step derivation

      1. lower-*.f64 36.7

         \[\leadsto x1 + \mathsf{fma}\left(x1 \cdot x1, x1, \color{blue}{-6 \cdot x2}\right) \]

   7. Applied rewrites 36.7%

      \[\leadsto x1 + \mathsf{fma}\left(x1 \cdot x1, x1, \color{blue}{-6 \cdot x2}\right) \]
   8. Step-by-step derivation

      1. lift-+.f64 N/A

         \[\leadsto \color{blue}{x1 + \mathsf{fma}\left(x1 \cdot x1, x1, -6 \cdot x2\right)} \]

   9. Applied rewrites 36.7%

      \[\leadsto \color{blue}{\mathsf{fma}\left(\mathsf{fma}\left(x1, x1, 1\right), x1, -6 \cdot x2\right)} \]

   10. Taylor expanded in x1 around -inf

       \[\leadsto \color{blue}{{x1}^{4} \cdot \left(6 + -1 \cdot \frac{3 + -1 \cdot \frac{9 + 4 \cdot \left(2 \cdot x2 - 3\right)}{x1}}{x1}\right)} \]
   11. Step-by-step derivation

       1. *-commutative N/A

          \[\leadsto \color{blue}{\left(6 + -1 \cdot \frac{3 + -1 \cdot \frac{9 + 4 \cdot \left(2 \cdot x2 - 3\right)}{x1}}{x1}\right) \cdot {x1}^{4}} \]

       2. lower-*.f64 N/A

          \[\leadsto \color{blue}{\left(6 + -1 \cdot \frac{3 + -1 \cdot \frac{9 + 4 \cdot \left(2 \cdot x2 - 3\right)}{x1}}{x1}\right) \cdot {x1}^{4}} \]

   12. Applied rewrites 97.5%

       \[\leadsto \color{blue}{\left(6 - \frac{3 - \frac{\mathsf{fma}\left(\mathsf{fma}\left(2, x2, -3\right), 4, 9\right)}{x1}}{x1}\right) \cdot {x1}^{4}} \]

   if -2.8e32 < x1 < 1.6499999999999999e40

   1. Initial program 99.3%

      \[x1 + \left(\left(\left(\left(\left(\left(\left(2 \cdot x1\right) \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 + 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \cdot \left(\frac{\left(\left(3 \cdot x1\right) \cdot x1 + 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1} - 3\right) + \left(x1 \cdot x1\right) \cdot \left(4 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 + 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1} - 6\right)\right) \cdot \left(x1 \cdot x1 + 1\right) + \left(\left(3 \cdot x1\right) \cdot x1\right) \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 + 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) + \left(x1 \cdot x1\right) \cdot x1\right) + x1\right) + 3 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 - 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \]
   2. Add Preprocessing

   3. Taylor expanded in x2 around inf

      \[\leadsto x1 + \left(\left(\color{blue}{8 \cdot \frac{x1 \cdot {x2}^{2}}{1 + {x1}^{2}}} + x1\right) + 3 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 - 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \]
   4. Step-by-step derivation

      1. associate-*r/ N/A

         \[\leadsto x1 + \left(\left(\color{blue}{\frac{8 \cdot \left(x1 \cdot {x2}^{2}\right)}{1 + {x1}^{2}}} + x1\right) + 3 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 - 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \]

      2. associate-*r* N/A

         \[\leadsto x1 + \left(\left(\frac{\color{blue}{\left(8 \cdot x1\right) \cdot {x2}^{2}}}{1 + {x1}^{2}} + x1\right) + 3 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 - 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \]

      3. associate-*l/ N/A

         \[\leadsto x1 + \left(\left(\color{blue}{\frac{8 \cdot x1}{1 + {x1}^{2}} \cdot {x2}^{2}} + x1\right) + 3 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 - 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \]

      4. associate-*r/ N/A

         \[\leadsto x1 + \left(\left(\color{blue}{\left(8 \cdot \frac{x1}{1 + {x1}^{2}}\right)} \cdot {x2}^{2} + x1\right) + 3 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 - 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \]

      5. unpow2 N/A

         \[\leadsto x1 + \left(\left(\left(8 \cdot \frac{x1}{1 + {x1}^{2}}\right) \cdot \color{blue}{\left(x2 \cdot x2\right)} + x1\right) + 3 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 - 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \]

      6. associate-*r* N/A

         \[\leadsto x1 + \left(\left(\color{blue}{\left(\left(8 \cdot \frac{x1}{1 + {x1}^{2}}\right) \cdot x2\right) \cdot x2} + x1\right) + 3 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 - 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \]

      7. lower-*.f64 N/A

         \[\leadsto x1 + \left(\left(\color{blue}{\left(\left(8 \cdot \frac{x1}{1 + {x1}^{2}}\right) \cdot x2\right) \cdot x2} + x1\right) + 3 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 - 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \]

      8. lower-*.f64 N/A

         \[\leadsto x1 + \left(\left(\color{blue}{\left(\left(8 \cdot \frac{x1}{1 + {x1}^{2}}\right) \cdot x2\right)} \cdot x2 + x1\right) + 3 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 - 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \]

      9. *-commutative N/A

         \[\leadsto x1 + \left(\left(\left(\color{blue}{\left(\frac{x1}{1 + {x1}^{2}} \cdot 8\right)} \cdot x2\right) \cdot x2 + x1\right) + 3 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 - 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \]

      10. lower-*.f64 N/A

          \[\leadsto x1 + \left(\left(\left(\color{blue}{\left(\frac{x1}{1 + {x1}^{2}} \cdot 8\right)} \cdot x2\right) \cdot x2 + x1\right) + 3 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 - 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \]

      11. lower-/.f64 N/A

          \[\leadsto x1 + \left(\left(\left(\left(\color{blue}{\frac{x1}{1 + {x1}^{2}}} \cdot 8\right) \cdot x2\right) \cdot x2 + x1\right) + 3 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 - 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \]

      12. +-commutative N/A

          \[\leadsto x1 + \left(\left(\left(\left(\frac{x1}{\color{blue}{{x1}^{2} + 1}} \cdot 8\right) \cdot x2\right) \cdot x2 + x1\right) + 3 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 - 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \]

      13. unpow2 N/A

          \[\leadsto x1 + \left(\left(\left(\left(\frac{x1}{\color{blue}{x1 \cdot x1} + 1} \cdot 8\right) \cdot x2\right) \cdot x2 + x1\right) + 3 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 - 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \]

      14. lower-fma.f64 97.6

          \[\leadsto x1 + \left(\left(\left(\left(\frac{x1}{\color{blue}{\mathsf{fma}\left(x1, x1, 1\right)}} \cdot 8\right) \cdot x2\right) \cdot x2 + x1\right) + 3 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 - 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \]

   5. Applied rewrites 97.6%

      \[\leadsto x1 + \left(\left(\color{blue}{\left(\left(\frac{x1}{\mathsf{fma}\left(x1, x1, 1\right)} \cdot 8\right) \cdot x2\right) \cdot x2} + x1\right) + 3 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 - 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \]
3. Recombined 2 regimes into one program.

4. Final simplification 97.6%

   \[\leadsto \begin{array}{l} \mathbf{if}\;x1 \leq -2.8 \cdot 10^{+32} \lor \neg \left(x1 \leq 1.65 \cdot 10^{+40}\right):\\ \;\;\;\;\left(6 - \frac{3 - \frac{\mathsf{fma}\left(\mathsf{fma}\left(2, x2, -3\right), 4, 9\right)}{x1}}{x1}\right) \cdot {x1}^{4}\\ \mathbf{else}:\\ \;\;\;\;x1 + \left(\left(\left(\left(\frac{x1}{\mathsf{fma}\left(x1, x1, 1\right)} \cdot 8\right) \cdot x2\right) \cdot x2 + x1\right) + 3 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 - 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right)\\ \end{array} \]
5. Add Preprocessing

Alternative 4: 66.5% accurate, 2.3× speedup

Math

\[\begin{array}{l} \\ \begin{array}{l} \mathbf{if}\;x1 \leq -7.5 \cdot 10^{+137}:\\ \;\;\;\;\left(\left(x2 \cdot x2\right) \cdot -36\right) \cdot {\left(x1 - -6 \cdot x2\right)}^{-1}\\ \mathbf{elif}\;x1 \leq 3.5 \cdot 10^{+88}:\\ \;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(2, x2, -3\right) \cdot x2, 4, -1\right), x1, -6 \cdot x2\right)\\ \mathbf{else}:\\ \;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(x1, x1, 1\right), x1, -6 \cdot x2\right)\\ \end{array} \end{array} \]

FPCore

(FPCore (x1 x2)
 :precision binary64
 (if (<= x1 -7.5e+137)
   (* (* (* x2 x2) -36.0) (pow (- x1 (* -6.0 x2)) -1.0))
   (if (<= x1 3.5e+88)
     (fma (fma (* (fma 2.0 x2 -3.0) x2) 4.0 -1.0) x1 (* -6.0 x2))
     (fma (fma x1 x1 1.0) x1 (* -6.0 x2)))))

C

double code(double x1, double x2) {
	double tmp;
	if (x1 <= -7.5e+137) {
		tmp = ((x2 * x2) * -36.0) * pow((x1 - (-6.0 * x2)), -1.0);
	} else if (x1 <= 3.5e+88) {
		tmp = fma(fma((fma(2.0, x2, -3.0) * x2), 4.0, -1.0), x1, (-6.0 * x2));
	} else {
		tmp = fma(fma(x1, x1, 1.0), x1, (-6.0 * x2));
	}
	return tmp;
}

Julia

function code(x1, x2)
	tmp = 0.0
	if (x1 <= -7.5e+137)
		tmp = Float64(Float64(Float64(x2 * x2) * -36.0) * (Float64(x1 - Float64(-6.0 * x2)) ^ -1.0));
	elseif (x1 <= 3.5e+88)
		tmp = fma(fma(Float64(fma(2.0, x2, -3.0) * x2), 4.0, -1.0), x1, Float64(-6.0 * x2));
	else
		tmp = fma(fma(x1, x1, 1.0), x1, Float64(-6.0 * x2));
	end
	return tmp
end

Wolfram

code[x1_, x2_] := If[LessEqual[x1, -7.5e+137], N[(N[(N[(x2 * x2), $MachinePrecision] * -36.0), $MachinePrecision] * N[Power[N[(x1 - N[(-6.0 * x2), $MachinePrecision]), $MachinePrecision], -1.0], $MachinePrecision]), $MachinePrecision], If[LessEqual[x1, 3.5e+88], N[(N[(N[(N[(2.0 * x2 + -3.0), $MachinePrecision] * x2), $MachinePrecision] * 4.0 + -1.0), $MachinePrecision] * x1 + N[(-6.0 * x2), $MachinePrecision]), $MachinePrecision], N[(N[(x1 * x1 + 1.0), $MachinePrecision] * x1 + N[(-6.0 * x2), $MachinePrecision]), $MachinePrecision]]]

Derivation

1. Split input into 3 regimes

2. if x1 < -7.50000000000000025e137

   1. Initial program 0.0%

      \[x1 + \left(\left(\left(\left(\left(\left(\left(2 \cdot x1\right) \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 + 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \cdot \left(\frac{\left(\left(3 \cdot x1\right) \cdot x1 + 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1} - 3\right) + \left(x1 \cdot x1\right) \cdot \left(4 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 + 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1} - 6\right)\right) \cdot \left(x1 \cdot x1 + 1\right) + \left(\left(3 \cdot x1\right) \cdot x1\right) \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 + 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) + \left(x1 \cdot x1\right) \cdot x1\right) + x1\right) + 3 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 - 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \]
   2. Add Preprocessing

   3. Taylor expanded in x1 around 0

      \[\leadsto x1 + \color{blue}{-6 \cdot x2} \]
   4. Step-by-step derivation

      1. lower-*.f64 0.5

         \[\leadsto x1 + \color{blue}{-6 \cdot x2} \]

   5. Applied rewrites 0.5%

      \[\leadsto x1 + \color{blue}{-6 \cdot x2} \]
   6. Step-by-step derivation

      1. lift-+.f64 N/A

         \[\leadsto \color{blue}{x1 + -6 \cdot x2} \]

   7. Applied rewrites 2.4%

      \[\leadsto \color{blue}{\left(x1 \cdot x1 - {\left(-6 \cdot x2\right)}^{2}\right) \cdot \frac{1}{x1 - -6 \cdot x2}} \]

   8. Taylor expanded in x1 around 0

      \[\leadsto \color{blue}{\left(-36 \cdot {x2}^{2}\right)} \cdot \frac{1}{x1 - -6 \cdot x2} \]
   9. Step-by-step derivation

      1. *-commutative N/A

         \[\leadsto \color{blue}{\left({x2}^{2} \cdot -36\right)} \cdot \frac{1}{x1 - -6 \cdot x2} \]

      2. lower-*.f64 N/A

         \[\leadsto \color{blue}{\left({x2}^{2} \cdot -36\right)} \cdot \frac{1}{x1 - -6 \cdot x2} \]

      3. unpow2 N/A

         \[\leadsto \left(\color{blue}{\left(x2 \cdot x2\right)} \cdot -36\right) \cdot \frac{1}{x1 - -6 \cdot x2} \]

      4. lower-*.f64 13.5

         \[\leadsto \left(\color{blue}{\left(x2 \cdot x2\right)} \cdot -36\right) \cdot \frac{1}{x1 - -6 \cdot x2} \]

   10. Applied rewrites 13.5%

       \[\leadsto \color{blue}{\left(\left(x2 \cdot x2\right) \cdot -36\right)} \cdot \frac{1}{x1 - -6 \cdot x2} \]

   if -7.50000000000000025e137 < x1 < 3.4999999999999998e88

   1. Initial program 94.0%

      \[x1 + \left(\left(\left(\left(\left(\left(\left(2 \cdot x1\right) \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 + 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \cdot \left(\frac{\left(\left(3 \cdot x1\right) \cdot x1 + 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1} - 3\right) + \left(x1 \cdot x1\right) \cdot \left(4 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 + 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1} - 6\right)\right) \cdot \left(x1 \cdot x1 + 1\right) + \left(\left(3 \cdot x1\right) \cdot x1\right) \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 + 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) + \left(x1 \cdot x1\right) \cdot x1\right) + x1\right) + 3 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 - 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \]
   2. Add Preprocessing

   4. Applied rewrites 99.7%

      \[\leadsto x1 + \color{blue}{\mathsf{fma}\left(x1 \cdot x1, x1, \mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(4, \frac{\mathsf{fma}\left(x2, 2, \left(3 \cdot x1\right) \cdot x1\right) - x1}{\mathsf{fma}\left(x1, x1, 1\right)}, -6\right), x1 \cdot x1, \left(\frac{\mathsf{fma}\left(x2, 2, \left(3 \cdot x1\right) \cdot x1\right) - x1}{\mathsf{fma}\left(x1, x1, 1\right)} - 3\right) \cdot \left(\frac{\mathsf{fma}\left(x2, 2, \left(3 \cdot x1\right) \cdot x1\right) - x1}{\mathsf{fma}\left(x1, x1, 1\right)} \cdot \left(x1 \cdot 2\right)\right)\right), \mathsf{fma}\left(x1, x1, 1\right), \frac{\mathsf{fma}\left(x2, 2, \left(3 \cdot x1\right) \cdot x1\right) - x1}{\mathsf{fma}\left(x1, x1, 1\right)} \cdot \left(\left(3 \cdot x1\right) \cdot x1\right)\right) + \mathsf{fma}\left(\frac{\mathsf{fma}\left(-2, x2, \left(3 \cdot x1\right) \cdot x1\right) - x1}{\mathsf{fma}\left(x1, x1, 1\right)}, 3, x1\right)\right)} \]

   5. Taylor expanded in x1 around 0

      \[\leadsto \color{blue}{-6 \cdot x2 + x1 \cdot \left(4 \cdot \left(x2 \cdot \left(2 \cdot x2 - 3\right)\right) - 1\right)} \]
   6. Step-by-step derivation

      1. +-commutative N/A

         \[\leadsto \color{blue}{x1 \cdot \left(4 \cdot \left(x2 \cdot \left(2 \cdot x2 - 3\right)\right) - 1\right) + -6 \cdot x2} \]

      2. *-commutative N/A

         \[\leadsto \color{blue}{\left(4 \cdot \left(x2 \cdot \left(2 \cdot x2 - 3\right)\right) - 1\right) \cdot x1} + -6 \cdot x2 \]

      3. lower-fma.f64 N/A

         \[\leadsto \color{blue}{\mathsf{fma}\left(4 \cdot \left(x2 \cdot \left(2 \cdot x2 - 3\right)\right) - 1, x1, -6 \cdot x2\right)} \]

      4. sub-neg N/A

         \[\leadsto \mathsf{fma}\left(\color{blue}{4 \cdot \left(x2 \cdot \left(2 \cdot x2 - 3\right)\right) + \left(\mathsf{neg}\left(1\right)\right)}, x1, -6 \cdot x2\right) \]

      5. *-commutative N/A

         \[\leadsto \mathsf{fma}\left(\color{blue}{\left(x2 \cdot \left(2 \cdot x2 - 3\right)\right) \cdot 4} + \left(\mathsf{neg}\left(1\right)\right), x1, -6 \cdot x2\right) \]

      6. metadata-eval N/A

         \[\leadsto \mathsf{fma}\left(\left(x2 \cdot \left(2 \cdot x2 - 3\right)\right) \cdot 4 + \color{blue}{-1}, x1, -6 \cdot x2\right) \]

      7. lower-fma.f64 N/A

         \[\leadsto \mathsf{fma}\left(\color{blue}{\mathsf{fma}\left(x2 \cdot \left(2 \cdot x2 - 3\right), 4, -1\right)}, x1, -6 \cdot x2\right) \]

      8. *-commutative N/A

         \[\leadsto \mathsf{fma}\left(\mathsf{fma}\left(\color{blue}{\left(2 \cdot x2 - 3\right) \cdot x2}, 4, -1\right), x1, -6 \cdot x2\right) \]

      9. lower-*.f64 N/A

         \[\leadsto \mathsf{fma}\left(\mathsf{fma}\left(\color{blue}{\left(2 \cdot x2 - 3\right) \cdot x2}, 4, -1\right), x1, -6 \cdot x2\right) \]

      10. sub-neg N/A

          \[\leadsto \mathsf{fma}\left(\mathsf{fma}\left(\color{blue}{\left(2 \cdot x2 + \left(\mathsf{neg}\left(3\right)\right)\right)} \cdot x2, 4, -1\right), x1, -6 \cdot x2\right) \]

      11. lower-fma.f64 N/A

          \[\leadsto \mathsf{fma}\left(\mathsf{fma}\left(\color{blue}{\mathsf{fma}\left(2, x2, \mathsf{neg}\left(3\right)\right)} \cdot x2, 4, -1\right), x1, -6 \cdot x2\right) \]

      12. metadata-eval N/A

          \[\leadsto \mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(2, x2, \color{blue}{-3}\right) \cdot x2, 4, -1\right), x1, -6 \cdot x2\right) \]

      13. lower-*.f64 72.3

          \[\leadsto \mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(2, x2, -3\right) \cdot x2, 4, -1\right), x1, \color{blue}{-6 \cdot x2}\right) \]

   7. Applied rewrites 72.3%

      \[\leadsto \color{blue}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(2, x2, -3\right) \cdot x2, 4, -1\right), x1, -6 \cdot x2\right)} \]

   if 3.4999999999999998e88 < x1

   1. Initial program 28.3%

      \[x1 + \left(\left(\left(\left(\left(\left(\left(2 \cdot x1\right) \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 + 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \cdot \left(\frac{\left(\left(3 \cdot x1\right) \cdot x1 + 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1} - 3\right) + \left(x1 \cdot x1\right) \cdot \left(4 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 + 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1} - 6\right)\right) \cdot \left(x1 \cdot x1 + 1\right) + \left(\left(3 \cdot x1\right) \cdot x1\right) \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 + 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) + \left(x1 \cdot x1\right) \cdot x1\right) + x1\right) + 3 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 - 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \]
   2. Add Preprocessing

   4. Applied rewrites 28.3%

      \[\leadsto x1 + \color{blue}{\mathsf{fma}\left(x1 \cdot x1, x1, \mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(4, \frac{\mathsf{fma}\left(x2, 2, \left(3 \cdot x1\right) \cdot x1\right) - x1}{\mathsf{fma}\left(x1, x1, 1\right)}, -6\right), x1 \cdot x1, \left(\frac{\mathsf{fma}\left(x2, 2, \left(3 \cdot x1\right) \cdot x1\right) - x1}{\mathsf{fma}\left(x1, x1, 1\right)} - 3\right) \cdot \left(\frac{\mathsf{fma}\left(x2, 2, \left(3 \cdot x1\right) \cdot x1\right) - x1}{\mathsf{fma}\left(x1, x1, 1\right)} \cdot \left(x1 \cdot 2\right)\right)\right), \mathsf{fma}\left(x1, x1, 1\right), \frac{\mathsf{fma}\left(x2, 2, \left(3 \cdot x1\right) \cdot x1\right) - x1}{\mathsf{fma}\left(x1, x1, 1\right)} \cdot \left(\left(3 \cdot x1\right) \cdot x1\right)\right) + \mathsf{fma}\left(\frac{\mathsf{fma}\left(-2, x2, \left(3 \cdot x1\right) \cdot x1\right) - x1}{\mathsf{fma}\left(x1, x1, 1\right)}, 3, x1\right)\right)} \]

   5. Taylor expanded in x1 around 0

      \[\leadsto x1 + \mathsf{fma}\left(x1 \cdot x1, x1, \color{blue}{-6 \cdot x2}\right) \]
   6. Step-by-step derivation

      1. lower-*.f64 92.3

         \[\leadsto x1 + \mathsf{fma}\left(x1 \cdot x1, x1, \color{blue}{-6 \cdot x2}\right) \]

   7. Applied rewrites 92.3%

      \[\leadsto x1 + \mathsf{fma}\left(x1 \cdot x1, x1, \color{blue}{-6 \cdot x2}\right) \]
   8. Step-by-step derivation

      1. lift-+.f64 N/A

         \[\leadsto \color{blue}{x1 + \mathsf{fma}\left(x1 \cdot x1, x1, -6 \cdot x2\right)} \]

   9. Applied rewrites 92.3%

      \[\leadsto \color{blue}{\mathsf{fma}\left(\mathsf{fma}\left(x1, x1, 1\right), x1, -6 \cdot x2\right)} \]
3. Recombined 3 regimes into one program.

4. Final simplification 66.2%

   \[\leadsto \begin{array}{l} \mathbf{if}\;x1 \leq -7.5 \cdot 10^{+137}:\\ \;\;\;\;\left(\left(x2 \cdot x2\right) \cdot -36\right) \cdot {\left(x1 - -6 \cdot x2\right)}^{-1}\\ \mathbf{elif}\;x1 \leq 3.5 \cdot 10^{+88}:\\ \;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(2, x2, -3\right) \cdot x2, 4, -1\right), x1, -6 \cdot x2\right)\\ \mathbf{else}:\\ \;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(x1, x1, 1\right), x1, -6 \cdot x2\right)\\ \end{array} \]
5. Add Preprocessing

Alternative 5: 92.5% accurate, 2.3× speedup

Math

\[\begin{array}{l} \\ \begin{array}{l} \mathbf{if}\;x1 \leq -1.2 \cdot 10^{+38} \lor \neg \left(x1 \leq 1.65 \cdot 10^{+40}\right):\\ \;\;\;\;x1 + \mathsf{fma}\left(6, x1, -3\right) \cdot {x1}^{3}\\ \mathbf{else}:\\ \;\;\;\;x1 + \left(\left(\left(\left(\frac{x1}{\mathsf{fma}\left(x1, x1, 1\right)} \cdot 8\right) \cdot x2\right) \cdot x2 + x1\right) + 3 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 - 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right)\\ \end{array} \end{array} \]

FPCore

(FPCore (x1 x2)
 :precision binary64
 (if (or (<= x1 -1.2e+38) (not (<= x1 1.65e+40)))
   (+ x1 (* (fma 6.0 x1 -3.0) (pow x1 3.0)))
   (+
    x1
    (+
     (+ (* (* (* (/ x1 (fma x1 x1 1.0)) 8.0) x2) x2) x1)
     (* 3.0 (/ (- (- (* (* 3.0 x1) x1) (* 2.0 x2)) x1) (+ (* x1 x1) 1.0)))))))

C

double code(double x1, double x2) {
	double tmp;
	if ((x1 <= -1.2e+38) || !(x1 <= 1.65e+40)) {
		tmp = x1 + (fma(6.0, x1, -3.0) * pow(x1, 3.0));
	} else {
		tmp = x1 + ((((((x1 / fma(x1, x1, 1.0)) * 8.0) * x2) * x2) + x1) + (3.0 * (((((3.0 * x1) * x1) - (2.0 * x2)) - x1) / ((x1 * x1) + 1.0))));
	}
	return tmp;
}

Julia

function code(x1, x2)
	tmp = 0.0
	if ((x1 <= -1.2e+38) || !(x1 <= 1.65e+40))
		tmp = Float64(x1 + Float64(fma(6.0, x1, -3.0) * (x1 ^ 3.0)));
	else
		tmp = Float64(x1 + Float64(Float64(Float64(Float64(Float64(Float64(x1 / fma(x1, x1, 1.0)) * 8.0) * x2) * x2) + x1) + Float64(3.0 * Float64(Float64(Float64(Float64(Float64(3.0 * x1) * x1) - Float64(2.0 * x2)) - x1) / Float64(Float64(x1 * x1) + 1.0)))));
	end
	return tmp
end

Wolfram

code[x1_, x2_] := If[Or[LessEqual[x1, -1.2e+38], N[Not[LessEqual[x1, 1.65e+40]], $MachinePrecision]], N[(x1 + N[(N[(6.0 * x1 + -3.0), $MachinePrecision] * N[Power[x1, 3.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x1 + N[(N[(N[(N[(N[(N[(x1 / N[(x1 * x1 + 1.0), $MachinePrecision]), $MachinePrecision] * 8.0), $MachinePrecision] * x2), $MachinePrecision] * x2), $MachinePrecision] + x1), $MachinePrecision] + N[(3.0 * N[(N[(N[(N[(N[(3.0 * x1), $MachinePrecision] * x1), $MachinePrecision] - N[(2.0 * x2), $MachinePrecision]), $MachinePrecision] - x1), $MachinePrecision] / N[(N[(x1 * x1), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]

Derivation

1. Split input into 2 regimes

2. if x1 < -1.20000000000000009e38 or 1.6499999999999999e40 < x1

   1. Initial program 28.7%

      \[x1 + \left(\left(\left(\left(\left(\left(\left(2 \cdot x1\right) \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 + 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \cdot \left(\frac{\left(\left(3 \cdot x1\right) \cdot x1 + 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1} - 3\right) + \left(x1 \cdot x1\right) \cdot \left(4 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 + 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1} - 6\right)\right) \cdot \left(x1 \cdot x1 + 1\right) + \left(\left(3 \cdot x1\right) \cdot x1\right) \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 + 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) + \left(x1 \cdot x1\right) \cdot x1\right) + x1\right) + 3 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 - 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \]
   2. Add Preprocessing

   4. Applied rewrites 37.2%

      \[\leadsto x1 + \color{blue}{\mathsf{fma}\left(x1 \cdot x1, x1, \mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(4, \frac{\mathsf{fma}\left(x2, 2, \left(3 \cdot x1\right) \cdot x1\right) - x1}{\mathsf{fma}\left(x1, x1, 1\right)}, -6\right), x1 \cdot x1, \left(\frac{\mathsf{fma}\left(x2, 2, \left(3 \cdot x1\right) \cdot x1\right) - x1}{\mathsf{fma}\left(x1, x1, 1\right)} - 3\right) \cdot \left(\frac{\mathsf{fma}\left(x2, 2, \left(3 \cdot x1\right) \cdot x1\right) - x1}{\mathsf{fma}\left(x1, x1, 1\right)} \cdot \left(x1 \cdot 2\right)\right)\right), \mathsf{fma}\left(x1, x1, 1\right), \frac{\mathsf{fma}\left(x2, 2, \left(3 \cdot x1\right) \cdot x1\right) - x1}{\mathsf{fma}\left(x1, x1, 1\right)} \cdot \left(\left(3 \cdot x1\right) \cdot x1\right)\right) + \mathsf{fma}\left(\frac{\mathsf{fma}\left(-2, x2, \left(3 \cdot x1\right) \cdot x1\right) - x1}{\mathsf{fma}\left(x1, x1, 1\right)}, 3, x1\right)\right)} \]

   5. Taylor expanded in x1 around inf

      \[\leadsto x1 + \color{blue}{{x1}^{4} \cdot \left(6 - 3 \cdot \frac{1}{x1}\right)} \]
   6. Step-by-step derivation

      1. *-commutative N/A

         \[\leadsto x1 + \color{blue}{\left(6 - 3 \cdot \frac{1}{x1}\right) \cdot {x1}^{4}} \]

      2. lower-*.f64 N/A

         \[\leadsto x1 + \color{blue}{\left(6 - 3 \cdot \frac{1}{x1}\right) \cdot {x1}^{4}} \]

      3. lower--.f64 N/A

         \[\leadsto x1 + \color{blue}{\left(6 - 3 \cdot \frac{1}{x1}\right)} \cdot {x1}^{4} \]

      4. associate-*r/ N/A

         \[\leadsto x1 + \left(6 - \color{blue}{\frac{3 \cdot 1}{x1}}\right) \cdot {x1}^{4} \]

      5. metadata-eval N/A

         \[\leadsto x1 + \left(6 - \frac{\color{blue}{3}}{x1}\right) \cdot {x1}^{4} \]

      6. lower-/.f64 N/A

         \[\leadsto x1 + \left(6 - \color{blue}{\frac{3}{x1}}\right) \cdot {x1}^{4} \]

      7. lower-pow.f64 96.3

         \[\leadsto x1 + \left(6 - \frac{3}{x1}\right) \cdot \color{blue}{{x1}^{4}} \]

   7. Applied rewrites 96.3%

      \[\leadsto x1 + \color{blue}{\left(6 - \frac{3}{x1}\right) \cdot {x1}^{4}} \]

   8. Taylor expanded in x1 around 0

      \[\leadsto x1 + {x1}^{3} \cdot \color{blue}{\left(6 \cdot x1 - 3\right)} \]
   9. Step-by-step derivation

   10. Applied rewrites 96.3%

       \[\leadsto x1 + \mathsf{fma}\left(6, x1, -3\right) \cdot \color{blue}{{x1}^{3}} \]

   if -1.20000000000000009e38 < x1 < 1.6499999999999999e40

   1. Initial program 99.3%

      \[x1 + \left(\left(\left(\left(\left(\left(\left(2 \cdot x1\right) \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 + 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \cdot \left(\frac{\left(\left(3 \cdot x1\right) \cdot x1 + 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1} - 3\right) + \left(x1 \cdot x1\right) \cdot \left(4 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 + 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1} - 6\right)\right) \cdot \left(x1 \cdot x1 + 1\right) + \left(\left(3 \cdot x1\right) \cdot x1\right) \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 + 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) + \left(x1 \cdot x1\right) \cdot x1\right) + x1\right) + 3 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 - 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \]
   2. Add Preprocessing

   3. Taylor expanded in x2 around inf

      \[\leadsto x1 + \left(\left(\color{blue}{8 \cdot \frac{x1 \cdot {x2}^{2}}{1 + {x1}^{2}}} + x1\right) + 3 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 - 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \]
   4. Step-by-step derivation

      1. associate-*r/ N/A

         \[\leadsto x1 + \left(\left(\color{blue}{\frac{8 \cdot \left(x1 \cdot {x2}^{2}\right)}{1 + {x1}^{2}}} + x1\right) + 3 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 - 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \]

      2. associate-*r* N/A

         \[\leadsto x1 + \left(\left(\frac{\color{blue}{\left(8 \cdot x1\right) \cdot {x2}^{2}}}{1 + {x1}^{2}} + x1\right) + 3 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 - 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \]

      3. associate-*l/ N/A

         \[\leadsto x1 + \left(\left(\color{blue}{\frac{8 \cdot x1}{1 + {x1}^{2}} \cdot {x2}^{2}} + x1\right) + 3 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 - 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \]

      4. associate-*r/ N/A

         \[\leadsto x1 + \left(\left(\color{blue}{\left(8 \cdot \frac{x1}{1 + {x1}^{2}}\right)} \cdot {x2}^{2} + x1\right) + 3 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 - 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \]

      5. unpow2 N/A

         \[\leadsto x1 + \left(\left(\left(8 \cdot \frac{x1}{1 + {x1}^{2}}\right) \cdot \color{blue}{\left(x2 \cdot x2\right)} + x1\right) + 3 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 - 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \]

      6. associate-*r* N/A

         \[\leadsto x1 + \left(\left(\color{blue}{\left(\left(8 \cdot \frac{x1}{1 + {x1}^{2}}\right) \cdot x2\right) \cdot x2} + x1\right) + 3 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 - 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \]

      7. lower-*.f64 N/A

         \[\leadsto x1 + \left(\left(\color{blue}{\left(\left(8 \cdot \frac{x1}{1 + {x1}^{2}}\right) \cdot x2\right) \cdot x2} + x1\right) + 3 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 - 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \]

      8. lower-*.f64 N/A

         \[\leadsto x1 + \left(\left(\color{blue}{\left(\left(8 \cdot \frac{x1}{1 + {x1}^{2}}\right) \cdot x2\right)} \cdot x2 + x1\right) + 3 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 - 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \]

      9. *-commutative N/A

         \[\leadsto x1 + \left(\left(\left(\color{blue}{\left(\frac{x1}{1 + {x1}^{2}} \cdot 8\right)} \cdot x2\right) \cdot x2 + x1\right) + 3 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 - 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \]

      10. lower-*.f64 N/A

          \[\leadsto x1 + \left(\left(\left(\color{blue}{\left(\frac{x1}{1 + {x1}^{2}} \cdot 8\right)} \cdot x2\right) \cdot x2 + x1\right) + 3 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 - 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \]

      11. lower-/.f64 N/A

          \[\leadsto x1 + \left(\left(\left(\left(\color{blue}{\frac{x1}{1 + {x1}^{2}}} \cdot 8\right) \cdot x2\right) \cdot x2 + x1\right) + 3 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 - 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \]

      12. +-commutative N/A

          \[\leadsto x1 + \left(\left(\left(\left(\frac{x1}{\color{blue}{{x1}^{2} + 1}} \cdot 8\right) \cdot x2\right) \cdot x2 + x1\right) + 3 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 - 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \]

      13. unpow2 N/A

          \[\leadsto x1 + \left(\left(\left(\left(\frac{x1}{\color{blue}{x1 \cdot x1} + 1} \cdot 8\right) \cdot x2\right) \cdot x2 + x1\right) + 3 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 - 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \]

      14. lower-fma.f64 96.9

          \[\leadsto x1 + \left(\left(\left(\left(\frac{x1}{\color{blue}{\mathsf{fma}\left(x1, x1, 1\right)}} \cdot 8\right) \cdot x2\right) \cdot x2 + x1\right) + 3 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 - 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \]

   5. Applied rewrites 96.9%

      \[\leadsto x1 + \left(\left(\color{blue}{\left(\left(\frac{x1}{\mathsf{fma}\left(x1, x1, 1\right)} \cdot 8\right) \cdot x2\right) \cdot x2} + x1\right) + 3 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 - 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \]
3. Recombined 2 regimes into one program.

4. Final simplification 96.7%

   \[\leadsto \begin{array}{l} \mathbf{if}\;x1 \leq -1.2 \cdot 10^{+38} \lor \neg \left(x1 \leq 1.65 \cdot 10^{+40}\right):\\ \;\;\;\;x1 + \mathsf{fma}\left(6, x1, -3\right) \cdot {x1}^{3}\\ \mathbf{else}:\\ \;\;\;\;x1 + \left(\left(\left(\left(\frac{x1}{\mathsf{fma}\left(x1, x1, 1\right)} \cdot 8\right) \cdot x2\right) \cdot x2 + x1\right) + 3 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 - 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right)\\ \end{array} \]
5. Add Preprocessing

Alternative 6: 92.5% accurate, 2.4× speedup

Math

\[\begin{array}{l} \\ \begin{array}{l} \mathbf{if}\;x1 \leq -1.2 \cdot 10^{+38}:\\ \;\;\;\;x1 + \left(\left(6 - \frac{3}{x1}\right) \cdot \left(x1 \cdot x1\right)\right) \cdot \left(x1 \cdot x1\right)\\ \mathbf{elif}\;x1 \leq 1.65 \cdot 10^{+40}:\\ \;\;\;\;x1 + \left(\left(\left(\left(\frac{x1}{\mathsf{fma}\left(x1, x1, 1\right)} \cdot 8\right) \cdot x2\right) \cdot x2 + x1\right) + 3 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 - 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right)\\ \mathbf{else}:\\ \;\;\;\;x1 + 6 \cdot {x1}^{4}\\ \end{array} \end{array} \]

FPCore

(FPCore (x1 x2)
 :precision binary64
 (if (<= x1 -1.2e+38)
   (+ x1 (* (* (- 6.0 (/ 3.0 x1)) (* x1 x1)) (* x1 x1)))
   (if (<= x1 1.65e+40)
     (+
      x1
      (+
       (+ (* (* (* (/ x1 (fma x1 x1 1.0)) 8.0) x2) x2) x1)
       (* 3.0 (/ (- (- (* (* 3.0 x1) x1) (* 2.0 x2)) x1) (+ (* x1 x1) 1.0)))))
     (+ x1 (* 6.0 (pow x1 4.0))))))

C

double code(double x1, double x2) {
	double tmp;
	if (x1 <= -1.2e+38) {
		tmp = x1 + (((6.0 - (3.0 / x1)) * (x1 * x1)) * (x1 * x1));
	} else if (x1 <= 1.65e+40) {
		tmp = x1 + ((((((x1 / fma(x1, x1, 1.0)) * 8.0) * x2) * x2) + x1) + (3.0 * (((((3.0 * x1) * x1) - (2.0 * x2)) - x1) / ((x1 * x1) + 1.0))));
	} else {
		tmp = x1 + (6.0 * pow(x1, 4.0));
	}
	return tmp;
}

Julia

function code(x1, x2)
	tmp = 0.0
	if (x1 <= -1.2e+38)
		tmp = Float64(x1 + Float64(Float64(Float64(6.0 - Float64(3.0 / x1)) * Float64(x1 * x1)) * Float64(x1 * x1)));
	elseif (x1 <= 1.65e+40)
		tmp = Float64(x1 + Float64(Float64(Float64(Float64(Float64(Float64(x1 / fma(x1, x1, 1.0)) * 8.0) * x2) * x2) + x1) + Float64(3.0 * Float64(Float64(Float64(Float64(Float64(3.0 * x1) * x1) - Float64(2.0 * x2)) - x1) / Float64(Float64(x1 * x1) + 1.0)))));
	else
		tmp = Float64(x1 + Float64(6.0 * (x1 ^ 4.0)));
	end
	return tmp
end

Wolfram

code[x1_, x2_] := If[LessEqual[x1, -1.2e+38], N[(x1 + N[(N[(N[(6.0 - N[(3.0 / x1), $MachinePrecision]), $MachinePrecision] * N[(x1 * x1), $MachinePrecision]), $MachinePrecision] * N[(x1 * x1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[x1, 1.65e+40], N[(x1 + N[(N[(N[(N[(N[(N[(x1 / N[(x1 * x1 + 1.0), $MachinePrecision]), $MachinePrecision] * 8.0), $MachinePrecision] * x2), $MachinePrecision] * x2), $MachinePrecision] + x1), $MachinePrecision] + N[(3.0 * N[(N[(N[(N[(N[(3.0 * x1), $MachinePrecision] * x1), $MachinePrecision] - N[(2.0 * x2), $MachinePrecision]), $MachinePrecision] - x1), $MachinePrecision] / N[(N[(x1 * x1), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x1 + N[(6.0 * N[Power[x1, 4.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]

Derivation

1. Split input into 3 regimes

2. if x1 < -1.20000000000000009e38

   1. Initial program 19.0%

      \[x1 + \left(\left(\left(\left(\left(\left(\left(2 \cdot x1\right) \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 + 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \cdot \left(\frac{\left(\left(3 \cdot x1\right) \cdot x1 + 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1} - 3\right) + \left(x1 \cdot x1\right) \cdot \left(4 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 + 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1} - 6\right)\right) \cdot \left(x1 \cdot x1 + 1\right) + \left(\left(3 \cdot x1\right) \cdot x1\right) \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 + 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) + \left(x1 \cdot x1\right) \cdot x1\right) + x1\right) + 3 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 - 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \]
   2. Add Preprocessing

   4. Applied rewrites 34.9%

      \[\leadsto x1 + \color{blue}{\mathsf{fma}\left(x1 \cdot x1, x1, \mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(4, \frac{\mathsf{fma}\left(x2, 2, \left(3 \cdot x1\right) \cdot x1\right) - x1}{\mathsf{fma}\left(x1, x1, 1\right)}, -6\right), x1 \cdot x1, \left(\frac{\mathsf{fma}\left(x2, 2, \left(3 \cdot x1\right) \cdot x1\right) - x1}{\mathsf{fma}\left(x1, x1, 1\right)} - 3\right) \cdot \left(\frac{\mathsf{fma}\left(x2, 2, \left(3 \cdot x1\right) \cdot x1\right) - x1}{\mathsf{fma}\left(x1, x1, 1\right)} \cdot \left(x1 \cdot 2\right)\right)\right), \mathsf{fma}\left(x1, x1, 1\right), \frac{\mathsf{fma}\left(x2, 2, \left(3 \cdot x1\right) \cdot x1\right) - x1}{\mathsf{fma}\left(x1, x1, 1\right)} \cdot \left(\left(3 \cdot x1\right) \cdot x1\right)\right) + \mathsf{fma}\left(\frac{\mathsf{fma}\left(-2, x2, \left(3 \cdot x1\right) \cdot x1\right) - x1}{\mathsf{fma}\left(x1, x1, 1\right)}, 3, x1\right)\right)} \]

   5. Taylor expanded in x1 around inf

      \[\leadsto x1 + \color{blue}{{x1}^{4} \cdot \left(6 - 3 \cdot \frac{1}{x1}\right)} \]
   6. Step-by-step derivation

      1. *-commutative N/A

         \[\leadsto x1 + \color{blue}{\left(6 - 3 \cdot \frac{1}{x1}\right) \cdot {x1}^{4}} \]

      2. lower-*.f64 N/A

         \[\leadsto x1 + \color{blue}{\left(6 - 3 \cdot \frac{1}{x1}\right) \cdot {x1}^{4}} \]

      3. lower--.f64 N/A

         \[\leadsto x1 + \color{blue}{\left(6 - 3 \cdot \frac{1}{x1}\right)} \cdot {x1}^{4} \]

      4. associate-*r/ N/A

         \[\leadsto x1 + \left(6 - \color{blue}{\frac{3 \cdot 1}{x1}}\right) \cdot {x1}^{4} \]

      5. metadata-eval N/A

         \[\leadsto x1 + \left(6 - \frac{\color{blue}{3}}{x1}\right) \cdot {x1}^{4} \]

      6. lower-/.f64 N/A

         \[\leadsto x1 + \left(6 - \color{blue}{\frac{3}{x1}}\right) \cdot {x1}^{4} \]

      7. lower-pow.f64 96.8

         \[\leadsto x1 + \left(6 - \frac{3}{x1}\right) \cdot \color{blue}{{x1}^{4}} \]

   7. Applied rewrites 96.8%

      \[\leadsto x1 + \color{blue}{\left(6 - \frac{3}{x1}\right) \cdot {x1}^{4}} \]
   8. Step-by-step derivation

   9. Applied rewrites 96.8%

      \[\leadsto x1 + \left(\left(6 - \frac{3}{x1}\right) \cdot \left(x1 \cdot x1\right)\right) \cdot \color{blue}{\left(x1 \cdot x1\right)} \]

   if -1.20000000000000009e38 < x1 < 1.6499999999999999e40

   1. Initial program 99.3%

      \[x1 + \left(\left(\left(\left(\left(\left(\left(2 \cdot x1\right) \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 + 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \cdot \left(\frac{\left(\left(3 \cdot x1\right) \cdot x1 + 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1} - 3\right) + \left(x1 \cdot x1\right) \cdot \left(4 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 + 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1} - 6\right)\right) \cdot \left(x1 \cdot x1 + 1\right) + \left(\left(3 \cdot x1\right) \cdot x1\right) \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 + 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) + \left(x1 \cdot x1\right) \cdot x1\right) + x1\right) + 3 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 - 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \]
   2. Add Preprocessing

   3. Taylor expanded in x2 around inf

      \[\leadsto x1 + \left(\left(\color{blue}{8 \cdot \frac{x1 \cdot {x2}^{2}}{1 + {x1}^{2}}} + x1\right) + 3 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 - 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \]
   4. Step-by-step derivation

      1. associate-*r/ N/A

         \[\leadsto x1 + \left(\left(\color{blue}{\frac{8 \cdot \left(x1 \cdot {x2}^{2}\right)}{1 + {x1}^{2}}} + x1\right) + 3 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 - 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \]

      2. associate-*r* N/A

         \[\leadsto x1 + \left(\left(\frac{\color{blue}{\left(8 \cdot x1\right) \cdot {x2}^{2}}}{1 + {x1}^{2}} + x1\right) + 3 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 - 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \]

      3. associate-*l/ N/A

         \[\leadsto x1 + \left(\left(\color{blue}{\frac{8 \cdot x1}{1 + {x1}^{2}} \cdot {x2}^{2}} + x1\right) + 3 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 - 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \]

      4. associate-*r/ N/A

         \[\leadsto x1 + \left(\left(\color{blue}{\left(8 \cdot \frac{x1}{1 + {x1}^{2}}\right)} \cdot {x2}^{2} + x1\right) + 3 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 - 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \]

      5. unpow2 N/A

         \[\leadsto x1 + \left(\left(\left(8 \cdot \frac{x1}{1 + {x1}^{2}}\right) \cdot \color{blue}{\left(x2 \cdot x2\right)} + x1\right) + 3 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 - 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \]

      6. associate-*r* N/A

         \[\leadsto x1 + \left(\left(\color{blue}{\left(\left(8 \cdot \frac{x1}{1 + {x1}^{2}}\right) \cdot x2\right) \cdot x2} + x1\right) + 3 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 - 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \]

      7. lower-*.f64 N/A

         \[\leadsto x1 + \left(\left(\color{blue}{\left(\left(8 \cdot \frac{x1}{1 + {x1}^{2}}\right) \cdot x2\right) \cdot x2} + x1\right) + 3 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 - 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \]

      8. lower-*.f64 N/A

         \[\leadsto x1 + \left(\left(\color{blue}{\left(\left(8 \cdot \frac{x1}{1 + {x1}^{2}}\right) \cdot x2\right)} \cdot x2 + x1\right) + 3 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 - 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \]

      9. *-commutative N/A

         \[\leadsto x1 + \left(\left(\left(\color{blue}{\left(\frac{x1}{1 + {x1}^{2}} \cdot 8\right)} \cdot x2\right) \cdot x2 + x1\right) + 3 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 - 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \]

      10. lower-*.f64 N/A

          \[\leadsto x1 + \left(\left(\left(\color{blue}{\left(\frac{x1}{1 + {x1}^{2}} \cdot 8\right)} \cdot x2\right) \cdot x2 + x1\right) + 3 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 - 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \]

      11. lower-/.f64 N/A

          \[\leadsto x1 + \left(\left(\left(\left(\color{blue}{\frac{x1}{1 + {x1}^{2}}} \cdot 8\right) \cdot x2\right) \cdot x2 + x1\right) + 3 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 - 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \]

      12. +-commutative N/A

          \[\leadsto x1 + \left(\left(\left(\left(\frac{x1}{\color{blue}{{x1}^{2} + 1}} \cdot 8\right) \cdot x2\right) \cdot x2 + x1\right) + 3 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 - 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \]

      13. unpow2 N/A

          \[\leadsto x1 + \left(\left(\left(\left(\frac{x1}{\color{blue}{x1 \cdot x1} + 1} \cdot 8\right) \cdot x2\right) \cdot x2 + x1\right) + 3 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 - 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \]

      14. lower-fma.f64 96.9

          \[\leadsto x1 + \left(\left(\left(\left(\frac{x1}{\color{blue}{\mathsf{fma}\left(x1, x1, 1\right)}} \cdot 8\right) \cdot x2\right) \cdot x2 + x1\right) + 3 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 - 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \]

   5. Applied rewrites 96.9%

      \[\leadsto x1 + \left(\left(\color{blue}{\left(\left(\frac{x1}{\mathsf{fma}\left(x1, x1, 1\right)} \cdot 8\right) \cdot x2\right) \cdot x2} + x1\right) + 3 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 - 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \]

   if 1.6499999999999999e40 < x1

   1. Initial program 39.9%

      \[x1 + \left(\left(\left(\left(\left(\left(\left(2 \cdot x1\right) \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 + 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \cdot \left(\frac{\left(\left(3 \cdot x1\right) \cdot x1 + 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1} - 3\right) + \left(x1 \cdot x1\right) \cdot \left(4 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 + 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1} - 6\right)\right) \cdot \left(x1 \cdot x1 + 1\right) + \left(\left(3 \cdot x1\right) \cdot x1\right) \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 + 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) + \left(x1 \cdot x1\right) \cdot x1\right) + x1\right) + 3 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 - 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \]
   2. Add Preprocessing

   4. Applied rewrites 39.9%

      \[\leadsto x1 + \color{blue}{\mathsf{fma}\left(x1 \cdot x1, x1, \mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(4, \frac{\mathsf{fma}\left(x2, 2, \left(3 \cdot x1\right) \cdot x1\right) - x1}{\mathsf{fma}\left(x1, x1, 1\right)}, -6\right), x1 \cdot x1, \left(\frac{\mathsf{fma}\left(x2, 2, \left(3 \cdot x1\right) \cdot x1\right) - x1}{\mathsf{fma}\left(x1, x1, 1\right)} - 3\right) \cdot \left(\frac{\mathsf{fma}\left(x2, 2, \left(3 \cdot x1\right) \cdot x1\right) - x1}{\mathsf{fma}\left(x1, x1, 1\right)} \cdot \left(x1 \cdot 2\right)\right)\right), \mathsf{fma}\left(x1, x1, 1\right), \frac{\mathsf{fma}\left(x2, 2, \left(3 \cdot x1\right) \cdot x1\right) - x1}{\mathsf{fma}\left(x1, x1, 1\right)} \cdot \left(\left(3 \cdot x1\right) \cdot x1\right)\right) + \mathsf{fma}\left(\frac{\mathsf{fma}\left(-2, x2, \left(3 \cdot x1\right) \cdot x1\right) - x1}{\mathsf{fma}\left(x1, x1, 1\right)}, 3, x1\right)\right)} \]

   5. Taylor expanded in x1 around inf

      \[\leadsto x1 + \color{blue}{{x1}^{4} \cdot \left(6 - 3 \cdot \frac{1}{x1}\right)} \]
   6. Step-by-step derivation

      1. *-commutative N/A

         \[\leadsto x1 + \color{blue}{\left(6 - 3 \cdot \frac{1}{x1}\right) \cdot {x1}^{4}} \]

      2. lower-*.f64 N/A

         \[\leadsto x1 + \color{blue}{\left(6 - 3 \cdot \frac{1}{x1}\right) \cdot {x1}^{4}} \]

      3. lower--.f64 N/A

         \[\leadsto x1 + \color{blue}{\left(6 - 3 \cdot \frac{1}{x1}\right)} \cdot {x1}^{4} \]

      4. associate-*r/ N/A

         \[\leadsto x1 + \left(6 - \color{blue}{\frac{3 \cdot 1}{x1}}\right) \cdot {x1}^{4} \]

      5. metadata-eval N/A

         \[\leadsto x1 + \left(6 - \frac{\color{blue}{3}}{x1}\right) \cdot {x1}^{4} \]

      6. lower-/.f64 N/A

         \[\leadsto x1 + \left(6 - \color{blue}{\frac{3}{x1}}\right) \cdot {x1}^{4} \]

      7. lower-pow.f64 95.8

         \[\leadsto x1 + \left(6 - \frac{3}{x1}\right) \cdot \color{blue}{{x1}^{4}} \]

   7. Applied rewrites 95.8%

      \[\leadsto x1 + \color{blue}{\left(6 - \frac{3}{x1}\right) \cdot {x1}^{4}} \]

   8. Taylor expanded in x1 around inf

      \[\leadsto x1 + 6 \cdot {\color{blue}{x1}}^{4} \]
   9. Step-by-step derivation

   10. Applied rewrites 95.8%

       \[\leadsto x1 + 6 \cdot {\color{blue}{x1}}^{4} \]
3. Recombined 3 regimes into one program.
4. Add Preprocessing

Alternative 7: 92.5% accurate, 3.0× speedup

Math

\[\begin{array}{l} \\ \begin{array}{l} \mathbf{if}\;x1 \leq -1.2 \cdot 10^{+38} \lor \neg \left(x1 \leq 1.65 \cdot 10^{+40}\right):\\ \;\;\;\;x1 + \left(\left(6 - \frac{3}{x1}\right) \cdot \left(x1 \cdot x1\right)\right) \cdot \left(x1 \cdot x1\right)\\ \mathbf{else}:\\ \;\;\;\;x1 + \left(\left(\left(\left(\frac{x1}{\mathsf{fma}\left(x1, x1, 1\right)} \cdot 8\right) \cdot x2\right) \cdot x2 + x1\right) + 3 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 - 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right)\\ \end{array} \end{array} \]

FPCore

(FPCore (x1 x2)
 :precision binary64
 (if (or (<= x1 -1.2e+38) (not (<= x1 1.65e+40)))
   (+ x1 (* (* (- 6.0 (/ 3.0 x1)) (* x1 x1)) (* x1 x1)))
   (+
    x1
    (+
     (+ (* (* (* (/ x1 (fma x1 x1 1.0)) 8.0) x2) x2) x1)
     (* 3.0 (/ (- (- (* (* 3.0 x1) x1) (* 2.0 x2)) x1) (+ (* x1 x1) 1.0)))))))

C

double code(double x1, double x2) {
	double tmp;
	if ((x1 <= -1.2e+38) || !(x1 <= 1.65e+40)) {
		tmp = x1 + (((6.0 - (3.0 / x1)) * (x1 * x1)) * (x1 * x1));
	} else {
		tmp = x1 + ((((((x1 / fma(x1, x1, 1.0)) * 8.0) * x2) * x2) + x1) + (3.0 * (((((3.0 * x1) * x1) - (2.0 * x2)) - x1) / ((x1 * x1) + 1.0))));
	}
	return tmp;
}

Julia

function code(x1, x2)
	tmp = 0.0
	if ((x1 <= -1.2e+38) || !(x1 <= 1.65e+40))
		tmp = Float64(x1 + Float64(Float64(Float64(6.0 - Float64(3.0 / x1)) * Float64(x1 * x1)) * Float64(x1 * x1)));
	else
		tmp = Float64(x1 + Float64(Float64(Float64(Float64(Float64(Float64(x1 / fma(x1, x1, 1.0)) * 8.0) * x2) * x2) + x1) + Float64(3.0 * Float64(Float64(Float64(Float64(Float64(3.0 * x1) * x1) - Float64(2.0 * x2)) - x1) / Float64(Float64(x1 * x1) + 1.0)))));
	end
	return tmp
end

Wolfram

code[x1_, x2_] := If[Or[LessEqual[x1, -1.2e+38], N[Not[LessEqual[x1, 1.65e+40]], $MachinePrecision]], N[(x1 + N[(N[(N[(6.0 - N[(3.0 / x1), $MachinePrecision]), $MachinePrecision] * N[(x1 * x1), $MachinePrecision]), $MachinePrecision] * N[(x1 * x1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x1 + N[(N[(N[(N[(N[(N[(x1 / N[(x1 * x1 + 1.0), $MachinePrecision]), $MachinePrecision] * 8.0), $MachinePrecision] * x2), $MachinePrecision] * x2), $MachinePrecision] + x1), $MachinePrecision] + N[(3.0 * N[(N[(N[(N[(N[(3.0 * x1), $MachinePrecision] * x1), $MachinePrecision] - N[(2.0 * x2), $MachinePrecision]), $MachinePrecision] - x1), $MachinePrecision] / N[(N[(x1 * x1), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]

Derivation

1. Split input into 2 regimes

2. if x1 < -1.20000000000000009e38 or 1.6499999999999999e40 < x1

   1. Initial program 28.7%

      \[x1 + \left(\left(\left(\left(\left(\left(\left(2 \cdot x1\right) \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 + 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \cdot \left(\frac{\left(\left(3 \cdot x1\right) \cdot x1 + 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1} - 3\right) + \left(x1 \cdot x1\right) \cdot \left(4 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 + 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1} - 6\right)\right) \cdot \left(x1 \cdot x1 + 1\right) + \left(\left(3 \cdot x1\right) \cdot x1\right) \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 + 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) + \left(x1 \cdot x1\right) \cdot x1\right) + x1\right) + 3 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 - 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \]
   2. Add Preprocessing

   4. Applied rewrites 37.2%

      \[\leadsto x1 + \color{blue}{\mathsf{fma}\left(x1 \cdot x1, x1, \mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(4, \frac{\mathsf{fma}\left(x2, 2, \left(3 \cdot x1\right) \cdot x1\right) - x1}{\mathsf{fma}\left(x1, x1, 1\right)}, -6\right), x1 \cdot x1, \left(\frac{\mathsf{fma}\left(x2, 2, \left(3 \cdot x1\right) \cdot x1\right) - x1}{\mathsf{fma}\left(x1, x1, 1\right)} - 3\right) \cdot \left(\frac{\mathsf{fma}\left(x2, 2, \left(3 \cdot x1\right) \cdot x1\right) - x1}{\mathsf{fma}\left(x1, x1, 1\right)} \cdot \left(x1 \cdot 2\right)\right)\right), \mathsf{fma}\left(x1, x1, 1\right), \frac{\mathsf{fma}\left(x2, 2, \left(3 \cdot x1\right) \cdot x1\right) - x1}{\mathsf{fma}\left(x1, x1, 1\right)} \cdot \left(\left(3 \cdot x1\right) \cdot x1\right)\right) + \mathsf{fma}\left(\frac{\mathsf{fma}\left(-2, x2, \left(3 \cdot x1\right) \cdot x1\right) - x1}{\mathsf{fma}\left(x1, x1, 1\right)}, 3, x1\right)\right)} \]

   5. Taylor expanded in x1 around inf

      \[\leadsto x1 + \color{blue}{{x1}^{4} \cdot \left(6 - 3 \cdot \frac{1}{x1}\right)} \]
   6. Step-by-step derivation

      1. *-commutative N/A

         \[\leadsto x1 + \color{blue}{\left(6 - 3 \cdot \frac{1}{x1}\right) \cdot {x1}^{4}} \]

      2. lower-*.f64 N/A

         \[\leadsto x1 + \color{blue}{\left(6 - 3 \cdot \frac{1}{x1}\right) \cdot {x1}^{4}} \]

      3. lower--.f64 N/A

         \[\leadsto x1 + \color{blue}{\left(6 - 3 \cdot \frac{1}{x1}\right)} \cdot {x1}^{4} \]

      4. associate-*r/ N/A

         \[\leadsto x1 + \left(6 - \color{blue}{\frac{3 \cdot 1}{x1}}\right) \cdot {x1}^{4} \]

      5. metadata-eval N/A

         \[\leadsto x1 + \left(6 - \frac{\color{blue}{3}}{x1}\right) \cdot {x1}^{4} \]

      6. lower-/.f64 N/A

         \[\leadsto x1 + \left(6 - \color{blue}{\frac{3}{x1}}\right) \cdot {x1}^{4} \]

      7. lower-pow.f64 96.3

         \[\leadsto x1 + \left(6 - \frac{3}{x1}\right) \cdot \color{blue}{{x1}^{4}} \]

   7. Applied rewrites 96.3%

      \[\leadsto x1 + \color{blue}{\left(6 - \frac{3}{x1}\right) \cdot {x1}^{4}} \]
   8. Step-by-step derivation

   9. Applied rewrites 96.3%

      \[\leadsto x1 + \left(\left(6 - \frac{3}{x1}\right) \cdot \left(x1 \cdot x1\right)\right) \cdot \color{blue}{\left(x1 \cdot x1\right)} \]

   if -1.20000000000000009e38 < x1 < 1.6499999999999999e40

   1. Initial program 99.3%

      \[x1 + \left(\left(\left(\left(\left(\left(\left(2 \cdot x1\right) \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 + 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \cdot \left(\frac{\left(\left(3 \cdot x1\right) \cdot x1 + 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1} - 3\right) + \left(x1 \cdot x1\right) \cdot \left(4 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 + 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1} - 6\right)\right) \cdot \left(x1 \cdot x1 + 1\right) + \left(\left(3 \cdot x1\right) \cdot x1\right) \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 + 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) + \left(x1 \cdot x1\right) \cdot x1\right) + x1\right) + 3 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 - 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \]
   2. Add Preprocessing

   3. Taylor expanded in x2 around inf

      \[\leadsto x1 + \left(\left(\color{blue}{8 \cdot \frac{x1 \cdot {x2}^{2}}{1 + {x1}^{2}}} + x1\right) + 3 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 - 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \]
   4. Step-by-step derivation

      1. associate-*r/ N/A

         \[\leadsto x1 + \left(\left(\color{blue}{\frac{8 \cdot \left(x1 \cdot {x2}^{2}\right)}{1 + {x1}^{2}}} + x1\right) + 3 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 - 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \]

      2. associate-*r* N/A

         \[\leadsto x1 + \left(\left(\frac{\color{blue}{\left(8 \cdot x1\right) \cdot {x2}^{2}}}{1 + {x1}^{2}} + x1\right) + 3 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 - 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \]

      3. associate-*l/ N/A

         \[\leadsto x1 + \left(\left(\color{blue}{\frac{8 \cdot x1}{1 + {x1}^{2}} \cdot {x2}^{2}} + x1\right) + 3 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 - 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \]

      4. associate-*r/ N/A

         \[\leadsto x1 + \left(\left(\color{blue}{\left(8 \cdot \frac{x1}{1 + {x1}^{2}}\right)} \cdot {x2}^{2} + x1\right) + 3 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 - 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \]

      5. unpow2 N/A

         \[\leadsto x1 + \left(\left(\left(8 \cdot \frac{x1}{1 + {x1}^{2}}\right) \cdot \color{blue}{\left(x2 \cdot x2\right)} + x1\right) + 3 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 - 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \]

      6. associate-*r* N/A

         \[\leadsto x1 + \left(\left(\color{blue}{\left(\left(8 \cdot \frac{x1}{1 + {x1}^{2}}\right) \cdot x2\right) \cdot x2} + x1\right) + 3 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 - 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \]

      7. lower-*.f64 N/A

         \[\leadsto x1 + \left(\left(\color{blue}{\left(\left(8 \cdot \frac{x1}{1 + {x1}^{2}}\right) \cdot x2\right) \cdot x2} + x1\right) + 3 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 - 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \]

      8. lower-*.f64 N/A

         \[\leadsto x1 + \left(\left(\color{blue}{\left(\left(8 \cdot \frac{x1}{1 + {x1}^{2}}\right) \cdot x2\right)} \cdot x2 + x1\right) + 3 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 - 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \]

      9. *-commutative N/A

         \[\leadsto x1 + \left(\left(\left(\color{blue}{\left(\frac{x1}{1 + {x1}^{2}} \cdot 8\right)} \cdot x2\right) \cdot x2 + x1\right) + 3 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 - 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \]

      10. lower-*.f64 N/A

          \[\leadsto x1 + \left(\left(\left(\color{blue}{\left(\frac{x1}{1 + {x1}^{2}} \cdot 8\right)} \cdot x2\right) \cdot x2 + x1\right) + 3 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 - 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \]

      11. lower-/.f64 N/A

          \[\leadsto x1 + \left(\left(\left(\left(\color{blue}{\frac{x1}{1 + {x1}^{2}}} \cdot 8\right) \cdot x2\right) \cdot x2 + x1\right) + 3 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 - 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \]

      12. +-commutative N/A

          \[\leadsto x1 + \left(\left(\left(\left(\frac{x1}{\color{blue}{{x1}^{2} + 1}} \cdot 8\right) \cdot x2\right) \cdot x2 + x1\right) + 3 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 - 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \]

      13. unpow2 N/A

          \[\leadsto x1 + \left(\left(\left(\left(\frac{x1}{\color{blue}{x1 \cdot x1} + 1} \cdot 8\right) \cdot x2\right) \cdot x2 + x1\right) + 3 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 - 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \]

      14. lower-fma.f64 96.9

          \[\leadsto x1 + \left(\left(\left(\left(\frac{x1}{\color{blue}{\mathsf{fma}\left(x1, x1, 1\right)}} \cdot 8\right) \cdot x2\right) \cdot x2 + x1\right) + 3 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 - 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \]

   5. Applied rewrites 96.9%

      \[\leadsto x1 + \left(\left(\color{blue}{\left(\left(\frac{x1}{\mathsf{fma}\left(x1, x1, 1\right)} \cdot 8\right) \cdot x2\right) \cdot x2} + x1\right) + 3 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 - 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \]
3. Recombined 2 regimes into one program.

4. Final simplification 96.6%

   \[\leadsto \begin{array}{l} \mathbf{if}\;x1 \leq -1.2 \cdot 10^{+38} \lor \neg \left(x1 \leq 1.65 \cdot 10^{+40}\right):\\ \;\;\;\;x1 + \left(\left(6 - \frac{3}{x1}\right) \cdot \left(x1 \cdot x1\right)\right) \cdot \left(x1 \cdot x1\right)\\ \mathbf{else}:\\ \;\;\;\;x1 + \left(\left(\left(\left(\frac{x1}{\mathsf{fma}\left(x1, x1, 1\right)} \cdot 8\right) \cdot x2\right) \cdot x2 + x1\right) + 3 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 - 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right)\\ \end{array} \]
5. Add Preprocessing

Alternative 8: 91.9% accurate, 3.4× speedup

Math

\[\begin{array}{l} \\ \begin{array}{l} \mathbf{if}\;x1 \leq -2.6 \cdot 10^{+37} \lor \neg \left(x1 \leq 1.75 \cdot 10^{+35}\right):\\ \;\;\;\;x1 + \left(\left(6 - \frac{3}{x1}\right) \cdot \left(x1 \cdot x1\right)\right) \cdot \left(x1 \cdot x1\right)\\ \mathbf{else}:\\ \;\;\;\;x1 + \mathsf{fma}\left(\frac{\mathsf{fma}\left(-2, x2, \left(x1 \cdot x1\right) \cdot 3\right) - x1}{\mathsf{fma}\left(x1, x1, 1\right)}, 3, \mathsf{fma}\left(\mathsf{fma}\left(6, x1, -12\right), x1, 8 \cdot \left(x2 \cdot x1\right)\right) \cdot x2 + x1\right)\\ \end{array} \end{array} \]

FPCore

(FPCore (x1 x2)
 :precision binary64
 (if (or (<= x1 -2.6e+37) (not (<= x1 1.75e+35)))
   (+ x1 (* (* (- 6.0 (/ 3.0 x1)) (* x1 x1)) (* x1 x1)))
   (+
    x1
    (fma
     (/ (- (fma -2.0 x2 (* (* x1 x1) 3.0)) x1) (fma x1 x1 1.0))
     3.0
     (+ (* (fma (fma 6.0 x1 -12.0) x1 (* 8.0 (* x2 x1))) x2) x1)))))

C

double code(double x1, double x2) {
	double tmp;
	if ((x1 <= -2.6e+37) || !(x1 <= 1.75e+35)) {
		tmp = x1 + (((6.0 - (3.0 / x1)) * (x1 * x1)) * (x1 * x1));
	} else {
		tmp = x1 + fma(((fma(-2.0, x2, ((x1 * x1) * 3.0)) - x1) / fma(x1, x1, 1.0)), 3.0, ((fma(fma(6.0, x1, -12.0), x1, (8.0 * (x2 * x1))) * x2) + x1));
	}
	return tmp;
}

Julia

function code(x1, x2)
	tmp = 0.0
	if ((x1 <= -2.6e+37) || !(x1 <= 1.75e+35))
		tmp = Float64(x1 + Float64(Float64(Float64(6.0 - Float64(3.0 / x1)) * Float64(x1 * x1)) * Float64(x1 * x1)));
	else
		tmp = Float64(x1 + fma(Float64(Float64(fma(-2.0, x2, Float64(Float64(x1 * x1) * 3.0)) - x1) / fma(x1, x1, 1.0)), 3.0, Float64(Float64(fma(fma(6.0, x1, -12.0), x1, Float64(8.0 * Float64(x2 * x1))) * x2) + x1)));
	end
	return tmp
end

Wolfram

code[x1_, x2_] := If[Or[LessEqual[x1, -2.6e+37], N[Not[LessEqual[x1, 1.75e+35]], $MachinePrecision]], N[(x1 + N[(N[(N[(6.0 - N[(3.0 / x1), $MachinePrecision]), $MachinePrecision] * N[(x1 * x1), $MachinePrecision]), $MachinePrecision] * N[(x1 * x1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x1 + N[(N[(N[(N[(-2.0 * x2 + N[(N[(x1 * x1), $MachinePrecision] * 3.0), $MachinePrecision]), $MachinePrecision] - x1), $MachinePrecision] / N[(x1 * x1 + 1.0), $MachinePrecision]), $MachinePrecision] * 3.0 + N[(N[(N[(N[(6.0 * x1 + -12.0), $MachinePrecision] * x1 + N[(8.0 * N[(x2 * x1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * x2), $MachinePrecision] + x1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]

Derivation

1. Split input into 2 regimes

2. if x1 < -2.5999999999999999e37 or 1.75e35 < x1

   1. Initial program 29.3%

      \[x1 + \left(\left(\left(\left(\left(\left(\left(2 \cdot x1\right) \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 + 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \cdot \left(\frac{\left(\left(3 \cdot x1\right) \cdot x1 + 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1} - 3\right) + \left(x1 \cdot x1\right) \cdot \left(4 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 + 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1} - 6\right)\right) \cdot \left(x1 \cdot x1 + 1\right) + \left(\left(3 \cdot x1\right) \cdot x1\right) \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 + 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) + \left(x1 \cdot x1\right) \cdot x1\right) + x1\right) + 3 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 - 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \]
   2. Add Preprocessing

   4. Applied rewrites 37.7%

      \[\leadsto x1 + \color{blue}{\mathsf{fma}\left(x1 \cdot x1, x1, \mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(4, \frac{\mathsf{fma}\left(x2, 2, \left(3 \cdot x1\right) \cdot x1\right) - x1}{\mathsf{fma}\left(x1, x1, 1\right)}, -6\right), x1 \cdot x1, \left(\frac{\mathsf{fma}\left(x2, 2, \left(3 \cdot x1\right) \cdot x1\right) - x1}{\mathsf{fma}\left(x1, x1, 1\right)} - 3\right) \cdot \left(\frac{\mathsf{fma}\left(x2, 2, \left(3 \cdot x1\right) \cdot x1\right) - x1}{\mathsf{fma}\left(x1, x1, 1\right)} \cdot \left(x1 \cdot 2\right)\right)\right), \mathsf{fma}\left(x1, x1, 1\right), \frac{\mathsf{fma}\left(x2, 2, \left(3 \cdot x1\right) \cdot x1\right) - x1}{\mathsf{fma}\left(x1, x1, 1\right)} \cdot \left(\left(3 \cdot x1\right) \cdot x1\right)\right) + \mathsf{fma}\left(\frac{\mathsf{fma}\left(-2, x2, \left(3 \cdot x1\right) \cdot x1\right) - x1}{\mathsf{fma}\left(x1, x1, 1\right)}, 3, x1\right)\right)} \]

   5. Taylor expanded in x1 around inf

      \[\leadsto x1 + \color{blue}{{x1}^{4} \cdot \left(6 - 3 \cdot \frac{1}{x1}\right)} \]
   6. Step-by-step derivation

      1. *-commutative N/A

         \[\leadsto x1 + \color{blue}{\left(6 - 3 \cdot \frac{1}{x1}\right) \cdot {x1}^{4}} \]

      2. lower-*.f64 N/A

         \[\leadsto x1 + \color{blue}{\left(6 - 3 \cdot \frac{1}{x1}\right) \cdot {x1}^{4}} \]

      3. lower--.f64 N/A

         \[\leadsto x1 + \color{blue}{\left(6 - 3 \cdot \frac{1}{x1}\right)} \cdot {x1}^{4} \]

      4. associate-*r/ N/A

         \[\leadsto x1 + \left(6 - \color{blue}{\frac{3 \cdot 1}{x1}}\right) \cdot {x1}^{4} \]

      5. metadata-eval N/A

         \[\leadsto x1 + \left(6 - \frac{\color{blue}{3}}{x1}\right) \cdot {x1}^{4} \]

      6. lower-/.f64 N/A

         \[\leadsto x1 + \left(6 - \color{blue}{\frac{3}{x1}}\right) \cdot {x1}^{4} \]

      7. lower-pow.f64 95.6

         \[\leadsto x1 + \left(6 - \frac{3}{x1}\right) \cdot \color{blue}{{x1}^{4}} \]

   7. Applied rewrites 95.6%

      \[\leadsto x1 + \color{blue}{\left(6 - \frac{3}{x1}\right) \cdot {x1}^{4}} \]
   8. Step-by-step derivation

   9. Applied rewrites 95.6%

      \[\leadsto x1 + \left(\left(6 - \frac{3}{x1}\right) \cdot \left(x1 \cdot x1\right)\right) \cdot \color{blue}{\left(x1 \cdot x1\right)} \]

   if -2.5999999999999999e37 < x1 < 1.75e35

   1. Initial program 99.3%

      \[x1 + \left(\left(\left(\left(\left(\left(\left(2 \cdot x1\right) \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 + 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \cdot \left(\frac{\left(\left(3 \cdot x1\right) \cdot x1 + 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1} - 3\right) + \left(x1 \cdot x1\right) \cdot \left(4 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 + 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1} - 6\right)\right) \cdot \left(x1 \cdot x1 + 1\right) + \left(\left(3 \cdot x1\right) \cdot x1\right) \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 + 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) + \left(x1 \cdot x1\right) \cdot x1\right) + x1\right) + 3 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 - 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \]
   2. Add Preprocessing

   3. Taylor expanded in x1 around 0

      \[\leadsto x1 + \left(\left(\color{blue}{x1 \cdot \left(4 \cdot \left(x2 \cdot \left(2 \cdot x2 - 3\right)\right) + x1 \cdot \left(\left(2 \cdot \left(-2 \cdot x2 + -1 \cdot \left(2 \cdot x2 - 3\right)\right) + \left(6 \cdot x2 + 8 \cdot x2\right)\right) - 6\right)\right)} + x1\right) + 3 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 - 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \]
   4. Step-by-step derivation

      1. *-commutative N/A

         \[\leadsto x1 + \left(\left(\color{blue}{\left(4 \cdot \left(x2 \cdot \left(2 \cdot x2 - 3\right)\right) + x1 \cdot \left(\left(2 \cdot \left(-2 \cdot x2 + -1 \cdot \left(2 \cdot x2 - 3\right)\right) + \left(6 \cdot x2 + 8 \cdot x2\right)\right) - 6\right)\right) \cdot x1} + x1\right) + 3 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 - 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \]

      2. lower-*.f64 N/A

         \[\leadsto x1 + \left(\left(\color{blue}{\left(4 \cdot \left(x2 \cdot \left(2 \cdot x2 - 3\right)\right) + x1 \cdot \left(\left(2 \cdot \left(-2 \cdot x2 + -1 \cdot \left(2 \cdot x2 - 3\right)\right) + \left(6 \cdot x2 + 8 \cdot x2\right)\right) - 6\right)\right) \cdot x1} + x1\right) + 3 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 - 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \]

   5. Applied rewrites 84.9%

      \[\leadsto x1 + \left(\left(\color{blue}{\mathsf{fma}\left(4 \cdot x2, \mathsf{fma}\left(2, x2, -3\right), \mathsf{fma}\left(-4, x2, \mathsf{fma}\left(\mathsf{fma}\left(-2, x2, 3\right), 2, \mathsf{fma}\left(14, x2, -6\right)\right)\right) \cdot x1\right) \cdot x1} + x1\right) + 3 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 - 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \]

   6. Taylor expanded in x2 around 0

      \[\leadsto x1 + \left(\left(\left(x2 \cdot \left(6 \cdot x1 - 12\right)\right) \cdot x1 + x1\right) + 3 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 - 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \]
   7. Step-by-step derivation

   8. Applied rewrites 73.8%

      \[\leadsto x1 + \left(\left(\left(\mathsf{fma}\left(6, x1, -12\right) \cdot x2\right) \cdot x1 + x1\right) + 3 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 - 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \]

   9. Taylor expanded in x2 around 0

      \[\leadsto x1 + \left(\left(x2 \cdot \color{blue}{\left(8 \cdot \left(x1 \cdot x2\right) + x1 \cdot \left(6 \cdot x1 - 12\right)\right)} + x1\right) + 3 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 - 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \]
   10. Step-by-step derivation

   11. Applied rewrites 94.0%

       \[\leadsto x1 + \left(\left(\mathsf{fma}\left(\mathsf{fma}\left(6, x1, -12\right), x1, \left(x2 \cdot x1\right) \cdot 8\right) \cdot \color{blue}{x2} + x1\right) + 3 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 - 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \]

   13. Applied rewrites 94.3%

       \[\leadsto x1 + \color{blue}{\mathsf{fma}\left(\frac{\mathsf{fma}\left(-2, x2, \left(x1 \cdot x1\right) \cdot 3\right) - x1}{\mathsf{fma}\left(x1, x1, 1\right)}, 3, \mathsf{fma}\left(\mathsf{fma}\left(6, x1, -12\right), x1, 8 \cdot \left(x2 \cdot x1\right)\right) \cdot x2 + x1\right)} \]
3. Recombined 2 regimes into one program.

4. Final simplification 94.9%

   \[\leadsto \begin{array}{l} \mathbf{if}\;x1 \leq -2.6 \cdot 10^{+37} \lor \neg \left(x1 \leq 1.75 \cdot 10^{+35}\right):\\ \;\;\;\;x1 + \left(\left(6 - \frac{3}{x1}\right) \cdot \left(x1 \cdot x1\right)\right) \cdot \left(x1 \cdot x1\right)\\ \mathbf{else}:\\ \;\;\;\;x1 + \mathsf{fma}\left(\frac{\mathsf{fma}\left(-2, x2, \left(x1 \cdot x1\right) \cdot 3\right) - x1}{\mathsf{fma}\left(x1, x1, 1\right)}, 3, \mathsf{fma}\left(\mathsf{fma}\left(6, x1, -12\right), x1, 8 \cdot \left(x2 \cdot x1\right)\right) \cdot x2 + x1\right)\\ \end{array} \]
5. Add Preprocessing

Alternative 9: 86.2% accurate, 3.7× speedup

Math

\[\begin{array}{l} \\ \begin{array}{l} \mathbf{if}\;x1 \leq -2.6 \cdot 10^{+37} \lor \neg \left(x1 \leq 1.75 \cdot 10^{+35}\right):\\ \;\;\;\;x1 + \left(\left(6 - \frac{3}{x1}\right) \cdot \left(x1 \cdot x1\right)\right) \cdot \left(x1 \cdot x1\right)\\ \mathbf{else}:\\ \;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(-4, x2, \mathsf{fma}\left(-4, x2, 6\right) + \mathsf{fma}\left(\mathsf{fma}\left(2, x2, 3\right), 3, \mathsf{fma}\left(14, x2, -6\right)\right)\right), x1, \mathsf{fma}\left(\mathsf{fma}\left(2, x2, -3\right) \cdot x2, 4, -1\right)\right), x1, -6 \cdot x2\right)\\ \end{array} \end{array} \]

FPCore

(FPCore (x1 x2)
 :precision binary64
 (if (or (<= x1 -2.6e+37) (not (<= x1 1.75e+35)))
   (+ x1 (* (* (- 6.0 (/ 3.0 x1)) (* x1 x1)) (* x1 x1)))
   (fma
    (fma
     (fma
      -4.0
      x2
      (+ (fma -4.0 x2 6.0) (fma (fma 2.0 x2 3.0) 3.0 (fma 14.0 x2 -6.0))))
     x1
     (fma (* (fma 2.0 x2 -3.0) x2) 4.0 -1.0))
    x1
    (* -6.0 x2))))

C

double code(double x1, double x2) {
	double tmp;
	if ((x1 <= -2.6e+37) || !(x1 <= 1.75e+35)) {
		tmp = x1 + (((6.0 - (3.0 / x1)) * (x1 * x1)) * (x1 * x1));
	} else {
		tmp = fma(fma(fma(-4.0, x2, (fma(-4.0, x2, 6.0) + fma(fma(2.0, x2, 3.0), 3.0, fma(14.0, x2, -6.0)))), x1, fma((fma(2.0, x2, -3.0) * x2), 4.0, -1.0)), x1, (-6.0 * x2));
	}
	return tmp;
}

Julia

function code(x1, x2)
	tmp = 0.0
	if ((x1 <= -2.6e+37) || !(x1 <= 1.75e+35))
		tmp = Float64(x1 + Float64(Float64(Float64(6.0 - Float64(3.0 / x1)) * Float64(x1 * x1)) * Float64(x1 * x1)));
	else
		tmp = fma(fma(fma(-4.0, x2, Float64(fma(-4.0, x2, 6.0) + fma(fma(2.0, x2, 3.0), 3.0, fma(14.0, x2, -6.0)))), x1, fma(Float64(fma(2.0, x2, -3.0) * x2), 4.0, -1.0)), x1, Float64(-6.0 * x2));
	end
	return tmp
end

Wolfram

code[x1_, x2_] := If[Or[LessEqual[x1, -2.6e+37], N[Not[LessEqual[x1, 1.75e+35]], $MachinePrecision]], N[(x1 + N[(N[(N[(6.0 - N[(3.0 / x1), $MachinePrecision]), $MachinePrecision] * N[(x1 * x1), $MachinePrecision]), $MachinePrecision] * N[(x1 * x1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(N[(-4.0 * x2 + N[(N[(-4.0 * x2 + 6.0), $MachinePrecision] + N[(N[(2.0 * x2 + 3.0), $MachinePrecision] * 3.0 + N[(14.0 * x2 + -6.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * x1 + N[(N[(N[(2.0 * x2 + -3.0), $MachinePrecision] * x2), $MachinePrecision] * 4.0 + -1.0), $MachinePrecision]), $MachinePrecision] * x1 + N[(-6.0 * x2), $MachinePrecision]), $MachinePrecision]]

Derivation

1. Split input into 2 regimes

2. if x1 < -2.5999999999999999e37 or 1.75e35 < x1

   1. Initial program 29.3%

      \[x1 + \left(\left(\left(\left(\left(\left(\left(2 \cdot x1\right) \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 + 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \cdot \left(\frac{\left(\left(3 \cdot x1\right) \cdot x1 + 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1} - 3\right) + \left(x1 \cdot x1\right) \cdot \left(4 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 + 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1} - 6\right)\right) \cdot \left(x1 \cdot x1 + 1\right) + \left(\left(3 \cdot x1\right) \cdot x1\right) \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 + 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) + \left(x1 \cdot x1\right) \cdot x1\right) + x1\right) + 3 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 - 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \]
   2. Add Preprocessing

   4. Applied rewrites 37.7%

      \[\leadsto x1 + \color{blue}{\mathsf{fma}\left(x1 \cdot x1, x1, \mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(4, \frac{\mathsf{fma}\left(x2, 2, \left(3 \cdot x1\right) \cdot x1\right) - x1}{\mathsf{fma}\left(x1, x1, 1\right)}, -6\right), x1 \cdot x1, \left(\frac{\mathsf{fma}\left(x2, 2, \left(3 \cdot x1\right) \cdot x1\right) - x1}{\mathsf{fma}\left(x1, x1, 1\right)} - 3\right) \cdot \left(\frac{\mathsf{fma}\left(x2, 2, \left(3 \cdot x1\right) \cdot x1\right) - x1}{\mathsf{fma}\left(x1, x1, 1\right)} \cdot \left(x1 \cdot 2\right)\right)\right), \mathsf{fma}\left(x1, x1, 1\right), \frac{\mathsf{fma}\left(x2, 2, \left(3 \cdot x1\right) \cdot x1\right) - x1}{\mathsf{fma}\left(x1, x1, 1\right)} \cdot \left(\left(3 \cdot x1\right) \cdot x1\right)\right) + \mathsf{fma}\left(\frac{\mathsf{fma}\left(-2, x2, \left(3 \cdot x1\right) \cdot x1\right) - x1}{\mathsf{fma}\left(x1, x1, 1\right)}, 3, x1\right)\right)} \]

   5. Taylor expanded in x1 around inf

      \[\leadsto x1 + \color{blue}{{x1}^{4} \cdot \left(6 - 3 \cdot \frac{1}{x1}\right)} \]
   6. Step-by-step derivation

      1. *-commutative N/A

         \[\leadsto x1 + \color{blue}{\left(6 - 3 \cdot \frac{1}{x1}\right) \cdot {x1}^{4}} \]

      2. lower-*.f64 N/A

         \[\leadsto x1 + \color{blue}{\left(6 - 3 \cdot \frac{1}{x1}\right) \cdot {x1}^{4}} \]

      3. lower--.f64 N/A

         \[\leadsto x1 + \color{blue}{\left(6 - 3 \cdot \frac{1}{x1}\right)} \cdot {x1}^{4} \]

      4. associate-*r/ N/A

         \[\leadsto x1 + \left(6 - \color{blue}{\frac{3 \cdot 1}{x1}}\right) \cdot {x1}^{4} \]

      5. metadata-eval N/A

         \[\leadsto x1 + \left(6 - \frac{\color{blue}{3}}{x1}\right) \cdot {x1}^{4} \]

      6. lower-/.f64 N/A

         \[\leadsto x1 + \left(6 - \color{blue}{\frac{3}{x1}}\right) \cdot {x1}^{4} \]

      7. lower-pow.f64 95.6

         \[\leadsto x1 + \left(6 - \frac{3}{x1}\right) \cdot \color{blue}{{x1}^{4}} \]

   7. Applied rewrites 95.6%

      \[\leadsto x1 + \color{blue}{\left(6 - \frac{3}{x1}\right) \cdot {x1}^{4}} \]
   8. Step-by-step derivation

   9. Applied rewrites 95.6%

      \[\leadsto x1 + \left(\left(6 - \frac{3}{x1}\right) \cdot \left(x1 \cdot x1\right)\right) \cdot \color{blue}{\left(x1 \cdot x1\right)} \]

   if -2.5999999999999999e37 < x1 < 1.75e35

   1. Initial program 99.3%

      \[x1 + \left(\left(\left(\left(\left(\left(\left(2 \cdot x1\right) \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 + 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \cdot \left(\frac{\left(\left(3 \cdot x1\right) \cdot x1 + 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1} - 3\right) + \left(x1 \cdot x1\right) \cdot \left(4 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 + 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1} - 6\right)\right) \cdot \left(x1 \cdot x1 + 1\right) + \left(\left(3 \cdot x1\right) \cdot x1\right) \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 + 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) + \left(x1 \cdot x1\right) \cdot x1\right) + x1\right) + 3 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 - 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \]
   2. Add Preprocessing

   4. Applied rewrites 99.7%

      \[\leadsto x1 + \color{blue}{\mathsf{fma}\left(x1 \cdot x1, x1, \mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(4, \frac{\mathsf{fma}\left(x2, 2, \left(3 \cdot x1\right) \cdot x1\right) - x1}{\mathsf{fma}\left(x1, x1, 1\right)}, -6\right), x1 \cdot x1, \left(\frac{\mathsf{fma}\left(x2, 2, \left(3 \cdot x1\right) \cdot x1\right) - x1}{\mathsf{fma}\left(x1, x1, 1\right)} - 3\right) \cdot \left(\frac{\mathsf{fma}\left(x2, 2, \left(3 \cdot x1\right) \cdot x1\right) - x1}{\mathsf{fma}\left(x1, x1, 1\right)} \cdot \left(x1 \cdot 2\right)\right)\right), \mathsf{fma}\left(x1, x1, 1\right), \frac{\mathsf{fma}\left(x2, 2, \left(3 \cdot x1\right) \cdot x1\right) - x1}{\mathsf{fma}\left(x1, x1, 1\right)} \cdot \left(\left(3 \cdot x1\right) \cdot x1\right)\right) + \mathsf{fma}\left(\frac{\mathsf{fma}\left(-2, x2, \left(3 \cdot x1\right) \cdot x1\right) - x1}{\mathsf{fma}\left(x1, x1, 1\right)}, 3, x1\right)\right)} \]

   6. Applied rewrites 99.7%

      \[\leadsto x1 + \mathsf{fma}\left(x1 \cdot x1, x1, \mathsf{fma}\left(\color{blue}{\mathsf{fma}\left(\mathsf{fma}\left(\frac{\mathsf{fma}\left(x1 \cdot x1, 3, 2 \cdot x2\right) - x1}{\mathsf{fma}\left(x1, x1, 1\right)}, 4, -6\right) \cdot x1, x1, \left(\frac{\mathsf{fma}\left(x1 \cdot x1, 3, 2 \cdot x2\right) - x1}{\mathsf{fma}\left(x1, x1, 1\right)} - 3\right) \cdot \left(\left(x1 \cdot 2\right) \cdot \frac{\mathsf{fma}\left(x1 \cdot x1, 3, 2 \cdot x2\right) - x1}{\mathsf{fma}\left(x1, x1, 1\right)}\right)\right)}, \mathsf{fma}\left(x1, x1, 1\right), \frac{\mathsf{fma}\left(x2, 2, \left(3 \cdot x1\right) \cdot x1\right) - x1}{\mathsf{fma}\left(x1, x1, 1\right)} \cdot \left(\left(3 \cdot x1\right) \cdot x1\right)\right) + \mathsf{fma}\left(\frac{\mathsf{fma}\left(-2, x2, \left(3 \cdot x1\right) \cdot x1\right) - x1}{\mathsf{fma}\left(x1, x1, 1\right)}, 3, x1\right)\right) \]

   7. Taylor expanded in x1 around 0

      \[\leadsto \color{blue}{-6 \cdot x2 + x1 \cdot \left(\left(4 \cdot \left(x2 \cdot \left(2 \cdot x2 - 3\right)\right) + x1 \cdot \left(\left(2 \cdot \left(-2 \cdot x2 + -1 \cdot \left(2 \cdot x2 - 3\right)\right) + \left(3 \cdot \left(3 - -2 \cdot x2\right) + \left(6 \cdot x2 + 8 \cdot x2\right)\right)\right) - 6\right)\right) - 1\right)} \]

   9. Applied rewrites 85.4%

      \[\leadsto \color{blue}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(-4, x2, \mathsf{fma}\left(-4, x2, 6\right) + \mathsf{fma}\left(\mathsf{fma}\left(2, x2, 3\right), 3, \mathsf{fma}\left(14, x2, -6\right)\right)\right), x1, \mathsf{fma}\left(\mathsf{fma}\left(2, x2, -3\right) \cdot x2, 4, -1\right)\right), x1, -6 \cdot x2\right)} \]
3. Recombined 2 regimes into one program.

4. Final simplification 90.1%

   \[\leadsto \begin{array}{l} \mathbf{if}\;x1 \leq -2.6 \cdot 10^{+37} \lor \neg \left(x1 \leq 1.75 \cdot 10^{+35}\right):\\ \;\;\;\;x1 + \left(\left(6 - \frac{3}{x1}\right) \cdot \left(x1 \cdot x1\right)\right) \cdot \left(x1 \cdot x1\right)\\ \mathbf{else}:\\ \;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(-4, x2, \mathsf{fma}\left(-4, x2, 6\right) + \mathsf{fma}\left(\mathsf{fma}\left(2, x2, 3\right), 3, \mathsf{fma}\left(14, x2, -6\right)\right)\right), x1, \mathsf{fma}\left(\mathsf{fma}\left(2, x2, -3\right) \cdot x2, 4, -1\right)\right), x1, -6 \cdot x2\right)\\ \end{array} \]
5. Add Preprocessing

Alternative 10: 85.9% accurate, 6.0× speedup

Math

\[\begin{array}{l} \\ \begin{array}{l} \mathbf{if}\;x1 \leq -2.6 \cdot 10^{+37} \lor \neg \left(x1 \leq 1.75 \cdot 10^{+35}\right):\\ \;\;\;\;x1 + \left(\left(6 - \frac{3}{x1}\right) \cdot \left(x1 \cdot x1\right)\right) \cdot \left(x1 \cdot x1\right)\\ \mathbf{else}:\\ \;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(2, x2, -3\right) \cdot x2, 4, -1\right), x1, -6 \cdot x2\right)\\ \end{array} \end{array} \]

FPCore

(FPCore (x1 x2)
 :precision binary64
 (if (or (<= x1 -2.6e+37) (not (<= x1 1.75e+35)))
   (+ x1 (* (* (- 6.0 (/ 3.0 x1)) (* x1 x1)) (* x1 x1)))
   (fma (fma (* (fma 2.0 x2 -3.0) x2) 4.0 -1.0) x1 (* -6.0 x2))))

C

double code(double x1, double x2) {
	double tmp;
	if ((x1 <= -2.6e+37) || !(x1 <= 1.75e+35)) {
		tmp = x1 + (((6.0 - (3.0 / x1)) * (x1 * x1)) * (x1 * x1));
	} else {
		tmp = fma(fma((fma(2.0, x2, -3.0) * x2), 4.0, -1.0), x1, (-6.0 * x2));
	}
	return tmp;
}

Julia

function code(x1, x2)
	tmp = 0.0
	if ((x1 <= -2.6e+37) || !(x1 <= 1.75e+35))
		tmp = Float64(x1 + Float64(Float64(Float64(6.0 - Float64(3.0 / x1)) * Float64(x1 * x1)) * Float64(x1 * x1)));
	else
		tmp = fma(fma(Float64(fma(2.0, x2, -3.0) * x2), 4.0, -1.0), x1, Float64(-6.0 * x2));
	end
	return tmp
end

Wolfram

code[x1_, x2_] := If[Or[LessEqual[x1, -2.6e+37], N[Not[LessEqual[x1, 1.75e+35]], $MachinePrecision]], N[(x1 + N[(N[(N[(6.0 - N[(3.0 / x1), $MachinePrecision]), $MachinePrecision] * N[(x1 * x1), $MachinePrecision]), $MachinePrecision] * N[(x1 * x1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(N[(N[(2.0 * x2 + -3.0), $MachinePrecision] * x2), $MachinePrecision] * 4.0 + -1.0), $MachinePrecision] * x1 + N[(-6.0 * x2), $MachinePrecision]), $MachinePrecision]]

Derivation

1. Split input into 2 regimes

2. if x1 < -2.5999999999999999e37 or 1.75e35 < x1

   1. Initial program 29.3%

      \[x1 + \left(\left(\left(\left(\left(\left(\left(2 \cdot x1\right) \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 + 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \cdot \left(\frac{\left(\left(3 \cdot x1\right) \cdot x1 + 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1} - 3\right) + \left(x1 \cdot x1\right) \cdot \left(4 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 + 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1} - 6\right)\right) \cdot \left(x1 \cdot x1 + 1\right) + \left(\left(3 \cdot x1\right) \cdot x1\right) \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 + 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) + \left(x1 \cdot x1\right) \cdot x1\right) + x1\right) + 3 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 - 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \]
   2. Add Preprocessing

   4. Applied rewrites 37.7%

      \[\leadsto x1 + \color{blue}{\mathsf{fma}\left(x1 \cdot x1, x1, \mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(4, \frac{\mathsf{fma}\left(x2, 2, \left(3 \cdot x1\right) \cdot x1\right) - x1}{\mathsf{fma}\left(x1, x1, 1\right)}, -6\right), x1 \cdot x1, \left(\frac{\mathsf{fma}\left(x2, 2, \left(3 \cdot x1\right) \cdot x1\right) - x1}{\mathsf{fma}\left(x1, x1, 1\right)} - 3\right) \cdot \left(\frac{\mathsf{fma}\left(x2, 2, \left(3 \cdot x1\right) \cdot x1\right) - x1}{\mathsf{fma}\left(x1, x1, 1\right)} \cdot \left(x1 \cdot 2\right)\right)\right), \mathsf{fma}\left(x1, x1, 1\right), \frac{\mathsf{fma}\left(x2, 2, \left(3 \cdot x1\right) \cdot x1\right) - x1}{\mathsf{fma}\left(x1, x1, 1\right)} \cdot \left(\left(3 \cdot x1\right) \cdot x1\right)\right) + \mathsf{fma}\left(\frac{\mathsf{fma}\left(-2, x2, \left(3 \cdot x1\right) \cdot x1\right) - x1}{\mathsf{fma}\left(x1, x1, 1\right)}, 3, x1\right)\right)} \]

   5. Taylor expanded in x1 around inf

      \[\leadsto x1 + \color{blue}{{x1}^{4} \cdot \left(6 - 3 \cdot \frac{1}{x1}\right)} \]
   6. Step-by-step derivation

      1. *-commutative N/A

         \[\leadsto x1 + \color{blue}{\left(6 - 3 \cdot \frac{1}{x1}\right) \cdot {x1}^{4}} \]

      2. lower-*.f64 N/A

         \[\leadsto x1 + \color{blue}{\left(6 - 3 \cdot \frac{1}{x1}\right) \cdot {x1}^{4}} \]

      3. lower--.f64 N/A

         \[\leadsto x1 + \color{blue}{\left(6 - 3 \cdot \frac{1}{x1}\right)} \cdot {x1}^{4} \]

      4. associate-*r/ N/A

         \[\leadsto x1 + \left(6 - \color{blue}{\frac{3 \cdot 1}{x1}}\right) \cdot {x1}^{4} \]

      5. metadata-eval N/A

         \[\leadsto x1 + \left(6 - \frac{\color{blue}{3}}{x1}\right) \cdot {x1}^{4} \]

      6. lower-/.f64 N/A

         \[\leadsto x1 + \left(6 - \color{blue}{\frac{3}{x1}}\right) \cdot {x1}^{4} \]

      7. lower-pow.f64 95.6

         \[\leadsto x1 + \left(6 - \frac{3}{x1}\right) \cdot \color{blue}{{x1}^{4}} \]

   7. Applied rewrites 95.6%

      \[\leadsto x1 + \color{blue}{\left(6 - \frac{3}{x1}\right) \cdot {x1}^{4}} \]
   8. Step-by-step derivation

   9. Applied rewrites 95.6%

      \[\leadsto x1 + \left(\left(6 - \frac{3}{x1}\right) \cdot \left(x1 \cdot x1\right)\right) \cdot \color{blue}{\left(x1 \cdot x1\right)} \]

   if -2.5999999999999999e37 < x1 < 1.75e35

   1. Initial program 99.3%

      \[x1 + \left(\left(\left(\left(\left(\left(\left(2 \cdot x1\right) \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 + 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \cdot \left(\frac{\left(\left(3 \cdot x1\right) \cdot x1 + 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1} - 3\right) + \left(x1 \cdot x1\right) \cdot \left(4 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 + 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1} - 6\right)\right) \cdot \left(x1 \cdot x1 + 1\right) + \left(\left(3 \cdot x1\right) \cdot x1\right) \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 + 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) + \left(x1 \cdot x1\right) \cdot x1\right) + x1\right) + 3 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 - 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \]
   2. Add Preprocessing

   4. Applied rewrites 99.7%

      \[\leadsto x1 + \color{blue}{\mathsf{fma}\left(x1 \cdot x1, x1, \mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(4, \frac{\mathsf{fma}\left(x2, 2, \left(3 \cdot x1\right) \cdot x1\right) - x1}{\mathsf{fma}\left(x1, x1, 1\right)}, -6\right), x1 \cdot x1, \left(\frac{\mathsf{fma}\left(x2, 2, \left(3 \cdot x1\right) \cdot x1\right) - x1}{\mathsf{fma}\left(x1, x1, 1\right)} - 3\right) \cdot \left(\frac{\mathsf{fma}\left(x2, 2, \left(3 \cdot x1\right) \cdot x1\right) - x1}{\mathsf{fma}\left(x1, x1, 1\right)} \cdot \left(x1 \cdot 2\right)\right)\right), \mathsf{fma}\left(x1, x1, 1\right), \frac{\mathsf{fma}\left(x2, 2, \left(3 \cdot x1\right) \cdot x1\right) - x1}{\mathsf{fma}\left(x1, x1, 1\right)} \cdot \left(\left(3 \cdot x1\right) \cdot x1\right)\right) + \mathsf{fma}\left(\frac{\mathsf{fma}\left(-2, x2, \left(3 \cdot x1\right) \cdot x1\right) - x1}{\mathsf{fma}\left(x1, x1, 1\right)}, 3, x1\right)\right)} \]

   5. Taylor expanded in x1 around 0

      \[\leadsto \color{blue}{-6 \cdot x2 + x1 \cdot \left(4 \cdot \left(x2 \cdot \left(2 \cdot x2 - 3\right)\right) - 1\right)} \]
   6. Step-by-step derivation

      1. +-commutative N/A

         \[\leadsto \color{blue}{x1 \cdot \left(4 \cdot \left(x2 \cdot \left(2 \cdot x2 - 3\right)\right) - 1\right) + -6 \cdot x2} \]

      2. *-commutative N/A

         \[\leadsto \color{blue}{\left(4 \cdot \left(x2 \cdot \left(2 \cdot x2 - 3\right)\right) - 1\right) \cdot x1} + -6 \cdot x2 \]

      3. lower-fma.f64 N/A

         \[\leadsto \color{blue}{\mathsf{fma}\left(4 \cdot \left(x2 \cdot \left(2 \cdot x2 - 3\right)\right) - 1, x1, -6 \cdot x2\right)} \]

      4. sub-neg N/A

         \[\leadsto \mathsf{fma}\left(\color{blue}{4 \cdot \left(x2 \cdot \left(2 \cdot x2 - 3\right)\right) + \left(\mathsf{neg}\left(1\right)\right)}, x1, -6 \cdot x2\right) \]

      5. *-commutative N/A

         \[\leadsto \mathsf{fma}\left(\color{blue}{\left(x2 \cdot \left(2 \cdot x2 - 3\right)\right) \cdot 4} + \left(\mathsf{neg}\left(1\right)\right), x1, -6 \cdot x2\right) \]

      6. metadata-eval N/A

         \[\leadsto \mathsf{fma}\left(\left(x2 \cdot \left(2 \cdot x2 - 3\right)\right) \cdot 4 + \color{blue}{-1}, x1, -6 \cdot x2\right) \]

      7. lower-fma.f64 N/A

         \[\leadsto \mathsf{fma}\left(\color{blue}{\mathsf{fma}\left(x2 \cdot \left(2 \cdot x2 - 3\right), 4, -1\right)}, x1, -6 \cdot x2\right) \]

      8. *-commutative N/A

         \[\leadsto \mathsf{fma}\left(\mathsf{fma}\left(\color{blue}{\left(2 \cdot x2 - 3\right) \cdot x2}, 4, -1\right), x1, -6 \cdot x2\right) \]

      9. lower-*.f64 N/A

         \[\leadsto \mathsf{fma}\left(\mathsf{fma}\left(\color{blue}{\left(2 \cdot x2 - 3\right) \cdot x2}, 4, -1\right), x1, -6 \cdot x2\right) \]

      10. sub-neg N/A

          \[\leadsto \mathsf{fma}\left(\mathsf{fma}\left(\color{blue}{\left(2 \cdot x2 + \left(\mathsf{neg}\left(3\right)\right)\right)} \cdot x2, 4, -1\right), x1, -6 \cdot x2\right) \]

      11. lower-fma.f64 N/A

          \[\leadsto \mathsf{fma}\left(\mathsf{fma}\left(\color{blue}{\mathsf{fma}\left(2, x2, \mathsf{neg}\left(3\right)\right)} \cdot x2, 4, -1\right), x1, -6 \cdot x2\right) \]

      12. metadata-eval N/A

          \[\leadsto \mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(2, x2, \color{blue}{-3}\right) \cdot x2, 4, -1\right), x1, -6 \cdot x2\right) \]

      13. lower-*.f64 85.1

          \[\leadsto \mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(2, x2, -3\right) \cdot x2, 4, -1\right), x1, \color{blue}{-6 \cdot x2}\right) \]

   7. Applied rewrites 85.1%

      \[\leadsto \color{blue}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(2, x2, -3\right) \cdot x2, 4, -1\right), x1, -6 \cdot x2\right)} \]
3. Recombined 2 regimes into one program.

4. Final simplification 90.0%

   \[\leadsto \begin{array}{l} \mathbf{if}\;x1 \leq -2.6 \cdot 10^{+37} \lor \neg \left(x1 \leq 1.75 \cdot 10^{+35}\right):\\ \;\;\;\;x1 + \left(\left(6 - \frac{3}{x1}\right) \cdot \left(x1 \cdot x1\right)\right) \cdot \left(x1 \cdot x1\right)\\ \mathbf{else}:\\ \;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(2, x2, -3\right) \cdot x2, 4, -1\right), x1, -6 \cdot x2\right)\\ \end{array} \]
5. Add Preprocessing

Alternative 11: 64.8% accurate, 8.5× speedup

Math

\[\begin{array}{l} \\ \begin{array}{l} \mathbf{if}\;x1 \leq 3.5 \cdot 10^{+88}:\\ \;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(2, x2, -3\right) \cdot x2, 4, -1\right), x1, -6 \cdot x2\right)\\ \mathbf{else}:\\ \;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(x1, x1, 1\right), x1, -6 \cdot x2\right)\\ \end{array} \end{array} \]

FPCore

(FPCore (x1 x2)
 :precision binary64
 (if (<= x1 3.5e+88)
   (fma (fma (* (fma 2.0 x2 -3.0) x2) 4.0 -1.0) x1 (* -6.0 x2))
   (fma (fma x1 x1 1.0) x1 (* -6.0 x2))))

C

double code(double x1, double x2) {
	double tmp;
	if (x1 <= 3.5e+88) {
		tmp = fma(fma((fma(2.0, x2, -3.0) * x2), 4.0, -1.0), x1, (-6.0 * x2));
	} else {
		tmp = fma(fma(x1, x1, 1.0), x1, (-6.0 * x2));
	}
	return tmp;
}

Julia

function code(x1, x2)
	tmp = 0.0
	if (x1 <= 3.5e+88)
		tmp = fma(fma(Float64(fma(2.0, x2, -3.0) * x2), 4.0, -1.0), x1, Float64(-6.0 * x2));
	else
		tmp = fma(fma(x1, x1, 1.0), x1, Float64(-6.0 * x2));
	end
	return tmp
end

Wolfram

code[x1_, x2_] := If[LessEqual[x1, 3.5e+88], N[(N[(N[(N[(2.0 * x2 + -3.0), $MachinePrecision] * x2), $MachinePrecision] * 4.0 + -1.0), $MachinePrecision] * x1 + N[(-6.0 * x2), $MachinePrecision]), $MachinePrecision], N[(N[(x1 * x1 + 1.0), $MachinePrecision] * x1 + N[(-6.0 * x2), $MachinePrecision]), $MachinePrecision]]

Derivation

1. Split input into 2 regimes

2. if x1 < 3.4999999999999998e88

   1. Initial program 75.2%

      \[x1 + \left(\left(\left(\left(\left(\left(\left(2 \cdot x1\right) \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 + 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \cdot \left(\frac{\left(\left(3 \cdot x1\right) \cdot x1 + 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1} - 3\right) + \left(x1 \cdot x1\right) \cdot \left(4 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 + 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1} - 6\right)\right) \cdot \left(x1 \cdot x1 + 1\right) + \left(\left(3 \cdot x1\right) \cdot x1\right) \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 + 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) + \left(x1 \cdot x1\right) \cdot x1\right) + x1\right) + 3 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 - 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \]
   2. Add Preprocessing

   4. Applied rewrites 80.2%

      \[\leadsto x1 + \color{blue}{\mathsf{fma}\left(x1 \cdot x1, x1, \mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(4, \frac{\mathsf{fma}\left(x2, 2, \left(3 \cdot x1\right) \cdot x1\right) - x1}{\mathsf{fma}\left(x1, x1, 1\right)}, -6\right), x1 \cdot x1, \left(\frac{\mathsf{fma}\left(x2, 2, \left(3 \cdot x1\right) \cdot x1\right) - x1}{\mathsf{fma}\left(x1, x1, 1\right)} - 3\right) \cdot \left(\frac{\mathsf{fma}\left(x2, 2, \left(3 \cdot x1\right) \cdot x1\right) - x1}{\mathsf{fma}\left(x1, x1, 1\right)} \cdot \left(x1 \cdot 2\right)\right)\right), \mathsf{fma}\left(x1, x1, 1\right), \frac{\mathsf{fma}\left(x2, 2, \left(3 \cdot x1\right) \cdot x1\right) - x1}{\mathsf{fma}\left(x1, x1, 1\right)} \cdot \left(\left(3 \cdot x1\right) \cdot x1\right)\right) + \mathsf{fma}\left(\frac{\mathsf{fma}\left(-2, x2, \left(3 \cdot x1\right) \cdot x1\right) - x1}{\mathsf{fma}\left(x1, x1, 1\right)}, 3, x1\right)\right)} \]

   5. Taylor expanded in x1 around 0

      \[\leadsto \color{blue}{-6 \cdot x2 + x1 \cdot \left(4 \cdot \left(x2 \cdot \left(2 \cdot x2 - 3\right)\right) - 1\right)} \]
   6. Step-by-step derivation

      1. +-commutative N/A

         \[\leadsto \color{blue}{x1 \cdot \left(4 \cdot \left(x2 \cdot \left(2 \cdot x2 - 3\right)\right) - 1\right) + -6 \cdot x2} \]

      2. *-commutative N/A

         \[\leadsto \color{blue}{\left(4 \cdot \left(x2 \cdot \left(2 \cdot x2 - 3\right)\right) - 1\right) \cdot x1} + -6 \cdot x2 \]

      3. lower-fma.f64 N/A

         \[\leadsto \color{blue}{\mathsf{fma}\left(4 \cdot \left(x2 \cdot \left(2 \cdot x2 - 3\right)\right) - 1, x1, -6 \cdot x2\right)} \]

      4. sub-neg N/A

         \[\leadsto \mathsf{fma}\left(\color{blue}{4 \cdot \left(x2 \cdot \left(2 \cdot x2 - 3\right)\right) + \left(\mathsf{neg}\left(1\right)\right)}, x1, -6 \cdot x2\right) \]

      5. *-commutative N/A

         \[\leadsto \mathsf{fma}\left(\color{blue}{\left(x2 \cdot \left(2 \cdot x2 - 3\right)\right) \cdot 4} + \left(\mathsf{neg}\left(1\right)\right), x1, -6 \cdot x2\right) \]

      6. metadata-eval N/A

         \[\leadsto \mathsf{fma}\left(\left(x2 \cdot \left(2 \cdot x2 - 3\right)\right) \cdot 4 + \color{blue}{-1}, x1, -6 \cdot x2\right) \]

      7. lower-fma.f64 N/A

         \[\leadsto \mathsf{fma}\left(\color{blue}{\mathsf{fma}\left(x2 \cdot \left(2 \cdot x2 - 3\right), 4, -1\right)}, x1, -6 \cdot x2\right) \]

      8. *-commutative N/A

         \[\leadsto \mathsf{fma}\left(\mathsf{fma}\left(\color{blue}{\left(2 \cdot x2 - 3\right) \cdot x2}, 4, -1\right), x1, -6 \cdot x2\right) \]

      9. lower-*.f64 N/A

         \[\leadsto \mathsf{fma}\left(\mathsf{fma}\left(\color{blue}{\left(2 \cdot x2 - 3\right) \cdot x2}, 4, -1\right), x1, -6 \cdot x2\right) \]

      10. sub-neg N/A

          \[\leadsto \mathsf{fma}\left(\mathsf{fma}\left(\color{blue}{\left(2 \cdot x2 + \left(\mathsf{neg}\left(3\right)\right)\right)} \cdot x2, 4, -1\right), x1, -6 \cdot x2\right) \]

      11. lower-fma.f64 N/A

          \[\leadsto \mathsf{fma}\left(\mathsf{fma}\left(\color{blue}{\mathsf{fma}\left(2, x2, \mathsf{neg}\left(3\right)\right)} \cdot x2, 4, -1\right), x1, -6 \cdot x2\right) \]

      12. metadata-eval N/A

          \[\leadsto \mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(2, x2, \color{blue}{-3}\right) \cdot x2, 4, -1\right), x1, -6 \cdot x2\right) \]

      13. lower-*.f64 58.4

          \[\leadsto \mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(2, x2, -3\right) \cdot x2, 4, -1\right), x1, \color{blue}{-6 \cdot x2}\right) \]

   7. Applied rewrites 58.4%

      \[\leadsto \color{blue}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(2, x2, -3\right) \cdot x2, 4, -1\right), x1, -6 \cdot x2\right)} \]

   if 3.4999999999999998e88 < x1

   1. Initial program 28.3%

      \[x1 + \left(\left(\left(\left(\left(\left(\left(2 \cdot x1\right) \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 + 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \cdot \left(\frac{\left(\left(3 \cdot x1\right) \cdot x1 + 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1} - 3\right) + \left(x1 \cdot x1\right) \cdot \left(4 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 + 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1} - 6\right)\right) \cdot \left(x1 \cdot x1 + 1\right) + \left(\left(3 \cdot x1\right) \cdot x1\right) \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 + 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) + \left(x1 \cdot x1\right) \cdot x1\right) + x1\right) + 3 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 - 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \]
   2. Add Preprocessing

   4. Applied rewrites 28.3%

      \[\leadsto x1 + \color{blue}{\mathsf{fma}\left(x1 \cdot x1, x1, \mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(4, \frac{\mathsf{fma}\left(x2, 2, \left(3 \cdot x1\right) \cdot x1\right) - x1}{\mathsf{fma}\left(x1, x1, 1\right)}, -6\right), x1 \cdot x1, \left(\frac{\mathsf{fma}\left(x2, 2, \left(3 \cdot x1\right) \cdot x1\right) - x1}{\mathsf{fma}\left(x1, x1, 1\right)} - 3\right) \cdot \left(\frac{\mathsf{fma}\left(x2, 2, \left(3 \cdot x1\right) \cdot x1\right) - x1}{\mathsf{fma}\left(x1, x1, 1\right)} \cdot \left(x1 \cdot 2\right)\right)\right), \mathsf{fma}\left(x1, x1, 1\right), \frac{\mathsf{fma}\left(x2, 2, \left(3 \cdot x1\right) \cdot x1\right) - x1}{\mathsf{fma}\left(x1, x1, 1\right)} \cdot \left(\left(3 \cdot x1\right) \cdot x1\right)\right) + \mathsf{fma}\left(\frac{\mathsf{fma}\left(-2, x2, \left(3 \cdot x1\right) \cdot x1\right) - x1}{\mathsf{fma}\left(x1, x1, 1\right)}, 3, x1\right)\right)} \]

   5. Taylor expanded in x1 around 0

      \[\leadsto x1 + \mathsf{fma}\left(x1 \cdot x1, x1, \color{blue}{-6 \cdot x2}\right) \]
   6. Step-by-step derivation

      1. lower-*.f64 92.3

         \[\leadsto x1 + \mathsf{fma}\left(x1 \cdot x1, x1, \color{blue}{-6 \cdot x2}\right) \]

   7. Applied rewrites 92.3%

      \[\leadsto x1 + \mathsf{fma}\left(x1 \cdot x1, x1, \color{blue}{-6 \cdot x2}\right) \]
   8. Step-by-step derivation

      1. lift-+.f64 N/A

         \[\leadsto \color{blue}{x1 + \mathsf{fma}\left(x1 \cdot x1, x1, -6 \cdot x2\right)} \]

   9. Applied rewrites 92.3%

      \[\leadsto \color{blue}{\mathsf{fma}\left(\mathsf{fma}\left(x1, x1, 1\right), x1, -6 \cdot x2\right)} \]
3. Recombined 2 regimes into one program.
4. Add Preprocessing

Alternative 12: 41.7% accurate, 16.6× speedup

Math

\[\begin{array}{l} \\ \mathsf{fma}\left(\mathsf{fma}\left(x1, x1, 1\right), x1, -6 \cdot x2\right) \end{array} \]

FPCore

(FPCore (x1 x2) :precision binary64 (fma (fma x1 x1 1.0) x1 (* -6.0 x2)))

C

double code(double x1, double x2) {
	return fma(fma(x1, x1, 1.0), x1, (-6.0 * x2));
}

Julia

function code(x1, x2)
	return fma(fma(x1, x1, 1.0), x1, Float64(-6.0 * x2))
end

Wolfram

code[x1_, x2_] := N[(N[(x1 * x1 + 1.0), $MachinePrecision] * x1 + N[(-6.0 * x2), $MachinePrecision]), $MachinePrecision]

Derivation

1. Initial program 66.8%

   \[x1 + \left(\left(\left(\left(\left(\left(\left(2 \cdot x1\right) \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 + 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \cdot \left(\frac{\left(\left(3 \cdot x1\right) \cdot x1 + 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1} - 3\right) + \left(x1 \cdot x1\right) \cdot \left(4 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 + 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1} - 6\right)\right) \cdot \left(x1 \cdot x1 + 1\right) + \left(\left(3 \cdot x1\right) \cdot x1\right) \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 + 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) + \left(x1 \cdot x1\right) \cdot x1\right) + x1\right) + 3 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 - 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \]
2. Add Preprocessing

4. Applied rewrites 70.9%

   \[\leadsto x1 + \color{blue}{\mathsf{fma}\left(x1 \cdot x1, x1, \mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(4, \frac{\mathsf{fma}\left(x2, 2, \left(3 \cdot x1\right) \cdot x1\right) - x1}{\mathsf{fma}\left(x1, x1, 1\right)}, -6\right), x1 \cdot x1, \left(\frac{\mathsf{fma}\left(x2, 2, \left(3 \cdot x1\right) \cdot x1\right) - x1}{\mathsf{fma}\left(x1, x1, 1\right)} - 3\right) \cdot \left(\frac{\mathsf{fma}\left(x2, 2, \left(3 \cdot x1\right) \cdot x1\right) - x1}{\mathsf{fma}\left(x1, x1, 1\right)} \cdot \left(x1 \cdot 2\right)\right)\right), \mathsf{fma}\left(x1, x1, 1\right), \frac{\mathsf{fma}\left(x2, 2, \left(3 \cdot x1\right) \cdot x1\right) - x1}{\mathsf{fma}\left(x1, x1, 1\right)} \cdot \left(\left(3 \cdot x1\right) \cdot x1\right)\right) + \mathsf{fma}\left(\frac{\mathsf{fma}\left(-2, x2, \left(3 \cdot x1\right) \cdot x1\right) - x1}{\mathsf{fma}\left(x1, x1, 1\right)}, 3, x1\right)\right)} \]

5. Taylor expanded in x1 around 0

   \[\leadsto x1 + \mathsf{fma}\left(x1 \cdot x1, x1, \color{blue}{-6 \cdot x2}\right) \]
6. Step-by-step derivation

   1. lower-*.f64 41.6

      \[\leadsto x1 + \mathsf{fma}\left(x1 \cdot x1, x1, \color{blue}{-6 \cdot x2}\right) \]

7. Applied rewrites 41.6%

   \[\leadsto x1 + \mathsf{fma}\left(x1 \cdot x1, x1, \color{blue}{-6 \cdot x2}\right) \]
8. Step-by-step derivation

   1. lift-+.f64 N/A

      \[\leadsto \color{blue}{x1 + \mathsf{fma}\left(x1 \cdot x1, x1, -6 \cdot x2\right)} \]

9. Applied rewrites 41.6%

   \[\leadsto \color{blue}{\mathsf{fma}\left(\mathsf{fma}\left(x1, x1, 1\right), x1, -6 \cdot x2\right)} \]
10. Add Preprocessing

Alternative 13: 26.2% accurate, 33.1× speedup

Math

\[\begin{array}{l} \\ x1 + -6 \cdot x2 \end{array} \]

FPCore

(FPCore (x1 x2) :precision binary64 (+ x1 (* -6.0 x2)))

C

double code(double x1, double x2) {
	return x1 + (-6.0 * x2);
}

Fortran

real(8) function code(x1, x2)
    real(8), intent (in) :: x1
    real(8), intent (in) :: x2
    code = x1 + ((-6.0d0) * x2)
end function

Java

public static double code(double x1, double x2) {
	return x1 + (-6.0 * x2);
}

Python

def code(x1, x2):
	return x1 + (-6.0 * x2)

Julia

function code(x1, x2)
	return Float64(x1 + Float64(-6.0 * x2))
end

MATLAB

function tmp = code(x1, x2)
	tmp = x1 + (-6.0 * x2);
end

Wolfram

code[x1_, x2_] := N[(x1 + N[(-6.0 * x2), $MachinePrecision]), $MachinePrecision]

Derivation

1. Initial program 66.8%

   \[x1 + \left(\left(\left(\left(\left(\left(\left(2 \cdot x1\right) \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 + 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \cdot \left(\frac{\left(\left(3 \cdot x1\right) \cdot x1 + 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1} - 3\right) + \left(x1 \cdot x1\right) \cdot \left(4 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 + 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1} - 6\right)\right) \cdot \left(x1 \cdot x1 + 1\right) + \left(\left(3 \cdot x1\right) \cdot x1\right) \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 + 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) + \left(x1 \cdot x1\right) \cdot x1\right) + x1\right) + 3 \cdot \frac{\left(\left(3 \cdot x1\right) \cdot x1 - 2 \cdot x2\right) - x1}{x1 \cdot x1 + 1}\right) \]
2. Add Preprocessing

3. Taylor expanded in x1 around 0

   \[\leadsto x1 + \color{blue}{-6 \cdot x2} \]
4. Step-by-step derivation

   1. lower-*.f64 25.8

      \[\leadsto x1 + \color{blue}{-6 \cdot x2} \]

5. Applied rewrites 25.8%

   \[\leadsto x1 + \color{blue}{-6 \cdot x2} \]
6. Add Preprocessing

Reproduce

herbie shell --seed 2024315 
(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))))))