Optimisation.CirclePacking:place from circle-packing-0.1.0.4, G

Percentage Accurate: 100.0% → 100.0%
Time: 4.0s
Alternatives: 6
Speedup: 1.0×

Specification

?
\[\begin{array}{l} \\ \left(x + y\right) \cdot \left(z + 1\right) \end{array} \]
(FPCore (x y z) :precision binary64 (* (+ x y) (+ z 1.0)))
double code(double x, double y, double z) {
	return (x + y) * (z + 1.0);
}
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) * (z + 1.0d0)
end function
public static double code(double x, double y, double z) {
	return (x + y) * (z + 1.0);
}
def code(x, y, z):
	return (x + y) * (z + 1.0)
function code(x, y, z)
	return Float64(Float64(x + y) * Float64(z + 1.0))
end
function tmp = code(x, y, z)
	tmp = (x + y) * (z + 1.0);
end
code[x_, y_, z_] := N[(N[(x + y), $MachinePrecision] * N[(z + 1.0), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}

\\
\left(x + y\right) \cdot \left(z + 1\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 6 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: 100.0% accurate, 1.0× speedup?

\[\begin{array}{l} \\ \left(x + y\right) \cdot \left(z + 1\right) \end{array} \]
(FPCore (x y z) :precision binary64 (* (+ x y) (+ z 1.0)))
double code(double x, double y, double z) {
	return (x + y) * (z + 1.0);
}
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) * (z + 1.0d0)
end function
public static double code(double x, double y, double z) {
	return (x + y) * (z + 1.0);
}
def code(x, y, z):
	return (x + y) * (z + 1.0)
function code(x, y, z)
	return Float64(Float64(x + y) * Float64(z + 1.0))
end
function tmp = code(x, y, z)
	tmp = (x + y) * (z + 1.0);
end
code[x_, y_, z_] := N[(N[(x + y), $MachinePrecision] * N[(z + 1.0), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}

\\
\left(x + y\right) \cdot \left(z + 1\right)
\end{array}

Alternative 1: 100.0% accurate, 1.0× speedup?

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

\\
\left(1 + z\right) \cdot \left(y + x\right)
\end{array}
Derivation
  1. Initial program 100.0%

    \[\left(x + y\right) \cdot \left(z + 1\right) \]
  2. Add Preprocessing
  3. Final simplification100.0%

    \[\leadsto \left(1 + z\right) \cdot \left(y + x\right) \]
  4. Add Preprocessing

Alternative 2: 41.3% accurate, 0.3× speedup?

\[\begin{array}{l} \\ \begin{array}{l} \mathbf{if}\;y + x \leq 5 \cdot 10^{-278}:\\ \;\;\;\;\mathsf{fma}\left(z, x, x\right)\\ \mathbf{elif}\;y + x \leq 2 \cdot 10^{-175}:\\ \;\;\;\;y + x\\ \mathbf{elif}\;y + x \leq 10^{+122}:\\ \;\;\;\;z \cdot y\\ \mathbf{elif}\;y + x \leq 2 \cdot 10^{+175}:\\ \;\;\;\;y + x\\ \mathbf{else}:\\ \;\;\;\;z \cdot y\\ \end{array} \end{array} \]
(FPCore (x y z)
 :precision binary64
 (if (<= (+ y x) 5e-278)
   (fma z x x)
   (if (<= (+ y x) 2e-175)
     (+ y x)
     (if (<= (+ y x) 1e+122)
       (* z y)
       (if (<= (+ y x) 2e+175) (+ y x) (* z y))))))
double code(double x, double y, double z) {
	double tmp;
	if ((y + x) <= 5e-278) {
		tmp = fma(z, x, x);
	} else if ((y + x) <= 2e-175) {
		tmp = y + x;
	} else if ((y + x) <= 1e+122) {
		tmp = z * y;
	} else if ((y + x) <= 2e+175) {
		tmp = y + x;
	} else {
		tmp = z * y;
	}
	return tmp;
}
function code(x, y, z)
	tmp = 0.0
	if (Float64(y + x) <= 5e-278)
		tmp = fma(z, x, x);
	elseif (Float64(y + x) <= 2e-175)
		tmp = Float64(y + x);
	elseif (Float64(y + x) <= 1e+122)
		tmp = Float64(z * y);
	elseif (Float64(y + x) <= 2e+175)
		tmp = Float64(y + x);
	else
		tmp = Float64(z * y);
	end
	return tmp
end
code[x_, y_, z_] := If[LessEqual[N[(y + x), $MachinePrecision], 5e-278], N[(z * x + x), $MachinePrecision], If[LessEqual[N[(y + x), $MachinePrecision], 2e-175], N[(y + x), $MachinePrecision], If[LessEqual[N[(y + x), $MachinePrecision], 1e+122], N[(z * y), $MachinePrecision], If[LessEqual[N[(y + x), $MachinePrecision], 2e+175], N[(y + x), $MachinePrecision], N[(z * y), $MachinePrecision]]]]]
\begin{array}{l}

\\
\begin{array}{l}
\mathbf{if}\;y + x \leq 5 \cdot 10^{-278}:\\
\;\;\;\;\mathsf{fma}\left(z, x, x\right)\\

\mathbf{elif}\;y + x \leq 2 \cdot 10^{-175}:\\
\;\;\;\;y + x\\

\mathbf{elif}\;y + x \leq 10^{+122}:\\
\;\;\;\;z \cdot y\\

\mathbf{elif}\;y + x \leq 2 \cdot 10^{+175}:\\
\;\;\;\;y + x\\

\mathbf{else}:\\
\;\;\;\;z \cdot y\\


\end{array}
\end{array}
Derivation
  1. Split input into 3 regimes
  2. if (+.f64 x y) < 4.99999999999999985e-278

    1. Initial program 100.0%

      \[\left(x + y\right) \cdot \left(z + 1\right) \]
    2. Add Preprocessing
    3. Taylor expanded in y around 0

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

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

        \[\leadsto \color{blue}{z \cdot x + 1 \cdot x} \]
      3. *-lft-identityN/A

        \[\leadsto z \cdot x + \color{blue}{x} \]
      4. lower-fma.f6447.0

        \[\leadsto \color{blue}{\mathsf{fma}\left(z, x, x\right)} \]
    5. Applied rewrites47.0%

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

    if 4.99999999999999985e-278 < (+.f64 x y) < 2e-175 or 1.00000000000000001e122 < (+.f64 x y) < 1.9999999999999999e175

    1. Initial program 100.0%

      \[\left(x + y\right) \cdot \left(z + 1\right) \]
    2. Add Preprocessing
    3. Taylor expanded in z around 0

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

        \[\leadsto \color{blue}{y + x} \]
      2. lower-+.f6479.8

        \[\leadsto \color{blue}{y + x} \]
    5. Applied rewrites79.8%

      \[\leadsto \color{blue}{y + x} \]

    if 2e-175 < (+.f64 x y) < 1.00000000000000001e122 or 1.9999999999999999e175 < (+.f64 x y)

    1. Initial program 100.0%

      \[\left(x + y\right) \cdot \left(z + 1\right) \]
    2. Add Preprocessing
    3. Taylor expanded in y around inf

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

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

        \[\leadsto \color{blue}{z \cdot y + 1 \cdot y} \]
      3. *-lft-identityN/A

        \[\leadsto z \cdot y + \color{blue}{y} \]
      4. lower-fma.f6448.1

        \[\leadsto \color{blue}{\mathsf{fma}\left(z, y, y\right)} \]
    5. Applied rewrites48.1%

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

      \[\leadsto y \cdot \color{blue}{z} \]
    7. Step-by-step derivation
      1. Applied rewrites30.3%

        \[\leadsto z \cdot \color{blue}{y} \]
    8. Recombined 3 regimes into one program.
    9. Final simplification42.7%

      \[\leadsto \begin{array}{l} \mathbf{if}\;y + x \leq 5 \cdot 10^{-278}:\\ \;\;\;\;\mathsf{fma}\left(z, x, x\right)\\ \mathbf{elif}\;y + x \leq 2 \cdot 10^{-175}:\\ \;\;\;\;y + x\\ \mathbf{elif}\;y + x \leq 10^{+122}:\\ \;\;\;\;z \cdot y\\ \mathbf{elif}\;y + x \leq 2 \cdot 10^{+175}:\\ \;\;\;\;y + x\\ \mathbf{else}:\\ \;\;\;\;z \cdot y\\ \end{array} \]
    10. Add Preprocessing

    Alternative 3: 74.5% accurate, 0.5× speedup?

    \[\begin{array}{l} \\ \begin{array}{l} \mathbf{if}\;1 + z \leq -5000:\\ \;\;\;\;z \cdot y\\ \mathbf{elif}\;1 + z \leq 200000:\\ \;\;\;\;y + x\\ \mathbf{else}:\\ \;\;\;\;z \cdot y\\ \end{array} \end{array} \]
    (FPCore (x y z)
     :precision binary64
     (if (<= (+ 1.0 z) -5000.0)
       (* z y)
       (if (<= (+ 1.0 z) 200000.0) (+ y x) (* z y))))
    double code(double x, double y, double z) {
    	double tmp;
    	if ((1.0 + z) <= -5000.0) {
    		tmp = z * y;
    	} else if ((1.0 + z) <= 200000.0) {
    		tmp = y + x;
    	} else {
    		tmp = z * y;
    	}
    	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 ((1.0d0 + z) <= (-5000.0d0)) then
            tmp = z * y
        else if ((1.0d0 + z) <= 200000.0d0) then
            tmp = y + x
        else
            tmp = z * y
        end if
        code = tmp
    end function
    
    public static double code(double x, double y, double z) {
    	double tmp;
    	if ((1.0 + z) <= -5000.0) {
    		tmp = z * y;
    	} else if ((1.0 + z) <= 200000.0) {
    		tmp = y + x;
    	} else {
    		tmp = z * y;
    	}
    	return tmp;
    }
    
    def code(x, y, z):
    	tmp = 0
    	if (1.0 + z) <= -5000.0:
    		tmp = z * y
    	elif (1.0 + z) <= 200000.0:
    		tmp = y + x
    	else:
    		tmp = z * y
    	return tmp
    
    function code(x, y, z)
    	tmp = 0.0
    	if (Float64(1.0 + z) <= -5000.0)
    		tmp = Float64(z * y);
    	elseif (Float64(1.0 + z) <= 200000.0)
    		tmp = Float64(y + x);
    	else
    		tmp = Float64(z * y);
    	end
    	return tmp
    end
    
    function tmp_2 = code(x, y, z)
    	tmp = 0.0;
    	if ((1.0 + z) <= -5000.0)
    		tmp = z * y;
    	elseif ((1.0 + z) <= 200000.0)
    		tmp = y + x;
    	else
    		tmp = z * y;
    	end
    	tmp_2 = tmp;
    end
    
    code[x_, y_, z_] := If[LessEqual[N[(1.0 + z), $MachinePrecision], -5000.0], N[(z * y), $MachinePrecision], If[LessEqual[N[(1.0 + z), $MachinePrecision], 200000.0], N[(y + x), $MachinePrecision], N[(z * y), $MachinePrecision]]]
    
    \begin{array}{l}
    
    \\
    \begin{array}{l}
    \mathbf{if}\;1 + z \leq -5000:\\
    \;\;\;\;z \cdot y\\
    
    \mathbf{elif}\;1 + z \leq 200000:\\
    \;\;\;\;y + x\\
    
    \mathbf{else}:\\
    \;\;\;\;z \cdot y\\
    
    
    \end{array}
    \end{array}
    
    Derivation
    1. Split input into 2 regimes
    2. if (+.f64 z #s(literal 1 binary64)) < -5e3 or 2e5 < (+.f64 z #s(literal 1 binary64))

      1. Initial program 100.0%

        \[\left(x + y\right) \cdot \left(z + 1\right) \]
      2. Add Preprocessing
      3. Taylor expanded in y around inf

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

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

          \[\leadsto \color{blue}{z \cdot y + 1 \cdot y} \]
        3. *-lft-identityN/A

          \[\leadsto z \cdot y + \color{blue}{y} \]
        4. lower-fma.f6454.8

          \[\leadsto \color{blue}{\mathsf{fma}\left(z, y, y\right)} \]
      5. Applied rewrites54.8%

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

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

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

        if -5e3 < (+.f64 z #s(literal 1 binary64)) < 2e5

        1. Initial program 100.0%

          \[\left(x + y\right) \cdot \left(z + 1\right) \]
        2. Add Preprocessing
        3. Taylor expanded in z around 0

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

            \[\leadsto \color{blue}{y + x} \]
          2. lower-+.f6496.5

            \[\leadsto \color{blue}{y + x} \]
        5. Applied rewrites96.5%

          \[\leadsto \color{blue}{y + x} \]
      8. Recombined 2 regimes into one program.
      9. Final simplification75.1%

        \[\leadsto \begin{array}{l} \mathbf{if}\;1 + z \leq -5000:\\ \;\;\;\;z \cdot y\\ \mathbf{elif}\;1 + z \leq 200000:\\ \;\;\;\;y + x\\ \mathbf{else}:\\ \;\;\;\;z \cdot y\\ \end{array} \]
      10. Add Preprocessing

      Alternative 4: 74.4% accurate, 0.5× speedup?

      \[\begin{array}{l} \\ \begin{array}{l} \mathbf{if}\;1 + z \leq -5000:\\ \;\;\;\;z \cdot x\\ \mathbf{elif}\;1 + z \leq 20:\\ \;\;\;\;y + x\\ \mathbf{else}:\\ \;\;\;\;z \cdot x\\ \end{array} \end{array} \]
      (FPCore (x y z)
       :precision binary64
       (if (<= (+ 1.0 z) -5000.0) (* z x) (if (<= (+ 1.0 z) 20.0) (+ y x) (* z x))))
      double code(double x, double y, double z) {
      	double tmp;
      	if ((1.0 + z) <= -5000.0) {
      		tmp = z * x;
      	} else if ((1.0 + z) <= 20.0) {
      		tmp = y + x;
      	} else {
      		tmp = z * 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 ((1.0d0 + z) <= (-5000.0d0)) then
              tmp = z * x
          else if ((1.0d0 + z) <= 20.0d0) then
              tmp = y + x
          else
              tmp = z * x
          end if
          code = tmp
      end function
      
      public static double code(double x, double y, double z) {
      	double tmp;
      	if ((1.0 + z) <= -5000.0) {
      		tmp = z * x;
      	} else if ((1.0 + z) <= 20.0) {
      		tmp = y + x;
      	} else {
      		tmp = z * x;
      	}
      	return tmp;
      }
      
      def code(x, y, z):
      	tmp = 0
      	if (1.0 + z) <= -5000.0:
      		tmp = z * x
      	elif (1.0 + z) <= 20.0:
      		tmp = y + x
      	else:
      		tmp = z * x
      	return tmp
      
      function code(x, y, z)
      	tmp = 0.0
      	if (Float64(1.0 + z) <= -5000.0)
      		tmp = Float64(z * x);
      	elseif (Float64(1.0 + z) <= 20.0)
      		tmp = Float64(y + x);
      	else
      		tmp = Float64(z * x);
      	end
      	return tmp
      end
      
      function tmp_2 = code(x, y, z)
      	tmp = 0.0;
      	if ((1.0 + z) <= -5000.0)
      		tmp = z * x;
      	elseif ((1.0 + z) <= 20.0)
      		tmp = y + x;
      	else
      		tmp = z * x;
      	end
      	tmp_2 = tmp;
      end
      
      code[x_, y_, z_] := If[LessEqual[N[(1.0 + z), $MachinePrecision], -5000.0], N[(z * x), $MachinePrecision], If[LessEqual[N[(1.0 + z), $MachinePrecision], 20.0], N[(y + x), $MachinePrecision], N[(z * x), $MachinePrecision]]]
      
      \begin{array}{l}
      
      \\
      \begin{array}{l}
      \mathbf{if}\;1 + z \leq -5000:\\
      \;\;\;\;z \cdot x\\
      
      \mathbf{elif}\;1 + z \leq 20:\\
      \;\;\;\;y + x\\
      
      \mathbf{else}:\\
      \;\;\;\;z \cdot x\\
      
      
      \end{array}
      \end{array}
      
      Derivation
      1. Split input into 2 regimes
      2. if (+.f64 z #s(literal 1 binary64)) < -5e3 or 20 < (+.f64 z #s(literal 1 binary64))

        1. Initial program 100.0%

          \[\left(x + y\right) \cdot \left(z + 1\right) \]
        2. Add Preprocessing
        3. Taylor expanded in y around 0

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

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

            \[\leadsto \color{blue}{z \cdot x + 1 \cdot x} \]
          3. *-lft-identityN/A

            \[\leadsto z \cdot x + \color{blue}{x} \]
          4. lower-fma.f6451.3

            \[\leadsto \color{blue}{\mathsf{fma}\left(z, x, x\right)} \]
        5. Applied rewrites51.3%

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

          \[\leadsto x \cdot \color{blue}{z} \]
        7. Step-by-step derivation
          1. Applied rewrites50.8%

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

          if -5e3 < (+.f64 z #s(literal 1 binary64)) < 20

          1. Initial program 100.0%

            \[\left(x + y\right) \cdot \left(z + 1\right) \]
          2. Add Preprocessing
          3. Taylor expanded in z around 0

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

              \[\leadsto \color{blue}{y + x} \]
            2. lower-+.f6497.2

              \[\leadsto \color{blue}{y + x} \]
          5. Applied rewrites97.2%

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

          \[\leadsto \begin{array}{l} \mathbf{if}\;1 + z \leq -5000:\\ \;\;\;\;z \cdot x\\ \mathbf{elif}\;1 + z \leq 20:\\ \;\;\;\;y + x\\ \mathbf{else}:\\ \;\;\;\;z \cdot x\\ \end{array} \]
        10. Add Preprocessing

        Alternative 5: 51.4% accurate, 0.7× speedup?

        \[\begin{array}{l} \\ \begin{array}{l} \mathbf{if}\;y + x \leq -5 \cdot 10^{-265}:\\ \;\;\;\;\mathsf{fma}\left(z, x, x\right)\\ \mathbf{else}:\\ \;\;\;\;\mathsf{fma}\left(z, y, y\right)\\ \end{array} \end{array} \]
        (FPCore (x y z)
         :precision binary64
         (if (<= (+ y x) -5e-265) (fma z x x) (fma z y y)))
        double code(double x, double y, double z) {
        	double tmp;
        	if ((y + x) <= -5e-265) {
        		tmp = fma(z, x, x);
        	} else {
        		tmp = fma(z, y, y);
        	}
        	return tmp;
        }
        
        function code(x, y, z)
        	tmp = 0.0
        	if (Float64(y + x) <= -5e-265)
        		tmp = fma(z, x, x);
        	else
        		tmp = fma(z, y, y);
        	end
        	return tmp
        end
        
        code[x_, y_, z_] := If[LessEqual[N[(y + x), $MachinePrecision], -5e-265], N[(z * x + x), $MachinePrecision], N[(z * y + y), $MachinePrecision]]
        
        \begin{array}{l}
        
        \\
        \begin{array}{l}
        \mathbf{if}\;y + x \leq -5 \cdot 10^{-265}:\\
        \;\;\;\;\mathsf{fma}\left(z, x, x\right)\\
        
        \mathbf{else}:\\
        \;\;\;\;\mathsf{fma}\left(z, y, y\right)\\
        
        
        \end{array}
        \end{array}
        
        Derivation
        1. Split input into 2 regimes
        2. if (+.f64 x y) < -5.0000000000000001e-265

          1. Initial program 100.0%

            \[\left(x + y\right) \cdot \left(z + 1\right) \]
          2. Add Preprocessing
          3. Taylor expanded in y around 0

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

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

              \[\leadsto \color{blue}{z \cdot x + 1 \cdot x} \]
            3. *-lft-identityN/A

              \[\leadsto z \cdot x + \color{blue}{x} \]
            4. lower-fma.f6447.2

              \[\leadsto \color{blue}{\mathsf{fma}\left(z, x, x\right)} \]
          5. Applied rewrites47.2%

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

          if -5.0000000000000001e-265 < (+.f64 x y)

          1. Initial program 100.0%

            \[\left(x + y\right) \cdot \left(z + 1\right) \]
          2. Add Preprocessing
          3. Taylor expanded in y around inf

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

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

              \[\leadsto \color{blue}{z \cdot y + 1 \cdot y} \]
            3. *-lft-identityN/A

              \[\leadsto z \cdot y + \color{blue}{y} \]
            4. lower-fma.f6447.9

              \[\leadsto \color{blue}{\mathsf{fma}\left(z, y, y\right)} \]
          5. Applied rewrites47.9%

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

          \[\leadsto \begin{array}{l} \mathbf{if}\;y + x \leq -5 \cdot 10^{-265}:\\ \;\;\;\;\mathsf{fma}\left(z, x, x\right)\\ \mathbf{else}:\\ \;\;\;\;\mathsf{fma}\left(z, y, y\right)\\ \end{array} \]
        5. Add Preprocessing

        Alternative 6: 50.2% accurate, 3.0× speedup?

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

          \[\left(x + y\right) \cdot \left(z + 1\right) \]
        2. Add Preprocessing
        3. Taylor expanded in z around 0

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

            \[\leadsto \color{blue}{y + x} \]
          2. lower-+.f6449.4

            \[\leadsto \color{blue}{y + x} \]
        5. Applied rewrites49.4%

          \[\leadsto \color{blue}{y + x} \]
        6. Add Preprocessing

        Reproduce

        ?
        herbie shell --seed 2024264 
        (FPCore (x y z)
          :name "Optimisation.CirclePacking:place from circle-packing-0.1.0.4, G"
          :precision binary64
          (* (+ x y) (+ z 1.0)))