Numeric.SpecFunctions:logGamma from math-functions-0.1.5.2, C

Percentage Accurate: 58.7% → 98.6%

Time: 15.8s

Alternatives: 14

Speedup: 4.4×

Specification

Math

\[\begin{array}{l} \\ \frac{\left(x - 2\right) \cdot \left(\left(\left(\left(x \cdot 4.16438922228 + 78.6994924154\right) \cdot x + 137.519416416\right) \cdot x + y\right) \cdot x + z\right)}{\left(\left(\left(x + 43.3400022514\right) \cdot x + 263.505074721\right) \cdot x + 313.399215894\right) \cdot x + 47.066876606} \end{array} \]

FPCore

(FPCore (x y z)
 :precision binary64
 (/
  (*
   (- x 2.0)
   (+
    (*
     (+ (* (+ (* (+ (* x 4.16438922228) 78.6994924154) x) 137.519416416) x) y)
     x)
    z))
  (+
   (* (+ (* (+ (* (+ x 43.3400022514) x) 263.505074721) x) 313.399215894) x)
   47.066876606)))

C

double code(double x, double y, double z) {
	return ((x - 2.0) * ((((((((x * 4.16438922228) + 78.6994924154) * x) + 137.519416416) * x) + y) * x) + z)) / (((((((x + 43.3400022514) * x) + 263.505074721) * x) + 313.399215894) * x) + 47.066876606);
}

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 - 2.0d0) * ((((((((x * 4.16438922228d0) + 78.6994924154d0) * x) + 137.519416416d0) * x) + y) * x) + z)) / (((((((x + 43.3400022514d0) * x) + 263.505074721d0) * x) + 313.399215894d0) * x) + 47.066876606d0)
end function

Java

public static double code(double x, double y, double z) {
	return ((x - 2.0) * ((((((((x * 4.16438922228) + 78.6994924154) * x) + 137.519416416) * x) + y) * x) + z)) / (((((((x + 43.3400022514) * x) + 263.505074721) * x) + 313.399215894) * x) + 47.066876606);
}

Python

def code(x, y, z):
	return ((x - 2.0) * ((((((((x * 4.16438922228) + 78.6994924154) * x) + 137.519416416) * x) + y) * x) + z)) / (((((((x + 43.3400022514) * x) + 263.505074721) * x) + 313.399215894) * x) + 47.066876606)

Julia

function code(x, y, z)
	return Float64(Float64(Float64(x - 2.0) * Float64(Float64(Float64(Float64(Float64(Float64(Float64(Float64(x * 4.16438922228) + 78.6994924154) * x) + 137.519416416) * x) + y) * x) + z)) / Float64(Float64(Float64(Float64(Float64(Float64(Float64(x + 43.3400022514) * x) + 263.505074721) * x) + 313.399215894) * x) + 47.066876606))
end

MATLAB

function tmp = code(x, y, z)
	tmp = ((x - 2.0) * ((((((((x * 4.16438922228) + 78.6994924154) * x) + 137.519416416) * x) + y) * x) + z)) / (((((((x + 43.3400022514) * x) + 263.505074721) * x) + 313.399215894) * x) + 47.066876606);
end

Wolfram

code[x_, y_, z_] := N[(N[(N[(x - 2.0), $MachinePrecision] * N[(N[(N[(N[(N[(N[(N[(N[(x * 4.16438922228), $MachinePrecision] + 78.6994924154), $MachinePrecision] * x), $MachinePrecision] + 137.519416416), $MachinePrecision] * x), $MachinePrecision] + y), $MachinePrecision] * x), $MachinePrecision] + z), $MachinePrecision]), $MachinePrecision] / N[(N[(N[(N[(N[(N[(N[(x + 43.3400022514), $MachinePrecision] * x), $MachinePrecision] + 263.505074721), $MachinePrecision] * x), $MachinePrecision] + 313.399215894), $MachinePrecision] * x), $MachinePrecision] + 47.066876606), $MachinePrecision]), $MachinePrecision]

Initial Program: 58.7% accurate, 1.0× speedup

Math

\[\begin{array}{l} \\ \frac{\left(x - 2\right) \cdot \left(\left(\left(\left(x \cdot 4.16438922228 + 78.6994924154\right) \cdot x + 137.519416416\right) \cdot x + y\right) \cdot x + z\right)}{\left(\left(\left(x + 43.3400022514\right) \cdot x + 263.505074721\right) \cdot x + 313.399215894\right) \cdot x + 47.066876606} \end{array} \]

FPCore

(FPCore (x y z)
 :precision binary64
 (/
  (*
   (- x 2.0)
   (+
    (*
     (+ (* (+ (* (+ (* x 4.16438922228) 78.6994924154) x) 137.519416416) x) y)
     x)
    z))
  (+
   (* (+ (* (+ (* (+ x 43.3400022514) x) 263.505074721) x) 313.399215894) x)
   47.066876606)))

C

double code(double x, double y, double z) {
	return ((x - 2.0) * ((((((((x * 4.16438922228) + 78.6994924154) * x) + 137.519416416) * x) + y) * x) + z)) / (((((((x + 43.3400022514) * x) + 263.505074721) * x) + 313.399215894) * x) + 47.066876606);
}

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 - 2.0d0) * ((((((((x * 4.16438922228d0) + 78.6994924154d0) * x) + 137.519416416d0) * x) + y) * x) + z)) / (((((((x + 43.3400022514d0) * x) + 263.505074721d0) * x) + 313.399215894d0) * x) + 47.066876606d0)
end function

Java

public static double code(double x, double y, double z) {
	return ((x - 2.0) * ((((((((x * 4.16438922228) + 78.6994924154) * x) + 137.519416416) * x) + y) * x) + z)) / (((((((x + 43.3400022514) * x) + 263.505074721) * x) + 313.399215894) * x) + 47.066876606);
}

Python

def code(x, y, z):
	return ((x - 2.0) * ((((((((x * 4.16438922228) + 78.6994924154) * x) + 137.519416416) * x) + y) * x) + z)) / (((((((x + 43.3400022514) * x) + 263.505074721) * x) + 313.399215894) * x) + 47.066876606)

Julia

function code(x, y, z)
	return Float64(Float64(Float64(x - 2.0) * Float64(Float64(Float64(Float64(Float64(Float64(Float64(Float64(x * 4.16438922228) + 78.6994924154) * x) + 137.519416416) * x) + y) * x) + z)) / Float64(Float64(Float64(Float64(Float64(Float64(Float64(x + 43.3400022514) * x) + 263.505074721) * x) + 313.399215894) * x) + 47.066876606))
end

MATLAB

function tmp = code(x, y, z)
	tmp = ((x - 2.0) * ((((((((x * 4.16438922228) + 78.6994924154) * x) + 137.519416416) * x) + y) * x) + z)) / (((((((x + 43.3400022514) * x) + 263.505074721) * x) + 313.399215894) * x) + 47.066876606);
end

Wolfram

code[x_, y_, z_] := N[(N[(N[(x - 2.0), $MachinePrecision] * N[(N[(N[(N[(N[(N[(N[(N[(x * 4.16438922228), $MachinePrecision] + 78.6994924154), $MachinePrecision] * x), $MachinePrecision] + 137.519416416), $MachinePrecision] * x), $MachinePrecision] + y), $MachinePrecision] * x), $MachinePrecision] + z), $MachinePrecision]), $MachinePrecision] / N[(N[(N[(N[(N[(N[(N[(x + 43.3400022514), $MachinePrecision] * x), $MachinePrecision] + 263.505074721), $MachinePrecision] * x), $MachinePrecision] + 313.399215894), $MachinePrecision] * x), $MachinePrecision] + 47.066876606), $MachinePrecision]), $MachinePrecision]

Alternative 1: 98.6% accurate, 0.5× speedup

Math

\[\begin{array}{l} \\ \begin{array}{l} \mathbf{if}\;\frac{\left(x - 2\right) \cdot \left(\left(\left(\left(x \cdot 4.16438922228 + 78.6994924154\right) \cdot x + 137.519416416\right) \cdot x + y\right) \cdot x + z\right)}{\left(\left(\left(x + 43.3400022514\right) \cdot x + 263.505074721\right) \cdot x + 313.399215894\right) \cdot x + 47.066876606} \leq \infty:\\ \;\;\;\;\frac{\mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, 4.16438922228, 78.6994924154\right), 137.519416416\right), y\right), z\right)}{\mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, x + 43.3400022514, 263.505074721\right), 313.399215894\right), 47.066876606\right)} \cdot \left(x + -2\right)\\ \mathbf{else}:\\ \;\;\;\;\left(4.16438922228 - \frac{101.7851458539211 - \frac{3451.550173699799 + \frac{y + -124074.40615218398}{x}}{x}}{x}\right) \cdot \left(x + -2\right)\\ \end{array} \end{array} \]

FPCore

(FPCore (x y z)
 :precision binary64
 (if (<=
      (/
       (*
        (- x 2.0)
        (+
         (*
          (+
           (* (+ (* (+ (* x 4.16438922228) 78.6994924154) x) 137.519416416) x)
           y)
          x)
         z))
       (+
        (*
         (+ (* (+ (* (+ x 43.3400022514) x) 263.505074721) x) 313.399215894)
         x)
        47.066876606))
      INFINITY)
   (*
    (/
     (fma
      x
      (fma x (fma x (fma x 4.16438922228 78.6994924154) 137.519416416) y)
      z)
     (fma
      x
      (fma x (fma x (+ x 43.3400022514) 263.505074721) 313.399215894)
      47.066876606))
    (+ x -2.0))
   (*
    (-
     4.16438922228
     (/
      (-
       101.7851458539211
       (/ (+ 3451.550173699799 (/ (+ y -124074.40615218398) x)) x))
      x))
    (+ x -2.0))))

C

double code(double x, double y, double z) {
	double tmp;
	if ((((x - 2.0) * ((((((((x * 4.16438922228) + 78.6994924154) * x) + 137.519416416) * x) + y) * x) + z)) / (((((((x + 43.3400022514) * x) + 263.505074721) * x) + 313.399215894) * x) + 47.066876606)) <= ((double) INFINITY)) {
		tmp = (fma(x, fma(x, fma(x, fma(x, 4.16438922228, 78.6994924154), 137.519416416), y), z) / fma(x, fma(x, fma(x, (x + 43.3400022514), 263.505074721), 313.399215894), 47.066876606)) * (x + -2.0);
	} else {
		tmp = (4.16438922228 - ((101.7851458539211 - ((3451.550173699799 + ((y + -124074.40615218398) / x)) / x)) / x)) * (x + -2.0);
	}
	return tmp;
}

Julia

function code(x, y, z)
	tmp = 0.0
	if (Float64(Float64(Float64(x - 2.0) * Float64(Float64(Float64(Float64(Float64(Float64(Float64(Float64(x * 4.16438922228) + 78.6994924154) * x) + 137.519416416) * x) + y) * x) + z)) / Float64(Float64(Float64(Float64(Float64(Float64(Float64(x + 43.3400022514) * x) + 263.505074721) * x) + 313.399215894) * x) + 47.066876606)) <= Inf)
		tmp = Float64(Float64(fma(x, fma(x, fma(x, fma(x, 4.16438922228, 78.6994924154), 137.519416416), y), z) / fma(x, fma(x, fma(x, Float64(x + 43.3400022514), 263.505074721), 313.399215894), 47.066876606)) * Float64(x + -2.0));
	else
		tmp = Float64(Float64(4.16438922228 - Float64(Float64(101.7851458539211 - Float64(Float64(3451.550173699799 + Float64(Float64(y + -124074.40615218398) / x)) / x)) / x)) * Float64(x + -2.0));
	end
	return tmp
end

Wolfram

code[x_, y_, z_] := If[LessEqual[N[(N[(N[(x - 2.0), $MachinePrecision] * N[(N[(N[(N[(N[(N[(N[(N[(x * 4.16438922228), $MachinePrecision] + 78.6994924154), $MachinePrecision] * x), $MachinePrecision] + 137.519416416), $MachinePrecision] * x), $MachinePrecision] + y), $MachinePrecision] * x), $MachinePrecision] + z), $MachinePrecision]), $MachinePrecision] / N[(N[(N[(N[(N[(N[(N[(x + 43.3400022514), $MachinePrecision] * x), $MachinePrecision] + 263.505074721), $MachinePrecision] * x), $MachinePrecision] + 313.399215894), $MachinePrecision] * x), $MachinePrecision] + 47.066876606), $MachinePrecision]), $MachinePrecision], Infinity], N[(N[(N[(x * N[(x * N[(x * N[(x * 4.16438922228 + 78.6994924154), $MachinePrecision] + 137.519416416), $MachinePrecision] + y), $MachinePrecision] + z), $MachinePrecision] / N[(x * N[(x * N[(x * N[(x + 43.3400022514), $MachinePrecision] + 263.505074721), $MachinePrecision] + 313.399215894), $MachinePrecision] + 47.066876606), $MachinePrecision]), $MachinePrecision] * N[(x + -2.0), $MachinePrecision]), $MachinePrecision], N[(N[(4.16438922228 - N[(N[(101.7851458539211 - N[(N[(3451.550173699799 + N[(N[(y + -124074.40615218398), $MachinePrecision] / x), $MachinePrecision]), $MachinePrecision] / x), $MachinePrecision]), $MachinePrecision] / x), $MachinePrecision]), $MachinePrecision] * N[(x + -2.0), $MachinePrecision]), $MachinePrecision]]

Derivation

1. Split input into 2 regimes

2. if (/.f64 (*.f64 (-.f64 x #s(literal 2 binary64)) (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 x #s(literal 104109730557/25000000000 binary64)) #s(literal 393497462077/5000000000 binary64)) x) #s(literal 4297481763/31250000 binary64)) x) y) x) z)) (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 x #s(literal 216700011257/5000000000 binary64)) x) #s(literal 263505074721/1000000000 binary64)) x) #s(literal 156699607947/500000000 binary64)) x) #s(literal 23533438303/500000000 binary64))) < +inf.0

   1. Initial program 93.4%

      \[\frac{\left(x - 2\right) \cdot \left(\left(\left(\left(x \cdot 4.16438922228 + 78.6994924154\right) \cdot x + 137.519416416\right) \cdot x + y\right) \cdot x + z\right)}{\left(\left(\left(x + 43.3400022514\right) \cdot x + 263.505074721\right) \cdot x + 313.399215894\right) \cdot x + 47.066876606} \]
   2. Add Preprocessing
   3. Step-by-step derivation

      1. lift-/.f64 N/A

         \[\leadsto \color{blue}{\frac{\left(x - 2\right) \cdot \left(\left(\left(\left(x \cdot \frac{104109730557}{25000000000} + \frac{393497462077}{5000000000}\right) \cdot x + \frac{4297481763}{31250000}\right) \cdot x + y\right) \cdot x + z\right)}{\left(\left(\left(x + \frac{216700011257}{5000000000}\right) \cdot x + \frac{263505074721}{1000000000}\right) \cdot x + \frac{156699607947}{500000000}\right) \cdot x + \frac{23533438303}{500000000}}} \]

      2. lift-*.f64 N/A

         \[\leadsto \frac{\color{blue}{\left(x - 2\right) \cdot \left(\left(\left(\left(x \cdot \frac{104109730557}{25000000000} + \frac{393497462077}{5000000000}\right) \cdot x + \frac{4297481763}{31250000}\right) \cdot x + y\right) \cdot x + z\right)}}{\left(\left(\left(x + \frac{216700011257}{5000000000}\right) \cdot x + \frac{263505074721}{1000000000}\right) \cdot x + \frac{156699607947}{500000000}\right) \cdot x + \frac{23533438303}{500000000}} \]

      3. associate-/l* N/A

         \[\leadsto \color{blue}{\left(x - 2\right) \cdot \frac{\left(\left(\left(x \cdot \frac{104109730557}{25000000000} + \frac{393497462077}{5000000000}\right) \cdot x + \frac{4297481763}{31250000}\right) \cdot x + y\right) \cdot x + z}{\left(\left(\left(x + \frac{216700011257}{5000000000}\right) \cdot x + \frac{263505074721}{1000000000}\right) \cdot x + \frac{156699607947}{500000000}\right) \cdot x + \frac{23533438303}{500000000}}} \]

      4. *-commutative N/A

         \[\leadsto \color{blue}{\frac{\left(\left(\left(x \cdot \frac{104109730557}{25000000000} + \frac{393497462077}{5000000000}\right) \cdot x + \frac{4297481763}{31250000}\right) \cdot x + y\right) \cdot x + z}{\left(\left(\left(x + \frac{216700011257}{5000000000}\right) \cdot x + \frac{263505074721}{1000000000}\right) \cdot x + \frac{156699607947}{500000000}\right) \cdot x + \frac{23533438303}{500000000}} \cdot \left(x - 2\right)} \]

      5. lower-*.f64 N/A

         \[\leadsto \color{blue}{\frac{\left(\left(\left(x \cdot \frac{104109730557}{25000000000} + \frac{393497462077}{5000000000}\right) \cdot x + \frac{4297481763}{31250000}\right) \cdot x + y\right) \cdot x + z}{\left(\left(\left(x + \frac{216700011257}{5000000000}\right) \cdot x + \frac{263505074721}{1000000000}\right) \cdot x + \frac{156699607947}{500000000}\right) \cdot x + \frac{23533438303}{500000000}} \cdot \left(x - 2\right)} \]

   4. Applied rewrites 99.1%

      \[\leadsto \color{blue}{\frac{\mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, 4.16438922228, 78.6994924154\right), 137.519416416\right), y\right), z\right)}{\mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, x + 43.3400022514, 263.505074721\right), 313.399215894\right), 47.066876606\right)} \cdot \left(x + -2\right)} \]

   if +inf.0 < (/.f64 (*.f64 (-.f64 x #s(literal 2 binary64)) (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 x #s(literal 104109730557/25000000000 binary64)) #s(literal 393497462077/5000000000 binary64)) x) #s(literal 4297481763/31250000 binary64)) x) y) x) z)) (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 x #s(literal 216700011257/5000000000 binary64)) x) #s(literal 263505074721/1000000000 binary64)) x) #s(literal 156699607947/500000000 binary64)) x) #s(literal 23533438303/500000000 binary64)))

   1. Initial program 0.0%

      \[\frac{\left(x - 2\right) \cdot \left(\left(\left(\left(x \cdot 4.16438922228 + 78.6994924154\right) \cdot x + 137.519416416\right) \cdot x + y\right) \cdot x + z\right)}{\left(\left(\left(x + 43.3400022514\right) \cdot x + 263.505074721\right) \cdot x + 313.399215894\right) \cdot x + 47.066876606} \]
   2. Add Preprocessing
   3. Step-by-step derivation

      1. lift-/.f64 N/A

         \[\leadsto \color{blue}{\frac{\left(x - 2\right) \cdot \left(\left(\left(\left(x \cdot \frac{104109730557}{25000000000} + \frac{393497462077}{5000000000}\right) \cdot x + \frac{4297481763}{31250000}\right) \cdot x + y\right) \cdot x + z\right)}{\left(\left(\left(x + \frac{216700011257}{5000000000}\right) \cdot x + \frac{263505074721}{1000000000}\right) \cdot x + \frac{156699607947}{500000000}\right) \cdot x + \frac{23533438303}{500000000}}} \]

      2. lift-*.f64 N/A

         \[\leadsto \frac{\color{blue}{\left(x - 2\right) \cdot \left(\left(\left(\left(x \cdot \frac{104109730557}{25000000000} + \frac{393497462077}{5000000000}\right) \cdot x + \frac{4297481763}{31250000}\right) \cdot x + y\right) \cdot x + z\right)}}{\left(\left(\left(x + \frac{216700011257}{5000000000}\right) \cdot x + \frac{263505074721}{1000000000}\right) \cdot x + \frac{156699607947}{500000000}\right) \cdot x + \frac{23533438303}{500000000}} \]

      3. associate-/l* N/A

         \[\leadsto \color{blue}{\left(x - 2\right) \cdot \frac{\left(\left(\left(x \cdot \frac{104109730557}{25000000000} + \frac{393497462077}{5000000000}\right) \cdot x + \frac{4297481763}{31250000}\right) \cdot x + y\right) \cdot x + z}{\left(\left(\left(x + \frac{216700011257}{5000000000}\right) \cdot x + \frac{263505074721}{1000000000}\right) \cdot x + \frac{156699607947}{500000000}\right) \cdot x + \frac{23533438303}{500000000}}} \]

      4. *-commutative N/A

         \[\leadsto \color{blue}{\frac{\left(\left(\left(x \cdot \frac{104109730557}{25000000000} + \frac{393497462077}{5000000000}\right) \cdot x + \frac{4297481763}{31250000}\right) \cdot x + y\right) \cdot x + z}{\left(\left(\left(x + \frac{216700011257}{5000000000}\right) \cdot x + \frac{263505074721}{1000000000}\right) \cdot x + \frac{156699607947}{500000000}\right) \cdot x + \frac{23533438303}{500000000}} \cdot \left(x - 2\right)} \]

      5. lower-*.f64 N/A

         \[\leadsto \color{blue}{\frac{\left(\left(\left(x \cdot \frac{104109730557}{25000000000} + \frac{393497462077}{5000000000}\right) \cdot x + \frac{4297481763}{31250000}\right) \cdot x + y\right) \cdot x + z}{\left(\left(\left(x + \frac{216700011257}{5000000000}\right) \cdot x + \frac{263505074721}{1000000000}\right) \cdot x + \frac{156699607947}{500000000}\right) \cdot x + \frac{23533438303}{500000000}} \cdot \left(x - 2\right)} \]

   4. Applied rewrites 0.0%

      \[\leadsto \color{blue}{\frac{\mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, 4.16438922228, 78.6994924154\right), 137.519416416\right), y\right), z\right)}{\mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, x + 43.3400022514, 263.505074721\right), 313.399215894\right), 47.066876606\right)} \cdot \left(x + -2\right)} \]

   5. Taylor expanded in x around -inf

      \[\leadsto \color{blue}{\left(\frac{104109730557}{25000000000} + -1 \cdot \frac{\frac{12723143231740136880149}{125000000000000000000} + -1 \cdot \frac{\frac{2157218858562374472887084159837293}{625000000000000000000000000000} + -1 \cdot \frac{\frac{387732519225574910908939577061312055388407301}{3125000000000000000000000000000000000000} + -1 \cdot y}{x}}{x}}{x}\right)} \cdot \left(x + -2\right) \]
   6. Step-by-step derivation

      1. mul-1-neg N/A

         \[\leadsto \left(\frac{104109730557}{25000000000} + \color{blue}{\left(\mathsf{neg}\left(\frac{\frac{12723143231740136880149}{125000000000000000000} + -1 \cdot \frac{\frac{2157218858562374472887084159837293}{625000000000000000000000000000} + -1 \cdot \frac{\frac{387732519225574910908939577061312055388407301}{3125000000000000000000000000000000000000} + -1 \cdot y}{x}}{x}}{x}\right)\right)}\right) \cdot \left(x + -2\right) \]

      2. unsub-neg N/A

         \[\leadsto \color{blue}{\left(\frac{104109730557}{25000000000} - \frac{\frac{12723143231740136880149}{125000000000000000000} + -1 \cdot \frac{\frac{2157218858562374472887084159837293}{625000000000000000000000000000} + -1 \cdot \frac{\frac{387732519225574910908939577061312055388407301}{3125000000000000000000000000000000000000} + -1 \cdot y}{x}}{x}}{x}\right)} \cdot \left(x + -2\right) \]

      3. lower--.f64 N/A

         \[\leadsto \color{blue}{\left(\frac{104109730557}{25000000000} - \frac{\frac{12723143231740136880149}{125000000000000000000} + -1 \cdot \frac{\frac{2157218858562374472887084159837293}{625000000000000000000000000000} + -1 \cdot \frac{\frac{387732519225574910908939577061312055388407301}{3125000000000000000000000000000000000000} + -1 \cdot y}{x}}{x}}{x}\right)} \cdot \left(x + -2\right) \]

      4. lower-/.f64 N/A

         \[\leadsto \left(\frac{104109730557}{25000000000} - \color{blue}{\frac{\frac{12723143231740136880149}{125000000000000000000} + -1 \cdot \frac{\frac{2157218858562374472887084159837293}{625000000000000000000000000000} + -1 \cdot \frac{\frac{387732519225574910908939577061312055388407301}{3125000000000000000000000000000000000000} + -1 \cdot y}{x}}{x}}{x}}\right) \cdot \left(x + -2\right) \]

   7. Applied rewrites 99.1%

      \[\leadsto \color{blue}{\left(4.16438922228 - \frac{101.7851458539211 - \frac{3451.550173699799 + \frac{y + -124074.40615218398}{x}}{x}}{x}\right)} \cdot \left(x + -2\right) \]
