Statistics.Distribution.Beta:$centropy from math-functions-0.1.5.2

Percentage Accurate: 95.2% → 97.5%

Time: 10.1s

Alternatives: 19

Speedup: 0.5×

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

Specification

Math

\[\begin{array}{l} \\ \left(\left(x - \left(y - 1\right) \cdot z\right) - \left(t - 1\right) \cdot a\right) + \left(\left(y + t\right) - 2\right) \cdot b \end{array} \]

FPCore

(FPCore (x y z t a b)
 :precision binary64
 (+ (- (- x (* (- y 1.0) z)) (* (- t 1.0) a)) (* (- (+ y t) 2.0) b)))

C

double code(double x, double y, double z, double t, double a, double b) {
	return ((x - ((y - 1.0) * z)) - ((t - 1.0) * a)) + (((y + t) - 2.0) * b);
}

Fortran

real(8) function code(x, y, z, t, a, b)
    real(8), intent (in) :: x
    real(8), intent (in) :: y
    real(8), intent (in) :: z
    real(8), intent (in) :: t
    real(8), intent (in) :: a
    real(8), intent (in) :: b
    code = ((x - ((y - 1.0d0) * z)) - ((t - 1.0d0) * a)) + (((y + t) - 2.0d0) * b)
end function

Java

public static double code(double x, double y, double z, double t, double a, double b) {
	return ((x - ((y - 1.0) * z)) - ((t - 1.0) * a)) + (((y + t) - 2.0) * b);
}

Python

def code(x, y, z, t, a, b):
	return ((x - ((y - 1.0) * z)) - ((t - 1.0) * a)) + (((y + t) - 2.0) * b)

Julia

function code(x, y, z, t, a, b)
	return Float64(Float64(Float64(x - Float64(Float64(y - 1.0) * z)) - Float64(Float64(t - 1.0) * a)) + Float64(Float64(Float64(y + t) - 2.0) * b))
end

MATLAB

function tmp = code(x, y, z, t, a, b)
	tmp = ((x - ((y - 1.0) * z)) - ((t - 1.0) * a)) + (((y + t) - 2.0) * b);
end

Wolfram

code[x_, y_, z_, t_, a_, b_] := N[(N[(N[(x - N[(N[(y - 1.0), $MachinePrecision] * z), $MachinePrecision]), $MachinePrecision] - N[(N[(t - 1.0), $MachinePrecision] * a), $MachinePrecision]), $MachinePrecision] + N[(N[(N[(y + t), $MachinePrecision] - 2.0), $MachinePrecision] * b), $MachinePrecision]), $MachinePrecision]

Initial Program: 95.2% accurate, 1.0× speedup

Math

\[\begin{array}{l} \\ \left(\left(x - \left(y - 1\right) \cdot z\right) - \left(t - 1\right) \cdot a\right) + \left(\left(y + t\right) - 2\right) \cdot b \end{array} \]

FPCore

(FPCore (x y z t a b)
 :precision binary64
 (+ (- (- x (* (- y 1.0) z)) (* (- t 1.0) a)) (* (- (+ y t) 2.0) b)))

C

double code(double x, double y, double z, double t, double a, double b) {
	return ((x - ((y - 1.0) * z)) - ((t - 1.0) * a)) + (((y + t) - 2.0) * b);
}

Fortran

real(8) function code(x, y, z, t, a, b)
    real(8), intent (in) :: x
    real(8), intent (in) :: y
    real(8), intent (in) :: z
    real(8), intent (in) :: t
    real(8), intent (in) :: a
    real(8), intent (in) :: b
    code = ((x - ((y - 1.0d0) * z)) - ((t - 1.0d0) * a)) + (((y + t) - 2.0d0) * b)
end function

Java

public static double code(double x, double y, double z, double t, double a, double b) {
	return ((x - ((y - 1.0) * z)) - ((t - 1.0) * a)) + (((y + t) - 2.0) * b);
}

Python

def code(x, y, z, t, a, b):
	return ((x - ((y - 1.0) * z)) - ((t - 1.0) * a)) + (((y + t) - 2.0) * b)

Julia

function code(x, y, z, t, a, b)
	return Float64(Float64(Float64(x - Float64(Float64(y - 1.0) * z)) - Float64(Float64(t - 1.0) * a)) + Float64(Float64(Float64(y + t) - 2.0) * b))
end

MATLAB

function tmp = code(x, y, z, t, a, b)
	tmp = ((x - ((y - 1.0) * z)) - ((t - 1.0) * a)) + (((y + t) - 2.0) * b);
end

Wolfram

code[x_, y_, z_, t_, a_, b_] := N[(N[(N[(x - N[(N[(y - 1.0), $MachinePrecision] * z), $MachinePrecision]), $MachinePrecision] - N[(N[(t - 1.0), $MachinePrecision] * a), $MachinePrecision]), $MachinePrecision] + N[(N[(N[(y + t), $MachinePrecision] - 2.0), $MachinePrecision] * b), $MachinePrecision]), $MachinePrecision]

Alternative 1: 97.5% accurate, 0.5× speedup

Math

\[\begin{array}{l} \\ \begin{array}{l} t_1 := \left(\left(x - \left(y - 1\right) \cdot z\right) - \left(t - 1\right) \cdot a\right) + \left(\left(y + t\right) - 2\right) \cdot b\\ \mathbf{if}\;t\_1 \leq \infty:\\ \;\;\;\;t\_1\\ \mathbf{else}:\\ \;\;\;\;\left(\left(t + y\right) - 2\right) \cdot b\\ \end{array} \end{array} \]

FPCore

(FPCore (x y z t a b)
 :precision binary64
 (let* ((t_1
         (+ (- (- x (* (- y 1.0) z)) (* (- t 1.0) a)) (* (- (+ y t) 2.0) b))))
   (if (<= t_1 INFINITY) t_1 (* (- (+ t y) 2.0) b))))

C

double code(double x, double y, double z, double t, double a, double b) {
	double t_1 = ((x - ((y - 1.0) * z)) - ((t - 1.0) * a)) + (((y + t) - 2.0) * b);
	double tmp;
	if (t_1 <= ((double) INFINITY)) {
		tmp = t_1;
	} else {
		tmp = ((t + y) - 2.0) * b;
	}
	return tmp;
}

Java

public static double code(double x, double y, double z, double t, double a, double b) {
	double t_1 = ((x - ((y - 1.0) * z)) - ((t - 1.0) * a)) + (((y + t) - 2.0) * b);
	double tmp;
	if (t_1 <= Double.POSITIVE_INFINITY) {
		tmp = t_1;
	} else {
		tmp = ((t + y) - 2.0) * b;
	}
	return tmp;
}

Python

def code(x, y, z, t, a, b):
	t_1 = ((x - ((y - 1.0) * z)) - ((t - 1.0) * a)) + (((y + t) - 2.0) * b)
	tmp = 0
	if t_1 <= math.inf:
		tmp = t_1
	else:
		tmp = ((t + y) - 2.0) * b
	return tmp

Julia

function code(x, y, z, t, a, b)
	t_1 = Float64(Float64(Float64(x - Float64(Float64(y - 1.0) * z)) - Float64(Float64(t - 1.0) * a)) + Float64(Float64(Float64(y + t) - 2.0) * b))
	tmp = 0.0
	if (t_1 <= Inf)
		tmp = t_1;
	else
		tmp = Float64(Float64(Float64(t + y) - 2.0) * b);
	end
	return tmp
end

MATLAB

function tmp_2 = code(x, y, z, t, a, b)
	t_1 = ((x - ((y - 1.0) * z)) - ((t - 1.0) * a)) + (((y + t) - 2.0) * b);
	tmp = 0.0;
	if (t_1 <= Inf)
		tmp = t_1;
	else
		tmp = ((t + y) - 2.0) * b;
	end
	tmp_2 = tmp;
end

Wolfram

code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(N[(N[(x - N[(N[(y - 1.0), $MachinePrecision] * z), $MachinePrecision]), $MachinePrecision] - N[(N[(t - 1.0), $MachinePrecision] * a), $MachinePrecision]), $MachinePrecision] + N[(N[(N[(y + t), $MachinePrecision] - 2.0), $MachinePrecision] * b), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$1, Infinity], t$95$1, N[(N[(N[(t + y), $MachinePrecision] - 2.0), $MachinePrecision] * b), $MachinePrecision]]]

Derivation

1. Split input into 2 regimes

