Data.Colour.RGBSpace.HSL:hsl from colour-2.3.3, D

Percentage Accurate: 99.5% → 99.8%
Time: 5.3s
Alternatives: 14
Speedup: 1.9×

Specification

?
\[\begin{array}{l} \\ x + \left(\left(y - x\right) \cdot 6\right) \cdot \left(\frac{2}{3} - z\right) \end{array} \]
(FPCore (x y z)
 :precision binary64
 (+ x (* (* (- y x) 6.0) (- (/ 2.0 3.0) z))))
double code(double x, double y, double z) {
	return x + (((y - x) * 6.0) * ((2.0 / 3.0) - z));
}
real(8) function code(x, y, z)
    real(8), intent (in) :: x
    real(8), intent (in) :: y
    real(8), intent (in) :: z
    code = x + (((y - x) * 6.0d0) * ((2.0d0 / 3.0d0) - z))
end function
public static double code(double x, double y, double z) {
	return x + (((y - x) * 6.0) * ((2.0 / 3.0) - z));
}
def code(x, y, z):
	return x + (((y - x) * 6.0) * ((2.0 / 3.0) - z))
function code(x, y, z)
	return Float64(x + Float64(Float64(Float64(y - x) * 6.0) * Float64(Float64(2.0 / 3.0) - z)))
end
function tmp = code(x, y, z)
	tmp = x + (((y - x) * 6.0) * ((2.0 / 3.0) - z));