3. Recombined 2 regimes into one program.
4. Add Preprocessing

Alternative 2: 95.4% accurate, 1.1× speedup

Math

\[\begin{array}{l} \\ \begin{array}{l} t_0 := \frac{x + -2}{0.24013125253755718}\\ t_1 := \mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, 4.16438922228, 78.6994924154\right), 137.519416416\right), y\right), z\right)\\ \mathbf{if}\;x \leq -1.3 \cdot 10^{+63}:\\ \;\;\;\;t\_0\\ \mathbf{elif}\;x \leq -36:\\ \;\;\;\;t\_1 \cdot \frac{1}{x \cdot \left(x \cdot x\right)}\\ \mathbf{elif}\;x \leq 1.96 \cdot 10^{-14}:\\ \;\;\;\;\frac{t\_1}{\mathsf{fma}\left(x, 313.399215894, 47.066876606\right)} \cdot \left(x + -2\right)\\ \mathbf{elif}\;x \leq 4.1 \cdot 10^{+49}:\\ \;\;\;\;\frac{\mathsf{fma}\left(x, y, z\right)}{\mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, x + 43.3400022514, 263.505074721\right), 313.399215894\right), 47.066876606\right)} \cdot \left(x + -2\right)\\ \mathbf{else}:\\ \;\;\;\;t\_0\\ \end{array} \end{array} \]

FPCore

(FPCore (x y z)
 :precision binary64
 (let* ((t_0 (/ (+ x -2.0) 0.24013125253755718))
        (t_1
         (fma
          x
          (fma x (fma x (fma x 4.16438922228 78.6994924154) 137.519416416) y)
          z)))
   (if (<= x -1.3e+63)
     t_0
     (if (<= x -36.0)
       (* t_1 (/ 1.0 (* x (* x x))))
       (if (<= x 1.96e-14)
         (* (/ t_1 (fma x 313.399215894 47.066876606)) (+ x -2.0))
         (if (<= x 4.1e+49)
           (*
            (/
             (fma x y z)
             (fma
              x
              (fma x (fma x (+ x 43.3400022514) 263.505074721) 313.399215894)
              47.066876606))
            (+ x -2.0))
           t_0))))))

C

double code(double x, double y, double z) {
	double t_0 = (x + -2.0) / 0.24013125253755718;
	double t_1 = fma(x, fma(x, fma(x, fma(x, 4.16438922228, 78.6994924154), 137.519416416), y), z);
	double tmp;
	if (x <= -1.3e+63) {
		tmp = t_0;
	} else if (x <= -36.0) {
		tmp = t_1 * (1.0 / (x * (x * x)));
	} else if (x <= 1.96e-14) {
		tmp = (t_1 / fma(x, 313.399215894, 47.066876606)) * (x + -2.0);
	} else if (x <= 4.1e+49) {
		tmp = (fma(x, y, z) / fma(x, fma(x, fma(x, (x + 43.3400022514), 263.505074721), 313.399215894), 47.066876606)) * (x + -2.0);
	} else {
		tmp = t_0;
	}
	return tmp;
}

Julia

function code(x, y, z)
	t_0 = Float64(Float64(x + -2.0) / 0.24013125253755718)
	t_1 = fma(x, fma(x, fma(x, fma(x, 4.16438922228, 78.6994924154), 137.519416416), y), z)
	tmp = 0.0
	if (x <= -1.3e+63)
		tmp = t_0;
	elseif (x <= -36.0)
		tmp = Float64(t_1 * Float64(1.0 / Float64(x * Float64(x * x))));
	elseif (x <= 1.96e-14)
		tmp = Float64(Float64(t_1 / fma(x, 313.399215894, 47.066876606)) * Float64(x + -2.0));
	elseif (x <= 4.1e+49)
		tmp = Float64(Float64(fma(x, y, z) / fma(x, fma(x, fma(x, Float64(x + 43.3400022514), 263.505074721), 313.399215894), 47.066876606)) * Float64(x + -2.0));
	else
		tmp = t_0;
	end
	return tmp
end

Wolfram

code[x_, y_, z_] := Block[{t$95$0 = N[(N[(x + -2.0), $MachinePrecision] / 0.24013125253755718), $MachinePrecision]}, Block[{t$95$1 = N[(x * N[(x * N[(x * N[(x * 4.16438922228 + 78.6994924154), $MachinePrecision] + 137.519416416), $MachinePrecision] + y), $MachinePrecision] + z), $MachinePrecision]}, If[LessEqual[x, -1.3e+63], t$95$0, If[LessEqual[x, -36.0], N[(t$95$1 * N[(1.0 / N[(x * N[(x * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[x, 1.96e-14], N[(N[(t$95$1 / N[(x * 313.399215894 + 47.066876606), $MachinePrecision]), $MachinePrecision] * N[(x + -2.0), $MachinePrecision]), $MachinePrecision], If[LessEqual[x, 4.1e+49], N[(N[(N[(x * y + z), $MachinePrecision] / N[(x * N[(x * N[(x * N[(x + 43.3400022514), $MachinePrecision] + 263.505074721), $MachinePrecision] + 313.399215894), $MachinePrecision] + 47.066876606), $MachinePrecision]), $MachinePrecision] * N[(x + -2.0), $MachinePrecision]), $MachinePrecision], t$95$0]]]]]]

Derivation

1. Split input into 4 regimes

2. if x < -1.3000000000000001e63 or 4.1e49 < x

   1. Initial program 1.2%

      \[\frac{\left(x - 2\right) \cdot \left(\left(\left(\left(x \cdot 4.16438922228 + 78.6994924154\right) \cdot x + 137.519416416\right) \cdot x + y\right) \cdot x + z\right)}{\left(\left(\left(x + 43.3400022514\right) \cdot x + 263.505074721\right) \cdot x + 313.399215894\right) \cdot x + 47.066876606} \]
   2. Add Preprocessing

   4. Applied rewrites 7.1%

      \[\leadsto \color{blue}{\frac{x + -2}{\frac{\mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, x + 43.3400022514, 263.505074721\right), 313.399215894\right), 47.066876606\right)}{\mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, 4.16438922228, 78.6994924154\right), 137.519416416\right), y\right), z\right)}}} \]

   5. Taylor expanded in x around inf

      \[\leadsto \frac{x + -2}{\color{blue}{\frac{25000000000}{104109730557}}} \]
   6. Step-by-step derivation

   7. Applied rewrites 98.6%

      \[\leadsto \frac{x + -2}{\color{blue}{0.24013125253755718}} \]

   if -1.3000000000000001e63 < x < -36

   1. Initial program 72.0%

      \[\frac{\left(x - 2\right) \cdot \left(\left(\left(\left(x \cdot 4.16438922228 + 78.6994924154\right) \cdot x + 137.519416416\right) \cdot x + y\right) \cdot x + z\right)}{\left(\left(\left(x + 43.3400022514\right) \cdot x + 263.505074721\right) \cdot x + 313.399215894\right) \cdot x + 47.066876606} \]
   2. Add Preprocessing
   3. Step-by-step derivation

      1. lift-/.f64 N/A

         \[\leadsto \color{blue}{\frac{\left(x - 2\right) \cdot \left(\left(\left(\left(x \cdot \frac{104109730557}{25000000000} + \frac{393497462077}{5000000000}\right) \cdot x + \frac{4297481763}{31250000}\right) \cdot x + y\right) \cdot x + z\right)}{\left(\left(\left(x + \frac{216700011257}{5000000000}\right) \cdot x + \frac{263505074721}{1000000000}\right) \cdot x + \frac{156699607947}{500000000}\right) \cdot x + \frac{23533438303}{500000000}}} \]

      2. lift-*.f64 N/A

         \[\leadsto \frac{\color{blue}{\left(x - 2\right) \cdot \left(\left(\left(\left(x \cdot \frac{104109730557}{25000000000} + \frac{393497462077}{5000000000}\right) \cdot x + \frac{4297481763}{31250000}\right) \cdot x + y\right) \cdot x + z\right)}}{\left(\left(\left(x + \frac{216700011257}{5000000000}\right) \cdot x + \frac{263505074721}{1000000000}\right) \cdot x + \frac{156699607947}{500000000}\right) \cdot x + \frac{23533438303}{500000000}} \]

      3. *-commutative N/A

         \[\leadsto \frac{\color{blue}{\left(\left(\left(\left(x \cdot \frac{104109730557}{25000000000} + \frac{393497462077}{5000000000}\right) \cdot x + \frac{4297481763}{31250000}\right) \cdot x + y\right) \cdot x + z\right) \cdot \left(x - 2\right)}}{\left(\left(\left(x + \frac{216700011257}{5000000000}\right) \cdot x + \frac{263505074721}{1000000000}\right) \cdot x + \frac{156699607947}{500000000}\right) \cdot x + \frac{23533438303}{500000000}} \]

      4. associate-/l* N/A

         \[\leadsto \color{blue}{\left(\left(\left(\left(x \cdot \frac{104109730557}{25000000000} + \frac{393497462077}{5000000000}\right) \cdot x + \frac{4297481763}{31250000}\right) \cdot x + y\right) \cdot x + z\right) \cdot \frac{x - 2}{\left(\left(\left(x + \frac{216700011257}{5000000000}\right) \cdot x + \frac{263505074721}{1000000000}\right) \cdot x + \frac{156699607947}{500000000}\right) \cdot x + \frac{23533438303}{500000000}}} \]

      5. lower-*.f64 N/A

         \[\leadsto \color{blue}{\left(\left(\left(\left(x \cdot \frac{104109730557}{25000000000} + \frac{393497462077}{5000000000}\right) \cdot x + \frac{4297481763}{31250000}\right) \cdot x + y\right) \cdot x + z\right) \cdot \frac{x - 2}{\left(\left(\left(x + \frac{216700011257}{5000000000}\right) \cdot x + \frac{263505074721}{1000000000}\right) \cdot x + \frac{156699607947}{500000000}\right) \cdot x + \frac{23533438303}{500000000}}} \]

   4. Applied rewrites 92.3%

      \[\leadsto \color{blue}{\mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, 4.16438922228, 78.6994924154\right), 137.519416416\right), y\right), z\right) \cdot \frac{x + -2}{\mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, x + 43.3400022514, 263.505074721\right), 313.399215894\right), 47.066876606\right)}} \]

   5. Taylor expanded in x around inf

      \[\leadsto \mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, \frac{104109730557}{25000000000}, \frac{393497462077}{5000000000}\right), \frac{4297481763}{31250000}\right), y\right), z\right) \cdot \color{blue}{\frac{1}{{x}^{3}}} \]
   6. Step-by-step derivation

      1. lower-/.f64 N/A

         \[\leadsto \mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, \frac{104109730557}{25000000000}, \frac{393497462077}{5000000000}\right), \frac{4297481763}{31250000}\right), y\right), z\right) \cdot \color{blue}{\frac{1}{{x}^{3}}} \]

      2. cube-mult N/A

         \[\leadsto \mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, \frac{104109730557}{25000000000}, \frac{393497462077}{5000000000}\right), \frac{4297481763}{31250000}\right), y\right), z\right) \cdot \frac{1}{\color{blue}{x \cdot \left(x \cdot x\right)}} \]

      3. unpow2 N/A

         \[\leadsto \mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, \frac{104109730557}{25000000000}, \frac{393497462077}{5000000000}\right), \frac{4297481763}{31250000}\right), y\right), z\right) \cdot \frac{1}{x \cdot \color{blue}{{x}^{2}}} \]

      4. lower-*.f64 N/A

         \[\leadsto \mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, \frac{104109730557}{25000000000}, \frac{393497462077}{5000000000}\right), \frac{4297481763}{31250000}\right), y\right), z\right) \cdot \frac{1}{\color{blue}{x \cdot {x}^{2}}} \]

      5. unpow2 N/A

         \[\leadsto \mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, \frac{104109730557}{25000000000}, \frac{393497462077}{5000000000}\right), \frac{4297481763}{31250000}\right), y\right), z\right) \cdot \frac{1}{x \cdot \color{blue}{\left(x \cdot x\right)}} \]

      6. lower-*.f64 90.0

         \[\leadsto \mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, 4.16438922228, 78.6994924154\right), 137.519416416\right), y\right), z\right) \cdot \frac{1}{x \cdot \color{blue}{\left(x \cdot x\right)}} \]

   7. Applied rewrites 90.0%

      \[\leadsto \mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, 4.16438922228, 78.6994924154\right), 137.519416416\right), y\right), z\right) \cdot \color{blue}{\frac{1}{x \cdot \left(x \cdot x\right)}} \]

   if -36 < x < 1.96e-14

   1. Initial program 99.8%

      \[\frac{\left(x - 2\right) \cdot \left(\left(\left(\left(x \cdot 4.16438922228 + 78.6994924154\right) \cdot x + 137.519416416\right) \cdot x + y\right) \cdot x + z\right)}{\left(\left(\left(x + 43.3400022514\right) \cdot x + 263.505074721\right) \cdot x + 313.399215894\right) \cdot x + 47.066876606} \]
   2. Add Preprocessing
   3. Step-by-step derivation

      1. lift-/.f64 N/A

         \[\leadsto \color{blue}{\frac{\left(x - 2\right) \cdot \left(\left(\left(\left(x \cdot \frac{104109730557}{25000000000} + \frac{393497462077}{5000000000}\right) \cdot x + \frac{4297481763}{31250000}\right) \cdot x + y\right) \cdot x + z\right)}{\left(\left(\left(x + \frac{216700011257}{5000000000}\right) \cdot x + \frac{263505074721}{1000000000}\right) \cdot x + \frac{156699607947}{500000000}\right) \cdot x + \frac{23533438303}{500000000}}} \]

      2. lift-*.f64 N/A

         \[\leadsto \frac{\color{blue}{\left(x - 2\right) \cdot \left(\left(\left(\left(x \cdot \frac{104109730557}{25000000000} + \frac{393497462077}{5000000000}\right) \cdot x + \frac{4297481763}{31250000}\right) \cdot x + y\right) \cdot x + z\right)}}{\left(\left(\left(x + \frac{216700011257}{5000000000}\right) \cdot x + \frac{263505074721}{1000000000}\right) \cdot x + \frac{156699607947}{500000000}\right) \cdot x + \frac{23533438303}{500000000}} \]

      3. associate-/l* N/A

         \[\leadsto \color{blue}{\left(x - 2\right) \cdot \frac{\left(\left(\left(x \cdot \frac{104109730557}{25000000000} + \frac{393497462077}{5000000000}\right) \cdot x + \frac{4297481763}{31250000}\right) \cdot x + y\right) \cdot x + z}{\left(\left(\left(x + \frac{216700011257}{5000000000}\right) \cdot x + \frac{263505074721}{1000000000}\right) \cdot x + \frac{156699607947}{500000000}\right) \cdot x + \frac{23533438303}{500000000}}} \]

      4. *-commutative N/A

         \[\leadsto \color{blue}{\frac{\left(\left(\left(x \cdot \frac{104109730557}{25000000000} + \frac{393497462077}{5000000000}\right) \cdot x + \frac{4297481763}{31250000}\right) \cdot x + y\right) \cdot x + z}{\left(\left(\left(x + \frac{216700011257}{5000000000}\right) \cdot x + \frac{263505074721}{1000000000}\right) \cdot x + \frac{156699607947}{500000000}\right) \cdot x + \frac{23533438303}{500000000}} \cdot \left(x - 2\right)} \]

      5. lower-*.f64 N/A

         \[\leadsto \color{blue}{\frac{\left(\left(\left(x \cdot \frac{104109730557}{25000000000} + \frac{393497462077}{5000000000}\right) \cdot x + \frac{4297481763}{31250000}\right) \cdot x + y\right) \cdot x + z}{\left(\left(\left(x + \frac{216700011257}{5000000000}\right) \cdot x + \frac{263505074721}{1000000000}\right) \cdot x + \frac{156699607947}{500000000}\right) \cdot x + \frac{23533438303}{500000000}} \cdot \left(x - 2\right)} \]

   4. Applied rewrites 99.8%

      \[\leadsto \color{blue}{\frac{\mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, 4.16438922228, 78.6994924154\right), 137.519416416\right), y\right), z\right)}{\mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, x + 43.3400022514, 263.505074721\right), 313.399215894\right), 47.066876606\right)} \cdot \left(x + -2\right)} \]

   5. Taylor expanded in x around 0

      \[\leadsto \frac{\mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, \frac{104109730557}{25000000000}, \frac{393497462077}{5000000000}\right), \frac{4297481763}{31250000}\right), y\right), z\right)}{\mathsf{fma}\left(x, \color{blue}{\frac{156699607947}{500000000}}, \frac{23533438303}{500000000}\right)} \cdot \left(x + -2\right) \]
   6. Step-by-step derivation

   7. Applied rewrites 99.6%

      \[\leadsto \frac{\mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, 4.16438922228, 78.6994924154\right), 137.519416416\right), y\right), z\right)}{\mathsf{fma}\left(x, \color{blue}{313.399215894}, 47.066876606\right)} \cdot \left(x + -2\right) \]

   if 1.96e-14 < x < 4.1e49

   1. Initial program 87.7%

      \[\frac{\left(x - 2\right) \cdot \left(\left(\left(\left(x \cdot 4.16438922228 + 78.6994924154\right) \cdot x + 137.519416416\right) \cdot x + y\right) \cdot x + z\right)}{\left(\left(\left(x + 43.3400022514\right) \cdot x + 263.505074721\right) \cdot x + 313.399215894\right) \cdot x + 47.066876606} \]
   2. Add Preprocessing
   3. Step-by-step derivation

      1. lift-/.f64 N/A

         \[\leadsto \color{blue}{\frac{\left(x - 2\right) \cdot \left(\left(\left(\left(x \cdot \frac{104109730557}{25000000000} + \frac{393497462077}{5000000000}\right) \cdot x + \frac{4297481763}{31250000}\right) \cdot x + y\right) \cdot x + z\right)}{\left(\left(\left(x + \frac{216700011257}{5000000000}\right) \cdot x + \frac{263505074721}{1000000000}\right) \cdot x + \frac{156699607947}{500000000}\right) \cdot x + \frac{23533438303}{500000000}}} \]

      2. lift-*.f64 N/A

         \[\leadsto \frac{\color{blue}{\left(x - 2\right) \cdot \left(\left(\left(\left(x \cdot \frac{104109730557}{25000000000} + \frac{393497462077}{5000000000}\right) \cdot x + \frac{4297481763}{31250000}\right) \cdot x + y\right) \cdot x + z\right)}}{\left(\left(\left(x + \frac{216700011257}{5000000000}\right) \cdot x + \frac{263505074721}{1000000000}\right) \cdot x + \frac{156699607947}{500000000}\right) \cdot x + \frac{23533438303}{500000000}} \]

      3. associate-/l* N/A

         \[\leadsto \color{blue}{\left(x - 2\right) \cdot \frac{\left(\left(\left(x \cdot \frac{104109730557}{25000000000} + \frac{393497462077}{5000000000}\right) \cdot x + \frac{4297481763}{31250000}\right) \cdot x + y\right) \cdot x + z}{\left(\left(\left(x + \frac{216700011257}{5000000000}\right) \cdot x + \frac{263505074721}{1000000000}\right) \cdot x + \frac{156699607947}{500000000}\right) \cdot x + \frac{23533438303}{500000000}}} \]

      4. *-commutative N/A

         \[\leadsto \color{blue}{\frac{\left(\left(\left(x \cdot \frac{104109730557}{25000000000} + \frac{393497462077}{5000000000}\right) \cdot x + \frac{4297481763}{31250000}\right) \cdot x + y\right) \cdot x + z}{\left(\left(\left(x + \frac{216700011257}{5000000000}\right) \cdot x + \frac{263505074721}{1000000000}\right) \cdot x + \frac{156699607947}{500000000}\right) \cdot x + \frac{23533438303}{500000000}} \cdot \left(x - 2\right)} \]

      5. lower-*.f64 N/A

         \[\leadsto \color{blue}{\frac{\left(\left(\left(x \cdot \frac{104109730557}{25000000000} + \frac{393497462077}{5000000000}\right) \cdot x + \frac{4297481763}{31250000}\right) \cdot x + y\right) \cdot x + z}{\left(\left(\left(x + \frac{216700011257}{5000000000}\right) \cdot x + \frac{263505074721}{1000000000}\right) \cdot x + \frac{156699607947}{500000000}\right) \cdot x + \frac{23533438303}{500000000}} \cdot \left(x - 2\right)} \]

   4. Applied rewrites 99.0%

      \[\leadsto \color{blue}{\frac{\mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, 4.16438922228, 78.6994924154\right), 137.519416416\right), y\right), z\right)}{\mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, x + 43.3400022514, 263.505074721\right), 313.399215894\right), 47.066876606\right)} \cdot \left(x + -2\right)} \]

   5. Taylor expanded in x around 0

      \[\leadsto \frac{\color{blue}{z + x \cdot y}}{\mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, x + \frac{216700011257}{5000000000}, \frac{263505074721}{1000000000}\right), \frac{156699607947}{500000000}\right), \frac{23533438303}{500000000}\right)} \cdot \left(x + -2\right) \]
   6. Step-by-step derivation

      1. +-commutative N/A

         \[\leadsto \frac{\color{blue}{x \cdot y + z}}{\mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, x + \frac{216700011257}{5000000000}, \frac{263505074721}{1000000000}\right), \frac{156699607947}{500000000}\right), \frac{23533438303}{500000000}\right)} \cdot \left(x + -2\right) \]

      2. lower-fma.f64 99.0

         \[\leadsto \frac{\color{blue}{\mathsf{fma}\left(x, y, z\right)}}{\mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, x + 43.3400022514, 263.505074721\right), 313.399215894\right), 47.066876606\right)} \cdot \left(x + -2\right) \]

   7. Applied rewrites 99.0%

      \[\leadsto \frac{\color{blue}{\mathsf{fma}\left(x, y, z\right)}}{\mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, x + 43.3400022514, 263.505074721\right), 313.399215894\right), 47.066876606\right)} \cdot \left(x + -2\right) \]
