Diagrams.TwoD.Apollonian:initialConfig from diagrams-contrib-1.3.0.5, A

Percentage Accurate: 68.7% → 99.9%

Time: 9.5s

Alternatives: 9

Speedup: 1.2×

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

Specification

Math

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

FPCore

(FPCore (x y z)
 :precision binary64
 (/ (- (+ (* x x) (* y y)) (* z z)) (* y 2.0)))

C

double code(double x, double y, double z) {
	return (((x * x) + (y * y)) - (z * z)) / (y * 2.0);
}

Fortran

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

Java

public static double code(double x, double y, double z) {
	return (((x * x) + (y * y)) - (z * z)) / (y * 2.0);
}

Python

def code(x, y, z):
	return (((x * x) + (y * y)) - (z * z)) / (y * 2.0)

Julia

function code(x, y, z)
	return Float64(Float64(Float64(Float64(x * x) + Float64(y * y)) - Float64(z * z)) / Float64(y * 2.0))
end

MATLAB

function tmp = code(x, y, z)
	tmp = (((x * x) + (y * y)) - (z * z)) / (y * 2.0);
end

Wolfram

code[x_, y_, z_] := N[(N[(N[(N[(x * x), $MachinePrecision] + N[(y * y), $MachinePrecision]), $MachinePrecision] - N[(z * z), $MachinePrecision]), $MachinePrecision] / N[(y * 2.0), $MachinePrecision]), $MachinePrecision]

Initial Program: 68.7% accurate, 1.0× speedup

Math

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

FPCore

(FPCore (x y z)
 :precision binary64
 (/ (- (+ (* x x) (* y y)) (* z z)) (* y 2.0)))

C

double code(double x, double y, double z) {
	return (((x * x) + (y * y)) - (z * z)) / (y * 2.0);
}

Fortran

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

Java

public static double code(double x, double y, double z) {
	return (((x * x) + (y * y)) - (z * z)) / (y * 2.0);
}

Python

def code(x, y, z):
	return (((x * x) + (y * y)) - (z * z)) / (y * 2.0)

Julia

function code(x, y, z)
	return Float64(Float64(Float64(Float64(x * x) + Float64(y * y)) - Float64(z * z)) / Float64(y * 2.0))
end

MATLAB

function tmp = code(x, y, z)
	tmp = (((x * x) + (y * y)) - (z * z)) / (y * 2.0);
end

Wolfram

code[x_, y_, z_] := N[(N[(N[(N[(x * x), $MachinePrecision] + N[(y * y), $MachinePrecision]), $MachinePrecision] - N[(z * z), $MachinePrecision]), $MachinePrecision] / N[(y * 2.0), $MachinePrecision]), $MachinePrecision]

Alternative 1: 99.9% accurate, 1.2× speedup

Math

\[\begin{array}{l} \\ 0.5 \cdot \left(y + \frac{x + z}{\frac{y}{x - z}}\right) \end{array} \]

FPCore

(FPCore (x y z) :precision binary64 (* 0.5 (+ y (/ (+ x z) (/ y (- x z))))))

C

double code(double x, double y, double z) {
	return 0.5 * (y + ((x + z) / (y / (x - z))));
}

Fortran

real(8) function code(x, y, z)
    real(8), intent (in) :: x
    real(8), intent (in) :: y
    real(8), intent (in) :: z
    code = 0.5d0 * (y + ((x + z) / (y / (x - z))))
end function

Java

public static double code(double x, double y, double z) {
	return 0.5 * (y + ((x + z) / (y / (x - z))));
}

Python

def code(x, y, z):
	return 0.5 * (y + ((x + z) / (y / (x - z))))

Julia

function code(x, y, z)
	return Float64(0.5 * Float64(y + Float64(Float64(x + z) / Float64(y / Float64(x - z)))))
end

MATLAB

function tmp = code(x, y, z)
	tmp = 0.5 * (y + ((x + z) / (y / (x - z))));
end

Wolfram