end
code[x_, y_, z_] := N[(x + N[(N[(N[(y - x), $MachinePrecision] * 6.0), $MachinePrecision] * N[(N[(2.0 / 3.0), $MachinePrecision] - z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}

\\
x + \left(\left(y - x\right) \cdot 6\right) \cdot \left(\frac{2}{3} - z\right)
\end{array}

Sampling outcomes in binary64 precision:

Local Percentage Accuracy vs ?

The average percentage accuracy by input value. Horizontal axis shows value of an input variable; the variable is choosen in the title. Vertical axis is accuracy; higher is better. Red represent the original program, while blue represents Herbie's suggestion. These can be toggled with buttons below the plot. The line is an average while dots represent individual samples.

Accuracy vs Speed?

Herbie found 14 alternatives:

AlternativeAccuracySpeedup
The accuracy (vertical axis) and speed (horizontal axis) of each alternatives. Up and to the right is better. The red square shows the initial program, and each blue circle shows an alternative.The line shows the best available speed-accuracy tradeoffs.

Initial Program: 99.5% accurate, 1.0× speedup?

\[\begin{array}{l} \\ x + \left(\left(y - x\right) \cdot 6\right) \cdot \left(\frac{2}{3} - z\right) \end{array} \]
(FPCore (x y z)
 :precision binary64
 (+ x (* (* (- y x) 6.0) (- (/ 2.0 3.0) z))))
double code(double x, double y, double z) {
	return x + (((y - x) * 6.0) * ((2.0 / 3.0) - z));
}
real(8) function code(x, y, z)
    real(8), intent (in) :: x
    real(8), intent (in) :: y
    real(8), intent (in) :: z
    code = x + (((y - x) * 6.0d0) * ((2.0d0 / 3.0d0) - z))
end function
public static double code(double x, double y, double z) {
	return x + (((y - x) * 6.0) * ((2.0 / 3.0) - z));
}
def code(x, y, z):
	return x + (((y - x) * 6.0) * ((2.0 / 3.0) - z))
function code(x, y, z)
	return Float64(x + Float64(Float64(Float64(y - x) * 6.0) * Float64(Float64(2.0 / 3.0) - z)))
end
function tmp = code(x, y, z)
	tmp = x + (((y - x) * 6.0) * ((2.0 / 3.0) - z));
end
code[x_, y_, z_] := N[(x + N[(N[(N[(y - x), $MachinePrecision] * 6.0), $MachinePrecision] * N[(N[(2.0 / 3.0), $MachinePrecision] - z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}

\\
x + \left(\left(y - x\right) \cdot 6\right) \cdot \left(\frac{2}{3} - z\right)
\end{array}

Alternative 1: 99.8% accurate, 1.9× speedup?

\[\begin{array}{l} \\ \mathsf{fma}\left(\mathsf{fma}\left(-6, z, 4\right), y - x, x\right) \end{array} \]
(FPCore (x y z) :precision binary64 (fma (fma -6.0 z 4.0) (- y x) x))
double code(double x, double y, double z) {
	return fma(fma(-6.0, z, 4.0), (y - x), x);
}
function code(x, y, z)
	return fma(fma(-6.0, z, 4.0), Float64(y - x), x)
end
code[x_, y_, z_] := N[(N[(-6.0 * z + 4.0), $MachinePrecision] * N[(y - x), $MachinePrecision] + x), $MachinePrecision]
\begin{array}{l}

\\
\mathsf{fma}\left(\mathsf{fma}\left(-6, z, 4\right), y - x, x\right)
\end{array}
Derivation
  1. Initial program 99.2%

    \[x + \left(\left(y - x\right) \cdot 6\right) \cdot \left(\frac{2}{3} - z\right) \]
  2. Add Preprocessing
  3. Step-by-step derivation
    1. lift-+.f64N/A

      \[\leadsto \color{blue}{x + \left(\left(y - x\right) \cdot 6\right) \cdot \left(\frac{2}{3} - z\right)} \]
    2. +-commutativeN/A

      \[\leadsto \color{blue}{\left(\left(y - x\right) \cdot 6\right) \cdot \left(\frac{2}{3} - z\right) + x} \]
    3. lift-*.f64N/A

      \[\leadsto \color{blue}{\left(\left(y - x\right) \cdot 6\right) \cdot \left(\frac{2}{3} - z\right)} + x \]
    4. lift-*.f64N/A

      \[\leadsto \color{blue}{\left(\left(y - x\right) \cdot 6\right)} \cdot \left(\frac{2}{3} - z\right) + x \]
    5. associate-*l*N/A

      \[\leadsto \color{blue}{\left(y - x\right) \cdot \left(6 \cdot \left(\frac{2}{3} - z\right)\right)} + x \]
    6. *-commutativeN/A

      \[\leadsto \color{blue}{\left(6 \cdot \left(\frac{2}{3} - z\right)\right) \cdot \left(y - x\right)} + x \]
    7. lower-fma.f64N/A

      \[\leadsto \color{blue}{\mathsf{fma}\left(6 \cdot \left(\frac{2}{3} - z\right), y - x, x\right)} \]
    8. lift--.f64N/A

      \[\leadsto \mathsf{fma}\left(6 \cdot \color{blue}{\left(\frac{2}{3} - z\right)}, y - x, x\right) \]
    9. sub-negN/A

      \[\leadsto \mathsf{fma}\left(6 \cdot \color{blue}{\left(\frac{2}{3} + \left(\mathsf{neg}\left(z\right)\right)\right)}, y - x, x\right) \]
    10. +-commutativeN/A

      \[\leadsto \mathsf{fma}\left(6 \cdot \color{blue}{\left(\left(\mathsf{neg}\left(z\right)\right) + \frac{2}{3}\right)}, y - x, x\right) \]
    11. distribute-lft-inN/A

      \[\leadsto \mathsf{fma}\left(\color{blue}{6 \cdot \left(\mathsf{neg}\left(z\right)\right) + 6 \cdot \frac{2}{3}}, y - x, x\right) \]
    12. neg-mul-1N/A

      \[\leadsto \mathsf{fma}\left(6 \cdot \color{blue}{\left(-1 \cdot z\right)} + 6 \cdot \frac{2}{3}, y - x, x\right) \]
    13. associate-*r*N/A

      \[\leadsto \mathsf{fma}\left(\color{blue}{\left(6 \cdot -1\right) \cdot z} + 6 \cdot \frac{2}{3}, y - x, x\right) \]
    14. metadata-evalN/A

      \[\leadsto \mathsf{fma}\left(\color{blue}{-6} \cdot z + 6 \cdot \frac{2}{3}, y - x, x\right) \]
    15. metadata-evalN/A

      \[\leadsto \mathsf{fma}\left(\color{blue}{\left(\mathsf{neg}\left(6\right)\right)} \cdot z + 6 \cdot \frac{2}{3}, y - x, x\right) \]
    16. lower-fma.f64N/A

      \[\leadsto \mathsf{fma}\left(\color{blue}{\mathsf{fma}\left(\mathsf{neg}\left(6\right), z, 6 \cdot \frac{2}{3}\right)}, y - x, x\right) \]
    17. metadata-evalN/A

      \[\leadsto \mathsf{fma}\left(\mathsf{fma}\left(\color{blue}{-6}, z, 6 \cdot \frac{2}{3}\right), y - x, x\right) \]
    18. lift-/.f64N/A

      \[\leadsto \mathsf{fma}\left(\mathsf{fma}\left(-6, z, 6 \cdot \color{blue}{\frac{2}{3}}\right), y - x, x\right) \]
    19. metadata-evalN/A

      \[\leadsto \mathsf{fma}\left(\mathsf{fma}\left(-6, z, 6 \cdot \color{blue}{\frac{2}{3}}\right), y - x, x\right) \]
    20. metadata-eval99.8

      \[\leadsto \mathsf{fma}\left(\mathsf{fma}\left(-6, z, \color{blue}{4}\right), y - x, x\right) \]
  4. Applied rewrites99.8%

    \[\leadsto \color{blue}{\mathsf{fma}\left(\mathsf{fma}\left(-6, z, 4\right), y - x, x\right)} \]
  5. Add Preprocessing

Alternative 2: 75.2% accurate, 0.4× speedup?

\[\begin{array}{l} \\ \begin{array}{l} t_0 := \frac{2}{3} - z\\ t_1 := \mathsf{fma}\left(6, z, -3\right) \cdot x\\ \mathbf{if}\;t\_0 \leq 0.66666666665:\\ \;\;\;\;t\_1\\ \mathbf{elif}\;t\_0 \leq 1:\\ \;\;\;\;\mathsf{fma}\left(-3, x, y \cdot 4\right)\\ \mathbf{elif}\;t\_0 \leq 10^{+156}:\\ \;\;\;\;y \cdot \mathsf{fma}\left(-6, z, 4\right)\\ \mathbf{else}:\\ \;\;\;\;t\_1\\ \end{array} \end{array} \]
(FPCore (x y z)
 :precision binary64
 (let* ((t_0 (- (/ 2.0 3.0) z)) (t_1 (* (fma 6.0 z -3.0) x)))
   (if (<= t_0 0.66666666665)
     t_1
     (if (<= t_0 1.0)
       (fma -3.0 x (* y 4.0))
       (if (<= t_0 1e+156) (* y (fma -6.0 z 4.0)) t_1)))))
double code(double x, double y, double z) {
	double t_0 = (2.0 / 3.0) - z;
	double t_1 = fma(6.0, z, -3.0) * x;
	double tmp;
	if (t_0 <= 0.66666666665) {
		tmp = t_1;
	} else if (t_0 <= 1.0) {
		tmp = fma(-3.0, x, (y * 4.0));
	} else if (t_0 <= 1e+156) {
		tmp = y * fma(-6.0, z, 4.0);
	} else {
		tmp = t_1;
	}
	return tmp;
}
function code(x, y, z)
	t_0 = Float64(Float64(2.0 / 3.0) - z)
	t_1 = Float64(fma(6.0, z, -3.0) * x)
	tmp = 0.0
	if (t_0 <= 0.66666666665)
		tmp = t_1;
	elseif (t_0 <= 1.0)
		tmp = fma(-3.0, x, Float64(y * 4.0));
	elseif (t_0 <= 1e+156)
		tmp = Float64(y * fma(-6.0, z, 4.0));
	else
		tmp = t_1;
	end
	return tmp
end
code[x_, y_, z_] := Block[{t$95$0 = N[(N[(2.0 / 3.0), $MachinePrecision] - z), $MachinePrecision]}, Block[{t$95$1 = N[(N[(6.0 * z + -3.0), $MachinePrecision] * x), $MachinePrecision]}, If[LessEqual[t$95$0, 0.66666666665], t$95$1, If[LessEqual[t$95$0, 1.0], N[(-3.0 * x + N[(y * 4.0), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$0, 1e+156], N[(y * N[(-6.0 * z + 4.0), $MachinePrecision]), $MachinePrecision], t$95$1]]]]]
\begin{array}{l}

\\
\begin{array}{l}
t_0 := \frac{2}{3} - z\\
t_1 := \mathsf{fma}\left(6, z, -3\right) \cdot x\\
\mathbf{if}\;t\_0 \leq 0.66666666665:\\
\;\;\;\;t\_1\\

\mathbf{elif}\;t\_0 \leq 1:\\
\;\;\;\;\mathsf{fma}\left(-3, x, y \cdot 4\right)\\

\mathbf{elif}\;t\_0 \leq 10^{+156}:\\
\;\;\;\;y \cdot \mathsf{fma}\left(-6, z, 4\right)\\

\mathbf{else}:\\
\;\;\;\;t\_1\\


\end{array}
\end{array}
Derivation
  1. Split input into 3 regimes
  2. if (-.f64 (/.f64 #s(literal 2 binary64) #s(literal 3 binary64)) z) < 0.666666666649999962 or 9.9999999999999998e155 < (-.f64 (/.f64 #s(literal 2 binary64) #s(literal 3 binary64)) z)

    1. Initial program 99.7%

      \[x + \left(\left(y - x\right) \cdot 6\right) \cdot \left(\frac{2}{3} - z\right) \]
    2. Add Preprocessing
    3. Step-by-step derivation
      1. lift-+.f64N/A

        \[\leadsto \color{blue}{x + \left(\left(y - x\right) \cdot 6\right) \cdot \left(\frac{2}{3} - z\right)} \]
      2. +-commutativeN/A

        \[\leadsto \color{blue}{\left(\left(y - x\right) \cdot 6\right) \cdot \left(\frac{2}{3} - z\right) + x} \]
      3. lift-*.f64N/A

        \[\leadsto \color{blue}{\left(\left(y - x\right) \cdot 6\right) \cdot \left(\frac{2}{3} - z\right)} + x \]
      4. lift-*.f64N/A

        \[\leadsto \color{blue}{\left(\left(y - x\right) \cdot 6\right)} \cdot \left(\frac{2}{3} - z\right) + x \]
      5. associate-*l*N/A

        \[\leadsto \color{blue}{\left(y - x\right) \cdot \left(6 \cdot \left(\frac{2}{3} - z\right)\right)} + x \]
      6. *-commutativeN/A

        \[\leadsto \color{blue}{\left(6 \cdot \left(\frac{2}{3} - z\right)\right) \cdot \left(y - x\right)} + x \]
      7. lower-fma.f64N/A

        \[\leadsto \color{blue}{\mathsf{fma}\left(6 \cdot \left(\frac{2}{3} - z\right), y - x, x\right)} \]
      8. lift--.f64N/A

        \[\leadsto \mathsf{fma}\left(6 \cdot \color{blue}{\left(\frac{2}{3} - z\right)}, y - x, x\right) \]
      9. sub-negN/A

        \[\leadsto \mathsf{fma}\left(6 \cdot \color{blue}{\left(\frac{2}{3} + \left(\mathsf{neg}\left(z\right)\right)\right)}, y - x, x\right) \]
      10. +-commutativeN/A

        \[\leadsto \mathsf{fma}\left(6 \cdot \color{blue}{\left(\left(\mathsf{neg}\left(z\right)\right) + \frac{2}{3}\right)}, y - x, x\right) \]
      11. distribute-lft-inN/A

        \[\leadsto \mathsf{fma}\left(\color{blue}{6 \cdot \left(\mathsf{neg}\left(z\right)\right) + 6 \cdot \frac{2}{3}}, y - x, x\right) \]
      12. neg-mul-1N/A

        \[\leadsto \mathsf{fma}\left(6 \cdot \color{blue}{\left(-1 \cdot z\right)} + 6 \cdot \frac{2}{3}, y - x, x\right) \]
      13. associate-*r*N/A

        \[\leadsto \mathsf{fma}\left(\color{blue}{\left(6 \cdot -1\right) \cdot z} + 6 \cdot \frac{2}{3}, y - x, x\right) \]
      14. metadata-evalN/A

        \[\leadsto \mathsf{fma}\left(\color{blue}{-6} \cdot z + 6 \cdot \frac{2}{3}, y - x, x\right) \]
      15. metadata-evalN/A

        \[\leadsto \mathsf{fma}\left(\color{blue}{\left(\mathsf{neg}\left(6\right)\right)} \cdot z + 6 \cdot \frac{2}{3}, y - x, x\right) \]
      16. lower-fma.f64N/A

        \[\leadsto \mathsf{fma}\left(\color{blue}{\mathsf{fma}\left(\mathsf{neg}\left(6\right), z, 6 \cdot \frac{2}{3}\right)}, y - x, x\right) \]
      17. metadata-evalN/A

        \[\leadsto \mathsf{fma}\left(\mathsf{fma}\left(\color{blue}{-6}, z, 6 \cdot \frac{2}{3}\right), y - x, x\right) \]
      18. lift-/.f64N/A

        \[\leadsto \mathsf{fma}\left(\mathsf{fma}\left(-6, z, 6 \cdot \color{blue}{\frac{2}{3}}\right), y - x, x\right) \]
      19. metadata-evalN/A

        \[\leadsto \mathsf{fma}\left(\mathsf{fma}\left(-6, z, 6 \cdot \color{blue}{\frac{2}{3}}\right), y - x, x\right) \]
      20. metadata-eval99.8

        \[\leadsto \mathsf{fma}\left(\mathsf{fma}\left(-6, z, \color{blue}{4}\right), y - x, x\right) \]
    4. Applied rewrites99.8%

      \[\leadsto \color{blue}{\mathsf{fma}\left(\mathsf{fma}\left(-6, z, 4\right), y - x, x\right)} \]
    5. Taylor expanded in x around 0

      \[\leadsto \color{blue}{y \cdot \left(4 + -6 \cdot z\right)} \]
    6. Step-by-step derivation
      1. *-commutativeN/A

        \[\leadsto \color{blue}{\left(4 + -6 \cdot z\right) \cdot y} \]
      2. lower-*.f64N/A

        \[\leadsto \color{blue}{\left(4 + -6 \cdot z\right) \cdot y} \]
      3. +-commutativeN/A

        \[\leadsto \color{blue}{\left(-6 \cdot z + 4\right)} \cdot y \]
      4. lower-fma.f6443.8

        \[\leadsto \color{blue}{\mathsf{fma}\left(-6, z, 4\right)} \cdot y \]
    7. Applied rewrites43.8%

      \[\leadsto \color{blue}{\mathsf{fma}\left(-6, z, 4\right) \cdot y} \]
    8. Taylor expanded in x around inf

      \[\leadsto \color{blue}{x \cdot \left(1 + -1 \cdot \left(4 + -6 \cdot z\right)\right)} \]
    9. Step-by-step derivation
      1. *-commutativeN/A

        \[\leadsto \color{blue}{\left(1 + -1 \cdot \left(4 + -6 \cdot z\right)\right) \cdot x} \]
      2. distribute-lft-inN/A

        \[\leadsto \left(1 + \color{blue}{\left(-1 \cdot 4 + -1 \cdot \left(-6 \cdot z\right)\right)}\right) \cdot x \]
      3. metadata-evalN/A

        \[\leadsto \left(1 + \left(\color{blue}{-4} + -1 \cdot \left(-6 \cdot z\right)\right)\right) \cdot x \]
      4. metadata-evalN/A

        \[\leadsto \left(1 + \left(\color{blue}{\frac{2}{3} \cdot -6} + -1 \cdot \left(-6 \cdot z\right)\right)\right) \cdot x \]
      5. *-commutativeN/A

        \[\leadsto \left(1 + \left(\frac{2}{3} \cdot -6 + -1 \cdot \color{blue}{\left(z \cdot -6\right)}\right)\right) \cdot x \]
      6. associate-*r*N/A

        \[\leadsto \left(1 + \left(\frac{2}{3} \cdot -6 + \color{blue}{\left(-1 \cdot z\right) \cdot -6}\right)\right) \cdot x \]
      7. distribute-rgt-inN/A

        \[\leadsto \left(1 + \color{blue}{-6 \cdot \left(\frac{2}{3} + -1 \cdot z\right)}\right) \cdot x \]
      8. mul-1-negN/A

        \[\leadsto \left(1 + -6 \cdot \left(\frac{2}{3} + \color{blue}{\left(\mathsf{neg}\left(z\right)\right)}\right)\right) \cdot x \]
      9. sub-negN/A

        \[\leadsto \left(1 + -6 \cdot \color{blue}{\left(\frac{2}{3} - z\right)}\right) \cdot x \]
      10. lower-*.f64N/A

        \[\leadsto \color{blue}{\left(1 + -6 \cdot \left(\frac{2}{3} - z\right)\right) \cdot x} \]
    10. Applied rewrites64.6%

      \[\leadsto \color{blue}{\mathsf{fma}\left(6, z, -3\right) \cdot x} \]

    if 0.666666666649999962 < (-.f64 (/.f64 #s(literal 2 binary64) #s(literal 3 binary64)) z) < 1

    1. Initial program 98.9%

      \[x + \left(\left(y - x\right) \cdot 6\right) \cdot \left(\frac{2}{3} - z\right) \]
    2. Add Preprocessing
    3. Taylor expanded in z around 0

      \[\leadsto \color{blue}{x + 4 \cdot \left(y - x\right)} \]
    4. Step-by-step derivation
      1. +-commutativeN/A

        \[\leadsto \color{blue}{4 \cdot \left(y - x\right) + x} \]
      2. lower-fma.f64N/A

        \[\leadsto \color{blue}{\mathsf{fma}\left(4, y - x, x\right)} \]
      3. lower--.f6499.8

        \[\leadsto \mathsf{fma}\left(4, \color{blue}{y - x}, x\right) \]
    5. Applied rewrites99.8%

      \[\leadsto \color{blue}{\mathsf{fma}\left(4, y - x, x\right)} \]
    6. Taylor expanded in x around 0

      \[\leadsto -3 \cdot x + \color{blue}{4 \cdot y} \]
    7. Step-by-step derivation
      1. Applied rewrites99.9%

        \[\leadsto \mathsf{fma}\left(-3, \color{blue}{x}, 4 \cdot y\right) \]

      if 1 < (-.f64 (/.f64 #s(literal 2 binary64) #s(literal 3 binary64)) z) < 9.9999999999999998e155

      1. Initial program 99.3%

        \[x + \left(\left(y - x\right) \cdot 6\right) \cdot \left(\frac{2}{3} - z\right) \]
      2. Add Preprocessing
      3. Taylor expanded in x around 0

        \[\leadsto \color{blue}{6 \cdot \left(y \cdot \left(\frac{2}{3} - z\right)\right)} \]
      4. Step-by-step derivation
        1. *-commutativeN/A

          \[\leadsto 6 \cdot \color{blue}{\left(\left(\frac{2}{3} - z\right) \cdot y\right)} \]
        2. associate-*r*N/A

          \[\leadsto \color{blue}{\left(6 \cdot \left(\frac{2}{3} - z\right)\right) \cdot y} \]
        3. lower-*.f64N/A

          \[\leadsto \color{blue}{\left(6 \cdot \left(\frac{2}{3} - z\right)\right) \cdot y} \]
        4. sub-negN/A

          \[\leadsto \left(6 \cdot \color{blue}{\left(\frac{2}{3} + \left(\mathsf{neg}\left(z\right)\right)\right)}\right) \cdot y \]
        5. neg-mul-1N/A

          \[\leadsto \left(6 \cdot \left(\frac{2}{3} + \color{blue}{-1 \cdot z}\right)\right) \cdot y \]
        6. +-commutativeN/A

          \[\leadsto \left(6 \cdot \color{blue}{\left(-1 \cdot z + \frac{2}{3}\right)}\right) \cdot y \]
        7. distribute-lft-inN/A

          \[\leadsto \color{blue}{\left(6 \cdot \left(-1 \cdot z\right) + 6 \cdot \frac{2}{3}\right)} \cdot y \]
        8. metadata-evalN/A

          \[\leadsto \left(6 \cdot \left(-1 \cdot z\right) + \color{blue}{4}\right) \cdot y \]
        9. associate-*r*N/A

          \[\leadsto \left(\color{blue}{\left(6 \cdot -1\right) \cdot z} + 4\right) \cdot y \]
        10. metadata-evalN/A

          \[\leadsto \left(\color{blue}{-6} \cdot z + 4\right) \cdot y \]
        11. lower-fma.f6460.5

          \[\leadsto \color{blue}{\mathsf{fma}\left(-6, z, 4\right)} \cdot y \]
      5. Applied rewrites60.5%

        \[\leadsto \color{blue}{\mathsf{fma}\left(-6, z, 4\right) \cdot y} \]
    8. Recombined 3 regimes into one program.
    9. Final simplification83.7%

      \[\leadsto \begin{array}{l} \mathbf{if}\;\frac{2}{3} - z \leq 0.66666666665:\\ \;\;\;\;\mathsf{fma}\left(6, z, -3\right) \cdot x\\ \mathbf{elif}\;\frac{2}{3} - z \leq 1:\\ \;\;\;\;\mathsf{fma}\left(-3, x, y \cdot 4\right)\\ \mathbf{elif}\;\frac{2}{3} - z \leq 10^{+156}:\\ \;\;\;\;y \cdot \mathsf{fma}\left(-6, z, 4\right)\\ \mathbf{else}:\\ \;\;\;\;\mathsf{fma}\left(6, z, -3\right) \cdot x\\ \end{array} \]
    10. Add Preprocessing

    Alternative 3: 74.7% accurate, 0.4× speedup?

    \[\begin{array}{l} \\ \begin{array}{l} t_0 := \frac{2}{3} - z\\ t_1 := \left(6 \cdot x\right) \cdot z\\ \mathbf{if}\;t\_0 \leq -10000000000:\\ \;\;\;\;t\_1\\ \mathbf{elif}\;t\_0 \leq 1:\\ \;\;\;\;\mathsf{fma}\left(-3, x, y \cdot 4\right)\\ \mathbf{elif}\;t\_0 \leq 10^{+156}:\\ \;\;\;\;y \cdot \mathsf{fma}\left(-6, z, 4\right)\\ \mathbf{else}:\\ \;\;\;\;t\_1\\ \end{array} \end{array} \]
    (FPCore (x y z)
     :precision binary64
     (let* ((t_0 (- (/ 2.0 3.0) z)) (t_1 (* (* 6.0 x) z)))
       (if (<= t_0 -10000000000.0)
         t_1
         (if (<= t_0 1.0)
           (fma -3.0 x (* y 4.0))
           (if (<= t_0 1e+156) (* y (fma -6.0 z 4.0)) t_1)))))
    double code(double x, double y, double z) {
    	double t_0 = (2.0 / 3.0) - z;
    	double t_1 = (6.0 * x) * z;
    	double tmp;
    	if (t_0 <= -10000000000.0) {
    		tmp = t_1;
    	} else if (t_0 <= 1.0) {
    		tmp = fma(-3.0, x, (y * 4.0));
    	} else if (t_0 <= 1e+156) {
    		tmp = y * fma(-6.0, z, 4.0);
    	} else {
    		tmp = t_1;
    	}
    	return tmp;
    }
    
    function code(x, y, z)
    	t_0 = Float64(Float64(2.0 / 3.0) - z)
    	t_1 = Float64(Float64(6.0 * x) * z)
    	tmp = 0.0
    	if (t_0 <= -10000000000.0)
    		tmp = t_1;
    	elseif (t_0 <= 1.0)
    		tmp = fma(-3.0, x, Float64(y * 4.0));
    	elseif (t_0 <= 1e+156)
    		tmp = Float64(y * fma(-6.0, z, 4.0));
    	else
    		tmp = t_1;
    	end
    	return tmp
    end
    
    code[x_, y_, z_] := Block[{t$95$0 = N[(N[(2.0 / 3.0), $MachinePrecision] - z), $MachinePrecision]}, Block[{t$95$1 = N[(N[(6.0 * x), $MachinePrecision] * z), $MachinePrecision]}, If[LessEqual[t$95$0, -10000000000.0], t$95$1, If[LessEqual[t$95$0, 1.0], N[(-3.0 * x + N[(y * 4.0), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$0, 1e+156], N[(y * N[(-6.0 * z + 4.0), $MachinePrecision]), $MachinePrecision], t$95$1]]]]]
    
    \begin{array}{l}
    
    \\
    \begin{array}{l}
    t_0 := \frac{2}{3} - z\\
    t_1 := \left(6 \cdot x\right) \cdot z\\
    \mathbf{if}\;t\_0 \leq -10000000000:\\
    \;\;\;\;t\_1\\
    
    \mathbf{elif}\;t\_0 \leq 1:\\
    \;\;\;\;\mathsf{fma}\left(-3, x, y \cdot 4\right)\\
    
    \mathbf{elif}\;t\_0 \leq 10^{+156}:\\
    \;\;\;\;y \cdot \mathsf{fma}\left(-6, z, 4\right)\\
    
    \mathbf{else}:\\
    \;\;\;\;t\_1\\
    
    
    \end{array}
    \end{array}
    
    Derivation
    1. Split input into 3 regimes
    2. if (-.f64 (/.f64 #s(literal 2 binary64) #s(literal 3 binary64)) z) < -1e10 or 9.9999999999999998e155 < (-.f64 (/.f64 #s(literal 2 binary64) #s(literal 3 binary64)) z)

      1. Initial program 99.7%

        \[x + \left(\left(y - x\right) \cdot 6\right) \cdot \left(\frac{2}{3} - z\right) \]
      2. Add Preprocessing
      3. Step-by-step derivation
        1. lift-+.f64N/A

          \[\leadsto \color{blue}{x + \left(\left(y - x\right) \cdot 6\right) \cdot \left(\frac{2}{3} - z\right)} \]
        2. +-commutativeN/A

          \[\leadsto \color{blue}{\left(\left(y - x\right) \cdot 6\right) \cdot \left(\frac{2}{3} - z\right) + x} \]
        3. lift-*.f64N/A

          \[\leadsto \color{blue}{\left(\left(y - x\right) \cdot 6\right) \cdot \left(\frac{2}{3} - z\right)} + x \]
        4. lift-*.f64N/A

          \[\leadsto \color{blue}{\left(\left(y - x\right) \cdot 6\right)} \cdot \left(\frac{2}{3} - z\right) + x \]
        5. associate-*l*N/A

          \[\leadsto \color{blue}{\left(y - x\right) \cdot \left(6 \cdot \left(\frac{2}{3} - z\right)\right)} + x \]
        6. *-commutativeN/A

          \[\leadsto \color{blue}{\left(6 \cdot \left(\frac{2}{3} - z\right)\right) \cdot \left(y - x\right)} + x \]
        7. lower-fma.f64N/A

          \[\leadsto \color{blue}{\mathsf{fma}\left(6 \cdot \left(\frac{2}{3} - z\right), y - x, x\right)} \]
        8. lift--.f64N/A

          \[\leadsto \mathsf{fma}\left(6 \cdot \color{blue}{\left(\frac{2}{3} - z\right)}, y - x, x\right) \]
        9. sub-negN/A

          \[\leadsto \mathsf{fma}\left(6 \cdot \color{blue}{\left(\frac{2}{3} + \left(\mathsf{neg}\left(z\right)\right)\right)}, y - x, x\right) \]
        10. +-commutativeN/A

          \[\leadsto \mathsf{fma}\left(6 \cdot \color{blue}{\left(\left(\mathsf{neg}\left(z\right)\right) + \frac{2}{3}\right)}, y - x, x\right) \]
        11. distribute-lft-inN/A

          \[\leadsto \mathsf{fma}\left(\color{blue}{6 \cdot \left(\mathsf{neg}\left(z\right)\right) + 6 \cdot \frac{2}{3}}, y - x, x\right) \]
        12. neg-mul-1N/A

          \[\leadsto \mathsf{fma}\left(6 \cdot \color{blue}{\left(-1 \cdot z\right)} + 6 \cdot \frac{2}{3}, y - x, x\right) \]
        13. associate-*r*N/A

          \[\leadsto \mathsf{fma}\left(\color{blue}{\left(6 \cdot -1\right) \cdot z} + 6 \cdot \frac{2}{3}, y - x, x\right) \]
        14. metadata-evalN/A

          \[\leadsto \mathsf{fma}\left(\color{blue}{-6} \cdot z + 6 \cdot \frac{2}{3}, y - x, x\right) \]
        15. metadata-evalN/A

          \[\leadsto \mathsf{fma}\left(\color{blue}{\left(\mathsf{neg}\left(6\right)\right)} \cdot z + 6 \cdot \frac{2}{3}, y - x, x\right) \]
        16. lower-fma.f64N/A

          \[\leadsto \mathsf{fma}\left(\color{blue}{\mathsf{fma}\left(\mathsf{neg}\left(6\right), z, 6 \cdot \frac{2}{3}\right)}, y - x, x\right) \]
        17. metadata-evalN/A

          \[\leadsto \mathsf{fma}\left(\mathsf{fma}\left(\color{blue}{-6}, z, 6 \cdot \frac{2}{3}\right), y - x, x\right) \]
        18. lift-/.f64N/A

          \[\leadsto \mathsf{fma}\left(\mathsf{fma}\left(-6, z, 6 \cdot \color{blue}{\frac{2}{3}}\right), y - x, x\right) \]
        19. metadata-evalN/A

          \[\leadsto \mathsf{fma}\left(\mathsf{fma}\left(-6, z, 6 \cdot \color{blue}{\frac{2}{3}}\right), y - x, x\right) \]
        20. metadata-eval99.8

          \[\leadsto \mathsf{fma}\left(\mathsf{fma}\left(-6, z, \color{blue}{4}\right), y - x, x\right) \]
      4. Applied rewrites99.8%

        \[\leadsto \color{blue}{\mathsf{fma}\left(\mathsf{fma}\left(-6, z, 4\right), y - x, x\right)} \]
      5. Taylor expanded in z around inf

        \[\leadsto \color{blue}{-6 \cdot \left(z \cdot \left(y - x\right)\right)} \]
      6. Step-by-step derivation
        1. *-commutativeN/A

          \[\leadsto -6 \cdot \color{blue}{\left(\left(y - x\right) \cdot z\right)} \]
        2. associate-*r*N/A

          \[\leadsto \color{blue}{\left(-6 \cdot \left(y - x\right)\right) \cdot z} \]
        3. lower-*.f64N/A

          \[\leadsto \color{blue}{\left(-6 \cdot \left(y - x\right)\right) \cdot z} \]
        4. lower-*.f64N/A

          \[\leadsto \color{blue}{\left(-6 \cdot \left(y - x\right)\right)} \cdot z \]
        5. lower--.f6499.1

          \[\leadsto \left(-6 \cdot \color{blue}{\left(y - x\right)}\right) \cdot z \]
      7. Applied rewrites99.1%

        \[\leadsto \color{blue}{\left(-6 \cdot \left(y - x\right)\right) \cdot z} \]
      8. Taylor expanded in x around inf

        \[\leadsto \left(6 \cdot x\right) \cdot z \]
      9. Step-by-step derivation
        1. Applied rewrites63.6%

          \[\leadsto \left(6 \cdot x\right) \cdot z \]

        if -1e10 < (-.f64 (/.f64 #s(literal 2 binary64) #s(literal 3 binary64)) z) < 1

        1. Initial program 98.8%

          \[x + \left(\left(y - x\right) \cdot 6\right) \cdot \left(\frac{2}{3} - z\right) \]
        2. Add Preprocessing
        3. Taylor expanded in z around 0

          \[\leadsto \color{blue}{x + 4 \cdot \left(y - x\right)} \]
        4. Step-by-step derivation
          1. +-commutativeN/A

            \[\leadsto \color{blue}{4 \cdot \left(y - x\right) + x} \]
          2. lower-fma.f64N/A

            \[\leadsto \color{blue}{\mathsf{fma}\left(4, y - x, x\right)} \]
          3. lower--.f6499.6

            \[\leadsto \mathsf{fma}\left(4, \color{blue}{y - x}, x\right) \]
        5. Applied rewrites99.6%

          \[\leadsto \color{blue}{\mathsf{fma}\left(4, y - x, x\right)} \]
        6. Taylor expanded in x around 0

          \[\leadsto -3 \cdot x + \color{blue}{4 \cdot y} \]
        7. Step-by-step derivation
          1. Applied rewrites99.7%

            \[\leadsto \mathsf{fma}\left(-3, \color{blue}{x}, 4 \cdot y\right) \]

          if 1 < (-.f64 (/.f64 #s(literal 2 binary64) #s(literal 3 binary64)) z) < 9.9999999999999998e155

          1. Initial program 99.3%

            \[x + \left(\left(y - x\right) \cdot 6\right) \cdot \left(\frac{2}{3} - z\right) \]
          2. Add Preprocessing
          3. Taylor expanded in x around 0

            \[\leadsto \color{blue}{6 \cdot \left(y \cdot \left(\frac{2}{3} - z\right)\right)} \]
          4. Step-by-step derivation
            1. *-commutativeN/A

              \[\leadsto 6 \cdot \color{blue}{\left(\left(\frac{2}{3} - z\right) \cdot y\right)} \]
            2. associate-*r*N/A

              \[\leadsto \color{blue}{\left(6 \cdot \left(\frac{2}{3} - z\right)\right) \cdot y} \]
            3. lower-*.f64N/A

              \[\leadsto \color{blue}{\left(6 \cdot \left(\frac{2}{3} - z\right)\right) \cdot y} \]
            4. sub-negN/A

              \[\leadsto \left(6 \cdot \color{blue}{\left(\frac{2}{3} + \left(\mathsf{neg}\left(z\right)\right)\right)}\right) \cdot y \]
            5. neg-mul-1N/A

              \[\leadsto \left(6 \cdot \left(\frac{2}{3} + \color{blue}{-1 \cdot z}\right)\right) \cdot y \]
            6. +-commutativeN/A

              \[\leadsto \left(6 \cdot \color{blue}{\left(-1 \cdot z + \frac{2}{3}\right)}\right) \cdot y \]
            7. distribute-lft-inN/A

              \[\leadsto \color{blue}{\left(6 \cdot \left(-1 \cdot z\right) + 6 \cdot \frac{2}{3}\right)} \cdot y \]
            8. metadata-evalN/A

              \[\leadsto \left(6 \cdot \left(-1 \cdot z\right) + \color{blue}{4}\right) \cdot y \]
            9. associate-*r*N/A

              \[\leadsto \left(\color{blue}{\left(6 \cdot -1\right) \cdot z} + 4\right) \cdot y \]
            10. metadata-evalN/A

              \[\leadsto \left(\color{blue}{-6} \cdot z + 4\right) \cdot y \]
            11. lower-fma.f6460.5

              \[\leadsto \color{blue}{\mathsf{fma}\left(-6, z, 4\right)} \cdot y \]
          5. Applied rewrites60.5%

            \[\leadsto \color{blue}{\mathsf{fma}\left(-6, z, 4\right) \cdot y} \]
        8. Recombined 3 regimes into one program.
        9. Final simplification83.4%

          \[\leadsto \begin{array}{l} \mathbf{if}\;\frac{2}{3} - z \leq -10000000000:\\ \;\;\;\;\left(6 \cdot x\right) \cdot z\\ \mathbf{elif}\;\frac{2}{3} - z \leq 1:\\ \;\;\;\;\mathsf{fma}\left(-3, x, y \cdot 4\right)\\ \mathbf{elif}\;\frac{2}{3} - z \leq 10^{+156}:\\ \;\;\;\;y \cdot \mathsf{fma}\left(-6, z, 4\right)\\ \mathbf{else}:\\ \;\;\;\;\left(6 \cdot x\right) \cdot z\\ \end{array} \]
        10. Add Preprocessing

        Alternative 4: 74.5% accurate, 0.4× speedup?

        \[\begin{array}{l} \\ \begin{array}{l} t_0 := \frac{2}{3} - z\\ t_1 := \left(6 \cdot x\right) \cdot z\\ \mathbf{if}\;t\_0 \leq -10000000000:\\ \;\;\;\;t\_1\\ \mathbf{elif}\;t\_0 \leq 1:\\ \;\;\;\;\mathsf{fma}\left(-3, x, y \cdot 4\right)\\ \mathbf{elif}\;t\_0 \leq 10^{+156}:\\ \;\;\;\;\left(z \cdot -6\right) \cdot y\\ \mathbf{else}:\\ \;\;\;\;t\_1\\ \end{array} \end{array} \]
        (FPCore (x y z)
         :precision binary64
         (let* ((t_0 (- (/ 2.0 3.0) z)) (t_1 (* (* 6.0 x) z)))
           (if (<= t_0 -10000000000.0)
             t_1
             (if (<= t_0 1.0)
               (fma -3.0 x (* y 4.0))
               (if (<= t_0 1e+156) (* (* z -6.0) y) t_1)))))
        double code(double x, double y, double z) {
        	double t_0 = (2.0 / 3.0) - z;
        	double t_1 = (6.0 * x) * z;
        	double tmp;
        	if (t_0 <= -10000000000.0) {
        		tmp = t_1;
        	} else if (t_0 <= 1.0) {
        		tmp = fma(-3.0, x, (y * 4.0));
        	} else if (t_0 <= 1e+156) {
        		tmp = (z * -6.0) * y;
        	} else {
        		tmp = t_1;
        	}
        	return tmp;
        }
        
        function code(x, y, z)
        	t_0 = Float64(Float64(2.0 / 3.0) - z)
        	t_1 = Float64(Float64(6.0 * x) * z)
        	tmp = 0.0
        	if (t_0 <= -10000000000.0)
        		tmp = t_1;
        	elseif (t_0 <= 1.0)
        		tmp = fma(-3.0, x, Float64(y * 4.0));
        	elseif (t_0 <= 1e+156)
        		tmp = Float64(Float64(z * -6.0) * y);
        	else
        		tmp = t_1;
        	end
        	return tmp
        end
        
        code[x_, y_, z_] := Block[{t$95$0 = N[(N[(2.0 / 3.0), $MachinePrecision] - z), $MachinePrecision]}, Block[{t$95$1 = N[(N[(6.0 * x), $MachinePrecision] * z), $MachinePrecision]}, If[LessEqual[t$95$0, -10000000000.0], t$95$1, If[LessEqual[t$95$0, 1.0], N[(-3.0 * x + N[(y * 4.0), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$0, 1e+156], N[(N[(z * -6.0), $MachinePrecision] * y), $MachinePrecision], t$95$1]]]]]
        
        \begin{array}{l}
        
        \\
        \begin{array}{l}
        t_0 := \frac{2}{3} - z\\
        t_1 := \left(6 \cdot x\right) \cdot z\\
        \mathbf{if}\;t\_0 \leq -10000000000:\\
        \;\;\;\;t\_1\\
        
        \mathbf{elif}\;t\_0 \leq 1:\\
        \;\;\;\;\mathsf{fma}\left(-3, x, y \cdot 4\right)\\
        
        \mathbf{elif}\;t\_0 \leq 10^{+156}:\\
        \;\;\;\;\left(z \cdot -6\right) \cdot y\\
        
        \mathbf{else}:\\
        \;\;\;\;t\_1\\
        
        
        \end{array}
        \end{array}
        
        Derivation
        1. Split input into 3 regimes
        2. if (-.f64 (/.f64 #s(literal 2 binary64) #s(literal 3 binary64)) z) < -1e10 or 9.9999999999999998e155 < (-.f64 (/.f64 #s(literal 2 binary64) #s(literal 3 binary64)) z)

          1. Initial program 99.7%

            \[x + \left(\left(y - x\right) \cdot 6\right) \cdot \left(\frac{2}{3} - z\right) \]
          2. Add Preprocessing
          3. Step-by-step derivation
            1. lift-+.f64N/A

              \[\leadsto \color{blue}{x + \left(\left(y - x\right) \cdot 6\right) \cdot \left(\frac{2}{3} - z\right)} \]
            2. +-commutativeN/A

              \[\leadsto \color{blue}{\left(\left(y - x\right) \cdot 6\right) \cdot \left(\frac{2}{3} - z\right) + x} \]
            3. lift-*.f64N/A

              \[\leadsto \color{blue}{\left(\left(y - x\right) \cdot 6\right) \cdot \left(\frac{2}{3} - z\right)} + x \]
            4. lift-*.f64N/A

              \[\leadsto \color{blue}{\left(\left(y - x\right) \cdot 6\right)} \cdot \left(\frac{2}{3} - z\right) + x \]
            5. associate-*l*N/A

              \[\leadsto \color{blue}{\left(y - x\right) \cdot \left(6 \cdot \left(\frac{2}{3} - z\right)\right)} + x \]
            6. *-commutativeN/A

              \[\leadsto \color{blue}{\left(6 \cdot \left(\frac{2}{3} - z\right)\right) \cdot \left(y - x\right)} + x \]
            7. lower-fma.f64N/A

              \[\leadsto \color{blue}{\mathsf{fma}\left(6 \cdot \left(\frac{2}{3} - z\right), y - x, x\right)} \]
            8. lift--.f64N/A

              \[\leadsto \mathsf{fma}\left(6 \cdot \color{blue}{\left(\frac{2}{3} - z\right)}, y - x, x\right) \]
            9. sub-negN/A

              \[\leadsto \mathsf{fma}\left(6 \cdot \color{blue}{\left(\frac{2}{3} + \left(\mathsf{neg}\left(z\right)\right)\right)}, y - x, x\right) \]
            10. +-commutativeN/A

              \[\leadsto \mathsf{fma}\left(6 \cdot \color{blue}{\left(\left(\mathsf{neg}\left(z\right)\right) + \frac{2}{3}\right)}, y - x, x\right) \]
            11. distribute-lft-inN/A

              \[\leadsto \mathsf{fma}\left(\color{blue}{6 \cdot \left(\mathsf{neg}\left(z\right)\right) + 6 \cdot \frac{2}{3}}, y - x, x\right) \]
            12. neg-mul-1N/A

              \[\leadsto \mathsf{fma}\left(6 \cdot \color{blue}{\left(-1 \cdot z\right)} + 6 \cdot \frac{2}{3}, y - x, x\right) \]
            13. associate-*r*N/A

              \[\leadsto \mathsf{fma}\left(\color{blue}{\left(6 \cdot -1\right) \cdot z} + 6 \cdot \frac{2}{3}, y - x, x\right) \]
            14. metadata-evalN/A

              \[\leadsto \mathsf{fma}\left(\color{blue}{-6} \cdot z + 6 \cdot \frac{2}{3}, y - x, x\right) \]
            15. metadata-evalN/A

              \[\leadsto \mathsf{fma}\left(\color{blue}{\left(\mathsf{neg}\left(6\right)\right)} \cdot z + 6 \cdot \frac{2}{3}, y - x, x\right) \]
            16. lower-fma.f64N/A

              \[\leadsto \mathsf{fma}\left(\color{blue}{\mathsf{fma}\left(\mathsf{neg}\left(6\right), z, 6 \cdot \frac{2}{3}\right)}, y - x, x\right) \]
            17. metadata-evalN/A

              \[\leadsto \mathsf{fma}\left(\mathsf{fma}\left(\color{blue}{-6}, z, 6 \cdot \frac{2}{3}\right), y - x, x\right) \]
            18. lift-/.f64N/A

              \[\leadsto \mathsf{fma}\left(\mathsf{fma}\left(-6, z, 6 \cdot \color{blue}{\frac{2}{3}}\right), y - x, x\right) \]
            19. metadata-evalN/A

              \[\leadsto \mathsf{fma}\left(\mathsf{fma}\left(-6, z, 6 \cdot \color{blue}{\frac{2}{3}}\right), y - x, x\right) \]
            20. metadata-eval99.8

              \[\leadsto \mathsf{fma}\left(\mathsf{fma}\left(-6, z, \color{blue}{4}\right), y - x, x\right) \]
          4. Applied rewrites99.8%

            \[\leadsto \color{blue}{\mathsf{fma}\left(\mathsf{fma}\left(-6, z, 4\right), y - x, x\right)} \]
          5. Taylor expanded in z around inf

            \[\leadsto \color{blue}{-6 \cdot \left(z \cdot \left(y - x\right)\right)} \]
          6. Step-by-step derivation
            1. *-commutativeN/A

              \[\leadsto -6 \cdot \color{blue}{\left(\left(y - x\right) \cdot z\right)} \]
            2. associate-*r*N/A

              \[\leadsto \color{blue}{\left(-6 \cdot \left(y - x\right)\right) \cdot z} \]
            3. lower-*.f64N/A

              \[\leadsto \color{blue}{\left(-6 \cdot \left(y - x\right)\right) \cdot z} \]
            4. lower-*.f64N/A

              \[\leadsto \color{blue}{\left(-6 \cdot \left(y - x\right)\right)} \cdot z \]
            5. lower--.f6499.1

              \[\leadsto \left(-6 \cdot \color{blue}{\left(y - x\right)}\right) \cdot z \]
          7. Applied rewrites99.1%

            \[\leadsto \color{blue}{\left(-6 \cdot \left(y - x\right)\right) \cdot z} \]
          8. Taylor expanded in x around inf

            \[\leadsto \left(6 \cdot x\right) \cdot z \]
          9. Step-by-step derivation
            1. Applied rewrites63.6%

              \[\leadsto \left(6 \cdot x\right) \cdot z \]

            if -1e10 < (-.f64 (/.f64 #s(literal 2 binary64) #s(literal 3 binary64)) z) < 1

            1. Initial program 98.8%

              \[x + \left(\left(y - x\right) \cdot 6\right) \cdot \left(\frac{2}{3} - z\right) \]
            2. Add Preprocessing
            3. Taylor expanded in z around 0

              \[\leadsto \color{blue}{x + 4 \cdot \left(y - x\right)} \]
            4. Step-by-step derivation
              1. +-commutativeN/A

                \[\leadsto \color{blue}{4 \cdot \left(y - x\right) + x} \]
              2. lower-fma.f64N/A

                \[\leadsto \color{blue}{\mathsf{fma}\left(4, y - x, x\right)} \]
              3. lower--.f6499.6

                \[\leadsto \mathsf{fma}\left(4, \color{blue}{y - x}, x\right) \]
            5. Applied rewrites99.6%

              \[\leadsto \color{blue}{\mathsf{fma}\left(4, y - x, x\right)} \]
            6. Taylor expanded in x around 0

              \[\leadsto -3 \cdot x + \color{blue}{4 \cdot y} \]
            7. Step-by-step derivation
              1. Applied rewrites99.7%

                \[\leadsto \mathsf{fma}\left(-3, \color{blue}{x}, 4 \cdot y\right) \]

              if 1 < (-.f64 (/.f64 #s(literal 2 binary64) #s(literal 3 binary64)) z) < 9.9999999999999998e155

              1. Initial program 99.3%

                \[x + \left(\left(y - x\right) \cdot 6\right) \cdot \left(\frac{2}{3} - z\right) \]
              2. Add Preprocessing
              3. Step-by-step derivation
                1. lift-+.f64N/A

                  \[\leadsto \color{blue}{x + \left(\left(y - x\right) \cdot 6\right) \cdot \left(\frac{2}{3} - z\right)} \]
                2. +-commutativeN/A

                  \[\leadsto \color{blue}{\left(\left(y - x\right) \cdot 6\right) \cdot \left(\frac{2}{3} - z\right) + x} \]
                3. lift-*.f64N/A

                  \[\leadsto \color{blue}{\left(\left(y - x\right) \cdot 6\right) \cdot \left(\frac{2}{3} - z\right)} + x \]
                4. lift-*.f64N/A

                  \[\leadsto \color{blue}{\left(\left(y - x\right) \cdot 6\right)} \cdot \left(\frac{2}{3} - z\right) + x \]
                5. associate-*l*N/A

                  \[\leadsto \color{blue}{\left(y - x\right) \cdot \left(6 \cdot \left(\frac{2}{3} - z\right)\right)} + x \]
                6. *-commutativeN/A

                  \[\leadsto \color{blue}{\left(6 \cdot \left(\frac{2}{3} - z\right)\right) \cdot \left(y - x\right)} + x \]
                7. lower-fma.f64N/A

                  \[\leadsto \color{blue}{\mathsf{fma}\left(6 \cdot \left(\frac{2}{3} - z\right), y - x, x\right)} \]
                8. lift--.f64N/A

                  \[\leadsto \mathsf{fma}\left(6 \cdot \color{blue}{\left(\frac{2}{3} - z\right)}, y - x, x\right) \]
                9. sub-negN/A

                  \[\leadsto \mathsf{fma}\left(6 \cdot \color{blue}{\left(\frac{2}{3} + \left(\mathsf{neg}\left(z\right)\right)\right)}, y - x, x\right) \]
                10. +-commutativeN/A

                  \[\leadsto \mathsf{fma}\left(6 \cdot \color{blue}{\left(\left(\mathsf{neg}\left(z\right)\right) + \frac{2}{3}\right)}, y - x, x\right) \]
                11. distribute-lft-inN/A

                  \[\leadsto \mathsf{fma}\left(\color{blue}{6 \cdot \left(\mathsf{neg}\left(z\right)\right) + 6 \cdot \frac{2}{3}}, y - x, x\right) \]
                12. neg-mul-1N/A

                  \[\leadsto \mathsf{fma}\left(6 \cdot \color{blue}{\left(-1 \cdot z\right)} + 6 \cdot \frac{2}{3}, y - x, x\right) \]
                13. associate-*r*N/A

                  \[\leadsto \mathsf{fma}\left(\color{blue}{\left(6 \cdot -1\right) \cdot z} + 6 \cdot \frac{2}{3}, y - x, x\right) \]
                14. metadata-evalN/A

                  \[\leadsto \mathsf{fma}\left(\color{blue}{-6} \cdot z + 6 \cdot \frac{2}{3}, y - x, x\right) \]
                15. metadata-evalN/A

                  \[\leadsto \mathsf{fma}\left(\color{blue}{\left(\mathsf{neg}\left(6\right)\right)} \cdot z + 6 \cdot \frac{2}{3}, y - x, x\right) \]
                16. lower-fma.f64N/A

                  \[\leadsto \mathsf{fma}\left(\color{blue}{\mathsf{fma}\left(\mathsf{neg}\left(6\right), z, 6 \cdot \frac{2}{3}\right)}, y - x, x\right) \]
                17. metadata-evalN/A

                  \[\leadsto \mathsf{fma}\left(\mathsf{fma}\left(\color{blue}{-6}, z, 6 \cdot \frac{2}{3}\right), y - x, x\right) \]
                18. lift-/.f64N/A

                  \[\leadsto \mathsf{fma}\left(\mathsf{fma}\left(-6, z, 6 \cdot \color{blue}{\frac{2}{3}}\right), y - x, x\right) \]
                19. metadata-evalN/A

                  \[\leadsto \mathsf{fma}\left(\mathsf{fma}\left(-6, z, 6 \cdot \color{blue}{\frac{2}{3}}\right), y - x, x\right) \]
                20. metadata-eval99.6

                  \[\leadsto \mathsf{fma}\left(\mathsf{fma}\left(-6, z, \color{blue}{4}\right), y - x, x\right) \]
              4. Applied rewrites99.6%

                \[\leadsto \color{blue}{\mathsf{fma}\left(\mathsf{fma}\left(-6, z, 4\right), y - x, x\right)} \]
              5. Taylor expanded in x around 0

                \[\leadsto \color{blue}{y \cdot \left(4 + -6 \cdot z\right)} \]
              6. Step-by-step derivation
                1. *-commutativeN/A

                  \[\leadsto \color{blue}{\left(4 + -6 \cdot z\right) \cdot y} \]
                2. lower-*.f64N/A

                  \[\leadsto \color{blue}{\left(4 + -6 \cdot z\right) \cdot y} \]
                3. +-commutativeN/A

                  \[\leadsto \color{blue}{\left(-6 \cdot z + 4\right)} \cdot y \]
                4. lower-fma.f6460.5

                  \[\leadsto \color{blue}{\mathsf{fma}\left(-6, z, 4\right)} \cdot y \]
              7. Applied rewrites60.5%

                \[\leadsto \color{blue}{\mathsf{fma}\left(-6, z, 4\right) \cdot y} \]
              8. Taylor expanded in z around inf

                \[\leadsto \left(-6 \cdot z\right) \cdot y \]
              9. Step-by-step derivation
                1. Applied rewrites58.9%

                  \[\leadsto \left(-6 \cdot z\right) \cdot y \]
              10. Recombined 3 regimes into one program.
              11. Final simplification83.3%

                \[\leadsto \begin{array}{l} \mathbf{if}\;\frac{2}{3} - z \leq -10000000000:\\ \;\;\;\;\left(6 \cdot x\right) \cdot z\\ \mathbf{elif}\;\frac{2}{3} - z \leq 1:\\ \;\;\;\;\mathsf{fma}\left(-3, x, y \cdot 4\right)\\ \mathbf{elif}\;\frac{2}{3} - z \leq 10^{+156}:\\ \;\;\;\;\left(z \cdot -6\right) \cdot y\\ \mathbf{else}:\\ \;\;\;\;\left(6 \cdot x\right) \cdot z\\ \end{array} \]
              12. Add Preprocessing

              Alternative 5: 74.4% accurate, 0.4× speedup?

              \[\begin{array}{l} \\ \begin{array}{l} t_0 := \frac{2}{3} - z\\ t_1 := \left(6 \cdot x\right) \cdot z\\ \mathbf{if}\;t\_0 \leq -10000000000:\\ \;\;\;\;t\_1\\ \mathbf{elif}\;t\_0 \leq 1:\\ \;\;\;\;\mathsf{fma}\left(4, y - x, x\right)\\ \mathbf{elif}\;t\_0 \leq 10^{+156}:\\ \;\;\;\;\left(z \cdot -6\right) \cdot y\\ \mathbf{else}:\\ \;\;\;\;t\_1\\ \end{array} \end{array} \]
              (FPCore (x y z)
               :precision binary64
               (let* ((t_0 (- (/ 2.0 3.0) z)) (t_1 (* (* 6.0 x) z)))
                 (if (<= t_0 -10000000000.0)
                   t_1
                   (if (<= t_0 1.0)
                     (fma 4.0 (- y x) x)
                     (if (<= t_0 1e+156) (* (* z -6.0) y) t_1)))))
              double code(double x, double y, double z) {
              	double t_0 = (2.0 / 3.0) - z;
              	double t_1 = (6.0 * x) * z;
              	double tmp;
              	if (t_0 <= -10000000000.0) {
              		tmp = t_1;
              	} else if (t_0 <= 1.0) {
              		tmp = fma(4.0, (y - x), x);
              	} else if (t_0 <= 1e+156) {
              		tmp = (z * -6.0) * y;
              	} else {
              		tmp = t_1;
              	}
              	return tmp;
              }
              
              function code(x, y, z)
              	t_0 = Float64(Float64(2.0 / 3.0) - z)
              	t_1 = Float64(Float64(6.0 * x) * z)
              	tmp = 0.0
              	if (t_0 <= -10000000000.0)
              		tmp = t_1;
              	elseif (t_0 <= 1.0)
              		tmp = fma(4.0, Float64(y - x), x);
              	elseif (t_0 <= 1e+156)
              		tmp = Float64(Float64(z * -6.0) * y);
              	else
              		tmp = t_1;
              	end
              	return tmp
              end
              
              code[x_, y_, z_] := Block[{t$95$0 = N[(N[(2.0 / 3.0), $MachinePrecision] - z), $MachinePrecision]}, Block[{t$95$1 = N[(N[(6.0 * x), $MachinePrecision] * z), $MachinePrecision]}, If[LessEqual[t$95$0, -10000000000.0], t$95$1, If[LessEqual[t$95$0, 1.0], N[(4.0 * N[(y - x), $MachinePrecision] + x), $MachinePrecision], If[LessEqual[t$95$0, 1e+156], N[(N[(z * -6.0), $MachinePrecision] * y), $MachinePrecision], t$95$1]]]]]
              
              \begin{array}{l}
              
              \\
              \begin{array}{l}
              t_0 := \frac{2}{3} - z\\
              t_1 := \left(6 \cdot x\right) \cdot z\\
              \mathbf{if}\;t\_0 \leq -10000000000:\\
              \;\;\;\;t\_1\\
              
              \mathbf{elif}\;t\_0 \leq 1:\\
              \;\;\;\;\mathsf{fma}\left(4, y - x, x\right)\\
              
              \mathbf{elif}\;t\_0 \leq 10^{+156}:\\
              \;\;\;\;\left(z \cdot -6\right) \cdot y\\
              
              \mathbf{else}:\\
              \;\;\;\;t\_1\\
              
              
              \end{array}
              \end{array}
              
              Derivation
              1. Split input into 3 regimes
              2. if (-.f64 (/.f64 #s(literal 2 binary64) #s(literal 3 binary64)) z) < -1e10 or 9.9999999999999998e155 < (-.f64 (/.f64 #s(literal 2 binary64) #s(literal 3 binary64)) z)

                1. Initial program 99.7%

                  \[x + \left(\left(y - x\right) \cdot 6\right) \cdot \left(\frac{2}{3} - z\right) \]
                2. Add Preprocessing
                3. Step-by-step derivation
                  1. lift-+.f64N/A

                    \[\leadsto \color{blue}{x + \left(\left(y - x\right) \cdot 6\right) \cdot \left(\frac{2}{3} - z\right)} \]
                  2. +-commutativeN/A

                    \[\leadsto \color{blue}{\left(\left(y - x\right) \cdot 6\right) \cdot \left(\frac{2}{3} - z\right) + x} \]
                  3. lift-*.f64N/A

                    \[\leadsto \color{blue}{\left(\left(y - x\right) \cdot 6\right) \cdot \left(\frac{2}{3} - z\right)} + x \]
                  4. lift-*.f64N/A

                    \[\leadsto \color{blue}{\left(\left(y - x\right) \cdot 6\right)} \cdot \left(\frac{2}{3} - z\right) + x \]
                  5. associate-*l*N/A

                    \[\leadsto \color{blue}{\left(y - x\right) \cdot \left(6 \cdot \left(\frac{2}{3} - z\right)\right)} + x \]
                  6. *-commutativeN/A

                    \[\leadsto \color{blue}{\left(6 \cdot \left(\frac{2}{3} - z\right)\right) \cdot \left(y - x\right)} + x \]
                  7. lower-fma.f64N/A

                    \[\leadsto \color{blue}{\mathsf{fma}\left(6 \cdot \left(\frac{2}{3} - z\right), y - x, x\right)} \]
                  8. lift--.f64N/A

                    \[\leadsto \mathsf{fma}\left(6 \cdot \color{blue}{\left(\frac{2}{3} - z\right)}, y - x, x\right) \]
                  9. sub-negN/A

                    \[\leadsto \mathsf{fma}\left(6 \cdot \color{blue}{\left(\frac{2}{3} + \left(\mathsf{neg}\left(z\right)\right)\right)}, y - x, x\right) \]
                  10. +-commutativeN/A

                    \[\leadsto \mathsf{fma}\left(6 \cdot \color{blue}{\left(\left(\mathsf{neg}\left(z\right)\right) + \frac{2}{3}\right)}, y - x, x\right) \]
                  11. distribute-lft-inN/A

                    \[\leadsto \mathsf{fma}\left(\color{blue}{6 \cdot \left(\mathsf{neg}\left(z\right)\right) + 6 \cdot \frac{2}{3}}, y - x, x\right) \]
                  12. neg-mul-1N/A

                    \[\leadsto \mathsf{fma}\left(6 \cdot \color{blue}{\left(-1 \cdot z\right)} + 6 \cdot \frac{2}{3}, y - x, x\right) \]
                  13. associate-*r*N/A

                    \[\leadsto \mathsf{fma}\left(\color{blue}{\left(6 \cdot -1\right) \cdot z} + 6 \cdot \frac{2}{3}, y - x, x\right) \]
                  14. metadata-evalN/A

                    \[\leadsto \mathsf{fma}\left(\color{blue}{-6} \cdot z + 6 \cdot \frac{2}{3}, y - x, x\right) \]
                  15. metadata-evalN/A

                    \[\leadsto \mathsf{fma}\left(\color{blue}{\left(\mathsf{neg}\left(6\right)\right)} \cdot z + 6 \cdot \frac{2}{3}, y - x, x\right) \]
                  16. lower-fma.f64N/A

                    \[\leadsto \mathsf{fma}\left(\color{blue}{\mathsf{fma}\left(\mathsf{neg}\left(6\right), z, 6 \cdot \frac{2}{3}\right)}, y - x, x\right) \]
                  17. metadata-evalN/A

                    \[\leadsto \mathsf{fma}\left(\mathsf{fma}\left(\color{blue}{-6}, z, 6 \cdot \frac{2}{3}\right), y - x, x\right) \]
                  18. lift-/.f64N/A

                    \[\leadsto \mathsf{fma}\left(\mathsf{fma}\left(-6, z, 6 \cdot \color{blue}{\frac{2}{3}}\right), y - x, x\right) \]
                  19. metadata-evalN/A

                    \[\leadsto \mathsf{fma}\left(\mathsf{fma}\left(-6, z, 6 \cdot \color{blue}{\frac{2}{3}}\right), y - x, x\right) \]
                  20. metadata-eval99.8

                    \[\leadsto \mathsf{fma}\left(\mathsf{fma}\left(-6, z, \color{blue}{4}\right), y - x, x\right) \]
                4. Applied rewrites99.8%

                  \[\leadsto \color{blue}{\mathsf{fma}\left(\mathsf{fma}\left(-6, z, 4\right), y - x, x\right)} \]
                5. Taylor expanded in z around inf

                  \[\leadsto \color{blue}{-6 \cdot \left(z \cdot \left(y - x\right)\right)} \]
                6. Step-by-step derivation
                  1. *-commutativeN/A

                    \[\leadsto -6 \cdot \color{blue}{\left(\left(y - x\right) \cdot z\right)} \]
                  2. associate-*r*N/A

                    \[\leadsto \color{blue}{\left(-6 \cdot \left(y - x\right)\right) \cdot z} \]
                  3. lower-*.f64N/A

                    \[\leadsto \color{blue}{\left(-6 \cdot \left(y - x\right)\right) \cdot z} \]
                  4. lower-*.f64N/A

                    \[\leadsto \color{blue}{\left(-6 \cdot \left(y - x\right)\right)} \cdot z \]
                  5. lower--.f6499.1

                    \[\leadsto \left(-6 \cdot \color{blue}{\left(y - x\right)}\right) \cdot z \]
                7. Applied rewrites99.1%

                  \[\leadsto \color{blue}{\left(-6 \cdot \left(y - x\right)\right) \cdot z} \]
                8. Taylor expanded in x around inf

                  \[\leadsto \left(6 \cdot x\right) \cdot z \]
                9. Step-by-step derivation
                  1. Applied rewrites63.6%

                    \[\leadsto \left(6 \cdot x\right) \cdot z \]

                  if -1e10 < (-.f64 (/.f64 #s(literal 2 binary64) #s(literal 3 binary64)) z) < 1

                  1. Initial program 98.8%

                    \[x + \left(\left(y - x\right) \cdot 6\right) \cdot \left(\frac{2}{3} - z\right) \]
                  2. Add Preprocessing
                  3. Taylor expanded in z around 0

                    \[\leadsto \color{blue}{x + 4 \cdot \left(y - x\right)} \]
                  4. Step-by-step derivation
                    1. +-commutativeN/A

                      \[\leadsto \color{blue}{4 \cdot \left(y - x\right) + x} \]
                    2. lower-fma.f64N/A

                      \[\leadsto \color{blue}{\mathsf{fma}\left(4, y - x, x\right)} \]
                    3. lower--.f6499.6

                      \[\leadsto \mathsf{fma}\left(4, \color{blue}{y - x}, x\right) \]
                  5. Applied rewrites99.6%

                    \[\leadsto \color{blue}{\mathsf{fma}\left(4, y - x, x\right)} \]

                  if 1 < (-.f64 (/.f64 #s(literal 2 binary64) #s(literal 3 binary64)) z) < 9.9999999999999998e155

                  1. Initial program 99.3%

                    \[x + \left(\left(y - x\right) \cdot 6\right) \cdot \left(\frac{2}{3} - z\right) \]
                  2. Add Preprocessing
                  3. Step-by-step derivation
                    1. lift-+.f64N/A

                      \[\leadsto \color{blue}{x + \left(\left(y - x\right) \cdot 6\right) \cdot \left(\frac{2}{3} - z\right)} \]
                    2. +-commutativeN/A

                      \[\leadsto \color{blue}{\left(\left(y - x\right) \cdot 6\right) \cdot \left(\frac{2}{3} - z\right) + x} \]
                    3. lift-*.f64N/A

                      \[\leadsto \color{blue}{\left(\left(y - x\right) \cdot 6\right) \cdot \left(\frac{2}{3} - z\right)} + x \]
                    4. lift-*.f64N/A

                      \[\leadsto \color{blue}{\left(\left(y - x\right) \cdot 6\right)} \cdot \left(\frac{2}{3} - z\right) + x \]
                    5. associate-*l*N/A

                      \[\leadsto \color{blue}{\left(y - x\right) \cdot \left(6 \cdot \left(\frac{2}{3} - z\right)\right)} + x \]
                    6. *-commutativeN/A

                      \[\leadsto \color{blue}{\left(6 \cdot \left(\frac{2}{3} - z\right)\right) \cdot \left(y - x\right)} + x \]
                    7. lower-fma.f64N/A

                      \[\leadsto \color{blue}{\mathsf{fma}\left(6 \cdot \left(\frac{2}{3} - z\right), y - x, x\right)} \]
                    8. lift--.f64N/A

                      \[\leadsto \mathsf{fma}\left(6 \cdot \color{blue}{\left(\frac{2}{3} - z\right)}, y - x, x\right) \]
                    9. sub-negN/A

                      \[\leadsto \mathsf{fma}\left(6 \cdot \color{blue}{\left(\frac{2}{3} + \left(\mathsf{neg}\left(z\right)\right)\right)}, y - x, x\right) \]
                    10. +-commutativeN/A

                      \[\leadsto \mathsf{fma}\left(6 \cdot \color{blue}{\left(\left(\mathsf{neg}\left(z\right)\right) + \frac{2}{3}\right)}, y - x, x\right) \]
                    11. distribute-lft-inN/A

                      \[\leadsto \mathsf{fma}\left(\color{blue}{6 \cdot \left(\mathsf{neg}\left(z\right)\right) + 6 \cdot \frac{2}{3}}, y - x, x\right) \]
                    12. neg-mul-1N/A

                      \[\leadsto \mathsf{fma}\left(6 \cdot \color{blue}{\left(-1 \cdot z\right)} + 6 \cdot \frac{2}{3}, y - x, x\right) \]
                    13. associate-*r*N/A

                      \[\leadsto \mathsf{fma}\left(\color{blue}{\left(6 \cdot -1\right) \cdot z} + 6 \cdot \frac{2}{3}, y - x, x\right) \]
                    14. metadata-evalN/A

                      \[\leadsto \mathsf{fma}\left(\color{blue}{-6} \cdot z + 6 \cdot \frac{2}{3}, y - x, x\right) \]
                    15. metadata-evalN/A

                      \[\leadsto \mathsf{fma}\left(\color{blue}{\left(\mathsf{neg}\left(6\right)\right)} \cdot z + 6 \cdot \frac{2}{3}, y - x, x\right) \]
                    16. lower-fma.f64N/A

                      \[\leadsto \mathsf{fma}\left(\color{blue}{\mathsf{fma}\left(\mathsf{neg}\left(6\right), z, 6 \cdot \frac{2}{3}\right)}, y - x, x\right) \]
                    17. metadata-evalN/A

                      \[\leadsto \mathsf{fma}\left(\mathsf{fma}\left(\color{blue}{-6}, z, 6 \cdot \frac{2}{3}\right), y - x, x\right) \]
                    18. lift-/.f64N/A

                      \[\leadsto \mathsf{fma}\left(\mathsf{fma}\left(-6, z, 6 \cdot \color{blue}{\frac{2}{3}}\right), y - x, x\right) \]
                    19. metadata-evalN/A

                      \[\leadsto \mathsf{fma}\left(\mathsf{fma}\left(-6, z, 6 \cdot \color{blue}{\frac{2}{3}}\right), y - x, x\right) \]
                    20. metadata-eval99.6

                      \[\leadsto \mathsf{fma}\left(\mathsf{fma}\left(-6, z, \color{blue}{4}\right), y - x, x\right) \]
                  4. Applied rewrites99.6%

                    \[\leadsto \color{blue}{\mathsf{fma}\left(\mathsf{fma}\left(-6, z, 4\right), y - x, x\right)} \]
                  5. Taylor expanded in x around 0

                    \[\leadsto \color{blue}{y \cdot \left(4 + -6 \cdot z\right)} \]
                  6. Step-by-step derivation
                    1. *-commutativeN/A

                      \[\leadsto \color{blue}{\left(4 + -6 \cdot z\right) \cdot y} \]
                    2. lower-*.f64N/A

                      \[\leadsto \color{blue}{\left(4 + -6 \cdot z\right) \cdot y} \]
                    3. +-commutativeN/A

                      \[\leadsto \color{blue}{\left(-6 \cdot z + 4\right)} \cdot y \]
                    4. lower-fma.f6460.5

                      \[\leadsto \color{blue}{\mathsf{fma}\left(-6, z, 4\right)} \cdot y \]
                  7. Applied rewrites60.5%

                    \[\leadsto \color{blue}{\mathsf{fma}\left(-6, z, 4\right) \cdot y} \]
                  8. Taylor expanded in z around inf

                    \[\leadsto \left(-6 \cdot z\right) \cdot y \]
                  9. Step-by-step derivation
                    1. Applied rewrites58.9%

                      \[\leadsto \left(-6 \cdot z\right) \cdot y \]
                  10. Recombined 3 regimes into one program.
                  11. Final simplification83.2%

                    \[\leadsto \begin{array}{l} \mathbf{if}\;\frac{2}{3} - z \leq -10000000000:\\ \;\;\;\;\left(6 \cdot x\right) \cdot z\\ \mathbf{elif}\;\frac{2}{3} - z \leq 1:\\ \;\;\;\;\mathsf{fma}\left(4, y - x, x\right)\\ \mathbf{elif}\;\frac{2}{3} - z \leq 10^{+156}:\\ \;\;\;\;\left(z \cdot -6\right) \cdot y\\ \mathbf{else}:\\ \;\;\;\;\left(6 \cdot x\right) \cdot z\\ \end{array} \]
                  12. Add Preprocessing

                  Alternative 6: 74.4% accurate, 0.4× speedup?

                  \[\begin{array}{l} \\ \begin{array}{l} t_0 := \frac{2}{3} - z\\ t_1 := \left(6 \cdot x\right) \cdot z\\ \mathbf{if}\;t\_0 \leq -10000000000:\\ \;\;\;\;t\_1\\ \mathbf{elif}\;t\_0 \leq 1:\\ \;\;\;\;\mathsf{fma}\left(4, y - x, x\right)\\ \mathbf{elif}\;t\_0 \leq 10^{+156}:\\ \;\;\;\;\left(y \cdot -6\right) \cdot z\\ \mathbf{else}:\\ \;\;\;\;t\_1\\ \end{array} \end{array} \]
                  (FPCore (x y z)
                   :precision binary64
                   (let* ((t_0 (- (/ 2.0 3.0) z)) (t_1 (* (* 6.0 x) z)))
                     (if (<= t_0 -10000000000.0)
                       t_1
                       (if (<= t_0 1.0)
                         (fma 4.0 (- y x) x)
                         (if (<= t_0 1e+156) (* (* y -6.0) z) t_1)))))
                  double code(double x, double y, double z) {
                  	double t_0 = (2.0 / 3.0) - z;
                  	double t_1 = (6.0 * x) * z;
                  	double tmp;
                  	if (t_0 <= -10000000000.0) {
                  		tmp = t_1;
                  	} else if (t_0 <= 1.0) {
                  		tmp = fma(4.0, (y - x), x);
                  	} else if (t_0 <= 1e+156) {
                  		tmp = (y * -6.0) * z;
                  	} else {
                  		tmp = t_1;
                  	}
                  	return tmp;
                  }
                  
                  function code(x, y, z)
                  	t_0 = Float64(Float64(2.0 / 3.0) - z)
                  	t_1 = Float64(Float64(6.0 * x) * z)
                  	tmp = 0.0
                  	if (t_0 <= -10000000000.0)
                  		tmp = t_1;
                  	elseif (t_0 <= 1.0)
                  		tmp = fma(4.0, Float64(y - x), x);
                  	elseif (t_0 <= 1e+156)
                  		tmp = Float64(Float64(y * -6.0) * z);
                  	else
                  		tmp = t_1;
                  	end
                  	return tmp
                  end
                  
                  code[x_, y_, z_] := Block[{t$95$0 = N[(N[(2.0 / 3.0), $MachinePrecision] - z), $MachinePrecision]}, Block[{t$95$1 = N[(N[(6.0 * x), $MachinePrecision] * z), $MachinePrecision]}, If[LessEqual[t$95$0, -10000000000.0], t$95$1, If[LessEqual[t$95$0, 1.0], N[(4.0 * N[(y - x), $MachinePrecision] + x), $MachinePrecision], If[LessEqual[t$95$0, 1e+156], N[(N[(y * -6.0), $MachinePrecision] * z), $MachinePrecision], t$95$1]]]]]
                  
                  \begin{array}{l}
                  
                  \\
                  \begin{array}{l}
                  t_0 := \frac{2}{3} - z\\
                  t_1 := \left(6 \cdot x\right) \cdot z\\
                  \mathbf{if}\;t\_0 \leq -10000000000:\\
                  \;\;\;\;t\_1\\
                  
                  \mathbf{elif}\;t\_0 \leq 1:\\
                  \;\;\;\;\mathsf{fma}\left(4, y - x, x\right)\\
                  
                  \mathbf{elif}\;t\_0 \leq 10^{+156}:\\
                  \;\;\;\;\left(y \cdot -6\right) \cdot z\\
                  
                  \mathbf{else}:\\
                  \;\;\;\;t\_1\\
                  
                  
                  \end{array}
                  \end{array}
                  
                  Derivation
                  1. Split input into 3 regimes
                  2. if (-.f64 (/.f64 #s(literal 2 binary64) #s(literal 3 binary64)) z) < -1e10 or 9.9999999999999998e155 < (-.f64 (/.f64 #s(literal 2 binary64) #s(literal 3 binary64)) z)

                    1. Initial program 99.7%

                      \[x + \left(\left(y - x\right) \cdot 6\right) \cdot \left(\frac{2}{3} - z\right) \]
                    2. Add Preprocessing
                    3. Step-by-step derivation
                      1. lift-+.f64N/A

                        \[\leadsto \color{blue}{x + \left(\left(y - x\right) \cdot 6\right) \cdot \left(\frac{2}{3} - z\right)} \]
                      2. +-commutativeN/A

                        \[\leadsto \color{blue}{\left(\left(y - x\right) \cdot 6\right) \cdot \left(\frac{2}{3} - z\right) + x} \]
                      3. lift-*.f64N/A

                        \[\leadsto \color{blue}{\left(\left(y - x\right) \cdot 6\right) \cdot \left(\frac{2}{3} - z\right)} + x \]
                      4. lift-*.f64N/A

                        \[\leadsto \color{blue}{\left(\left(y - x\right) \cdot 6\right)} \cdot \left(\frac{2}{3} - z\right) + x \]
                      5. associate-*l*N/A

                        \[\leadsto \color{blue}{\left(y - x\right) \cdot \left(6 \cdot \left(\frac{2}{3} - z\right)\right)} + x \]
                      6. *-commutativeN/A

                        \[\leadsto \color{blue}{\left(6 \cdot \left(\frac{2}{3} - z\right)\right) \cdot \left(y - x\right)} + x \]
                      7. lower-fma.f64N/A

                        \[\leadsto \color{blue}{\mathsf{fma}\left(6 \cdot \left(\frac{2}{3} - z\right), y - x, x\right)} \]
                      8. lift--.f64N/A

                        \[\leadsto \mathsf{fma}\left(6 \cdot \color{blue}{\left(\frac{2}{3} - z\right)}, y - x, x\right) \]
                      9. sub-negN/A

                        \[\leadsto \mathsf{fma}\left(6 \cdot \color{blue}{\left(\frac{2}{3} + \left(\mathsf{neg}\left(z\right)\right)\right)}, y - x, x\right) \]
                      10. +-commutativeN/A

                        \[\leadsto \mathsf{fma}\left(6 \cdot \color{blue}{\left(\left(\mathsf{neg}\left(z\right)\right) + \frac{2}{3}\right)}, y - x, x\right) \]
                      11. distribute-lft-inN/A

                        \[\leadsto \mathsf{fma}\left(\color{blue}{6 \cdot \left(\mathsf{neg}\left(z\right)\right) + 6 \cdot \frac{2}{3}}, y - x, x\right) \]
                      12. neg-mul-1N/A

                        \[\leadsto \mathsf{fma}\left(6 \cdot \color{blue}{\left(-1 \cdot z\right)} + 6 \cdot \frac{2}{3}, y - x, x\right) \]
                      13. associate-*r*N/A

                        \[\leadsto \mathsf{fma}\left(\color{blue}{\left(6 \cdot -1\right) \cdot z} + 6 \cdot \frac{2}{3}, y - x, x\right) \]
                      14. metadata-evalN/A

                        \[\leadsto \mathsf{fma}\left(\color{blue}{-6} \cdot z + 6 \cdot \frac{2}{3}, y - x, x\right) \]
                      15. metadata-evalN/A

                        \[\leadsto \mathsf{fma}\left(\color{blue}{\left(\mathsf{neg}\left(6\right)\right)} \cdot z + 6 \cdot \frac{2}{3}, y - x, x\right) \]
                      16. lower-fma.f64N/A

                        \[\leadsto \mathsf{fma}\left(\color{blue}{\mathsf{fma}\left(\mathsf{neg}\left(6\right), z, 6 \cdot \frac{2}{3}\right)}, y - x, x\right) \]
                      17. metadata-evalN/A

                        \[\leadsto \mathsf{fma}\left(\mathsf{fma}\left(\color{blue}{-6}, z, 6 \cdot \frac{2}{3}\right), y - x, x\right) \]
                      18. lift-/.f64N/A

                        \[\leadsto \mathsf{fma}\left(\mathsf{fma}\left(-6, z, 6 \cdot \color{blue}{\frac{2}{3}}\right), y - x, x\right) \]
                      19. metadata-evalN/A

                        \[\leadsto \mathsf{fma}\left(\mathsf{fma}\left(-6, z, 6 \cdot \color{blue}{\frac{2}{3}}\right), y - x, x\right) \]
                      20. metadata-eval99.8

                        \[\leadsto \mathsf{fma}\left(\mathsf{fma}\left(-6, z, \color{blue}{4}\right), y - x, x\right) \]
                    4. Applied rewrites99.8%

                      \[\leadsto \color{blue}{\mathsf{fma}\left(\mathsf{fma}\left(-6, z, 4\right), y - x, x\right)} \]
                    5. Taylor expanded in z around inf

                      \[\leadsto \color{blue}{-6 \cdot \left(z \cdot \left(y - x\right)\right)} \]
                    6. Step-by-step derivation
                      1. *-commutativeN/A

                        \[\leadsto -6 \cdot \color{blue}{\left(\left(y - x\right) \cdot z\right)} \]
                      2. associate-*r*N/A

                        \[\leadsto \color{blue}{\left(-6 \cdot \left(y - x\right)\right) \cdot z} \]
                      3. lower-*.f64N/A

                        \[\leadsto \color{blue}{\left(-6 \cdot \left(y - x\right)\right) \cdot z} \]
                      4. lower-*.f64N/A

                        \[\leadsto \color{blue}{\left(-6 \cdot \left(y - x\right)\right)} \cdot z \]
                      5. lower--.f6499.1

                        \[\leadsto \left(-6 \cdot \color{blue}{\left(y - x\right)}\right) \cdot z \]
                    7. Applied rewrites99.1%

                      \[\leadsto \color{blue}{\left(-6 \cdot \left(y - x\right)\right) \cdot z} \]
                    8. Taylor expanded in x around inf

                      \[\leadsto \left(6 \cdot x\right) \cdot z \]
                    9. Step-by-step derivation
                      1. Applied rewrites63.6%

                        \[\leadsto \left(6 \cdot x\right) \cdot z \]

                      if -1e10 < (-.f64 (/.f64 #s(literal 2 binary64) #s(literal 3 binary64)) z) < 1

                      1. Initial program 98.8%

                        \[x + \left(\left(y - x\right) \cdot 6\right) \cdot \left(\frac{2}{3} - z\right) \]
                      2. Add Preprocessing
                      3. Taylor expanded in z around 0

                        \[\leadsto \color{blue}{x + 4 \cdot \left(y - x\right)} \]
                      4. Step-by-step derivation
                        1. +-commutativeN/A

                          \[\leadsto \color{blue}{4 \cdot \left(y - x\right) + x} \]
                        2. lower-fma.f64N/A

                          \[\leadsto \color{blue}{\mathsf{fma}\left(4, y - x, x\right)} \]
                        3. lower--.f6499.6

                          \[\leadsto \mathsf{fma}\left(4, \color{blue}{y - x}, x\right) \]
                      5. Applied rewrites99.6%

                        \[\leadsto \color{blue}{\mathsf{fma}\left(4, y - x, x\right)} \]

                      if 1 < (-.f64 (/.f64 #s(literal 2 binary64) #s(literal 3 binary64)) z) < 9.9999999999999998e155

                      1. Initial program 99.3%

                        \[x + \left(\left(y - x\right) \cdot 6\right) \cdot \left(\frac{2}{3} - z\right) \]
                      2. Add Preprocessing
                      3. Step-by-step derivation
                        1. lift-+.f64N/A

                          \[\leadsto \color{blue}{x + \left(\left(y - x\right) \cdot 6\right) \cdot \left(\frac{2}{3} - z\right)} \]
                        2. +-commutativeN/A

                          \[\leadsto \color{blue}{\left(\left(y - x\right) \cdot 6\right) \cdot \left(\frac{2}{3} - z\right) + x} \]
                        3. lift-*.f64N/A

                          \[\leadsto \color{blue}{\left(\left(y - x\right) \cdot 6\right) \cdot \left(\frac{2}{3} - z\right)} + x \]
                        4. lift-*.f64N/A

                          \[\leadsto \color{blue}{\left(\left(y - x\right) \cdot 6\right)} \cdot \left(\frac{2}{3} - z\right) + x \]
                        5. associate-*l*N/A

                          \[\leadsto \color{blue}{\left(y - x\right) \cdot \left(6 \cdot \left(\frac{2}{3} - z\right)\right)} + x \]
                        6. *-commutativeN/A

                          \[\leadsto \color{blue}{\left(6 \cdot \left(\frac{2}{3} - z\right)\right) \cdot \left(y - x\right)} + x \]
                        7. lower-fma.f64N/A

                          \[\leadsto \color{blue}{\mathsf{fma}\left(6 \cdot \left(\frac{2}{3} - z\right), y - x, x\right)} \]
                        8. lift--.f64N/A

                          \[\leadsto \mathsf{fma}\left(6 \cdot \color{blue}{\left(\frac{2}{3} - z\right)}, y - x, x\right) \]
                        9. sub-negN/A

                          \[\leadsto \mathsf{fma}\left(6 \cdot \color{blue}{\left(\frac{2}{3} + \left(\mathsf{neg}\left(z\right)\right)\right)}, y - x, x\right) \]
                        10. +-commutativeN/A

                          \[\leadsto \mathsf{fma}\left(6 \cdot \color{blue}{\left(\left(\mathsf{neg}\left(z\right)\right) + \frac{2}{3}\right)}, y - x, x\right) \]
                        11. distribute-lft-inN/A

                          \[\leadsto \mathsf{fma}\left(\color{blue}{6 \cdot \left(\mathsf{neg}\left(z\right)\right) + 6 \cdot \frac{2}{3}}, y - x, x\right) \]
                        12. neg-mul-1N/A

                          \[\leadsto \mathsf{fma}\left(6 \cdot \color{blue}{\left(-1 \cdot z\right)} + 6 \cdot \frac{2}{3}, y - x, x\right) \]
                        13. associate-*r*N/A

                          \[\leadsto \mathsf{fma}\left(\color{blue}{\left(6 \cdot -1\right) \cdot z} + 6 \cdot \frac{2}{3}, y - x, x\right) \]
                        14. metadata-evalN/A

                          \[\leadsto \mathsf{fma}\left(\color{blue}{-6} \cdot z + 6 \cdot \frac{2}{3}, y - x, x\right) \]
                        15. metadata-evalN/A

                          \[\leadsto \mathsf{fma}\left(\color{blue}{\left(\mathsf{neg}\left(6\right)\right)} \cdot z + 6 \cdot \frac{2}{3}, y - x, x\right) \]
                        16. lower-fma.f64N/A

                          \[\leadsto \mathsf{fma}\left(\color{blue}{\mathsf{fma}\left(\mathsf{neg}\left(6\right), z, 6 \cdot \frac{2}{3}\right)}, y - x, x\right) \]
                        17. metadata-evalN/A

                          \[\leadsto \mathsf{fma}\left(\mathsf{fma}\left(\color{blue}{-6}, z, 6 \cdot \frac{2}{3}\right), y - x, x\right) \]
                        18. lift-/.f64N/A

                          \[\leadsto \mathsf{fma}\left(\mathsf{fma}\left(-6, z, 6 \cdot \color{blue}{\frac{2}{3}}\right), y - x, x\right) \]
                        19. metadata-evalN/A

                          \[\leadsto \mathsf{fma}\left(\mathsf{fma}\left(-6, z, 6 \cdot \color{blue}{\frac{2}{3}}\right), y - x, x\right) \]
                        20. metadata-eval99.6

                          \[\leadsto \mathsf{fma}\left(\mathsf{fma}\left(-6, z, \color{blue}{4}\right), y - x, x\right) \]
                      4. Applied rewrites99.6%

                        \[\leadsto \color{blue}{\mathsf{fma}\left(\mathsf{fma}\left(-6, z, 4\right), y - x, x\right)} \]
                      5. Taylor expanded in z around inf

                        \[\leadsto \color{blue}{-6 \cdot \left(z \cdot \left(y - x\right)\right)} \]
                      6. Step-by-step derivation
                        1. *-commutativeN/A

                          \[\leadsto -6 \cdot \color{blue}{\left(\left(y - x\right) \cdot z\right)} \]
                        2. associate-*r*N/A

                          \[\leadsto \color{blue}{\left(-6 \cdot \left(y - x\right)\right) \cdot z} \]
                        3. lower-*.f64N/A

                          \[\leadsto \color{blue}{\left(-6 \cdot \left(y - x\right)\right) \cdot z} \]
                        4. lower-*.f64N/A

                          \[\leadsto \color{blue}{\left(-6 \cdot \left(y - x\right)\right)} \cdot z \]
                        5. lower--.f6497.8

                          \[\leadsto \left(-6 \cdot \color{blue}{\left(y - x\right)}\right) \cdot z \]
                      7. Applied rewrites97.8%

                        \[\leadsto \color{blue}{\left(-6 \cdot \left(y - x\right)\right) \cdot z} \]
                      8. Taylor expanded in x around 0

                        \[\leadsto \left(-6 \cdot y\right) \cdot z \]
                      9. Step-by-step derivation
                        1. Applied rewrites58.8%

                          \[\leadsto \left(-6 \cdot y\right) \cdot z \]
                      10. Recombined 3 regimes into one program.
                      11. Final simplification83.2%

                        \[\leadsto \begin{array}{l} \mathbf{if}\;\frac{2}{3} - z \leq -10000000000:\\ \;\;\;\;\left(6 \cdot x\right) \cdot z\\ \mathbf{elif}\;\frac{2}{3} - z \leq 1:\\ \;\;\;\;\mathsf{fma}\left(4, y - x, x\right)\\ \mathbf{elif}\;\frac{2}{3} - z \leq 10^{+156}:\\ \;\;\;\;\left(y \cdot -6\right) \cdot z\\ \mathbf{else}:\\ \;\;\;\;\left(6 \cdot x\right) \cdot z\\ \end{array} \]
                      12. Add Preprocessing

                      Alternative 7: 74.4% accurate, 0.4× speedup?

                      \[\begin{array}{l} \\ \begin{array}{l} t_0 := \frac{2}{3} - z\\ t_1 := \left(6 \cdot z\right) \cdot x\\ \mathbf{if}\;t\_0 \leq -10000000000:\\ \;\;\;\;t\_1\\ \mathbf{elif}\;t\_0 \leq 1:\\ \;\;\;\;\mathsf{fma}\left(4, y - x, x\right)\\ \mathbf{elif}\;t\_0 \leq 10^{+156}:\\ \;\;\;\;\left(y \cdot -6\right) \cdot z\\ \mathbf{else}:\\ \;\;\;\;t\_1\\ \end{array} \end{array} \]
                      (FPCore (x y z)
                       :precision binary64
                       (let* ((t_0 (- (/ 2.0 3.0) z)) (t_1 (* (* 6.0 z) x)))
                         (if (<= t_0 -10000000000.0)
                           t_1
                           (if (<= t_0 1.0)
                             (fma 4.0 (- y x) x)
                             (if (<= t_0 1e+156) (* (* y -6.0) z) t_1)))))
                      double code(double x, double y, double z) {
                      	double t_0 = (2.0 / 3.0) - z;
                      	double t_1 = (6.0 * z) * x;
                      	double tmp;
                      	if (t_0 <= -10000000000.0) {
                      		tmp = t_1;
                      	} else if (t_0 <= 1.0) {
                      		tmp = fma(4.0, (y - x), x);
                      	} else if (t_0 <= 1e+156) {
                      		tmp = (y * -6.0) * z;
                      	} else {
                      		tmp = t_1;
                      	}
                      	return tmp;
                      }
                      
                      function code(x, y, z)
                      	t_0 = Float64(Float64(2.0 / 3.0) - z)
                      	t_1 = Float64(Float64(6.0 * z) * x)
                      	tmp = 0.0
                      	if (t_0 <= -10000000000.0)
                      		tmp = t_1;
                      	elseif (t_0 <= 1.0)
                      		tmp = fma(4.0, Float64(y - x), x);
                      	elseif (t_0 <= 1e+156)
                      		tmp = Float64(Float64(y * -6.0) * z);
                      	else
                      		tmp = t_1;
                      	end
                      	return tmp
                      end
                      
                      code[x_, y_, z_] := Block[{t$95$0 = N[(N[(2.0 / 3.0), $MachinePrecision] - z), $MachinePrecision]}, Block[{t$95$1 = N[(N[(6.0 * z), $MachinePrecision] * x), $MachinePrecision]}, If[LessEqual[t$95$0, -10000000000.0], t$95$1, If[LessEqual[t$95$0, 1.0], N[(4.0 * N[(y - x), $MachinePrecision] + x), $MachinePrecision], If[LessEqual[t$95$0, 1e+156], N[(N[(y * -6.0), $MachinePrecision] * z), $MachinePrecision], t$95$1]]]]]
                      
                      \begin{array}{l}
                      
                      \\
                      \begin{array}{l}
                      t_0 := \frac{2}{3} - z\\
                      t_1 := \left(6 \cdot z\right) \cdot x\\
                      \mathbf{if}\;t\_0 \leq -10000000000:\\
                      \;\;\;\;t\_1\\
                      
                      \mathbf{elif}\;t\_0 \leq 1:\\
                      \;\;\;\;\mathsf{fma}\left(4, y - x, x\right)\\
                      
                      \mathbf{elif}\;t\_0 \leq 10^{+156}:\\
                      \;\;\;\;\left(y \cdot -6\right) \cdot z\\
                      
                      \mathbf{else}:\\
                      \;\;\;\;t\_1\\
                      
                      
                      \end{array}
                      \end{array}
                      
                      Derivation
                      1. Split input into 3 regimes
                      2. if (-.f64 (/.f64 #s(literal 2 binary64) #s(literal 3 binary64)) z) < -1e10 or 9.9999999999999998e155 < (-.f64 (/.f64 #s(literal 2 binary64) #s(literal 3 binary64)) z)

                        1. Initial program 99.7%

                          \[x + \left(\left(y - x\right) \cdot 6\right) \cdot \left(\frac{2}{3} - z\right) \]
                        2. Add Preprocessing
                        3. Step-by-step derivation
                          1. lift-+.f64N/A

                            \[\leadsto \color{blue}{x + \left(\left(y - x\right) \cdot 6\right) \cdot \left(\frac{2}{3} - z\right)} \]
                          2. +-commutativeN/A

                            \[\leadsto \color{blue}{\left(\left(y - x\right) \cdot 6\right) \cdot \left(\frac{2}{3} - z\right) + x} \]
                          3. lift-*.f64N/A

                            \[\leadsto \color{blue}{\left(\left(y - x\right) \cdot 6\right) \cdot \left(\frac{2}{3} - z\right)} + x \]
                          4. lift-*.f64N/A

                            \[\leadsto \color{blue}{\left(\left(y - x\right) \cdot 6\right)} \cdot \left(\frac{2}{3} - z\right) + x \]
                          5. associate-*l*N/A

                            \[\leadsto \color{blue}{\left(y - x\right) \cdot \left(6 \cdot \left(\frac{2}{3} - z\right)\right)} + x \]
                          6. *-commutativeN/A

                            \[\leadsto \color{blue}{\left(6 \cdot \left(\frac{2}{3} - z\right)\right) \cdot \left(y - x\right)} + x \]
                          7. lower-fma.f64N/A

                            \[\leadsto \color{blue}{\mathsf{fma}\left(6 \cdot \left(\frac{2}{3} - z\right), y - x, x\right)} \]
                          8. lift--.f64N/A

                            \[\leadsto \mathsf{fma}\left(6 \cdot \color{blue}{\left(\frac{2}{3} - z\right)}, y - x, x\right) \]
                          9. sub-negN/A

                            \[\leadsto \mathsf{fma}\left(6 \cdot \color{blue}{\left(\frac{2}{3} + \left(\mathsf{neg}\left(z\right)\right)\right)}, y - x, x\right) \]
                          10. +-commutativeN/A

                            \[\leadsto \mathsf{fma}\left(6 \cdot \color{blue}{\left(\left(\mathsf{neg}\left(z\right)\right) + \frac{2}{3}\right)}, y - x, x\right) \]
                          11. distribute-lft-inN/A

                            \[\leadsto \mathsf{fma}\left(\color{blue}{6 \cdot \left(\mathsf{neg}\left(z\right)\right) + 6 \cdot \frac{2}{3}}, y - x, x\right) \]
                          12. neg-mul-1N/A

                            \[\leadsto \mathsf{fma}\left(6 \cdot \color{blue}{\left(-1 \cdot z\right)} + 6 \cdot \frac{2}{3}, y - x, x\right) \]
                          13. associate-*r*N/A

                            \[\leadsto \mathsf{fma}\left(\color{blue}{\left(6 \cdot -1\right) \cdot z} + 6 \cdot \frac{2}{3}, y - x, x\right) \]
                          14. metadata-evalN/A

                            \[\leadsto \mathsf{fma}\left(\color{blue}{-6} \cdot z + 6 \cdot \frac{2}{3}, y - x, x\right) \]
                          15. metadata-evalN/A

                            \[\leadsto \mathsf{fma}\left(\color{blue}{\left(\mathsf{neg}\left(6\right)\right)} \cdot z + 6 \cdot \frac{2}{3}, y - x, x\right) \]
                          16. lower-fma.f64N/A

                            \[\leadsto \mathsf{fma}\left(\color{blue}{\mathsf{fma}\left(\mathsf{neg}\left(6\right), z, 6 \cdot \frac{2}{3}\right)}, y - x, x\right) \]
                          17. metadata-evalN/A

                            \[\leadsto \mathsf{fma}\left(\mathsf{fma}\left(\color{blue}{-6}, z, 6 \cdot \frac{2}{3}\right), y - x, x\right) \]
                          18. lift-/.f64N/A

                            \[\leadsto \mathsf{fma}\left(\mathsf{fma}\left(-6, z, 6 \cdot \color{blue}{\frac{2}{3}}\right), y - x, x\right) \]
                          19. metadata-evalN/A

                            \[\leadsto \mathsf{fma}\left(\mathsf{fma}\left(-6, z, 6 \cdot \color{blue}{\frac{2}{3}}\right), y - x, x\right) \]
                          20. metadata-eval99.8

                            \[\leadsto \mathsf{fma}\left(\mathsf{fma}\left(-6, z, \color{blue}{4}\right), y - x, x\right) \]
                        4. Applied rewrites99.8%

                          \[\leadsto \color{blue}{\mathsf{fma}\left(\mathsf{fma}\left(-6, z, 4\right), y - x, x\right)} \]
                        5. Taylor expanded in z around inf

                          \[\leadsto \color{blue}{-6 \cdot \left(z \cdot \left(y - x\right)\right)} \]
                        6. Step-by-step derivation
                          1. *-commutativeN/A

                            \[\leadsto -6 \cdot \color{blue}{\left(\left(y - x\right) \cdot z\right)} \]
                          2. associate-*r*N/A

                            \[\leadsto \color{blue}{\left(-6 \cdot \left(y - x\right)\right) \cdot z} \]
                          3. lower-*.f64N/A

                            \[\leadsto \color{blue}{\left(-6 \cdot \left(y - x\right)\right) \cdot z} \]
                          4. lower-*.f64N/A

                            \[\leadsto \color{blue}{\left(-6 \cdot \left(y - x\right)\right)} \cdot z \]
                          5. lower--.f6499.1

                            \[\leadsto \left(-6 \cdot \color{blue}{\left(y - x\right)}\right) \cdot z \]
                        7. Applied rewrites99.1%

                          \[\leadsto \color{blue}{\left(-6 \cdot \left(y - x\right)\right) \cdot z} \]
                        8. Taylor expanded in x around inf

                          \[\leadsto 6 \cdot \color{blue}{\left(x \cdot z\right)} \]
                        9. Step-by-step derivation
                          1. Applied rewrites63.5%

                            \[\leadsto \left(6 \cdot z\right) \cdot \color{blue}{x} \]

                          if -1e10 < (-.f64 (/.f64 #s(literal 2 binary64) #s(literal 3 binary64)) z) < 1

                          1. Initial program 98.8%

                            \[x + \left(\left(y - x\right) \cdot 6\right) \cdot \left(\frac{2}{3} - z\right) \]
                          2. Add Preprocessing
                          3. Taylor expanded in z around 0

                            \[\leadsto \color{blue}{x + 4 \cdot \left(y - x\right)} \]
                          4. Step-by-step derivation
                            1. +-commutativeN/A

                              \[\leadsto \color{blue}{4 \cdot \left(y - x\right) + x} \]
                            2. lower-fma.f64N/A

                              \[\leadsto \color{blue}{\mathsf{fma}\left(4, y - x, x\right)} \]
                            3. lower--.f6499.6

                              \[\leadsto \mathsf{fma}\left(4, \color{blue}{y - x}, x\right) \]
                          5. Applied rewrites99.6%

                            \[\leadsto \color{blue}{\mathsf{fma}\left(4, y - x, x\right)} \]

                          if 1 < (-.f64 (/.f64 #s(literal 2 binary64) #s(literal 3 binary64)) z) < 9.9999999999999998e155

                          1. Initial program 99.3%

                            \[x + \left(\left(y - x\right) \cdot 6\right) \cdot \left(\frac{2}{3} - z\right) \]
                          2. Add Preprocessing
                          3. Step-by-step derivation
                            1. lift-+.f64N/A

                              \[\leadsto \color{blue}{x + \left(\left(y - x\right) \cdot 6\right) \cdot \left(\frac{2}{3} - z\right)} \]
                            2. +-commutativeN/A

                              \[\leadsto \color{blue}{\left(\left(y - x\right) \cdot 6\right) \cdot \left(\frac{2}{3} - z\right) + x} \]
                            3. lift-*.f64N/A

                              \[\leadsto \color{blue}{\left(\left(y - x\right) \cdot 6\right) \cdot \left(\frac{2}{3} - z\right)} + x \]
                            4. lift-*.f64N/A

                              \[\leadsto \color{blue}{\left(\left(y - x\right) \cdot 6\right)} \cdot \left(\frac{2}{3} - z\right) + x \]
                            5. associate-*l*N/A

                              \[\leadsto \color{blue}{\left(y - x\right) \cdot \left(6 \cdot \left(\frac{2}{3} - z\right)\right)} + x \]
                            6. *-commutativeN/A

                              \[\leadsto \color{blue}{\left(6 \cdot \left(\frac{2}{3} - z\right)\right) \cdot \left(y - x\right)} + x \]
                            7. lower-fma.f64N/A

                              \[\leadsto \color{blue}{\mathsf{fma}\left(6 \cdot \left(\frac{2}{3} - z\right), y - x, x\right)} \]
                            8. lift--.f64N/A

                              \[\leadsto \mathsf{fma}\left(6 \cdot \color{blue}{\left(\frac{2}{3} - z\right)}, y - x, x\right) \]
                            9. sub-negN/A

                              \[\leadsto \mathsf{fma}\left(6 \cdot \color{blue}{\left(\frac{2}{3} + \left(\mathsf{neg}\left(z\right)\right)\right)}, y - x, x\right) \]
                            10. +-commutativeN/A

                              \[\leadsto \mathsf{fma}\left(6 \cdot \color{blue}{\left(\left(\mathsf{neg}\left(z\right)\right) + \frac{2}{3}\right)}, y - x, x\right) \]
                            11. distribute-lft-inN/A

                              \[\leadsto \mathsf{fma}\left(\color{blue}{6 \cdot \left(\mathsf{neg}\left(z\right)\right) + 6 \cdot \frac{2}{3}}, y - x, x\right) \]
                            12. neg-mul-1N/A

                              \[\leadsto \mathsf{fma}\left(6 \cdot \color{blue}{\left(-1 \cdot z\right)} + 6 \cdot \frac{2}{3}, y - x, x\right) \]
                            13. associate-*r*N/A

                              \[\leadsto \mathsf{fma}\left(\color{blue}{\left(6 \cdot -1\right) \cdot z} + 6 \cdot \frac{2}{3}, y - x, x\right) \]
                            14. metadata-evalN/A

                              \[\leadsto \mathsf{fma}\left(\color{blue}{-6} \cdot z + 6 \cdot \frac{2}{3}, y - x, x\right) \]
                            15. metadata-evalN/A

                              \[\leadsto \mathsf{fma}\left(\color{blue}{\left(\mathsf{neg}\left(6\right)\right)} \cdot z + 6 \cdot \frac{2}{3}, y - x, x\right) \]
                            16. lower-fma.f64N/A

                              \[\leadsto \mathsf{fma}\left(\color{blue}{\mathsf{fma}\left(\mathsf{neg}\left(6\right), z, 6 \cdot \frac{2}{3}\right)}, y - x, x\right) \]
                            17. metadata-evalN/A

                              \[\leadsto \mathsf{fma}\left(\mathsf{fma}\left(\color{blue}{-6}, z, 6 \cdot \frac{2}{3}\right), y - x, x\right) \]
                            18. lift-/.f64N/A

                              \[\leadsto \mathsf{fma}\left(\mathsf{fma}\left(-6, z, 6 \cdot \color{blue}{\frac{2}{3}}\right), y - x, x\right) \]
                            19. metadata-evalN/A

                              \[\leadsto \mathsf{fma}\left(\mathsf{fma}\left(-6, z, 6 \cdot \color{blue}{\frac{2}{3}}\right), y - x, x\right) \]
                            20. metadata-eval99.6

                              \[\leadsto \mathsf{fma}\left(\mathsf{fma}\left(-6, z, \color{blue}{4}\right), y - x, x\right) \]
                          4. Applied rewrites99.6%

                            \[\leadsto \color{blue}{\mathsf{fma}\left(\mathsf{fma}\left(-6, z, 4\right), y - x, x\right)} \]
                          5. Taylor expanded in z around inf

                            \[\leadsto \color{blue}{-6 \cdot \left(z \cdot \left(y - x\right)\right)} \]
                          6. Step-by-step derivation
                            1. *-commutativeN/A

                              \[\leadsto -6 \cdot \color{blue}{\left(\left(y - x\right) \cdot z\right)} \]
                            2. associate-*r*N/A

                              \[\leadsto \color{blue}{\left(-6 \cdot \left(y - x\right)\right) \cdot z} \]
                            3. lower-*.f64N/A

                              \[\leadsto \color{blue}{\left(-6 \cdot \left(y - x\right)\right) \cdot z} \]
                            4. lower-*.f64N/A

                              \[\leadsto \color{blue}{\left(-6 \cdot \left(y - x\right)\right)} \cdot z \]
                            5. lower--.f6497.8

                              \[\leadsto \left(-6 \cdot \color{blue}{\left(y - x\right)}\right) \cdot z \]
                          7. Applied rewrites97.8%

                            \[\leadsto \color{blue}{\left(-6 \cdot \left(y - x\right)\right) \cdot z} \]
                          8. Taylor expanded in x around 0

                            \[\leadsto \left(-6 \cdot y\right) \cdot z \]
                          9. Step-by-step derivation
                            1. Applied rewrites58.8%

                              \[\leadsto \left(-6 \cdot y\right) \cdot z \]
                          10. Recombined 3 regimes into one program.
                          11. Final simplification83.2%

                            \[\leadsto \begin{array}{l} \mathbf{if}\;\frac{2}{3} - z \leq -10000000000:\\ \;\;\;\;\left(6 \cdot z\right) \cdot x\\ \mathbf{elif}\;\frac{2}{3} - z \leq 1:\\ \;\;\;\;\mathsf{fma}\left(4, y - x, x\right)\\ \mathbf{elif}\;\frac{2}{3} - z \leq 10^{+156}:\\ \;\;\;\;\left(y \cdot -6\right) \cdot z\\ \mathbf{else}:\\ \;\;\;\;\left(6 \cdot z\right) \cdot x\\ \end{array} \]
                          12. Add Preprocessing

                          Alternative 8: 97.4% accurate, 0.6× speedup?

                          \[\begin{array}{l} \\ \begin{array}{l} t_0 := \frac{2}{3} - z\\ t_1 := \left(z \cdot -6\right) \cdot \left(y - x\right)\\ \mathbf{if}\;t\_0 \leq -10000000000:\\ \;\;\;\;t\_1\\ \mathbf{elif}\;t\_0 \leq 1:\\ \;\;\;\;\mathsf{fma}\left(-3, x, y \cdot 4\right)\\ \mathbf{else}:\\ \;\;\;\;t\_1\\ \end{array} \end{array} \]
                          (FPCore (x y z)
                           :precision binary64
                           (let* ((t_0 (- (/ 2.0 3.0) z)) (t_1 (* (* z -6.0) (- y x))))
                             (if (<= t_0 -10000000000.0)
                               t_1
                               (if (<= t_0 1.0) (fma -3.0 x (* y 4.0)) t_1))))
                          double code(double x, double y, double z) {
                          	double t_0 = (2.0 / 3.0) - z;
                          	double t_1 = (z * -6.0) * (y - x);
                          	double tmp;
                          	if (t_0 <= -10000000000.0) {
                          		tmp = t_1;
                          	} else if (t_0 <= 1.0) {
                          		tmp = fma(-3.0, x, (y * 4.0));
                          	} else {
                          		tmp = t_1;
                          	}
                          	return tmp;
                          }
                          
                          function code(x, y, z)
                          	t_0 = Float64(Float64(2.0 / 3.0) - z)
                          	t_1 = Float64(Float64(z * -6.0) * Float64(y - x))
                          	tmp = 0.0
                          	if (t_0 <= -10000000000.0)
                          		tmp = t_1;
                          	elseif (t_0 <= 1.0)
                          		tmp = fma(-3.0, x, Float64(y * 4.0));
                          	else
                          		tmp = t_1;
                          	end
                          	return tmp
                          end
                          
                          code[x_, y_, z_] := Block[{t$95$0 = N[(N[(2.0 / 3.0), $MachinePrecision] - z), $MachinePrecision]}, Block[{t$95$1 = N[(N[(z * -6.0), $MachinePrecision] * N[(y - x), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$0, -10000000000.0], t$95$1, If[LessEqual[t$95$0, 1.0], N[(-3.0 * x + N[(y * 4.0), $MachinePrecision]), $MachinePrecision], t$95$1]]]]
                          
                          \begin{array}{l}
                          
                          \\
                          \begin{array}{l}
                          t_0 := \frac{2}{3} - z\\
                          t_1 := \left(z \cdot -6\right) \cdot \left(y - x\right)\\
                          \mathbf{if}\;t\_0 \leq -10000000000:\\
                          \;\;\;\;t\_1\\
                          
                          \mathbf{elif}\;t\_0 \leq 1:\\
                          \;\;\;\;\mathsf{fma}\left(-3, x, y \cdot 4\right)\\
                          
                          \mathbf{else}:\\
                          \;\;\;\;t\_1\\
                          
                          
                          \end{array}
                          \end{array}
                          
                          Derivation
                          1. Split input into 2 regimes
                          2. if (-.f64 (/.f64 #s(literal 2 binary64) #s(literal 3 binary64)) z) < -1e10 or 1 < (-.f64 (/.f64 #s(literal 2 binary64) #s(literal 3 binary64)) z)

                            1. Initial program 99.7%

                              \[x + \left(\left(y - x\right) \cdot 6\right) \cdot \left(\frac{2}{3} - z\right) \]
                            2. Add Preprocessing
                            3. Step-by-step derivation
                              1. lift-+.f64N/A

                                \[\leadsto \color{blue}{x + \left(\left(y - x\right) \cdot 6\right) \cdot \left(\frac{2}{3} - z\right)} \]
                              2. +-commutativeN/A

                                \[\leadsto \color{blue}{\left(\left(y - x\right) \cdot 6\right) \cdot \left(\frac{2}{3} - z\right) + x} \]
                              3. lift-*.f64N/A

                                \[\leadsto \color{blue}{\left(\left(y - x\right) \cdot 6\right) \cdot \left(\frac{2}{3} - z\right)} + x \]
                              4. lift-*.f64N/A

                                \[\leadsto \color{blue}{\left(\left(y - x\right) \cdot 6\right)} \cdot \left(\frac{2}{3} - z\right) + x \]
                              5. associate-*l*N/A

                                \[\leadsto \color{blue}{\left(y - x\right) \cdot \left(6 \cdot \left(\frac{2}{3} - z\right)\right)} + x \]
                              6. *-commutativeN/A

                                \[\leadsto \color{blue}{\left(6 \cdot \left(\frac{2}{3} - z\right)\right) \cdot \left(y - x\right)} + x \]
                              7. lower-fma.f64N/A

                                \[\leadsto \color{blue}{\mathsf{fma}\left(6 \cdot \left(\frac{2}{3} - z\right), y - x, x\right)} \]
                              8. lift--.f64N/A

                                \[\leadsto \mathsf{fma}\left(6 \cdot \color{blue}{\left(\frac{2}{3} - z\right)}, y - x, x\right) \]
                              9. sub-negN/A

                                \[\leadsto \mathsf{fma}\left(6 \cdot \color{blue}{\left(\frac{2}{3} + \left(\mathsf{neg}\left(z\right)\right)\right)}, y - x, x\right) \]
                              10. +-commutativeN/A

                                \[\leadsto \mathsf{fma}\left(6 \cdot \color{blue}{\left(\left(\mathsf{neg}\left(z\right)\right) + \frac{2}{3}\right)}, y - x, x\right) \]
                              11. distribute-lft-inN/A

                                \[\leadsto \mathsf{fma}\left(\color{blue}{6 \cdot \left(\mathsf{neg}\left(z\right)\right) + 6 \cdot \frac{2}{3}}, y - x, x\right) \]
                              12. neg-mul-1N/A

                                \[\leadsto \mathsf{fma}\left(6 \cdot \color{blue}{\left(-1 \cdot z\right)} + 6 \cdot \frac{2}{3}, y - x, x\right) \]
                              13. associate-*r*N/A

                                \[\leadsto \mathsf{fma}\left(\color{blue}{\left(6 \cdot -1\right) \cdot z} + 6 \cdot \frac{2}{3}, y - x, x\right) \]
                              14. metadata-evalN/A

                                \[\leadsto \mathsf{fma}\left(\color{blue}{-6} \cdot z + 6 \cdot \frac{2}{3}, y - x, x\right) \]
                              15. metadata-evalN/A

                                \[\leadsto \mathsf{fma}\left(\color{blue}{\left(\mathsf{neg}\left(6\right)\right)} \cdot z + 6 \cdot \frac{2}{3}, y - x, x\right) \]
                              16. lower-fma.f64N/A

                                \[\leadsto \mathsf{fma}\left(\color{blue}{\mathsf{fma}\left(\mathsf{neg}\left(6\right), z, 6 \cdot \frac{2}{3}\right)}, y - x, x\right) \]
                              17. metadata-evalN/A

                                \[\leadsto \mathsf{fma}\left(\mathsf{fma}\left(\color{blue}{-6}, z, 6 \cdot \frac{2}{3}\right), y - x, x\right) \]
                              18. lift-/.f64N/A

                                \[\leadsto \mathsf{fma}\left(\mathsf{fma}\left(-6, z, 6 \cdot \color{blue}{\frac{2}{3}}\right), y - x, x\right) \]
                              19. metadata-evalN/A

                                \[\leadsto \mathsf{fma}\left(\mathsf{fma}\left(-6, z, 6 \cdot \color{blue}{\frac{2}{3}}\right), y - x, x\right) \]
                              20. metadata-eval99.8

                                \[\leadsto \mathsf{fma}\left(\mathsf{fma}\left(-6, z, \color{blue}{4}\right), y - x, x\right) \]
                            4. Applied rewrites99.8%

                              \[\leadsto \color{blue}{\mathsf{fma}\left(\mathsf{fma}\left(-6, z, 4\right), y - x, x\right)} \]
                            5. Taylor expanded in z around inf

                              \[\leadsto \color{blue}{-6 \cdot \left(z \cdot \left(y - x\right)\right)} \]
                            6. Step-by-step derivation
                              1. *-commutativeN/A

                                \[\leadsto -6 \cdot \color{blue}{\left(\left(y - x\right) \cdot z\right)} \]
                              2. associate-*r*N/A

                                \[\leadsto \color{blue}{\left(-6 \cdot \left(y - x\right)\right) \cdot z} \]
                              3. lower-*.f64N/A

                                \[\leadsto \color{blue}{\left(-6 \cdot \left(y - x\right)\right) \cdot z} \]
                              4. lower-*.f64N/A

                                \[\leadsto \color{blue}{\left(-6 \cdot \left(y - x\right)\right)} \cdot z \]
                              5. lower--.f6498.9

                                \[\leadsto \left(-6 \cdot \color{blue}{\left(y - x\right)}\right) \cdot z \]
                            7. Applied rewrites98.9%

                              \[\leadsto \color{blue}{\left(-6 \cdot \left(y - x\right)\right) \cdot z} \]
                            8. Step-by-step derivation
                              1. Applied rewrites99.0%

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

                              if -1e10 < (-.f64 (/.f64 #s(literal 2 binary64) #s(literal 3 binary64)) z) < 1

                              1. Initial program 98.8%

                                \[x + \left(\left(y - x\right) \cdot 6\right) \cdot \left(\frac{2}{3} - z\right) \]
                              2. Add Preprocessing
                              3. Taylor expanded in z around 0

                                \[\leadsto \color{blue}{x + 4 \cdot \left(y - x\right)} \]
                              4. Step-by-step derivation
                                1. +-commutativeN/A

                                  \[\leadsto \color{blue}{4 \cdot \left(y - x\right) + x} \]
                                2. lower-fma.f64N/A

                                  \[\leadsto \color{blue}{\mathsf{fma}\left(4, y - x, x\right)} \]
                                3. lower--.f6499.6

                                  \[\leadsto \mathsf{fma}\left(4, \color{blue}{y - x}, x\right) \]
                              5. Applied rewrites99.6%

                                \[\leadsto \color{blue}{\mathsf{fma}\left(4, y - x, x\right)} \]
                              6. Taylor expanded in x around 0

                                \[\leadsto -3 \cdot x + \color{blue}{4 \cdot y} \]
                              7. Step-by-step derivation
                                1. Applied rewrites99.7%

                                  \[\leadsto \mathsf{fma}\left(-3, \color{blue}{x}, 4 \cdot y\right) \]
                              8. Recombined 2 regimes into one program.
                              9. Final simplification99.4%

                                \[\leadsto \begin{array}{l} \mathbf{if}\;\frac{2}{3} - z \leq -10000000000:\\ \;\;\;\;\left(z \cdot -6\right) \cdot \left(y - x\right)\\ \mathbf{elif}\;\frac{2}{3} - z \leq 1:\\ \;\;\;\;\mathsf{fma}\left(-3, x, y \cdot 4\right)\\ \mathbf{else}:\\ \;\;\;\;\left(z \cdot -6\right) \cdot \left(y - x\right)\\ \end{array} \]
                              10. Add Preprocessing

                              Alternative 9: 97.4% accurate, 0.6× speedup?

                              \[\begin{array}{l} \\ \begin{array}{l} t_0 := \frac{2}{3} - z\\ \mathbf{if}\;t\_0 \leq -10000000000:\\ \;\;\;\;\left(\left(y - x\right) \cdot -6\right) \cdot z\\ \mathbf{elif}\;t\_0 \leq 1:\\ \;\;\;\;\mathsf{fma}\left(-3, x, y \cdot 4\right)\\ \mathbf{else}:\\ \;\;\;\;\left(\left(y - x\right) \cdot z\right) \cdot -6\\ \end{array} \end{array} \]
                              (FPCore (x y z)
                               :precision binary64
                               (let* ((t_0 (- (/ 2.0 3.0) z)))
                                 (if (<= t_0 -10000000000.0)
                                   (* (* (- y x) -6.0) z)
                                   (if (<= t_0 1.0) (fma -3.0 x (* y 4.0)) (* (* (- y x) z) -6.0)))))
                              double code(double x, double y, double z) {
                              	double t_0 = (2.0 / 3.0) - z;
                              	double tmp;
                              	if (t_0 <= -10000000000.0) {
                              		tmp = ((y - x) * -6.0) * z;
                              	} else if (t_0 <= 1.0) {
                              		tmp = fma(-3.0, x, (y * 4.0));
                              	} else {
                              		tmp = ((y - x) * z) * -6.0;
                              	}
                              	return tmp;
                              }
                              
                              function code(x, y, z)
                              	t_0 = Float64(Float64(2.0 / 3.0) - z)
                              	tmp = 0.0
                              	if (t_0 <= -10000000000.0)
                              		tmp = Float64(Float64(Float64(y - x) * -6.0) * z);
                              	elseif (t_0 <= 1.0)
                              		tmp = fma(-3.0, x, Float64(y * 4.0));
                              	else
                              		tmp = Float64(Float64(Float64(y - x) * z) * -6.0);
                              	end
                              	return tmp
                              end
                              
                              code[x_, y_, z_] := Block[{t$95$0 = N[(N[(2.0 / 3.0), $MachinePrecision] - z), $MachinePrecision]}, If[LessEqual[t$95$0, -10000000000.0], N[(N[(N[(y - x), $MachinePrecision] * -6.0), $MachinePrecision] * z), $MachinePrecision], If[LessEqual[t$95$0, 1.0], N[(-3.0 * x + N[(y * 4.0), $MachinePrecision]), $MachinePrecision], N[(N[(N[(y - x), $MachinePrecision] * z), $MachinePrecision] * -6.0), $MachinePrecision]]]]
                              
                              \begin{array}{l}
                              
                              \\
                              \begin{array}{l}
                              t_0 := \frac{2}{3} - z\\
                              \mathbf{if}\;t\_0 \leq -10000000000:\\
                              \;\;\;\;\left(\left(y - x\right) \cdot -6\right) \cdot z\\
                              
                              \mathbf{elif}\;t\_0 \leq 1:\\
                              \;\;\;\;\mathsf{fma}\left(-3, x, y \cdot 4\right)\\
                              
                              \mathbf{else}:\\
                              \;\;\;\;\left(\left(y - x\right) \cdot z\right) \cdot -6\\
                              
                              
                              \end{array}
                              \end{array}
                              
                              Derivation
                              1. Split input into 3 regimes
                              2. if (-.f64 (/.f64 #s(literal 2 binary64) #s(literal 3 binary64)) z) < -1e10

                                1. Initial program 99.7%

                                  \[x + \left(\left(y - x\right) \cdot 6\right) \cdot \left(\frac{2}{3} - z\right) \]
                                2. Add Preprocessing
                                3. Step-by-step derivation
                                  1. lift-+.f64N/A

                                    \[\leadsto \color{blue}{x + \left(\left(y - x\right) \cdot 6\right) \cdot \left(\frac{2}{3} - z\right)} \]
                                  2. +-commutativeN/A

                                    \[\leadsto \color{blue}{\left(\left(y - x\right) \cdot 6\right) \cdot \left(\frac{2}{3} - z\right) + x} \]
                                  3. lift-*.f64N/A

                                    \[\leadsto \color{blue}{\left(\left(y - x\right) \cdot 6\right) \cdot \left(\frac{2}{3} - z\right)} + x \]
                                  4. lift-*.f64N/A

                                    \[\leadsto \color{blue}{\left(\left(y - x\right) \cdot 6\right)} \cdot \left(\frac{2}{3} - z\right) + x \]
                                  5. associate-*l*N/A

                                    \[\leadsto \color{blue}{\left(y - x\right) \cdot \left(6 \cdot \left(\frac{2}{3} - z\right)\right)} + x \]
                                  6. *-commutativeN/A

                                    \[\leadsto \color{blue}{\left(6 \cdot \left(\frac{2}{3} - z\right)\right) \cdot \left(y - x\right)} + x \]
                                  7. lower-fma.f64N/A

                                    \[\leadsto \color{blue}{\mathsf{fma}\left(6 \cdot \left(\frac{2}{3} - z\right), y - x, x\right)} \]
                                  8. lift--.f64N/A

                                    \[\leadsto \mathsf{fma}\left(6 \cdot \color{blue}{\left(\frac{2}{3} - z\right)}, y - x, x\right) \]
                                  9. sub-negN/A

                                    \[\leadsto \mathsf{fma}\left(6 \cdot \color{blue}{\left(\frac{2}{3} + \left(\mathsf{neg}\left(z\right)\right)\right)}, y - x, x\right) \]
                                  10. +-commutativeN/A

                                    \[\leadsto \mathsf{fma}\left(6 \cdot \color{blue}{\left(\left(\mathsf{neg}\left(z\right)\right) + \frac{2}{3}\right)}, y - x, x\right) \]
                                  11. distribute-lft-inN/A

                                    \[\leadsto \mathsf{fma}\left(\color{blue}{6 \cdot \left(\mathsf{neg}\left(z\right)\right) + 6 \cdot \frac{2}{3}}, y - x, x\right) \]
                                  12. neg-mul-1N/A

                                    \[\leadsto \mathsf{fma}\left(6 \cdot \color{blue}{\left(-1 \cdot z\right)} + 6 \cdot \frac{2}{3}, y - x, x\right) \]
                                  13. associate-*r*N/A

                                    \[\leadsto \mathsf{fma}\left(\color{blue}{\left(6 \cdot -1\right) \cdot z} + 6 \cdot \frac{2}{3}, y - x, x\right) \]
                                  14. metadata-evalN/A

                                    \[\leadsto \mathsf{fma}\left(\color{blue}{-6} \cdot z + 6 \cdot \frac{2}{3}, y - x, x\right) \]
                                  15. metadata-evalN/A

                                    \[\leadsto \mathsf{fma}\left(\color{blue}{\left(\mathsf{neg}\left(6\right)\right)} \cdot z + 6 \cdot \frac{2}{3}, y - x, x\right) \]
                                  16. lower-fma.f64N/A

                                    \[\leadsto \mathsf{fma}\left(\color{blue}{\mathsf{fma}\left(\mathsf{neg}\left(6\right), z, 6 \cdot \frac{2}{3}\right)}, y - x, x\right) \]
                                  17. metadata-evalN/A

                                    \[\leadsto \mathsf{fma}\left(\mathsf{fma}\left(\color{blue}{-6}, z, 6 \cdot \frac{2}{3}\right), y - x, x\right) \]
                                  18. lift-/.f64N/A

                                    \[\leadsto \mathsf{fma}\left(\mathsf{fma}\left(-6, z, 6 \cdot \color{blue}{\frac{2}{3}}\right), y - x, x\right) \]
                                  19. metadata-evalN/A

                                    \[\leadsto \mathsf{fma}\left(\mathsf{fma}\left(-6, z, 6 \cdot \color{blue}{\frac{2}{3}}\right), y - x, x\right) \]
                                  20. metadata-eval99.7

                                    \[\leadsto \mathsf{fma}\left(\mathsf{fma}\left(-6, z, \color{blue}{4}\right), y - x, x\right) \]
                                4. Applied rewrites99.7%

                                  \[\leadsto \color{blue}{\mathsf{fma}\left(\mathsf{fma}\left(-6, z, 4\right), y - x, x\right)} \]
                                5. Taylor expanded in z around inf

                                  \[\leadsto \color{blue}{-6 \cdot \left(z \cdot \left(y - x\right)\right)} \]
                                6. Step-by-step derivation
                                  1. *-commutativeN/A

                                    \[\leadsto -6 \cdot \color{blue}{\left(\left(y - x\right) \cdot z\right)} \]
                                  2. associate-*r*N/A

                                    \[\leadsto \color{blue}{\left(-6 \cdot \left(y - x\right)\right) \cdot z} \]
                                  3. lower-*.f64N/A

                                    \[\leadsto \color{blue}{\left(-6 \cdot \left(y - x\right)\right) \cdot z} \]
                                  4. lower-*.f64N/A

                                    \[\leadsto \color{blue}{\left(-6 \cdot \left(y - x\right)\right)} \cdot z \]
                                  5. lower--.f6498.8

                                    \[\leadsto \left(-6 \cdot \color{blue}{\left(y - x\right)}\right) \cdot z \]
                                7. Applied rewrites98.8%

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

                                if -1e10 < (-.f64 (/.f64 #s(literal 2 binary64) #s(literal 3 binary64)) z) < 1

                                1. Initial program 98.8%

                                  \[x + \left(\left(y - x\right) \cdot 6\right) \cdot \left(\frac{2}{3} - z\right) \]
                                2. Add Preprocessing
                                3. Taylor expanded in z around 0

                                  \[\leadsto \color{blue}{x + 4 \cdot \left(y - x\right)} \]
                                4. Step-by-step derivation
                                  1. +-commutativeN/A

                                    \[\leadsto \color{blue}{4 \cdot \left(y - x\right) + x} \]
                                  2. lower-fma.f64N/A

                                    \[\leadsto \color{blue}{\mathsf{fma}\left(4, y - x, x\right)} \]
                                  3. lower--.f6499.6

                                    \[\leadsto \mathsf{fma}\left(4, \color{blue}{y - x}, x\right) \]
                                5. Applied rewrites99.6%

                                  \[\leadsto \color{blue}{\mathsf{fma}\left(4, y - x, x\right)} \]
                                6. Taylor expanded in x around 0

                                  \[\leadsto -3 \cdot x + \color{blue}{4 \cdot y} \]
                                7. Step-by-step derivation
                                  1. Applied rewrites99.7%

                                    \[\leadsto \mathsf{fma}\left(-3, \color{blue}{x}, 4 \cdot y\right) \]

                                  if 1 < (-.f64 (/.f64 #s(literal 2 binary64) #s(literal 3 binary64)) z)

                                  1. Initial program 99.6%

                                    \[x + \left(\left(y - x\right) \cdot 6\right) \cdot \left(\frac{2}{3} - z\right) \]
                                  2. Add Preprocessing
                                  3. Taylor expanded in z around inf

                                    \[\leadsto \color{blue}{-6 \cdot \left(z \cdot \left(y - x\right)\right)} \]
                                  4. Step-by-step derivation
                                    1. *-commutativeN/A

                                      \[\leadsto \color{blue}{\left(z \cdot \left(y - x\right)\right) \cdot -6} \]
                                    2. lower-*.f64N/A

                                      \[\leadsto \color{blue}{\left(z \cdot \left(y - x\right)\right) \cdot -6} \]
                                    3. *-commutativeN/A

                                      \[\leadsto \color{blue}{\left(\left(y - x\right) \cdot z\right)} \cdot -6 \]
                                    4. lower-*.f64N/A

                                      \[\leadsto \color{blue}{\left(\left(y - x\right) \cdot z\right)} \cdot -6 \]
                                    5. lower--.f6499.0

                                      \[\leadsto \left(\color{blue}{\left(y - x\right)} \cdot z\right) \cdot -6 \]
                                  5. Applied rewrites99.0%

                                    \[\leadsto \color{blue}{\left(\left(y - x\right) \cdot z\right) \cdot -6} \]
                                8. Recombined 3 regimes into one program.
                                9. Final simplification99.4%

                                  \[\leadsto \begin{array}{l} \mathbf{if}\;\frac{2}{3} - z \leq -10000000000:\\ \;\;\;\;\left(\left(y - x\right) \cdot -6\right) \cdot z\\ \mathbf{elif}\;\frac{2}{3} - z \leq 1:\\ \;\;\;\;\mathsf{fma}\left(-3, x, y \cdot 4\right)\\ \mathbf{else}:\\ \;\;\;\;\left(\left(y - x\right) \cdot z\right) \cdot -6\\ \end{array} \]
                                10. Add Preprocessing

                                Alternative 10: 97.4% accurate, 0.6× speedup?

                                \[\begin{array}{l} \\ \begin{array}{l} t_0 := \frac{2}{3} - z\\ t_1 := \left(\left(y - x\right) \cdot z\right) \cdot -6\\ \mathbf{if}\;t\_0 \leq -10000000000:\\ \;\;\;\;t\_1\\ \mathbf{elif}\;t\_0 \leq 1:\\ \;\;\;\;\mathsf{fma}\left(-3, x, y \cdot 4\right)\\ \mathbf{else}:\\ \;\;\;\;t\_1\\ \end{array} \end{array} \]
                                (FPCore (x y z)
                                 :precision binary64
                                 (let* ((t_0 (- (/ 2.0 3.0) z)) (t_1 (* (* (- y x) z) -6.0)))
                                   (if (<= t_0 -10000000000.0)
                                     t_1
                                     (if (<= t_0 1.0) (fma -3.0 x (* y 4.0)) t_1))))
                                double code(double x, double y, double z) {
                                	double t_0 = (2.0 / 3.0) - z;
                                	double t_1 = ((y - x) * z) * -6.0;
                                	double tmp;
                                	if (t_0 <= -10000000000.0) {
                                		tmp = t_1;
                                	} else if (t_0 <= 1.0) {
                                		tmp = fma(-3.0, x, (y * 4.0));
                                	} else {
                                		tmp = t_1;
                                	}
                                	return tmp;
                                }
                                
                                function code(x, y, z)
                                	t_0 = Float64(Float64(2.0 / 3.0) - z)
                                	t_1 = Float64(Float64(Float64(y - x) * z) * -6.0)
                                	tmp = 0.0
                                	if (t_0 <= -10000000000.0)
                                		tmp = t_1;
                                	elseif (t_0 <= 1.0)
                                		tmp = fma(-3.0, x, Float64(y * 4.0));
                                	else
                                		tmp = t_1;
                                	end
                                	return tmp
                                end
                                
                                code[x_, y_, z_] := Block[{t$95$0 = N[(N[(2.0 / 3.0), $MachinePrecision] - z), $MachinePrecision]}, Block[{t$95$1 = N[(N[(N[(y - x), $MachinePrecision] * z), $MachinePrecision] * -6.0), $MachinePrecision]}, If[LessEqual[t$95$0, -10000000000.0], t$95$1, If[LessEqual[t$95$0, 1.0], N[(-3.0 * x + N[(y * 4.0), $MachinePrecision]), $MachinePrecision], t$95$1]]]]
                                
                                \begin{array}{l}
                                
                                \\
                                \begin{array}{l}
                                t_0 := \frac{2}{3} - z\\
                                t_1 := \left(\left(y - x\right) \cdot z\right) \cdot -6\\
                                \mathbf{if}\;t\_0 \leq -10000000000:\\
                                \;\;\;\;t\_1\\
                                
                                \mathbf{elif}\;t\_0 \leq 1:\\
                                \;\;\;\;\mathsf{fma}\left(-3, x, y \cdot 4\right)\\
                                
                                \mathbf{else}:\\
                                \;\;\;\;t\_1\\
                                
                                
                                \end{array}
                                \end{array}
                                
                                Derivation
                                1. Split input into 2 regimes
                                2. if (-.f64 (/.f64 #s(literal 2 binary64) #s(literal 3 binary64)) z) < -1e10 or 1 < (-.f64 (/.f64 #s(literal 2 binary64) #s(literal 3 binary64)) z)

                                  1. Initial program 99.7%

                                    \[x + \left(\left(y - x\right) \cdot 6\right) \cdot \left(\frac{2}{3} - z\right) \]
                                  2. Add Preprocessing
                                  3. Taylor expanded in z around inf

                                    \[\leadsto \color{blue}{-6 \cdot \left(z \cdot \left(y - x\right)\right)} \]
                                  4. Step-by-step derivation
                                    1. *-commutativeN/A

                                      \[\leadsto \color{blue}{\left(z \cdot \left(y - x\right)\right) \cdot -6} \]
                                    2. lower-*.f64N/A

                                      \[\leadsto \color{blue}{\left(z \cdot \left(y - x\right)\right) \cdot -6} \]
                                    3. *-commutativeN/A

                                      \[\leadsto \color{blue}{\left(\left(y - x\right) \cdot z\right)} \cdot -6 \]
                                    4. lower-*.f64N/A

                                      \[\leadsto \color{blue}{\left(\left(y - x\right) \cdot z\right)} \cdot -6 \]
                                    5. lower--.f6498.9

                                      \[\leadsto \left(\color{blue}{\left(y - x\right)} \cdot z\right) \cdot -6 \]
                                  5. Applied rewrites98.9%

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

                                  if -1e10 < (-.f64 (/.f64 #s(literal 2 binary64) #s(literal 3 binary64)) z) < 1

                                  1. Initial program 98.8%

                                    \[x + \left(\left(y - x\right) \cdot 6\right) \cdot \left(\frac{2}{3} - z\right) \]
                                  2. Add Preprocessing
                                  3. Taylor expanded in z around 0

                                    \[\leadsto \color{blue}{x + 4 \cdot \left(y - x\right)} \]
                                  4. Step-by-step derivation
                                    1. +-commutativeN/A

                                      \[\leadsto \color{blue}{4 \cdot \left(y - x\right) + x} \]
                                    2. lower-fma.f64N/A

                                      \[\leadsto \color{blue}{\mathsf{fma}\left(4, y - x, x\right)} \]
                                    3. lower--.f6499.6

                                      \[\leadsto \mathsf{fma}\left(4, \color{blue}{y - x}, x\right) \]
                                  5. Applied rewrites99.6%

                                    \[\leadsto \color{blue}{\mathsf{fma}\left(4, y - x, x\right)} \]
                                  6. Taylor expanded in x around 0

                                    \[\leadsto -3 \cdot x + \color{blue}{4 \cdot y} \]
                                  7. Step-by-step derivation
                                    1. Applied rewrites99.7%

                                      \[\leadsto \mathsf{fma}\left(-3, \color{blue}{x}, 4 \cdot y\right) \]
                                  8. Recombined 2 regimes into one program.
                                  9. Final simplification99.3%

                                    \[\leadsto \begin{array}{l} \mathbf{if}\;\frac{2}{3} - z \leq -10000000000:\\ \;\;\;\;\left(\left(y - x\right) \cdot z\right) \cdot -6\\ \mathbf{elif}\;\frac{2}{3} - z \leq 1:\\ \;\;\;\;\mathsf{fma}\left(-3, x, y \cdot 4\right)\\ \mathbf{else}:\\ \;\;\;\;\left(\left(y - x\right) \cdot z\right) \cdot -6\\ \end{array} \]
                                  10. Add Preprocessing

                                  Alternative 11: 75.0% accurate, 0.6× speedup?

                                  \[\begin{array}{l} \\ \begin{array}{l} t_0 := \frac{2}{3} - z\\ t_1 := \left(6 \cdot z\right) \cdot x\\ \mathbf{if}\;t\_0 \leq -10000000000:\\ \;\;\;\;t\_1\\ \mathbf{elif}\;t\_0 \leq 2000000000:\\ \;\;\;\;\mathsf{fma}\left(4, y - x, x\right)\\ \mathbf{else}:\\ \;\;\;\;t\_1\\ \end{array} \end{array} \]
                                  (FPCore (x y z)
                                   :precision binary64
                                   (let* ((t_0 (- (/ 2.0 3.0) z)) (t_1 (* (* 6.0 z) x)))
                                     (if (<= t_0 -10000000000.0)
                                       t_1
                                       (if (<= t_0 2000000000.0) (fma 4.0 (- y x) x) t_1))))
                                  double code(double x, double y, double z) {
                                  	double t_0 = (2.0 / 3.0) - z;
                                  	double t_1 = (6.0 * z) * x;
                                  	double tmp;
                                  	if (t_0 <= -10000000000.0) {
                                  		tmp = t_1;
                                  	} else if (t_0 <= 2000000000.0) {
                                  		tmp = fma(4.0, (y - x), x);
                                  	} else {
                                  		tmp = t_1;
                                  	}
                                  	return tmp;
                                  }
                                  
                                  function code(x, y, z)
                                  	t_0 = Float64(Float64(2.0 / 3.0) - z)
                                  	t_1 = Float64(Float64(6.0 * z) * x)
                                  	tmp = 0.0
                                  	if (t_0 <= -10000000000.0)
                                  		tmp = t_1;
                                  	elseif (t_0 <= 2000000000.0)
                                  		tmp = fma(4.0, Float64(y - x), x);
                                  	else
                                  		tmp = t_1;
                                  	end
                                  	return tmp
                                  end
                                  
                                  code[x_, y_, z_] := Block[{t$95$0 = N[(N[(2.0 / 3.0), $MachinePrecision] - z), $MachinePrecision]}, Block[{t$95$1 = N[(N[(6.0 * z), $MachinePrecision] * x), $MachinePrecision]}, If[LessEqual[t$95$0, -10000000000.0], t$95$1, If[LessEqual[t$95$0, 2000000000.0], N[(4.0 * N[(y - x), $MachinePrecision] + x), $MachinePrecision], t$95$1]]]]
                                  
                                  \begin{array}{l}
                                  
                                  \\
                                  \begin{array}{l}
                                  t_0 := \frac{2}{3} - z\\
                                  t_1 := \left(6 \cdot z\right) \cdot x\\
                                  \mathbf{if}\;t\_0 \leq -10000000000:\\
                                  \;\;\;\;t\_1\\
                                  
                                  \mathbf{elif}\;t\_0 \leq 2000000000:\\
                                  \;\;\;\;\mathsf{fma}\left(4, y - x, x\right)\\
                                  
                                  \mathbf{else}:\\
                                  \;\;\;\;t\_1\\
                                  
                                  
                                  \end{array}
                                  \end{array}
                                  
                                  Derivation
                                  1. Split input into 2 regimes
                                  2. if (-.f64 (/.f64 #s(literal 2 binary64) #s(literal 3 binary64)) z) < -1e10 or 2e9 < (-.f64 (/.f64 #s(literal 2 binary64) #s(literal 3 binary64)) z)

                                    1. Initial program 99.7%

                                      \[x + \left(\left(y - x\right) \cdot 6\right) \cdot \left(\frac{2}{3} - z\right) \]
                                    2. Add Preprocessing
                                    3. Step-by-step derivation
                                      1. lift-+.f64N/A

                                        \[\leadsto \color{blue}{x + \left(\left(y - x\right) \cdot 6\right) \cdot \left(\frac{2}{3} - z\right)} \]
                                      2. +-commutativeN/A

                                        \[\leadsto \color{blue}{\left(\left(y - x\right) \cdot 6\right) \cdot \left(\frac{2}{3} - z\right) + x} \]
                                      3. lift-*.f64N/A

                                        \[\leadsto \color{blue}{\left(\left(y - x\right) \cdot 6\right) \cdot \left(\frac{2}{3} - z\right)} + x \]
                                      4. lift-*.f64N/A

                                        \[\leadsto \color{blue}{\left(\left(y - x\right) \cdot 6\right)} \cdot \left(\frac{2}{3} - z\right) + x \]
                                      5. associate-*l*N/A

                                        \[\leadsto \color{blue}{\left(y - x\right) \cdot \left(6 \cdot \left(\frac{2}{3} - z\right)\right)} + x \]
                                      6. *-commutativeN/A

                                        \[\leadsto \color{blue}{\left(6 \cdot \left(\frac{2}{3} - z\right)\right) \cdot \left(y - x\right)} + x \]
                                      7. lower-fma.f64N/A

                                        \[\leadsto \color{blue}{\mathsf{fma}\left(6 \cdot \left(\frac{2}{3} - z\right), y - x, x\right)} \]
                                      8. lift--.f64N/A

                                        \[\leadsto \mathsf{fma}\left(6 \cdot \color{blue}{\left(\frac{2}{3} - z\right)}, y - x, x\right) \]
                                      9. sub-negN/A

                                        \[\leadsto \mathsf{fma}\left(6 \cdot \color{blue}{\left(\frac{2}{3} + \left(\mathsf{neg}\left(z\right)\right)\right)}, y - x, x\right) \]
                                      10. +-commutativeN/A

                                        \[\leadsto \mathsf{fma}\left(6 \cdot \color{blue}{\left(\left(\mathsf{neg}\left(z\right)\right) + \frac{2}{3}\right)}, y - x, x\right) \]
                                      11. distribute-lft-inN/A

                                        \[\leadsto \mathsf{fma}\left(\color{blue}{6 \cdot \left(\mathsf{neg}\left(z\right)\right) + 6 \cdot \frac{2}{3}}, y - x, x\right) \]
                                      12. neg-mul-1N/A

                                        \[\leadsto \mathsf{fma}\left(6 \cdot \color{blue}{\left(-1 \cdot z\right)} + 6 \cdot \frac{2}{3}, y - x, x\right) \]
                                      13. associate-*r*N/A

                                        \[\leadsto \mathsf{fma}\left(\color{blue}{\left(6 \cdot -1\right) \cdot z} + 6 \cdot \frac{2}{3}, y - x, x\right) \]
                                      14. metadata-evalN/A

                                        \[\leadsto \mathsf{fma}\left(\color{blue}{-6} \cdot z + 6 \cdot \frac{2}{3}, y - x, x\right) \]
                                      15. metadata-evalN/A

                                        \[\leadsto \mathsf{fma}\left(\color{blue}{\left(\mathsf{neg}\left(6\right)\right)} \cdot z + 6 \cdot \frac{2}{3}, y - x, x\right) \]
                                      16. lower-fma.f64N/A

                                        \[\leadsto \mathsf{fma}\left(\color{blue}{\mathsf{fma}\left(\mathsf{neg}\left(6\right), z, 6 \cdot \frac{2}{3}\right)}, y - x, x\right) \]
                                      17. metadata-evalN/A

                                        \[\leadsto \mathsf{fma}\left(\mathsf{fma}\left(\color{blue}{-6}, z, 6 \cdot \frac{2}{3}\right), y - x, x\right) \]
                                      18. lift-/.f64N/A

                                        \[\leadsto \mathsf{fma}\left(\mathsf{fma}\left(-6, z, 6 \cdot \color{blue}{\frac{2}{3}}\right), y - x, x\right) \]
                                      19. metadata-evalN/A

                                        \[\leadsto \mathsf{fma}\left(\mathsf{fma}\left(-6, z, 6 \cdot \color{blue}{\frac{2}{3}}\right), y - x, x\right) \]
                                      20. metadata-eval99.8

                                        \[\leadsto \mathsf{fma}\left(\mathsf{fma}\left(-6, z, \color{blue}{4}\right), y - x, x\right) \]
                                    4. Applied rewrites99.8%

                                      \[\leadsto \color{blue}{\mathsf{fma}\left(\mathsf{fma}\left(-6, z, 4\right), y - x, x\right)} \]
                                    5. Taylor expanded in z around inf

                                      \[\leadsto \color{blue}{-6 \cdot \left(z \cdot \left(y - x\right)\right)} \]
                                    6. Step-by-step derivation
                                      1. *-commutativeN/A

                                        \[\leadsto -6 \cdot \color{blue}{\left(\left(y - x\right) \cdot z\right)} \]
                                      2. associate-*r*N/A

                                        \[\leadsto \color{blue}{\left(-6 \cdot \left(y - x\right)\right) \cdot z} \]
                                      3. lower-*.f64N/A

                                        \[\leadsto \color{blue}{\left(-6 \cdot \left(y - x\right)\right) \cdot z} \]
                                      4. lower-*.f64N/A

                                        \[\leadsto \color{blue}{\left(-6 \cdot \left(y - x\right)\right)} \cdot z \]
                                      5. lower--.f6499.2

                                        \[\leadsto \left(-6 \cdot \color{blue}{\left(y - x\right)}\right) \cdot z \]
                                    7. Applied rewrites99.2%

                                      \[\leadsto \color{blue}{\left(-6 \cdot \left(y - x\right)\right) \cdot z} \]
                                    8. Taylor expanded in x around inf

                                      \[\leadsto 6 \cdot \color{blue}{\left(x \cdot z\right)} \]
                                    9. Step-by-step derivation
                                      1. Applied rewrites59.5%

                                        \[\leadsto \left(6 \cdot z\right) \cdot \color{blue}{x} \]

                                      if -1e10 < (-.f64 (/.f64 #s(literal 2 binary64) #s(literal 3 binary64)) z) < 2e9

                                      1. Initial program 98.9%

                                        \[x + \left(\left(y - x\right) \cdot 6\right) \cdot \left(\frac{2}{3} - z\right) \]
                                      2. Add Preprocessing
                                      3. Taylor expanded in z around 0

                                        \[\leadsto \color{blue}{x + 4 \cdot \left(y - x\right)} \]
                                      4. Step-by-step derivation
                                        1. +-commutativeN/A

                                          \[\leadsto \color{blue}{4 \cdot \left(y - x\right) + x} \]
                                        2. lower-fma.f64N/A

                                          \[\leadsto \color{blue}{\mathsf{fma}\left(4, y - x, x\right)} \]
                                        3. lower--.f6499.0

                                          \[\leadsto \mathsf{fma}\left(4, \color{blue}{y - x}, x\right) \]
                                      5. Applied rewrites99.0%

                                        \[\leadsto \color{blue}{\mathsf{fma}\left(4, y - x, x\right)} \]
                                    10. Recombined 2 regimes into one program.
                                    11. Add Preprocessing

                                    Alternative 12: 38.8% accurate, 1.7× speedup?

                                    \[\begin{array}{l} \\ \begin{array}{l} \mathbf{if}\;x \leq -1.08 \cdot 10^{+21}:\\ \;\;\;\;-3 \cdot x\\ \mathbf{elif}\;x \leq 3.5 \cdot 10^{-67}:\\ \;\;\;\;y \cdot 4\\ \mathbf{else}:\\ \;\;\;\;-3 \cdot x\\ \end{array} \end{array} \]
                                    (FPCore (x y z)
                                     :precision binary64
                                     (if (<= x -1.08e+21) (* -3.0 x) (if (<= x 3.5e-67) (* y 4.0) (* -3.0 x))))
                                    double code(double x, double y, double z) {
                                    	double tmp;
                                    	if (x <= -1.08e+21) {
                                    		tmp = -3.0 * x;
                                    	} else if (x <= 3.5e-67) {
                                    		tmp = y * 4.0;
                                    	} else {
                                    		tmp = -3.0 * x;
                                    	}
                                    	return tmp;
                                    }
                                    
                                    real(8) function code(x, y, z)
                                        real(8), intent (in) :: x
                                        real(8), intent (in) :: y
                                        real(8), intent (in) :: z
                                        real(8) :: tmp
                                        if (x <= (-1.08d+21)) then
                                            tmp = (-3.0d0) * x
                                        else if (x <= 3.5d-67) then
                                            tmp = y * 4.0d0
                                        else
                                            tmp = (-3.0d0) * x
                                        end if
                                        code = tmp
                                    end function
                                    
                                    public static double code(double x, double y, double z) {
                                    	double tmp;
                                    	if (x <= -1.08e+21) {
                                    		tmp = -3.0 * x;
                                    	} else if (x <= 3.5e-67) {
                                    		tmp = y * 4.0;
                                    	} else {
                                    		tmp = -3.0 * x;
                                    	}
                                    	return tmp;
                                    }
                                    
                                    def code(x, y, z):
                                    	tmp = 0
                                    	if x <= -1.08e+21:
                                    		tmp = -3.0 * x
                                    	elif x <= 3.5e-67:
                                    		tmp = y * 4.0
                                    	else:
                                    		tmp = -3.0 * x
                                    	return tmp
                                    
                                    function code(x, y, z)
                                    	tmp = 0.0
                                    	if (x <= -1.08e+21)
                                    		tmp = Float64(-3.0 * x);
                                    	elseif (x <= 3.5e-67)
                                    		tmp = Float64(y * 4.0);
                                    	else
                                    		tmp = Float64(-3.0 * x);
                                    	end
                                    	return tmp
                                    end
                                    
                                    function tmp_2 = code(x, y, z)
                                    	tmp = 0.0;
                                    	if (x <= -1.08e+21)
                                    		tmp = -3.0 * x;
                                    	elseif (x <= 3.5e-67)
                                    		tmp = y * 4.0;
                                    	else
                                    		tmp = -3.0 * x;
                                    	end
                                    	tmp_2 = tmp;
                                    end
                                    
                                    code[x_, y_, z_] := If[LessEqual[x, -1.08e+21], N[(-3.0 * x), $MachinePrecision], If[LessEqual[x, 3.5e-67], N[(y * 4.0), $MachinePrecision], N[(-3.0 * x), $MachinePrecision]]]
                                    
                                    \begin{array}{l}
                                    
                                    \\
                                    \begin{array}{l}
                                    \mathbf{if}\;x \leq -1.08 \cdot 10^{+21}:\\
                                    \;\;\;\;-3 \cdot x\\
                                    
                                    \mathbf{elif}\;x \leq 3.5 \cdot 10^{-67}:\\
                                    \;\;\;\;y \cdot 4\\
                                    
                                    \mathbf{else}:\\
                                    \;\;\;\;-3 \cdot x\\
                                    
                                    
                                    \end{array}
                                    \end{array}
                                    
                                    Derivation
                                    1. Split input into 2 regimes
                                    2. if x < -1.08e21 or 3.5e-67 < x

                                      1. Initial program 99.6%

                                        \[x + \left(\left(y - x\right) \cdot 6\right) \cdot \left(\frac{2}{3} - z\right) \]
                                      2. Add Preprocessing
                                      3. Taylor expanded in z around 0

                                        \[\leadsto \color{blue}{x + 4 \cdot \left(y - x\right)} \]
                                      4. Step-by-step derivation
                                        1. +-commutativeN/A

                                          \[\leadsto \color{blue}{4 \cdot \left(y - x\right) + x} \]
                                        2. lower-fma.f64N/A

                                          \[\leadsto \color{blue}{\mathsf{fma}\left(4, y - x, x\right)} \]
                                        3. lower--.f6457.2

                                          \[\leadsto \mathsf{fma}\left(4, \color{blue}{y - x}, x\right) \]
                                      5. Applied rewrites57.2%

                                        \[\leadsto \color{blue}{\mathsf{fma}\left(4, y - x, x\right)} \]
                                      6. Taylor expanded in x around inf

                                        \[\leadsto -3 \cdot \color{blue}{x} \]
                                      7. Step-by-step derivation
                                        1. Applied rewrites44.4%

                                          \[\leadsto -3 \cdot \color{blue}{x} \]

                                        if -1.08e21 < x < 3.5e-67

                                        1. Initial program 98.8%

                                          \[x + \left(\left(y - x\right) \cdot 6\right) \cdot \left(\frac{2}{3} - z\right) \]
                                        2. Add Preprocessing
                                        3. Taylor expanded in z around 0

                                          \[\leadsto \color{blue}{x + 4 \cdot \left(y - x\right)} \]
                                        4. Step-by-step derivation
                                          1. +-commutativeN/A

                                            \[\leadsto \color{blue}{4 \cdot \left(y - x\right) + x} \]
                                          2. lower-fma.f64N/A

                                            \[\leadsto \color{blue}{\mathsf{fma}\left(4, y - x, x\right)} \]
                                          3. lower--.f6456.1

                                            \[\leadsto \mathsf{fma}\left(4, \color{blue}{y - x}, x\right) \]
                                        5. Applied rewrites56.1%

                                          \[\leadsto \color{blue}{\mathsf{fma}\left(4, y - x, x\right)} \]
                                        6. Taylor expanded in x around 0

                                          \[\leadsto 4 \cdot \color{blue}{y} \]
                                        7. Step-by-step derivation
                                          1. Applied rewrites45.6%

                                            \[\leadsto 4 \cdot \color{blue}{y} \]
                                        8. Recombined 2 regimes into one program.
                                        9. Final simplification44.9%

                                          \[\leadsto \begin{array}{l} \mathbf{if}\;x \leq -1.08 \cdot 10^{+21}:\\ \;\;\;\;-3 \cdot x\\ \mathbf{elif}\;x \leq 3.5 \cdot 10^{-67}:\\ \;\;\;\;y \cdot 4\\ \mathbf{else}:\\ \;\;\;\;-3 \cdot x\\ \end{array} \]
                                        10. Add Preprocessing

                                        Alternative 13: 51.5% accurate, 3.1× speedup?

                                        \[\begin{array}{l} \\ \mathsf{fma}\left(4, y - x, x\right) \end{array} \]
                                        (FPCore (x y z) :precision binary64 (fma 4.0 (- y x) x))
                                        double code(double x, double y, double z) {
                                        	return fma(4.0, (y - x), x);
                                        }
                                        
                                        function code(x, y, z)
                                        	return fma(4.0, Float64(y - x), x)
                                        end
                                        
                                        code[x_, y_, z_] := N[(4.0 * N[(y - x), $MachinePrecision] + x), $MachinePrecision]
                                        
                                        \begin{array}{l}
                                        
                                        \\
                                        \mathsf{fma}\left(4, y - x, x\right)
                                        \end{array}
                                        
                                        Derivation
                                        1. Initial program 99.2%

                                          \[x + \left(\left(y - x\right) \cdot 6\right) \cdot \left(\frac{2}{3} - z\right) \]
                                        2. Add Preprocessing
                                        3. Taylor expanded in z around 0

                                          \[\leadsto \color{blue}{x + 4 \cdot \left(y - x\right)} \]
                                        4. Step-by-step derivation
                                          1. +-commutativeN/A

                                            \[\leadsto \color{blue}{4 \cdot \left(y - x\right) + x} \]
                                          2. lower-fma.f64N/A

                                            \[\leadsto \color{blue}{\mathsf{fma}\left(4, y - x, x\right)} \]
                                          3. lower--.f6456.7

                                            \[\leadsto \mathsf{fma}\left(4, \color{blue}{y - x}, x\right) \]
                                        5. Applied rewrites56.7%

                                          \[\leadsto \color{blue}{\mathsf{fma}\left(4, y - x, x\right)} \]
                                        6. Add Preprocessing

                                        Alternative 14: 26.5% accurate, 5.2× speedup?

                                        \[\begin{array}{l} \\ -3 \cdot x \end{array} \]
                                        (FPCore (x y z) :precision binary64 (* -3.0 x))
                                        double code(double x, double y, double z) {
                                        	return -3.0 * x;
                                        }
                                        
                                        real(8) function code(x, y, z)
                                            real(8), intent (in) :: x
                                            real(8), intent (in) :: y
                                            real(8), intent (in) :: z
                                            code = (-3.0d0) * x
                                        end function
                                        
                                        public static double code(double x, double y, double z) {
                                        	return -3.0 * x;
                                        }
                                        
                                        def code(x, y, z):
                                        	return -3.0 * x
                                        
                                        function code(x, y, z)
                                        	return Float64(-3.0 * x)
                                        end
                                        
                                        function tmp = code(x, y, z)
                                        	tmp = -3.0 * x;
                                        end
                                        
                                        code[x_, y_, z_] := N[(-3.0 * x), $MachinePrecision]
                                        
                                        \begin{array}{l}
                                        
                                        \\
                                        -3 \cdot x
                                        \end{array}
                                        
                                        Derivation
                                        1. Initial program 99.2%

                                          \[x + \left(\left(y - x\right) \cdot 6\right) \cdot \left(\frac{2}{3} - z\right) \]
                                        2. Add Preprocessing
                                        3. Taylor expanded in z around 0

                                          \[\leadsto \color{blue}{x + 4 \cdot \left(y - x\right)} \]
                                        4. Step-by-step derivation
                                          1. +-commutativeN/A

                                            \[\leadsto \color{blue}{4 \cdot \left(y - x\right) + x} \]
                                          2. lower-fma.f64N/A

                                            \[\leadsto \color{blue}{\mathsf{fma}\left(4, y - x, x\right)} \]
                                          3. lower--.f6456.7

                                            \[\leadsto \mathsf{fma}\left(4, \color{blue}{y - x}, x\right) \]
                                        5. Applied rewrites56.7%

                                          \[\leadsto \color{blue}{\mathsf{fma}\left(4, y - x, x\right)} \]
                                        6. Taylor expanded in x around inf

                                          \[\leadsto -3 \cdot \color{blue}{x} \]
                                        7. Step-by-step derivation
                                          1. Applied rewrites29.5%

                                            \[\leadsto -3 \cdot \color{blue}{x} \]
                                          2. Add Preprocessing

                                          Reproduce

                                          ?
                                          herbie shell --seed 2024304 
                                          (FPCore (x y z)
                                            :name "Data.Colour.RGBSpace.HSL:hsl from colour-2.3.3, D"
                                            :precision binary64
                                            (+ x (* (* (- y x) 6.0) (- (/ 2.0 3.0) z))))