3. Recombined 4 regimes into one program.
4. Add Preprocessing

Alternative 3: 95.4% accurate, 1.1× speedup

Math

\[\begin{array}{l} \\ \begin{array}{l} \mathbf{if}\;x \leq -6.8 \cdot 10^{+28}:\\ \;\;\;\;\left(4.16438922228 - \frac{101.7851458539211 - \frac{3451.550173699799 + \frac{y + -124074.40615218398}{x}}{x}}{x}\right) \cdot \left(x + -2\right)\\ \mathbf{elif}\;x \leq 1.45 \cdot 10^{+37}:\\ \;\;\;\;\frac{\left(x - 2\right) \cdot \mathsf{fma}\left(x, \mathsf{fma}\left(x, 137.519416416, y\right), z\right)}{\left(\left(\left(x + 43.3400022514\right) \cdot x + 263.505074721\right) \cdot x + 313.399215894\right) \cdot x + 47.066876606}\\ \mathbf{else}:\\ \;\;\;\;\frac{x + -2}{0.24013125253755718}\\ \end{array} \end{array} \]

FPCore

(FPCore (x y z)
 :precision binary64
 (if (<= x -6.8e+28)
   (*
    (-
     4.16438922228
     (/
      (-
       101.7851458539211
       (/ (+ 3451.550173699799 (/ (+ y -124074.40615218398) x)) x))
      x))
    (+ x -2.0))
   (if (<= x 1.45e+37)
     (/
      (* (- x 2.0) (fma x (fma x 137.519416416 y) z))
      (+
       (*
        (+ (* (+ (* (+ x 43.3400022514) x) 263.505074721) x) 313.399215894)
        x)
       47.066876606))
     (/ (+ x -2.0) 0.24013125253755718))))

C

double code(double x, double y, double z) {
	double tmp;
	if (x <= -6.8e+28) {
		tmp = (4.16438922228 - ((101.7851458539211 - ((3451.550173699799 + ((y + -124074.40615218398) / x)) / x)) / x)) * (x + -2.0);
	} else if (x <= 1.45e+37) {
		tmp = ((x - 2.0) * fma(x, fma(x, 137.519416416, y), z)) / (((((((x + 43.3400022514) * x) + 263.505074721) * x) + 313.399215894) * x) + 47.066876606);
	} else {
		tmp = (x + -2.0) / 0.24013125253755718;
	}
	return tmp;
}

Julia

function code(x, y, z)
	tmp = 0.0
	if (x <= -6.8e+28)
		tmp = Float64(Float64(4.16438922228 - Float64(Float64(101.7851458539211 - Float64(Float64(3451.550173699799 + Float64(Float64(y + -124074.40615218398) / x)) / x)) / x)) * Float64(x + -2.0));
	elseif (x <= 1.45e+37)
		tmp = Float64(Float64(Float64(x - 2.0) * fma(x, fma(x, 137.519416416, y), z)) / Float64(Float64(Float64(Float64(Float64(Float64(Float64(x + 43.3400022514) * x) + 263.505074721) * x) + 313.399215894) * x) + 47.066876606));
	else
		tmp = Float64(Float64(x + -2.0) / 0.24013125253755718);
	end
	return tmp
end

Wolfram

code[x_, y_, z_] := If[LessEqual[x, -6.8e+28], N[(N[(4.16438922228 - N[(N[(101.7851458539211 - N[(N[(3451.550173699799 + N[(N[(y + -124074.40615218398), $MachinePrecision] / x), $MachinePrecision]), $MachinePrecision] / x), $MachinePrecision]), $MachinePrecision] / x), $MachinePrecision]), $MachinePrecision] * N[(x + -2.0), $MachinePrecision]), $MachinePrecision], If[LessEqual[x, 1.45e+37], N[(N[(N[(x - 2.0), $MachinePrecision] * N[(x * N[(x * 137.519416416 + y), $MachinePrecision] + z), $MachinePrecision]), $MachinePrecision] / N[(N[(N[(N[(N[(N[(N[(x + 43.3400022514), $MachinePrecision] * x), $MachinePrecision] + 263.505074721), $MachinePrecision] * x), $MachinePrecision] + 313.399215894), $MachinePrecision] * x), $MachinePrecision] + 47.066876606), $MachinePrecision]), $MachinePrecision], N[(N[(x + -2.0), $MachinePrecision] / 0.24013125253755718), $MachinePrecision]]]

Derivation

1. Split input into 3 regimes

2. if x < -6.8e28

   1. Initial program 13.9%

      \[\frac{\left(x - 2\right) \cdot \left(\left(\left(\left(x \cdot 4.16438922228 + 78.6994924154\right) \cdot x + 137.519416416\right) \cdot x + y\right) \cdot x + z\right)}{\left(\left(\left(x + 43.3400022514\right) \cdot x + 263.505074721\right) \cdot x + 313.399215894\right) \cdot x + 47.066876606} \]
   2. Add Preprocessing
   3. Step-by-step derivation

      1. lift-/.f64 N/A

         \[\leadsto \color{blue}{\frac{\left(x - 2\right) \cdot \left(\left(\left(\left(x \cdot \frac{104109730557}{25000000000} + \frac{393497462077}{5000000000}\right) \cdot x + \frac{4297481763}{31250000}\right) \cdot x + y\right) \cdot x + z\right)}{\left(\left(\left(x + \frac{216700011257}{5000000000}\right) \cdot x + \frac{263505074721}{1000000000}\right) \cdot x + \frac{156699607947}{500000000}\right) \cdot x + \frac{23533438303}{500000000}}} \]

      2. lift-*.f64 N/A

         \[\leadsto \frac{\color{blue}{\left(x - 2\right) \cdot \left(\left(\left(\left(x \cdot \frac{104109730557}{25000000000} + \frac{393497462077}{5000000000}\right) \cdot x + \frac{4297481763}{31250000}\right) \cdot x + y\right) \cdot x + z\right)}}{\left(\left(\left(x + \frac{216700011257}{5000000000}\right) \cdot x + \frac{263505074721}{1000000000}\right) \cdot x + \frac{156699607947}{500000000}\right) \cdot x + \frac{23533438303}{500000000}} \]

      3. associate-/l* N/A

         \[\leadsto \color{blue}{\left(x - 2\right) \cdot \frac{\left(\left(\left(x \cdot \frac{104109730557}{25000000000} + \frac{393497462077}{5000000000}\right) \cdot x + \frac{4297481763}{31250000}\right) \cdot x + y\right) \cdot x + z}{\left(\left(\left(x + \frac{216700011257}{5000000000}\right) \cdot x + \frac{263505074721}{1000000000}\right) \cdot x + \frac{156699607947}{500000000}\right) \cdot x + \frac{23533438303}{500000000}}} \]

      4. *-commutative N/A

         \[\leadsto \color{blue}{\frac{\left(\left(\left(x \cdot \frac{104109730557}{25000000000} + \frac{393497462077}{5000000000}\right) \cdot x + \frac{4297481763}{31250000}\right) \cdot x + y\right) \cdot x + z}{\left(\left(\left(x + \frac{216700011257}{5000000000}\right) \cdot x + \frac{263505074721}{1000000000}\right) \cdot x + \frac{156699607947}{500000000}\right) \cdot x + \frac{23533438303}{500000000}} \cdot \left(x - 2\right)} \]

      5. lower-*.f64 N/A

         \[\leadsto \color{blue}{\frac{\left(\left(\left(x \cdot \frac{104109730557}{25000000000} + \frac{393497462077}{5000000000}\right) \cdot x + \frac{4297481763}{31250000}\right) \cdot x + y\right) \cdot x + z}{\left(\left(\left(x + \frac{216700011257}{5000000000}\right) \cdot x + \frac{263505074721}{1000000000}\right) \cdot x + \frac{156699607947}{500000000}\right) \cdot x + \frac{23533438303}{500000000}} \cdot \left(x - 2\right)} \]

   4. Applied rewrites 23.5%

      \[\leadsto \color{blue}{\frac{\mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, 4.16438922228, 78.6994924154\right), 137.519416416\right), y\right), z\right)}{\mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, x + 43.3400022514, 263.505074721\right), 313.399215894\right), 47.066876606\right)} \cdot \left(x + -2\right)} \]

   5. Taylor expanded in x around -inf

      \[\leadsto \color{blue}{\left(\frac{104109730557}{25000000000} + -1 \cdot \frac{\frac{12723143231740136880149}{125000000000000000000} + -1 \cdot \frac{\frac{2157218858562374472887084159837293}{625000000000000000000000000000} + -1 \cdot \frac{\frac{387732519225574910908939577061312055388407301}{3125000000000000000000000000000000000000} + -1 \cdot y}{x}}{x}}{x}\right)} \cdot \left(x + -2\right) \]
   6. Step-by-step derivation

      1. mul-1-neg N/A

         \[\leadsto \left(\frac{104109730557}{25000000000} + \color{blue}{\left(\mathsf{neg}\left(\frac{\frac{12723143231740136880149}{125000000000000000000} + -1 \cdot \frac{\frac{2157218858562374472887084159837293}{625000000000000000000000000000} + -1 \cdot \frac{\frac{387732519225574910908939577061312055388407301}{3125000000000000000000000000000000000000} + -1 \cdot y}{x}}{x}}{x}\right)\right)}\right) \cdot \left(x + -2\right) \]

      2. unsub-neg N/A

         \[\leadsto \color{blue}{\left(\frac{104109730557}{25000000000} - \frac{\frac{12723143231740136880149}{125000000000000000000} + -1 \cdot \frac{\frac{2157218858562374472887084159837293}{625000000000000000000000000000} + -1 \cdot \frac{\frac{387732519225574910908939577061312055388407301}{3125000000000000000000000000000000000000} + -1 \cdot y}{x}}{x}}{x}\right)} \cdot \left(x + -2\right) \]

      3. lower--.f64 N/A

         \[\leadsto \color{blue}{\left(\frac{104109730557}{25000000000} - \frac{\frac{12723143231740136880149}{125000000000000000000} + -1 \cdot \frac{\frac{2157218858562374472887084159837293}{625000000000000000000000000000} + -1 \cdot \frac{\frac{387732519225574910908939577061312055388407301}{3125000000000000000000000000000000000000} + -1 \cdot y}{x}}{x}}{x}\right)} \cdot \left(x + -2\right) \]

      4. lower-/.f64 N/A

         \[\leadsto \left(\frac{104109730557}{25000000000} - \color{blue}{\frac{\frac{12723143231740136880149}{125000000000000000000} + -1 \cdot \frac{\frac{2157218858562374472887084159837293}{625000000000000000000000000000} + -1 \cdot \frac{\frac{387732519225574910908939577061312055388407301}{3125000000000000000000000000000000000000} + -1 \cdot y}{x}}{x}}{x}}\right) \cdot \left(x + -2\right) \]

   7. Applied rewrites 96.4%

      \[\leadsto \color{blue}{\left(4.16438922228 - \frac{101.7851458539211 - \frac{3451.550173699799 + \frac{y + -124074.40615218398}{x}}{x}}{x}\right)} \cdot \left(x + -2\right) \]

   if -6.8e28 < x < 1.44999999999999989e37

   1. Initial program 99.8%

      \[\frac{\left(x - 2\right) \cdot \left(\left(\left(\left(x \cdot 4.16438922228 + 78.6994924154\right) \cdot x + 137.519416416\right) \cdot x + y\right) \cdot x + z\right)}{\left(\left(\left(x + 43.3400022514\right) \cdot x + 263.505074721\right) \cdot x + 313.399215894\right) \cdot x + 47.066876606} \]
   2. Add Preprocessing

   3. Taylor expanded in x around 0

      \[\leadsto \frac{\left(x - 2\right) \cdot \color{blue}{\left(z + x \cdot \left(y + \frac{4297481763}{31250000} \cdot x\right)\right)}}{\left(\left(\left(x + \frac{216700011257}{5000000000}\right) \cdot x + \frac{263505074721}{1000000000}\right) \cdot x + \frac{156699607947}{500000000}\right) \cdot x + \frac{23533438303}{500000000}} \]
   4. Step-by-step derivation

      1. +-commutative N/A

         \[\leadsto \frac{\left(x - 2\right) \cdot \color{blue}{\left(x \cdot \left(y + \frac{4297481763}{31250000} \cdot x\right) + z\right)}}{\left(\left(\left(x + \frac{216700011257}{5000000000}\right) \cdot x + \frac{263505074721}{1000000000}\right) \cdot x + \frac{156699607947}{500000000}\right) \cdot x + \frac{23533438303}{500000000}} \]

      2. lower-fma.f64 N/A

         \[\leadsto \frac{\left(x - 2\right) \cdot \color{blue}{\mathsf{fma}\left(x, y + \frac{4297481763}{31250000} \cdot x, z\right)}}{\left(\left(\left(x + \frac{216700011257}{5000000000}\right) \cdot x + \frac{263505074721}{1000000000}\right) \cdot x + \frac{156699607947}{500000000}\right) \cdot x + \frac{23533438303}{500000000}} \]

      3. +-commutative N/A

         \[\leadsto \frac{\left(x - 2\right) \cdot \mathsf{fma}\left(x, \color{blue}{\frac{4297481763}{31250000} \cdot x + y}, z\right)}{\left(\left(\left(x + \frac{216700011257}{5000000000}\right) \cdot x + \frac{263505074721}{1000000000}\right) \cdot x + \frac{156699607947}{500000000}\right) \cdot x + \frac{23533438303}{500000000}} \]

      4. *-commutative N/A

         \[\leadsto \frac{\left(x - 2\right) \cdot \mathsf{fma}\left(x, \color{blue}{x \cdot \frac{4297481763}{31250000}} + y, z\right)}{\left(\left(\left(x + \frac{216700011257}{5000000000}\right) \cdot x + \frac{263505074721}{1000000000}\right) \cdot x + \frac{156699607947}{500000000}\right) \cdot x + \frac{23533438303}{500000000}} \]

      5. lower-fma.f64 99.7

         \[\leadsto \frac{\left(x - 2\right) \cdot \mathsf{fma}\left(x, \color{blue}{\mathsf{fma}\left(x, 137.519416416, y\right)}, z\right)}{\left(\left(\left(x + 43.3400022514\right) \cdot x + 263.505074721\right) \cdot x + 313.399215894\right) \cdot x + 47.066876606} \]

   5. Applied rewrites 99.7%

      \[\leadsto \frac{\left(x - 2\right) \cdot \color{blue}{\mathsf{fma}\left(x, \mathsf{fma}\left(x, 137.519416416, y\right), z\right)}}{\left(\left(\left(x + 43.3400022514\right) \cdot x + 263.505074721\right) \cdot x + 313.399215894\right) \cdot x + 47.066876606} \]

   if 1.44999999999999989e37 < x

   1. Initial program 2.2%

      \[\frac{\left(x - 2\right) \cdot \left(\left(\left(\left(x \cdot 4.16438922228 + 78.6994924154\right) \cdot x + 137.519416416\right) \cdot x + y\right) \cdot x + z\right)}{\left(\left(\left(x + 43.3400022514\right) \cdot x + 263.505074721\right) \cdot x + 313.399215894\right) \cdot x + 47.066876606} \]
   2. Add Preprocessing

   4. Applied rewrites 9.7%

      \[\leadsto \color{blue}{\frac{x + -2}{\frac{\mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, x + 43.3400022514, 263.505074721\right), 313.399215894\right), 47.066876606\right)}{\mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, 4.16438922228, 78.6994924154\right), 137.519416416\right), y\right), z\right)}}} \]

   5. Taylor expanded in x around inf

      \[\leadsto \frac{x + -2}{\color{blue}{\frac{25000000000}{104109730557}}} \]
   6. Step-by-step derivation

   7. Applied rewrites 97.6%

      \[\leadsto \frac{x + -2}{\color{blue}{0.24013125253755718}} \]
3. Recombined 3 regimes into one program.
4. Add Preprocessing

Alternative 4: 96.0% accurate, 1.2× speedup

Math

\[\begin{array}{l} \\ \begin{array}{l} t_0 := \left(4.16438922228 - \frac{101.7851458539211 - \frac{3451.550173699799 + \frac{y + -124074.40615218398}{x}}{x}}{x}\right) \cdot \left(x + -2\right)\\ \mathbf{if}\;x \leq -36:\\ \;\;\;\;t\_0\\ \mathbf{elif}\;x \leq 48:\\ \;\;\;\;\frac{\mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, 4.16438922228, 78.6994924154\right), 137.519416416\right), y\right), z\right)}{\mathsf{fma}\left(x, 313.399215894, 47.066876606\right)} \cdot \left(x + -2\right)\\ \mathbf{else}:\\ \;\;\;\;t\_0\\ \end{array} \end{array} \]