code[x_, y_, z_] := N[(0.5 * N[(y + N[(N[(x + z), $MachinePrecision] / N[(y / N[(x - z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]

Derivation

1. Initial program 67.4%

   \[\frac{\left(x \cdot x + y \cdot y\right) - z \cdot z}{y \cdot 2} \]
2. Step-by-step derivation

   1. remove-double-neg 67.4%

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

   2. distribute-lft-neg-out 67.4%

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

   3. distribute-frac-neg2 67.4%

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

   4. distribute-frac-neg 67.4%

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

   5. neg-mul-1 67.4%

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

   6. distribute-lft-neg-out 67.4%

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

   7. *-commutative 67.4%

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

   8. distribute-lft-neg-in 67.4%

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

   9. times-frac 67.4%

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

   10. metadata-eval 67.4%

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

   11. metadata-eval 67.4%

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

   12. associate--l+ 67.4%

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

   13. fma-define 70.1%

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

3. Simplified 70.1%

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

5. Taylor expanded in x around 0 78.8%

   \[\leadsto 0.5 \cdot \color{blue}{\left(\left(y + \frac{{x}^{2}}{y}\right) - \frac{{z}^{2}}{y}\right)} \]
6. Step-by-step derivation

   1. associate--l+ 78.8%

      \[\leadsto 0.5 \cdot \color{blue}{\left(y + \left(\frac{{x}^{2}}{y} - \frac{{z}^{2}}{y}\right)\right)} \]

   2. div-sub 83.8%

      \[\leadsto 0.5 \cdot \left(y + \color{blue}{\frac{{x}^{2} - {z}^{2}}{y}}\right) \]

7. Simplified 83.8%

   \[\leadsto 0.5 \cdot \color{blue}{\left(y + \frac{{x}^{2} - {z}^{2}}{y}\right)} \]
8. Step-by-step derivation

   1. pow2 83.8%

      \[\leadsto 0.5 \cdot \left(y + \frac{\color{blue}{x \cdot x} - {z}^{2}}{y}\right) \]

   2. pow2 83.8%

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

   3. difference-of-squares 89.1%

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

9. Applied egg-rr 89.1%

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

    1. associate-/l* 99.9%

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

    2. *-commutative 99.9%

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

11. Applied egg-rr 99.9%

    \[\leadsto 0.5 \cdot \left(y + \color{blue}{\frac{x - z}{y} \cdot \left(x + z\right)}\right) \]
12. Step-by-step derivation

    1. *-commutative 99.9%

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

    2. clear-num 99.9%

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

    3. un-div-inv 99.9%

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

13. Applied egg-rr 99.9%

    \[\leadsto 0.5 \cdot \left(y + \color{blue}{\frac{x + z}{\frac{y}{x - z}}}\right) \]
14. Add Preprocessing

Alternative 2: 80.3% accurate, 0.9× speedup

Math

\[\begin{array}{l} \\ \begin{array}{l} \mathbf{if}\;z \leq 10^{+29}:\\ \;\;\;\;0.5 \cdot \left(y + \left(x + z\right) \cdot \frac{x}{y}\right)\\ \mathbf{else}:\\ \;\;\;\;0.5 \cdot \left(y + z \cdot \frac{x - z}{y}\right)\\ \end{array} \end{array} \]

FPCore

(FPCore (x y z)
 :precision binary64
 (if (<= z 1e+29)
   (* 0.5 (+ y (* (+ x z) (/ x y))))
   (* 0.5 (+ y (* z (/ (- x z) y))))))

C

double code(double x, double y, double z) {
	double tmp;
	if (z <= 1e+29) {
		tmp = 0.5 * (y + ((x + z) * (x / y)));
	} else {
		tmp = 0.5 * (y + (z * ((x - z) / y)));
	}
	return tmp;
}

Fortran

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 (z <= 1d+29) then
        tmp = 0.5d0 * (y + ((x + z) * (x / y)))
    else
        tmp = 0.5d0 * (y + (z * ((x - z) / y)))
    end if
    code = tmp
end function

Java

public static double code(double x, double y, double z) {
	double tmp;
	if (z <= 1e+29) {
		tmp = 0.5 * (y + ((x + z) * (x / y)));
	} else {
		tmp = 0.5 * (y + (z * ((x - z) / y)));
	}
	return tmp;
}

Python

def code(x, y, z):
	tmp = 0
	if z <= 1e+29:
		tmp = 0.5 * (y + ((x + z) * (x / y)))
	else:
		tmp = 0.5 * (y + (z * ((x - z) / y)))
	return tmp

Julia

function code(x, y, z)
	tmp = 0.0
	if (z <= 1e+29)
		tmp = Float64(0.5 * Float64(y + Float64(Float64(x + z) * Float64(x / y))));
	else
		tmp = Float64(0.5 * Float64(y + Float64(z * Float64(Float64(x - z) / y))));
	end
	return tmp
end

MATLAB

function tmp_2 = code(x, y, z)
	tmp = 0.0;
	if (z <= 1e+29)
		tmp = 0.5 * (y + ((x + z) * (x / y)));
	else
		tmp = 0.5 * (y + (z * ((x - z) / y)));
	end
	tmp_2 = tmp;
end

Wolfram

code[x_, y_, z_] := If[LessEqual[z, 1e+29], N[(0.5 * N[(y + N[(N[(x + z), $MachinePrecision] * N[(x / y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(0.5 * N[(y + N[(z * N[(N[(x - z), $MachinePrecision] / y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]

Derivation

1. Split input into 2 regimes

2. if z < 9.99999999999999914e28

   1. Initial program 69.0%

      \[\frac{\left(x \cdot x + y \cdot y\right) - z \cdot z}{y \cdot 2} \]
   2. Step-by-step derivation

      1. remove-double-neg 69.0%

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

      2. distribute-lft-neg-out 69.0%

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

      3. distribute-frac-neg2 69.0%

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

      4. distribute-frac-neg 69.0%

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

      5. neg-mul-1 69.0%

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

      6. distribute-lft-neg-out 69.0%

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

      7. *-commutative 69.0%

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

      8. distribute-lft-neg-in 69.0%

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

      9. times-frac 69.0%

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

      10. metadata-eval 69.0%

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

      11. metadata-eval 69.0%

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

      12. associate--l+ 69.0%

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

      13. fma-define 71.0%

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

   3. Simplified 71.0%

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

   5. Taylor expanded in x around 0 83.1%

      \[\leadsto 0.5 \cdot \color{blue}{\left(\left(y + \frac{{x}^{2}}{y}\right) - \frac{{z}^{2}}{y}\right)} \]
   6. Step-by-step derivation

      1. associate--l+ 83.1%

         \[\leadsto 0.5 \cdot \color{blue}{\left(y + \left(\frac{{x}^{2}}{y} - \frac{{z}^{2}}{y}\right)\right)} \]

      2. div-sub 87.0%

         \[\leadsto 0.5 \cdot \left(y + \color{blue}{\frac{{x}^{2} - {z}^{2}}{y}}\right) \]

   7. Simplified 87.0%

      \[\leadsto 0.5 \cdot \color{blue}{\left(y + \frac{{x}^{2} - {z}^{2}}{y}\right)} \]
   8. Step-by-step derivation

      1. pow2 87.0%

         \[\leadsto 0.5 \cdot \left(y + \frac{\color{blue}{x \cdot x} - {z}^{2}}{y}\right) \]

      2. pow2 87.0%

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

      3. difference-of-squares 91.0%

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

   9. Applied egg-rr 91.0%

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

       1. associate-/l* 99.9%

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

       2. *-commutative 99.9%

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

   11. Applied egg-rr 99.9%

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

   12. Taylor expanded in x around inf 79.3%

       \[\leadsto 0.5 \cdot \left(y + \color{blue}{\frac{x}{y}} \cdot \left(x + z\right)\right) \]

   if 9.99999999999999914e28 < z

   1. Initial program 60.7%

      \[\frac{\left(x \cdot x + y \cdot y\right) - z \cdot z}{y \cdot 2} \]
   2. Step-by-step derivation

      1. remove-double-neg 60.7%

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

      2. distribute-lft-neg-out 60.7%

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

      3. distribute-frac-neg2 60.7%

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

      4. distribute-frac-neg 60.7%

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

      5. neg-mul-1 60.7%

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

      6. distribute-lft-neg-out 60.7%

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

      7. *-commutative 60.7%

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

      8. distribute-lft-neg-in 60.7%

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

      9. times-frac 60.7%

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

      10. metadata-eval 60.7%

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

      11. metadata-eval 60.7%

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

      12. associate--l+ 60.7%

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

      13. fma-define 66.7%

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

   3. Simplified 66.7%

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

   5. Taylor expanded in x around 0 60.9%

      \[\leadsto 0.5 \cdot \color{blue}{\left(\left(y + \frac{{x}^{2}}{y}\right) - \frac{{z}^{2}}{y}\right)} \]
   6. Step-by-step derivation

      1. associate--l+ 60.9%

         \[\leadsto 0.5 \cdot \color{blue}{\left(y + \left(\frac{{x}^{2}}{y} - \frac{{z}^{2}}{y}\right)\right)} \]

      2. div-sub 70.9%

         \[\leadsto 0.5 \cdot \left(y + \color{blue}{\frac{{x}^{2} - {z}^{2}}{y}}\right) \]

   7. Simplified 70.9%

      \[\leadsto 0.5 \cdot \color{blue}{\left(y + \frac{{x}^{2} - {z}^{2}}{y}\right)} \]
   8. Step-by-step derivation

      1. pow2 70.9%

         \[\leadsto 0.5 \cdot \left(y + \frac{\color{blue}{x \cdot x} - {z}^{2}}{y}\right) \]

      2. pow2 70.9%

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

      3. difference-of-squares 81.3%

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

   9. Applied egg-rr 81.3%

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

       1. associate-/l* 99.8%

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

       2. *-commutative 99.8%

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

   11. Applied egg-rr 99.8%

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

   12. Taylor expanded in x around 0 88.5%

       \[\leadsto 0.5 \cdot \left(y + \frac{x - z}{y} \cdot \color{blue}{z}\right) \]
3. Recombined 2 regimes into one program.

4. Final simplification 81.1%

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

Alternative 3: 76.1% accurate, 0.9× speedup

Math

\[\begin{array}{l} \\ \begin{array}{l} \mathbf{if}\;y \leq 2.2 \cdot 10^{+88}:\\ \;\;\;\;0.5 \cdot \left(\left(x + z\right) \cdot \frac{x - z}{y}\right)\\ \mathbf{else}:\\ \;\;\;\;0.5 \cdot \left(y + \left(x + z\right) \cdot \frac{x}{y}\right)\\ \end{array} \end{array} \]

FPCore

(FPCore (x y z)
 :precision binary64
 (if (<= y 2.2e+88)
   (* 0.5 (* (+ x z) (/ (- x z) y)))
   (* 0.5 (+ y (* (+ x z) (/ x y))))))

C

double code(double x, double y, double z) {
	double tmp;
	if (y <= 2.2e+88) {
		tmp = 0.5 * ((x + z) * ((x - z) / y));
	} else {
		tmp = 0.5 * (y + ((x + z) * (x / y)));
	}
	return tmp;
}

Fortran

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 (y <= 2.2d+88) then
        tmp = 0.5d0 * ((x + z) * ((x - z) / y))
    else
        tmp = 0.5d0 * (y + ((x + z) * (x / y)))
    end if
    code = tmp
end function

Java

public static double code(double x, double y, double z) {
	double tmp;
	if (y <= 2.2e+88) {
		tmp = 0.5 * ((x + z) * ((x - z) / y));
	} else {
		tmp = 0.5 * (y + ((x + z) * (x / y)));
	}
	return tmp;
}

Python

def code(x, y, z):
	tmp = 0
	if y <= 2.2e+88:
		tmp = 0.5 * ((x + z) * ((x - z) / y))
	else:
		tmp = 0.5 * (y + ((x + z) * (x / y)))
	return tmp

Julia

function code(x, y, z)
	tmp = 0.0
	if (y <= 2.2e+88)
		tmp = Float64(0.5 * Float64(Float64(x + z) * Float64(Float64(x - z) / y)));
	else
		tmp = Float64(0.5 * Float64(y + Float64(Float64(x + z) * Float64(x / y))));
	end
	return tmp
end

MATLAB

function tmp_2 = code(x, y, z)
	tmp = 0.0;
	if (y <= 2.2e+88)
		tmp = 0.5 * ((x + z) * ((x - z) / y));
	else
		tmp = 0.5 * (y + ((x + z) * (x / y)));
	end
	tmp_2 = tmp;
end

Wolfram

code[x_, y_, z_] := If[LessEqual[y, 2.2e+88], N[(0.5 * N[(N[(x + z), $MachinePrecision] * N[(N[(x - z), $MachinePrecision] / y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(0.5 * N[(y + N[(N[(x + z), $MachinePrecision] * N[(x / y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]

Derivation

1. Split input into 2 regimes

2. if y < 2.20000000000000009e88

   1. Initial program 78.0%

      \[\frac{\left(x \cdot x + y \cdot y\right) - z \cdot z}{y \cdot 2} \]
   2. Step-by-step derivation

      1. remove-double-neg 78.0%

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

      2. distribute-lft-neg-out 78.0%

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

      3. distribute-frac-neg2 78.0%

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

      4. distribute-frac-neg 78.0%

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

      5. neg-mul-1 78.0%

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

      6. distribute-lft-neg-out 78.0%

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

      7. *-commutative 78.0%

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

      8. distribute-lft-neg-in 78.0%

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

      9. times-frac 78.0%

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

      10. metadata-eval 78.0%

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

      11. metadata-eval 78.0%

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

      12. associate--l+ 78.0%

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

      13. fma-define 81.4%

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

   3. Simplified 81.4%

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

   5. Taylor expanded in x around 0 81.1%

      \[\leadsto 0.5 \cdot \color{blue}{\left(\left(y + \frac{{x}^{2}}{y}\right) - \frac{{z}^{2}}{y}\right)} \]
   6. Step-by-step derivation

      1. associate--l+ 81.1%

         \[\leadsto 0.5 \cdot \color{blue}{\left(y + \left(\frac{{x}^{2}}{y} - \frac{{z}^{2}}{y}\right)\right)} \]

      2. div-sub 87.4%

         \[\leadsto 0.5 \cdot \left(y + \color{blue}{\frac{{x}^{2} - {z}^{2}}{y}}\right) \]

   7. Simplified 87.4%

      \[\leadsto 0.5 \cdot \color{blue}{\left(y + \frac{{x}^{2} - {z}^{2}}{y}\right)} \]
   8. Step-by-step derivation

      1. pow2 87.4%

         \[\leadsto 0.5 \cdot \left(y + \frac{\color{blue}{x \cdot x} - {z}^{2}}{y}\right) \]

      2. pow2 87.4%

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

      3. difference-of-squares 93.4%

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

   9. Applied egg-rr 93.4%

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

   10. Taylor expanded in y around 0 70.3%

       \[\leadsto \color{blue}{0.5 \cdot \frac{\left(x + z\right) \cdot \left(x - z\right)}{y}} \]
   11. Step-by-step derivation

       1. associate-*r/ 75.2%

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

       2. +-commutative 75.2%

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

   12. Simplified 75.2%

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

   if 2.20000000000000009e88 < y

   1. Initial program 23.6%

      \[\frac{\left(x \cdot x + y \cdot y\right) - z \cdot z}{y \cdot 2} \]
   2. Step-by-step derivation

      1. remove-double-neg 23.6%

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

      2. distribute-lft-neg-out 23.6%

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

      3. distribute-frac-neg2 23.6%

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

      4. distribute-frac-neg 23.6%

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

      5. neg-mul-1 23.6%

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

      6. distribute-lft-neg-out 23.6%

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

      7. *-commutative 23.6%

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

      8. distribute-lft-neg-in 23.6%

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

      9. times-frac 23.6%

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

      10. metadata-eval 23.6%

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

      11. metadata-eval 23.6%

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

      12. associate--l+ 23.6%

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

      13. fma-define 23.6%

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

   3. Simplified 23.6%

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

   5. Taylor expanded in x around 0 69.0%

      \[\leadsto 0.5 \cdot \color{blue}{\left(\left(y + \frac{{x}^{2}}{y}\right) - \frac{{z}^{2}}{y}\right)} \]
   6. Step-by-step derivation

      1. associate--l+ 69.0%

         \[\leadsto 0.5 \cdot \color{blue}{\left(y + \left(\frac{{x}^{2}}{y} - \frac{{z}^{2}}{y}\right)\right)} \]

      2. div-sub 69.0%

         \[\leadsto 0.5 \cdot \left(y + \color{blue}{\frac{{x}^{2} - {z}^{2}}{y}}\right) \]

   7. Simplified 69.0%

      \[\leadsto 0.5 \cdot \color{blue}{\left(y + \frac{{x}^{2} - {z}^{2}}{y}\right)} \]
   8. Step-by-step derivation

      1. pow2 69.0%

         \[\leadsto 0.5 \cdot \left(y + \frac{\color{blue}{x \cdot x} - {z}^{2}}{y}\right) \]

      2. pow2 69.0%

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

      3. difference-of-squares 71.3%

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

   9. Applied egg-rr 71.3%

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

       1. associate-/l* 100.0%

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

       2. *-commutative 100.0%

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

   11. Applied egg-rr 100.0%

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

   12. Taylor expanded in x around inf 86.7%

       \[\leadsto 0.5 \cdot \left(y + \color{blue}{\frac{x}{y}} \cdot \left(x + z\right)\right) \]
3. Recombined 2 regimes into one program.

4. Final simplification 77.4%

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

Alternative 4: 73.8% accurate, 0.9× speedup

Math

\[\begin{array}{l} \\ \begin{array}{l} \mathbf{if}\;y \leq 2.7 \cdot 10^{+117}:\\ \;\;\;\;0.5 \cdot \left(\left(x + z\right) \cdot \frac{x - z}{y}\right)\\ \mathbf{else}:\\ \;\;\;\;0.5 \cdot y\\ \end{array} \end{array} \]

FPCore

(FPCore (x y z)
 :precision binary64
 (if (<= y 2.7e+117) (* 0.5 (* (+ x z) (/ (- x z) y))) (* 0.5 y)))

C

double code(double x, double y, double z) {
	double tmp;
	if (y <= 2.7e+117) {
		tmp = 0.5 * ((x + z) * ((x - z) / y));
	} else {
		tmp = 0.5 * y;
	}
	return tmp;
}

Fortran

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 (y <= 2.7d+117) then
        tmp = 0.5d0 * ((x + z) * ((x - z) / y))
    else
        tmp = 0.5d0 * y
    end if
    code = tmp
end function

Java

public static double code(double x, double y, double z) {
	double tmp;
	if (y <= 2.7e+117) {
		tmp = 0.5 * ((x + z) * ((x - z) / y));
	} else {
		tmp = 0.5 * y;
	}
	return tmp;
}

Python

def code(x, y, z):
	tmp = 0
	if y <= 2.7e+117:
		tmp = 0.5 * ((x + z) * ((x - z) / y))
	else:
		tmp = 0.5 * y
	return tmp

Julia

function code(x, y, z)
	tmp = 0.0
	if (y <= 2.7e+117)
		tmp = Float64(0.5 * Float64(Float64(x + z) * Float64(Float64(x - z) / y)));
	else
		tmp = Float64(0.5 * y);
	end
	return tmp
end

MATLAB

function tmp_2 = code(x, y, z)
	tmp = 0.0;
	if (y <= 2.7e+117)
		tmp = 0.5 * ((x + z) * ((x - z) / y));
	else
		tmp = 0.5 * y;
	end
	tmp_2 = tmp;
end

Wolfram

code[x_, y_, z_] := If[LessEqual[y, 2.7e+117], N[(0.5 * N[(N[(x + z), $MachinePrecision] * N[(N[(x - z), $MachinePrecision] / y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(0.5 * y), $MachinePrecision]]

Derivation

1. Split input into 2 regimes

2. if y < 2.7000000000000002e117

   1. Initial program 78.1%

      \[\frac{\left(x \cdot x + y \cdot y\right) - z \cdot z}{y \cdot 2} \]
   2. Step-by-step derivation

      1. remove-double-neg 78.1%

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

      2. distribute-lft-neg-out 78.1%

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

      3. distribute-frac-neg2 78.1%

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

      4. distribute-frac-neg 78.1%

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

      5. neg-mul-1 78.1%

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

      6. distribute-lft-neg-out 78.1%

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

      7. *-commutative 78.1%

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

      8. distribute-lft-neg-in 78.1%

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

      9. times-frac 78.1%

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

      10. metadata-eval 78.1%

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

      11. metadata-eval 78.1%

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

      12. associate--l+ 78.1%

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

      13. fma-define 81.4%

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

   3. Simplified 81.4%

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

   5. Taylor expanded in x around 0 81.1%

      \[\leadsto 0.5 \cdot \color{blue}{\left(\left(y + \frac{{x}^{2}}{y}\right) - \frac{{z}^{2}}{y}\right)} \]
   6. Step-by-step derivation

      1. associate--l+ 81.1%

         \[\leadsto 0.5 \cdot \color{blue}{\left(y + \left(\frac{{x}^{2}}{y} - \frac{{z}^{2}}{y}\right)\right)} \]

      2. div-sub 87.3%

         \[\leadsto 0.5 \cdot \left(y + \color{blue}{\frac{{x}^{2} - {z}^{2}}{y}}\right) \]

   7. Simplified 87.3%

      \[\leadsto 0.5 \cdot \color{blue}{\left(y + \frac{{x}^{2} - {z}^{2}}{y}\right)} \]
   8. Step-by-step derivation

      1. pow2 87.3%

         \[\leadsto 0.5 \cdot \left(y + \frac{\color{blue}{x \cdot x} - {z}^{2}}{y}\right) \]

      2. pow2 87.3%

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

      3. difference-of-squares 93.1%

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

   9. Applied egg-rr 93.1%

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

   10. Taylor expanded in y around 0 69.6%

       \[\leadsto \color{blue}{0.5 \cdot \frac{\left(x + z\right) \cdot \left(x - z\right)}{y}} \]
   11. Step-by-step derivation

       1. associate-*r/ 74.8%

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

       2. +-commutative 74.8%

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

   12. Simplified 74.8%

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

   if 2.7000000000000002e117 < y

   1. Initial program 17.2%

      \[\frac{\left(x \cdot x + y \cdot y\right) - z \cdot z}{y \cdot 2} \]
   2. Step-by-step derivation

      1. remove-double-neg 17.2%

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

      2. distribute-lft-neg-out 17.2%

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

      3. distribute-frac-neg2 17.2%

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

      4. distribute-frac-neg 17.2%

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

      5. neg-mul-1 17.2%

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

      6. distribute-lft-neg-out 17.2%

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

      7. *-commutative 17.2%

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

      8. distribute-lft-neg-in 17.2%

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

      9. times-frac 17.2%

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

      10. metadata-eval 17.2%

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

      11. metadata-eval 17.2%

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

      12. associate--l+ 17.2%

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

      13. fma-define 17.2%

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

   3. Simplified 17.2%

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

   5. Taylor expanded in y around inf 76.9%

      \[\leadsto \color{blue}{0.5 \cdot y} \]
3. Recombined 2 regimes into one program.

4. Final simplification 75.2%

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

Alternative 5: 49.2% accurate, 1.1× speedup

Math

\[\begin{array}{l} \\ \begin{array}{l} \mathbf{if}\;y \leq 9.5 \cdot 10^{+90}:\\ \;\;\;\;0.5 \cdot \left(z \cdot \frac{x - z}{y}\right)\\ \mathbf{else}:\\ \;\;\;\;0.5 \cdot y\\ \end{array} \end{array} \]

FPCore

(FPCore (x y z)
 :precision binary64
 (if (<= y 9.5e+90) (* 0.5 (* z (/ (- x z) y))) (* 0.5 y)))

C

double code(double x, double y, double z) {
	double tmp;
	if (y <= 9.5e+90) {
		tmp = 0.5 * (z * ((x - z) / y));
	} else {
		tmp = 0.5 * y;
	}
	return tmp;
}

Fortran

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 (y <= 9.5d+90) then
        tmp = 0.5d0 * (z * ((x - z) / y))
    else
        tmp = 0.5d0 * y
    end if
    code = tmp
end function

Java

public static double code(double x, double y, double z) {
	double tmp;
	if (y <= 9.5e+90) {
		tmp = 0.5 * (z * ((x - z) / y));
	} else {
		tmp = 0.5 * y;
	}
	return tmp;
}

Python

def code(x, y, z):
	tmp = 0
	if y <= 9.5e+90:
		tmp = 0.5 * (z * ((x - z) / y))
	else:
		tmp = 0.5 * y
	return tmp

Julia

function code(x, y, z)
	tmp = 0.0
	if (y <= 9.5e+90)
		tmp = Float64(0.5 * Float64(z * Float64(Float64(x - z) / y)));
	else
		tmp = Float64(0.5 * y);
	end
	return tmp
end

MATLAB

function tmp_2 = code(x, y, z)
	tmp = 0.0;
	if (y <= 9.5e+90)
		tmp = 0.5 * (z * ((x - z) / y));
	else
		tmp = 0.5 * y;
	end
	tmp_2 = tmp;
end

Wolfram

code[x_, y_, z_] := If[LessEqual[y, 9.5e+90], N[(0.5 * N[(z * N[(N[(x - z), $MachinePrecision] / y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(0.5 * y), $MachinePrecision]]

Derivation

1. Split input into 2 regimes

2. if y < 9.4999999999999994e90

   1. Initial program 78.0%

      \[\frac{\left(x \cdot x + y \cdot y\right) - z \cdot z}{y \cdot 2} \]
   2. Step-by-step derivation

      1. remove-double-neg 78.0%

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

      2. distribute-lft-neg-out 78.0%

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

      3. distribute-frac-neg2 78.0%

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

      4. distribute-frac-neg 78.0%

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

      5. neg-mul-1 78.0%

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

      6. distribute-lft-neg-out 78.0%

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

      7. *-commutative 78.0%

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

      8. distribute-lft-neg-in 78.0%

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

      9. times-frac 78.0%

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

      10. metadata-eval 78.0%

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

      11. metadata-eval 78.0%

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

      12. associate--l+ 78.0%

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

      13. fma-define 81.4%

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

   3. Simplified 81.4%

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

   5. Taylor expanded in x around 0 81.1%

      \[\leadsto 0.5 \cdot \color{blue}{\left(\left(y + \frac{{x}^{2}}{y}\right) - \frac{{z}^{2}}{y}\right)} \]
   6. Step-by-step derivation

      1. associate--l+ 81.1%

         \[\leadsto 0.5 \cdot \color{blue}{\left(y + \left(\frac{{x}^{2}}{y} - \frac{{z}^{2}}{y}\right)\right)} \]

      2. div-sub 87.4%

         \[\leadsto 0.5 \cdot \left(y + \color{blue}{\frac{{x}^{2} - {z}^{2}}{y}}\right) \]

   7. Simplified 87.4%

      \[\leadsto 0.5 \cdot \color{blue}{\left(y + \frac{{x}^{2} - {z}^{2}}{y}\right)} \]
   8. Step-by-step derivation

      1. pow2 87.4%

         \[\leadsto 0.5 \cdot \left(y + \frac{\color{blue}{x \cdot x} - {z}^{2}}{y}\right) \]

      2. pow2 87.4%

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

      3. difference-of-squares 93.4%

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

   9. Applied egg-rr 93.4%

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

   10. Taylor expanded in y around 0 70.3%

       \[\leadsto \color{blue}{0.5 \cdot \frac{\left(x + z\right) \cdot \left(x - z\right)}{y}} \]
   11. Step-by-step derivation

       1. associate-*r/ 75.2%

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

       2. +-commutative 75.2%

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

   12. Simplified 75.2%

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

   13. Taylor expanded in z around inf 44.9%

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

   if 9.4999999999999994e90 < y

   1. Initial program 23.6%

      \[\frac{\left(x \cdot x + y \cdot y\right) - z \cdot z}{y \cdot 2} \]
   2. Step-by-step derivation

      1. remove-double-neg 23.6%

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

      2. distribute-lft-neg-out 23.6%

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

      3. distribute-frac-neg2 23.6%

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

      4. distribute-frac-neg 23.6%

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

      5. neg-mul-1 23.6%

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

      6. distribute-lft-neg-out 23.6%

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

      7. *-commutative 23.6%

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

      8. distribute-lft-neg-in 23.6%

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

      9. times-frac 23.6%

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

      10. metadata-eval 23.6%

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

      11. metadata-eval 23.6%

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

      12. associate--l+ 23.6%

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

      13. fma-define 23.6%

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

   3. Simplified 23.6%

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

   5. Taylor expanded in y around inf 73.5%

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

Alternative 6: 99.9% accurate, 1.2× speedup

Math

\[\begin{array}{l} \\ 0.5 \cdot \left(y + \left(x + z\right) \cdot \frac{x - z}{y}\right) \end{array} \]

FPCore

(FPCore (x y z) :precision binary64 (* 0.5 (+ y (* (+ x z) (/ (- x z) y)))))

C

double code(double x, double y, double z) {
	return 0.5 * (y + ((x + z) * ((x - z) / y)));
}

Fortran

real(8) function code(x, y, z)
    real(8), intent (in) :: x
    real(8), intent (in) :: y
    real(8), intent (in) :: z
    code = 0.5d0 * (y + ((x + z) * ((x - z) / y)))
end function

Java

public static double code(double x, double y, double z) {
	return 0.5 * (y + ((x + z) * ((x - z) / y)));
}

Python

def code(x, y, z):
	return 0.5 * (y + ((x + z) * ((x - z) / y)))

Julia

function code(x, y, z)
	return Float64(0.5 * Float64(y + Float64(Float64(x + z) * Float64(Float64(x - z) / y))))
end

MATLAB

function tmp = code(x, y, z)
	tmp = 0.5 * (y + ((x + z) * ((x - z) / y)));
end

Wolfram

code[x_, y_, z_] := N[(0.5 * N[(y + N[(N[(x + z), $MachinePrecision] * N[(N[(x - z), $MachinePrecision] / y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]

Derivation

1. Initial program 67.4%

   \[\frac{\left(x \cdot x + y \cdot y\right) - z \cdot z}{y \cdot 2} \]
2. Step-by-step derivation

   1. remove-double-neg 67.4%

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

   2. distribute-lft-neg-out 67.4%

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

   3. distribute-frac-neg2 67.4%

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

   4. distribute-frac-neg 67.4%

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

   5. neg-mul-1 67.4%

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

   6. distribute-lft-neg-out 67.4%

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

   7. *-commutative 67.4%

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

   8. distribute-lft-neg-in 67.4%

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

   9. times-frac 67.4%

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

   10. metadata-eval 67.4%

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

   11. metadata-eval 67.4%

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

   12. associate--l+ 67.4%

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

   13. fma-define 70.1%

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

3. Simplified 70.1%

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

5. Taylor expanded in x around 0 78.8%

   \[\leadsto 0.5 \cdot \color{blue}{\left(\left(y + \frac{{x}^{2}}{y}\right) - \frac{{z}^{2}}{y}\right)} \]
6. Step-by-step derivation

   1. associate--l+ 78.8%

      \[\leadsto 0.5 \cdot \color{blue}{\left(y + \left(\frac{{x}^{2}}{y} - \frac{{z}^{2}}{y}\right)\right)} \]

   2. div-sub 83.8%

      \[\leadsto 0.5 \cdot \left(y + \color{blue}{\frac{{x}^{2} - {z}^{2}}{y}}\right) \]

7. Simplified 83.8%

   \[\leadsto 0.5 \cdot \color{blue}{\left(y + \frac{{x}^{2} - {z}^{2}}{y}\right)} \]
8. Step-by-step derivation

   1. pow2 83.8%

      \[\leadsto 0.5 \cdot \left(y + \frac{\color{blue}{x \cdot x} - {z}^{2}}{y}\right) \]

   2. pow2 83.8%

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

   3. difference-of-squares 89.1%

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

9. Applied egg-rr 89.1%

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

    1. associate-/l* 99.9%

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

    2. *-commutative 99.9%

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

11. Applied egg-rr 99.9%

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

12. Final simplification 99.9%

    \[\leadsto 0.5 \cdot \left(y + \left(x + z\right) \cdot \frac{x - z}{y}\right) \]
13. Add Preprocessing

Alternative 7: 43.9% accurate, 1.2× speedup

Math

\[\begin{array}{l} \\ \begin{array}{l} \mathbf{if}\;x \leq 1.8 \cdot 10^{+64}:\\ \;\;\;\;0.5 \cdot y\\ \mathbf{else}:\\ \;\;\;\;x \cdot \frac{x}{y \cdot 2}\\ \end{array} \end{array} \]

FPCore

(FPCore (x y z)
 :precision binary64
 (if (<= x 1.8e+64) (* 0.5 y) (* x (/ x (* y 2.0)))))

C

double code(double x, double y, double z) {
	double tmp;
	if (x <= 1.8e+64) {
		tmp = 0.5 * y;
	} else {
		tmp = x * (x / (y * 2.0));
	}
	return tmp;
}

Fortran

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.8d+64) then
        tmp = 0.5d0 * y
    else
        tmp = x * (x / (y * 2.0d0))
    end if
    code = tmp
end function

Java

public static double code(double x, double y, double z) {
	double tmp;
	if (x <= 1.8e+64) {
		tmp = 0.5 * y;
	} else {
		tmp = x * (x / (y * 2.0));
	}
	return tmp;
}

Python

def code(x, y, z):
	tmp = 0
	if x <= 1.8e+64:
		tmp = 0.5 * y
	else:
		tmp = x * (x / (y * 2.0))
	return tmp

Julia

function code(x, y, z)
	tmp = 0.0
	if (x <= 1.8e+64)
		tmp = Float64(0.5 * y);
	else
		tmp = Float64(x * Float64(x / Float64(y * 2.0)));
	end
	return tmp
end

MATLAB

function tmp_2 = code(x, y, z)
	tmp = 0.0;
	if (x <= 1.8e+64)
		tmp = 0.5 * y;
	else
		tmp = x * (x / (y * 2.0));
	end
	tmp_2 = tmp;
end

Wolfram

code[x_, y_, z_] := If[LessEqual[x, 1.8e+64], N[(0.5 * y), $MachinePrecision], N[(x * N[(x / N[(y * 2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]

Derivation

1. Split input into 2 regimes

2. if x < 1.80000000000000007e64

   1. Initial program 70.1%

      \[\frac{\left(x \cdot x + y \cdot y\right) - z \cdot z}{y \cdot 2} \]
   2. Step-by-step derivation

      1. remove-double-neg 70.1%

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

      2. distribute-lft-neg-out 70.1%

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

      3. distribute-frac-neg2 70.1%

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

      4. distribute-frac-neg 70.1%

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

      5. neg-mul-1 70.1%

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

      6. distribute-lft-neg-out 70.1%

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

      7. *-commutative 70.1%

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

      8. distribute-lft-neg-in 70.1%

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

      9. times-frac 70.1%

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

      10. metadata-eval 70.1%

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

      11. metadata-eval 70.1%

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

      12. associate--l+ 70.1%

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

      13. fma-define 73.0%

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

   3. Simplified 73.0%

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

   5. Taylor expanded in y around inf 38.4%

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

   if 1.80000000000000007e64 < x

   1. Initial program 55.8%

      \[\frac{\left(x \cdot x + y \cdot y\right) - z \cdot z}{y \cdot 2} \]
   2. Step-by-step derivation

      1. remove-double-neg 55.8%

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

      2. distribute-lft-neg-out 55.8%

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

      3. distribute-frac-neg2 55.8%

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

      4. distribute-frac-neg 55.8%

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

      5. neg-mul-1 55.8%

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

      6. distribute-lft-neg-out 55.8%

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

      7. *-commutative 55.8%

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

      8. distribute-lft-neg-in 55.8%

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

      9. times-frac 55.8%

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

      10. metadata-eval 55.8%

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

      11. metadata-eval 55.8%

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

      12. associate--l+ 55.8%

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

      13. fma-define 57.9%

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

   3. Simplified 57.9%

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

   5. Taylor expanded in x around inf 54.5%

      \[\leadsto \color{blue}{0.5 \cdot \frac{{x}^{2}}{y}} \]
   6. Step-by-step derivation

      1. *-commutative 54.5%

         \[\leadsto \color{blue}{\frac{{x}^{2}}{y} \cdot 0.5} \]

      2. associate-*l/ 54.5%

         \[\leadsto \color{blue}{\frac{{x}^{2} \cdot 0.5}{y}} \]

      3. associate-*r/ 54.4%

         \[\leadsto \color{blue}{{x}^{2} \cdot \frac{0.5}{y}} \]

   7. Simplified 54.4%

      \[\leadsto \color{blue}{{x}^{2} \cdot \frac{0.5}{y}} \]
   8. Step-by-step derivation

      1. associate-*r/ 54.5%

         \[\leadsto \color{blue}{\frac{{x}^{2} \cdot 0.5}{y}} \]

      2. clear-num 54.4%

         \[\leadsto \color{blue}{\frac{1}{\frac{y}{{x}^{2} \cdot 0.5}}} \]

      3. *-commutative 54.4%

         \[\leadsto \frac{1}{\frac{y}{\color{blue}{0.5 \cdot {x}^{2}}}} \]

   9. Applied egg-rr 54.4%

      \[\leadsto \color{blue}{\frac{1}{\frac{y}{0.5 \cdot {x}^{2}}}} \]
   10. Step-by-step derivation

       1. clear-num 54.5%

          \[\leadsto \color{blue}{\frac{0.5 \cdot {x}^{2}}{y}} \]

       2. associate-*l/ 54.4%

          \[\leadsto \color{blue}{\frac{0.5}{y} \cdot {x}^{2}} \]

       3. pow2 54.4%

          \[\leadsto \frac{0.5}{y} \cdot \color{blue}{\left(x \cdot x\right)} \]

       4. associate-*r* 69.0%

          \[\leadsto \color{blue}{\left(\frac{0.5}{y} \cdot x\right) \cdot x} \]

       5. *-commutative 69.0%

          \[\leadsto \color{blue}{\left(x \cdot \frac{0.5}{y}\right)} \cdot x \]

       6. clear-num 69.0%

          \[\leadsto \left(x \cdot \color{blue}{\frac{1}{\frac{y}{0.5}}}\right) \cdot x \]

       7. un-div-inv 69.1%

          \[\leadsto \color{blue}{\frac{x}{\frac{y}{0.5}}} \cdot x \]

       8. div-inv 69.1%

          \[\leadsto \frac{x}{\color{blue}{y \cdot \frac{1}{0.5}}} \cdot x \]

       9. metadata-eval 69.1%

          \[\leadsto \frac{x}{y \cdot \color{blue}{2}} \cdot x \]

   11. Applied egg-rr 69.1%

       \[\leadsto \color{blue}{\frac{x}{y \cdot 2} \cdot x} \]
3. Recombined 2 regimes into one program.

4. Final simplification 44.1%

   \[\leadsto \begin{array}{l} \mathbf{if}\;x \leq 1.8 \cdot 10^{+64}:\\ \;\;\;\;0.5 \cdot y\\ \mathbf{else}:\\ \;\;\;\;x \cdot \frac{x}{y \cdot 2}\\ \end{array} \]
5. Add Preprocessing

Alternative 8: 41.5% accurate, 1.2× speedup

Math

\[\begin{array}{l} \\ \begin{array}{l} \mathbf{if}\;y \leq 4.8 \cdot 10^{+112}:\\ \;\;\;\;\left(x \cdot x\right) \cdot \frac{0.5}{y}\\ \mathbf{else}:\\ \;\;\;\;0.5 \cdot y\\ \end{array} \end{array} \]

FPCore

(FPCore (x y z)
 :precision binary64
 (if (<= y 4.8e+112) (* (* x x) (/ 0.5 y)) (* 0.5 y)))

C

double code(double x, double y, double z) {
	double tmp;
	if (y <= 4.8e+112) {
		tmp = (x * x) * (0.5 / y);
	} else {
		tmp = 0.5 * y;
	}
	return tmp;
}

Fortran

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 (y <= 4.8d+112) then
        tmp = (x * x) * (0.5d0 / y)
    else
        tmp = 0.5d0 * y
    end if
    code = tmp
end function

Java

public static double code(double x, double y, double z) {
	double tmp;
	if (y <= 4.8e+112) {
		tmp = (x * x) * (0.5 / y);
	} else {
		tmp = 0.5 * y;
	}
	return tmp;
}

Python

def code(x, y, z):
	tmp = 0
	if y <= 4.8e+112:
		tmp = (x * x) * (0.5 / y)
	else:
		tmp = 0.5 * y
	return tmp

Julia

function code(x, y, z)
	tmp = 0.0
	if (y <= 4.8e+112)
		tmp = Float64(Float64(x * x) * Float64(0.5 / y));
	else
		tmp = Float64(0.5 * y);
	end
	return tmp
end

MATLAB

function tmp_2 = code(x, y, z)
	tmp = 0.0;
	if (y <= 4.8e+112)
		tmp = (x * x) * (0.5 / y);
	else
		tmp = 0.5 * y;
	end
	tmp_2 = tmp;
end

Wolfram

code[x_, y_, z_] := If[LessEqual[y, 4.8e+112], N[(N[(x * x), $MachinePrecision] * N[(0.5 / y), $MachinePrecision]), $MachinePrecision], N[(0.5 * y), $MachinePrecision]]

Derivation

1. Split input into 2 regimes

2. if y < 4.8e112

   1. Initial program 78.1%

      \[\frac{\left(x \cdot x + y \cdot y\right) - z \cdot z}{y \cdot 2} \]
   2. Step-by-step derivation

      1. remove-double-neg 78.1%

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

      2. distribute-lft-neg-out 78.1%

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

      3. distribute-frac-neg2 78.1%

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

      4. distribute-frac-neg 78.1%

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

      5. neg-mul-1 78.1%

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

      6. distribute-lft-neg-out 78.1%

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

      7. *-commutative 78.1%

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

      8. distribute-lft-neg-in 78.1%

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

      9. times-frac 78.1%

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

      10. metadata-eval 78.1%

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

      11. metadata-eval 78.1%

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

      12. associate--l+ 78.1%

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

      13. fma-define 81.4%

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

   3. Simplified 81.4%

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

   5. Taylor expanded in x around inf 36.7%

      \[\leadsto \color{blue}{0.5 \cdot \frac{{x}^{2}}{y}} \]
   6. Step-by-step derivation

      1. *-commutative 36.7%

         \[\leadsto \color{blue}{\frac{{x}^{2}}{y} \cdot 0.5} \]

      2. associate-*l/ 36.7%

         \[\leadsto \color{blue}{\frac{{x}^{2} \cdot 0.5}{y}} \]

      3. associate-*r/ 36.6%

         \[\leadsto \color{blue}{{x}^{2} \cdot \frac{0.5}{y}} \]

   7. Simplified 36.6%

      \[\leadsto \color{blue}{{x}^{2} \cdot \frac{0.5}{y}} \]
   8. Step-by-step derivation

      1. pow2 36.6%

         \[\leadsto \color{blue}{\left(x \cdot x\right)} \cdot \frac{0.5}{y} \]

   9. Applied egg-rr 36.6%

      \[\leadsto \color{blue}{\left(x \cdot x\right)} \cdot \frac{0.5}{y} \]

   if 4.8e112 < y

   1. Initial program 17.2%

      \[\frac{\left(x \cdot x + y \cdot y\right) - z \cdot z}{y \cdot 2} \]
   2. Step-by-step derivation

      1. remove-double-neg 17.2%

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

      2. distribute-lft-neg-out 17.2%

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

      3. distribute-frac-neg2 17.2%

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

      4. distribute-frac-neg 17.2%

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

      5. neg-mul-1 17.2%

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

      6. distribute-lft-neg-out 17.2%

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

      7. *-commutative 17.2%

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

      8. distribute-lft-neg-in 17.2%

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

      9. times-frac 17.2%

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

      10. metadata-eval 17.2%

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

      11. metadata-eval 17.2%

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

      12. associate--l+ 17.2%

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

      13. fma-define 17.2%

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

   3. Simplified 17.2%

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

   5. Taylor expanded in y around inf 76.9%

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

Alternative 9: 34.4% accurate, 5.0× speedup

Math

\[\begin{array}{l} \\ 0.5 \cdot y \end{array} \]

FPCore

(FPCore (x y z) :precision binary64 (* 0.5 y))

C

double code(double x, double y, double z) {
	return 0.5 * y;
}

Fortran

real(8) function code(x, y, z)
    real(8), intent (in) :: x
    real(8), intent (in) :: y
    real(8), intent (in) :: z
    code = 0.5d0 * y
end function

Java

public static double code(double x, double y, double z) {
	return 0.5 * y;
}

Python

def code(x, y, z):
	return 0.5 * y

Julia

function code(x, y, z)
	return Float64(0.5 * y)
end

MATLAB

function tmp = code(x, y, z)
	tmp = 0.5 * y;
end

Wolfram

code[x_, y_, z_] := N[(0.5 * y), $MachinePrecision]

Derivation

1. Initial program 67.4%

   \[\frac{\left(x \cdot x + y \cdot y\right) - z \cdot z}{y \cdot 2} \]
2. Step-by-step derivation

   1. remove-double-neg 67.4%

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

   2. distribute-lft-neg-out 67.4%

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

   3. distribute-frac-neg2 67.4%

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

   4. distribute-frac-neg 67.4%

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

   5. neg-mul-1 67.4%

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

   6. distribute-lft-neg-out 67.4%

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

   7. *-commutative 67.4%

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

   8. distribute-lft-neg-in 67.4%

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

   9. times-frac 67.4%

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

   10. metadata-eval 67.4%

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

   11. metadata-eval 67.4%

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

   12. associate--l+ 67.4%

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

   13. fma-define 70.1%

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

3. Simplified 70.1%

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

5. Taylor expanded in y around inf 35.6%

   \[\leadsto \color{blue}{0.5 \cdot y} \]
6. Add Preprocessing

Developer Target 1: 99.8% accurate, 1.0× speedup

Math

\[\begin{array}{l} \\ y \cdot 0.5 - \left(\frac{0.5}{y} \cdot \left(z + x\right)\right) \cdot \left(z - x\right) \end{array} \]

FPCore

(FPCore (x y z)
 :precision binary64
 (- (* y 0.5) (* (* (/ 0.5 y) (+ z x)) (- z x))))

C

double code(double x, double y, double z) {
	return (y * 0.5) - (((0.5 / y) * (z + x)) * (z - x));
}

Fortran

real(8) function code(x, y, z)
    real(8), intent (in) :: x
    real(8), intent (in) :: y
    real(8), intent (in) :: z
    code = (y * 0.5d0) - (((0.5d0 / y) * (z + x)) * (z - x))
end function

Java

public static double code(double x, double y, double z) {
	return (y * 0.5) - (((0.5 / y) * (z + x)) * (z - x));
}

Python

def code(x, y, z):
	return (y * 0.5) - (((0.5 / y) * (z + x)) * (z - x))

Julia

function code(x, y, z)
	return Float64(Float64(y * 0.5) - Float64(Float64(Float64(0.5 / y) * Float64(z + x)) * Float64(z - x)))
end

MATLAB

function tmp = code(x, y, z)
	tmp = (y * 0.5) - (((0.5 / y) * (z + x)) * (z - x));
end

Wolfram

code[x_, y_, z_] := N[(N[(y * 0.5), $MachinePrecision] - N[(N[(N[(0.5 / y), $MachinePrecision] * N[(z + x), $MachinePrecision]), $MachinePrecision] * N[(z - x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]

Reproduce

herbie shell --seed 2024165 
(FPCore (x y z)
  :name "Diagrams.TwoD.Apollonian:initialConfig from diagrams-contrib-1.3.0.5, A"
  :precision binary64

  :alt
  (! :herbie-platform default (- (* y 1/2) (* (* (/ 1/2 y) (+ z x)) (- z x))))

  (/ (- (+ (* x x) (* y y)) (* z z)) (* y 2.0)))