2. if (+.f64 (-.f64 (-.f64 x (*.f64 (-.f64 y #s(literal 1 binary64)) z)) (*.f64 (-.f64 t #s(literal 1 binary64)) a)) (*.f64 (-.f64 (+.f64 y t) #s(literal 2 binary64)) b)) < +inf.0

   1. Initial program 100.0%

      \[\left(\left(x - \left(y - 1\right) \cdot z\right) - \left(t - 1\right) \cdot a\right) + \left(\left(y + t\right) - 2\right) \cdot b \]
   2. Add Preprocessing

   if +inf.0 < (+.f64 (-.f64 (-.f64 x (*.f64 (-.f64 y #s(literal 1 binary64)) z)) (*.f64 (-.f64 t #s(literal 1 binary64)) a)) (*.f64 (-.f64 (+.f64 y t) #s(literal 2 binary64)) b))

   1. Initial program 0.0%

      \[\left(\left(x - \left(y - 1\right) \cdot z\right) - \left(t - 1\right) \cdot a\right) + \left(\left(y + t\right) - 2\right) \cdot b \]
   2. Add Preprocessing

   3. Taylor expanded in y around 0

      \[\leadsto \color{blue}{\left(x + b \cdot \left(t - 2\right)\right) - \left(-1 \cdot z + a \cdot \left(t - 1\right)\right)} \]
   4. Step-by-step derivation

      1. sub-neg N/A

         \[\leadsto \color{blue}{\left(x + b \cdot \left(t - 2\right)\right) + \left(\mathsf{neg}\left(\left(-1 \cdot z + a \cdot \left(t - 1\right)\right)\right)\right)} \]

      2. +-commutative N/A

         \[\leadsto \color{blue}{\left(\mathsf{neg}\left(\left(-1 \cdot z + a \cdot \left(t - 1\right)\right)\right)\right) + \left(x + b \cdot \left(t - 2\right)\right)} \]

      3. +-commutative N/A

         \[\leadsto \left(\mathsf{neg}\left(\color{blue}{\left(a \cdot \left(t - 1\right) + -1 \cdot z\right)}\right)\right) + \left(x + b \cdot \left(t - 2\right)\right) \]

      4. distribute-neg-in N/A

         \[\leadsto \color{blue}{\left(\left(\mathsf{neg}\left(a \cdot \left(t - 1\right)\right)\right) + \left(\mathsf{neg}\left(-1 \cdot z\right)\right)\right)} + \left(x + b \cdot \left(t - 2\right)\right) \]

      5. mul-1-neg N/A

         \[\leadsto \left(\left(\mathsf{neg}\left(a \cdot \left(t - 1\right)\right)\right) + \left(\mathsf{neg}\left(\color{blue}{\left(\mathsf{neg}\left(z\right)\right)}\right)\right)\right) + \left(x + b \cdot \left(t - 2\right)\right) \]

      6. remove-double-neg N/A

         \[\leadsto \left(\left(\mathsf{neg}\left(a \cdot \left(t - 1\right)\right)\right) + \color{blue}{z}\right) + \left(x + b \cdot \left(t - 2\right)\right) \]

      7. associate-+l+ N/A

         \[\leadsto \color{blue}{\left(\mathsf{neg}\left(a \cdot \left(t - 1\right)\right)\right) + \left(z + \left(x + b \cdot \left(t - 2\right)\right)\right)} \]

      8. *-commutative N/A

         \[\leadsto \left(\mathsf{neg}\left(\color{blue}{\left(t - 1\right) \cdot a}\right)\right) + \left(z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      9. distribute-lft-neg-in N/A

         \[\leadsto \color{blue}{\left(\mathsf{neg}\left(\left(t - 1\right)\right)\right) \cdot a} + \left(z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      10. mul-1-neg N/A

          \[\leadsto \color{blue}{\left(-1 \cdot \left(t - 1\right)\right)} \cdot a + \left(z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      11. lower-fma.f64 N/A

          \[\leadsto \color{blue}{\mathsf{fma}\left(-1 \cdot \left(t - 1\right), a, z + \left(x + b \cdot \left(t - 2\right)\right)\right)} \]

      12. sub-neg N/A

          \[\leadsto \mathsf{fma}\left(-1 \cdot \color{blue}{\left(t + \left(\mathsf{neg}\left(1\right)\right)\right)}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      13. metadata-eval N/A

          \[\leadsto \mathsf{fma}\left(-1 \cdot \left(t + \color{blue}{-1}\right), a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      14. distribute-lft-in N/A

          \[\leadsto \mathsf{fma}\left(\color{blue}{-1 \cdot t + -1 \cdot -1}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      15. metadata-eval N/A

          \[\leadsto \mathsf{fma}\left(-1 \cdot t + \color{blue}{1}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      16. +-commutative N/A

          \[\leadsto \mathsf{fma}\left(\color{blue}{1 + -1 \cdot t}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      17. neg-mul-1 N/A

          \[\leadsto \mathsf{fma}\left(1 + \color{blue}{\left(\mathsf{neg}\left(t\right)\right)}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      18. sub-neg N/A

          \[\leadsto \mathsf{fma}\left(\color{blue}{1 - t}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      19. lower--.f64 N/A

          \[\leadsto \mathsf{fma}\left(\color{blue}{1 - t}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      20. lower-+.f64 N/A

          \[\leadsto \mathsf{fma}\left(1 - t, a, \color{blue}{z + \left(x + b \cdot \left(t - 2\right)\right)}\right) \]

   5. Applied rewrites 67.2%

      \[\leadsto \color{blue}{\mathsf{fma}\left(1 - t, a, z + \mathsf{fma}\left(t - 2, b, x\right)\right)} \]

   6. Taylor expanded in b around inf

      \[\leadsto \color{blue}{b \cdot \left(\left(t + y\right) - 2\right)} \]
   7. Step-by-step derivation

      1. *-commutative N/A

         \[\leadsto \color{blue}{\left(\left(t + y\right) - 2\right) \cdot b} \]

      2. lower-*.f64 N/A

         \[\leadsto \color{blue}{\left(\left(t + y\right) - 2\right) \cdot b} \]

      3. lower--.f64 N/A

         \[\leadsto \color{blue}{\left(\left(t + y\right) - 2\right)} \cdot b \]

      4. lower-+.f64 75.4

         \[\leadsto \left(\color{blue}{\left(t + y\right)} - 2\right) \cdot b \]

   8. Applied rewrites 75.4%

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

Alternative 2: 84.4% accurate, 0.9× speedup

Math

\[\begin{array}{l} \\ \begin{array}{l} t_1 := \left(b - z\right) \cdot y\\ \mathbf{if}\;y \leq -3.45 \cdot 10^{+124}:\\ \;\;\;\;t\_1\\ \mathbf{elif}\;y \leq -2.1 \cdot 10^{+30}:\\ \;\;\;\;x - \mathsf{fma}\left(a, t - 1, z \cdot \left(y - 1\right)\right)\\ \mathbf{elif}\;y \leq 2.95 \cdot 10^{+99}:\\ \;\;\;\;\mathsf{fma}\left(1 - t, a, z + \mathsf{fma}\left(t - 2, b, x\right)\right)\\ \mathbf{else}:\\ \;\;\;\;t\_1\\ \end{array} \end{array} \]

FPCore

(FPCore (x y z t a b)
 :precision binary64
 (let* ((t_1 (* (- b z) y)))
   (if (<= y -3.45e+124)
     t_1
     (if (<= y -2.1e+30)
       (- x (fma a (- t 1.0) (* z (- y 1.0))))
       (if (<= y 2.95e+99) (fma (- 1.0 t) a (+ z (fma (- t 2.0) b x))) t_1)))))

C

double code(double x, double y, double z, double t, double a, double b) {
	double t_1 = (b - z) * y;
	double tmp;
	if (y <= -3.45e+124) {
		tmp = t_1;
	} else if (y <= -2.1e+30) {
		tmp = x - fma(a, (t - 1.0), (z * (y - 1.0)));
	} else if (y <= 2.95e+99) {
		tmp = fma((1.0 - t), a, (z + fma((t - 2.0), b, x)));
	} else {
		tmp = t_1;
	}
	return tmp;
}

Julia

function code(x, y, z, t, a, b)
	t_1 = Float64(Float64(b - z) * y)
	tmp = 0.0
	if (y <= -3.45e+124)
		tmp = t_1;
	elseif (y <= -2.1e+30)
		tmp = Float64(x - fma(a, Float64(t - 1.0), Float64(z * Float64(y - 1.0))));
	elseif (y <= 2.95e+99)
		tmp = fma(Float64(1.0 - t), a, Float64(z + fma(Float64(t - 2.0), b, x)));
	else
		tmp = t_1;
	end
	return tmp
end

Wolfram

code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(N[(b - z), $MachinePrecision] * y), $MachinePrecision]}, If[LessEqual[y, -3.45e+124], t$95$1, If[LessEqual[y, -2.1e+30], N[(x - N[(a * N[(t - 1.0), $MachinePrecision] + N[(z * N[(y - 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 2.95e+99], N[(N[(1.0 - t), $MachinePrecision] * a + N[(z + N[(N[(t - 2.0), $MachinePrecision] * b + x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$1]]]]

Derivation

1. Split input into 3 regimes

2. if y < -3.45e124 or 2.9499999999999999e99 < y

   1. Initial program 90.8%

      \[\left(\left(x - \left(y - 1\right) \cdot z\right) - \left(t - 1\right) \cdot a\right) + \left(\left(y + t\right) - 2\right) \cdot b \]
   2. Add Preprocessing

   3. Taylor expanded in y around inf

      \[\leadsto \color{blue}{y \cdot \left(b - z\right)} \]
   4. Step-by-step derivation

      1. *-commutative N/A

         \[\leadsto \color{blue}{\left(b - z\right) \cdot y} \]

      2. lower-*.f64 N/A

         \[\leadsto \color{blue}{\left(b - z\right) \cdot y} \]

      3. lower--.f64 75.5

         \[\leadsto \color{blue}{\left(b - z\right)} \cdot y \]

   5. Applied rewrites 75.5%

      \[\leadsto \color{blue}{\left(b - z\right) \cdot y} \]

   if -3.45e124 < y < -2.1e30

   1. Initial program 100.0%

      \[\left(\left(x - \left(y - 1\right) \cdot z\right) - \left(t - 1\right) \cdot a\right) + \left(\left(y + t\right) - 2\right) \cdot b \]
   2. Add Preprocessing

   3. Taylor expanded in y around 0

      \[\leadsto \color{blue}{\left(x + b \cdot \left(t - 2\right)\right) - \left(-1 \cdot z + a \cdot \left(t - 1\right)\right)} \]
   4. Step-by-step derivation

      1. sub-neg N/A

         \[\leadsto \color{blue}{\left(x + b \cdot \left(t - 2\right)\right) + \left(\mathsf{neg}\left(\left(-1 \cdot z + a \cdot \left(t - 1\right)\right)\right)\right)} \]

      2. +-commutative N/A

         \[\leadsto \color{blue}{\left(\mathsf{neg}\left(\left(-1 \cdot z + a \cdot \left(t - 1\right)\right)\right)\right) + \left(x + b \cdot \left(t - 2\right)\right)} \]

      3. +-commutative N/A

         \[\leadsto \left(\mathsf{neg}\left(\color{blue}{\left(a \cdot \left(t - 1\right) + -1 \cdot z\right)}\right)\right) + \left(x + b \cdot \left(t - 2\right)\right) \]

      4. distribute-neg-in N/A

         \[\leadsto \color{blue}{\left(\left(\mathsf{neg}\left(a \cdot \left(t - 1\right)\right)\right) + \left(\mathsf{neg}\left(-1 \cdot z\right)\right)\right)} + \left(x + b \cdot \left(t - 2\right)\right) \]

      5. mul-1-neg N/A

         \[\leadsto \left(\left(\mathsf{neg}\left(a \cdot \left(t - 1\right)\right)\right) + \left(\mathsf{neg}\left(\color{blue}{\left(\mathsf{neg}\left(z\right)\right)}\right)\right)\right) + \left(x + b \cdot \left(t - 2\right)\right) \]

      6. remove-double-neg N/A

         \[\leadsto \left(\left(\mathsf{neg}\left(a \cdot \left(t - 1\right)\right)\right) + \color{blue}{z}\right) + \left(x + b \cdot \left(t - 2\right)\right) \]

      7. associate-+l+ N/A

         \[\leadsto \color{blue}{\left(\mathsf{neg}\left(a \cdot \left(t - 1\right)\right)\right) + \left(z + \left(x + b \cdot \left(t - 2\right)\right)\right)} \]

      8. *-commutative N/A

         \[\leadsto \left(\mathsf{neg}\left(\color{blue}{\left(t - 1\right) \cdot a}\right)\right) + \left(z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      9. distribute-lft-neg-in N/A

         \[\leadsto \color{blue}{\left(\mathsf{neg}\left(\left(t - 1\right)\right)\right) \cdot a} + \left(z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      10. mul-1-neg N/A

          \[\leadsto \color{blue}{\left(-1 \cdot \left(t - 1\right)\right)} \cdot a + \left(z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      11. lower-fma.f64 N/A

          \[\leadsto \color{blue}{\mathsf{fma}\left(-1 \cdot \left(t - 1\right), a, z + \left(x + b \cdot \left(t - 2\right)\right)\right)} \]

      12. sub-neg N/A

          \[\leadsto \mathsf{fma}\left(-1 \cdot \color{blue}{\left(t + \left(\mathsf{neg}\left(1\right)\right)\right)}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      13. metadata-eval N/A

          \[\leadsto \mathsf{fma}\left(-1 \cdot \left(t + \color{blue}{-1}\right), a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      14. distribute-lft-in N/A

          \[\leadsto \mathsf{fma}\left(\color{blue}{-1 \cdot t + -1 \cdot -1}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      15. metadata-eval N/A

          \[\leadsto \mathsf{fma}\left(-1 \cdot t + \color{blue}{1}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      16. +-commutative N/A

          \[\leadsto \mathsf{fma}\left(\color{blue}{1 + -1 \cdot t}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      17. neg-mul-1 N/A

          \[\leadsto \mathsf{fma}\left(1 + \color{blue}{\left(\mathsf{neg}\left(t\right)\right)}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      18. sub-neg N/A

          \[\leadsto \mathsf{fma}\left(\color{blue}{1 - t}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      19. lower--.f64 N/A

          \[\leadsto \mathsf{fma}\left(\color{blue}{1 - t}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      20. lower-+.f64 N/A

          \[\leadsto \mathsf{fma}\left(1 - t, a, \color{blue}{z + \left(x + b \cdot \left(t - 2\right)\right)}\right) \]

   5. Applied rewrites 70.8%

      \[\leadsto \color{blue}{\mathsf{fma}\left(1 - t, a, z + \mathsf{fma}\left(t - 2, b, x\right)\right)} \]

   6. Taylor expanded in b around 0

      \[\leadsto \color{blue}{x - \left(a \cdot \left(t - 1\right) + z \cdot \left(y - 1\right)\right)} \]
   7. Step-by-step derivation

      1. lower--.f64 N/A

         \[\leadsto \color{blue}{x - \left(a \cdot \left(t - 1\right) + z \cdot \left(y - 1\right)\right)} \]

      2. lower-fma.f64 N/A

         \[\leadsto x - \color{blue}{\mathsf{fma}\left(a, t - 1, z \cdot \left(y - 1\right)\right)} \]

      3. lower--.f64 N/A

         \[\leadsto x - \mathsf{fma}\left(a, \color{blue}{t - 1}, z \cdot \left(y - 1\right)\right) \]

      4. lower-*.f64 N/A

         \[\leadsto x - \mathsf{fma}\left(a, t - 1, \color{blue}{z \cdot \left(y - 1\right)}\right) \]

      5. lower--.f64 91.4

         \[\leadsto x - \mathsf{fma}\left(a, t - 1, z \cdot \color{blue}{\left(y - 1\right)}\right) \]

   8. Applied rewrites 91.4%

      \[\leadsto \color{blue}{x - \mathsf{fma}\left(a, t - 1, z \cdot \left(y - 1\right)\right)} \]

   if -2.1e30 < y < 2.9499999999999999e99

   1. Initial program 96.8%

      \[\left(\left(x - \left(y - 1\right) \cdot z\right) - \left(t - 1\right) \cdot a\right) + \left(\left(y + t\right) - 2\right) \cdot b \]
   2. Add Preprocessing

   3. Taylor expanded in y around 0

      \[\leadsto \color{blue}{\left(x + b \cdot \left(t - 2\right)\right) - \left(-1 \cdot z + a \cdot \left(t - 1\right)\right)} \]
   4. Step-by-step derivation

      1. sub-neg N/A

         \[\leadsto \color{blue}{\left(x + b \cdot \left(t - 2\right)\right) + \left(\mathsf{neg}\left(\left(-1 \cdot z + a \cdot \left(t - 1\right)\right)\right)\right)} \]

      2. +-commutative N/A

         \[\leadsto \color{blue}{\left(\mathsf{neg}\left(\left(-1 \cdot z + a \cdot \left(t - 1\right)\right)\right)\right) + \left(x + b \cdot \left(t - 2\right)\right)} \]

      3. +-commutative N/A

         \[\leadsto \left(\mathsf{neg}\left(\color{blue}{\left(a \cdot \left(t - 1\right) + -1 \cdot z\right)}\right)\right) + \left(x + b \cdot \left(t - 2\right)\right) \]

      4. distribute-neg-in N/A

         \[\leadsto \color{blue}{\left(\left(\mathsf{neg}\left(a \cdot \left(t - 1\right)\right)\right) + \left(\mathsf{neg}\left(-1 \cdot z\right)\right)\right)} + \left(x + b \cdot \left(t - 2\right)\right) \]

      5. mul-1-neg N/A

         \[\leadsto \left(\left(\mathsf{neg}\left(a \cdot \left(t - 1\right)\right)\right) + \left(\mathsf{neg}\left(\color{blue}{\left(\mathsf{neg}\left(z\right)\right)}\right)\right)\right) + \left(x + b \cdot \left(t - 2\right)\right) \]

      6. remove-double-neg N/A

         \[\leadsto \left(\left(\mathsf{neg}\left(a \cdot \left(t - 1\right)\right)\right) + \color{blue}{z}\right) + \left(x + b \cdot \left(t - 2\right)\right) \]

      7. associate-+l+ N/A

         \[\leadsto \color{blue}{\left(\mathsf{neg}\left(a \cdot \left(t - 1\right)\right)\right) + \left(z + \left(x + b \cdot \left(t - 2\right)\right)\right)} \]

      8. *-commutative N/A

         \[\leadsto \left(\mathsf{neg}\left(\color{blue}{\left(t - 1\right) \cdot a}\right)\right) + \left(z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      9. distribute-lft-neg-in N/A

         \[\leadsto \color{blue}{\left(\mathsf{neg}\left(\left(t - 1\right)\right)\right) \cdot a} + \left(z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      10. mul-1-neg N/A

          \[\leadsto \color{blue}{\left(-1 \cdot \left(t - 1\right)\right)} \cdot a + \left(z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      11. lower-fma.f64 N/A

          \[\leadsto \color{blue}{\mathsf{fma}\left(-1 \cdot \left(t - 1\right), a, z + \left(x + b \cdot \left(t - 2\right)\right)\right)} \]

      12. sub-neg N/A

          \[\leadsto \mathsf{fma}\left(-1 \cdot \color{blue}{\left(t + \left(\mathsf{neg}\left(1\right)\right)\right)}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      13. metadata-eval N/A

          \[\leadsto \mathsf{fma}\left(-1 \cdot \left(t + \color{blue}{-1}\right), a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      14. distribute-lft-in N/A

          \[\leadsto \mathsf{fma}\left(\color{blue}{-1 \cdot t + -1 \cdot -1}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      15. metadata-eval N/A

          \[\leadsto \mathsf{fma}\left(-1 \cdot t + \color{blue}{1}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      16. +-commutative N/A

          \[\leadsto \mathsf{fma}\left(\color{blue}{1 + -1 \cdot t}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      17. neg-mul-1 N/A

          \[\leadsto \mathsf{fma}\left(1 + \color{blue}{\left(\mathsf{neg}\left(t\right)\right)}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      18. sub-neg N/A

          \[\leadsto \mathsf{fma}\left(\color{blue}{1 - t}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      19. lower--.f64 N/A

          \[\leadsto \mathsf{fma}\left(\color{blue}{1 - t}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      20. lower-+.f64 N/A

          \[\leadsto \mathsf{fma}\left(1 - t, a, \color{blue}{z + \left(x + b \cdot \left(t - 2\right)\right)}\right) \]

   5. Applied rewrites 96.8%

      \[\leadsto \color{blue}{\mathsf{fma}\left(1 - t, a, z + \mathsf{fma}\left(t - 2, b, x\right)\right)} \]
3. Recombined 3 regimes into one program.
4. Add Preprocessing

Alternative 3: 87.9% accurate, 1.1× speedup

Math

\[\begin{array}{l} \\ \begin{array}{l} \mathbf{if}\;t \leq -7800000 \lor \neg \left(t \leq 4.3 \cdot 10^{+42}\right):\\ \;\;\;\;\mathsf{fma}\left(1 - t, a, z + \mathsf{fma}\left(t - 2, b, x\right)\right)\\ \mathbf{else}:\\ \;\;\;\;\mathsf{fma}\left(1 - y, z, a + \mathsf{fma}\left(y - 2, b, x\right)\right)\\ \end{array} \end{array} \]

FPCore

(FPCore (x y z t a b)
 :precision binary64
 (if (or (<= t -7800000.0) (not (<= t 4.3e+42)))
   (fma (- 1.0 t) a (+ z (fma (- t 2.0) b x)))
   (fma (- 1.0 y) z (+ a (fma (- y 2.0) b x)))))

C

double code(double x, double y, double z, double t, double a, double b) {
	double tmp;
	if ((t <= -7800000.0) || !(t <= 4.3e+42)) {
		tmp = fma((1.0 - t), a, (z + fma((t - 2.0), b, x)));
	} else {
		tmp = fma((1.0 - y), z, (a + fma((y - 2.0), b, x)));
	}
	return tmp;
}

Julia

function code(x, y, z, t, a, b)
	tmp = 0.0
	if ((t <= -7800000.0) || !(t <= 4.3e+42))
		tmp = fma(Float64(1.0 - t), a, Float64(z + fma(Float64(t - 2.0), b, x)));
	else
		tmp = fma(Float64(1.0 - y), z, Float64(a + fma(Float64(y - 2.0), b, x)));
	end
	return tmp
end

Wolfram

code[x_, y_, z_, t_, a_, b_] := If[Or[LessEqual[t, -7800000.0], N[Not[LessEqual[t, 4.3e+42]], $MachinePrecision]], N[(N[(1.0 - t), $MachinePrecision] * a + N[(z + N[(N[(t - 2.0), $MachinePrecision] * b + x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(1.0 - y), $MachinePrecision] * z + N[(a + N[(N[(y - 2.0), $MachinePrecision] * b + x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]

Derivation

1. Split input into 2 regimes

2. if t < -7.8e6 or 4.2999999999999998e42 < t

   1. Initial program 92.8%

      \[\left(\left(x - \left(y - 1\right) \cdot z\right) - \left(t - 1\right) \cdot a\right) + \left(\left(y + t\right) - 2\right) \cdot b \]
   2. Add Preprocessing

   3. Taylor expanded in y around 0

      \[\leadsto \color{blue}{\left(x + b \cdot \left(t - 2\right)\right) - \left(-1 \cdot z + a \cdot \left(t - 1\right)\right)} \]
   4. Step-by-step derivation

      1. sub-neg N/A

         \[\leadsto \color{blue}{\left(x + b \cdot \left(t - 2\right)\right) + \left(\mathsf{neg}\left(\left(-1 \cdot z + a \cdot \left(t - 1\right)\right)\right)\right)} \]

      2. +-commutative N/A

         \[\leadsto \color{blue}{\left(\mathsf{neg}\left(\left(-1 \cdot z + a \cdot \left(t - 1\right)\right)\right)\right) + \left(x + b \cdot \left(t - 2\right)\right)} \]

      3. +-commutative N/A

         \[\leadsto \left(\mathsf{neg}\left(\color{blue}{\left(a \cdot \left(t - 1\right) + -1 \cdot z\right)}\right)\right) + \left(x + b \cdot \left(t - 2\right)\right) \]

      4. distribute-neg-in N/A

         \[\leadsto \color{blue}{\left(\left(\mathsf{neg}\left(a \cdot \left(t - 1\right)\right)\right) + \left(\mathsf{neg}\left(-1 \cdot z\right)\right)\right)} + \left(x + b \cdot \left(t - 2\right)\right) \]

      5. mul-1-neg N/A

         \[\leadsto \left(\left(\mathsf{neg}\left(a \cdot \left(t - 1\right)\right)\right) + \left(\mathsf{neg}\left(\color{blue}{\left(\mathsf{neg}\left(z\right)\right)}\right)\right)\right) + \left(x + b \cdot \left(t - 2\right)\right) \]

      6. remove-double-neg N/A

         \[\leadsto \left(\left(\mathsf{neg}\left(a \cdot \left(t - 1\right)\right)\right) + \color{blue}{z}\right) + \left(x + b \cdot \left(t - 2\right)\right) \]

      7. associate-+l+ N/A

         \[\leadsto \color{blue}{\left(\mathsf{neg}\left(a \cdot \left(t - 1\right)\right)\right) + \left(z + \left(x + b \cdot \left(t - 2\right)\right)\right)} \]

      8. *-commutative N/A

         \[\leadsto \left(\mathsf{neg}\left(\color{blue}{\left(t - 1\right) \cdot a}\right)\right) + \left(z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      9. distribute-lft-neg-in N/A

         \[\leadsto \color{blue}{\left(\mathsf{neg}\left(\left(t - 1\right)\right)\right) \cdot a} + \left(z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      10. mul-1-neg N/A

          \[\leadsto \color{blue}{\left(-1 \cdot \left(t - 1\right)\right)} \cdot a + \left(z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      11. lower-fma.f64 N/A

          \[\leadsto \color{blue}{\mathsf{fma}\left(-1 \cdot \left(t - 1\right), a, z + \left(x + b \cdot \left(t - 2\right)\right)\right)} \]

      12. sub-neg N/A

          \[\leadsto \mathsf{fma}\left(-1 \cdot \color{blue}{\left(t + \left(\mathsf{neg}\left(1\right)\right)\right)}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      13. metadata-eval N/A

          \[\leadsto \mathsf{fma}\left(-1 \cdot \left(t + \color{blue}{-1}\right), a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      14. distribute-lft-in N/A

          \[\leadsto \mathsf{fma}\left(\color{blue}{-1 \cdot t + -1 \cdot -1}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      15. metadata-eval N/A

          \[\leadsto \mathsf{fma}\left(-1 \cdot t + \color{blue}{1}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      16. +-commutative N/A

          \[\leadsto \mathsf{fma}\left(\color{blue}{1 + -1 \cdot t}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      17. neg-mul-1 N/A

          \[\leadsto \mathsf{fma}\left(1 + \color{blue}{\left(\mathsf{neg}\left(t\right)\right)}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      18. sub-neg N/A

          \[\leadsto \mathsf{fma}\left(\color{blue}{1 - t}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      19. lower--.f64 N/A

          \[\leadsto \mathsf{fma}\left(\color{blue}{1 - t}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      20. lower-+.f64 N/A

          \[\leadsto \mathsf{fma}\left(1 - t, a, \color{blue}{z + \left(x + b \cdot \left(t - 2\right)\right)}\right) \]

   5. Applied rewrites 88.8%

      \[\leadsto \color{blue}{\mathsf{fma}\left(1 - t, a, z + \mathsf{fma}\left(t - 2, b, x\right)\right)} \]

   if -7.8e6 < t < 4.2999999999999998e42

   1. Initial program 97.7%

      \[\left(\left(x - \left(y - 1\right) \cdot z\right) - \left(t - 1\right) \cdot a\right) + \left(\left(y + t\right) - 2\right) \cdot b \]
   2. Add Preprocessing

   3. Taylor expanded in t around 0

      \[\leadsto \color{blue}{\left(x + b \cdot \left(y - 2\right)\right) - \left(-1 \cdot a + z \cdot \left(y - 1\right)\right)} \]
   4. Step-by-step derivation

      1. sub-neg N/A

         \[\leadsto \color{blue}{\left(x + b \cdot \left(y - 2\right)\right) + \left(\mathsf{neg}\left(\left(-1 \cdot a + z \cdot \left(y - 1\right)\right)\right)\right)} \]

      2. +-commutative N/A

         \[\leadsto \color{blue}{\left(\mathsf{neg}\left(\left(-1 \cdot a + z \cdot \left(y - 1\right)\right)\right)\right) + \left(x + b \cdot \left(y - 2\right)\right)} \]

      3. +-commutative N/A

         \[\leadsto \left(\mathsf{neg}\left(\color{blue}{\left(z \cdot \left(y - 1\right) + -1 \cdot a\right)}\right)\right) + \left(x + b \cdot \left(y - 2\right)\right) \]

      4. distribute-neg-in N/A

         \[\leadsto \color{blue}{\left(\left(\mathsf{neg}\left(z \cdot \left(y - 1\right)\right)\right) + \left(\mathsf{neg}\left(-1 \cdot a\right)\right)\right)} + \left(x + b \cdot \left(y - 2\right)\right) \]

      5. mul-1-neg N/A

         \[\leadsto \left(\left(\mathsf{neg}\left(z \cdot \left(y - 1\right)\right)\right) + \left(\mathsf{neg}\left(\color{blue}{\left(\mathsf{neg}\left(a\right)\right)}\right)\right)\right) + \left(x + b \cdot \left(y - 2\right)\right) \]

      6. remove-double-neg N/A

         \[\leadsto \left(\left(\mathsf{neg}\left(z \cdot \left(y - 1\right)\right)\right) + \color{blue}{a}\right) + \left(x + b \cdot \left(y - 2\right)\right) \]

      7. associate-+l+ N/A

         \[\leadsto \color{blue}{\left(\mathsf{neg}\left(z \cdot \left(y - 1\right)\right)\right) + \left(a + \left(x + b \cdot \left(y - 2\right)\right)\right)} \]

      8. *-commutative N/A

         \[\leadsto \left(\mathsf{neg}\left(\color{blue}{\left(y - 1\right) \cdot z}\right)\right) + \left(a + \left(x + b \cdot \left(y - 2\right)\right)\right) \]

      9. distribute-lft-neg-in N/A

         \[\leadsto \color{blue}{\left(\mathsf{neg}\left(\left(y - 1\right)\right)\right) \cdot z} + \left(a + \left(x + b \cdot \left(y - 2\right)\right)\right) \]

      10. mul-1-neg N/A

          \[\leadsto \color{blue}{\left(-1 \cdot \left(y - 1\right)\right)} \cdot z + \left(a + \left(x + b \cdot \left(y - 2\right)\right)\right) \]

      11. lower-fma.f64 N/A

          \[\leadsto \color{blue}{\mathsf{fma}\left(-1 \cdot \left(y - 1\right), z, a + \left(x + b \cdot \left(y - 2\right)\right)\right)} \]

      12. sub-neg N/A

          \[\leadsto \mathsf{fma}\left(-1 \cdot \color{blue}{\left(y + \left(\mathsf{neg}\left(1\right)\right)\right)}, z, a + \left(x + b \cdot \left(y - 2\right)\right)\right) \]

      13. metadata-eval N/A

          \[\leadsto \mathsf{fma}\left(-1 \cdot \left(y + \color{blue}{-1}\right), z, a + \left(x + b \cdot \left(y - 2\right)\right)\right) \]

      14. distribute-lft-in N/A

          \[\leadsto \mathsf{fma}\left(\color{blue}{-1 \cdot y + -1 \cdot -1}, z, a + \left(x + b \cdot \left(y - 2\right)\right)\right) \]

      15. metadata-eval N/A

          \[\leadsto \mathsf{fma}\left(-1 \cdot y + \color{blue}{1}, z, a + \left(x + b \cdot \left(y - 2\right)\right)\right) \]

      16. +-commutative N/A

          \[\leadsto \mathsf{fma}\left(\color{blue}{1 + -1 \cdot y}, z, a + \left(x + b \cdot \left(y - 2\right)\right)\right) \]

      17. neg-mul-1 N/A

          \[\leadsto \mathsf{fma}\left(1 + \color{blue}{\left(\mathsf{neg}\left(y\right)\right)}, z, a + \left(x + b \cdot \left(y - 2\right)\right)\right) \]

      18. sub-neg N/A

          \[\leadsto \mathsf{fma}\left(\color{blue}{1 - y}, z, a + \left(x + b \cdot \left(y - 2\right)\right)\right) \]

      19. lower--.f64 N/A

          \[\leadsto \mathsf{fma}\left(\color{blue}{1 - y}, z, a + \left(x + b \cdot \left(y - 2\right)\right)\right) \]

      20. lower-+.f64 N/A

          \[\leadsto \mathsf{fma}\left(1 - y, z, \color{blue}{a + \left(x + b \cdot \left(y - 2\right)\right)}\right) \]

   5. Applied rewrites 98.4%

      \[\leadsto \color{blue}{\mathsf{fma}\left(1 - y, z, a + \mathsf{fma}\left(y - 2, b, x\right)\right)} \]
3. Recombined 2 regimes into one program.

4. Final simplification 93.7%

   \[\leadsto \begin{array}{l} \mathbf{if}\;t \leq -7800000 \lor \neg \left(t \leq 4.3 \cdot 10^{+42}\right):\\ \;\;\;\;\mathsf{fma}\left(1 - t, a, z + \mathsf{fma}\left(t - 2, b, x\right)\right)\\ \mathbf{else}:\\ \;\;\;\;\mathsf{fma}\left(1 - y, z, a + \mathsf{fma}\left(y - 2, b, x\right)\right)\\ \end{array} \]
5. Add Preprocessing

Alternative 4: 81.9% accurate, 1.1× speedup

Math

\[\begin{array}{l} \\ \begin{array}{l} \mathbf{if}\;b \leq -1.8 \cdot 10^{+85} \lor \neg \left(b \leq 4.4 \cdot 10^{+109}\right):\\ \;\;\;\;\left(\left(t + y\right) - 2\right) \cdot b\\ \mathbf{else}:\\ \;\;\;\;x - \mathsf{fma}\left(a, t - 1, z \cdot \left(y - 1\right)\right)\\ \end{array} \end{array} \]

FPCore

(FPCore (x y z t a b)
 :precision binary64
 (if (or (<= b -1.8e+85) (not (<= b 4.4e+109)))
   (* (- (+ t y) 2.0) b)
   (- x (fma a (- t 1.0) (* z (- y 1.0))))))

C

double code(double x, double y, double z, double t, double a, double b) {
	double tmp;
	if ((b <= -1.8e+85) || !(b <= 4.4e+109)) {
		tmp = ((t + y) - 2.0) * b;
	} else {
		tmp = x - fma(a, (t - 1.0), (z * (y - 1.0)));
	}
	return tmp;
}

Julia

function code(x, y, z, t, a, b)
	tmp = 0.0
	if ((b <= -1.8e+85) || !(b <= 4.4e+109))
		tmp = Float64(Float64(Float64(t + y) - 2.0) * b);
	else
		tmp = Float64(x - fma(a, Float64(t - 1.0), Float64(z * Float64(y - 1.0))));
	end
	return tmp
end

Wolfram

code[x_, y_, z_, t_, a_, b_] := If[Or[LessEqual[b, -1.8e+85], N[Not[LessEqual[b, 4.4e+109]], $MachinePrecision]], N[(N[(N[(t + y), $MachinePrecision] - 2.0), $MachinePrecision] * b), $MachinePrecision], N[(x - N[(a * N[(t - 1.0), $MachinePrecision] + N[(z * N[(y - 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]

Derivation

1. Split input into 2 regimes

2. if b < -1.7999999999999999e85 or 4.3999999999999998e109 < b

   1. Initial program 87.8%

      \[\left(\left(x - \left(y - 1\right) \cdot z\right) - \left(t - 1\right) \cdot a\right) + \left(\left(y + t\right) - 2\right) \cdot b \]
   2. Add Preprocessing

   3. Taylor expanded in y around 0

      \[\leadsto \color{blue}{\left(x + b \cdot \left(t - 2\right)\right) - \left(-1 \cdot z + a \cdot \left(t - 1\right)\right)} \]
   4. Step-by-step derivation

      1. sub-neg N/A

         \[\leadsto \color{blue}{\left(x + b \cdot \left(t - 2\right)\right) + \left(\mathsf{neg}\left(\left(-1 \cdot z + a \cdot \left(t - 1\right)\right)\right)\right)} \]

      2. +-commutative N/A

         \[\leadsto \color{blue}{\left(\mathsf{neg}\left(\left(-1 \cdot z + a \cdot \left(t - 1\right)\right)\right)\right) + \left(x + b \cdot \left(t - 2\right)\right)} \]

      3. +-commutative N/A

         \[\leadsto \left(\mathsf{neg}\left(\color{blue}{\left(a \cdot \left(t - 1\right) + -1 \cdot z\right)}\right)\right) + \left(x + b \cdot \left(t - 2\right)\right) \]

      4. distribute-neg-in N/A

         \[\leadsto \color{blue}{\left(\left(\mathsf{neg}\left(a \cdot \left(t - 1\right)\right)\right) + \left(\mathsf{neg}\left(-1 \cdot z\right)\right)\right)} + \left(x + b \cdot \left(t - 2\right)\right) \]

      5. mul-1-neg N/A

         \[\leadsto \left(\left(\mathsf{neg}\left(a \cdot \left(t - 1\right)\right)\right) + \left(\mathsf{neg}\left(\color{blue}{\left(\mathsf{neg}\left(z\right)\right)}\right)\right)\right) + \left(x + b \cdot \left(t - 2\right)\right) \]

      6. remove-double-neg N/A

         \[\leadsto \left(\left(\mathsf{neg}\left(a \cdot \left(t - 1\right)\right)\right) + \color{blue}{z}\right) + \left(x + b \cdot \left(t - 2\right)\right) \]

      7. associate-+l+ N/A

         \[\leadsto \color{blue}{\left(\mathsf{neg}\left(a \cdot \left(t - 1\right)\right)\right) + \left(z + \left(x + b \cdot \left(t - 2\right)\right)\right)} \]

      8. *-commutative N/A

         \[\leadsto \left(\mathsf{neg}\left(\color{blue}{\left(t - 1\right) \cdot a}\right)\right) + \left(z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      9. distribute-lft-neg-in N/A

         \[\leadsto \color{blue}{\left(\mathsf{neg}\left(\left(t - 1\right)\right)\right) \cdot a} + \left(z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      10. mul-1-neg N/A

          \[\leadsto \color{blue}{\left(-1 \cdot \left(t - 1\right)\right)} \cdot a + \left(z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      11. lower-fma.f64 N/A

          \[\leadsto \color{blue}{\mathsf{fma}\left(-1 \cdot \left(t - 1\right), a, z + \left(x + b \cdot \left(t - 2\right)\right)\right)} \]

      12. sub-neg N/A

          \[\leadsto \mathsf{fma}\left(-1 \cdot \color{blue}{\left(t + \left(\mathsf{neg}\left(1\right)\right)\right)}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      13. metadata-eval N/A

          \[\leadsto \mathsf{fma}\left(-1 \cdot \left(t + \color{blue}{-1}\right), a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      14. distribute-lft-in N/A

          \[\leadsto \mathsf{fma}\left(\color{blue}{-1 \cdot t + -1 \cdot -1}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      15. metadata-eval N/A

          \[\leadsto \mathsf{fma}\left(-1 \cdot t + \color{blue}{1}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      16. +-commutative N/A

          \[\leadsto \mathsf{fma}\left(\color{blue}{1 + -1 \cdot t}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      17. neg-mul-1 N/A

          \[\leadsto \mathsf{fma}\left(1 + \color{blue}{\left(\mathsf{neg}\left(t\right)\right)}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      18. sub-neg N/A

          \[\leadsto \mathsf{fma}\left(\color{blue}{1 - t}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      19. lower--.f64 N/A

          \[\leadsto \mathsf{fma}\left(\color{blue}{1 - t}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      20. lower-+.f64 N/A

          \[\leadsto \mathsf{fma}\left(1 - t, a, \color{blue}{z + \left(x + b \cdot \left(t - 2\right)\right)}\right) \]

   5. Applied rewrites 69.3%

      \[\leadsto \color{blue}{\mathsf{fma}\left(1 - t, a, z + \mathsf{fma}\left(t - 2, b, x\right)\right)} \]

   6. Taylor expanded in b around inf

      \[\leadsto \color{blue}{b \cdot \left(\left(t + y\right) - 2\right)} \]
   7. Step-by-step derivation

      1. *-commutative N/A

         \[\leadsto \color{blue}{\left(\left(t + y\right) - 2\right) \cdot b} \]

      2. lower-*.f64 N/A

         \[\leadsto \color{blue}{\left(\left(t + y\right) - 2\right) \cdot b} \]

      3. lower--.f64 N/A

         \[\leadsto \color{blue}{\left(\left(t + y\right) - 2\right)} \cdot b \]

      4. lower-+.f64 79.1

         \[\leadsto \left(\color{blue}{\left(t + y\right)} - 2\right) \cdot b \]

   8. Applied rewrites 79.1%

      \[\leadsto \color{blue}{\left(\left(t + y\right) - 2\right) \cdot b} \]

   if -1.7999999999999999e85 < b < 4.3999999999999998e109

   1. Initial program 98.8%

      \[\left(\left(x - \left(y - 1\right) \cdot z\right) - \left(t - 1\right) \cdot a\right) + \left(\left(y + t\right) - 2\right) \cdot b \]
   2. Add Preprocessing

   3. Taylor expanded in y around 0

      \[\leadsto \color{blue}{\left(x + b \cdot \left(t - 2\right)\right) - \left(-1 \cdot z + a \cdot \left(t - 1\right)\right)} \]
   4. Step-by-step derivation

      1. sub-neg N/A

         \[\leadsto \color{blue}{\left(x + b \cdot \left(t - 2\right)\right) + \left(\mathsf{neg}\left(\left(-1 \cdot z + a \cdot \left(t - 1\right)\right)\right)\right)} \]

      2. +-commutative N/A

         \[\leadsto \color{blue}{\left(\mathsf{neg}\left(\left(-1 \cdot z + a \cdot \left(t - 1\right)\right)\right)\right) + \left(x + b \cdot \left(t - 2\right)\right)} \]

      3. +-commutative N/A

         \[\leadsto \left(\mathsf{neg}\left(\color{blue}{\left(a \cdot \left(t - 1\right) + -1 \cdot z\right)}\right)\right) + \left(x + b \cdot \left(t - 2\right)\right) \]

      4. distribute-neg-in N/A

         \[\leadsto \color{blue}{\left(\left(\mathsf{neg}\left(a \cdot \left(t - 1\right)\right)\right) + \left(\mathsf{neg}\left(-1 \cdot z\right)\right)\right)} + \left(x + b \cdot \left(t - 2\right)\right) \]

      5. mul-1-neg N/A

         \[\leadsto \left(\left(\mathsf{neg}\left(a \cdot \left(t - 1\right)\right)\right) + \left(\mathsf{neg}\left(\color{blue}{\left(\mathsf{neg}\left(z\right)\right)}\right)\right)\right) + \left(x + b \cdot \left(t - 2\right)\right) \]

      6. remove-double-neg N/A

         \[\leadsto \left(\left(\mathsf{neg}\left(a \cdot \left(t - 1\right)\right)\right) + \color{blue}{z}\right) + \left(x + b \cdot \left(t - 2\right)\right) \]

      7. associate-+l+ N/A

         \[\leadsto \color{blue}{\left(\mathsf{neg}\left(a \cdot \left(t - 1\right)\right)\right) + \left(z + \left(x + b \cdot \left(t - 2\right)\right)\right)} \]

      8. *-commutative N/A

         \[\leadsto \left(\mathsf{neg}\left(\color{blue}{\left(t - 1\right) \cdot a}\right)\right) + \left(z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      9. distribute-lft-neg-in N/A

         \[\leadsto \color{blue}{\left(\mathsf{neg}\left(\left(t - 1\right)\right)\right) \cdot a} + \left(z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      10. mul-1-neg N/A

          \[\leadsto \color{blue}{\left(-1 \cdot \left(t - 1\right)\right)} \cdot a + \left(z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      11. lower-fma.f64 N/A

          \[\leadsto \color{blue}{\mathsf{fma}\left(-1 \cdot \left(t - 1\right), a, z + \left(x + b \cdot \left(t - 2\right)\right)\right)} \]

      12. sub-neg N/A

          \[\leadsto \mathsf{fma}\left(-1 \cdot \color{blue}{\left(t + \left(\mathsf{neg}\left(1\right)\right)\right)}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      13. metadata-eval N/A

          \[\leadsto \mathsf{fma}\left(-1 \cdot \left(t + \color{blue}{-1}\right), a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      14. distribute-lft-in N/A

          \[\leadsto \mathsf{fma}\left(\color{blue}{-1 \cdot t + -1 \cdot -1}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      15. metadata-eval N/A

          \[\leadsto \mathsf{fma}\left(-1 \cdot t + \color{blue}{1}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      16. +-commutative N/A

          \[\leadsto \mathsf{fma}\left(\color{blue}{1 + -1 \cdot t}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      17. neg-mul-1 N/A

          \[\leadsto \mathsf{fma}\left(1 + \color{blue}{\left(\mathsf{neg}\left(t\right)\right)}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      18. sub-neg N/A

          \[\leadsto \mathsf{fma}\left(\color{blue}{1 - t}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      19. lower--.f64 N/A

          \[\leadsto \mathsf{fma}\left(\color{blue}{1 - t}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      20. lower-+.f64 N/A

          \[\leadsto \mathsf{fma}\left(1 - t, a, \color{blue}{z + \left(x + b \cdot \left(t - 2\right)\right)}\right) \]

   5. Applied rewrites 78.5%

      \[\leadsto \color{blue}{\mathsf{fma}\left(1 - t, a, z + \mathsf{fma}\left(t - 2, b, x\right)\right)} \]

   6. Taylor expanded in b around 0

      \[\leadsto \color{blue}{x - \left(a \cdot \left(t - 1\right) + z \cdot \left(y - 1\right)\right)} \]
   7. Step-by-step derivation

      1. lower--.f64 N/A

         \[\leadsto \color{blue}{x - \left(a \cdot \left(t - 1\right) + z \cdot \left(y - 1\right)\right)} \]

      2. lower-fma.f64 N/A

         \[\leadsto x - \color{blue}{\mathsf{fma}\left(a, t - 1, z \cdot \left(y - 1\right)\right)} \]

      3. lower--.f64 N/A

         \[\leadsto x - \mathsf{fma}\left(a, \color{blue}{t - 1}, z \cdot \left(y - 1\right)\right) \]

      4. lower-*.f64 N/A

         \[\leadsto x - \mathsf{fma}\left(a, t - 1, \color{blue}{z \cdot \left(y - 1\right)}\right) \]

      5. lower--.f64 88.7

         \[\leadsto x - \mathsf{fma}\left(a, t - 1, z \cdot \color{blue}{\left(y - 1\right)}\right) \]

   8. Applied rewrites 88.7%

      \[\leadsto \color{blue}{x - \mathsf{fma}\left(a, t - 1, z \cdot \left(y - 1\right)\right)} \]
3. Recombined 2 regimes into one program.

4. Final simplification 85.6%

   \[\leadsto \begin{array}{l} \mathbf{if}\;b \leq -1.8 \cdot 10^{+85} \lor \neg \left(b \leq 4.4 \cdot 10^{+109}\right):\\ \;\;\;\;\left(\left(t + y\right) - 2\right) \cdot b\\ \mathbf{else}:\\ \;\;\;\;x - \mathsf{fma}\left(a, t - 1, z \cdot \left(y - 1\right)\right)\\ \end{array} \]
5. Add Preprocessing

Alternative 5: 74.3% accurate, 1.2× speedup

Math

\[\begin{array}{l} \\ \begin{array}{l} \mathbf{if}\;y \leq -1.15 \cdot 10^{+146} \lor \neg \left(y \leq 2.95 \cdot 10^{+99}\right):\\ \;\;\;\;\left(b - z\right) \cdot y\\ \mathbf{else}:\\ \;\;\;\;\mathsf{fma}\left(1 - t, a, \mathsf{fma}\left(t - 2, b, x\right)\right)\\ \end{array} \end{array} \]

FPCore

(FPCore (x y z t a b)
 :precision binary64
 (if (or (<= y -1.15e+146) (not (<= y 2.95e+99)))
   (* (- b z) y)
   (fma (- 1.0 t) a (fma (- t 2.0) b x))))

C

double code(double x, double y, double z, double t, double a, double b) {
	double tmp;
	if ((y <= -1.15e+146) || !(y <= 2.95e+99)) {
		tmp = (b - z) * y;
	} else {
		tmp = fma((1.0 - t), a, fma((t - 2.0), b, x));
	}
	return tmp;
}

Julia

function code(x, y, z, t, a, b)
	tmp = 0.0
	if ((y <= -1.15e+146) || !(y <= 2.95e+99))
		tmp = Float64(Float64(b - z) * y);
	else
		tmp = fma(Float64(1.0 - t), a, fma(Float64(t - 2.0), b, x));
	end
	return tmp
end

Wolfram

code[x_, y_, z_, t_, a_, b_] := If[Or[LessEqual[y, -1.15e+146], N[Not[LessEqual[y, 2.95e+99]], $MachinePrecision]], N[(N[(b - z), $MachinePrecision] * y), $MachinePrecision], N[(N[(1.0 - t), $MachinePrecision] * a + N[(N[(t - 2.0), $MachinePrecision] * b + x), $MachinePrecision]), $MachinePrecision]]

Derivation

1. Split input into 2 regimes

2. if y < -1.15e146 or 2.9499999999999999e99 < y

   1. Initial program 91.5%

      \[\left(\left(x - \left(y - 1\right) \cdot z\right) - \left(t - 1\right) \cdot a\right) + \left(\left(y + t\right) - 2\right) \cdot b \]
   2. Add Preprocessing

   3. Taylor expanded in y around inf

      \[\leadsto \color{blue}{y \cdot \left(b - z\right)} \]
   4. Step-by-step derivation

      1. *-commutative N/A

         \[\leadsto \color{blue}{\left(b - z\right) \cdot y} \]

      2. lower-*.f64 N/A

         \[\leadsto \color{blue}{\left(b - z\right) \cdot y} \]

      3. lower--.f64 76.5

         \[\leadsto \color{blue}{\left(b - z\right)} \cdot y \]

   5. Applied rewrites 76.5%

      \[\leadsto \color{blue}{\left(b - z\right) \cdot y} \]

   if -1.15e146 < y < 2.9499999999999999e99

   1. Initial program 96.7%

      \[\left(\left(x - \left(y - 1\right) \cdot z\right) - \left(t - 1\right) \cdot a\right) + \left(\left(y + t\right) - 2\right) \cdot b \]
   2. Add Preprocessing

   3. Taylor expanded in y around 0

      \[\leadsto \color{blue}{\left(x + b \cdot \left(t - 2\right)\right) - \left(-1 \cdot z + a \cdot \left(t - 1\right)\right)} \]
   4. Step-by-step derivation

      1. sub-neg N/A

         \[\leadsto \color{blue}{\left(x + b \cdot \left(t - 2\right)\right) + \left(\mathsf{neg}\left(\left(-1 \cdot z + a \cdot \left(t - 1\right)\right)\right)\right)} \]

      2. +-commutative N/A

         \[\leadsto \color{blue}{\left(\mathsf{neg}\left(\left(-1 \cdot z + a \cdot \left(t - 1\right)\right)\right)\right) + \left(x + b \cdot \left(t - 2\right)\right)} \]

      3. +-commutative N/A

         \[\leadsto \left(\mathsf{neg}\left(\color{blue}{\left(a \cdot \left(t - 1\right) + -1 \cdot z\right)}\right)\right) + \left(x + b \cdot \left(t - 2\right)\right) \]

      4. distribute-neg-in N/A

         \[\leadsto \color{blue}{\left(\left(\mathsf{neg}\left(a \cdot \left(t - 1\right)\right)\right) + \left(\mathsf{neg}\left(-1 \cdot z\right)\right)\right)} + \left(x + b \cdot \left(t - 2\right)\right) \]

      5. mul-1-neg N/A

         \[\leadsto \left(\left(\mathsf{neg}\left(a \cdot \left(t - 1\right)\right)\right) + \left(\mathsf{neg}\left(\color{blue}{\left(\mathsf{neg}\left(z\right)\right)}\right)\right)\right) + \left(x + b \cdot \left(t - 2\right)\right) \]

      6. remove-double-neg N/A

         \[\leadsto \left(\left(\mathsf{neg}\left(a \cdot \left(t - 1\right)\right)\right) + \color{blue}{z}\right) + \left(x + b \cdot \left(t - 2\right)\right) \]

      7. associate-+l+ N/A

         \[\leadsto \color{blue}{\left(\mathsf{neg}\left(a \cdot \left(t - 1\right)\right)\right) + \left(z + \left(x + b \cdot \left(t - 2\right)\right)\right)} \]

      8. *-commutative N/A

         \[\leadsto \left(\mathsf{neg}\left(\color{blue}{\left(t - 1\right) \cdot a}\right)\right) + \left(z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      9. distribute-lft-neg-in N/A

         \[\leadsto \color{blue}{\left(\mathsf{neg}\left(\left(t - 1\right)\right)\right) \cdot a} + \left(z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      10. mul-1-neg N/A

          \[\leadsto \color{blue}{\left(-1 \cdot \left(t - 1\right)\right)} \cdot a + \left(z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      11. lower-fma.f64 N/A

          \[\leadsto \color{blue}{\mathsf{fma}\left(-1 \cdot \left(t - 1\right), a, z + \left(x + b \cdot \left(t - 2\right)\right)\right)} \]

      12. sub-neg N/A

          \[\leadsto \mathsf{fma}\left(-1 \cdot \color{blue}{\left(t + \left(\mathsf{neg}\left(1\right)\right)\right)}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      13. metadata-eval N/A

          \[\leadsto \mathsf{fma}\left(-1 \cdot \left(t + \color{blue}{-1}\right), a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      14. distribute-lft-in N/A

          \[\leadsto \mathsf{fma}\left(\color{blue}{-1 \cdot t + -1 \cdot -1}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      15. metadata-eval N/A

          \[\leadsto \mathsf{fma}\left(-1 \cdot t + \color{blue}{1}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      16. +-commutative N/A

          \[\leadsto \mathsf{fma}\left(\color{blue}{1 + -1 \cdot t}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      17. neg-mul-1 N/A

          \[\leadsto \mathsf{fma}\left(1 + \color{blue}{\left(\mathsf{neg}\left(t\right)\right)}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      18. sub-neg N/A

          \[\leadsto \mathsf{fma}\left(\color{blue}{1 - t}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      19. lower--.f64 N/A

          \[\leadsto \mathsf{fma}\left(\color{blue}{1 - t}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      20. lower-+.f64 N/A

          \[\leadsto \mathsf{fma}\left(1 - t, a, \color{blue}{z + \left(x + b \cdot \left(t - 2\right)\right)}\right) \]

   5. Applied rewrites 92.8%

      \[\leadsto \color{blue}{\mathsf{fma}\left(1 - t, a, z + \mathsf{fma}\left(t - 2, b, x\right)\right)} \]

   6. Taylor expanded in z around 0

      \[\leadsto x + \color{blue}{\left(a \cdot \left(1 - t\right) + b \cdot \left(t - 2\right)\right)} \]
   7. Step-by-step derivation

   8. Applied rewrites 78.5%

      \[\leadsto \mathsf{fma}\left(1 - t, \color{blue}{a}, \mathsf{fma}\left(t - 2, b, x\right)\right) \]
3. Recombined 2 regimes into one program.

4. Final simplification 77.9%

   \[\leadsto \begin{array}{l} \mathbf{if}\;y \leq -1.15 \cdot 10^{+146} \lor \neg \left(y \leq 2.95 \cdot 10^{+99}\right):\\ \;\;\;\;\left(b - z\right) \cdot y\\ \mathbf{else}:\\ \;\;\;\;\mathsf{fma}\left(1 - t, a, \mathsf{fma}\left(t - 2, b, x\right)\right)\\ \end{array} \]
5. Add Preprocessing

Alternative 6: 65.9% accurate, 1.2× speedup

Math

\[\begin{array}{l} \\ \begin{array}{l} t_1 := \left(b - a\right) \cdot t\\ \mathbf{if}\;t \leq -2.9 \cdot 10^{+203}:\\ \;\;\;\;t\_1\\ \mathbf{elif}\;t \leq -145:\\ \;\;\;\;\mathsf{fma}\left(-t, a, z + x\right)\\ \mathbf{elif}\;t \leq 1.75 \cdot 10^{+47}:\\ \;\;\;\;\mathsf{fma}\left(1 - y, z, a\right) + x\\ \mathbf{else}:\\ \;\;\;\;t\_1\\ \end{array} \end{array} \]

FPCore

(FPCore (x y z t a b)
 :precision binary64
 (let* ((t_1 (* (- b a) t)))
   (if (<= t -2.9e+203)
     t_1
     (if (<= t -145.0)
       (fma (- t) a (+ z x))
       (if (<= t 1.75e+47) (+ (fma (- 1.0 y) z a) x) t_1)))))

C

double code(double x, double y, double z, double t, double a, double b) {
	double t_1 = (b - a) * t;
	double tmp;
	if (t <= -2.9e+203) {
		tmp = t_1;
	} else if (t <= -145.0) {
		tmp = fma(-t, a, (z + x));
	} else if (t <= 1.75e+47) {
		tmp = fma((1.0 - y), z, a) + x;
	} else {
		tmp = t_1;
	}
	return tmp;
}

Julia

function code(x, y, z, t, a, b)
	t_1 = Float64(Float64(b - a) * t)
	tmp = 0.0
	if (t <= -2.9e+203)
		tmp = t_1;
	elseif (t <= -145.0)
		tmp = fma(Float64(-t), a, Float64(z + x));
	elseif (t <= 1.75e+47)
		tmp = Float64(fma(Float64(1.0 - y), z, a) + x);
	else
		tmp = t_1;
	end
	return tmp
end

Wolfram

code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(N[(b - a), $MachinePrecision] * t), $MachinePrecision]}, If[LessEqual[t, -2.9e+203], t$95$1, If[LessEqual[t, -145.0], N[((-t) * a + N[(z + x), $MachinePrecision]), $MachinePrecision], If[LessEqual[t, 1.75e+47], N[(N[(N[(1.0 - y), $MachinePrecision] * z + a), $MachinePrecision] + x), $MachinePrecision], t$95$1]]]]

Derivation

1. Split input into 3 regimes

2. if t < -2.90000000000000011e203 or 1.75000000000000008e47 < t

   1. Initial program 92.0%

      \[\left(\left(x - \left(y - 1\right) \cdot z\right) - \left(t - 1\right) \cdot a\right) + \left(\left(y + t\right) - 2\right) \cdot b \]
   2. Add Preprocessing

   3. Taylor expanded in t around inf

      \[\leadsto \color{blue}{t \cdot \left(b - a\right)} \]
   4. Step-by-step derivation

      1. *-commutative N/A

         \[\leadsto \color{blue}{\left(b - a\right) \cdot t} \]

      2. lower-*.f64 N/A

         \[\leadsto \color{blue}{\left(b - a\right) \cdot t} \]

      3. lower--.f64 84.8

         \[\leadsto \color{blue}{\left(b - a\right)} \cdot t \]

   5. Applied rewrites 84.8%

      \[\leadsto \color{blue}{\left(b - a\right) \cdot t} \]

   if -2.90000000000000011e203 < t < -145

   1. Initial program 92.6%

      \[\left(\left(x - \left(y - 1\right) \cdot z\right) - \left(t - 1\right) \cdot a\right) + \left(\left(y + t\right) - 2\right) \cdot b \]
   2. Add Preprocessing

   3. Taylor expanded in y around 0

      \[\leadsto \color{blue}{\left(x + b \cdot \left(t - 2\right)\right) - \left(-1 \cdot z + a \cdot \left(t - 1\right)\right)} \]
   4. Step-by-step derivation

      1. sub-neg N/A

         \[\leadsto \color{blue}{\left(x + b \cdot \left(t - 2\right)\right) + \left(\mathsf{neg}\left(\left(-1 \cdot z + a \cdot \left(t - 1\right)\right)\right)\right)} \]

      2. +-commutative N/A

         \[\leadsto \color{blue}{\left(\mathsf{neg}\left(\left(-1 \cdot z + a \cdot \left(t - 1\right)\right)\right)\right) + \left(x + b \cdot \left(t - 2\right)\right)} \]

      3. +-commutative N/A

         \[\leadsto \left(\mathsf{neg}\left(\color{blue}{\left(a \cdot \left(t - 1\right) + -1 \cdot z\right)}\right)\right) + \left(x + b \cdot \left(t - 2\right)\right) \]

      4. distribute-neg-in N/A

         \[\leadsto \color{blue}{\left(\left(\mathsf{neg}\left(a \cdot \left(t - 1\right)\right)\right) + \left(\mathsf{neg}\left(-1 \cdot z\right)\right)\right)} + \left(x + b \cdot \left(t - 2\right)\right) \]

      5. mul-1-neg N/A

         \[\leadsto \left(\left(\mathsf{neg}\left(a \cdot \left(t - 1\right)\right)\right) + \left(\mathsf{neg}\left(\color{blue}{\left(\mathsf{neg}\left(z\right)\right)}\right)\right)\right) + \left(x + b \cdot \left(t - 2\right)\right) \]

      6. remove-double-neg N/A

         \[\leadsto \left(\left(\mathsf{neg}\left(a \cdot \left(t - 1\right)\right)\right) + \color{blue}{z}\right) + \left(x + b \cdot \left(t - 2\right)\right) \]

      7. associate-+l+ N/A

         \[\leadsto \color{blue}{\left(\mathsf{neg}\left(a \cdot \left(t - 1\right)\right)\right) + \left(z + \left(x + b \cdot \left(t - 2\right)\right)\right)} \]

      8. *-commutative N/A

         \[\leadsto \left(\mathsf{neg}\left(\color{blue}{\left(t - 1\right) \cdot a}\right)\right) + \left(z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      9. distribute-lft-neg-in N/A

         \[\leadsto \color{blue}{\left(\mathsf{neg}\left(\left(t - 1\right)\right)\right) \cdot a} + \left(z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      10. mul-1-neg N/A

          \[\leadsto \color{blue}{\left(-1 \cdot \left(t - 1\right)\right)} \cdot a + \left(z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      11. lower-fma.f64 N/A

          \[\leadsto \color{blue}{\mathsf{fma}\left(-1 \cdot \left(t - 1\right), a, z + \left(x + b \cdot \left(t - 2\right)\right)\right)} \]

      12. sub-neg N/A

          \[\leadsto \mathsf{fma}\left(-1 \cdot \color{blue}{\left(t + \left(\mathsf{neg}\left(1\right)\right)\right)}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      13. metadata-eval N/A

          \[\leadsto \mathsf{fma}\left(-1 \cdot \left(t + \color{blue}{-1}\right), a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      14. distribute-lft-in N/A

          \[\leadsto \mathsf{fma}\left(\color{blue}{-1 \cdot t + -1 \cdot -1}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      15. metadata-eval N/A

          \[\leadsto \mathsf{fma}\left(-1 \cdot t + \color{blue}{1}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      16. +-commutative N/A

          \[\leadsto \mathsf{fma}\left(\color{blue}{1 + -1 \cdot t}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      17. neg-mul-1 N/A

          \[\leadsto \mathsf{fma}\left(1 + \color{blue}{\left(\mathsf{neg}\left(t\right)\right)}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      18. sub-neg N/A

          \[\leadsto \mathsf{fma}\left(\color{blue}{1 - t}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      19. lower--.f64 N/A

          \[\leadsto \mathsf{fma}\left(\color{blue}{1 - t}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      20. lower-+.f64 N/A

          \[\leadsto \mathsf{fma}\left(1 - t, a, \color{blue}{z + \left(x + b \cdot \left(t - 2\right)\right)}\right) \]

   5. Applied rewrites 80.7%

      \[\leadsto \color{blue}{\mathsf{fma}\left(1 - t, a, z + \mathsf{fma}\left(t - 2, b, x\right)\right)} \]

   6. Taylor expanded in b around 0

      \[\leadsto \mathsf{fma}\left(1 - t, a, x + z\right) \]
   7. Step-by-step derivation

   8. Applied rewrites 71.6%

      \[\leadsto \mathsf{fma}\left(1 - t, a, z + x\right) \]

   9. Taylor expanded in t around inf

      \[\leadsto \mathsf{fma}\left(-1 \cdot t, a, z + x\right) \]
   10. Step-by-step derivation

   11. Applied rewrites 68.9%

       \[\leadsto \mathsf{fma}\left(-t, a, z + x\right) \]

   if -145 < t < 1.75000000000000008e47

   1. Initial program 98.4%

      \[\left(\left(x - \left(y - 1\right) \cdot z\right) - \left(t - 1\right) \cdot a\right) + \left(\left(y + t\right) - 2\right) \cdot b \]
   2. Add Preprocessing

   3. Taylor expanded in y around 0

      \[\leadsto \color{blue}{\left(x + b \cdot \left(t - 2\right)\right) - \left(-1 \cdot z + a \cdot \left(t - 1\right)\right)} \]
   4. Step-by-step derivation

      1. sub-neg N/A

         \[\leadsto \color{blue}{\left(x + b \cdot \left(t - 2\right)\right) + \left(\mathsf{neg}\left(\left(-1 \cdot z + a \cdot \left(t - 1\right)\right)\right)\right)} \]

      2. +-commutative N/A

         \[\leadsto \color{blue}{\left(\mathsf{neg}\left(\left(-1 \cdot z + a \cdot \left(t - 1\right)\right)\right)\right) + \left(x + b \cdot \left(t - 2\right)\right)} \]

      3. +-commutative N/A

         \[\leadsto \left(\mathsf{neg}\left(\color{blue}{\left(a \cdot \left(t - 1\right) + -1 \cdot z\right)}\right)\right) + \left(x + b \cdot \left(t - 2\right)\right) \]

      4. distribute-neg-in N/A

         \[\leadsto \color{blue}{\left(\left(\mathsf{neg}\left(a \cdot \left(t - 1\right)\right)\right) + \left(\mathsf{neg}\left(-1 \cdot z\right)\right)\right)} + \left(x + b \cdot \left(t - 2\right)\right) \]

      5. mul-1-neg N/A

         \[\leadsto \left(\left(\mathsf{neg}\left(a \cdot \left(t - 1\right)\right)\right) + \left(\mathsf{neg}\left(\color{blue}{\left(\mathsf{neg}\left(z\right)\right)}\right)\right)\right) + \left(x + b \cdot \left(t - 2\right)\right) \]

      6. remove-double-neg N/A

         \[\leadsto \left(\left(\mathsf{neg}\left(a \cdot \left(t - 1\right)\right)\right) + \color{blue}{z}\right) + \left(x + b \cdot \left(t - 2\right)\right) \]

      7. associate-+l+ N/A

         \[\leadsto \color{blue}{\left(\mathsf{neg}\left(a \cdot \left(t - 1\right)\right)\right) + \left(z + \left(x + b \cdot \left(t - 2\right)\right)\right)} \]

      8. *-commutative N/A

         \[\leadsto \left(\mathsf{neg}\left(\color{blue}{\left(t - 1\right) \cdot a}\right)\right) + \left(z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      9. distribute-lft-neg-in N/A

         \[\leadsto \color{blue}{\left(\mathsf{neg}\left(\left(t - 1\right)\right)\right) \cdot a} + \left(z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      10. mul-1-neg N/A

          \[\leadsto \color{blue}{\left(-1 \cdot \left(t - 1\right)\right)} \cdot a + \left(z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      11. lower-fma.f64 N/A

          \[\leadsto \color{blue}{\mathsf{fma}\left(-1 \cdot \left(t - 1\right), a, z + \left(x + b \cdot \left(t - 2\right)\right)\right)} \]

      12. sub-neg N/A

          \[\leadsto \mathsf{fma}\left(-1 \cdot \color{blue}{\left(t + \left(\mathsf{neg}\left(1\right)\right)\right)}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      13. metadata-eval N/A

          \[\leadsto \mathsf{fma}\left(-1 \cdot \left(t + \color{blue}{-1}\right), a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      14. distribute-lft-in N/A

          \[\leadsto \mathsf{fma}\left(\color{blue}{-1 \cdot t + -1 \cdot -1}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      15. metadata-eval N/A

          \[\leadsto \mathsf{fma}\left(-1 \cdot t + \color{blue}{1}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      16. +-commutative N/A

          \[\leadsto \mathsf{fma}\left(\color{blue}{1 + -1 \cdot t}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      17. neg-mul-1 N/A

          \[\leadsto \mathsf{fma}\left(1 + \color{blue}{\left(\mathsf{neg}\left(t\right)\right)}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      18. sub-neg N/A

          \[\leadsto \mathsf{fma}\left(\color{blue}{1 - t}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      19. lower--.f64 N/A

          \[\leadsto \mathsf{fma}\left(\color{blue}{1 - t}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      20. lower-+.f64 N/A

          \[\leadsto \mathsf{fma}\left(1 - t, a, \color{blue}{z + \left(x + b \cdot \left(t - 2\right)\right)}\right) \]

   5. Applied rewrites 62.6%

      \[\leadsto \color{blue}{\mathsf{fma}\left(1 - t, a, z + \mathsf{fma}\left(t - 2, b, x\right)\right)} \]

   6. Taylor expanded in b around 0

      \[\leadsto \color{blue}{x - \left(a \cdot \left(t - 1\right) + z \cdot \left(y - 1\right)\right)} \]
   7. Step-by-step derivation

      1. lower--.f64 N/A

         \[\leadsto \color{blue}{x - \left(a \cdot \left(t - 1\right) + z \cdot \left(y - 1\right)\right)} \]

      2. lower-fma.f64 N/A

         \[\leadsto x - \color{blue}{\mathsf{fma}\left(a, t - 1, z \cdot \left(y - 1\right)\right)} \]

      3. lower--.f64 N/A

         \[\leadsto x - \mathsf{fma}\left(a, \color{blue}{t - 1}, z \cdot \left(y - 1\right)\right) \]

      4. lower-*.f64 N/A

         \[\leadsto x - \mathsf{fma}\left(a, t - 1, \color{blue}{z \cdot \left(y - 1\right)}\right) \]

      5. lower--.f64 71.0

         \[\leadsto x - \mathsf{fma}\left(a, t - 1, z \cdot \color{blue}{\left(y - 1\right)}\right) \]

   8. Applied rewrites 71.0%

      \[\leadsto \color{blue}{x - \mathsf{fma}\left(a, t - 1, z \cdot \left(y - 1\right)\right)} \]

   9. Taylor expanded in t around 0

      \[\leadsto x - \color{blue}{\left(-1 \cdot a + z \cdot \left(y - 1\right)\right)} \]
   10. Step-by-step derivation

   11. Applied rewrites 70.9%

       \[\leadsto \mathsf{fma}\left(1 - y, z, a\right) + \color{blue}{x} \]
3. Recombined 3 regimes into one program.
4. Add Preprocessing

Alternative 7: 50.3% accurate, 1.4× speedup

Math

\[\begin{array}{l} \\ \begin{array}{l} t_1 := \left(b - a\right) \cdot t\\ \mathbf{if}\;t \leq -145:\\ \;\;\;\;t\_1\\ \mathbf{elif}\;t \leq -5.2 \cdot 10^{-215}:\\ \;\;\;\;x - y \cdot z\\ \mathbf{elif}\;t \leq 4.3 \cdot 10^{+42}:\\ \;\;\;\;\left(b - z\right) \cdot y\\ \mathbf{else}:\\ \;\;\;\;t\_1\\ \end{array} \end{array} \]

FPCore

(FPCore (x y z t a b)
 :precision binary64
 (let* ((t_1 (* (- b a) t)))
   (if (<= t -145.0)
     t_1
     (if (<= t -5.2e-215)
       (- x (* y z))
       (if (<= t 4.3e+42) (* (- b z) y) t_1)))))

C

double code(double x, double y, double z, double t, double a, double b) {
	double t_1 = (b - a) * t;
	double tmp;
	if (t <= -145.0) {
		tmp = t_1;
	} else if (t <= -5.2e-215) {
		tmp = x - (y * z);
	} else if (t <= 4.3e+42) {
		tmp = (b - z) * y;
	} else {
		tmp = t_1;
	}
	return tmp;
}

Fortran

real(8) function code(x, y, z, t, a, b)
    real(8), intent (in) :: x
    real(8), intent (in) :: y
    real(8), intent (in) :: z
    real(8), intent (in) :: t
    real(8), intent (in) :: a
    real(8), intent (in) :: b
    real(8) :: t_1
    real(8) :: tmp
    t_1 = (b - a) * t
    if (t <= (-145.0d0)) then
        tmp = t_1
    else if (t <= (-5.2d-215)) then
        tmp = x - (y * z)
    else if (t <= 4.3d+42) then
        tmp = (b - z) * y
    else
        tmp = t_1
    end if
    code = tmp
end function

Java

public static double code(double x, double y, double z, double t, double a, double b) {
	double t_1 = (b - a) * t;
	double tmp;
	if (t <= -145.0) {
		tmp = t_1;
	} else if (t <= -5.2e-215) {
		tmp = x - (y * z);
	} else if (t <= 4.3e+42) {
		tmp = (b - z) * y;
	} else {
		tmp = t_1;
	}
	return tmp;
}

Python

def code(x, y, z, t, a, b):
	t_1 = (b - a) * t
	tmp = 0
	if t <= -145.0:
		tmp = t_1
	elif t <= -5.2e-215:
		tmp = x - (y * z)
	elif t <= 4.3e+42:
		tmp = (b - z) * y
	else:
		tmp = t_1
	return tmp

Julia

function code(x, y, z, t, a, b)
	t_1 = Float64(Float64(b - a) * t)
	tmp = 0.0
	if (t <= -145.0)
		tmp = t_1;
	elseif (t <= -5.2e-215)
		tmp = Float64(x - Float64(y * z));
	elseif (t <= 4.3e+42)
		tmp = Float64(Float64(b - z) * y);
	else
		tmp = t_1;
	end
	return tmp
end

MATLAB

function tmp_2 = code(x, y, z, t, a, b)
	t_1 = (b - a) * t;
	tmp = 0.0;
	if (t <= -145.0)
		tmp = t_1;
	elseif (t <= -5.2e-215)
		tmp = x - (y * z);
	elseif (t <= 4.3e+42)
		tmp = (b - z) * y;
	else
		tmp = t_1;
	end
	tmp_2 = tmp;
end

Wolfram

code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(N[(b - a), $MachinePrecision] * t), $MachinePrecision]}, If[LessEqual[t, -145.0], t$95$1, If[LessEqual[t, -5.2e-215], N[(x - N[(y * z), $MachinePrecision]), $MachinePrecision], If[LessEqual[t, 4.3e+42], N[(N[(b - z), $MachinePrecision] * y), $MachinePrecision], t$95$1]]]]

Derivation

1. Split input into 3 regimes

2. if t < -145 or 4.2999999999999998e42 < t

   1. Initial program 92.2%

      \[\left(\left(x - \left(y - 1\right) \cdot z\right) - \left(t - 1\right) \cdot a\right) + \left(\left(y + t\right) - 2\right) \cdot b \]
   2. Add Preprocessing

   3. Taylor expanded in t around inf

      \[\leadsto \color{blue}{t \cdot \left(b - a\right)} \]
   4. Step-by-step derivation

      1. *-commutative N/A

         \[\leadsto \color{blue}{\left(b - a\right) \cdot t} \]

      2. lower-*.f64 N/A

         \[\leadsto \color{blue}{\left(b - a\right) \cdot t} \]

      3. lower--.f64 72.6

         \[\leadsto \color{blue}{\left(b - a\right)} \cdot t \]

   5. Applied rewrites 72.6%

      \[\leadsto \color{blue}{\left(b - a\right) \cdot t} \]

   if -145 < t < -5.2e-215

   1. Initial program 98.1%

      \[\left(\left(x - \left(y - 1\right) \cdot z\right) - \left(t - 1\right) \cdot a\right) + \left(\left(y + t\right) - 2\right) \cdot b \]
   2. Add Preprocessing

   3. Taylor expanded in y around 0

      \[\leadsto \color{blue}{\left(x + b \cdot \left(t - 2\right)\right) - \left(-1 \cdot z + a \cdot \left(t - 1\right)\right)} \]
   4. Step-by-step derivation

      1. sub-neg N/A

         \[\leadsto \color{blue}{\left(x + b \cdot \left(t - 2\right)\right) + \left(\mathsf{neg}\left(\left(-1 \cdot z + a \cdot \left(t - 1\right)\right)\right)\right)} \]

      2. +-commutative N/A

         \[\leadsto \color{blue}{\left(\mathsf{neg}\left(\left(-1 \cdot z + a \cdot \left(t - 1\right)\right)\right)\right) + \left(x + b \cdot \left(t - 2\right)\right)} \]

      3. +-commutative N/A

         \[\leadsto \left(\mathsf{neg}\left(\color{blue}{\left(a \cdot \left(t - 1\right) + -1 \cdot z\right)}\right)\right) + \left(x + b \cdot \left(t - 2\right)\right) \]

      4. distribute-neg-in N/A

         \[\leadsto \color{blue}{\left(\left(\mathsf{neg}\left(a \cdot \left(t - 1\right)\right)\right) + \left(\mathsf{neg}\left(-1 \cdot z\right)\right)\right)} + \left(x + b \cdot \left(t - 2\right)\right) \]

      5. mul-1-neg N/A

         \[\leadsto \left(\left(\mathsf{neg}\left(a \cdot \left(t - 1\right)\right)\right) + \left(\mathsf{neg}\left(\color{blue}{\left(\mathsf{neg}\left(z\right)\right)}\right)\right)\right) + \left(x + b \cdot \left(t - 2\right)\right) \]

      6. remove-double-neg N/A

         \[\leadsto \left(\left(\mathsf{neg}\left(a \cdot \left(t - 1\right)\right)\right) + \color{blue}{z}\right) + \left(x + b \cdot \left(t - 2\right)\right) \]

      7. associate-+l+ N/A

         \[\leadsto \color{blue}{\left(\mathsf{neg}\left(a \cdot \left(t - 1\right)\right)\right) + \left(z + \left(x + b \cdot \left(t - 2\right)\right)\right)} \]

      8. *-commutative N/A

         \[\leadsto \left(\mathsf{neg}\left(\color{blue}{\left(t - 1\right) \cdot a}\right)\right) + \left(z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      9. distribute-lft-neg-in N/A

         \[\leadsto \color{blue}{\left(\mathsf{neg}\left(\left(t - 1\right)\right)\right) \cdot a} + \left(z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      10. mul-1-neg N/A

          \[\leadsto \color{blue}{\left(-1 \cdot \left(t - 1\right)\right)} \cdot a + \left(z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      11. lower-fma.f64 N/A

          \[\leadsto \color{blue}{\mathsf{fma}\left(-1 \cdot \left(t - 1\right), a, z + \left(x + b \cdot \left(t - 2\right)\right)\right)} \]

      12. sub-neg N/A

          \[\leadsto \mathsf{fma}\left(-1 \cdot \color{blue}{\left(t + \left(\mathsf{neg}\left(1\right)\right)\right)}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      13. metadata-eval N/A

          \[\leadsto \mathsf{fma}\left(-1 \cdot \left(t + \color{blue}{-1}\right), a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      14. distribute-lft-in N/A

          \[\leadsto \mathsf{fma}\left(\color{blue}{-1 \cdot t + -1 \cdot -1}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      15. metadata-eval N/A

          \[\leadsto \mathsf{fma}\left(-1 \cdot t + \color{blue}{1}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      16. +-commutative N/A

          \[\leadsto \mathsf{fma}\left(\color{blue}{1 + -1 \cdot t}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      17. neg-mul-1 N/A

          \[\leadsto \mathsf{fma}\left(1 + \color{blue}{\left(\mathsf{neg}\left(t\right)\right)}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      18. sub-neg N/A

          \[\leadsto \mathsf{fma}\left(\color{blue}{1 - t}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      19. lower--.f64 N/A

          \[\leadsto \mathsf{fma}\left(\color{blue}{1 - t}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      20. lower-+.f64 N/A

          \[\leadsto \mathsf{fma}\left(1 - t, a, \color{blue}{z + \left(x + b \cdot \left(t - 2\right)\right)}\right) \]

   5. Applied rewrites 65.1%

      \[\leadsto \color{blue}{\mathsf{fma}\left(1 - t, a, z + \mathsf{fma}\left(t - 2, b, x\right)\right)} \]

   6. Taylor expanded in b around 0

      \[\leadsto \color{blue}{x - \left(a \cdot \left(t - 1\right) + z \cdot \left(y - 1\right)\right)} \]
   7. Step-by-step derivation

      1. lower--.f64 N/A

         \[\leadsto \color{blue}{x - \left(a \cdot \left(t - 1\right) + z \cdot \left(y - 1\right)\right)} \]

      2. lower-fma.f64 N/A

         \[\leadsto x - \color{blue}{\mathsf{fma}\left(a, t - 1, z \cdot \left(y - 1\right)\right)} \]

      3. lower--.f64 N/A

         \[\leadsto x - \mathsf{fma}\left(a, \color{blue}{t - 1}, z \cdot \left(y - 1\right)\right) \]

      4. lower-*.f64 N/A

         \[\leadsto x - \mathsf{fma}\left(a, t - 1, \color{blue}{z \cdot \left(y - 1\right)}\right) \]

      5. lower--.f64 77.7

         \[\leadsto x - \mathsf{fma}\left(a, t - 1, z \cdot \color{blue}{\left(y - 1\right)}\right) \]

   8. Applied rewrites 77.7%

      \[\leadsto \color{blue}{x - \mathsf{fma}\left(a, t - 1, z \cdot \left(y - 1\right)\right)} \]

   9. Taylor expanded in y around inf

      \[\leadsto x - y \cdot \color{blue}{z} \]
   10. Step-by-step derivation

   11. Applied rewrites 48.6%

       \[\leadsto x - y \cdot \color{blue}{z} \]

   if -5.2e-215 < t < 4.2999999999999998e42

   1. Initial program 98.6%

      \[\left(\left(x - \left(y - 1\right) \cdot z\right) - \left(t - 1\right) \cdot a\right) + \left(\left(y + t\right) - 2\right) \cdot b \]
   2. Add Preprocessing

   3. Taylor expanded in y around inf

      \[\leadsto \color{blue}{y \cdot \left(b - z\right)} \]
   4. Step-by-step derivation

      1. *-commutative N/A

         \[\leadsto \color{blue}{\left(b - z\right) \cdot y} \]

      2. lower-*.f64 N/A

         \[\leadsto \color{blue}{\left(b - z\right) \cdot y} \]

      3. lower--.f64 41.8

         \[\leadsto \color{blue}{\left(b - z\right)} \cdot y \]

   5. Applied rewrites 41.8%

      \[\leadsto \color{blue}{\left(b - z\right) \cdot y} \]
3. Recombined 3 regimes into one program.
4. Add Preprocessing

Alternative 8: 43.5% accurate, 1.4× speedup

Math

\[\begin{array}{l} \\ \begin{array}{l} t_1 := \left(b - a\right) \cdot t\\ \mathbf{if}\;t \leq -26500000000:\\ \;\;\;\;t\_1\\ \mathbf{elif}\;t \leq -1.75 \cdot 10^{-242}:\\ \;\;\;\;\left(1 - y\right) \cdot z\\ \mathbf{elif}\;t \leq 1.35 \cdot 10^{+41}:\\ \;\;\;\;b \cdot y\\ \mathbf{else}:\\ \;\;\;\;t\_1\\ \end{array} \end{array} \]

FPCore

(FPCore (x y z t a b)
 :precision binary64
 (let* ((t_1 (* (- b a) t)))
   (if (<= t -26500000000.0)
     t_1
     (if (<= t -1.75e-242) (* (- 1.0 y) z) (if (<= t 1.35e+41) (* b y) t_1)))))

C

double code(double x, double y, double z, double t, double a, double b) {
	double t_1 = (b - a) * t;
	double tmp;
	if (t <= -26500000000.0) {
		tmp = t_1;
	} else if (t <= -1.75e-242) {
		tmp = (1.0 - y) * z;
	} else if (t <= 1.35e+41) {
		tmp = b * y;
	} else {
		tmp = t_1;
	}
	return tmp;
}

Fortran

real(8) function code(x, y, z, t, a, b)
    real(8), intent (in) :: x
    real(8), intent (in) :: y
    real(8), intent (in) :: z
    real(8), intent (in) :: t
    real(8), intent (in) :: a
    real(8), intent (in) :: b
    real(8) :: t_1
    real(8) :: tmp
    t_1 = (b - a) * t
    if (t <= (-26500000000.0d0)) then
        tmp = t_1
    else if (t <= (-1.75d-242)) then
        tmp = (1.0d0 - y) * z
    else if (t <= 1.35d+41) then
        tmp = b * y
    else
        tmp = t_1
    end if
    code = tmp
end function

Java

public static double code(double x, double y, double z, double t, double a, double b) {
	double t_1 = (b - a) * t;
	double tmp;
	if (t <= -26500000000.0) {
		tmp = t_1;
	} else if (t <= -1.75e-242) {
		tmp = (1.0 - y) * z;
	} else if (t <= 1.35e+41) {
		tmp = b * y;
	} else {
		tmp = t_1;
	}
	return tmp;
}

Python

def code(x, y, z, t, a, b):
	t_1 = (b - a) * t
	tmp = 0
	if t <= -26500000000.0:
		tmp = t_1
	elif t <= -1.75e-242:
		tmp = (1.0 - y) * z
	elif t <= 1.35e+41:
		tmp = b * y
	else:
		tmp = t_1
	return tmp

Julia

function code(x, y, z, t, a, b)
	t_1 = Float64(Float64(b - a) * t)
	tmp = 0.0
	if (t <= -26500000000.0)
		tmp = t_1;
	elseif (t <= -1.75e-242)
		tmp = Float64(Float64(1.0 - y) * z);
	elseif (t <= 1.35e+41)
		tmp = Float64(b * y);
	else
		tmp = t_1;
	end
	return tmp
end

MATLAB

function tmp_2 = code(x, y, z, t, a, b)
	t_1 = (b - a) * t;
	tmp = 0.0;
	if (t <= -26500000000.0)
		tmp = t_1;
	elseif (t <= -1.75e-242)
		tmp = (1.0 - y) * z;
	elseif (t <= 1.35e+41)
		tmp = b * y;
	else
		tmp = t_1;
	end
	tmp_2 = tmp;
end

Wolfram

code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(N[(b - a), $MachinePrecision] * t), $MachinePrecision]}, If[LessEqual[t, -26500000000.0], t$95$1, If[LessEqual[t, -1.75e-242], N[(N[(1.0 - y), $MachinePrecision] * z), $MachinePrecision], If[LessEqual[t, 1.35e+41], N[(b * y), $MachinePrecision], t$95$1]]]]

Derivation

1. Split input into 3 regimes

2. if t < -2.65e10 or 1.35e41 < t

   1. Initial program 92.8%

      \[\left(\left(x - \left(y - 1\right) \cdot z\right) - \left(t - 1\right) \cdot a\right) + \left(\left(y + t\right) - 2\right) \cdot b \]
   2. Add Preprocessing

   3. Taylor expanded in t around inf

      \[\leadsto \color{blue}{t \cdot \left(b - a\right)} \]
   4. Step-by-step derivation

      1. *-commutative N/A

         \[\leadsto \color{blue}{\left(b - a\right) \cdot t} \]

      2. lower-*.f64 N/A

         \[\leadsto \color{blue}{\left(b - a\right) \cdot t} \]

      3. lower--.f64 74.2

         \[\leadsto \color{blue}{\left(b - a\right)} \cdot t \]

   5. Applied rewrites 74.2%

      \[\leadsto \color{blue}{\left(b - a\right) \cdot t} \]

   if -2.65e10 < t < -1.7499999999999999e-242

   1. Initial program 96.9%

      \[\left(\left(x - \left(y - 1\right) \cdot z\right) - \left(t - 1\right) \cdot a\right) + \left(\left(y + t\right) - 2\right) \cdot b \]
   2. Add Preprocessing

   3. Taylor expanded in z around inf

      \[\leadsto \color{blue}{z \cdot \left(1 - y\right)} \]
   4. Step-by-step derivation

      1. *-commutative N/A

         \[\leadsto \color{blue}{\left(1 - y\right) \cdot z} \]

      2. sub-neg N/A

         \[\leadsto \color{blue}{\left(1 + \left(\mathsf{neg}\left(y\right)\right)\right)} \cdot z \]

      3. neg-mul-1 N/A

         \[\leadsto \left(1 + \color{blue}{-1 \cdot y}\right) \cdot z \]

      4. +-commutative N/A

         \[\leadsto \color{blue}{\left(-1 \cdot y + 1\right)} \cdot z \]

      5. metadata-eval N/A

         \[\leadsto \left(-1 \cdot y + \color{blue}{-1 \cdot -1}\right) \cdot z \]

      6. distribute-lft-in N/A

         \[\leadsto \color{blue}{\left(-1 \cdot \left(y + -1\right)\right)} \cdot z \]

      7. metadata-eval N/A

         \[\leadsto \left(-1 \cdot \left(y + \color{blue}{\left(\mathsf{neg}\left(1\right)\right)}\right)\right) \cdot z \]

      8. sub-neg N/A

         \[\leadsto \left(-1 \cdot \color{blue}{\left(y - 1\right)}\right) \cdot z \]

      9. lower-*.f64 N/A

         \[\leadsto \color{blue}{\left(-1 \cdot \left(y - 1\right)\right) \cdot z} \]

      10. sub-neg N/A

          \[\leadsto \left(-1 \cdot \color{blue}{\left(y + \left(\mathsf{neg}\left(1\right)\right)\right)}\right) \cdot z \]

      11. metadata-eval N/A

          \[\leadsto \left(-1 \cdot \left(y + \color{blue}{-1}\right)\right) \cdot z \]

      12. distribute-lft-in N/A

          \[\leadsto \color{blue}{\left(-1 \cdot y + -1 \cdot -1\right)} \cdot z \]

      13. metadata-eval N/A

          \[\leadsto \left(-1 \cdot y + \color{blue}{1}\right) \cdot z \]

      14. +-commutative N/A

          \[\leadsto \color{blue}{\left(1 + -1 \cdot y\right)} \cdot z \]

      15. neg-mul-1 N/A

          \[\leadsto \left(1 + \color{blue}{\left(\mathsf{neg}\left(y\right)\right)}\right) \cdot z \]

      16. sub-neg N/A

          \[\leadsto \color{blue}{\left(1 - y\right)} \cdot z \]

      17. lower--.f64 40.9

          \[\leadsto \color{blue}{\left(1 - y\right)} \cdot z \]

   5. Applied rewrites 40.9%

      \[\leadsto \color{blue}{\left(1 - y\right) \cdot z} \]

   if -1.7499999999999999e-242 < t < 1.35e41

   1. Initial program 98.5%

      \[\left(\left(x - \left(y - 1\right) \cdot z\right) - \left(t - 1\right) \cdot a\right) + \left(\left(y + t\right) - 2\right) \cdot b \]
   2. Add Preprocessing

   3. Taylor expanded in y around inf

      \[\leadsto \color{blue}{y \cdot \left(b - z\right)} \]
   4. Step-by-step derivation

      1. *-commutative N/A

         \[\leadsto \color{blue}{\left(b - z\right) \cdot y} \]

      2. lower-*.f64 N/A

         \[\leadsto \color{blue}{\left(b - z\right) \cdot y} \]

      3. lower--.f64 41.3

         \[\leadsto \color{blue}{\left(b - z\right)} \cdot y \]

   5. Applied rewrites 41.3%

      \[\leadsto \color{blue}{\left(b - z\right) \cdot y} \]

   6. Taylor expanded in z around 0

      \[\leadsto b \cdot \color{blue}{y} \]
   7. Step-by-step derivation

   8. Applied rewrites 30.7%

      \[\leadsto b \cdot \color{blue}{y} \]
3. Recombined 3 regimes into one program.
4. Add Preprocessing

Alternative 9: 39.2% accurate, 1.4× speedup

Math

\[\begin{array}{l} \\ \begin{array}{l} t_1 := \left(t - 2\right) \cdot b\\ \mathbf{if}\;b \leq -5.6 \cdot 10^{+142}:\\ \;\;\;\;t\_1\\ \mathbf{elif}\;b \leq -8.5 \cdot 10^{+30}:\\ \;\;\;\;b \cdot y\\ \mathbf{elif}\;b \leq 1.2 \cdot 10^{+111}:\\ \;\;\;\;\left(1 - t\right) \cdot a\\ \mathbf{else}:\\ \;\;\;\;t\_1\\ \end{array} \end{array} \]

FPCore

(FPCore (x y z t a b)
 :precision binary64
 (let* ((t_1 (* (- t 2.0) b)))
   (if (<= b -5.6e+142)
     t_1
     (if (<= b -8.5e+30) (* b y) (if (<= b 1.2e+111) (* (- 1.0 t) a) t_1)))))

C

double code(double x, double y, double z, double t, double a, double b) {
	double t_1 = (t - 2.0) * b;
	double tmp;
	if (b <= -5.6e+142) {
		tmp = t_1;
	} else if (b <= -8.5e+30) {
		tmp = b * y;
	} else if (b <= 1.2e+111) {
		tmp = (1.0 - t) * a;
	} else {
		tmp = t_1;
	}
	return tmp;
}

Fortran

real(8) function code(x, y, z, t, a, b)
    real(8), intent (in) :: x
    real(8), intent (in) :: y
    real(8), intent (in) :: z
    real(8), intent (in) :: t
    real(8), intent (in) :: a
    real(8), intent (in) :: b
    real(8) :: t_1
    real(8) :: tmp
    t_1 = (t - 2.0d0) * b
    if (b <= (-5.6d+142)) then
        tmp = t_1
    else if (b <= (-8.5d+30)) then
        tmp = b * y
    else if (b <= 1.2d+111) then
        tmp = (1.0d0 - t) * a
    else
        tmp = t_1
    end if
    code = tmp
end function

Java

public static double code(double x, double y, double z, double t, double a, double b) {
	double t_1 = (t - 2.0) * b;
	double tmp;
	if (b <= -5.6e+142) {
		tmp = t_1;
	} else if (b <= -8.5e+30) {
		tmp = b * y;
	} else if (b <= 1.2e+111) {
		tmp = (1.0 - t) * a;
	} else {
		tmp = t_1;
	}
	return tmp;
}

Python

def code(x, y, z, t, a, b):
	t_1 = (t - 2.0) * b
	tmp = 0
	if b <= -5.6e+142:
		tmp = t_1
	elif b <= -8.5e+30:
		tmp = b * y
	elif b <= 1.2e+111:
		tmp = (1.0 - t) * a
	else:
		tmp = t_1
	return tmp

Julia

function code(x, y, z, t, a, b)
	t_1 = Float64(Float64(t - 2.0) * b)
	tmp = 0.0
	if (b <= -5.6e+142)
		tmp = t_1;
	elseif (b <= -8.5e+30)
		tmp = Float64(b * y);
	elseif (b <= 1.2e+111)
		tmp = Float64(Float64(1.0 - t) * a);
	else
		tmp = t_1;
	end
	return tmp
end

MATLAB

function tmp_2 = code(x, y, z, t, a, b)
	t_1 = (t - 2.0) * b;
	tmp = 0.0;
	if (b <= -5.6e+142)
		tmp = t_1;
	elseif (b <= -8.5e+30)
		tmp = b * y;
	elseif (b <= 1.2e+111)
		tmp = (1.0 - t) * a;
	else
		tmp = t_1;
	end
	tmp_2 = tmp;
end

Wolfram

code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(N[(t - 2.0), $MachinePrecision] * b), $MachinePrecision]}, If[LessEqual[b, -5.6e+142], t$95$1, If[LessEqual[b, -8.5e+30], N[(b * y), $MachinePrecision], If[LessEqual[b, 1.2e+111], N[(N[(1.0 - t), $MachinePrecision] * a), $MachinePrecision], t$95$1]]]]

Derivation

1. Split input into 3 regimes

2. if b < -5.6e142 or 1.20000000000000003e111 < b

   1. Initial program 86.5%

      \[\left(\left(x - \left(y - 1\right) \cdot z\right) - \left(t - 1\right) \cdot a\right) + \left(\left(y + t\right) - 2\right) \cdot b \]
   2. Add Preprocessing

   3. Taylor expanded in y around 0

      \[\leadsto \color{blue}{\left(x + b \cdot \left(t - 2\right)\right) - \left(-1 \cdot z + a \cdot \left(t - 1\right)\right)} \]
   4. Step-by-step derivation

      1. sub-neg N/A

         \[\leadsto \color{blue}{\left(x + b \cdot \left(t - 2\right)\right) + \left(\mathsf{neg}\left(\left(-1 \cdot z + a \cdot \left(t - 1\right)\right)\right)\right)} \]

      2. +-commutative N/A

         \[\leadsto \color{blue}{\left(\mathsf{neg}\left(\left(-1 \cdot z + a \cdot \left(t - 1\right)\right)\right)\right) + \left(x + b \cdot \left(t - 2\right)\right)} \]

      3. +-commutative N/A

         \[\leadsto \left(\mathsf{neg}\left(\color{blue}{\left(a \cdot \left(t - 1\right) + -1 \cdot z\right)}\right)\right) + \left(x + b \cdot \left(t - 2\right)\right) \]

      4. distribute-neg-in N/A

         \[\leadsto \color{blue}{\left(\left(\mathsf{neg}\left(a \cdot \left(t - 1\right)\right)\right) + \left(\mathsf{neg}\left(-1 \cdot z\right)\right)\right)} + \left(x + b \cdot \left(t - 2\right)\right) \]

      5. mul-1-neg N/A

         \[\leadsto \left(\left(\mathsf{neg}\left(a \cdot \left(t - 1\right)\right)\right) + \left(\mathsf{neg}\left(\color{blue}{\left(\mathsf{neg}\left(z\right)\right)}\right)\right)\right) + \left(x + b \cdot \left(t - 2\right)\right) \]

      6. remove-double-neg N/A

         \[\leadsto \left(\left(\mathsf{neg}\left(a \cdot \left(t - 1\right)\right)\right) + \color{blue}{z}\right) + \left(x + b \cdot \left(t - 2\right)\right) \]

      7. associate-+l+ N/A

         \[\leadsto \color{blue}{\left(\mathsf{neg}\left(a \cdot \left(t - 1\right)\right)\right) + \left(z + \left(x + b \cdot \left(t - 2\right)\right)\right)} \]

      8. *-commutative N/A

         \[\leadsto \left(\mathsf{neg}\left(\color{blue}{\left(t - 1\right) \cdot a}\right)\right) + \left(z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      9. distribute-lft-neg-in N/A

         \[\leadsto \color{blue}{\left(\mathsf{neg}\left(\left(t - 1\right)\right)\right) \cdot a} + \left(z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      10. mul-1-neg N/A

          \[\leadsto \color{blue}{\left(-1 \cdot \left(t - 1\right)\right)} \cdot a + \left(z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      11. lower-fma.f64 N/A

          \[\leadsto \color{blue}{\mathsf{fma}\left(-1 \cdot \left(t - 1\right), a, z + \left(x + b \cdot \left(t - 2\right)\right)\right)} \]

      12. sub-neg N/A

          \[\leadsto \mathsf{fma}\left(-1 \cdot \color{blue}{\left(t + \left(\mathsf{neg}\left(1\right)\right)\right)}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      13. metadata-eval N/A

          \[\leadsto \mathsf{fma}\left(-1 \cdot \left(t + \color{blue}{-1}\right), a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      14. distribute-lft-in N/A

          \[\leadsto \mathsf{fma}\left(\color{blue}{-1 \cdot t + -1 \cdot -1}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      15. metadata-eval N/A

          \[\leadsto \mathsf{fma}\left(-1 \cdot t + \color{blue}{1}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      16. +-commutative N/A

          \[\leadsto \mathsf{fma}\left(\color{blue}{1 + -1 \cdot t}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      17. neg-mul-1 N/A

          \[\leadsto \mathsf{fma}\left(1 + \color{blue}{\left(\mathsf{neg}\left(t\right)\right)}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      18. sub-neg N/A

          \[\leadsto \mathsf{fma}\left(\color{blue}{1 - t}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      19. lower--.f64 N/A

          \[\leadsto \mathsf{fma}\left(\color{blue}{1 - t}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      20. lower-+.f64 N/A

          \[\leadsto \mathsf{fma}\left(1 - t, a, \color{blue}{z + \left(x + b \cdot \left(t - 2\right)\right)}\right) \]

   5. Applied rewrites 71.3%

      \[\leadsto \color{blue}{\mathsf{fma}\left(1 - t, a, z + \mathsf{fma}\left(t - 2, b, x\right)\right)} \]

   6. Taylor expanded in b around inf

      \[\leadsto \color{blue}{b \cdot \left(\left(t + y\right) - 2\right)} \]
   7. Step-by-step derivation

      1. *-commutative N/A

         \[\leadsto \color{blue}{\left(\left(t + y\right) - 2\right) \cdot b} \]

      2. lower-*.f64 N/A

         \[\leadsto \color{blue}{\left(\left(t + y\right) - 2\right) \cdot b} \]

      3. lower--.f64 N/A

         \[\leadsto \color{blue}{\left(\left(t + y\right) - 2\right)} \cdot b \]

      4. lower-+.f64 79.5

         \[\leadsto \left(\color{blue}{\left(t + y\right)} - 2\right) \cdot b \]

   8. Applied rewrites 79.5%

      \[\leadsto \color{blue}{\left(\left(t + y\right) - 2\right) \cdot b} \]

   9. Taylor expanded in y around 0

      \[\leadsto \left(t - 2\right) \cdot b \]
   10. Step-by-step derivation

   11. Applied rewrites 53.7%

       \[\leadsto \left(t - 2\right) \cdot b \]

   if -5.6e142 < b < -8.4999999999999995e30

   1. Initial program 100.0%

      \[\left(\left(x - \left(y - 1\right) \cdot z\right) - \left(t - 1\right) \cdot a\right) + \left(\left(y + t\right) - 2\right) \cdot b \]
   2. Add Preprocessing

   3. Taylor expanded in y around inf

      \[\leadsto \color{blue}{y \cdot \left(b - z\right)} \]
   4. Step-by-step derivation

      1. *-commutative N/A

         \[\leadsto \color{blue}{\left(b - z\right) \cdot y} \]

      2. lower-*.f64 N/A

         \[\leadsto \color{blue}{\left(b - z\right) \cdot y} \]

      3. lower--.f64 59.3

         \[\leadsto \color{blue}{\left(b - z\right)} \cdot y \]

   5. Applied rewrites 59.3%

      \[\leadsto \color{blue}{\left(b - z\right) \cdot y} \]

   6. Taylor expanded in z around 0

      \[\leadsto b \cdot \color{blue}{y} \]
   7. Step-by-step derivation

   8. Applied rewrites 49.4%

      \[\leadsto b \cdot \color{blue}{y} \]

   if -8.4999999999999995e30 < b < 1.20000000000000003e111

   1. Initial program 98.8%

      \[\left(\left(x - \left(y - 1\right) \cdot z\right) - \left(t - 1\right) \cdot a\right) + \left(\left(y + t\right) - 2\right) \cdot b \]
   2. Add Preprocessing

   3. Taylor expanded in a around inf

      \[\leadsto \color{blue}{a \cdot \left(1 - t\right)} \]
   4. Step-by-step derivation

      1. *-commutative N/A

         \[\leadsto \color{blue}{\left(1 - t\right) \cdot a} \]

      2. sub-neg N/A

         \[\leadsto \color{blue}{\left(1 + \left(\mathsf{neg}\left(t\right)\right)\right)} \cdot a \]

      3. neg-mul-1 N/A

         \[\leadsto \left(1 + \color{blue}{-1 \cdot t}\right) \cdot a \]

      4. +-commutative N/A

         \[\leadsto \color{blue}{\left(-1 \cdot t + 1\right)} \cdot a \]

      5. metadata-eval N/A

         \[\leadsto \left(-1 \cdot t + \color{blue}{-1 \cdot -1}\right) \cdot a \]

      6. distribute-lft-in N/A

         \[\leadsto \color{blue}{\left(-1 \cdot \left(t + -1\right)\right)} \cdot a \]

      7. metadata-eval N/A

         \[\leadsto \left(-1 \cdot \left(t + \color{blue}{\left(\mathsf{neg}\left(1\right)\right)}\right)\right) \cdot a \]

      8. sub-neg N/A

         \[\leadsto \left(-1 \cdot \color{blue}{\left(t - 1\right)}\right) \cdot a \]

      9. lower-*.f64 N/A

         \[\leadsto \color{blue}{\left(-1 \cdot \left(t - 1\right)\right) \cdot a} \]

      10. sub-neg N/A

          \[\leadsto \left(-1 \cdot \color{blue}{\left(t + \left(\mathsf{neg}\left(1\right)\right)\right)}\right) \cdot a \]

      11. metadata-eval N/A

          \[\leadsto \left(-1 \cdot \left(t + \color{blue}{-1}\right)\right) \cdot a \]

      12. distribute-lft-in N/A

          \[\leadsto \color{blue}{\left(-1 \cdot t + -1 \cdot -1\right)} \cdot a \]

      13. metadata-eval N/A

          \[\leadsto \left(-1 \cdot t + \color{blue}{1}\right) \cdot a \]

      14. +-commutative N/A

          \[\leadsto \color{blue}{\left(1 + -1 \cdot t\right)} \cdot a \]

      15. neg-mul-1 N/A

          \[\leadsto \left(1 + \color{blue}{\left(\mathsf{neg}\left(t\right)\right)}\right) \cdot a \]

      16. sub-neg N/A

          \[\leadsto \color{blue}{\left(1 - t\right)} \cdot a \]

      17. lower--.f64 46.2

          \[\leadsto \color{blue}{\left(1 - t\right)} \cdot a \]

   5. Applied rewrites 46.2%

      \[\leadsto \color{blue}{\left(1 - t\right) \cdot a} \]
3. Recombined 3 regimes into one program.
4. Add Preprocessing

Alternative 10: 31.6% accurate, 1.4× speedup

Math

\[\begin{array}{l} \\ \begin{array}{l} t_1 := \left(t - 2\right) \cdot b\\ \mathbf{if}\;b \leq -5.6 \cdot 10^{+142}:\\ \;\;\;\;t\_1\\ \mathbf{elif}\;b \leq -9.4 \cdot 10^{+26}:\\ \;\;\;\;b \cdot y\\ \mathbf{elif}\;b \leq 2.6 \cdot 10^{+110}:\\ \;\;\;\;\left(-a\right) \cdot t\\ \mathbf{else}:\\ \;\;\;\;t\_1\\ \end{array} \end{array} \]

FPCore

(FPCore (x y z t a b)
 :precision binary64
 (let* ((t_1 (* (- t 2.0) b)))
   (if (<= b -5.6e+142)
     t_1
     (if (<= b -9.4e+26) (* b y) (if (<= b 2.6e+110) (* (- a) t) t_1)))))

C

double code(double x, double y, double z, double t, double a, double b) {
	double t_1 = (t - 2.0) * b;
	double tmp;
	if (b <= -5.6e+142) {
		tmp = t_1;
	} else if (b <= -9.4e+26) {
		tmp = b * y;
	} else if (b <= 2.6e+110) {
		tmp = -a * t;
	} else {
		tmp = t_1;
	}
	return tmp;
}

Fortran

real(8) function code(x, y, z, t, a, b)
    real(8), intent (in) :: x
    real(8), intent (in) :: y
    real(8), intent (in) :: z
    real(8), intent (in) :: t
    real(8), intent (in) :: a
    real(8), intent (in) :: b
    real(8) :: t_1
    real(8) :: tmp
    t_1 = (t - 2.0d0) * b
    if (b <= (-5.6d+142)) then
        tmp = t_1
    else if (b <= (-9.4d+26)) then
        tmp = b * y
    else if (b <= 2.6d+110) then
        tmp = -a * t
    else
        tmp = t_1
    end if
    code = tmp
end function

Java

public static double code(double x, double y, double z, double t, double a, double b) {
	double t_1 = (t - 2.0) * b;
	double tmp;
	if (b <= -5.6e+142) {
		tmp = t_1;
	} else if (b <= -9.4e+26) {
		tmp = b * y;
	} else if (b <= 2.6e+110) {
		tmp = -a * t;
	} else {
		tmp = t_1;
	}
	return tmp;
}

Python

def code(x, y, z, t, a, b):
	t_1 = (t - 2.0) * b
	tmp = 0
	if b <= -5.6e+142:
		tmp = t_1
	elif b <= -9.4e+26:
		tmp = b * y
	elif b <= 2.6e+110:
		tmp = -a * t
	else:
		tmp = t_1
	return tmp

Julia

function code(x, y, z, t, a, b)
	t_1 = Float64(Float64(t - 2.0) * b)
	tmp = 0.0
	if (b <= -5.6e+142)
		tmp = t_1;
	elseif (b <= -9.4e+26)
		tmp = Float64(b * y);
	elseif (b <= 2.6e+110)
		tmp = Float64(Float64(-a) * t);
	else
		tmp = t_1;
	end
	return tmp
end

MATLAB

function tmp_2 = code(x, y, z, t, a, b)
	t_1 = (t - 2.0) * b;
	tmp = 0.0;
	if (b <= -5.6e+142)
		tmp = t_1;
	elseif (b <= -9.4e+26)
		tmp = b * y;
	elseif (b <= 2.6e+110)
		tmp = -a * t;
	else
		tmp = t_1;
	end
	tmp_2 = tmp;
end

Wolfram

code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(N[(t - 2.0), $MachinePrecision] * b), $MachinePrecision]}, If[LessEqual[b, -5.6e+142], t$95$1, If[LessEqual[b, -9.4e+26], N[(b * y), $MachinePrecision], If[LessEqual[b, 2.6e+110], N[((-a) * t), $MachinePrecision], t$95$1]]]]

Derivation

1. Split input into 3 regimes

2. if b < -5.6e142 or 2.6e110 < b

   1. Initial program 86.5%

      \[\left(\left(x - \left(y - 1\right) \cdot z\right) - \left(t - 1\right) \cdot a\right) + \left(\left(y + t\right) - 2\right) \cdot b \]
   2. Add Preprocessing

   3. Taylor expanded in y around 0

      \[\leadsto \color{blue}{\left(x + b \cdot \left(t - 2\right)\right) - \left(-1 \cdot z + a \cdot \left(t - 1\right)\right)} \]
   4. Step-by-step derivation

      1. sub-neg N/A

         \[\leadsto \color{blue}{\left(x + b \cdot \left(t - 2\right)\right) + \left(\mathsf{neg}\left(\left(-1 \cdot z + a \cdot \left(t - 1\right)\right)\right)\right)} \]

      2. +-commutative N/A

         \[\leadsto \color{blue}{\left(\mathsf{neg}\left(\left(-1 \cdot z + a \cdot \left(t - 1\right)\right)\right)\right) + \left(x + b \cdot \left(t - 2\right)\right)} \]

      3. +-commutative N/A

         \[\leadsto \left(\mathsf{neg}\left(\color{blue}{\left(a \cdot \left(t - 1\right) + -1 \cdot z\right)}\right)\right) + \left(x + b \cdot \left(t - 2\right)\right) \]

      4. distribute-neg-in N/A

         \[\leadsto \color{blue}{\left(\left(\mathsf{neg}\left(a \cdot \left(t - 1\right)\right)\right) + \left(\mathsf{neg}\left(-1 \cdot z\right)\right)\right)} + \left(x + b \cdot \left(t - 2\right)\right) \]

      5. mul-1-neg N/A

         \[\leadsto \left(\left(\mathsf{neg}\left(a \cdot \left(t - 1\right)\right)\right) + \left(\mathsf{neg}\left(\color{blue}{\left(\mathsf{neg}\left(z\right)\right)}\right)\right)\right) + \left(x + b \cdot \left(t - 2\right)\right) \]

      6. remove-double-neg N/A

         \[\leadsto \left(\left(\mathsf{neg}\left(a \cdot \left(t - 1\right)\right)\right) + \color{blue}{z}\right) + \left(x + b \cdot \left(t - 2\right)\right) \]

      7. associate-+l+ N/A

         \[\leadsto \color{blue}{\left(\mathsf{neg}\left(a \cdot \left(t - 1\right)\right)\right) + \left(z + \left(x + b \cdot \left(t - 2\right)\right)\right)} \]

      8. *-commutative N/A

         \[\leadsto \left(\mathsf{neg}\left(\color{blue}{\left(t - 1\right) \cdot a}\right)\right) + \left(z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      9. distribute-lft-neg-in N/A

         \[\leadsto \color{blue}{\left(\mathsf{neg}\left(\left(t - 1\right)\right)\right) \cdot a} + \left(z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      10. mul-1-neg N/A

          \[\leadsto \color{blue}{\left(-1 \cdot \left(t - 1\right)\right)} \cdot a + \left(z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      11. lower-fma.f64 N/A

          \[\leadsto \color{blue}{\mathsf{fma}\left(-1 \cdot \left(t - 1\right), a, z + \left(x + b \cdot \left(t - 2\right)\right)\right)} \]

      12. sub-neg N/A

          \[\leadsto \mathsf{fma}\left(-1 \cdot \color{blue}{\left(t + \left(\mathsf{neg}\left(1\right)\right)\right)}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      13. metadata-eval N/A

          \[\leadsto \mathsf{fma}\left(-1 \cdot \left(t + \color{blue}{-1}\right), a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      14. distribute-lft-in N/A

          \[\leadsto \mathsf{fma}\left(\color{blue}{-1 \cdot t + -1 \cdot -1}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      15. metadata-eval N/A

          \[\leadsto \mathsf{fma}\left(-1 \cdot t + \color{blue}{1}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      16. +-commutative N/A

          \[\leadsto \mathsf{fma}\left(\color{blue}{1 + -1 \cdot t}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      17. neg-mul-1 N/A

          \[\leadsto \mathsf{fma}\left(1 + \color{blue}{\left(\mathsf{neg}\left(t\right)\right)}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      18. sub-neg N/A

          \[\leadsto \mathsf{fma}\left(\color{blue}{1 - t}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      19. lower--.f64 N/A

          \[\leadsto \mathsf{fma}\left(\color{blue}{1 - t}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      20. lower-+.f64 N/A

          \[\leadsto \mathsf{fma}\left(1 - t, a, \color{blue}{z + \left(x + b \cdot \left(t - 2\right)\right)}\right) \]

   5. Applied rewrites 71.3%

      \[\leadsto \color{blue}{\mathsf{fma}\left(1 - t, a, z + \mathsf{fma}\left(t - 2, b, x\right)\right)} \]

   6. Taylor expanded in b around inf

      \[\leadsto \color{blue}{b \cdot \left(\left(t + y\right) - 2\right)} \]
   7. Step-by-step derivation

      1. *-commutative N/A

         \[\leadsto \color{blue}{\left(\left(t + y\right) - 2\right) \cdot b} \]

      2. lower-*.f64 N/A

         \[\leadsto \color{blue}{\left(\left(t + y\right) - 2\right) \cdot b} \]

      3. lower--.f64 N/A

         \[\leadsto \color{blue}{\left(\left(t + y\right) - 2\right)} \cdot b \]

      4. lower-+.f64 79.5

         \[\leadsto \left(\color{blue}{\left(t + y\right)} - 2\right) \cdot b \]

   8. Applied rewrites 79.5%

      \[\leadsto \color{blue}{\left(\left(t + y\right) - 2\right) \cdot b} \]

   9. Taylor expanded in y around 0

      \[\leadsto \left(t - 2\right) \cdot b \]
   10. Step-by-step derivation

   11. Applied rewrites 53.7%

       \[\leadsto \left(t - 2\right) \cdot b \]

   if -5.6e142 < b < -9.3999999999999995e26

   1. Initial program 100.0%

      \[\left(\left(x - \left(y - 1\right) \cdot z\right) - \left(t - 1\right) \cdot a\right) + \left(\left(y + t\right) - 2\right) \cdot b \]
   2. Add Preprocessing

   3. Taylor expanded in y around inf

      \[\leadsto \color{blue}{y \cdot \left(b - z\right)} \]
   4. Step-by-step derivation

      1. *-commutative N/A

         \[\leadsto \color{blue}{\left(b - z\right) \cdot y} \]

      2. lower-*.f64 N/A

         \[\leadsto \color{blue}{\left(b - z\right) \cdot y} \]

      3. lower--.f64 59.3

         \[\leadsto \color{blue}{\left(b - z\right)} \cdot y \]

   5. Applied rewrites 59.3%

      \[\leadsto \color{blue}{\left(b - z\right) \cdot y} \]

   6. Taylor expanded in z around 0

      \[\leadsto b \cdot \color{blue}{y} \]
   7. Step-by-step derivation

   8. Applied rewrites 49.4%

      \[\leadsto b \cdot \color{blue}{y} \]

   if -9.3999999999999995e26 < b < 2.6e110

   1. Initial program 98.8%

      \[\left(\left(x - \left(y - 1\right) \cdot z\right) - \left(t - 1\right) \cdot a\right) + \left(\left(y + t\right) - 2\right) \cdot b \]
   2. Add Preprocessing

   3. Taylor expanded in t around inf

      \[\leadsto \color{blue}{t \cdot \left(b - a\right)} \]
   4. Step-by-step derivation

      1. *-commutative N/A

         \[\leadsto \color{blue}{\left(b - a\right) \cdot t} \]

      2. lower-*.f64 N/A

         \[\leadsto \color{blue}{\left(b - a\right) \cdot t} \]

      3. lower--.f64 39.0

         \[\leadsto \color{blue}{\left(b - a\right)} \cdot t \]

   5. Applied rewrites 39.0%

      \[\leadsto \color{blue}{\left(b - a\right) \cdot t} \]

   6. Taylor expanded in a around inf

      \[\leadsto \left(-1 \cdot a\right) \cdot t \]
   7. Step-by-step derivation

   8. Applied rewrites 35.8%

      \[\leadsto \left(-a\right) \cdot t \]
3. Recombined 3 regimes into one program.
4. Add Preprocessing

Alternative 11: 68.9% accurate, 1.5× speedup

Math

\[\begin{array}{l} \\ \begin{array}{l} \mathbf{if}\;b \leq -4 \cdot 10^{+36} \lor \neg \left(b \leq 4 \cdot 10^{+109}\right):\\ \;\;\;\;\left(\left(t + y\right) - 2\right) \cdot b\\ \mathbf{else}:\\ \;\;\;\;\mathsf{fma}\left(1 - t, a, z + x\right)\\ \end{array} \end{array} \]

FPCore

(FPCore (x y z t a b)
 :precision binary64
 (if (or (<= b -4e+36) (not (<= b 4e+109)))
   (* (- (+ t y) 2.0) b)
   (fma (- 1.0 t) a (+ z x))))

C

double code(double x, double y, double z, double t, double a, double b) {
	double tmp;
	if ((b <= -4e+36) || !(b <= 4e+109)) {
		tmp = ((t + y) - 2.0) * b;
	} else {
		tmp = fma((1.0 - t), a, (z + x));
	}
	return tmp;
}

Julia

function code(x, y, z, t, a, b)
	tmp = 0.0
	if ((b <= -4e+36) || !(b <= 4e+109))
		tmp = Float64(Float64(Float64(t + y) - 2.0) * b);
	else
		tmp = fma(Float64(1.0 - t), a, Float64(z + x));
	end
	return tmp
end

Wolfram

code[x_, y_, z_, t_, a_, b_] := If[Or[LessEqual[b, -4e+36], N[Not[LessEqual[b, 4e+109]], $MachinePrecision]], N[(N[(N[(t + y), $MachinePrecision] - 2.0), $MachinePrecision] * b), $MachinePrecision], N[(N[(1.0 - t), $MachinePrecision] * a + N[(z + x), $MachinePrecision]), $MachinePrecision]]

Derivation

1. Split input into 2 regimes

2. if b < -4.00000000000000017e36 or 3.99999999999999993e109 < b

   1. Initial program 89.2%

      \[\left(\left(x - \left(y - 1\right) \cdot z\right) - \left(t - 1\right) \cdot a\right) + \left(\left(y + t\right) - 2\right) \cdot b \]
   2. Add Preprocessing

   3. Taylor expanded in y around 0

      \[\leadsto \color{blue}{\left(x + b \cdot \left(t - 2\right)\right) - \left(-1 \cdot z + a \cdot \left(t - 1\right)\right)} \]
   4. Step-by-step derivation

      1. sub-neg N/A

         \[\leadsto \color{blue}{\left(x + b \cdot \left(t - 2\right)\right) + \left(\mathsf{neg}\left(\left(-1 \cdot z + a \cdot \left(t - 1\right)\right)\right)\right)} \]

      2. +-commutative N/A

         \[\leadsto \color{blue}{\left(\mathsf{neg}\left(\left(-1 \cdot z + a \cdot \left(t - 1\right)\right)\right)\right) + \left(x + b \cdot \left(t - 2\right)\right)} \]

      3. +-commutative N/A

         \[\leadsto \left(\mathsf{neg}\left(\color{blue}{\left(a \cdot \left(t - 1\right) + -1 \cdot z\right)}\right)\right) + \left(x + b \cdot \left(t - 2\right)\right) \]

      4. distribute-neg-in N/A

         \[\leadsto \color{blue}{\left(\left(\mathsf{neg}\left(a \cdot \left(t - 1\right)\right)\right) + \left(\mathsf{neg}\left(-1 \cdot z\right)\right)\right)} + \left(x + b \cdot \left(t - 2\right)\right) \]

      5. mul-1-neg N/A

         \[\leadsto \left(\left(\mathsf{neg}\left(a \cdot \left(t - 1\right)\right)\right) + \left(\mathsf{neg}\left(\color{blue}{\left(\mathsf{neg}\left(z\right)\right)}\right)\right)\right) + \left(x + b \cdot \left(t - 2\right)\right) \]

      6. remove-double-neg N/A

         \[\leadsto \left(\left(\mathsf{neg}\left(a \cdot \left(t - 1\right)\right)\right) + \color{blue}{z}\right) + \left(x + b \cdot \left(t - 2\right)\right) \]

      7. associate-+l+ N/A

         \[\leadsto \color{blue}{\left(\mathsf{neg}\left(a \cdot \left(t - 1\right)\right)\right) + \left(z + \left(x + b \cdot \left(t - 2\right)\right)\right)} \]

      8. *-commutative N/A

         \[\leadsto \left(\mathsf{neg}\left(\color{blue}{\left(t - 1\right) \cdot a}\right)\right) + \left(z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      9. distribute-lft-neg-in N/A

         \[\leadsto \color{blue}{\left(\mathsf{neg}\left(\left(t - 1\right)\right)\right) \cdot a} + \left(z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      10. mul-1-neg N/A

          \[\leadsto \color{blue}{\left(-1 \cdot \left(t - 1\right)\right)} \cdot a + \left(z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      11. lower-fma.f64 N/A

          \[\leadsto \color{blue}{\mathsf{fma}\left(-1 \cdot \left(t - 1\right), a, z + \left(x + b \cdot \left(t - 2\right)\right)\right)} \]

      12. sub-neg N/A

          \[\leadsto \mathsf{fma}\left(-1 \cdot \color{blue}{\left(t + \left(\mathsf{neg}\left(1\right)\right)\right)}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      13. metadata-eval N/A

          \[\leadsto \mathsf{fma}\left(-1 \cdot \left(t + \color{blue}{-1}\right), a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      14. distribute-lft-in N/A

          \[\leadsto \mathsf{fma}\left(\color{blue}{-1 \cdot t + -1 \cdot -1}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      15. metadata-eval N/A

          \[\leadsto \mathsf{fma}\left(-1 \cdot t + \color{blue}{1}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      16. +-commutative N/A

          \[\leadsto \mathsf{fma}\left(\color{blue}{1 + -1 \cdot t}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      17. neg-mul-1 N/A

          \[\leadsto \mathsf{fma}\left(1 + \color{blue}{\left(\mathsf{neg}\left(t\right)\right)}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      18. sub-neg N/A

          \[\leadsto \mathsf{fma}\left(\color{blue}{1 - t}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      19. lower--.f64 N/A

          \[\leadsto \mathsf{fma}\left(\color{blue}{1 - t}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      20. lower-+.f64 N/A

          \[\leadsto \mathsf{fma}\left(1 - t, a, \color{blue}{z + \left(x + b \cdot \left(t - 2\right)\right)}\right) \]

   5. Applied rewrites 65.8%

      \[\leadsto \color{blue}{\mathsf{fma}\left(1 - t, a, z + \mathsf{fma}\left(t - 2, b, x\right)\right)} \]

   6. Taylor expanded in b around inf

      \[\leadsto \color{blue}{b \cdot \left(\left(t + y\right) - 2\right)} \]
   7. Step-by-step derivation

      1. *-commutative N/A

         \[\leadsto \color{blue}{\left(\left(t + y\right) - 2\right) \cdot b} \]

      2. lower-*.f64 N/A

         \[\leadsto \color{blue}{\left(\left(t + y\right) - 2\right) \cdot b} \]

      3. lower--.f64 N/A

         \[\leadsto \color{blue}{\left(\left(t + y\right) - 2\right)} \cdot b \]

      4. lower-+.f64 75.4

         \[\leadsto \left(\color{blue}{\left(t + y\right)} - 2\right) \cdot b \]

   8. Applied rewrites 75.4%

      \[\leadsto \color{blue}{\left(\left(t + y\right) - 2\right) \cdot b} \]

   if -4.00000000000000017e36 < b < 3.99999999999999993e109

   1. Initial program 98.8%

      \[\left(\left(x - \left(y - 1\right) \cdot z\right) - \left(t - 1\right) \cdot a\right) + \left(\left(y + t\right) - 2\right) \cdot b \]
   2. Add Preprocessing

   3. Taylor expanded in y around 0

      \[\leadsto \color{blue}{\left(x + b \cdot \left(t - 2\right)\right) - \left(-1 \cdot z + a \cdot \left(t - 1\right)\right)} \]
   4. Step-by-step derivation

      1. sub-neg N/A

         \[\leadsto \color{blue}{\left(x + b \cdot \left(t - 2\right)\right) + \left(\mathsf{neg}\left(\left(-1 \cdot z + a \cdot \left(t - 1\right)\right)\right)\right)} \]

      2. +-commutative N/A

         \[\leadsto \color{blue}{\left(\mathsf{neg}\left(\left(-1 \cdot z + a \cdot \left(t - 1\right)\right)\right)\right) + \left(x + b \cdot \left(t - 2\right)\right)} \]

      3. +-commutative N/A

         \[\leadsto \left(\mathsf{neg}\left(\color{blue}{\left(a \cdot \left(t - 1\right) + -1 \cdot z\right)}\right)\right) + \left(x + b \cdot \left(t - 2\right)\right) \]

      4. distribute-neg-in N/A

         \[\leadsto \color{blue}{\left(\left(\mathsf{neg}\left(a \cdot \left(t - 1\right)\right)\right) + \left(\mathsf{neg}\left(-1 \cdot z\right)\right)\right)} + \left(x + b \cdot \left(t - 2\right)\right) \]

      5. mul-1-neg N/A

         \[\leadsto \left(\left(\mathsf{neg}\left(a \cdot \left(t - 1\right)\right)\right) + \left(\mathsf{neg}\left(\color{blue}{\left(\mathsf{neg}\left(z\right)\right)}\right)\right)\right) + \left(x + b \cdot \left(t - 2\right)\right) \]

      6. remove-double-neg N/A

         \[\leadsto \left(\left(\mathsf{neg}\left(a \cdot \left(t - 1\right)\right)\right) + \color{blue}{z}\right) + \left(x + b \cdot \left(t - 2\right)\right) \]

      7. associate-+l+ N/A

         \[\leadsto \color{blue}{\left(\mathsf{neg}\left(a \cdot \left(t - 1\right)\right)\right) + \left(z + \left(x + b \cdot \left(t - 2\right)\right)\right)} \]

      8. *-commutative N/A

         \[\leadsto \left(\mathsf{neg}\left(\color{blue}{\left(t - 1\right) \cdot a}\right)\right) + \left(z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      9. distribute-lft-neg-in N/A

         \[\leadsto \color{blue}{\left(\mathsf{neg}\left(\left(t - 1\right)\right)\right) \cdot a} + \left(z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      10. mul-1-neg N/A

          \[\leadsto \color{blue}{\left(-1 \cdot \left(t - 1\right)\right)} \cdot a + \left(z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      11. lower-fma.f64 N/A

          \[\leadsto \color{blue}{\mathsf{fma}\left(-1 \cdot \left(t - 1\right), a, z + \left(x + b \cdot \left(t - 2\right)\right)\right)} \]

      12. sub-neg N/A

          \[\leadsto \mathsf{fma}\left(-1 \cdot \color{blue}{\left(t + \left(\mathsf{neg}\left(1\right)\right)\right)}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      13. metadata-eval N/A

          \[\leadsto \mathsf{fma}\left(-1 \cdot \left(t + \color{blue}{-1}\right), a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      14. distribute-lft-in N/A

          \[\leadsto \mathsf{fma}\left(\color{blue}{-1 \cdot t + -1 \cdot -1}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      15. metadata-eval N/A

          \[\leadsto \mathsf{fma}\left(-1 \cdot t + \color{blue}{1}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      16. +-commutative N/A

          \[\leadsto \mathsf{fma}\left(\color{blue}{1 + -1 \cdot t}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      17. neg-mul-1 N/A

          \[\leadsto \mathsf{fma}\left(1 + \color{blue}{\left(\mathsf{neg}\left(t\right)\right)}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      18. sub-neg N/A

          \[\leadsto \mathsf{fma}\left(\color{blue}{1 - t}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      19. lower--.f64 N/A

          \[\leadsto \mathsf{fma}\left(\color{blue}{1 - t}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      20. lower-+.f64 N/A

          \[\leadsto \mathsf{fma}\left(1 - t, a, \color{blue}{z + \left(x + b \cdot \left(t - 2\right)\right)}\right) \]

   5. Applied rewrites 81.1%

      \[\leadsto \color{blue}{\mathsf{fma}\left(1 - t, a, z + \mathsf{fma}\left(t - 2, b, x\right)\right)} \]

   6. Taylor expanded in b around 0

      \[\leadsto \mathsf{fma}\left(1 - t, a, x + z\right) \]
   7. Step-by-step derivation

   8. Applied rewrites 76.9%

      \[\leadsto \mathsf{fma}\left(1 - t, a, z + x\right) \]
3. Recombined 2 regimes into one program.

4. Final simplification 76.4%

   \[\leadsto \begin{array}{l} \mathbf{if}\;b \leq -4 \cdot 10^{+36} \lor \neg \left(b \leq 4 \cdot 10^{+109}\right):\\ \;\;\;\;\left(\left(t + y\right) - 2\right) \cdot b\\ \mathbf{else}:\\ \;\;\;\;\mathsf{fma}\left(1 - t, a, z + x\right)\\ \end{array} \]
5. Add Preprocessing

Alternative 12: 61.2% accurate, 1.5× speedup

Math

\[\begin{array}{l} \\ \begin{array}{l} \mathbf{if}\;b \leq -4 \cdot 10^{+36} \lor \neg \left(b \leq 3.8 \cdot 10^{+109}\right):\\ \;\;\;\;\left(\left(t + y\right) - 2\right) \cdot b\\ \mathbf{else}:\\ \;\;\;\;\mathsf{fma}\left(-t, a, z + x\right)\\ \end{array} \end{array} \]

FPCore

(FPCore (x y z t a b)
 :precision binary64
 (if (or (<= b -4e+36) (not (<= b 3.8e+109)))
   (* (- (+ t y) 2.0) b)
   (fma (- t) a (+ z x))))

C

double code(double x, double y, double z, double t, double a, double b) {
	double tmp;
	if ((b <= -4e+36) || !(b <= 3.8e+109)) {
		tmp = ((t + y) - 2.0) * b;
	} else {
		tmp = fma(-t, a, (z + x));
	}
	return tmp;
}

Julia

function code(x, y, z, t, a, b)
	tmp = 0.0
	if ((b <= -4e+36) || !(b <= 3.8e+109))
		tmp = Float64(Float64(Float64(t + y) - 2.0) * b);
	else
		tmp = fma(Float64(-t), a, Float64(z + x));
	end
	return tmp
end

Wolfram

code[x_, y_, z_, t_, a_, b_] := If[Or[LessEqual[b, -4e+36], N[Not[LessEqual[b, 3.8e+109]], $MachinePrecision]], N[(N[(N[(t + y), $MachinePrecision] - 2.0), $MachinePrecision] * b), $MachinePrecision], N[((-t) * a + N[(z + x), $MachinePrecision]), $MachinePrecision]]

Derivation

1. Split input into 2 regimes

2. if b < -4.00000000000000017e36 or 3.80000000000000039e109 < b

   1. Initial program 89.2%

      \[\left(\left(x - \left(y - 1\right) \cdot z\right) - \left(t - 1\right) \cdot a\right) + \left(\left(y + t\right) - 2\right) \cdot b \]
   2. Add Preprocessing

   3. Taylor expanded in y around 0

      \[\leadsto \color{blue}{\left(x + b \cdot \left(t - 2\right)\right) - \left(-1 \cdot z + a \cdot \left(t - 1\right)\right)} \]
   4. Step-by-step derivation

      1. sub-neg N/A

         \[\leadsto \color{blue}{\left(x + b \cdot \left(t - 2\right)\right) + \left(\mathsf{neg}\left(\left(-1 \cdot z + a \cdot \left(t - 1\right)\right)\right)\right)} \]

      2. +-commutative N/A

         \[\leadsto \color{blue}{\left(\mathsf{neg}\left(\left(-1 \cdot z + a \cdot \left(t - 1\right)\right)\right)\right) + \left(x + b \cdot \left(t - 2\right)\right)} \]

      3. +-commutative N/A

         \[\leadsto \left(\mathsf{neg}\left(\color{blue}{\left(a \cdot \left(t - 1\right) + -1 \cdot z\right)}\right)\right) + \left(x + b \cdot \left(t - 2\right)\right) \]

      4. distribute-neg-in N/A

         \[\leadsto \color{blue}{\left(\left(\mathsf{neg}\left(a \cdot \left(t - 1\right)\right)\right) + \left(\mathsf{neg}\left(-1 \cdot z\right)\right)\right)} + \left(x + b \cdot \left(t - 2\right)\right) \]

      5. mul-1-neg N/A

         \[\leadsto \left(\left(\mathsf{neg}\left(a \cdot \left(t - 1\right)\right)\right) + \left(\mathsf{neg}\left(\color{blue}{\left(\mathsf{neg}\left(z\right)\right)}\right)\right)\right) + \left(x + b \cdot \left(t - 2\right)\right) \]

      6. remove-double-neg N/A

         \[\leadsto \left(\left(\mathsf{neg}\left(a \cdot \left(t - 1\right)\right)\right) + \color{blue}{z}\right) + \left(x + b \cdot \left(t - 2\right)\right) \]

      7. associate-+l+ N/A

         \[\leadsto \color{blue}{\left(\mathsf{neg}\left(a \cdot \left(t - 1\right)\right)\right) + \left(z + \left(x + b \cdot \left(t - 2\right)\right)\right)} \]

      8. *-commutative N/A

         \[\leadsto \left(\mathsf{neg}\left(\color{blue}{\left(t - 1\right) \cdot a}\right)\right) + \left(z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      9. distribute-lft-neg-in N/A

         \[\leadsto \color{blue}{\left(\mathsf{neg}\left(\left(t - 1\right)\right)\right) \cdot a} + \left(z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      10. mul-1-neg N/A

          \[\leadsto \color{blue}{\left(-1 \cdot \left(t - 1\right)\right)} \cdot a + \left(z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      11. lower-fma.f64 N/A

          \[\leadsto \color{blue}{\mathsf{fma}\left(-1 \cdot \left(t - 1\right), a, z + \left(x + b \cdot \left(t - 2\right)\right)\right)} \]

      12. sub-neg N/A

          \[\leadsto \mathsf{fma}\left(-1 \cdot \color{blue}{\left(t + \left(\mathsf{neg}\left(1\right)\right)\right)}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      13. metadata-eval N/A

          \[\leadsto \mathsf{fma}\left(-1 \cdot \left(t + \color{blue}{-1}\right), a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      14. distribute-lft-in N/A

          \[\leadsto \mathsf{fma}\left(\color{blue}{-1 \cdot t + -1 \cdot -1}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      15. metadata-eval N/A

          \[\leadsto \mathsf{fma}\left(-1 \cdot t + \color{blue}{1}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      16. +-commutative N/A

          \[\leadsto \mathsf{fma}\left(\color{blue}{1 + -1 \cdot t}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      17. neg-mul-1 N/A

          \[\leadsto \mathsf{fma}\left(1 + \color{blue}{\left(\mathsf{neg}\left(t\right)\right)}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      18. sub-neg N/A

          \[\leadsto \mathsf{fma}\left(\color{blue}{1 - t}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      19. lower--.f64 N/A

          \[\leadsto \mathsf{fma}\left(\color{blue}{1 - t}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      20. lower-+.f64 N/A

          \[\leadsto \mathsf{fma}\left(1 - t, a, \color{blue}{z + \left(x + b \cdot \left(t - 2\right)\right)}\right) \]

   5. Applied rewrites 65.8%

      \[\leadsto \color{blue}{\mathsf{fma}\left(1 - t, a, z + \mathsf{fma}\left(t - 2, b, x\right)\right)} \]

   6. Taylor expanded in b around inf

      \[\leadsto \color{blue}{b \cdot \left(\left(t + y\right) - 2\right)} \]
   7. Step-by-step derivation

      1. *-commutative N/A

         \[\leadsto \color{blue}{\left(\left(t + y\right) - 2\right) \cdot b} \]

      2. lower-*.f64 N/A

         \[\leadsto \color{blue}{\left(\left(t + y\right) - 2\right) \cdot b} \]

      3. lower--.f64 N/A

         \[\leadsto \color{blue}{\left(\left(t + y\right) - 2\right)} \cdot b \]

      4. lower-+.f64 75.4

         \[\leadsto \left(\color{blue}{\left(t + y\right)} - 2\right) \cdot b \]

   8. Applied rewrites 75.4%

      \[\leadsto \color{blue}{\left(\left(t + y\right) - 2\right) \cdot b} \]

   if -4.00000000000000017e36 < b < 3.80000000000000039e109

   1. Initial program 98.8%

      \[\left(\left(x - \left(y - 1\right) \cdot z\right) - \left(t - 1\right) \cdot a\right) + \left(\left(y + t\right) - 2\right) \cdot b \]
   2. Add Preprocessing

   3. Taylor expanded in y around 0

      \[\leadsto \color{blue}{\left(x + b \cdot \left(t - 2\right)\right) - \left(-1 \cdot z + a \cdot \left(t - 1\right)\right)} \]
   4. Step-by-step derivation

      1. sub-neg N/A

         \[\leadsto \color{blue}{\left(x + b \cdot \left(t - 2\right)\right) + \left(\mathsf{neg}\left(\left(-1 \cdot z + a \cdot \left(t - 1\right)\right)\right)\right)} \]

      2. +-commutative N/A

         \[\leadsto \color{blue}{\left(\mathsf{neg}\left(\left(-1 \cdot z + a \cdot \left(t - 1\right)\right)\right)\right) + \left(x + b \cdot \left(t - 2\right)\right)} \]

      3. +-commutative N/A

         \[\leadsto \left(\mathsf{neg}\left(\color{blue}{\left(a \cdot \left(t - 1\right) + -1 \cdot z\right)}\right)\right) + \left(x + b \cdot \left(t - 2\right)\right) \]

      4. distribute-neg-in N/A

         \[\leadsto \color{blue}{\left(\left(\mathsf{neg}\left(a \cdot \left(t - 1\right)\right)\right) + \left(\mathsf{neg}\left(-1 \cdot z\right)\right)\right)} + \left(x + b \cdot \left(t - 2\right)\right) \]

      5. mul-1-neg N/A

         \[\leadsto \left(\left(\mathsf{neg}\left(a \cdot \left(t - 1\right)\right)\right) + \left(\mathsf{neg}\left(\color{blue}{\left(\mathsf{neg}\left(z\right)\right)}\right)\right)\right) + \left(x + b \cdot \left(t - 2\right)\right) \]

      6. remove-double-neg N/A

         \[\leadsto \left(\left(\mathsf{neg}\left(a \cdot \left(t - 1\right)\right)\right) + \color{blue}{z}\right) + \left(x + b \cdot \left(t - 2\right)\right) \]

      7. associate-+l+ N/A

         \[\leadsto \color{blue}{\left(\mathsf{neg}\left(a \cdot \left(t - 1\right)\right)\right) + \left(z + \left(x + b \cdot \left(t - 2\right)\right)\right)} \]

      8. *-commutative N/A

         \[\leadsto \left(\mathsf{neg}\left(\color{blue}{\left(t - 1\right) \cdot a}\right)\right) + \left(z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      9. distribute-lft-neg-in N/A

         \[\leadsto \color{blue}{\left(\mathsf{neg}\left(\left(t - 1\right)\right)\right) \cdot a} + \left(z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      10. mul-1-neg N/A

          \[\leadsto \color{blue}{\left(-1 \cdot \left(t - 1\right)\right)} \cdot a + \left(z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      11. lower-fma.f64 N/A

          \[\leadsto \color{blue}{\mathsf{fma}\left(-1 \cdot \left(t - 1\right), a, z + \left(x + b \cdot \left(t - 2\right)\right)\right)} \]

      12. sub-neg N/A

          \[\leadsto \mathsf{fma}\left(-1 \cdot \color{blue}{\left(t + \left(\mathsf{neg}\left(1\right)\right)\right)}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      13. metadata-eval N/A

          \[\leadsto \mathsf{fma}\left(-1 \cdot \left(t + \color{blue}{-1}\right), a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      14. distribute-lft-in N/A

          \[\leadsto \mathsf{fma}\left(\color{blue}{-1 \cdot t + -1 \cdot -1}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      15. metadata-eval N/A

          \[\leadsto \mathsf{fma}\left(-1 \cdot t + \color{blue}{1}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      16. +-commutative N/A

          \[\leadsto \mathsf{fma}\left(\color{blue}{1 + -1 \cdot t}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      17. neg-mul-1 N/A

          \[\leadsto \mathsf{fma}\left(1 + \color{blue}{\left(\mathsf{neg}\left(t\right)\right)}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      18. sub-neg N/A

          \[\leadsto \mathsf{fma}\left(\color{blue}{1 - t}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      19. lower--.f64 N/A

          \[\leadsto \mathsf{fma}\left(\color{blue}{1 - t}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      20. lower-+.f64 N/A

          \[\leadsto \mathsf{fma}\left(1 - t, a, \color{blue}{z + \left(x + b \cdot \left(t - 2\right)\right)}\right) \]

   5. Applied rewrites 81.1%

      \[\leadsto \color{blue}{\mathsf{fma}\left(1 - t, a, z + \mathsf{fma}\left(t - 2, b, x\right)\right)} \]

   6. Taylor expanded in b around 0

      \[\leadsto \mathsf{fma}\left(1 - t, a, x + z\right) \]
   7. Step-by-step derivation

   8. Applied rewrites 76.9%

      \[\leadsto \mathsf{fma}\left(1 - t, a, z + x\right) \]

   9. Taylor expanded in t around inf

      \[\leadsto \mathsf{fma}\left(-1 \cdot t, a, z + x\right) \]
   10. Step-by-step derivation

   11. Applied rewrites 66.7%

       \[\leadsto \mathsf{fma}\left(-t, a, z + x\right) \]
3. Recombined 2 regimes into one program.

4. Final simplification 69.8%

   \[\leadsto \begin{array}{l} \mathbf{if}\;b \leq -4 \cdot 10^{+36} \lor \neg \left(b \leq 3.8 \cdot 10^{+109}\right):\\ \;\;\;\;\left(\left(t + y\right) - 2\right) \cdot b\\ \mathbf{else}:\\ \;\;\;\;\mathsf{fma}\left(-t, a, z + x\right)\\ \end{array} \]
5. Add Preprocessing

Alternative 13: 26.8% accurate, 1.5× speedup

Math

\[\begin{array}{l} \\ \begin{array}{l} \mathbf{if}\;y \leq -720:\\ \;\;\;\;b \cdot y\\ \mathbf{elif}\;y \leq 7.5 \cdot 10^{-93}:\\ \;\;\;\;1 \cdot z\\ \mathbf{elif}\;y \leq 4.4 \cdot 10^{+60}:\\ \;\;\;\;b \cdot t\\ \mathbf{else}:\\ \;\;\;\;b \cdot y\\ \end{array} \end{array} \]

FPCore

(FPCore (x y z t a b)
 :precision binary64
 (if (<= y -720.0)
   (* b y)
   (if (<= y 7.5e-93) (* 1.0 z) (if (<= y 4.4e+60) (* b t) (* b y)))))

C

double code(double x, double y, double z, double t, double a, double b) {
	double tmp;
	if (y <= -720.0) {
		tmp = b * y;
	} else if (y <= 7.5e-93) {
		tmp = 1.0 * z;
	} else if (y <= 4.4e+60) {
		tmp = b * t;
	} else {
		tmp = b * y;
	}
	return tmp;
}

Fortran

real(8) function code(x, y, z, t, a, b)
    real(8), intent (in) :: x
    real(8), intent (in) :: y
    real(8), intent (in) :: z
    real(8), intent (in) :: t
    real(8), intent (in) :: a
    real(8), intent (in) :: b
    real(8) :: tmp
    if (y <= (-720.0d0)) then
        tmp = b * y
    else if (y <= 7.5d-93) then
        tmp = 1.0d0 * z
    else if (y <= 4.4d+60) then
        tmp = b * t
    else
        tmp = b * y
    end if
    code = tmp
end function

Java

public static double code(double x, double y, double z, double t, double a, double b) {
	double tmp;
	if (y <= -720.0) {
		tmp = b * y;
	} else if (y <= 7.5e-93) {
		tmp = 1.0 * z;
	} else if (y <= 4.4e+60) {
		tmp = b * t;
	} else {
		tmp = b * y;
	}
	return tmp;
}

Python

def code(x, y, z, t, a, b):
	tmp = 0
	if y <= -720.0:
		tmp = b * y
	elif y <= 7.5e-93:
		tmp = 1.0 * z
	elif y <= 4.4e+60:
		tmp = b * t
	else:
		tmp = b * y
	return tmp

Julia

function code(x, y, z, t, a, b)
	tmp = 0.0
	if (y <= -720.0)
		tmp = Float64(b * y);
	elseif (y <= 7.5e-93)
		tmp = Float64(1.0 * z);
	elseif (y <= 4.4e+60)
		tmp = Float64(b * t);
	else
		tmp = Float64(b * y);
	end
	return tmp
end

MATLAB

function tmp_2 = code(x, y, z, t, a, b)
	tmp = 0.0;
	if (y <= -720.0)
		tmp = b * y;
	elseif (y <= 7.5e-93)
		tmp = 1.0 * z;
	elseif (y <= 4.4e+60)
		tmp = b * t;
	else
		tmp = b * y;
	end
	tmp_2 = tmp;
end

Wolfram

code[x_, y_, z_, t_, a_, b_] := If[LessEqual[y, -720.0], N[(b * y), $MachinePrecision], If[LessEqual[y, 7.5e-93], N[(1.0 * z), $MachinePrecision], If[LessEqual[y, 4.4e+60], N[(b * t), $MachinePrecision], N[(b * y), $MachinePrecision]]]]

Derivation

1. Split input into 3 regimes

2. if y < -720 or 4.39999999999999992e60 < y

   1. Initial program 92.5%

      \[\left(\left(x - \left(y - 1\right) \cdot z\right) - \left(t - 1\right) \cdot a\right) + \left(\left(y + t\right) - 2\right) \cdot b \]
   2. Add Preprocessing

   3. Taylor expanded in y around inf

      \[\leadsto \color{blue}{y \cdot \left(b - z\right)} \]
   4. Step-by-step derivation

      1. *-commutative N/A

         \[\leadsto \color{blue}{\left(b - z\right) \cdot y} \]

      2. lower-*.f64 N/A

         \[\leadsto \color{blue}{\left(b - z\right) \cdot y} \]

      3. lower--.f64 63.5

         \[\leadsto \color{blue}{\left(b - z\right)} \cdot y \]

   5. Applied rewrites 63.5%

      \[\leadsto \color{blue}{\left(b - z\right) \cdot y} \]

   6. Taylor expanded in z around 0

      \[\leadsto b \cdot \color{blue}{y} \]
   7. Step-by-step derivation

   8. Applied rewrites 41.3%

      \[\leadsto b \cdot \color{blue}{y} \]

   if -720 < y < 7.50000000000000034e-93

   1. Initial program 97.3%

      \[\left(\left(x - \left(y - 1\right) \cdot z\right) - \left(t - 1\right) \cdot a\right) + \left(\left(y + t\right) - 2\right) \cdot b \]
   2. Add Preprocessing

   3. Taylor expanded in z around inf

      \[\leadsto \color{blue}{z \cdot \left(1 - y\right)} \]
   4. Step-by-step derivation

      1. *-commutative N/A

         \[\leadsto \color{blue}{\left(1 - y\right) \cdot z} \]

      2. sub-neg N/A

         \[\leadsto \color{blue}{\left(1 + \left(\mathsf{neg}\left(y\right)\right)\right)} \cdot z \]

      3. neg-mul-1 N/A

         \[\leadsto \left(1 + \color{blue}{-1 \cdot y}\right) \cdot z \]

      4. +-commutative N/A

         \[\leadsto \color{blue}{\left(-1 \cdot y + 1\right)} \cdot z \]

      5. metadata-eval N/A

         \[\leadsto \left(-1 \cdot y + \color{blue}{-1 \cdot -1}\right) \cdot z \]

      6. distribute-lft-in N/A

         \[\leadsto \color{blue}{\left(-1 \cdot \left(y + -1\right)\right)} \cdot z \]

      7. metadata-eval N/A

         \[\leadsto \left(-1 \cdot \left(y + \color{blue}{\left(\mathsf{neg}\left(1\right)\right)}\right)\right) \cdot z \]

      8. sub-neg N/A

         \[\leadsto \left(-1 \cdot \color{blue}{\left(y - 1\right)}\right) \cdot z \]

      9. lower-*.f64 N/A

         \[\leadsto \color{blue}{\left(-1 \cdot \left(y - 1\right)\right) \cdot z} \]

      10. sub-neg N/A

          \[\leadsto \left(-1 \cdot \color{blue}{\left(y + \left(\mathsf{neg}\left(1\right)\right)\right)}\right) \cdot z \]

      11. metadata-eval N/A

          \[\leadsto \left(-1 \cdot \left(y + \color{blue}{-1}\right)\right) \cdot z \]

      12. distribute-lft-in N/A

          \[\leadsto \color{blue}{\left(-1 \cdot y + -1 \cdot -1\right)} \cdot z \]

      13. metadata-eval N/A

          \[\leadsto \left(-1 \cdot y + \color{blue}{1}\right) \cdot z \]

      14. +-commutative N/A

          \[\leadsto \color{blue}{\left(1 + -1 \cdot y\right)} \cdot z \]

      15. neg-mul-1 N/A

          \[\leadsto \left(1 + \color{blue}{\left(\mathsf{neg}\left(y\right)\right)}\right) \cdot z \]

      16. sub-neg N/A

          \[\leadsto \color{blue}{\left(1 - y\right)} \cdot z \]

      17. lower--.f64 21.8

          \[\leadsto \color{blue}{\left(1 - y\right)} \cdot z \]

   5. Applied rewrites 21.8%

      \[\leadsto \color{blue}{\left(1 - y\right) \cdot z} \]

   6. Taylor expanded in y around 0

      \[\leadsto 1 \cdot z \]
   7. Step-by-step derivation

   8. Applied rewrites 21.8%

      \[\leadsto 1 \cdot z \]

   if 7.50000000000000034e-93 < y < 4.39999999999999992e60

   1. Initial program 97.2%

      \[\left(\left(x - \left(y - 1\right) \cdot z\right) - \left(t - 1\right) \cdot a\right) + \left(\left(y + t\right) - 2\right) \cdot b \]
   2. Add Preprocessing

   3. Taylor expanded in y around 0

      \[\leadsto \color{blue}{\left(x + b \cdot \left(t - 2\right)\right) - \left(-1 \cdot z + a \cdot \left(t - 1\right)\right)} \]
   4. Step-by-step derivation

      1. sub-neg N/A

         \[\leadsto \color{blue}{\left(x + b \cdot \left(t - 2\right)\right) + \left(\mathsf{neg}\left(\left(-1 \cdot z + a \cdot \left(t - 1\right)\right)\right)\right)} \]

      2. +-commutative N/A

         \[\leadsto \color{blue}{\left(\mathsf{neg}\left(\left(-1 \cdot z + a \cdot \left(t - 1\right)\right)\right)\right) + \left(x + b \cdot \left(t - 2\right)\right)} \]

      3. +-commutative N/A

         \[\leadsto \left(\mathsf{neg}\left(\color{blue}{\left(a \cdot \left(t - 1\right) + -1 \cdot z\right)}\right)\right) + \left(x + b \cdot \left(t - 2\right)\right) \]

      4. distribute-neg-in N/A

         \[\leadsto \color{blue}{\left(\left(\mathsf{neg}\left(a \cdot \left(t - 1\right)\right)\right) + \left(\mathsf{neg}\left(-1 \cdot z\right)\right)\right)} + \left(x + b \cdot \left(t - 2\right)\right) \]

      5. mul-1-neg N/A

         \[\leadsto \left(\left(\mathsf{neg}\left(a \cdot \left(t - 1\right)\right)\right) + \left(\mathsf{neg}\left(\color{blue}{\left(\mathsf{neg}\left(z\right)\right)}\right)\right)\right) + \left(x + b \cdot \left(t - 2\right)\right) \]

      6. remove-double-neg N/A

         \[\leadsto \left(\left(\mathsf{neg}\left(a \cdot \left(t - 1\right)\right)\right) + \color{blue}{z}\right) + \left(x + b \cdot \left(t - 2\right)\right) \]

      7. associate-+l+ N/A

         \[\leadsto \color{blue}{\left(\mathsf{neg}\left(a \cdot \left(t - 1\right)\right)\right) + \left(z + \left(x + b \cdot \left(t - 2\right)\right)\right)} \]

      8. *-commutative N/A

         \[\leadsto \left(\mathsf{neg}\left(\color{blue}{\left(t - 1\right) \cdot a}\right)\right) + \left(z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      9. distribute-lft-neg-in N/A

         \[\leadsto \color{blue}{\left(\mathsf{neg}\left(\left(t - 1\right)\right)\right) \cdot a} + \left(z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      10. mul-1-neg N/A

          \[\leadsto \color{blue}{\left(-1 \cdot \left(t - 1\right)\right)} \cdot a + \left(z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      11. lower-fma.f64 N/A

          \[\leadsto \color{blue}{\mathsf{fma}\left(-1 \cdot \left(t - 1\right), a, z + \left(x + b \cdot \left(t - 2\right)\right)\right)} \]

      12. sub-neg N/A

          \[\leadsto \mathsf{fma}\left(-1 \cdot \color{blue}{\left(t + \left(\mathsf{neg}\left(1\right)\right)\right)}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      13. metadata-eval N/A

          \[\leadsto \mathsf{fma}\left(-1 \cdot \left(t + \color{blue}{-1}\right), a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      14. distribute-lft-in N/A

          \[\leadsto \mathsf{fma}\left(\color{blue}{-1 \cdot t + -1 \cdot -1}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      15. metadata-eval N/A

          \[\leadsto \mathsf{fma}\left(-1 \cdot t + \color{blue}{1}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      16. +-commutative N/A

          \[\leadsto \mathsf{fma}\left(\color{blue}{1 + -1 \cdot t}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      17. neg-mul-1 N/A

          \[\leadsto \mathsf{fma}\left(1 + \color{blue}{\left(\mathsf{neg}\left(t\right)\right)}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      18. sub-neg N/A

          \[\leadsto \mathsf{fma}\left(\color{blue}{1 - t}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      19. lower--.f64 N/A

          \[\leadsto \mathsf{fma}\left(\color{blue}{1 - t}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      20. lower-+.f64 N/A

          \[\leadsto \mathsf{fma}\left(1 - t, a, \color{blue}{z + \left(x + b \cdot \left(t - 2\right)\right)}\right) \]

   5. Applied rewrites 94.2%

      \[\leadsto \color{blue}{\mathsf{fma}\left(1 - t, a, z + \mathsf{fma}\left(t - 2, b, x\right)\right)} \]

   6. Taylor expanded in b around inf

      \[\leadsto \color{blue}{b \cdot \left(\left(t + y\right) - 2\right)} \]
   7. Step-by-step derivation

      1. *-commutative N/A

         \[\leadsto \color{blue}{\left(\left(t + y\right) - 2\right) \cdot b} \]

      2. lower-*.f64 N/A

         \[\leadsto \color{blue}{\left(\left(t + y\right) - 2\right) \cdot b} \]

      3. lower--.f64 N/A

         \[\leadsto \color{blue}{\left(\left(t + y\right) - 2\right)} \cdot b \]

      4. lower-+.f64 44.0

         \[\leadsto \left(\color{blue}{\left(t + y\right)} - 2\right) \cdot b \]

   8. Applied rewrites 44.0%

      \[\leadsto \color{blue}{\left(\left(t + y\right) - 2\right) \cdot b} \]

   9. Taylor expanded in t around inf

      \[\leadsto b \cdot \color{blue}{t} \]
   10. Step-by-step derivation

   11. Applied rewrites 35.4%

       \[\leadsto b \cdot \color{blue}{t} \]
3. Recombined 3 regimes into one program.
4. Add Preprocessing

Alternative 14: 56.4% accurate, 1.7× speedup

Math

\[\begin{array}{l} \\ \begin{array}{l} \mathbf{if}\;t \leq -2.6 \cdot 10^{+41} \lor \neg \left(t \leq 7600000000\right):\\ \;\;\;\;\left(b - a\right) \cdot t\\ \mathbf{else}:\\ \;\;\;\;\mathsf{fma}\left(1, a, z + x\right)\\ \end{array} \end{array} \]

FPCore

(FPCore (x y z t a b)
 :precision binary64
 (if (or (<= t -2.6e+41) (not (<= t 7600000000.0)))
   (* (- b a) t)
   (fma 1.0 a (+ z x))))

C

double code(double x, double y, double z, double t, double a, double b) {
	double tmp;
	if ((t <= -2.6e+41) || !(t <= 7600000000.0)) {
		tmp = (b - a) * t;
	} else {
		tmp = fma(1.0, a, (z + x));
	}
	return tmp;
}

Julia

function code(x, y, z, t, a, b)
	tmp = 0.0
	if ((t <= -2.6e+41) || !(t <= 7600000000.0))
		tmp = Float64(Float64(b - a) * t);
	else
		tmp = fma(1.0, a, Float64(z + x));
	end
	return tmp
end

Wolfram

code[x_, y_, z_, t_, a_, b_] := If[Or[LessEqual[t, -2.6e+41], N[Not[LessEqual[t, 7600000000.0]], $MachinePrecision]], N[(N[(b - a), $MachinePrecision] * t), $MachinePrecision], N[(1.0 * a + N[(z + x), $MachinePrecision]), $MachinePrecision]]

Derivation

1. Split input into 2 regimes

2. if t < -2.6000000000000001e41 or 7.6e9 < t

   1. Initial program 92.3%

      \[\left(\left(x - \left(y - 1\right) \cdot z\right) - \left(t - 1\right) \cdot a\right) + \left(\left(y + t\right) - 2\right) \cdot b \]
   2. Add Preprocessing

   3. Taylor expanded in t around inf

      \[\leadsto \color{blue}{t \cdot \left(b - a\right)} \]
   4. Step-by-step derivation

      1. *-commutative N/A

         \[\leadsto \color{blue}{\left(b - a\right) \cdot t} \]

      2. lower-*.f64 N/A

         \[\leadsto \color{blue}{\left(b - a\right) \cdot t} \]

      3. lower--.f64 76.2

         \[\leadsto \color{blue}{\left(b - a\right)} \cdot t \]

   5. Applied rewrites 76.2%

      \[\leadsto \color{blue}{\left(b - a\right) \cdot t} \]

   if -2.6000000000000001e41 < t < 7.6e9

   1. Initial program 97.8%

      \[\left(\left(x - \left(y - 1\right) \cdot z\right) - \left(t - 1\right) \cdot a\right) + \left(\left(y + t\right) - 2\right) \cdot b \]
   2. Add Preprocessing

   3. Taylor expanded in y around 0

      \[\leadsto \color{blue}{\left(x + b \cdot \left(t - 2\right)\right) - \left(-1 \cdot z + a \cdot \left(t - 1\right)\right)} \]
   4. Step-by-step derivation

      1. sub-neg N/A

         \[\leadsto \color{blue}{\left(x + b \cdot \left(t - 2\right)\right) + \left(\mathsf{neg}\left(\left(-1 \cdot z + a \cdot \left(t - 1\right)\right)\right)\right)} \]

      2. +-commutative N/A

         \[\leadsto \color{blue}{\left(\mathsf{neg}\left(\left(-1 \cdot z + a \cdot \left(t - 1\right)\right)\right)\right) + \left(x + b \cdot \left(t - 2\right)\right)} \]

      3. +-commutative N/A

         \[\leadsto \left(\mathsf{neg}\left(\color{blue}{\left(a \cdot \left(t - 1\right) + -1 \cdot z\right)}\right)\right) + \left(x + b \cdot \left(t - 2\right)\right) \]

      4. distribute-neg-in N/A

         \[\leadsto \color{blue}{\left(\left(\mathsf{neg}\left(a \cdot \left(t - 1\right)\right)\right) + \left(\mathsf{neg}\left(-1 \cdot z\right)\right)\right)} + \left(x + b \cdot \left(t - 2\right)\right) \]

      5. mul-1-neg N/A

         \[\leadsto \left(\left(\mathsf{neg}\left(a \cdot \left(t - 1\right)\right)\right) + \left(\mathsf{neg}\left(\color{blue}{\left(\mathsf{neg}\left(z\right)\right)}\right)\right)\right) + \left(x + b \cdot \left(t - 2\right)\right) \]

      6. remove-double-neg N/A

         \[\leadsto \left(\left(\mathsf{neg}\left(a \cdot \left(t - 1\right)\right)\right) + \color{blue}{z}\right) + \left(x + b \cdot \left(t - 2\right)\right) \]

      7. associate-+l+ N/A

         \[\leadsto \color{blue}{\left(\mathsf{neg}\left(a \cdot \left(t - 1\right)\right)\right) + \left(z + \left(x + b \cdot \left(t - 2\right)\right)\right)} \]

      8. *-commutative N/A

         \[\leadsto \left(\mathsf{neg}\left(\color{blue}{\left(t - 1\right) \cdot a}\right)\right) + \left(z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      9. distribute-lft-neg-in N/A

         \[\leadsto \color{blue}{\left(\mathsf{neg}\left(\left(t - 1\right)\right)\right) \cdot a} + \left(z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      10. mul-1-neg N/A

          \[\leadsto \color{blue}{\left(-1 \cdot \left(t - 1\right)\right)} \cdot a + \left(z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      11. lower-fma.f64 N/A

          \[\leadsto \color{blue}{\mathsf{fma}\left(-1 \cdot \left(t - 1\right), a, z + \left(x + b \cdot \left(t - 2\right)\right)\right)} \]

      12. sub-neg N/A

          \[\leadsto \mathsf{fma}\left(-1 \cdot \color{blue}{\left(t + \left(\mathsf{neg}\left(1\right)\right)\right)}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      13. metadata-eval N/A

          \[\leadsto \mathsf{fma}\left(-1 \cdot \left(t + \color{blue}{-1}\right), a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      14. distribute-lft-in N/A

          \[\leadsto \mathsf{fma}\left(\color{blue}{-1 \cdot t + -1 \cdot -1}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      15. metadata-eval N/A

          \[\leadsto \mathsf{fma}\left(-1 \cdot t + \color{blue}{1}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      16. +-commutative N/A

          \[\leadsto \mathsf{fma}\left(\color{blue}{1 + -1 \cdot t}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      17. neg-mul-1 N/A

          \[\leadsto \mathsf{fma}\left(1 + \color{blue}{\left(\mathsf{neg}\left(t\right)\right)}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      18. sub-neg N/A

          \[\leadsto \mathsf{fma}\left(\color{blue}{1 - t}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      19. lower--.f64 N/A

          \[\leadsto \mathsf{fma}\left(\color{blue}{1 - t}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      20. lower-+.f64 N/A

          \[\leadsto \mathsf{fma}\left(1 - t, a, \color{blue}{z + \left(x + b \cdot \left(t - 2\right)\right)}\right) \]

   5. Applied rewrites 65.8%

      \[\leadsto \color{blue}{\mathsf{fma}\left(1 - t, a, z + \mathsf{fma}\left(t - 2, b, x\right)\right)} \]

   6. Taylor expanded in b around 0

      \[\leadsto \mathsf{fma}\left(1 - t, a, x + z\right) \]
   7. Step-by-step derivation

   8. Applied rewrites 56.7%

      \[\leadsto \mathsf{fma}\left(1 - t, a, z + x\right) \]

   9. Taylor expanded in t around 0

      \[\leadsto \mathsf{fma}\left(1, a, z + x\right) \]
   10. Step-by-step derivation

   11. Applied rewrites 54.1%

       \[\leadsto \mathsf{fma}\left(1, a, z + x\right) \]
3. Recombined 2 regimes into one program.

4. Final simplification 64.2%

   \[\leadsto \begin{array}{l} \mathbf{if}\;t \leq -2.6 \cdot 10^{+41} \lor \neg \left(t \leq 7600000000\right):\\ \;\;\;\;\left(b - a\right) \cdot t\\ \mathbf{else}:\\ \;\;\;\;\mathsf{fma}\left(1, a, z + x\right)\\ \end{array} \]
5. Add Preprocessing

Alternative 15: 50.3% accurate, 1.8× speedup

Math

\[\begin{array}{l} \\ \begin{array}{l} \mathbf{if}\;t \leq -8000000 \lor \neg \left(t \leq 4.3 \cdot 10^{+42}\right):\\ \;\;\;\;\left(b - a\right) \cdot t\\ \mathbf{else}:\\ \;\;\;\;\left(b - z\right) \cdot y\\ \end{array} \end{array} \]

FPCore

(FPCore (x y z t a b)
 :precision binary64
 (if (or (<= t -8000000.0) (not (<= t 4.3e+42))) (* (- b a) t) (* (- b z) y)))

C

double code(double x, double y, double z, double t, double a, double b) {
	double tmp;
	if ((t <= -8000000.0) || !(t <= 4.3e+42)) {
		tmp = (b - a) * t;
	} else {
		tmp = (b - z) * y;
	}
	return tmp;
}

Fortran

real(8) function code(x, y, z, t, a, b)
    real(8), intent (in) :: x
    real(8), intent (in) :: y
    real(8), intent (in) :: z
    real(8), intent (in) :: t
    real(8), intent (in) :: a
    real(8), intent (in) :: b
    real(8) :: tmp
    if ((t <= (-8000000.0d0)) .or. (.not. (t <= 4.3d+42))) then
        tmp = (b - a) * t
    else
        tmp = (b - z) * y
    end if
    code = tmp
end function

Java

public static double code(double x, double y, double z, double t, double a, double b) {
	double tmp;
	if ((t <= -8000000.0) || !(t <= 4.3e+42)) {
		tmp = (b - a) * t;
	} else {
		tmp = (b - z) * y;
	}
	return tmp;
}

Python

def code(x, y, z, t, a, b):
	tmp = 0
	if (t <= -8000000.0) or not (t <= 4.3e+42):
		tmp = (b - a) * t
	else:
		tmp = (b - z) * y
	return tmp

Julia

function code(x, y, z, t, a, b)
	tmp = 0.0
	if ((t <= -8000000.0) || !(t <= 4.3e+42))
		tmp = Float64(Float64(b - a) * t);
	else
		tmp = Float64(Float64(b - z) * y);
	end
	return tmp
end

MATLAB

function tmp_2 = code(x, y, z, t, a, b)
	tmp = 0.0;
	if ((t <= -8000000.0) || ~((t <= 4.3e+42)))
		tmp = (b - a) * t;
	else
		tmp = (b - z) * y;
	end
	tmp_2 = tmp;
end

Wolfram

code[x_, y_, z_, t_, a_, b_] := If[Or[LessEqual[t, -8000000.0], N[Not[LessEqual[t, 4.3e+42]], $MachinePrecision]], N[(N[(b - a), $MachinePrecision] * t), $MachinePrecision], N[(N[(b - z), $MachinePrecision] * y), $MachinePrecision]]

Derivation

1. Split input into 2 regimes

2. if t < -8e6 or 4.2999999999999998e42 < t

   1. Initial program 92.8%

      \[\left(\left(x - \left(y - 1\right) \cdot z\right) - \left(t - 1\right) \cdot a\right) + \left(\left(y + t\right) - 2\right) \cdot b \]
   2. Add Preprocessing

   3. Taylor expanded in t around inf

      \[\leadsto \color{blue}{t \cdot \left(b - a\right)} \]
   4. Step-by-step derivation

      1. *-commutative N/A

         \[\leadsto \color{blue}{\left(b - a\right) \cdot t} \]

      2. lower-*.f64 N/A

         \[\leadsto \color{blue}{\left(b - a\right) \cdot t} \]

      3. lower--.f64 74.0

         \[\leadsto \color{blue}{\left(b - a\right)} \cdot t \]

   5. Applied rewrites 74.0%

      \[\leadsto \color{blue}{\left(b - a\right) \cdot t} \]

   if -8e6 < t < 4.2999999999999998e42

   1. Initial program 97.7%

      \[\left(\left(x - \left(y - 1\right) \cdot z\right) - \left(t - 1\right) \cdot a\right) + \left(\left(y + t\right) - 2\right) \cdot b \]
   2. Add Preprocessing

   3. Taylor expanded in y around inf

      \[\leadsto \color{blue}{y \cdot \left(b - z\right)} \]
   4. Step-by-step derivation

      1. *-commutative N/A

         \[\leadsto \color{blue}{\left(b - z\right) \cdot y} \]

      2. lower-*.f64 N/A

         \[\leadsto \color{blue}{\left(b - z\right) \cdot y} \]

      3. lower--.f64 39.3

         \[\leadsto \color{blue}{\left(b - z\right)} \cdot y \]

   5. Applied rewrites 39.3%

      \[\leadsto \color{blue}{\left(b - z\right) \cdot y} \]
3. Recombined 2 regimes into one program.

4. Final simplification 56.4%

   \[\leadsto \begin{array}{l} \mathbf{if}\;t \leq -8000000 \lor \neg \left(t \leq 4.3 \cdot 10^{+42}\right):\\ \;\;\;\;\left(b - a\right) \cdot t\\ \mathbf{else}:\\ \;\;\;\;\left(b - z\right) \cdot y\\ \end{array} \]
5. Add Preprocessing

Alternative 16: 42.1% accurate, 1.8× speedup

Math

\[\begin{array}{l} \\ \begin{array}{l} \mathbf{if}\;a \leq -34000000 \lor \neg \left(a \leq 2.9 \cdot 10^{-5}\right):\\ \;\;\;\;\left(1 - t\right) \cdot a\\ \mathbf{else}:\\ \;\;\;\;\left(1 - y\right) \cdot z\\ \end{array} \end{array} \]

FPCore

(FPCore (x y z t a b)
 :precision binary64
 (if (or (<= a -34000000.0) (not (<= a 2.9e-5)))
   (* (- 1.0 t) a)
   (* (- 1.0 y) z)))

C

double code(double x, double y, double z, double t, double a, double b) {
	double tmp;
	if ((a <= -34000000.0) || !(a <= 2.9e-5)) {
		tmp = (1.0 - t) * a;
	} else {
		tmp = (1.0 - y) * z;
	}
	return tmp;
}

Fortran

real(8) function code(x, y, z, t, a, b)
    real(8), intent (in) :: x
    real(8), intent (in) :: y
    real(8), intent (in) :: z
    real(8), intent (in) :: t
    real(8), intent (in) :: a
    real(8), intent (in) :: b
    real(8) :: tmp
    if ((a <= (-34000000.0d0)) .or. (.not. (a <= 2.9d-5))) then
        tmp = (1.0d0 - t) * a
    else
        tmp = (1.0d0 - y) * z
    end if
    code = tmp
end function

Java

public static double code(double x, double y, double z, double t, double a, double b) {
	double tmp;
	if ((a <= -34000000.0) || !(a <= 2.9e-5)) {
		tmp = (1.0 - t) * a;
	} else {
		tmp = (1.0 - y) * z;
	}
	return tmp;
}

Python

def code(x, y, z, t, a, b):
	tmp = 0
	if (a <= -34000000.0) or not (a <= 2.9e-5):
		tmp = (1.0 - t) * a
	else:
		tmp = (1.0 - y) * z
	return tmp

Julia

function code(x, y, z, t, a, b)
	tmp = 0.0
	if ((a <= -34000000.0) || !(a <= 2.9e-5))
		tmp = Float64(Float64(1.0 - t) * a);
	else
		tmp = Float64(Float64(1.0 - y) * z);
	end
	return tmp
end

MATLAB

function tmp_2 = code(x, y, z, t, a, b)
	tmp = 0.0;
	if ((a <= -34000000.0) || ~((a <= 2.9e-5)))
		tmp = (1.0 - t) * a;
	else
		tmp = (1.0 - y) * z;
	end
	tmp_2 = tmp;
end

Wolfram

code[x_, y_, z_, t_, a_, b_] := If[Or[LessEqual[a, -34000000.0], N[Not[LessEqual[a, 2.9e-5]], $MachinePrecision]], N[(N[(1.0 - t), $MachinePrecision] * a), $MachinePrecision], N[(N[(1.0 - y), $MachinePrecision] * z), $MachinePrecision]]

Derivation

1. Split input into 2 regimes

2. if a < -3.4e7 or 2.9e-5 < a

   1. Initial program 91.7%

      \[\left(\left(x - \left(y - 1\right) \cdot z\right) - \left(t - 1\right) \cdot a\right) + \left(\left(y + t\right) - 2\right) \cdot b \]
   2. Add Preprocessing

   3. Taylor expanded in a around inf

      \[\leadsto \color{blue}{a \cdot \left(1 - t\right)} \]
   4. Step-by-step derivation

      1. *-commutative N/A

         \[\leadsto \color{blue}{\left(1 - t\right) \cdot a} \]

      2. sub-neg N/A

         \[\leadsto \color{blue}{\left(1 + \left(\mathsf{neg}\left(t\right)\right)\right)} \cdot a \]

      3. neg-mul-1 N/A

         \[\leadsto \left(1 + \color{blue}{-1 \cdot t}\right) \cdot a \]

      4. +-commutative N/A

         \[\leadsto \color{blue}{\left(-1 \cdot t + 1\right)} \cdot a \]

      5. metadata-eval N/A

         \[\leadsto \left(-1 \cdot t + \color{blue}{-1 \cdot -1}\right) \cdot a \]

      6. distribute-lft-in N/A

         \[\leadsto \color{blue}{\left(-1 \cdot \left(t + -1\right)\right)} \cdot a \]

      7. metadata-eval N/A

         \[\leadsto \left(-1 \cdot \left(t + \color{blue}{\left(\mathsf{neg}\left(1\right)\right)}\right)\right) \cdot a \]

      8. sub-neg N/A

         \[\leadsto \left(-1 \cdot \color{blue}{\left(t - 1\right)}\right) \cdot a \]

      9. lower-*.f64 N/A

         \[\leadsto \color{blue}{\left(-1 \cdot \left(t - 1\right)\right) \cdot a} \]

      10. sub-neg N/A

          \[\leadsto \left(-1 \cdot \color{blue}{\left(t + \left(\mathsf{neg}\left(1\right)\right)\right)}\right) \cdot a \]

      11. metadata-eval N/A

          \[\leadsto \left(-1 \cdot \left(t + \color{blue}{-1}\right)\right) \cdot a \]

      12. distribute-lft-in N/A

          \[\leadsto \color{blue}{\left(-1 \cdot t + -1 \cdot -1\right)} \cdot a \]

      13. metadata-eval N/A

          \[\leadsto \left(-1 \cdot t + \color{blue}{1}\right) \cdot a \]

      14. +-commutative N/A

          \[\leadsto \color{blue}{\left(1 + -1 \cdot t\right)} \cdot a \]

      15. neg-mul-1 N/A

          \[\leadsto \left(1 + \color{blue}{\left(\mathsf{neg}\left(t\right)\right)}\right) \cdot a \]

      16. sub-neg N/A

          \[\leadsto \color{blue}{\left(1 - t\right)} \cdot a \]

      17. lower--.f64 58.7

          \[\leadsto \color{blue}{\left(1 - t\right)} \cdot a \]

   5. Applied rewrites 58.7%

      \[\leadsto \color{blue}{\left(1 - t\right) \cdot a} \]

   if -3.4e7 < a < 2.9e-5

   1. Initial program 99.2%

      \[\left(\left(x - \left(y - 1\right) \cdot z\right) - \left(t - 1\right) \cdot a\right) + \left(\left(y + t\right) - 2\right) \cdot b \]
   2. Add Preprocessing

   3. Taylor expanded in z around inf

      \[\leadsto \color{blue}{z \cdot \left(1 - y\right)} \]
   4. Step-by-step derivation

      1. *-commutative N/A

         \[\leadsto \color{blue}{\left(1 - y\right) \cdot z} \]

      2. sub-neg N/A

         \[\leadsto \color{blue}{\left(1 + \left(\mathsf{neg}\left(y\right)\right)\right)} \cdot z \]

      3. neg-mul-1 N/A

         \[\leadsto \left(1 + \color{blue}{-1 \cdot y}\right) \cdot z \]

      4. +-commutative N/A

         \[\leadsto \color{blue}{\left(-1 \cdot y + 1\right)} \cdot z \]

      5. metadata-eval N/A

         \[\leadsto \left(-1 \cdot y + \color{blue}{-1 \cdot -1}\right) \cdot z \]

      6. distribute-lft-in N/A

         \[\leadsto \color{blue}{\left(-1 \cdot \left(y + -1\right)\right)} \cdot z \]

      7. metadata-eval N/A

         \[\leadsto \left(-1 \cdot \left(y + \color{blue}{\left(\mathsf{neg}\left(1\right)\right)}\right)\right) \cdot z \]

      8. sub-neg N/A

         \[\leadsto \left(-1 \cdot \color{blue}{\left(y - 1\right)}\right) \cdot z \]

      9. lower-*.f64 N/A

         \[\leadsto \color{blue}{\left(-1 \cdot \left(y - 1\right)\right) \cdot z} \]

      10. sub-neg N/A

          \[\leadsto \left(-1 \cdot \color{blue}{\left(y + \left(\mathsf{neg}\left(1\right)\right)\right)}\right) \cdot z \]

      11. metadata-eval N/A

          \[\leadsto \left(-1 \cdot \left(y + \color{blue}{-1}\right)\right) \cdot z \]

      12. distribute-lft-in N/A

          \[\leadsto \color{blue}{\left(-1 \cdot y + -1 \cdot -1\right)} \cdot z \]

      13. metadata-eval N/A

          \[\leadsto \left(-1 \cdot y + \color{blue}{1}\right) \cdot z \]

      14. +-commutative N/A

          \[\leadsto \color{blue}{\left(1 + -1 \cdot y\right)} \cdot z \]

      15. neg-mul-1 N/A

          \[\leadsto \left(1 + \color{blue}{\left(\mathsf{neg}\left(y\right)\right)}\right) \cdot z \]

      16. sub-neg N/A

          \[\leadsto \color{blue}{\left(1 - y\right)} \cdot z \]

      17. lower--.f64 37.8

          \[\leadsto \color{blue}{\left(1 - y\right)} \cdot z \]

   5. Applied rewrites 37.8%

      \[\leadsto \color{blue}{\left(1 - y\right) \cdot z} \]
3. Recombined 2 regimes into one program.

4. Final simplification 48.6%

   \[\leadsto \begin{array}{l} \mathbf{if}\;a \leq -34000000 \lor \neg \left(a \leq 2.9 \cdot 10^{-5}\right):\\ \;\;\;\;\left(1 - t\right) \cdot a\\ \mathbf{else}:\\ \;\;\;\;\left(1 - y\right) \cdot z\\ \end{array} \]
5. Add Preprocessing

Alternative 17: 27.0% accurate, 1.8× speedup

Math

\[\begin{array}{l} \\ \begin{array}{l} \mathbf{if}\;b \leq -9.4 \cdot 10^{+26}:\\ \;\;\;\;b \cdot y\\ \mathbf{elif}\;b \leq 2.85 \cdot 10^{+110}:\\ \;\;\;\;\left(-a\right) \cdot t\\ \mathbf{else}:\\ \;\;\;\;b \cdot t\\ \end{array} \end{array} \]

FPCore

(FPCore (x y z t a b)
 :precision binary64
 (if (<= b -9.4e+26) (* b y) (if (<= b 2.85e+110) (* (- a) t) (* b t))))

C

double code(double x, double y, double z, double t, double a, double b) {
	double tmp;
	if (b <= -9.4e+26) {
		tmp = b * y;
	} else if (b <= 2.85e+110) {
		tmp = -a * t;
	} else {
		tmp = b * t;
	}
	return tmp;
}

Fortran

real(8) function code(x, y, z, t, a, b)
    real(8), intent (in) :: x
    real(8), intent (in) :: y
    real(8), intent (in) :: z
    real(8), intent (in) :: t
    real(8), intent (in) :: a
    real(8), intent (in) :: b
    real(8) :: tmp
    if (b <= (-9.4d+26)) then
        tmp = b * y
    else if (b <= 2.85d+110) then
        tmp = -a * t
    else
        tmp = b * t
    end if
    code = tmp
end function

Java

public static double code(double x, double y, double z, double t, double a, double b) {
	double tmp;
	if (b <= -9.4e+26) {
		tmp = b * y;
	} else if (b <= 2.85e+110) {
		tmp = -a * t;
	} else {
		tmp = b * t;
	}
	return tmp;
}

Python

def code(x, y, z, t, a, b):
	tmp = 0
	if b <= -9.4e+26:
		tmp = b * y
	elif b <= 2.85e+110:
		tmp = -a * t
	else:
		tmp = b * t
	return tmp

Julia

function code(x, y, z, t, a, b)
	tmp = 0.0
	if (b <= -9.4e+26)
		tmp = Float64(b * y);
	elseif (b <= 2.85e+110)
		tmp = Float64(Float64(-a) * t);
	else
		tmp = Float64(b * t);
	end
	return tmp
end

MATLAB

function tmp_2 = code(x, y, z, t, a, b)
	tmp = 0.0;
	if (b <= -9.4e+26)
		tmp = b * y;
	elseif (b <= 2.85e+110)
		tmp = -a * t;
	else
		tmp = b * t;
	end
	tmp_2 = tmp;
end

Wolfram

code[x_, y_, z_, t_, a_, b_] := If[LessEqual[b, -9.4e+26], N[(b * y), $MachinePrecision], If[LessEqual[b, 2.85e+110], N[((-a) * t), $MachinePrecision], N[(b * t), $MachinePrecision]]]

Derivation

1. Split input into 3 regimes

2. if b < -9.3999999999999995e26

   1. Initial program 91.5%

      \[\left(\left(x - \left(y - 1\right) \cdot z\right) - \left(t - 1\right) \cdot a\right) + \left(\left(y + t\right) - 2\right) \cdot b \]
   2. Add Preprocessing

   3. Taylor expanded in y around inf

      \[\leadsto \color{blue}{y \cdot \left(b - z\right)} \]
   4. Step-by-step derivation

      1. *-commutative N/A

         \[\leadsto \color{blue}{\left(b - z\right) \cdot y} \]

      2. lower-*.f64 N/A

         \[\leadsto \color{blue}{\left(b - z\right) \cdot y} \]

      3. lower--.f64 48.5

         \[\leadsto \color{blue}{\left(b - z\right)} \cdot y \]

   5. Applied rewrites 48.5%

      \[\leadsto \color{blue}{\left(b - z\right) \cdot y} \]

   6. Taylor expanded in z around 0

      \[\leadsto b \cdot \color{blue}{y} \]
   7. Step-by-step derivation

   8. Applied rewrites 42.2%

      \[\leadsto b \cdot \color{blue}{y} \]

   if -9.3999999999999995e26 < b < 2.8500000000000001e110

   1. Initial program 98.8%

      \[\left(\left(x - \left(y - 1\right) \cdot z\right) - \left(t - 1\right) \cdot a\right) + \left(\left(y + t\right) - 2\right) \cdot b \]
   2. Add Preprocessing

   3. Taylor expanded in t around inf

      \[\leadsto \color{blue}{t \cdot \left(b - a\right)} \]
   4. Step-by-step derivation

      1. *-commutative N/A

         \[\leadsto \color{blue}{\left(b - a\right) \cdot t} \]

      2. lower-*.f64 N/A

         \[\leadsto \color{blue}{\left(b - a\right) \cdot t} \]

      3. lower--.f64 39.0

         \[\leadsto \color{blue}{\left(b - a\right)} \cdot t \]

   5. Applied rewrites 39.0%

      \[\leadsto \color{blue}{\left(b - a\right) \cdot t} \]

   6. Taylor expanded in a around inf

      \[\leadsto \left(-1 \cdot a\right) \cdot t \]
   7. Step-by-step derivation

   8. Applied rewrites 35.8%

      \[\leadsto \left(-a\right) \cdot t \]

   if 2.8500000000000001e110 < b

   1. Initial program 87.0%

      \[\left(\left(x - \left(y - 1\right) \cdot z\right) - \left(t - 1\right) \cdot a\right) + \left(\left(y + t\right) - 2\right) \cdot b \]
   2. Add Preprocessing

   3. Taylor expanded in y around 0

      \[\leadsto \color{blue}{\left(x + b \cdot \left(t - 2\right)\right) - \left(-1 \cdot z + a \cdot \left(t - 1\right)\right)} \]
   4. Step-by-step derivation

      1. sub-neg N/A

         \[\leadsto \color{blue}{\left(x + b \cdot \left(t - 2\right)\right) + \left(\mathsf{neg}\left(\left(-1 \cdot z + a \cdot \left(t - 1\right)\right)\right)\right)} \]

      2. +-commutative N/A

         \[\leadsto \color{blue}{\left(\mathsf{neg}\left(\left(-1 \cdot z + a \cdot \left(t - 1\right)\right)\right)\right) + \left(x + b \cdot \left(t - 2\right)\right)} \]

      3. +-commutative N/A

         \[\leadsto \left(\mathsf{neg}\left(\color{blue}{\left(a \cdot \left(t - 1\right) + -1 \cdot z\right)}\right)\right) + \left(x + b \cdot \left(t - 2\right)\right) \]

      4. distribute-neg-in N/A

         \[\leadsto \color{blue}{\left(\left(\mathsf{neg}\left(a \cdot \left(t - 1\right)\right)\right) + \left(\mathsf{neg}\left(-1 \cdot z\right)\right)\right)} + \left(x + b \cdot \left(t - 2\right)\right) \]

      5. mul-1-neg N/A

         \[\leadsto \left(\left(\mathsf{neg}\left(a \cdot \left(t - 1\right)\right)\right) + \left(\mathsf{neg}\left(\color{blue}{\left(\mathsf{neg}\left(z\right)\right)}\right)\right)\right) + \left(x + b \cdot \left(t - 2\right)\right) \]

      6. remove-double-neg N/A

         \[\leadsto \left(\left(\mathsf{neg}\left(a \cdot \left(t - 1\right)\right)\right) + \color{blue}{z}\right) + \left(x + b \cdot \left(t - 2\right)\right) \]

      7. associate-+l+ N/A

         \[\leadsto \color{blue}{\left(\mathsf{neg}\left(a \cdot \left(t - 1\right)\right)\right) + \left(z + \left(x + b \cdot \left(t - 2\right)\right)\right)} \]

      8. *-commutative N/A

         \[\leadsto \left(\mathsf{neg}\left(\color{blue}{\left(t - 1\right) \cdot a}\right)\right) + \left(z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      9. distribute-lft-neg-in N/A

         \[\leadsto \color{blue}{\left(\mathsf{neg}\left(\left(t - 1\right)\right)\right) \cdot a} + \left(z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      10. mul-1-neg N/A

          \[\leadsto \color{blue}{\left(-1 \cdot \left(t - 1\right)\right)} \cdot a + \left(z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      11. lower-fma.f64 N/A

          \[\leadsto \color{blue}{\mathsf{fma}\left(-1 \cdot \left(t - 1\right), a, z + \left(x + b \cdot \left(t - 2\right)\right)\right)} \]

      12. sub-neg N/A

          \[\leadsto \mathsf{fma}\left(-1 \cdot \color{blue}{\left(t + \left(\mathsf{neg}\left(1\right)\right)\right)}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      13. metadata-eval N/A

          \[\leadsto \mathsf{fma}\left(-1 \cdot \left(t + \color{blue}{-1}\right), a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      14. distribute-lft-in N/A

          \[\leadsto \mathsf{fma}\left(\color{blue}{-1 \cdot t + -1 \cdot -1}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      15. metadata-eval N/A

          \[\leadsto \mathsf{fma}\left(-1 \cdot t + \color{blue}{1}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      16. +-commutative N/A

          \[\leadsto \mathsf{fma}\left(\color{blue}{1 + -1 \cdot t}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      17. neg-mul-1 N/A

          \[\leadsto \mathsf{fma}\left(1 + \color{blue}{\left(\mathsf{neg}\left(t\right)\right)}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      18. sub-neg N/A

          \[\leadsto \mathsf{fma}\left(\color{blue}{1 - t}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      19. lower--.f64 N/A

          \[\leadsto \mathsf{fma}\left(\color{blue}{1 - t}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      20. lower-+.f64 N/A

          \[\leadsto \mathsf{fma}\left(1 - t, a, \color{blue}{z + \left(x + b \cdot \left(t - 2\right)\right)}\right) \]

   5. Applied rewrites 72.7%

      \[\leadsto \color{blue}{\mathsf{fma}\left(1 - t, a, z + \mathsf{fma}\left(t - 2, b, x\right)\right)} \]

   6. Taylor expanded in b around inf

      \[\leadsto \color{blue}{b \cdot \left(\left(t + y\right) - 2\right)} \]
   7. Step-by-step derivation

      1. *-commutative N/A

         \[\leadsto \color{blue}{\left(\left(t + y\right) - 2\right) \cdot b} \]

      2. lower-*.f64 N/A

         \[\leadsto \color{blue}{\left(\left(t + y\right) - 2\right) \cdot b} \]

      3. lower--.f64 N/A

         \[\leadsto \color{blue}{\left(\left(t + y\right) - 2\right)} \cdot b \]

      4. lower-+.f64 79.1

         \[\leadsto \left(\color{blue}{\left(t + y\right)} - 2\right) \cdot b \]

   8. Applied rewrites 79.1%

      \[\leadsto \color{blue}{\left(\left(t + y\right) - 2\right) \cdot b} \]

   9. Taylor expanded in t around inf

      \[\leadsto b \cdot \color{blue}{t} \]
   10. Step-by-step derivation

   11. Applied rewrites 39.1%

       \[\leadsto b \cdot \color{blue}{t} \]
3. Recombined 3 regimes into one program.
4. Add Preprocessing

Alternative 18: 27.2% accurate, 2.1× speedup

Math

\[\begin{array}{l} \\ \begin{array}{l} \mathbf{if}\;y \leq -3.8 \cdot 10^{-17} \lor \neg \left(y \leq 4.4 \cdot 10^{+60}\right):\\ \;\;\;\;b \cdot y\\ \mathbf{else}:\\ \;\;\;\;b \cdot t\\ \end{array} \end{array} \]

FPCore

(FPCore (x y z t a b)
 :precision binary64
 (if (or (<= y -3.8e-17) (not (<= y 4.4e+60))) (* b y) (* b t)))

C

double code(double x, double y, double z, double t, double a, double b) {
	double tmp;
	if ((y <= -3.8e-17) || !(y <= 4.4e+60)) {
		tmp = b * y;
	} else {
		tmp = b * t;
	}
	return tmp;
}

Fortran

real(8) function code(x, y, z, t, a, b)
    real(8), intent (in) :: x
    real(8), intent (in) :: y
    real(8), intent (in) :: z
    real(8), intent (in) :: t
    real(8), intent (in) :: a
    real(8), intent (in) :: b
    real(8) :: tmp
    if ((y <= (-3.8d-17)) .or. (.not. (y <= 4.4d+60))) then
        tmp = b * y
    else
        tmp = b * t
    end if
    code = tmp
end function

Java

public static double code(double x, double y, double z, double t, double a, double b) {
	double tmp;
	if ((y <= -3.8e-17) || !(y <= 4.4e+60)) {
		tmp = b * y;
	} else {
		tmp = b * t;
	}
	return tmp;
}

Python

def code(x, y, z, t, a, b):
	tmp = 0
	if (y <= -3.8e-17) or not (y <= 4.4e+60):
		tmp = b * y
	else:
		tmp = b * t
	return tmp

Julia

function code(x, y, z, t, a, b)
	tmp = 0.0
	if ((y <= -3.8e-17) || !(y <= 4.4e+60))
		tmp = Float64(b * y);
	else
		tmp = Float64(b * t);
	end
	return tmp
end

MATLAB

function tmp_2 = code(x, y, z, t, a, b)
	tmp = 0.0;
	if ((y <= -3.8e-17) || ~((y <= 4.4e+60)))
		tmp = b * y;
	else
		tmp = b * t;
	end
	tmp_2 = tmp;
end

Wolfram

code[x_, y_, z_, t_, a_, b_] := If[Or[LessEqual[y, -3.8e-17], N[Not[LessEqual[y, 4.4e+60]], $MachinePrecision]], N[(b * y), $MachinePrecision], N[(b * t), $MachinePrecision]]

Derivation

1. Split input into 2 regimes

2. if y < -3.8000000000000001e-17 or 4.39999999999999992e60 < y

   1. Initial program 92.8%

      \[\left(\left(x - \left(y - 1\right) \cdot z\right) - \left(t - 1\right) \cdot a\right) + \left(\left(y + t\right) - 2\right) \cdot b \]
   2. Add Preprocessing

   3. Taylor expanded in y around inf

      \[\leadsto \color{blue}{y \cdot \left(b - z\right)} \]
   4. Step-by-step derivation

      1. *-commutative N/A

         \[\leadsto \color{blue}{\left(b - z\right) \cdot y} \]

      2. lower-*.f64 N/A

         \[\leadsto \color{blue}{\left(b - z\right) \cdot y} \]

      3. lower--.f64 61.3

         \[\leadsto \color{blue}{\left(b - z\right)} \cdot y \]

   5. Applied rewrites 61.3%

      \[\leadsto \color{blue}{\left(b - z\right) \cdot y} \]

   6. Taylor expanded in z around 0

      \[\leadsto b \cdot \color{blue}{y} \]
   7. Step-by-step derivation

   8. Applied rewrites 39.9%

      \[\leadsto b \cdot \color{blue}{y} \]

   if -3.8000000000000001e-17 < y < 4.39999999999999992e60

   1. Initial program 97.2%

      \[\left(\left(x - \left(y - 1\right) \cdot z\right) - \left(t - 1\right) \cdot a\right) + \left(\left(y + t\right) - 2\right) \cdot b \]
   2. Add Preprocessing

   3. Taylor expanded in y around 0

      \[\leadsto \color{blue}{\left(x + b \cdot \left(t - 2\right)\right) - \left(-1 \cdot z + a \cdot \left(t - 1\right)\right)} \]
   4. Step-by-step derivation

      1. sub-neg N/A

         \[\leadsto \color{blue}{\left(x + b \cdot \left(t - 2\right)\right) + \left(\mathsf{neg}\left(\left(-1 \cdot z + a \cdot \left(t - 1\right)\right)\right)\right)} \]

      2. +-commutative N/A

         \[\leadsto \color{blue}{\left(\mathsf{neg}\left(\left(-1 \cdot z + a \cdot \left(t - 1\right)\right)\right)\right) + \left(x + b \cdot \left(t - 2\right)\right)} \]

      3. +-commutative N/A

         \[\leadsto \left(\mathsf{neg}\left(\color{blue}{\left(a \cdot \left(t - 1\right) + -1 \cdot z\right)}\right)\right) + \left(x + b \cdot \left(t - 2\right)\right) \]

      4. distribute-neg-in N/A

         \[\leadsto \color{blue}{\left(\left(\mathsf{neg}\left(a \cdot \left(t - 1\right)\right)\right) + \left(\mathsf{neg}\left(-1 \cdot z\right)\right)\right)} + \left(x + b \cdot \left(t - 2\right)\right) \]

      5. mul-1-neg N/A

         \[\leadsto \left(\left(\mathsf{neg}\left(a \cdot \left(t - 1\right)\right)\right) + \left(\mathsf{neg}\left(\color{blue}{\left(\mathsf{neg}\left(z\right)\right)}\right)\right)\right) + \left(x + b \cdot \left(t - 2\right)\right) \]

      6. remove-double-neg N/A

         \[\leadsto \left(\left(\mathsf{neg}\left(a \cdot \left(t - 1\right)\right)\right) + \color{blue}{z}\right) + \left(x + b \cdot \left(t - 2\right)\right) \]

      7. associate-+l+ N/A

         \[\leadsto \color{blue}{\left(\mathsf{neg}\left(a \cdot \left(t - 1\right)\right)\right) + \left(z + \left(x + b \cdot \left(t - 2\right)\right)\right)} \]

      8. *-commutative N/A

         \[\leadsto \left(\mathsf{neg}\left(\color{blue}{\left(t - 1\right) \cdot a}\right)\right) + \left(z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      9. distribute-lft-neg-in N/A

         \[\leadsto \color{blue}{\left(\mathsf{neg}\left(\left(t - 1\right)\right)\right) \cdot a} + \left(z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      10. mul-1-neg N/A

          \[\leadsto \color{blue}{\left(-1 \cdot \left(t - 1\right)\right)} \cdot a + \left(z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      11. lower-fma.f64 N/A

          \[\leadsto \color{blue}{\mathsf{fma}\left(-1 \cdot \left(t - 1\right), a, z + \left(x + b \cdot \left(t - 2\right)\right)\right)} \]

      12. sub-neg N/A

          \[\leadsto \mathsf{fma}\left(-1 \cdot \color{blue}{\left(t + \left(\mathsf{neg}\left(1\right)\right)\right)}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      13. metadata-eval N/A

          \[\leadsto \mathsf{fma}\left(-1 \cdot \left(t + \color{blue}{-1}\right), a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      14. distribute-lft-in N/A

          \[\leadsto \mathsf{fma}\left(\color{blue}{-1 \cdot t + -1 \cdot -1}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      15. metadata-eval N/A

          \[\leadsto \mathsf{fma}\left(-1 \cdot t + \color{blue}{1}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      16. +-commutative N/A

          \[\leadsto \mathsf{fma}\left(\color{blue}{1 + -1 \cdot t}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      17. neg-mul-1 N/A

          \[\leadsto \mathsf{fma}\left(1 + \color{blue}{\left(\mathsf{neg}\left(t\right)\right)}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      18. sub-neg N/A

          \[\leadsto \mathsf{fma}\left(\color{blue}{1 - t}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      19. lower--.f64 N/A

          \[\leadsto \mathsf{fma}\left(\color{blue}{1 - t}, a, z + \left(x + b \cdot \left(t - 2\right)\right)\right) \]

      20. lower-+.f64 N/A

          \[\leadsto \mathsf{fma}\left(1 - t, a, \color{blue}{z + \left(x + b \cdot \left(t - 2\right)\right)}\right) \]

   5. Applied rewrites 98.6%

      \[\leadsto \color{blue}{\mathsf{fma}\left(1 - t, a, z + \mathsf{fma}\left(t - 2, b, x\right)\right)} \]

   6. Taylor expanded in b around inf

      \[\leadsto \color{blue}{b \cdot \left(\left(t + y\right) - 2\right)} \]
   7. Step-by-step derivation

      1. *-commutative N/A

         \[\leadsto \color{blue}{\left(\left(t + y\right) - 2\right) \cdot b} \]

      2. lower-*.f64 N/A

         \[\leadsto \color{blue}{\left(\left(t + y\right) - 2\right) \cdot b} \]

      3. lower--.f64 N/A

         \[\leadsto \color{blue}{\left(\left(t + y\right) - 2\right)} \cdot b \]

      4. lower-+.f64 28.1

         \[\leadsto \left(\color{blue}{\left(t + y\right)} - 2\right) \cdot b \]

   8. Applied rewrites 28.1%

      \[\leadsto \color{blue}{\left(\left(t + y\right) - 2\right) \cdot b} \]

   9. Taylor expanded in t around inf

      \[\leadsto b \cdot \color{blue}{t} \]
   10. Step-by-step derivation

   11. Applied rewrites 20.8%

       \[\leadsto b \cdot \color{blue}{t} \]
3. Recombined 2 regimes into one program.

4. Final simplification 29.0%

   \[\leadsto \begin{array}{l} \mathbf{if}\;y \leq -3.8 \cdot 10^{-17} \lor \neg \left(y \leq 4.4 \cdot 10^{+60}\right):\\ \;\;\;\;b \cdot y\\ \mathbf{else}:\\ \;\;\;\;b \cdot t\\ \end{array} \]
5. Add Preprocessing

Alternative 19: 17.7% accurate, 6.2× speedup

Math

\[\begin{array}{l} \\ b \cdot y \end{array} \]

FPCore

(FPCore (x y z t a b) :precision binary64 (* b y))

C

double code(double x, double y, double z, double t, double a, double b) {
	return b * y;
}

Fortran

real(8) function code(x, y, z, t, a, b)
    real(8), intent (in) :: x
    real(8), intent (in) :: y
    real(8), intent (in) :: z
    real(8), intent (in) :: t
    real(8), intent (in) :: a
    real(8), intent (in) :: b
    code = b * y
end function

Java

public static double code(double x, double y, double z, double t, double a, double b) {
	return b * y;
}

Python

def code(x, y, z, t, a, b):
	return b * y

Julia

function code(x, y, z, t, a, b)
	return Float64(b * y)
end

MATLAB

function tmp = code(x, y, z, t, a, b)
	tmp = b * y;
end

Wolfram

code[x_, y_, z_, t_, a_, b_] := N[(b * y), $MachinePrecision]

Derivation

1. Initial program 95.3%

   \[\left(\left(x - \left(y - 1\right) \cdot z\right) - \left(t - 1\right) \cdot a\right) + \left(\left(y + t\right) - 2\right) \cdot b \]
2. Add Preprocessing

3. Taylor expanded in y around inf

   \[\leadsto \color{blue}{y \cdot \left(b - z\right)} \]
4. Step-by-step derivation

   1. *-commutative N/A

      \[\leadsto \color{blue}{\left(b - z\right) \cdot y} \]

   2. lower-*.f64 N/A

      \[\leadsto \color{blue}{\left(b - z\right) \cdot y} \]

   3. lower--.f64 28.5

      \[\leadsto \color{blue}{\left(b - z\right)} \cdot y \]

5. Applied rewrites 28.5%

   \[\leadsto \color{blue}{\left(b - z\right) \cdot y} \]

6. Taylor expanded in z around 0

   \[\leadsto b \cdot \color{blue}{y} \]
7. Step-by-step derivation

8. Applied rewrites 18.6%

   \[\leadsto b \cdot \color{blue}{y} \]
9. Add Preprocessing

Reproduce

herbie shell --seed 2024318 
(FPCore (x y z t a b)
  :name "Statistics.Distribution.Beta:$centropy from math-functions-0.1.5.2"
  :precision binary64
  (+ (- (- x (* (- y 1.0) z)) (* (- t 1.0) a)) (* (- (+ y t) 2.0) b)))