FPCore

(FPCore (x y z)
 :precision binary64
 (let* ((t_0
         (*
          (-
           4.16438922228
           (/
            (-
             101.7851458539211
             (/ (+ 3451.550173699799 (/ (+ y -124074.40615218398) x)) x))
            x))
          (+ x -2.0))))
   (if (<= x -36.0)
     t_0
     (if (<= x 48.0)
       (*
        (/
         (fma
          x
          (fma x (fma x (fma x 4.16438922228 78.6994924154) 137.519416416) y)
          z)
         (fma x 313.399215894 47.066876606))
        (+ x -2.0))
       t_0))))

C

double code(double x, double y, double z) {
	double t_0 = (4.16438922228 - ((101.7851458539211 - ((3451.550173699799 + ((y + -124074.40615218398) / x)) / x)) / x)) * (x + -2.0);
	double tmp;
	if (x <= -36.0) {
		tmp = t_0;
	} else if (x <= 48.0) {
		tmp = (fma(x, fma(x, fma(x, fma(x, 4.16438922228, 78.6994924154), 137.519416416), y), z) / fma(x, 313.399215894, 47.066876606)) * (x + -2.0);
	} else {
		tmp = t_0;
	}
	return tmp;
}

Julia

function code(x, y, z)
	t_0 = Float64(Float64(4.16438922228 - Float64(Float64(101.7851458539211 - Float64(Float64(3451.550173699799 + Float64(Float64(y + -124074.40615218398) / x)) / x)) / x)) * Float64(x + -2.0))
	tmp = 0.0
	if (x <= -36.0)
		tmp = t_0;
	elseif (x <= 48.0)
		tmp = Float64(Float64(fma(x, fma(x, fma(x, fma(x, 4.16438922228, 78.6994924154), 137.519416416), y), z) / fma(x, 313.399215894, 47.066876606)) * Float64(x + -2.0));
	else
		tmp = t_0;
	end
	return tmp
end

Wolfram

code[x_, y_, z_] := Block[{t$95$0 = N[(N[(4.16438922228 - N[(N[(101.7851458539211 - N[(N[(3451.550173699799 + N[(N[(y + -124074.40615218398), $MachinePrecision] / x), $MachinePrecision]), $MachinePrecision] / x), $MachinePrecision]), $MachinePrecision] / x), $MachinePrecision]), $MachinePrecision] * N[(x + -2.0), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[x, -36.0], t$95$0, If[LessEqual[x, 48.0], N[(N[(N[(x * N[(x * N[(x * N[(x * 4.16438922228 + 78.6994924154), $MachinePrecision] + 137.519416416), $MachinePrecision] + y), $MachinePrecision] + z), $MachinePrecision] / N[(x * 313.399215894 + 47.066876606), $MachinePrecision]), $MachinePrecision] * N[(x + -2.0), $MachinePrecision]), $MachinePrecision], t$95$0]]]

Derivation

1. Split input into 2 regimes

2. if x < -36 or 48 < x

   1. Initial program 13.2%

      \[\frac{\left(x - 2\right) \cdot \left(\left(\left(\left(x \cdot 4.16438922228 + 78.6994924154\right) \cdot x + 137.519416416\right) \cdot x + y\right) \cdot x + z\right)}{\left(\left(\left(x + 43.3400022514\right) \cdot x + 263.505074721\right) \cdot x + 313.399215894\right) \cdot x + 47.066876606} \]
   2. Add Preprocessing
   3. Step-by-step derivation

      1. lift-/.f64 N/A

         \[\leadsto \color{blue}{\frac{\left(x - 2\right) \cdot \left(\left(\left(\left(x \cdot \frac{104109730557}{25000000000} + \frac{393497462077}{5000000000}\right) \cdot x + \frac{4297481763}{31250000}\right) \cdot x + y\right) \cdot x + z\right)}{\left(\left(\left(x + \frac{216700011257}{5000000000}\right) \cdot x + \frac{263505074721}{1000000000}\right) \cdot x + \frac{156699607947}{500000000}\right) \cdot x + \frac{23533438303}{500000000}}} \]

      2. lift-*.f64 N/A

         \[\leadsto \frac{\color{blue}{\left(x - 2\right) \cdot \left(\left(\left(\left(x \cdot \frac{104109730557}{25000000000} + \frac{393497462077}{5000000000}\right) \cdot x + \frac{4297481763}{31250000}\right) \cdot x + y\right) \cdot x + z\right)}}{\left(\left(\left(x + \frac{216700011257}{5000000000}\right) \cdot x + \frac{263505074721}{1000000000}\right) \cdot x + \frac{156699607947}{500000000}\right) \cdot x + \frac{23533438303}{500000000}} \]

      3. associate-/l* N/A

         \[\leadsto \color{blue}{\left(x - 2\right) \cdot \frac{\left(\left(\left(x \cdot \frac{104109730557}{25000000000} + \frac{393497462077}{5000000000}\right) \cdot x + \frac{4297481763}{31250000}\right) \cdot x + y\right) \cdot x + z}{\left(\left(\left(x + \frac{216700011257}{5000000000}\right) \cdot x + \frac{263505074721}{1000000000}\right) \cdot x + \frac{156699607947}{500000000}\right) \cdot x + \frac{23533438303}{500000000}}} \]

      4. *-commutative N/A

         \[\leadsto \color{blue}{\frac{\left(\left(\left(x \cdot \frac{104109730557}{25000000000} + \frac{393497462077}{5000000000}\right) \cdot x + \frac{4297481763}{31250000}\right) \cdot x + y\right) \cdot x + z}{\left(\left(\left(x + \frac{216700011257}{5000000000}\right) \cdot x + \frac{263505074721}{1000000000}\right) \cdot x + \frac{156699607947}{500000000}\right) \cdot x + \frac{23533438303}{500000000}} \cdot \left(x - 2\right)} \]

      5. lower-*.f64 N/A

         \[\leadsto \color{blue}{\frac{\left(\left(\left(x \cdot \frac{104109730557}{25000000000} + \frac{393497462077}{5000000000}\right) \cdot x + \frac{4297481763}{31250000}\right) \cdot x + y\right) \cdot x + z}{\left(\left(\left(x + \frac{216700011257}{5000000000}\right) \cdot x + \frac{263505074721}{1000000000}\right) \cdot x + \frac{156699607947}{500000000}\right) \cdot x + \frac{23533438303}{500000000}} \cdot \left(x - 2\right)} \]

   4. Applied rewrites 21.3%

      \[\leadsto \color{blue}{\frac{\mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, 4.16438922228, 78.6994924154\right), 137.519416416\right), y\right), z\right)}{\mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, x + 43.3400022514, 263.505074721\right), 313.399215894\right), 47.066876606\right)} \cdot \left(x + -2\right)} \]

   5. Taylor expanded in x around -inf

      \[\leadsto \color{blue}{\left(\frac{104109730557}{25000000000} + -1 \cdot \frac{\frac{12723143231740136880149}{125000000000000000000} + -1 \cdot \frac{\frac{2157218858562374472887084159837293}{625000000000000000000000000000} + -1 \cdot \frac{\frac{387732519225574910908939577061312055388407301}{3125000000000000000000000000000000000000} + -1 \cdot y}{x}}{x}}{x}\right)} \cdot \left(x + -2\right) \]
   6. Step-by-step derivation

      1. mul-1-neg N/A

         \[\leadsto \left(\frac{104109730557}{25000000000} + \color{blue}{\left(\mathsf{neg}\left(\frac{\frac{12723143231740136880149}{125000000000000000000} + -1 \cdot \frac{\frac{2157218858562374472887084159837293}{625000000000000000000000000000} + -1 \cdot \frac{\frac{387732519225574910908939577061312055388407301}{3125000000000000000000000000000000000000} + -1 \cdot y}{x}}{x}}{x}\right)\right)}\right) \cdot \left(x + -2\right) \]

      2. unsub-neg N/A

         \[\leadsto \color{blue}{\left(\frac{104109730557}{25000000000} - \frac{\frac{12723143231740136880149}{125000000000000000000} + -1 \cdot \frac{\frac{2157218858562374472887084159837293}{625000000000000000000000000000} + -1 \cdot \frac{\frac{387732519225574910908939577061312055388407301}{3125000000000000000000000000000000000000} + -1 \cdot y}{x}}{x}}{x}\right)} \cdot \left(x + -2\right) \]

      3. lower--.f64 N/A

         \[\leadsto \color{blue}{\left(\frac{104109730557}{25000000000} - \frac{\frac{12723143231740136880149}{125000000000000000000} + -1 \cdot \frac{\frac{2157218858562374472887084159837293}{625000000000000000000000000000} + -1 \cdot \frac{\frac{387732519225574910908939577061312055388407301}{3125000000000000000000000000000000000000} + -1 \cdot y}{x}}{x}}{x}\right)} \cdot \left(x + -2\right) \]

      4. lower-/.f64 N/A

         \[\leadsto \left(\frac{104109730557}{25000000000} - \color{blue}{\frac{\frac{12723143231740136880149}{125000000000000000000} + -1 \cdot \frac{\frac{2157218858562374472887084159837293}{625000000000000000000000000000} + -1 \cdot \frac{\frac{387732519225574910908939577061312055388407301}{3125000000000000000000000000000000000000} + -1 \cdot y}{x}}{x}}{x}}\right) \cdot \left(x + -2\right) \]

   7. Applied rewrites 94.1%

      \[\leadsto \color{blue}{\left(4.16438922228 - \frac{101.7851458539211 - \frac{3451.550173699799 + \frac{y + -124074.40615218398}{x}}{x}}{x}\right)} \cdot \left(x + -2\right) \]

   if -36 < x < 48

   1. Initial program 99.8%

      \[\frac{\left(x - 2\right) \cdot \left(\left(\left(\left(x \cdot 4.16438922228 + 78.6994924154\right) \cdot x + 137.519416416\right) \cdot x + y\right) \cdot x + z\right)}{\left(\left(\left(x + 43.3400022514\right) \cdot x + 263.505074721\right) \cdot x + 313.399215894\right) \cdot x + 47.066876606} \]
   2. Add Preprocessing
   3. Step-by-step derivation

      1. lift-/.f64 N/A

         \[\leadsto \color{blue}{\frac{\left(x - 2\right) \cdot \left(\left(\left(\left(x \cdot \frac{104109730557}{25000000000} + \frac{393497462077}{5000000000}\right) \cdot x + \frac{4297481763}{31250000}\right) \cdot x + y\right) \cdot x + z\right)}{\left(\left(\left(x + \frac{216700011257}{5000000000}\right) \cdot x + \frac{263505074721}{1000000000}\right) \cdot x + \frac{156699607947}{500000000}\right) \cdot x + \frac{23533438303}{500000000}}} \]

      2. lift-*.f64 N/A

         \[\leadsto \frac{\color{blue}{\left(x - 2\right) \cdot \left(\left(\left(\left(x \cdot \frac{104109730557}{25000000000} + \frac{393497462077}{5000000000}\right) \cdot x + \frac{4297481763}{31250000}\right) \cdot x + y\right) \cdot x + z\right)}}{\left(\left(\left(x + \frac{216700011257}{5000000000}\right) \cdot x + \frac{263505074721}{1000000000}\right) \cdot x + \frac{156699607947}{500000000}\right) \cdot x + \frac{23533438303}{500000000}} \]

      3. associate-/l* N/A

         \[\leadsto \color{blue}{\left(x - 2\right) \cdot \frac{\left(\left(\left(x \cdot \frac{104109730557}{25000000000} + \frac{393497462077}{5000000000}\right) \cdot x + \frac{4297481763}{31250000}\right) \cdot x + y\right) \cdot x + z}{\left(\left(\left(x + \frac{216700011257}{5000000000}\right) \cdot x + \frac{263505074721}{1000000000}\right) \cdot x + \frac{156699607947}{500000000}\right) \cdot x + \frac{23533438303}{500000000}}} \]

      4. *-commutative N/A

         \[\leadsto \color{blue}{\frac{\left(\left(\left(x \cdot \frac{104109730557}{25000000000} + \frac{393497462077}{5000000000}\right) \cdot x + \frac{4297481763}{31250000}\right) \cdot x + y\right) \cdot x + z}{\left(\left(\left(x + \frac{216700011257}{5000000000}\right) \cdot x + \frac{263505074721}{1000000000}\right) \cdot x + \frac{156699607947}{500000000}\right) \cdot x + \frac{23533438303}{500000000}} \cdot \left(x - 2\right)} \]

      5. lower-*.f64 N/A

         \[\leadsto \color{blue}{\frac{\left(\left(\left(x \cdot \frac{104109730557}{25000000000} + \frac{393497462077}{5000000000}\right) \cdot x + \frac{4297481763}{31250000}\right) \cdot x + y\right) \cdot x + z}{\left(\left(\left(x + \frac{216700011257}{5000000000}\right) \cdot x + \frac{263505074721}{1000000000}\right) \cdot x + \frac{156699607947}{500000000}\right) \cdot x + \frac{23533438303}{500000000}} \cdot \left(x - 2\right)} \]

   4. Applied rewrites 99.8%

      \[\leadsto \color{blue}{\frac{\mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, 4.16438922228, 78.6994924154\right), 137.519416416\right), y\right), z\right)}{\mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, x + 43.3400022514, 263.505074721\right), 313.399215894\right), 47.066876606\right)} \cdot \left(x + -2\right)} \]

   5. Taylor expanded in x around 0

      \[\leadsto \frac{\mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, \frac{104109730557}{25000000000}, \frac{393497462077}{5000000000}\right), \frac{4297481763}{31250000}\right), y\right), z\right)}{\mathsf{fma}\left(x, \color{blue}{\frac{156699607947}{500000000}}, \frac{23533438303}{500000000}\right)} \cdot \left(x + -2\right) \]
   6. Step-by-step derivation

   7. Applied rewrites 98.9%

      \[\leadsto \frac{\mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, 4.16438922228, 78.6994924154\right), 137.519416416\right), y\right), z\right)}{\mathsf{fma}\left(x, \color{blue}{313.399215894}, 47.066876606\right)} \cdot \left(x + -2\right) \]
3. Recombined 2 regimes into one program.
4. Add Preprocessing

Alternative 5: 92.1% accurate, 1.3× speedup

Math

\[\begin{array}{l} \\ \begin{array}{l} t_0 := \frac{x + -2}{0.24013125253755718}\\ \mathbf{if}\;x \leq -1.5 \cdot 10^{+39}:\\ \;\;\;\;t\_0\\ \mathbf{elif}\;x \leq 4.1 \cdot 10^{+49}:\\ \;\;\;\;\frac{\mathsf{fma}\left(x, y, z\right)}{\mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, x + 43.3400022514, 263.505074721\right), 313.399215894\right), 47.066876606\right)} \cdot \left(x + -2\right)\\ \mathbf{else}:\\ \;\;\;\;t\_0\\ \end{array} \end{array} \]

FPCore

(FPCore (x y z)
 :precision binary64
 (let* ((t_0 (/ (+ x -2.0) 0.24013125253755718)))
   (if (<= x -1.5e+39)
     t_0
     (if (<= x 4.1e+49)
       (*
        (/
         (fma x y z)
         (fma
          x
          (fma x (fma x (+ x 43.3400022514) 263.505074721) 313.399215894)
          47.066876606))
        (+ x -2.0))
       t_0))))

C

double code(double x, double y, double z) {
	double t_0 = (x + -2.0) / 0.24013125253755718;
	double tmp;
	if (x <= -1.5e+39) {
		tmp = t_0;
	} else if (x <= 4.1e+49) {
		tmp = (fma(x, y, z) / fma(x, fma(x, fma(x, (x + 43.3400022514), 263.505074721), 313.399215894), 47.066876606)) * (x + -2.0);
	} else {
		tmp = t_0;
	}
	return tmp;
}

Julia

function code(x, y, z)
	t_0 = Float64(Float64(x + -2.0) / 0.24013125253755718)
	tmp = 0.0
	if (x <= -1.5e+39)
		tmp = t_0;
	elseif (x <= 4.1e+49)
		tmp = Float64(Float64(fma(x, y, z) / fma(x, fma(x, fma(x, Float64(x + 43.3400022514), 263.505074721), 313.399215894), 47.066876606)) * Float64(x + -2.0));
	else
		tmp = t_0;
	end
	return tmp
end

Wolfram

code[x_, y_, z_] := Block[{t$95$0 = N[(N[(x + -2.0), $MachinePrecision] / 0.24013125253755718), $MachinePrecision]}, If[LessEqual[x, -1.5e+39], t$95$0, If[LessEqual[x, 4.1e+49], N[(N[(N[(x * y + z), $MachinePrecision] / N[(x * N[(x * N[(x * N[(x + 43.3400022514), $MachinePrecision] + 263.505074721), $MachinePrecision] + 313.399215894), $MachinePrecision] + 47.066876606), $MachinePrecision]), $MachinePrecision] * N[(x + -2.0), $MachinePrecision]), $MachinePrecision], t$95$0]]]

Derivation

1. Split input into 2 regimes

2. if x < -1.5e39 or 4.1e49 < x

   1. Initial program 6.9%

      \[\frac{\left(x - 2\right) \cdot \left(\left(\left(\left(x \cdot 4.16438922228 + 78.6994924154\right) \cdot x + 137.519416416\right) \cdot x + y\right) \cdot x + z\right)}{\left(\left(\left(x + 43.3400022514\right) \cdot x + 263.505074721\right) \cdot x + 313.399215894\right) \cdot x + 47.066876606} \]
   2. Add Preprocessing

   4. Applied rewrites 13.3%

      \[\leadsto \color{blue}{\frac{x + -2}{\frac{\mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, x + 43.3400022514, 263.505074721\right), 313.399215894\right), 47.066876606\right)}{\mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, 4.16438922228, 78.6994924154\right), 137.519416416\right), y\right), z\right)}}} \]

   5. Taylor expanded in x around inf

      \[\leadsto \frac{x + -2}{\color{blue}{\frac{25000000000}{104109730557}}} \]
   6. Step-by-step derivation

   7. Applied rewrites 97.0%

      \[\leadsto \frac{x + -2}{\color{blue}{0.24013125253755718}} \]

   if -1.5e39 < x < 4.1e49

   1. Initial program 97.3%

      \[\frac{\left(x - 2\right) \cdot \left(\left(\left(\left(x \cdot 4.16438922228 + 78.6994924154\right) \cdot x + 137.519416416\right) \cdot x + y\right) \cdot x + z\right)}{\left(\left(\left(x + 43.3400022514\right) \cdot x + 263.505074721\right) \cdot x + 313.399215894\right) \cdot x + 47.066876606} \]
   2. Add Preprocessing
   3. Step-by-step derivation

      1. lift-/.f64 N/A

         \[\leadsto \color{blue}{\frac{\left(x - 2\right) \cdot \left(\left(\left(\left(x \cdot \frac{104109730557}{25000000000} + \frac{393497462077}{5000000000}\right) \cdot x + \frac{4297481763}{31250000}\right) \cdot x + y\right) \cdot x + z\right)}{\left(\left(\left(x + \frac{216700011257}{5000000000}\right) \cdot x + \frac{263505074721}{1000000000}\right) \cdot x + \frac{156699607947}{500000000}\right) \cdot x + \frac{23533438303}{500000000}}} \]

      2. lift-*.f64 N/A

         \[\leadsto \frac{\color{blue}{\left(x - 2\right) \cdot \left(\left(\left(\left(x \cdot \frac{104109730557}{25000000000} + \frac{393497462077}{5000000000}\right) \cdot x + \frac{4297481763}{31250000}\right) \cdot x + y\right) \cdot x + z\right)}}{\left(\left(\left(x + \frac{216700011257}{5000000000}\right) \cdot x + \frac{263505074721}{1000000000}\right) \cdot x + \frac{156699607947}{500000000}\right) \cdot x + \frac{23533438303}{500000000}} \]

      3. associate-/l* N/A

         \[\leadsto \color{blue}{\left(x - 2\right) \cdot \frac{\left(\left(\left(x \cdot \frac{104109730557}{25000000000} + \frac{393497462077}{5000000000}\right) \cdot x + \frac{4297481763}{31250000}\right) \cdot x + y\right) \cdot x + z}{\left(\left(\left(x + \frac{216700011257}{5000000000}\right) \cdot x + \frac{263505074721}{1000000000}\right) \cdot x + \frac{156699607947}{500000000}\right) \cdot x + \frac{23533438303}{500000000}}} \]

      4. *-commutative N/A

         \[\leadsto \color{blue}{\frac{\left(\left(\left(x \cdot \frac{104109730557}{25000000000} + \frac{393497462077}{5000000000}\right) \cdot x + \frac{4297481763}{31250000}\right) \cdot x + y\right) \cdot x + z}{\left(\left(\left(x + \frac{216700011257}{5000000000}\right) \cdot x + \frac{263505074721}{1000000000}\right) \cdot x + \frac{156699607947}{500000000}\right) \cdot x + \frac{23533438303}{500000000}} \cdot \left(x - 2\right)} \]

      5. lower-*.f64 N/A

         \[\leadsto \color{blue}{\frac{\left(\left(\left(x \cdot \frac{104109730557}{25000000000} + \frac{393497462077}{5000000000}\right) \cdot x + \frac{4297481763}{31250000}\right) \cdot x + y\right) \cdot x + z}{\left(\left(\left(x + \frac{216700011257}{5000000000}\right) \cdot x + \frac{263505074721}{1000000000}\right) \cdot x + \frac{156699607947}{500000000}\right) \cdot x + \frac{23533438303}{500000000}} \cdot \left(x - 2\right)} \]

   4. Applied rewrites 99.1%

      \[\leadsto \color{blue}{\frac{\mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, 4.16438922228, 78.6994924154\right), 137.519416416\right), y\right), z\right)}{\mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, x + 43.3400022514, 263.505074721\right), 313.399215894\right), 47.066876606\right)} \cdot \left(x + -2\right)} \]

   5. Taylor expanded in x around 0

      \[\leadsto \frac{\color{blue}{z + x \cdot y}}{\mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, x + \frac{216700011257}{5000000000}, \frac{263505074721}{1000000000}\right), \frac{156699607947}{500000000}\right), \frac{23533438303}{500000000}\right)} \cdot \left(x + -2\right) \]
   6. Step-by-step derivation

      1. +-commutative N/A

         \[\leadsto \frac{\color{blue}{x \cdot y + z}}{\mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, x + \frac{216700011257}{5000000000}, \frac{263505074721}{1000000000}\right), \frac{156699607947}{500000000}\right), \frac{23533438303}{500000000}\right)} \cdot \left(x + -2\right) \]

      2. lower-fma.f64 94.4

         \[\leadsto \frac{\color{blue}{\mathsf{fma}\left(x, y, z\right)}}{\mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, x + 43.3400022514, 263.505074721\right), 313.399215894\right), 47.066876606\right)} \cdot \left(x + -2\right) \]

   7. Applied rewrites 94.4%

      \[\leadsto \frac{\color{blue}{\mathsf{fma}\left(x, y, z\right)}}{\mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, x + 43.3400022514, 263.505074721\right), 313.399215894\right), 47.066876606\right)} \cdot \left(x + -2\right) \]
3. Recombined 2 regimes into one program.
4. Add Preprocessing

Alternative 6: 91.2% accurate, 1.6× speedup

Math

\[\begin{array}{l} \\ \begin{array}{l} t_0 := \frac{x + -2}{0.24013125253755718}\\ \mathbf{if}\;x \leq -1.8 \cdot 10^{+31}:\\ \;\;\;\;t\_0\\ \mathbf{elif}\;x \leq 1.3 \cdot 10^{+37}:\\ \;\;\;\;\mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, 4.16438922228, 78.6994924154\right), 137.519416416\right), y\right), z\right) \cdot \mathsf{fma}\left(x, 0.3041881842569256, -0.0424927283095952\right)\\ \mathbf{else}:\\ \;\;\;\;t\_0\\ \end{array} \end{array} \]

FPCore

(FPCore (x y z)
 :precision binary64
 (let* ((t_0 (/ (+ x -2.0) 0.24013125253755718)))
   (if (<= x -1.8e+31)
     t_0
     (if (<= x 1.3e+37)
       (*
        (fma
         x
         (fma x (fma x (fma x 4.16438922228 78.6994924154) 137.519416416) y)
         z)
        (fma x 0.3041881842569256 -0.0424927283095952))
       t_0))))

C

double code(double x, double y, double z) {
	double t_0 = (x + -2.0) / 0.24013125253755718;
	double tmp;
	if (x <= -1.8e+31) {
		tmp = t_0;
	} else if (x <= 1.3e+37) {
		tmp = fma(x, fma(x, fma(x, fma(x, 4.16438922228, 78.6994924154), 137.519416416), y), z) * fma(x, 0.3041881842569256, -0.0424927283095952);
	} else {
		tmp = t_0;
	}
	return tmp;
}

Julia

function code(x, y, z)
	t_0 = Float64(Float64(x + -2.0) / 0.24013125253755718)
	tmp = 0.0
	if (x <= -1.8e+31)
		tmp = t_0;
	elseif (x <= 1.3e+37)
		tmp = Float64(fma(x, fma(x, fma(x, fma(x, 4.16438922228, 78.6994924154), 137.519416416), y), z) * fma(x, 0.3041881842569256, -0.0424927283095952));
	else
		tmp = t_0;
	end
	return tmp
end

Wolfram

code[x_, y_, z_] := Block[{t$95$0 = N[(N[(x + -2.0), $MachinePrecision] / 0.24013125253755718), $MachinePrecision]}, If[LessEqual[x, -1.8e+31], t$95$0, If[LessEqual[x, 1.3e+37], N[(N[(x * N[(x * N[(x * N[(x * 4.16438922228 + 78.6994924154), $MachinePrecision] + 137.519416416), $MachinePrecision] + y), $MachinePrecision] + z), $MachinePrecision] * N[(x * 0.3041881842569256 + -0.0424927283095952), $MachinePrecision]), $MachinePrecision], t$95$0]]]

Derivation

1. Split input into 2 regimes

2. if x < -1.79999999999999998e31 or 1.3e37 < x

   1. Initial program 8.5%

      \[\frac{\left(x - 2\right) \cdot \left(\left(\left(\left(x \cdot 4.16438922228 + 78.6994924154\right) \cdot x + 137.519416416\right) \cdot x + y\right) \cdot x + z\right)}{\left(\left(\left(x + 43.3400022514\right) \cdot x + 263.505074721\right) \cdot x + 313.399215894\right) \cdot x + 47.066876606} \]
   2. Add Preprocessing

   4. Applied rewrites 16.4%

      \[\leadsto \color{blue}{\frac{x + -2}{\frac{\mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, x + 43.3400022514, 263.505074721\right), 313.399215894\right), 47.066876606\right)}{\mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, 4.16438922228, 78.6994924154\right), 137.519416416\right), y\right), z\right)}}} \]

   5. Taylor expanded in x around inf

      \[\leadsto \frac{x + -2}{\color{blue}{\frac{25000000000}{104109730557}}} \]
   6. Step-by-step derivation

   7. Applied rewrites 94.4%

      \[\leadsto \frac{x + -2}{\color{blue}{0.24013125253755718}} \]

   if -1.79999999999999998e31 < x < 1.3e37

   1. Initial program 98.5%

      \[\frac{\left(x - 2\right) \cdot \left(\left(\left(\left(x \cdot 4.16438922228 + 78.6994924154\right) \cdot x + 137.519416416\right) \cdot x + y\right) \cdot x + z\right)}{\left(\left(\left(x + 43.3400022514\right) \cdot x + 263.505074721\right) \cdot x + 313.399215894\right) \cdot x + 47.066876606} \]
   2. Add Preprocessing
   3. Step-by-step derivation

      1. lift-/.f64 N/A

         \[\leadsto \color{blue}{\frac{\left(x - 2\right) \cdot \left(\left(\left(\left(x \cdot \frac{104109730557}{25000000000} + \frac{393497462077}{5000000000}\right) \cdot x + \frac{4297481763}{31250000}\right) \cdot x + y\right) \cdot x + z\right)}{\left(\left(\left(x + \frac{216700011257}{5000000000}\right) \cdot x + \frac{263505074721}{1000000000}\right) \cdot x + \frac{156699607947}{500000000}\right) \cdot x + \frac{23533438303}{500000000}}} \]

      2. lift-*.f64 N/A

         \[\leadsto \frac{\color{blue}{\left(x - 2\right) \cdot \left(\left(\left(\left(x \cdot \frac{104109730557}{25000000000} + \frac{393497462077}{5000000000}\right) \cdot x + \frac{4297481763}{31250000}\right) \cdot x + y\right) \cdot x + z\right)}}{\left(\left(\left(x + \frac{216700011257}{5000000000}\right) \cdot x + \frac{263505074721}{1000000000}\right) \cdot x + \frac{156699607947}{500000000}\right) \cdot x + \frac{23533438303}{500000000}} \]

      3. *-commutative N/A

         \[\leadsto \frac{\color{blue}{\left(\left(\left(\left(x \cdot \frac{104109730557}{25000000000} + \frac{393497462077}{5000000000}\right) \cdot x + \frac{4297481763}{31250000}\right) \cdot x + y\right) \cdot x + z\right) \cdot \left(x - 2\right)}}{\left(\left(\left(x + \frac{216700011257}{5000000000}\right) \cdot x + \frac{263505074721}{1000000000}\right) \cdot x + \frac{156699607947}{500000000}\right) \cdot x + \frac{23533438303}{500000000}} \]

      4. associate-/l* N/A

         \[\leadsto \color{blue}{\left(\left(\left(\left(x \cdot \frac{104109730557}{25000000000} + \frac{393497462077}{5000000000}\right) \cdot x + \frac{4297481763}{31250000}\right) \cdot x + y\right) \cdot x + z\right) \cdot \frac{x - 2}{\left(\left(\left(x + \frac{216700011257}{5000000000}\right) \cdot x + \frac{263505074721}{1000000000}\right) \cdot x + \frac{156699607947}{500000000}\right) \cdot x + \frac{23533438303}{500000000}}} \]

      5. lower-*.f64 N/A

         \[\leadsto \color{blue}{\left(\left(\left(\left(x \cdot \frac{104109730557}{25000000000} + \frac{393497462077}{5000000000}\right) \cdot x + \frac{4297481763}{31250000}\right) \cdot x + y\right) \cdot x + z\right) \cdot \frac{x - 2}{\left(\left(\left(x + \frac{216700011257}{5000000000}\right) \cdot x + \frac{263505074721}{1000000000}\right) \cdot x + \frac{156699607947}{500000000}\right) \cdot x + \frac{23533438303}{500000000}}} \]

   4. Applied rewrites 98.7%

      \[\leadsto \color{blue}{\mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, 4.16438922228, 78.6994924154\right), 137.519416416\right), y\right), z\right) \cdot \frac{x + -2}{\mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, x + 43.3400022514, 263.505074721\right), 313.399215894\right), 47.066876606\right)}} \]

   5. Taylor expanded in x around 0

      \[\leadsto \mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, \frac{104109730557}{25000000000}, \frac{393497462077}{5000000000}\right), \frac{4297481763}{31250000}\right), y\right), z\right) \cdot \color{blue}{\left(\frac{168466327098500000000}{553822718361107519809} \cdot x - \frac{1000000000}{23533438303}\right)} \]
   6. Step-by-step derivation

      1. sub-neg N/A

         \[\leadsto \mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, \frac{104109730557}{25000000000}, \frac{393497462077}{5000000000}\right), \frac{4297481763}{31250000}\right), y\right), z\right) \cdot \color{blue}{\left(\frac{168466327098500000000}{553822718361107519809} \cdot x + \left(\mathsf{neg}\left(\frac{1000000000}{23533438303}\right)\right)\right)} \]

      2. *-commutative N/A

         \[\leadsto \mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, \frac{104109730557}{25000000000}, \frac{393497462077}{5000000000}\right), \frac{4297481763}{31250000}\right), y\right), z\right) \cdot \left(\color{blue}{x \cdot \frac{168466327098500000000}{553822718361107519809}} + \left(\mathsf{neg}\left(\frac{1000000000}{23533438303}\right)\right)\right) \]

      3. metadata-eval N/A

         \[\leadsto \mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, \frac{104109730557}{25000000000}, \frac{393497462077}{5000000000}\right), \frac{4297481763}{31250000}\right), y\right), z\right) \cdot \left(x \cdot \frac{168466327098500000000}{553822718361107519809} + \color{blue}{\frac{-1000000000}{23533438303}}\right) \]

      4. lower-fma.f64 93.4

         \[\leadsto \mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, 4.16438922228, 78.6994924154\right), 137.519416416\right), y\right), z\right) \cdot \color{blue}{\mathsf{fma}\left(x, 0.3041881842569256, -0.0424927283095952\right)} \]

   7. Applied rewrites 93.4%

      \[\leadsto \mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, 4.16438922228, 78.6994924154\right), 137.519416416\right), y\right), z\right) \cdot \color{blue}{\mathsf{fma}\left(x, 0.3041881842569256, -0.0424927283095952\right)} \]
3. Recombined 2 regimes into one program.
4. Add Preprocessing

Alternative 7: 91.4% accurate, 1.9× speedup

Math

\[\begin{array}{l} \\ \begin{array}{l} t_0 := \frac{x + -2}{0.24013125253755718}\\ \mathbf{if}\;x \leq -1.8 \cdot 10^{+31}:\\ \;\;\;\;t\_0\\ \mathbf{elif}\;x \leq 2:\\ \;\;\;\;\mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, 4.16438922228, 78.6994924154\right), 137.519416416\right), y\right), z\right) \cdot -0.0424927283095952\\ \mathbf{else}:\\ \;\;\;\;t\_0\\ \end{array} \end{array} \]

FPCore

(FPCore (x y z)
 :precision binary64
 (let* ((t_0 (/ (+ x -2.0) 0.24013125253755718)))
   (if (<= x -1.8e+31)
     t_0
     (if (<= x 2.0)
       (*
        (fma
         x
         (fma x (fma x (fma x 4.16438922228 78.6994924154) 137.519416416) y)
         z)
        -0.0424927283095952)
       t_0))))

C

double code(double x, double y, double z) {
	double t_0 = (x + -2.0) / 0.24013125253755718;
	double tmp;
	if (x <= -1.8e+31) {
		tmp = t_0;
	} else if (x <= 2.0) {
		tmp = fma(x, fma(x, fma(x, fma(x, 4.16438922228, 78.6994924154), 137.519416416), y), z) * -0.0424927283095952;
	} else {
		tmp = t_0;
	}
	return tmp;
}

Julia

function code(x, y, z)
	t_0 = Float64(Float64(x + -2.0) / 0.24013125253755718)
	tmp = 0.0
	if (x <= -1.8e+31)
		tmp = t_0;
	elseif (x <= 2.0)
		tmp = Float64(fma(x, fma(x, fma(x, fma(x, 4.16438922228, 78.6994924154), 137.519416416), y), z) * -0.0424927283095952);
	else
		tmp = t_0;
	end
	return tmp
end

Wolfram

code[x_, y_, z_] := Block[{t$95$0 = N[(N[(x + -2.0), $MachinePrecision] / 0.24013125253755718), $MachinePrecision]}, If[LessEqual[x, -1.8e+31], t$95$0, If[LessEqual[x, 2.0], N[(N[(x * N[(x * N[(x * N[(x * 4.16438922228 + 78.6994924154), $MachinePrecision] + 137.519416416), $MachinePrecision] + y), $MachinePrecision] + z), $MachinePrecision] * -0.0424927283095952), $MachinePrecision], t$95$0]]]

Derivation

1. Split input into 2 regimes

2. if x < -1.79999999999999998e31 or 2 < x

   1. Initial program 11.8%

      \[\frac{\left(x - 2\right) \cdot \left(\left(\left(\left(x \cdot 4.16438922228 + 78.6994924154\right) \cdot x + 137.519416416\right) \cdot x + y\right) \cdot x + z\right)}{\left(\left(\left(x + 43.3400022514\right) \cdot x + 263.505074721\right) \cdot x + 313.399215894\right) \cdot x + 47.066876606} \]
   2. Add Preprocessing

   4. Applied rewrites 19.4%

      \[\leadsto \color{blue}{\frac{x + -2}{\frac{\mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, x + 43.3400022514, 263.505074721\right), 313.399215894\right), 47.066876606\right)}{\mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, 4.16438922228, 78.6994924154\right), 137.519416416\right), y\right), z\right)}}} \]

   5. Taylor expanded in x around inf

      \[\leadsto \frac{x + -2}{\color{blue}{\frac{25000000000}{104109730557}}} \]
   6. Step-by-step derivation

   7. Applied rewrites 91.1%

      \[\leadsto \frac{x + -2}{\color{blue}{0.24013125253755718}} \]

   if -1.79999999999999998e31 < x < 2

   1. Initial program 98.5%

      \[\frac{\left(x - 2\right) \cdot \left(\left(\left(\left(x \cdot 4.16438922228 + 78.6994924154\right) \cdot x + 137.519416416\right) \cdot x + y\right) \cdot x + z\right)}{\left(\left(\left(x + 43.3400022514\right) \cdot x + 263.505074721\right) \cdot x + 313.399215894\right) \cdot x + 47.066876606} \]
   2. Add Preprocessing
   3. Step-by-step derivation

      1. lift-/.f64 N/A

         \[\leadsto \color{blue}{\frac{\left(x - 2\right) \cdot \left(\left(\left(\left(x \cdot \frac{104109730557}{25000000000} + \frac{393497462077}{5000000000}\right) \cdot x + \frac{4297481763}{31250000}\right) \cdot x + y\right) \cdot x + z\right)}{\left(\left(\left(x + \frac{216700011257}{5000000000}\right) \cdot x + \frac{263505074721}{1000000000}\right) \cdot x + \frac{156699607947}{500000000}\right) \cdot x + \frac{23533438303}{500000000}}} \]

      2. lift-*.f64 N/A

         \[\leadsto \frac{\color{blue}{\left(x - 2\right) \cdot \left(\left(\left(\left(x \cdot \frac{104109730557}{25000000000} + \frac{393497462077}{5000000000}\right) \cdot x + \frac{4297481763}{31250000}\right) \cdot x + y\right) \cdot x + z\right)}}{\left(\left(\left(x + \frac{216700011257}{5000000000}\right) \cdot x + \frac{263505074721}{1000000000}\right) \cdot x + \frac{156699607947}{500000000}\right) \cdot x + \frac{23533438303}{500000000}} \]

      3. *-commutative N/A

         \[\leadsto \frac{\color{blue}{\left(\left(\left(\left(x \cdot \frac{104109730557}{25000000000} + \frac{393497462077}{5000000000}\right) \cdot x + \frac{4297481763}{31250000}\right) \cdot x + y\right) \cdot x + z\right) \cdot \left(x - 2\right)}}{\left(\left(\left(x + \frac{216700011257}{5000000000}\right) \cdot x + \frac{263505074721}{1000000000}\right) \cdot x + \frac{156699607947}{500000000}\right) \cdot x + \frac{23533438303}{500000000}} \]

      4. associate-/l* N/A

         \[\leadsto \color{blue}{\left(\left(\left(\left(x \cdot \frac{104109730557}{25000000000} + \frac{393497462077}{5000000000}\right) \cdot x + \frac{4297481763}{31250000}\right) \cdot x + y\right) \cdot x + z\right) \cdot \frac{x - 2}{\left(\left(\left(x + \frac{216700011257}{5000000000}\right) \cdot x + \frac{263505074721}{1000000000}\right) \cdot x + \frac{156699607947}{500000000}\right) \cdot x + \frac{23533438303}{500000000}}} \]

      5. lower-*.f64 N/A

         \[\leadsto \color{blue}{\left(\left(\left(\left(x \cdot \frac{104109730557}{25000000000} + \frac{393497462077}{5000000000}\right) \cdot x + \frac{4297481763}{31250000}\right) \cdot x + y\right) \cdot x + z\right) \cdot \frac{x - 2}{\left(\left(\left(x + \frac{216700011257}{5000000000}\right) \cdot x + \frac{263505074721}{1000000000}\right) \cdot x + \frac{156699607947}{500000000}\right) \cdot x + \frac{23533438303}{500000000}}} \]

   4. Applied rewrites 98.7%

      \[\leadsto \color{blue}{\mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, 4.16438922228, 78.6994924154\right), 137.519416416\right), y\right), z\right) \cdot \frac{x + -2}{\mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, x + 43.3400022514, 263.505074721\right), 313.399215894\right), 47.066876606\right)}} \]

   5. Taylor expanded in x around 0

      \[\leadsto \mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, \frac{104109730557}{25000000000}, \frac{393497462077}{5000000000}\right), \frac{4297481763}{31250000}\right), y\right), z\right) \cdot \color{blue}{\frac{-1000000000}{23533438303}} \]
   6. Step-by-step derivation

   7. Applied rewrites 95.0%

      \[\leadsto \mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, 4.16438922228, 78.6994924154\right), 137.519416416\right), y\right), z\right) \cdot \color{blue}{-0.0424927283095952} \]
3. Recombined 2 regimes into one program.
4. Add Preprocessing

Alternative 8: 89.3% accurate, 2.3× speedup

Math

\[\begin{array}{l} \\ \begin{array}{l} t_0 := \frac{x + -2}{0.24013125253755718}\\ \mathbf{if}\;x \leq -7.6 \cdot 10^{+28}:\\ \;\;\;\;t\_0\\ \mathbf{elif}\;x \leq 1.5:\\ \;\;\;\;\mathsf{fma}\left(x, \mathsf{fma}\left(-0.0424927283095952, y, z \cdot 0.3041881842569256\right), -0.0424927283095952 \cdot z\right)\\ \mathbf{else}:\\ \;\;\;\;t\_0\\ \end{array} \end{array} \]

FPCore

(FPCore (x y z)
 :precision binary64
 (let* ((t_0 (/ (+ x -2.0) 0.24013125253755718)))
   (if (<= x -7.6e+28)
     t_0
     (if (<= x 1.5)
       (fma
        x
        (fma -0.0424927283095952 y (* z 0.3041881842569256))
        (* -0.0424927283095952 z))
       t_0))))

C

double code(double x, double y, double z) {
	double t_0 = (x + -2.0) / 0.24013125253755718;
	double tmp;
	if (x <= -7.6e+28) {
		tmp = t_0;
	} else if (x <= 1.5) {
		tmp = fma(x, fma(-0.0424927283095952, y, (z * 0.3041881842569256)), (-0.0424927283095952 * z));
	} else {
		tmp = t_0;
	}
	return tmp;
}

Julia

function code(x, y, z)
	t_0 = Float64(Float64(x + -2.0) / 0.24013125253755718)
	tmp = 0.0
	if (x <= -7.6e+28)
		tmp = t_0;
	elseif (x <= 1.5)
		tmp = fma(x, fma(-0.0424927283095952, y, Float64(z * 0.3041881842569256)), Float64(-0.0424927283095952 * z));
	else
		tmp = t_0;
	end
	return tmp
end

Wolfram

code[x_, y_, z_] := Block[{t$95$0 = N[(N[(x + -2.0), $MachinePrecision] / 0.24013125253755718), $MachinePrecision]}, If[LessEqual[x, -7.6e+28], t$95$0, If[LessEqual[x, 1.5], N[(x * N[(-0.0424927283095952 * y + N[(z * 0.3041881842569256), $MachinePrecision]), $MachinePrecision] + N[(-0.0424927283095952 * z), $MachinePrecision]), $MachinePrecision], t$95$0]]]

Derivation

1. Split input into 2 regimes

2. if x < -7.5999999999999998e28 or 1.5 < x

   1. Initial program 11.7%

      \[\frac{\left(x - 2\right) \cdot \left(\left(\left(\left(x \cdot 4.16438922228 + 78.6994924154\right) \cdot x + 137.519416416\right) \cdot x + y\right) \cdot x + z\right)}{\left(\left(\left(x + 43.3400022514\right) \cdot x + 263.505074721\right) \cdot x + 313.399215894\right) \cdot x + 47.066876606} \]
   2. Add Preprocessing

   4. Applied rewrites 20.1%

      \[\leadsto \color{blue}{\frac{x + -2}{\frac{\mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, x + 43.3400022514, 263.505074721\right), 313.399215894\right), 47.066876606\right)}{\mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, 4.16438922228, 78.6994924154\right), 137.519416416\right), y\right), z\right)}}} \]

   5. Taylor expanded in x around inf

      \[\leadsto \frac{x + -2}{\color{blue}{\frac{25000000000}{104109730557}}} \]
   6. Step-by-step derivation

   7. Applied rewrites 90.3%

      \[\leadsto \frac{x + -2}{\color{blue}{0.24013125253755718}} \]

   if -7.5999999999999998e28 < x < 1.5

   1. Initial program 99.1%

      \[\frac{\left(x - 2\right) \cdot \left(\left(\left(\left(x \cdot 4.16438922228 + 78.6994924154\right) \cdot x + 137.519416416\right) \cdot x + y\right) \cdot x + z\right)}{\left(\left(\left(x + 43.3400022514\right) \cdot x + 263.505074721\right) \cdot x + 313.399215894\right) \cdot x + 47.066876606} \]
   2. Add Preprocessing

   4. Applied rewrites 99.1%

      \[\leadsto \color{blue}{\frac{\mathsf{fma}\left(x, x, -4\right) \cdot \frac{\mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, 4.16438922228, 78.6994924154\right), 137.519416416\right), y\right), z\right)}{\mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, x + 43.3400022514, 263.505074721\right), 313.399215894\right), 47.066876606\right)}}{x + 2}} \]

   5. Taylor expanded in x around 0

      \[\leadsto \color{blue}{\frac{-1000000000}{23533438303} \cdot z + x \cdot \left(\frac{-1000000000}{23533438303} \cdot y - \frac{-168466327098500000000}{553822718361107519809} \cdot z\right)} \]
   6. Step-by-step derivation

      1. +-commutative N/A

         \[\leadsto \color{blue}{x \cdot \left(\frac{-1000000000}{23533438303} \cdot y - \frac{-168466327098500000000}{553822718361107519809} \cdot z\right) + \frac{-1000000000}{23533438303} \cdot z} \]

      2. lower-fma.f64 N/A

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

      3. sub-neg N/A

         \[\leadsto \mathsf{fma}\left(x, \color{blue}{\frac{-1000000000}{23533438303} \cdot y + \left(\mathsf{neg}\left(\frac{-168466327098500000000}{553822718361107519809} \cdot z\right)\right)}, \frac{-1000000000}{23533438303} \cdot z\right) \]

      4. lower-fma.f64 N/A

         \[\leadsto \mathsf{fma}\left(x, \color{blue}{\mathsf{fma}\left(\frac{-1000000000}{23533438303}, y, \mathsf{neg}\left(\frac{-168466327098500000000}{553822718361107519809} \cdot z\right)\right)}, \frac{-1000000000}{23533438303} \cdot z\right) \]

      5. *-commutative N/A

         \[\leadsto \mathsf{fma}\left(x, \mathsf{fma}\left(\frac{-1000000000}{23533438303}, y, \mathsf{neg}\left(\color{blue}{z \cdot \frac{-168466327098500000000}{553822718361107519809}}\right)\right), \frac{-1000000000}{23533438303} \cdot z\right) \]

      6. distribute-rgt-neg-in N/A

         \[\leadsto \mathsf{fma}\left(x, \mathsf{fma}\left(\frac{-1000000000}{23533438303}, y, \color{blue}{z \cdot \left(\mathsf{neg}\left(\frac{-168466327098500000000}{553822718361107519809}\right)\right)}\right), \frac{-1000000000}{23533438303} \cdot z\right) \]

      7. lower-*.f64 N/A

         \[\leadsto \mathsf{fma}\left(x, \mathsf{fma}\left(\frac{-1000000000}{23533438303}, y, \color{blue}{z \cdot \left(\mathsf{neg}\left(\frac{-168466327098500000000}{553822718361107519809}\right)\right)}\right), \frac{-1000000000}{23533438303} \cdot z\right) \]

      8. metadata-eval N/A

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

      9. lower-*.f64 91.7

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

   7. Applied rewrites 91.7%

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

Alternative 9: 78.9% accurate, 2.7× speedup

Math

\[\begin{array}{l} \\ \begin{array}{l} t_0 := \frac{x + -2}{0.24013125253755718}\\ \mathbf{if}\;x \leq -6 \cdot 10^{+28}:\\ \;\;\;\;t\_0\\ \mathbf{elif}\;x \leq 17500000000:\\ \;\;\;\;\mathsf{fma}\left(x, x \cdot -5.843575199059173, -0.0424927283095952 \cdot z\right)\\ \mathbf{else}:\\ \;\;\;\;t\_0\\ \end{array} \end{array} \]

FPCore

(FPCore (x y z)
 :precision binary64
 (let* ((t_0 (/ (+ x -2.0) 0.24013125253755718)))
   (if (<= x -6e+28)
     t_0
     (if (<= x 17500000000.0)
       (fma x (* x -5.843575199059173) (* -0.0424927283095952 z))
       t_0))))

C

double code(double x, double y, double z) {
	double t_0 = (x + -2.0) / 0.24013125253755718;
	double tmp;
	if (x <= -6e+28) {
		tmp = t_0;
	} else if (x <= 17500000000.0) {
		tmp = fma(x, (x * -5.843575199059173), (-0.0424927283095952 * z));
	} else {
		tmp = t_0;
	}
	return tmp;
}

Julia

function code(x, y, z)
	t_0 = Float64(Float64(x + -2.0) / 0.24013125253755718)
	tmp = 0.0
	if (x <= -6e+28)
		tmp = t_0;
	elseif (x <= 17500000000.0)
		tmp = fma(x, Float64(x * -5.843575199059173), Float64(-0.0424927283095952 * z));
	else
		tmp = t_0;
	end
	return tmp
end

Wolfram

code[x_, y_, z_] := Block[{t$95$0 = N[(N[(x + -2.0), $MachinePrecision] / 0.24013125253755718), $MachinePrecision]}, If[LessEqual[x, -6e+28], t$95$0, If[LessEqual[x, 17500000000.0], N[(x * N[(x * -5.843575199059173), $MachinePrecision] + N[(-0.0424927283095952 * z), $MachinePrecision]), $MachinePrecision], t$95$0]]]

Derivation

1. Split input into 2 regimes

2. if x < -6.0000000000000002e28 or 1.75e10 < x

   1. Initial program 10.9%

      \[\frac{\left(x - 2\right) \cdot \left(\left(\left(\left(x \cdot 4.16438922228 + 78.6994924154\right) \cdot x + 137.519416416\right) \cdot x + y\right) \cdot x + z\right)}{\left(\left(\left(x + 43.3400022514\right) \cdot x + 263.505074721\right) \cdot x + 313.399215894\right) \cdot x + 47.066876606} \]
   2. Add Preprocessing

   4. Applied rewrites 19.3%

      \[\leadsto \color{blue}{\frac{x + -2}{\frac{\mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, x + 43.3400022514, 263.505074721\right), 313.399215894\right), 47.066876606\right)}{\mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, 4.16438922228, 78.6994924154\right), 137.519416416\right), y\right), z\right)}}} \]

   5. Taylor expanded in x around inf

      \[\leadsto \frac{x + -2}{\color{blue}{\frac{25000000000}{104109730557}}} \]
   6. Step-by-step derivation

   7. Applied rewrites 90.3%

      \[\leadsto \frac{x + -2}{\color{blue}{0.24013125253755718}} \]

   if -6.0000000000000002e28 < x < 1.75e10

   1. Initial program 99.8%

      \[\frac{\left(x - 2\right) \cdot \left(\left(\left(\left(x \cdot 4.16438922228 + 78.6994924154\right) \cdot x + 137.519416416\right) \cdot x + y\right) \cdot x + z\right)}{\left(\left(\left(x + 43.3400022514\right) \cdot x + 263.505074721\right) \cdot x + 313.399215894\right) \cdot x + 47.066876606} \]
   2. Add Preprocessing

   3. Taylor expanded in y around 0

      \[\leadsto \color{blue}{\frac{\left(z + {x}^{2} \cdot \left(\frac{4297481763}{31250000} + x \cdot \left(\frac{393497462077}{5000000000} + \frac{104109730557}{25000000000} \cdot x\right)\right)\right) \cdot \left(x - 2\right)}{\frac{23533438303}{500000000} + x \cdot \left(\frac{156699607947}{500000000} + x \cdot \left(\frac{263505074721}{1000000000} + x \cdot \left(\frac{216700011257}{5000000000} + x\right)\right)\right)}} \]
   4. Step-by-step derivation

      1. lower-/.f64 N/A

         \[\leadsto \color{blue}{\frac{\left(z + {x}^{2} \cdot \left(\frac{4297481763}{31250000} + x \cdot \left(\frac{393497462077}{5000000000} + \frac{104109730557}{25000000000} \cdot x\right)\right)\right) \cdot \left(x - 2\right)}{\frac{23533438303}{500000000} + x \cdot \left(\frac{156699607947}{500000000} + x \cdot \left(\frac{263505074721}{1000000000} + x \cdot \left(\frac{216700011257}{5000000000} + x\right)\right)\right)}} \]

   5. Applied rewrites 71.7%

      \[\leadsto \color{blue}{\frac{\mathsf{fma}\left(x \cdot x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, 4.16438922228, 78.6994924154\right), 137.519416416\right), z\right) \cdot \left(x + -2\right)}{\mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, x + 43.3400022514, 263.505074721\right), 313.399215894\right), 47.066876606\right)}} \]

   6. Taylor expanded in x around 0

      \[\leadsto \frac{-1000000000}{23533438303} \cdot z + \color{blue}{x \cdot \left(\left(-1 \cdot \left(x \cdot \left(\frac{137519416416}{23533438303} + \left(\frac{-131752537360500000000}{553822718361107519809} \cdot z + \frac{156699607947}{23533438303} \cdot \left(\frac{500000000}{23533438303} \cdot z - \frac{-156699607947000000000}{553822718361107519809} \cdot z\right)\right)\right)\right) + \frac{500000000}{23533438303} \cdot z\right) - \frac{-156699607947000000000}{553822718361107519809} \cdot z\right)} \]
   7. Step-by-step derivation

   8. Applied rewrites 69.5%

      \[\leadsto \mathsf{fma}\left(x, \color{blue}{\mathsf{fma}\left(-x, \mathsf{fma}\left(z, -0.23789659216289816, \mathsf{fma}\left(z, 2.025465578019411, 5.843575199059173\right)\right), z \cdot 0.3041881842569256\right)}, -0.0424927283095952 \cdot z\right) \]

   9. Taylor expanded in z around 0

      \[\leadsto \mathsf{fma}\left(x, \frac{-137519416416}{23533438303} \cdot x, \frac{-1000000000}{23533438303} \cdot z\right) \]
   10. Step-by-step derivation

   11. Applied rewrites 69.0%

       \[\leadsto \mathsf{fma}\left(x, x \cdot -5.843575199059173, -0.0424927283095952 \cdot z\right) \]
3. Recombined 2 regimes into one program.
4. Add Preprocessing

Alternative 10: 76.1% accurate, 2.7× speedup

Math

\[\begin{array}{l} \\ \begin{array}{l} t_0 := \frac{x + -2}{0.24013125253755718}\\ \mathbf{if}\;x \leq -6 \cdot 10^{+28}:\\ \;\;\;\;t\_0\\ \mathbf{elif}\;x \leq 1.3 \cdot 10^{+37}:\\ \;\;\;\;\mathsf{fma}\left(-0.0424927283095952, z, \left(z \cdot x\right) \cdot 0.3041881842569256\right)\\ \mathbf{else}:\\ \;\;\;\;t\_0\\ \end{array} \end{array} \]

FPCore

(FPCore (x y z)
 :precision binary64
 (let* ((t_0 (/ (+ x -2.0) 0.24013125253755718)))
   (if (<= x -6e+28)
     t_0
     (if (<= x 1.3e+37)
       (fma -0.0424927283095952 z (* (* z x) 0.3041881842569256))
       t_0))))

C

double code(double x, double y, double z) {
	double t_0 = (x + -2.0) / 0.24013125253755718;
	double tmp;
	if (x <= -6e+28) {
		tmp = t_0;
	} else if (x <= 1.3e+37) {
		tmp = fma(-0.0424927283095952, z, ((z * x) * 0.3041881842569256));
	} else {
		tmp = t_0;
	}
	return tmp;
}

Julia

function code(x, y, z)
	t_0 = Float64(Float64(x + -2.0) / 0.24013125253755718)
	tmp = 0.0
	if (x <= -6e+28)
		tmp = t_0;
	elseif (x <= 1.3e+37)
		tmp = fma(-0.0424927283095952, z, Float64(Float64(z * x) * 0.3041881842569256));
	else
		tmp = t_0;
	end
	return tmp
end

Wolfram

code[x_, y_, z_] := Block[{t$95$0 = N[(N[(x + -2.0), $MachinePrecision] / 0.24013125253755718), $MachinePrecision]}, If[LessEqual[x, -6e+28], t$95$0, If[LessEqual[x, 1.3e+37], N[(-0.0424927283095952 * z + N[(N[(z * x), $MachinePrecision] * 0.3041881842569256), $MachinePrecision]), $MachinePrecision], t$95$0]]]

Derivation

1. Split input into 2 regimes

2. if x < -6.0000000000000002e28 or 1.3e37 < x

   1. Initial program 8.5%

      \[\frac{\left(x - 2\right) \cdot \left(\left(\left(\left(x \cdot 4.16438922228 + 78.6994924154\right) \cdot x + 137.519416416\right) \cdot x + y\right) \cdot x + z\right)}{\left(\left(\left(x + 43.3400022514\right) \cdot x + 263.505074721\right) \cdot x + 313.399215894\right) \cdot x + 47.066876606} \]
   2. Add Preprocessing

   4. Applied rewrites 17.1%

      \[\leadsto \color{blue}{\frac{x + -2}{\frac{\mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, x + 43.3400022514, 263.505074721\right), 313.399215894\right), 47.066876606\right)}{\mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, 4.16438922228, 78.6994924154\right), 137.519416416\right), y\right), z\right)}}} \]

   5. Taylor expanded in x around inf

      \[\leadsto \frac{x + -2}{\color{blue}{\frac{25000000000}{104109730557}}} \]
   6. Step-by-step derivation

   7. Applied rewrites 92.7%

      \[\leadsto \frac{x + -2}{\color{blue}{0.24013125253755718}} \]

   if -6.0000000000000002e28 < x < 1.3e37

   1. Initial program 99.8%

      \[\frac{\left(x - 2\right) \cdot \left(\left(\left(\left(x \cdot 4.16438922228 + 78.6994924154\right) \cdot x + 137.519416416\right) \cdot x + y\right) \cdot x + z\right)}{\left(\left(\left(x + 43.3400022514\right) \cdot x + 263.505074721\right) \cdot x + 313.399215894\right) \cdot x + 47.066876606} \]
   2. Add Preprocessing

   3. Taylor expanded in y around 0

      \[\leadsto \color{blue}{\frac{\left(z + {x}^{2} \cdot \left(\frac{4297481763}{31250000} + x \cdot \left(\frac{393497462077}{5000000000} + \frac{104109730557}{25000000000} \cdot x\right)\right)\right) \cdot \left(x - 2\right)}{\frac{23533438303}{500000000} + x \cdot \left(\frac{156699607947}{500000000} + x \cdot \left(\frac{263505074721}{1000000000} + x \cdot \left(\frac{216700011257}{5000000000} + x\right)\right)\right)}} \]
   4. Step-by-step derivation

      1. lower-/.f64 N/A

         \[\leadsto \color{blue}{\frac{\left(z + {x}^{2} \cdot \left(\frac{4297481763}{31250000} + x \cdot \left(\frac{393497462077}{5000000000} + \frac{104109730557}{25000000000} \cdot x\right)\right)\right) \cdot \left(x - 2\right)}{\frac{23533438303}{500000000} + x \cdot \left(\frac{156699607947}{500000000} + x \cdot \left(\frac{263505074721}{1000000000} + x \cdot \left(\frac{216700011257}{5000000000} + x\right)\right)\right)}} \]

   5. Applied rewrites 70.9%

      \[\leadsto \color{blue}{\frac{\mathsf{fma}\left(x \cdot x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, 4.16438922228, 78.6994924154\right), 137.519416416\right), z\right) \cdot \left(x + -2\right)}{\mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, x + 43.3400022514, 263.505074721\right), 313.399215894\right), 47.066876606\right)}} \]

   6. Taylor expanded in x around 0

      \[\leadsto \frac{-1000000000}{23533438303} \cdot z + \color{blue}{x \cdot \left(\frac{500000000}{23533438303} \cdot z - \frac{-156699607947000000000}{553822718361107519809} \cdot z\right)} \]
   7. Step-by-step derivation

   8. Applied rewrites 64.3%

      \[\leadsto \mathsf{fma}\left(-0.0424927283095952, \color{blue}{z}, \left(z \cdot x\right) \cdot 0.3041881842569256\right) \]
3. Recombined 2 regimes into one program.
4. Add Preprocessing

Alternative 11: 76.1% accurate, 2.9× speedup

Math

\[\begin{array}{l} \\ \begin{array}{l} t_0 := \frac{x + -2}{0.24013125253755718}\\ \mathbf{if}\;x \leq -6 \cdot 10^{+28}:\\ \;\;\;\;t\_0\\ \mathbf{elif}\;x \leq 1.3 \cdot 10^{+37}:\\ \;\;\;\;z \cdot \mathsf{fma}\left(0.3041881842569256, x, -0.0424927283095952\right)\\ \mathbf{else}:\\ \;\;\;\;t\_0\\ \end{array} \end{array} \]

FPCore

(FPCore (x y z)
 :precision binary64
 (let* ((t_0 (/ (+ x -2.0) 0.24013125253755718)))
   (if (<= x -6e+28)
     t_0
     (if (<= x 1.3e+37)
       (* z (fma 0.3041881842569256 x -0.0424927283095952))
       t_0))))

C

double code(double x, double y, double z) {
	double t_0 = (x + -2.0) / 0.24013125253755718;
	double tmp;
	if (x <= -6e+28) {
		tmp = t_0;
	} else if (x <= 1.3e+37) {
		tmp = z * fma(0.3041881842569256, x, -0.0424927283095952);
	} else {
		tmp = t_0;
	}
	return tmp;
}

Julia

function code(x, y, z)
	t_0 = Float64(Float64(x + -2.0) / 0.24013125253755718)
	tmp = 0.0
	if (x <= -6e+28)
		tmp = t_0;
	elseif (x <= 1.3e+37)
		tmp = Float64(z * fma(0.3041881842569256, x, -0.0424927283095952));
	else
		tmp = t_0;
	end
	return tmp
end

Wolfram

code[x_, y_, z_] := Block[{t$95$0 = N[(N[(x + -2.0), $MachinePrecision] / 0.24013125253755718), $MachinePrecision]}, If[LessEqual[x, -6e+28], t$95$0, If[LessEqual[x, 1.3e+37], N[(z * N[(0.3041881842569256 * x + -0.0424927283095952), $MachinePrecision]), $MachinePrecision], t$95$0]]]

Derivation

1. Split input into 2 regimes

2. if x < -6.0000000000000002e28 or 1.3e37 < x

   1. Initial program 8.5%

      \[\frac{\left(x - 2\right) \cdot \left(\left(\left(\left(x \cdot 4.16438922228 + 78.6994924154\right) \cdot x + 137.519416416\right) \cdot x + y\right) \cdot x + z\right)}{\left(\left(\left(x + 43.3400022514\right) \cdot x + 263.505074721\right) \cdot x + 313.399215894\right) \cdot x + 47.066876606} \]
   2. Add Preprocessing

   4. Applied rewrites 17.1%

      \[\leadsto \color{blue}{\frac{x + -2}{\frac{\mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, x + 43.3400022514, 263.505074721\right), 313.399215894\right), 47.066876606\right)}{\mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, 4.16438922228, 78.6994924154\right), 137.519416416\right), y\right), z\right)}}} \]

   5. Taylor expanded in x around inf

      \[\leadsto \frac{x + -2}{\color{blue}{\frac{25000000000}{104109730557}}} \]
   6. Step-by-step derivation

   7. Applied rewrites 92.7%

      \[\leadsto \frac{x + -2}{\color{blue}{0.24013125253755718}} \]

   if -6.0000000000000002e28 < x < 1.3e37

   1. Initial program 99.8%

      \[\frac{\left(x - 2\right) \cdot \left(\left(\left(\left(x \cdot 4.16438922228 + 78.6994924154\right) \cdot x + 137.519416416\right) \cdot x + y\right) \cdot x + z\right)}{\left(\left(\left(x + 43.3400022514\right) \cdot x + 263.505074721\right) \cdot x + 313.399215894\right) \cdot x + 47.066876606} \]
   2. Add Preprocessing

   4. Applied rewrites 99.7%

      \[\leadsto \color{blue}{\frac{\mathsf{fma}\left(x, x, -4\right) \cdot \frac{\mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, 4.16438922228, 78.6994924154\right), 137.519416416\right), y\right), z\right)}{\mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, x + 43.3400022514, 263.505074721\right), 313.399215894\right), 47.066876606\right)}}{x + 2}} \]

   5. Taylor expanded in z around inf

      \[\leadsto \color{blue}{\frac{z \cdot \left({x}^{2} - 4\right)}{\left(2 + x\right) \cdot \left(\frac{23533438303}{500000000} + x \cdot \left(\frac{156699607947}{500000000} + x \cdot \left(\frac{263505074721}{1000000000} + x \cdot \left(\frac{216700011257}{5000000000} + x\right)\right)\right)\right)}} \]
   6. Step-by-step derivation

      1. lower-/.f64 N/A

         \[\leadsto \color{blue}{\frac{z \cdot \left({x}^{2} - 4\right)}{\left(2 + x\right) \cdot \left(\frac{23533438303}{500000000} + x \cdot \left(\frac{156699607947}{500000000} + x \cdot \left(\frac{263505074721}{1000000000} + x \cdot \left(\frac{216700011257}{5000000000} + x\right)\right)\right)\right)}} \]

      2. lower-*.f64 N/A

         \[\leadsto \frac{\color{blue}{z \cdot \left({x}^{2} - 4\right)}}{\left(2 + x\right) \cdot \left(\frac{23533438303}{500000000} + x \cdot \left(\frac{156699607947}{500000000} + x \cdot \left(\frac{263505074721}{1000000000} + x \cdot \left(\frac{216700011257}{5000000000} + x\right)\right)\right)\right)} \]

      3. sub-neg N/A

         \[\leadsto \frac{z \cdot \color{blue}{\left({x}^{2} + \left(\mathsf{neg}\left(4\right)\right)\right)}}{\left(2 + x\right) \cdot \left(\frac{23533438303}{500000000} + x \cdot \left(\frac{156699607947}{500000000} + x \cdot \left(\frac{263505074721}{1000000000} + x \cdot \left(\frac{216700011257}{5000000000} + x\right)\right)\right)\right)} \]

      4. unpow2 N/A

         \[\leadsto \frac{z \cdot \left(\color{blue}{x \cdot x} + \left(\mathsf{neg}\left(4\right)\right)\right)}{\left(2 + x\right) \cdot \left(\frac{23533438303}{500000000} + x \cdot \left(\frac{156699607947}{500000000} + x \cdot \left(\frac{263505074721}{1000000000} + x \cdot \left(\frac{216700011257}{5000000000} + x\right)\right)\right)\right)} \]

      5. lower-fma.f64 N/A

         \[\leadsto \frac{z \cdot \color{blue}{\mathsf{fma}\left(x, x, \mathsf{neg}\left(4\right)\right)}}{\left(2 + x\right) \cdot \left(\frac{23533438303}{500000000} + x \cdot \left(\frac{156699607947}{500000000} + x \cdot \left(\frac{263505074721}{1000000000} + x \cdot \left(\frac{216700011257}{5000000000} + x\right)\right)\right)\right)} \]

      6. metadata-eval N/A

         \[\leadsto \frac{z \cdot \mathsf{fma}\left(x, x, \color{blue}{-4}\right)}{\left(2 + x\right) \cdot \left(\frac{23533438303}{500000000} + x \cdot \left(\frac{156699607947}{500000000} + x \cdot \left(\frac{263505074721}{1000000000} + x \cdot \left(\frac{216700011257}{5000000000} + x\right)\right)\right)\right)} \]

      7. *-commutative N/A

         \[\leadsto \frac{z \cdot \mathsf{fma}\left(x, x, -4\right)}{\color{blue}{\left(\frac{23533438303}{500000000} + x \cdot \left(\frac{156699607947}{500000000} + x \cdot \left(\frac{263505074721}{1000000000} + x \cdot \left(\frac{216700011257}{5000000000} + x\right)\right)\right)\right) \cdot \left(2 + x\right)}} \]

      8. lower-*.f64 N/A

         \[\leadsto \frac{z \cdot \mathsf{fma}\left(x, x, -4\right)}{\color{blue}{\left(\frac{23533438303}{500000000} + x \cdot \left(\frac{156699607947}{500000000} + x \cdot \left(\frac{263505074721}{1000000000} + x \cdot \left(\frac{216700011257}{5000000000} + x\right)\right)\right)\right) \cdot \left(2 + x\right)}} \]

      9. +-commutative N/A

         \[\leadsto \frac{z \cdot \mathsf{fma}\left(x, x, -4\right)}{\color{blue}{\left(x \cdot \left(\frac{156699607947}{500000000} + x \cdot \left(\frac{263505074721}{1000000000} + x \cdot \left(\frac{216700011257}{5000000000} + x\right)\right)\right) + \frac{23533438303}{500000000}\right)} \cdot \left(2 + x\right)} \]

      10. lower-fma.f64 N/A

          \[\leadsto \frac{z \cdot \mathsf{fma}\left(x, x, -4\right)}{\color{blue}{\mathsf{fma}\left(x, \frac{156699607947}{500000000} + x \cdot \left(\frac{263505074721}{1000000000} + x \cdot \left(\frac{216700011257}{5000000000} + x\right)\right), \frac{23533438303}{500000000}\right)} \cdot \left(2 + x\right)} \]

      11. +-commutative N/A

          \[\leadsto \frac{z \cdot \mathsf{fma}\left(x, x, -4\right)}{\mathsf{fma}\left(x, \color{blue}{x \cdot \left(\frac{263505074721}{1000000000} + x \cdot \left(\frac{216700011257}{5000000000} + x\right)\right) + \frac{156699607947}{500000000}}, \frac{23533438303}{500000000}\right) \cdot \left(2 + x\right)} \]

      12. lower-fma.f64 N/A

          \[\leadsto \frac{z \cdot \mathsf{fma}\left(x, x, -4\right)}{\mathsf{fma}\left(x, \color{blue}{\mathsf{fma}\left(x, \frac{263505074721}{1000000000} + x \cdot \left(\frac{216700011257}{5000000000} + x\right), \frac{156699607947}{500000000}\right)}, \frac{23533438303}{500000000}\right) \cdot \left(2 + x\right)} \]

      13. +-commutative N/A

          \[\leadsto \frac{z \cdot \mathsf{fma}\left(x, x, -4\right)}{\mathsf{fma}\left(x, \mathsf{fma}\left(x, \color{blue}{x \cdot \left(\frac{216700011257}{5000000000} + x\right) + \frac{263505074721}{1000000000}}, \frac{156699607947}{500000000}\right), \frac{23533438303}{500000000}\right) \cdot \left(2 + x\right)} \]

      14. lower-fma.f64 N/A

          \[\leadsto \frac{z \cdot \mathsf{fma}\left(x, x, -4\right)}{\mathsf{fma}\left(x, \mathsf{fma}\left(x, \color{blue}{\mathsf{fma}\left(x, \frac{216700011257}{5000000000} + x, \frac{263505074721}{1000000000}\right)}, \frac{156699607947}{500000000}\right), \frac{23533438303}{500000000}\right) \cdot \left(2 + x\right)} \]

      15. +-commutative N/A

          \[\leadsto \frac{z \cdot \mathsf{fma}\left(x, x, -4\right)}{\mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, \color{blue}{x + \frac{216700011257}{5000000000}}, \frac{263505074721}{1000000000}\right), \frac{156699607947}{500000000}\right), \frac{23533438303}{500000000}\right) \cdot \left(2 + x\right)} \]

      16. lower-+.f64 N/A

          \[\leadsto \frac{z \cdot \mathsf{fma}\left(x, x, -4\right)}{\mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, \color{blue}{x + \frac{216700011257}{5000000000}}, \frac{263505074721}{1000000000}\right), \frac{156699607947}{500000000}\right), \frac{23533438303}{500000000}\right) \cdot \left(2 + x\right)} \]

      17. +-commutative N/A

          \[\leadsto \frac{z \cdot \mathsf{fma}\left(x, x, -4\right)}{\mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, x + \frac{216700011257}{5000000000}, \frac{263505074721}{1000000000}\right), \frac{156699607947}{500000000}\right), \frac{23533438303}{500000000}\right) \cdot \color{blue}{\left(x + 2\right)}} \]

      18. lower-+.f64 66.5

          \[\leadsto \frac{z \cdot \mathsf{fma}\left(x, x, -4\right)}{\mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, x + 43.3400022514, 263.505074721\right), 313.399215894\right), 47.066876606\right) \cdot \color{blue}{\left(x + 2\right)}} \]

   7. Applied rewrites 66.5%

      \[\leadsto \color{blue}{\frac{z \cdot \mathsf{fma}\left(x, x, -4\right)}{\mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, x + 43.3400022514, 263.505074721\right), 313.399215894\right), 47.066876606\right) \cdot \left(x + 2\right)}} \]

   8. Taylor expanded in x around 0

      \[\leadsto \frac{-1000000000}{23533438303} \cdot z + \color{blue}{\frac{168466327098500000000}{553822718361107519809} \cdot \left(x \cdot z\right)} \]
   9. Step-by-step derivation

   10. Applied rewrites 64.3%

       \[\leadsto z \cdot \color{blue}{\mathsf{fma}\left(0.3041881842569256, x, -0.0424927283095952\right)} \]
3. Recombined 2 regimes into one program.
4. Add Preprocessing

Alternative 12: 75.8% accurate, 3.3× speedup

Math

\[\begin{array}{l} \\ \begin{array}{l} \mathbf{if}\;x \leq -6 \cdot 10^{+28}:\\ \;\;\;\;x \cdot 4.16438922228\\ \mathbf{elif}\;x \leq 1.3 \cdot 10^{+37}:\\ \;\;\;\;z \cdot \mathsf{fma}\left(0.3041881842569256, x, -0.0424927283095952\right)\\ \mathbf{else}:\\ \;\;\;\;4.16438922228 \cdot \left(x + -2\right)\\ \end{array} \end{array} \]

FPCore

(FPCore (x y z)
 :precision binary64
 (if (<= x -6e+28)
   (* x 4.16438922228)
   (if (<= x 1.3e+37)
     (* z (fma 0.3041881842569256 x -0.0424927283095952))
     (* 4.16438922228 (+ x -2.0)))))

C

double code(double x, double y, double z) {
	double tmp;
	if (x <= -6e+28) {
		tmp = x * 4.16438922228;
	} else if (x <= 1.3e+37) {
		tmp = z * fma(0.3041881842569256, x, -0.0424927283095952);
	} else {
		tmp = 4.16438922228 * (x + -2.0);
	}
	return tmp;
}

Julia

function code(x, y, z)
	tmp = 0.0
	if (x <= -6e+28)
		tmp = Float64(x * 4.16438922228);
	elseif (x <= 1.3e+37)
		tmp = Float64(z * fma(0.3041881842569256, x, -0.0424927283095952));
	else
		tmp = Float64(4.16438922228 * Float64(x + -2.0));
	end
	return tmp
end

Wolfram

code[x_, y_, z_] := If[LessEqual[x, -6e+28], N[(x * 4.16438922228), $MachinePrecision], If[LessEqual[x, 1.3e+37], N[(z * N[(0.3041881842569256 * x + -0.0424927283095952), $MachinePrecision]), $MachinePrecision], N[(4.16438922228 * N[(x + -2.0), $MachinePrecision]), $MachinePrecision]]]

Derivation

1. Split input into 3 regimes

2. if x < -6.0000000000000002e28

   1. Initial program 13.9%

      \[\frac{\left(x - 2\right) \cdot \left(\left(\left(\left(x \cdot 4.16438922228 + 78.6994924154\right) \cdot x + 137.519416416\right) \cdot x + y\right) \cdot x + z\right)}{\left(\left(\left(x + 43.3400022514\right) \cdot x + 263.505074721\right) \cdot x + 313.399215894\right) \cdot x + 47.066876606} \]
   2. Add Preprocessing

   3. Taylor expanded in x around inf

      \[\leadsto \color{blue}{\frac{104109730557}{25000000000} \cdot x} \]
   4. Step-by-step derivation

      1. *-commutative N/A

         \[\leadsto \color{blue}{x \cdot \frac{104109730557}{25000000000}} \]

      2. lower-*.f64 88.1

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

   5. Applied rewrites 88.1%

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

   if -6.0000000000000002e28 < x < 1.3e37

   1. Initial program 99.8%

      \[\frac{\left(x - 2\right) \cdot \left(\left(\left(\left(x \cdot 4.16438922228 + 78.6994924154\right) \cdot x + 137.519416416\right) \cdot x + y\right) \cdot x + z\right)}{\left(\left(\left(x + 43.3400022514\right) \cdot x + 263.505074721\right) \cdot x + 313.399215894\right) \cdot x + 47.066876606} \]
   2. Add Preprocessing

   4. Applied rewrites 99.7%

      \[\leadsto \color{blue}{\frac{\mathsf{fma}\left(x, x, -4\right) \cdot \frac{\mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, 4.16438922228, 78.6994924154\right), 137.519416416\right), y\right), z\right)}{\mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, x + 43.3400022514, 263.505074721\right), 313.399215894\right), 47.066876606\right)}}{x + 2}} \]

   5. Taylor expanded in z around inf

      \[\leadsto \color{blue}{\frac{z \cdot \left({x}^{2} - 4\right)}{\left(2 + x\right) \cdot \left(\frac{23533438303}{500000000} + x \cdot \left(\frac{156699607947}{500000000} + x \cdot \left(\frac{263505074721}{1000000000} + x \cdot \left(\frac{216700011257}{5000000000} + x\right)\right)\right)\right)}} \]
   6. Step-by-step derivation

      1. lower-/.f64 N/A

         \[\leadsto \color{blue}{\frac{z \cdot \left({x}^{2} - 4\right)}{\left(2 + x\right) \cdot \left(\frac{23533438303}{500000000} + x \cdot \left(\frac{156699607947}{500000000} + x \cdot \left(\frac{263505074721}{1000000000} + x \cdot \left(\frac{216700011257}{5000000000} + x\right)\right)\right)\right)}} \]

      2. lower-*.f64 N/A

         \[\leadsto \frac{\color{blue}{z \cdot \left({x}^{2} - 4\right)}}{\left(2 + x\right) \cdot \left(\frac{23533438303}{500000000} + x \cdot \left(\frac{156699607947}{500000000} + x \cdot \left(\frac{263505074721}{1000000000} + x \cdot \left(\frac{216700011257}{5000000000} + x\right)\right)\right)\right)} \]

      3. sub-neg N/A

         \[\leadsto \frac{z \cdot \color{blue}{\left({x}^{2} + \left(\mathsf{neg}\left(4\right)\right)\right)}}{\left(2 + x\right) \cdot \left(\frac{23533438303}{500000000} + x \cdot \left(\frac{156699607947}{500000000} + x \cdot \left(\frac{263505074721}{1000000000} + x \cdot \left(\frac{216700011257}{5000000000} + x\right)\right)\right)\right)} \]

      4. unpow2 N/A

         \[\leadsto \frac{z \cdot \left(\color{blue}{x \cdot x} + \left(\mathsf{neg}\left(4\right)\right)\right)}{\left(2 + x\right) \cdot \left(\frac{23533438303}{500000000} + x \cdot \left(\frac{156699607947}{500000000} + x \cdot \left(\frac{263505074721}{1000000000} + x \cdot \left(\frac{216700011257}{5000000000} + x\right)\right)\right)\right)} \]

      5. lower-fma.f64 N/A

         \[\leadsto \frac{z \cdot \color{blue}{\mathsf{fma}\left(x, x, \mathsf{neg}\left(4\right)\right)}}{\left(2 + x\right) \cdot \left(\frac{23533438303}{500000000} + x \cdot \left(\frac{156699607947}{500000000} + x \cdot \left(\frac{263505074721}{1000000000} + x \cdot \left(\frac{216700011257}{5000000000} + x\right)\right)\right)\right)} \]

      6. metadata-eval N/A

         \[\leadsto \frac{z \cdot \mathsf{fma}\left(x, x, \color{blue}{-4}\right)}{\left(2 + x\right) \cdot \left(\frac{23533438303}{500000000} + x \cdot \left(\frac{156699607947}{500000000} + x \cdot \left(\frac{263505074721}{1000000000} + x \cdot \left(\frac{216700011257}{5000000000} + x\right)\right)\right)\right)} \]

      7. *-commutative N/A

         \[\leadsto \frac{z \cdot \mathsf{fma}\left(x, x, -4\right)}{\color{blue}{\left(\frac{23533438303}{500000000} + x \cdot \left(\frac{156699607947}{500000000} + x \cdot \left(\frac{263505074721}{1000000000} + x \cdot \left(\frac{216700011257}{5000000000} + x\right)\right)\right)\right) \cdot \left(2 + x\right)}} \]

      8. lower-*.f64 N/A

         \[\leadsto \frac{z \cdot \mathsf{fma}\left(x, x, -4\right)}{\color{blue}{\left(\frac{23533438303}{500000000} + x \cdot \left(\frac{156699607947}{500000000} + x \cdot \left(\frac{263505074721}{1000000000} + x \cdot \left(\frac{216700011257}{5000000000} + x\right)\right)\right)\right) \cdot \left(2 + x\right)}} \]

      9. +-commutative N/A

         \[\leadsto \frac{z \cdot \mathsf{fma}\left(x, x, -4\right)}{\color{blue}{\left(x \cdot \left(\frac{156699607947}{500000000} + x \cdot \left(\frac{263505074721}{1000000000} + x \cdot \left(\frac{216700011257}{5000000000} + x\right)\right)\right) + \frac{23533438303}{500000000}\right)} \cdot \left(2 + x\right)} \]

      10. lower-fma.f64 N/A

          \[\leadsto \frac{z \cdot \mathsf{fma}\left(x, x, -4\right)}{\color{blue}{\mathsf{fma}\left(x, \frac{156699607947}{500000000} + x \cdot \left(\frac{263505074721}{1000000000} + x \cdot \left(\frac{216700011257}{5000000000} + x\right)\right), \frac{23533438303}{500000000}\right)} \cdot \left(2 + x\right)} \]

      11. +-commutative N/A

          \[\leadsto \frac{z \cdot \mathsf{fma}\left(x, x, -4\right)}{\mathsf{fma}\left(x, \color{blue}{x \cdot \left(\frac{263505074721}{1000000000} + x \cdot \left(\frac{216700011257}{5000000000} + x\right)\right) + \frac{156699607947}{500000000}}, \frac{23533438303}{500000000}\right) \cdot \left(2 + x\right)} \]

      12. lower-fma.f64 N/A

          \[\leadsto \frac{z \cdot \mathsf{fma}\left(x, x, -4\right)}{\mathsf{fma}\left(x, \color{blue}{\mathsf{fma}\left(x, \frac{263505074721}{1000000000} + x \cdot \left(\frac{216700011257}{5000000000} + x\right), \frac{156699607947}{500000000}\right)}, \frac{23533438303}{500000000}\right) \cdot \left(2 + x\right)} \]

      13. +-commutative N/A

          \[\leadsto \frac{z \cdot \mathsf{fma}\left(x, x, -4\right)}{\mathsf{fma}\left(x, \mathsf{fma}\left(x, \color{blue}{x \cdot \left(\frac{216700011257}{5000000000} + x\right) + \frac{263505074721}{1000000000}}, \frac{156699607947}{500000000}\right), \frac{23533438303}{500000000}\right) \cdot \left(2 + x\right)} \]

      14. lower-fma.f64 N/A

          \[\leadsto \frac{z \cdot \mathsf{fma}\left(x, x, -4\right)}{\mathsf{fma}\left(x, \mathsf{fma}\left(x, \color{blue}{\mathsf{fma}\left(x, \frac{216700011257}{5000000000} + x, \frac{263505074721}{1000000000}\right)}, \frac{156699607947}{500000000}\right), \frac{23533438303}{500000000}\right) \cdot \left(2 + x\right)} \]

      15. +-commutative N/A

          \[\leadsto \frac{z \cdot \mathsf{fma}\left(x, x, -4\right)}{\mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, \color{blue}{x + \frac{216700011257}{5000000000}}, \frac{263505074721}{1000000000}\right), \frac{156699607947}{500000000}\right), \frac{23533438303}{500000000}\right) \cdot \left(2 + x\right)} \]

      16. lower-+.f64 N/A

          \[\leadsto \frac{z \cdot \mathsf{fma}\left(x, x, -4\right)}{\mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, \color{blue}{x + \frac{216700011257}{5000000000}}, \frac{263505074721}{1000000000}\right), \frac{156699607947}{500000000}\right), \frac{23533438303}{500000000}\right) \cdot \left(2 + x\right)} \]

      17. +-commutative N/A

          \[\leadsto \frac{z \cdot \mathsf{fma}\left(x, x, -4\right)}{\mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, x + \frac{216700011257}{5000000000}, \frac{263505074721}{1000000000}\right), \frac{156699607947}{500000000}\right), \frac{23533438303}{500000000}\right) \cdot \color{blue}{\left(x + 2\right)}} \]

      18. lower-+.f64 66.5

          \[\leadsto \frac{z \cdot \mathsf{fma}\left(x, x, -4\right)}{\mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, x + 43.3400022514, 263.505074721\right), 313.399215894\right), 47.066876606\right) \cdot \color{blue}{\left(x + 2\right)}} \]

   7. Applied rewrites 66.5%

      \[\leadsto \color{blue}{\frac{z \cdot \mathsf{fma}\left(x, x, -4\right)}{\mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, x + 43.3400022514, 263.505074721\right), 313.399215894\right), 47.066876606\right) \cdot \left(x + 2\right)}} \]

   8. Taylor expanded in x around 0

      \[\leadsto \frac{-1000000000}{23533438303} \cdot z + \color{blue}{\frac{168466327098500000000}{553822718361107519809} \cdot \left(x \cdot z\right)} \]
   9. Step-by-step derivation

   10. Applied rewrites 64.3%

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

   if 1.3e37 < x

   1. Initial program 2.2%

      \[\frac{\left(x - 2\right) \cdot \left(\left(\left(\left(x \cdot 4.16438922228 + 78.6994924154\right) \cdot x + 137.519416416\right) \cdot x + y\right) \cdot x + z\right)}{\left(\left(\left(x + 43.3400022514\right) \cdot x + 263.505074721\right) \cdot x + 313.399215894\right) \cdot x + 47.066876606} \]
   2. Add Preprocessing
   3. Step-by-step derivation

      1. lift-/.f64 N/A

         \[\leadsto \color{blue}{\frac{\left(x - 2\right) \cdot \left(\left(\left(\left(x \cdot \frac{104109730557}{25000000000} + \frac{393497462077}{5000000000}\right) \cdot x + \frac{4297481763}{31250000}\right) \cdot x + y\right) \cdot x + z\right)}{\left(\left(\left(x + \frac{216700011257}{5000000000}\right) \cdot x + \frac{263505074721}{1000000000}\right) \cdot x + \frac{156699607947}{500000000}\right) \cdot x + \frac{23533438303}{500000000}}} \]

      2. lift-*.f64 N/A

         \[\leadsto \frac{\color{blue}{\left(x - 2\right) \cdot \left(\left(\left(\left(x \cdot \frac{104109730557}{25000000000} + \frac{393497462077}{5000000000}\right) \cdot x + \frac{4297481763}{31250000}\right) \cdot x + y\right) \cdot x + z\right)}}{\left(\left(\left(x + \frac{216700011257}{5000000000}\right) \cdot x + \frac{263505074721}{1000000000}\right) \cdot x + \frac{156699607947}{500000000}\right) \cdot x + \frac{23533438303}{500000000}} \]

      3. associate-/l* N/A

         \[\leadsto \color{blue}{\left(x - 2\right) \cdot \frac{\left(\left(\left(x \cdot \frac{104109730557}{25000000000} + \frac{393497462077}{5000000000}\right) \cdot x + \frac{4297481763}{31250000}\right) \cdot x + y\right) \cdot x + z}{\left(\left(\left(x + \frac{216700011257}{5000000000}\right) \cdot x + \frac{263505074721}{1000000000}\right) \cdot x + \frac{156699607947}{500000000}\right) \cdot x + \frac{23533438303}{500000000}}} \]

      4. *-commutative N/A

         \[\leadsto \color{blue}{\frac{\left(\left(\left(x \cdot \frac{104109730557}{25000000000} + \frac{393497462077}{5000000000}\right) \cdot x + \frac{4297481763}{31250000}\right) \cdot x + y\right) \cdot x + z}{\left(\left(\left(x + \frac{216700011257}{5000000000}\right) \cdot x + \frac{263505074721}{1000000000}\right) \cdot x + \frac{156699607947}{500000000}\right) \cdot x + \frac{23533438303}{500000000}} \cdot \left(x - 2\right)} \]

      5. lower-*.f64 N/A

         \[\leadsto \color{blue}{\frac{\left(\left(\left(x \cdot \frac{104109730557}{25000000000} + \frac{393497462077}{5000000000}\right) \cdot x + \frac{4297481763}{31250000}\right) \cdot x + y\right) \cdot x + z}{\left(\left(\left(x + \frac{216700011257}{5000000000}\right) \cdot x + \frac{263505074721}{1000000000}\right) \cdot x + \frac{156699607947}{500000000}\right) \cdot x + \frac{23533438303}{500000000}} \cdot \left(x - 2\right)} \]

   4. Applied rewrites 9.7%

      \[\leadsto \color{blue}{\frac{\mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, 4.16438922228, 78.6994924154\right), 137.519416416\right), y\right), z\right)}{\mathsf{fma}\left(x, \mathsf{fma}\left(x, \mathsf{fma}\left(x, x + 43.3400022514, 263.505074721\right), 313.399215894\right), 47.066876606\right)} \cdot \left(x + -2\right)} \]

   5. Taylor expanded in x around inf

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

   7. Applied rewrites 96.9%

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

Alternative 13: 76.3% accurate, 4.4× speedup

Math

\[\begin{array}{l} \\ \begin{array}{l} \mathbf{if}\;x \leq -6 \cdot 10^{+28}:\\ \;\;\;\;x \cdot 4.16438922228\\ \mathbf{elif}\;x \leq 17500000000:\\ \;\;\;\;z \cdot -0.0424927283095952\\ \mathbf{else}:\\ \;\;\;\;x \cdot 4.16438922228\\ \end{array} \end{array} \]

FPCore

(FPCore (x y z)
 :precision binary64
 (if (<= x -6e+28)
   (* x 4.16438922228)
   (if (<= x 17500000000.0) (* z -0.0424927283095952) (* x 4.16438922228))))

C

double code(double x, double y, double z) {
	double tmp;
	if (x <= -6e+28) {
		tmp = x * 4.16438922228;
	} else if (x <= 17500000000.0) {
		tmp = z * -0.0424927283095952;
	} else {
		tmp = x * 4.16438922228;
	}
	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 <= (-6d+28)) then
        tmp = x * 4.16438922228d0
    else if (x <= 17500000000.0d0) then
        tmp = z * (-0.0424927283095952d0)
    else
        tmp = x * 4.16438922228d0
    end if
    code = tmp
end function

Java

public static double code(double x, double y, double z) {
	double tmp;
	if (x <= -6e+28) {
		tmp = x * 4.16438922228;
	} else if (x <= 17500000000.0) {
		tmp = z * -0.0424927283095952;
	} else {
		tmp = x * 4.16438922228;
	}
	return tmp;
}

Python

def code(x, y, z):
	tmp = 0
	if x <= -6e+28:
		tmp = x * 4.16438922228
	elif x <= 17500000000.0:
		tmp = z * -0.0424927283095952
	else:
		tmp = x * 4.16438922228
	return tmp

Julia

function code(x, y, z)
	tmp = 0.0
	if (x <= -6e+28)
		tmp = Float64(x * 4.16438922228);
	elseif (x <= 17500000000.0)
		tmp = Float64(z * -0.0424927283095952);
	else
		tmp = Float64(x * 4.16438922228);
	end
	return tmp
end

MATLAB

function tmp_2 = code(x, y, z)
	tmp = 0.0;
	if (x <= -6e+28)
		tmp = x * 4.16438922228;
	elseif (x <= 17500000000.0)
		tmp = z * -0.0424927283095952;
	else
		tmp = x * 4.16438922228;
	end
	tmp_2 = tmp;
end

Wolfram

code[x_, y_, z_] := If[LessEqual[x, -6e+28], N[(x * 4.16438922228), $MachinePrecision], If[LessEqual[x, 17500000000.0], N[(z * -0.0424927283095952), $MachinePrecision], N[(x * 4.16438922228), $MachinePrecision]]]

Derivation

1. Split input into 2 regimes

2. if x < -6.0000000000000002e28 or 1.75e10 < x

   1. Initial program 10.9%

      \[\frac{\left(x - 2\right) \cdot \left(\left(\left(\left(x \cdot 4.16438922228 + 78.6994924154\right) \cdot x + 137.519416416\right) \cdot x + y\right) \cdot x + z\right)}{\left(\left(\left(x + 43.3400022514\right) \cdot x + 263.505074721\right) \cdot x + 313.399215894\right) \cdot x + 47.066876606} \]
   2. Add Preprocessing

   3. Taylor expanded in x around inf

      \[\leadsto \color{blue}{\frac{104109730557}{25000000000} \cdot x} \]
   4. Step-by-step derivation

      1. *-commutative N/A

         \[\leadsto \color{blue}{x \cdot \frac{104109730557}{25000000000}} \]

      2. lower-*.f64 89.8

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

   5. Applied rewrites 89.8%

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

   if -6.0000000000000002e28 < x < 1.75e10

   1. Initial program 99.8%

      \[\frac{\left(x - 2\right) \cdot \left(\left(\left(\left(x \cdot 4.16438922228 + 78.6994924154\right) \cdot x + 137.519416416\right) \cdot x + y\right) \cdot x + z\right)}{\left(\left(\left(x + 43.3400022514\right) \cdot x + 263.505074721\right) \cdot x + 313.399215894\right) \cdot x + 47.066876606} \]
   2. Add Preprocessing

   3. Taylor expanded in x around 0

      \[\leadsto \color{blue}{\frac{-1000000000}{23533438303} \cdot z} \]
   4. Step-by-step derivation

      1. *-commutative N/A

         \[\leadsto \color{blue}{z \cdot \frac{-1000000000}{23533438303}} \]

      2. lower-*.f64 65.3

         \[\leadsto \color{blue}{z \cdot -0.0424927283095952} \]

   5. Applied rewrites 65.3%

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

Alternative 14: 44.6% accurate, 13.2× speedup

Math

\[\begin{array}{l} \\ x \cdot 4.16438922228 \end{array} \]

FPCore

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

C

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

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 * 4.16438922228d0
end function

Java

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

Python

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

Julia

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

MATLAB

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

Wolfram

code[x_, y_, z_] := N[(x * 4.16438922228), $MachinePrecision]

Derivation

1. Initial program 60.5%

   \[\frac{\left(x - 2\right) \cdot \left(\left(\left(\left(x \cdot 4.16438922228 + 78.6994924154\right) \cdot x + 137.519416416\right) \cdot x + y\right) \cdot x + z\right)}{\left(\left(\left(x + 43.3400022514\right) \cdot x + 263.505074721\right) \cdot x + 313.399215894\right) \cdot x + 47.066876606} \]
2. Add Preprocessing

3. Taylor expanded in x around inf

   \[\leadsto \color{blue}{\frac{104109730557}{25000000000} \cdot x} \]
4. Step-by-step derivation

   1. *-commutative N/A

      \[\leadsto \color{blue}{x \cdot \frac{104109730557}{25000000000}} \]

   2. lower-*.f64 41.5

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

5. Applied rewrites 41.5%

   \[\leadsto \color{blue}{x \cdot 4.16438922228} \]
6. Add Preprocessing

Developer Target 1: 98.8% accurate, 0.7× speedup

Math

\[\begin{array}{l} \\ \begin{array}{l} t_0 := \left(\frac{y}{x \cdot x} + 4.16438922228 \cdot x\right) - 110.1139242984811\\ \mathbf{if}\;x < -3.326128725870005 \cdot 10^{+62}:\\ \;\;\;\;t\_0\\ \mathbf{elif}\;x < 9.429991714554673 \cdot 10^{+55}:\\ \;\;\;\;\frac{x - 2}{1} \cdot \frac{\left(\left(\left(x \cdot 4.16438922228 + 78.6994924154\right) \cdot x + 137.519416416\right) \cdot x + y\right) \cdot x + z}{\left(\left(263.505074721 \cdot x + \left(43.3400022514 \cdot \left(x \cdot x\right) + x \cdot \left(x \cdot x\right)\right)\right) + 313.399215894\right) \cdot x + 47.066876606}\\ \mathbf{else}:\\ \;\;\;\;t\_0\\ \end{array} \end{array} \]

FPCore

(FPCore (x y z)
 :precision binary64
 (let* ((t_0 (- (+ (/ y (* x x)) (* 4.16438922228 x)) 110.1139242984811)))
   (if (< x -3.326128725870005e+62)
     t_0
     (if (< x 9.429991714554673e+55)
       (*
        (/ (- x 2.0) 1.0)
        (/
         (+
          (*
           (+
            (* (+ (* (+ (* x 4.16438922228) 78.6994924154) x) 137.519416416) x)
            y)
           x)
          z)
         (+
          (*
           (+
            (+ (* 263.505074721 x) (+ (* 43.3400022514 (* x x)) (* x (* x x))))
            313.399215894)
           x)
          47.066876606)))
       t_0))))

C

double code(double x, double y, double z) {
	double t_0 = ((y / (x * x)) + (4.16438922228 * x)) - 110.1139242984811;
	double tmp;
	if (x < -3.326128725870005e+62) {
		tmp = t_0;
	} else if (x < 9.429991714554673e+55) {
		tmp = ((x - 2.0) / 1.0) * (((((((((x * 4.16438922228) + 78.6994924154) * x) + 137.519416416) * x) + y) * x) + z) / (((((263.505074721 * x) + ((43.3400022514 * (x * x)) + (x * (x * x)))) + 313.399215894) * x) + 47.066876606));
	} else {
		tmp = t_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) :: t_0
    real(8) :: tmp
    t_0 = ((y / (x * x)) + (4.16438922228d0 * x)) - 110.1139242984811d0
    if (x < (-3.326128725870005d+62)) then
        tmp = t_0
    else if (x < 9.429991714554673d+55) then
        tmp = ((x - 2.0d0) / 1.0d0) * (((((((((x * 4.16438922228d0) + 78.6994924154d0) * x) + 137.519416416d0) * x) + y) * x) + z) / (((((263.505074721d0 * x) + ((43.3400022514d0 * (x * x)) + (x * (x * x)))) + 313.399215894d0) * x) + 47.066876606d0))
    else
        tmp = t_0
    end if
    code = tmp
end function

Java

public static double code(double x, double y, double z) {
	double t_0 = ((y / (x * x)) + (4.16438922228 * x)) - 110.1139242984811;
	double tmp;
	if (x < -3.326128725870005e+62) {
		tmp = t_0;
	} else if (x < 9.429991714554673e+55) {
		tmp = ((x - 2.0) / 1.0) * (((((((((x * 4.16438922228) + 78.6994924154) * x) + 137.519416416) * x) + y) * x) + z) / (((((263.505074721 * x) + ((43.3400022514 * (x * x)) + (x * (x * x)))) + 313.399215894) * x) + 47.066876606));
	} else {
		tmp = t_0;
	}
	return tmp;
}

Python

def code(x, y, z):
	t_0 = ((y / (x * x)) + (4.16438922228 * x)) - 110.1139242984811
	tmp = 0
	if x < -3.326128725870005e+62:
		tmp = t_0
	elif x < 9.429991714554673e+55:
		tmp = ((x - 2.0) / 1.0) * (((((((((x * 4.16438922228) + 78.6994924154) * x) + 137.519416416) * x) + y) * x) + z) / (((((263.505074721 * x) + ((43.3400022514 * (x * x)) + (x * (x * x)))) + 313.399215894) * x) + 47.066876606))
	else:
		tmp = t_0
	return tmp

Julia

function code(x, y, z)
	t_0 = Float64(Float64(Float64(y / Float64(x * x)) + Float64(4.16438922228 * x)) - 110.1139242984811)
	tmp = 0.0
	if (x < -3.326128725870005e+62)
		tmp = t_0;
	elseif (x < 9.429991714554673e+55)
		tmp = Float64(Float64(Float64(x - 2.0) / 1.0) * Float64(Float64(Float64(Float64(Float64(Float64(Float64(Float64(Float64(x * 4.16438922228) + 78.6994924154) * x) + 137.519416416) * x) + y) * x) + z) / Float64(Float64(Float64(Float64(Float64(263.505074721 * x) + Float64(Float64(43.3400022514 * Float64(x * x)) + Float64(x * Float64(x * x)))) + 313.399215894) * x) + 47.066876606)));
	else
		tmp = t_0;
	end
	return tmp
end

MATLAB

function tmp_2 = code(x, y, z)
	t_0 = ((y / (x * x)) + (4.16438922228 * x)) - 110.1139242984811;
	tmp = 0.0;
	if (x < -3.326128725870005e+62)
		tmp = t_0;
	elseif (x < 9.429991714554673e+55)
		tmp = ((x - 2.0) / 1.0) * (((((((((x * 4.16438922228) + 78.6994924154) * x) + 137.519416416) * x) + y) * x) + z) / (((((263.505074721 * x) + ((43.3400022514 * (x * x)) + (x * (x * x)))) + 313.399215894) * x) + 47.066876606));
	else
		tmp = t_0;
	end
	tmp_2 = tmp;
end

Wolfram

code[x_, y_, z_] := Block[{t$95$0 = N[(N[(N[(y / N[(x * x), $MachinePrecision]), $MachinePrecision] + N[(4.16438922228 * x), $MachinePrecision]), $MachinePrecision] - 110.1139242984811), $MachinePrecision]}, If[Less[x, -3.326128725870005e+62], t$95$0, If[Less[x, 9.429991714554673e+55], N[(N[(N[(x - 2.0), $MachinePrecision] / 1.0), $MachinePrecision] * N[(N[(N[(N[(N[(N[(N[(N[(N[(x * 4.16438922228), $MachinePrecision] + 78.6994924154), $MachinePrecision] * x), $MachinePrecision] + 137.519416416), $MachinePrecision] * x), $MachinePrecision] + y), $MachinePrecision] * x), $MachinePrecision] + z), $MachinePrecision] / N[(N[(N[(N[(N[(263.505074721 * x), $MachinePrecision] + N[(N[(43.3400022514 * N[(x * x), $MachinePrecision]), $MachinePrecision] + N[(x * N[(x * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + 313.399215894), $MachinePrecision] * x), $MachinePrecision] + 47.066876606), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$0]]]

Reproduce

herbie shell --seed 2024254 
(FPCore (x y z)
  :name "Numeric.SpecFunctions:logGamma from math-functions-0.1.5.2, C"
  :precision binary64

  :alt
  (! :herbie-platform default (if (< x -332612872587000500000000000000000000000000000000000000000000000) (- (+ (/ y (* x x)) (* 104109730557/25000000000 x)) 1101139242984811/10000000000000) (if (< x 94299917145546730000000000000000000000000000000000000000) (* (/ (- x 2) 1) (/ (+ (* (+ (* (+ (* (+ (* x 104109730557/25000000000) 393497462077/5000000000) x) 4297481763/31250000) x) y) x) z) (+ (* (+ (+ (* 263505074721/1000000000 x) (+ (* 216700011257/5000000000 (* x x)) (* x (* x x)))) 156699607947/500000000) x) 23533438303/500000000))) (- (+ (/ y (* x x)) (* 104109730557/25000000000 x)) 1101139242984811/10000000000000))))

  (/ (* (- x 2.0) (+ (* (+ (* (+ (* (+ (* x 4.16438922228) 78.6994924154) x) 137.519416416) x) y) x) z)) (+ (* (+ (* (+ (* (+ x 43.3400022514) x) 263.505074721) x) 313.399215894) x) 47.066876606)))