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

Percentage Accurate: 58.7% → 98.7%

Time: 8.0s

Alternatives: 19

Speedup: 8.8×

Specification

Math

\[\begin{array}{l} \\ x + \frac{y \cdot \left(\left(\left(\left(z \cdot 3.13060547623 + 11.1667541262\right) \cdot z + t\right) \cdot z + a\right) \cdot z + b\right)}{\left(\left(\left(z + 15.234687407\right) \cdot z + 31.4690115749\right) \cdot z + 11.9400905721\right) \cdot z + 0.607771387771} \end{array} \]

FPCore

(FPCore (x y z t a b)
 :precision binary64
 (+
  x
  (/
   (*
    y
    (+ (* (+ (* (+ (* (+ (* z 3.13060547623) 11.1667541262) z) t) z) a) z) b))
   (+
    (* (+ (* (+ (* (+ z 15.234687407) z) 31.4690115749) z) 11.9400905721) z)
    0.607771387771))))

C

double code(double x, double y, double z, double t, double a, double b) {
	return x + ((y * ((((((((z * 3.13060547623) + 11.1667541262) * z) + t) * z) + a) * z) + b)) / (((((((z + 15.234687407) * z) + 31.4690115749) * z) + 11.9400905721) * z) + 0.607771387771));
}

Fortran

module fmin_fmax_functions
    implicit none
    private
    public fmax
    public fmin

    interface fmax
        module procedure fmax88
        module procedure fmax44
        module procedure fmax84
        module procedure fmax48
    end interface
    interface fmin
        module procedure fmin88
        module procedure fmin44
        module procedure fmin84
        module procedure fmin48
    end interface
contains
    real(8) function fmax88(x, y) result (res)
        real(8), intent (in) :: x
        real(8), intent (in) :: y
        res = merge(y, merge(x, max(x, y), y /= y), x /= x)
    end function
    real(4) function fmax44(x, y) result (res)
        real(4), intent (in) :: x
        real(4), intent (in) :: y
        res = merge(y, merge(x, max(x, y), y /= y), x /= x)
    end function
    real(8) function fmax84(x, y) result(res)
        real(8), intent (in) :: x
        real(4), intent (in) :: y
        res = merge(dble(y), merge(x, max(x, dble(y)), y /= y), x /= x)
    end function
    real(8) function fmax48(x, y) result(res)
        real(4), intent (in) :: x
        real(8), intent (in) :: y
        res = merge(y, merge(dble(x), max(dble(x), y), y /= y), x /= x)
    end function
    real(8) function fmin88(x, y) result (res)
        real(8), intent (in) :: x
        real(8), intent (in) :: y
        res = merge(y, merge(x, min(x, y), y /= y), x /= x)
    end function
    real(4) function fmin44(x, y) result (res)
        real(4), intent (in) :: x
        real(4), intent (in) :: y
        res = merge(y, merge(x, min(x, y), y /= y), x /= x)
    end function
    real(8) function fmin84(x, y) result(res)
        real(8), intent (in) :: x
        real(4), intent (in) :: y
        res = merge(dble(y), merge(x, min(x, dble(y)), y /= y), x /= x)
    end function
    real(8) function fmin48(x, y) result(res)
        real(4), intent (in) :: x
        real(8), intent (in) :: y
        res = merge(y, merge(dble(x), min(dble(x), y), y /= y), x /= x)
    end function
end module

real(8) function code(x, y, z, t, a, b)
use fmin_fmax_functions
    real(8), intent (in) :: x
    real(8), intent (in) :: y
    real(8), intent (in) :: z
    real(8), intent (in) :: t
    real(8), intent (in) :: a
    real(8), intent (in) :: b
    code = x + ((y * ((((((((z * 3.13060547623d0) + 11.1667541262d0) * z) + t) * z) + a) * z) + b)) / (((((((z + 15.234687407d0) * z) + 31.4690115749d0) * z) + 11.9400905721d0) * z) + 0.607771387771d0))
end function

Java

public static double code(double x, double y, double z, double t, double a, double b) {
	return x + ((y * ((((((((z * 3.13060547623) + 11.1667541262) * z) + t) * z) + a) * z) + b)) / (((((((z + 15.234687407) * z) + 31.4690115749) * z) + 11.9400905721) * z) + 0.607771387771));
}

Python

def code(x, y, z, t, a, b):
	return x + ((y * ((((((((z * 3.13060547623) + 11.1667541262) * z) + t) * z) + a) * z) + b)) / (((((((z + 15.234687407) * z) + 31.4690115749) * z) + 11.9400905721) * z) + 0.607771387771))

Julia

function code(x, y, z, t, a, b)
	return Float64(x + Float64(Float64(y * Float64(Float64(Float64(Float64(Float64(Float64(Float64(Float64(z * 3.13060547623) + 11.1667541262) * z) + t) * z) + a) * z) + b)) / Float64(Float64(Float64(Float64(Float64(Float64(Float64(z + 15.234687407) * z) + 31.4690115749) * z) + 11.9400905721) * z) + 0.607771387771)))
end

MATLAB

function tmp = code(x, y, z, t, a, b)
	tmp = x + ((y * ((((((((z * 3.13060547623) + 11.1667541262) * z) + t) * z) + a) * z) + b)) / (((((((z + 15.234687407) * z) + 31.4690115749) * z) + 11.9400905721) * z) + 0.607771387771));
end

Wolfram

code[x_, y_, z_, t_, a_, b_] := N[(x + N[(N[(y * N[(N[(N[(N[(N[(N[(N[(N[(z * 3.13060547623), $MachinePrecision] + 11.1667541262), $MachinePrecision] * z), $MachinePrecision] + t), $MachinePrecision] * z), $MachinePrecision] + a), $MachinePrecision] * z), $MachinePrecision] + b), $MachinePrecision]), $MachinePrecision] / N[(N[(N[(N[(N[(N[(N[(z + 15.234687407), $MachinePrecision] * z), $MachinePrecision] + 31.4690115749), $MachinePrecision] * z), $MachinePrecision] + 11.9400905721), $MachinePrecision] * z), $MachinePrecision] + 0.607771387771), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]

Initial Program: 58.7% accurate, 1.0× speedup

Math

\[\begin{array}{l} \\ x + \frac{y \cdot \left(\left(\left(\left(z \cdot 3.13060547623 + 11.1667541262\right) \cdot z + t\right) \cdot z + a\right) \cdot z + b\right)}{\left(\left(\left(z + 15.234687407\right) \cdot z + 31.4690115749\right) \cdot z + 11.9400905721\right) \cdot z + 0.607771387771} \end{array} \]

FPCore

(FPCore (x y z t a b)
 :precision binary64
 (+
  x
  (/
   (*
    y
    (+ (* (+ (* (+ (* (+ (* z 3.13060547623) 11.1667541262) z) t) z) a) z) b))
   (+
    (* (+ (* (+ (* (+ z 15.234687407) z) 31.4690115749) z) 11.9400905721) z)
    0.607771387771))))

C

double code(double x, double y, double z, double t, double a, double b) {
	return x + ((y * ((((((((z * 3.13060547623) + 11.1667541262) * z) + t) * z) + a) * z) + b)) / (((((((z + 15.234687407) * z) + 31.4690115749) * z) + 11.9400905721) * z) + 0.607771387771));
}

Fortran

module fmin_fmax_functions
    implicit none
    private
    public fmax
    public fmin

    interface fmax
        module procedure fmax88
        module procedure fmax44
        module procedure fmax84
        module procedure fmax48
    end interface
    interface fmin
        module procedure fmin88
        module procedure fmin44
        module procedure fmin84
        module procedure fmin48
    end interface
contains
    real(8) function fmax88(x, y) result (res)
        real(8), intent (in) :: x
        real(8), intent (in) :: y
        res = merge(y, merge(x, max(x, y), y /= y), x /= x)
    end function
    real(4) function fmax44(x, y) result (res)
        real(4), intent (in) :: x
        real(4), intent (in) :: y
        res = merge(y, merge(x, max(x, y), y /= y), x /= x)
    end function
    real(8) function fmax84(x, y) result(res)
        real(8), intent (in) :: x
        real(4), intent (in) :: y
        res = merge(dble(y), merge(x, max(x, dble(y)), y /= y), x /= x)
    end function
    real(8) function fmax48(x, y) result(res)
        real(4), intent (in) :: x
        real(8), intent (in) :: y
        res = merge(y, merge(dble(x), max(dble(x), y), y /= y), x /= x)
    end function
    real(8) function fmin88(x, y) result (res)
        real(8), intent (in) :: x
        real(8), intent (in) :: y
        res = merge(y, merge(x, min(x, y), y /= y), x /= x)
    end function
    real(4) function fmin44(x, y) result (res)
        real(4), intent (in) :: x
        real(4), intent (in) :: y
        res = merge(y, merge(x, min(x, y), y /= y), x /= x)
    end function
    real(8) function fmin84(x, y) result(res)
        real(8), intent (in) :: x
        real(4), intent (in) :: y
        res = merge(dble(y), merge(x, min(x, dble(y)), y /= y), x /= x)
    end function
    real(8) function fmin48(x, y) result(res)
        real(4), intent (in) :: x
        real(8), intent (in) :: y
        res = merge(y, merge(dble(x), min(dble(x), y), y /= y), x /= x)
    end function
end module

real(8) function code(x, y, z, t, a, b)
use fmin_fmax_functions
    real(8), intent (in) :: x
    real(8), intent (in) :: y
    real(8), intent (in) :: z
    real(8), intent (in) :: t
    real(8), intent (in) :: a
    real(8), intent (in) :: b
    code = x + ((y * ((((((((z * 3.13060547623d0) + 11.1667541262d0) * z) + t) * z) + a) * z) + b)) / (((((((z + 15.234687407d0) * z) + 31.4690115749d0) * z) + 11.9400905721d0) * z) + 0.607771387771d0))
end function

Java

public static double code(double x, double y, double z, double t, double a, double b) {
	return x + ((y * ((((((((z * 3.13060547623) + 11.1667541262) * z) + t) * z) + a) * z) + b)) / (((((((z + 15.234687407) * z) + 31.4690115749) * z) + 11.9400905721) * z) + 0.607771387771));
}

Python

def code(x, y, z, t, a, b):
	return x + ((y * ((((((((z * 3.13060547623) + 11.1667541262) * z) + t) * z) + a) * z) + b)) / (((((((z + 15.234687407) * z) + 31.4690115749) * z) + 11.9400905721) * z) + 0.607771387771))

Julia

function code(x, y, z, t, a, b)
	return Float64(x + Float64(Float64(y * Float64(Float64(Float64(Float64(Float64(Float64(Float64(Float64(z * 3.13060547623) + 11.1667541262) * z) + t) * z) + a) * z) + b)) / Float64(Float64(Float64(Float64(Float64(Float64(Float64(z + 15.234687407) * z) + 31.4690115749) * z) + 11.9400905721) * z) + 0.607771387771)))
end

MATLAB

function tmp = code(x, y, z, t, a, b)
	tmp = x + ((y * ((((((((z * 3.13060547623) + 11.1667541262) * z) + t) * z) + a) * z) + b)) / (((((((z + 15.234687407) * z) + 31.4690115749) * z) + 11.9400905721) * z) + 0.607771387771));
end

Wolfram

code[x_, y_, z_, t_, a_, b_] := N[(x + N[(N[(y * N[(N[(N[(N[(N[(N[(N[(N[(z * 3.13060547623), $MachinePrecision] + 11.1667541262), $MachinePrecision] * z), $MachinePrecision] + t), $MachinePrecision] * z), $MachinePrecision] + a), $MachinePrecision] * z), $MachinePrecision] + b), $MachinePrecision]), $MachinePrecision] / N[(N[(N[(N[(N[(N[(N[(z + 15.234687407), $MachinePrecision] * z), $MachinePrecision] + 31.4690115749), $MachinePrecision] * z), $MachinePrecision] + 11.9400905721), $MachinePrecision] * z), $MachinePrecision] + 0.607771387771), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]

Alternative 1: 98.7% accurate, 0.4× speedup

Math

\[\begin{array}{l} \\ \begin{array}{l} t_1 := \mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(z - -15.234687407, z, 31.4690115749\right), z, 11.9400905721\right), z, 0.607771387771\right)\\ \mathbf{if}\;x + \frac{y \cdot \left(\left(\left(\left(z \cdot 3.13060547623 + 11.1667541262\right) \cdot z + t\right) \cdot z + a\right) \cdot z + b\right)}{\left(\left(\left(z + 15.234687407\right) \cdot z + 31.4690115749\right) \cdot z + 11.9400905721\right) \cdot z + 0.607771387771} \leq \infty:\\ \;\;\;\;\mathsf{fma}\left(y, z \cdot \frac{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(3.13060547623, z, 11.1667541262\right), z, t\right), z, a\right)}{t\_1}, \mathsf{fma}\left(b, \frac{y}{t\_1}, x\right)\right)\\ \mathbf{else}:\\ \;\;\;\;\mathsf{fma}\left(3.13060547623 + -1 \cdot \frac{36.52704169880642 + -1 \cdot \frac{457.9610022158428 + t}{z}}{z}, y, x\right)\\ \end{array} \end{array} \]

FPCore

(FPCore (x y z t a b)
 :precision binary64
 (let* ((t_1
         (fma
          (fma (fma (- z -15.234687407) z 31.4690115749) z 11.9400905721)
          z
          0.607771387771)))
   (if (<=
        (+
         x
         (/
          (*
           y
           (+
            (* (+ (* (+ (* (+ (* z 3.13060547623) 11.1667541262) z) t) z) a) z)
            b))
          (+
           (*
            (+ (* (+ (* (+ z 15.234687407) z) 31.4690115749) z) 11.9400905721)
            z)
           0.607771387771)))
        INFINITY)
     (fma
      y
      (* z (/ (fma (fma (fma 3.13060547623 z 11.1667541262) z t) z a) t_1))
      (fma b (/ y t_1) x))
     (fma
      (+
       3.13060547623
       (*
        -1.0
        (/ (+ 36.52704169880642 (* -1.0 (/ (+ 457.9610022158428 t) z))) z)))
      y
      x))))

C

double code(double x, double y, double z, double t, double a, double b) {
	double t_1 = fma(fma(fma((z - -15.234687407), z, 31.4690115749), z, 11.9400905721), z, 0.607771387771);
	double tmp;
	if ((x + ((y * ((((((((z * 3.13060547623) + 11.1667541262) * z) + t) * z) + a) * z) + b)) / (((((((z + 15.234687407) * z) + 31.4690115749) * z) + 11.9400905721) * z) + 0.607771387771))) <= ((double) INFINITY)) {
		tmp = fma(y, (z * (fma(fma(fma(3.13060547623, z, 11.1667541262), z, t), z, a) / t_1)), fma(b, (y / t_1), x));
	} else {
		tmp = fma((3.13060547623 + (-1.0 * ((36.52704169880642 + (-1.0 * ((457.9610022158428 + t) / z))) / z))), y, x);
	}
	return tmp;
}

Julia

function code(x, y, z, t, a, b)
	t_1 = fma(fma(fma(Float64(z - -15.234687407), z, 31.4690115749), z, 11.9400905721), z, 0.607771387771)
	tmp = 0.0
	if (Float64(x + Float64(Float64(y * Float64(Float64(Float64(Float64(Float64(Float64(Float64(Float64(z * 3.13060547623) + 11.1667541262) * z) + t) * z) + a) * z) + b)) / Float64(Float64(Float64(Float64(Float64(Float64(Float64(z + 15.234687407) * z) + 31.4690115749) * z) + 11.9400905721) * z) + 0.607771387771))) <= Inf)
		tmp = fma(y, Float64(z * Float64(fma(fma(fma(3.13060547623, z, 11.1667541262), z, t), z, a) / t_1)), fma(b, Float64(y / t_1), x));
	else
		tmp = fma(Float64(3.13060547623 + Float64(-1.0 * Float64(Float64(36.52704169880642 + Float64(-1.0 * Float64(Float64(457.9610022158428 + t) / z))) / z))), y, x);
	end
	return tmp
end

Wolfram

code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(N[(N[(N[(z - -15.234687407), $MachinePrecision] * z + 31.4690115749), $MachinePrecision] * z + 11.9400905721), $MachinePrecision] * z + 0.607771387771), $MachinePrecision]}, If[LessEqual[N[(x + N[(N[(y * N[(N[(N[(N[(N[(N[(N[(N[(z * 3.13060547623), $MachinePrecision] + 11.1667541262), $MachinePrecision] * z), $MachinePrecision] + t), $MachinePrecision] * z), $MachinePrecision] + a), $MachinePrecision] * z), $MachinePrecision] + b), $MachinePrecision]), $MachinePrecision] / N[(N[(N[(N[(N[(N[(N[(z + 15.234687407), $MachinePrecision] * z), $MachinePrecision] + 31.4690115749), $MachinePrecision] * z), $MachinePrecision] + 11.9400905721), $MachinePrecision] * z), $MachinePrecision] + 0.607771387771), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], Infinity], N[(y * N[(z * N[(N[(N[(N[(3.13060547623 * z + 11.1667541262), $MachinePrecision] * z + t), $MachinePrecision] * z + a), $MachinePrecision] / t$95$1), $MachinePrecision]), $MachinePrecision] + N[(b * N[(y / t$95$1), $MachinePrecision] + x), $MachinePrecision]), $MachinePrecision], N[(N[(3.13060547623 + N[(-1.0 * N[(N[(36.52704169880642 + N[(-1.0 * N[(N[(457.9610022158428 + t), $MachinePrecision] / z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * y + x), $MachinePrecision]]]

Derivation

1. Split input into 2 regimes

2. if (+.f64 x (/.f64 (*.f64 y (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 z #s(literal 313060547623/100000000000 binary64)) #s(literal 55833770631/5000000000 binary64)) z) t) z) a) z) b)) (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 z #s(literal 15234687407/1000000000 binary64)) z) #s(literal 314690115749/10000000000 binary64)) z) #s(literal 119400905721/10000000000 binary64)) z) #s(literal 607771387771/1000000000000 binary64)))) < +inf.0

   1. Initial program 58.7%

      \[x + \frac{y \cdot \left(\left(\left(\left(z \cdot 3.13060547623 + 11.1667541262\right) \cdot z + t\right) \cdot z + a\right) \cdot z + b\right)}{\left(\left(\left(z + 15.234687407\right) \cdot z + 31.4690115749\right) \cdot z + 11.9400905721\right) \cdot z + 0.607771387771} \]

   3. Applied rewrites 63.3%

      \[\leadsto \color{blue}{\mathsf{fma}\left(y, z \cdot \frac{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(3.13060547623, z, 11.1667541262\right), z, t\right), z, a\right)}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(z - -15.234687407, z, 31.4690115749\right), z, 11.9400905721\right), z, 0.607771387771\right)}, \mathsf{fma}\left(b, \frac{y}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(z - -15.234687407, z, 31.4690115749\right), z, 11.9400905721\right), z, 0.607771387771\right)}, x\right)\right)} \]

   if +inf.0 < (+.f64 x (/.f64 (*.f64 y (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 z #s(literal 313060547623/100000000000 binary64)) #s(literal 55833770631/5000000000 binary64)) z) t) z) a) z) b)) (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 z #s(literal 15234687407/1000000000 binary64)) z) #s(literal 314690115749/10000000000 binary64)) z) #s(literal 119400905721/10000000000 binary64)) z) #s(literal 607771387771/1000000000000 binary64))))

   1. Initial program 58.7%

      \[x + \frac{y \cdot \left(\left(\left(\left(z \cdot 3.13060547623 + 11.1667541262\right) \cdot z + t\right) \cdot z + a\right) \cdot z + b\right)}{\left(\left(\left(z + 15.234687407\right) \cdot z + 31.4690115749\right) \cdot z + 11.9400905721\right) \cdot z + 0.607771387771} \]

   3. Applied rewrites 60.9%

      \[\leadsto \color{blue}{\mathsf{fma}\left(\frac{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(3.13060547623, z, 11.1667541262\right), z, t\right), z, a\right), z, b\right)}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(z - -15.234687407, z, 31.4690115749\right), z, 11.9400905721\right), z, 0.607771387771\right)}, y, x\right)} \]

   4. Taylor expanded in z around -inf

      \[\leadsto \mathsf{fma}\left(\color{blue}{\frac{313060547623}{100000000000} + -1 \cdot \frac{\frac{3652704169880641883561}{100000000000000000000} + -1 \cdot \frac{\frac{45796100221584283915100827016327}{100000000000000000000000000000} + t}{z}}{z}}, y, x\right) \]
   5. Step-by-step derivation

      1. lower-+.f64 N/A

         \[\leadsto \mathsf{fma}\left(\frac{313060547623}{100000000000} + \color{blue}{-1 \cdot \frac{\frac{3652704169880641883561}{100000000000000000000} + -1 \cdot \frac{\frac{45796100221584283915100827016327}{100000000000000000000000000000} + t}{z}}{z}}, y, x\right) \]

      2. lower-*.f64 N/A

         \[\leadsto \mathsf{fma}\left(\frac{313060547623}{100000000000} + -1 \cdot \color{blue}{\frac{\frac{3652704169880641883561}{100000000000000000000} + -1 \cdot \frac{\frac{45796100221584283915100827016327}{100000000000000000000000000000} + t}{z}}{z}}, y, x\right) \]

      3. lower-/.f64 N/A

         \[\leadsto \mathsf{fma}\left(\frac{313060547623}{100000000000} + -1 \cdot \frac{\frac{3652704169880641883561}{100000000000000000000} + -1 \cdot \frac{\frac{45796100221584283915100827016327}{100000000000000000000000000000} + t}{z}}{\color{blue}{z}}, y, x\right) \]

      4. lower-+.f64 N/A

         \[\leadsto \mathsf{fma}\left(\frac{313060547623}{100000000000} + -1 \cdot \frac{\frac{3652704169880641883561}{100000000000000000000} + -1 \cdot \frac{\frac{45796100221584283915100827016327}{100000000000000000000000000000} + t}{z}}{z}, y, x\right) \]

      5. lower-*.f64 N/A

         \[\leadsto \mathsf{fma}\left(\frac{313060547623}{100000000000} + -1 \cdot \frac{\frac{3652704169880641883561}{100000000000000000000} + -1 \cdot \frac{\frac{45796100221584283915100827016327}{100000000000000000000000000000} + t}{z}}{z}, y, x\right) \]

      6. lower-/.f64 N/A

         \[\leadsto \mathsf{fma}\left(\frac{313060547623}{100000000000} + -1 \cdot \frac{\frac{3652704169880641883561}{100000000000000000000} + -1 \cdot \frac{\frac{45796100221584283915100827016327}{100000000000000000000000000000} + t}{z}}{z}, y, x\right) \]

      7. lower-+.f64 55.8

         \[\leadsto \mathsf{fma}\left(3.13060547623 + -1 \cdot \frac{36.52704169880642 + -1 \cdot \frac{457.9610022158428 + t}{z}}{z}, y, x\right) \]

   6. Applied rewrites 55.8%

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

Alternative 2: 98.0% accurate, 0.5× speedup

Math

\[\begin{array}{l} \\ \begin{array}{l} \mathbf{if}\;x + \frac{y \cdot \left(\left(\left(\left(z \cdot 3.13060547623 + 11.1667541262\right) \cdot z + t\right) \cdot z + a\right) \cdot z + b\right)}{\left(\left(\left(z + 15.234687407\right) \cdot z + 31.4690115749\right) \cdot z + 11.9400905721\right) \cdot z + 0.607771387771} \leq \infty:\\ \;\;\;\;\mathsf{fma}\left(\frac{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(3.13060547623, z, 11.1667541262\right), z, t\right), z, a\right), z, b\right)}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(z - -15.234687407, z, 31.4690115749\right), z, 11.9400905721\right), z, 0.607771387771\right)}, y, x\right)\\ \mathbf{else}:\\ \;\;\;\;\mathsf{fma}\left(3.13060547623 + -1 \cdot \frac{36.52704169880642 + -1 \cdot \frac{457.9610022158428 + t}{z}}{z}, y, x\right)\\ \end{array} \end{array} \]

FPCore

(FPCore (x y z t a b)
 :precision binary64
 (if (<=
      (+
       x
       (/
        (*
         y
         (+
          (* (+ (* (+ (* (+ (* z 3.13060547623) 11.1667541262) z) t) z) a) z)
          b))
        (+
         (*
          (+ (* (+ (* (+ z 15.234687407) z) 31.4690115749) z) 11.9400905721)
          z)
         0.607771387771)))
      INFINITY)
   (fma
    (/
     (fma (fma (fma (fma 3.13060547623 z 11.1667541262) z t) z a) z b)
     (fma
      (fma (fma (- z -15.234687407) z 31.4690115749) z 11.9400905721)
      z
      0.607771387771))
    y
    x)
   (fma
    (+
     3.13060547623
     (*
      -1.0
      (/ (+ 36.52704169880642 (* -1.0 (/ (+ 457.9610022158428 t) z))) z)))
    y
    x)))

C

double code(double x, double y, double z, double t, double a, double b) {
	double tmp;
	if ((x + ((y * ((((((((z * 3.13060547623) + 11.1667541262) * z) + t) * z) + a) * z) + b)) / (((((((z + 15.234687407) * z) + 31.4690115749) * z) + 11.9400905721) * z) + 0.607771387771))) <= ((double) INFINITY)) {
		tmp = fma((fma(fma(fma(fma(3.13060547623, z, 11.1667541262), z, t), z, a), z, b) / fma(fma(fma((z - -15.234687407), z, 31.4690115749), z, 11.9400905721), z, 0.607771387771)), y, x);
	} else {
		tmp = fma((3.13060547623 + (-1.0 * ((36.52704169880642 + (-1.0 * ((457.9610022158428 + t) / z))) / z))), y, x);
	}
	return tmp;
}

Julia

function code(x, y, z, t, a, b)
	tmp = 0.0
	if (Float64(x + Float64(Float64(y * Float64(Float64(Float64(Float64(Float64(Float64(Float64(Float64(z * 3.13060547623) + 11.1667541262) * z) + t) * z) + a) * z) + b)) / Float64(Float64(Float64(Float64(Float64(Float64(Float64(z + 15.234687407) * z) + 31.4690115749) * z) + 11.9400905721) * z) + 0.607771387771))) <= Inf)
		tmp = fma(Float64(fma(fma(fma(fma(3.13060547623, z, 11.1667541262), z, t), z, a), z, b) / fma(fma(fma(Float64(z - -15.234687407), z, 31.4690115749), z, 11.9400905721), z, 0.607771387771)), y, x);
	else
		tmp = fma(Float64(3.13060547623 + Float64(-1.0 * Float64(Float64(36.52704169880642 + Float64(-1.0 * Float64(Float64(457.9610022158428 + t) / z))) / z))), y, x);
	end
	return tmp
end

Wolfram

code[x_, y_, z_, t_, a_, b_] := If[LessEqual[N[(x + N[(N[(y * N[(N[(N[(N[(N[(N[(N[(N[(z * 3.13060547623), $MachinePrecision] + 11.1667541262), $MachinePrecision] * z), $MachinePrecision] + t), $MachinePrecision] * z), $MachinePrecision] + a), $MachinePrecision] * z), $MachinePrecision] + b), $MachinePrecision]), $MachinePrecision] / N[(N[(N[(N[(N[(N[(N[(z + 15.234687407), $MachinePrecision] * z), $MachinePrecision] + 31.4690115749), $MachinePrecision] * z), $MachinePrecision] + 11.9400905721), $MachinePrecision] * z), $MachinePrecision] + 0.607771387771), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], Infinity], N[(N[(N[(N[(N[(N[(3.13060547623 * z + 11.1667541262), $MachinePrecision] * z + t), $MachinePrecision] * z + a), $MachinePrecision] * z + b), $MachinePrecision] / N[(N[(N[(N[(z - -15.234687407), $MachinePrecision] * z + 31.4690115749), $MachinePrecision] * z + 11.9400905721), $MachinePrecision] * z + 0.607771387771), $MachinePrecision]), $MachinePrecision] * y + x), $MachinePrecision], N[(N[(3.13060547623 + N[(-1.0 * N[(N[(36.52704169880642 + N[(-1.0 * N[(N[(457.9610022158428 + t), $MachinePrecision] / z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * y + x), $MachinePrecision]]

Derivation

1. Split input into 2 regimes

2. if (+.f64 x (/.f64 (*.f64 y (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 z #s(literal 313060547623/100000000000 binary64)) #s(literal 55833770631/5000000000 binary64)) z) t) z) a) z) b)) (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 z #s(literal 15234687407/1000000000 binary64)) z) #s(literal 314690115749/10000000000 binary64)) z) #s(literal 119400905721/10000000000 binary64)) z) #s(literal 607771387771/1000000000000 binary64)))) < +inf.0

   1. Initial program 58.7%

      \[x + \frac{y \cdot \left(\left(\left(\left(z \cdot 3.13060547623 + 11.1667541262\right) \cdot z + t\right) \cdot z + a\right) \cdot z + b\right)}{\left(\left(\left(z + 15.234687407\right) \cdot z + 31.4690115749\right) \cdot z + 11.9400905721\right) \cdot z + 0.607771387771} \]

   3. Applied rewrites 60.9%

      \[\leadsto \color{blue}{\mathsf{fma}\left(\frac{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(3.13060547623, z, 11.1667541262\right), z, t\right), z, a\right), z, b\right)}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(z - -15.234687407, z, 31.4690115749\right), z, 11.9400905721\right), z, 0.607771387771\right)}, y, x\right)} \]

   if +inf.0 < (+.f64 x (/.f64 (*.f64 y (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 z #s(literal 313060547623/100000000000 binary64)) #s(literal 55833770631/5000000000 binary64)) z) t) z) a) z) b)) (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 z #s(literal 15234687407/1000000000 binary64)) z) #s(literal 314690115749/10000000000 binary64)) z) #s(literal 119400905721/10000000000 binary64)) z) #s(literal 607771387771/1000000000000 binary64))))

   1. Initial program 58.7%

      \[x + \frac{y \cdot \left(\left(\left(\left(z \cdot 3.13060547623 + 11.1667541262\right) \cdot z + t\right) \cdot z + a\right) \cdot z + b\right)}{\left(\left(\left(z + 15.234687407\right) \cdot z + 31.4690115749\right) \cdot z + 11.9400905721\right) \cdot z + 0.607771387771} \]

   3. Applied rewrites 60.9%

      \[\leadsto \color{blue}{\mathsf{fma}\left(\frac{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(3.13060547623, z, 11.1667541262\right), z, t\right), z, a\right), z, b\right)}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(z - -15.234687407, z, 31.4690115749\right), z, 11.9400905721\right), z, 0.607771387771\right)}, y, x\right)} \]

   4. Taylor expanded in z around -inf

      \[\leadsto \mathsf{fma}\left(\color{blue}{\frac{313060547623}{100000000000} + -1 \cdot \frac{\frac{3652704169880641883561}{100000000000000000000} + -1 \cdot \frac{\frac{45796100221584283915100827016327}{100000000000000000000000000000} + t}{z}}{z}}, y, x\right) \]
   5. Step-by-step derivation

      1. lower-+.f64 N/A

         \[\leadsto \mathsf{fma}\left(\frac{313060547623}{100000000000} + \color{blue}{-1 \cdot \frac{\frac{3652704169880641883561}{100000000000000000000} + -1 \cdot \frac{\frac{45796100221584283915100827016327}{100000000000000000000000000000} + t}{z}}{z}}, y, x\right) \]

      2. lower-*.f64 N/A

         \[\leadsto \mathsf{fma}\left(\frac{313060547623}{100000000000} + -1 \cdot \color{blue}{\frac{\frac{3652704169880641883561}{100000000000000000000} + -1 \cdot \frac{\frac{45796100221584283915100827016327}{100000000000000000000000000000} + t}{z}}{z}}, y, x\right) \]

      3. lower-/.f64 N/A

         \[\leadsto \mathsf{fma}\left(\frac{313060547623}{100000000000} + -1 \cdot \frac{\frac{3652704169880641883561}{100000000000000000000} + -1 \cdot \frac{\frac{45796100221584283915100827016327}{100000000000000000000000000000} + t}{z}}{\color{blue}{z}}, y, x\right) \]

      4. lower-+.f64 N/A

         \[\leadsto \mathsf{fma}\left(\frac{313060547623}{100000000000} + -1 \cdot \frac{\frac{3652704169880641883561}{100000000000000000000} + -1 \cdot \frac{\frac{45796100221584283915100827016327}{100000000000000000000000000000} + t}{z}}{z}, y, x\right) \]

      5. lower-*.f64 N/A

         \[\leadsto \mathsf{fma}\left(\frac{313060547623}{100000000000} + -1 \cdot \frac{\frac{3652704169880641883561}{100000000000000000000} + -1 \cdot \frac{\frac{45796100221584283915100827016327}{100000000000000000000000000000} + t}{z}}{z}, y, x\right) \]

      6. lower-/.f64 N/A

         \[\leadsto \mathsf{fma}\left(\frac{313060547623}{100000000000} + -1 \cdot \frac{\frac{3652704169880641883561}{100000000000000000000} + -1 \cdot \frac{\frac{45796100221584283915100827016327}{100000000000000000000000000000} + t}{z}}{z}, y, x\right) \]

      7. lower-+.f64 55.8

         \[\leadsto \mathsf{fma}\left(3.13060547623 + -1 \cdot \frac{36.52704169880642 + -1 \cdot \frac{457.9610022158428 + t}{z}}{z}, y, x\right) \]

   6. Applied rewrites 55.8%

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

Alternative 3: 97.6% accurate, 0.5× speedup

Math

\[\begin{array}{l} \\ \begin{array}{l} \mathbf{if}\;x + \frac{y \cdot \left(\left(\left(\left(z \cdot 3.13060547623 + 11.1667541262\right) \cdot z + t\right) \cdot z + a\right) \cdot z + b\right)}{\left(\left(\left(z + 15.234687407\right) \cdot z + 31.4690115749\right) \cdot z + 11.9400905721\right) \cdot z + 0.607771387771} \leq \infty:\\ \;\;\;\;\mathsf{fma}\left(\frac{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(3.13060547623 \cdot z, z, t\right), z, a\right), z, b\right)}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(z - -15.234687407, z, 31.4690115749\right), z, 11.9400905721\right), z, 0.607771387771\right)}, y, x\right)\\ \mathbf{else}:\\ \;\;\;\;\mathsf{fma}\left(3.13060547623 + -1 \cdot \frac{36.52704169880642 + -1 \cdot \frac{457.9610022158428 + t}{z}}{z}, y, x\right)\\ \end{array} \end{array} \]

FPCore

(FPCore (x y z t a b)
 :precision binary64
 (if (<=
      (+
       x
       (/
        (*
         y
         (+
          (* (+ (* (+ (* (+ (* z 3.13060547623) 11.1667541262) z) t) z) a) z)
          b))
        (+
         (*
          (+ (* (+ (* (+ z 15.234687407) z) 31.4690115749) z) 11.9400905721)
          z)
         0.607771387771)))
      INFINITY)
   (fma
    (/
     (fma (fma (fma (* 3.13060547623 z) z t) z a) z b)
     (fma
      (fma (fma (- z -15.234687407) z 31.4690115749) z 11.9400905721)
      z
      0.607771387771))
    y
    x)
   (fma
    (+
     3.13060547623
     (*
      -1.0
      (/ (+ 36.52704169880642 (* -1.0 (/ (+ 457.9610022158428 t) z))) z)))
    y
    x)))

C

double code(double x, double y, double z, double t, double a, double b) {
	double tmp;
	if ((x + ((y * ((((((((z * 3.13060547623) + 11.1667541262) * z) + t) * z) + a) * z) + b)) / (((((((z + 15.234687407) * z) + 31.4690115749) * z) + 11.9400905721) * z) + 0.607771387771))) <= ((double) INFINITY)) {
		tmp = fma((fma(fma(fma((3.13060547623 * z), z, t), z, a), z, b) / fma(fma(fma((z - -15.234687407), z, 31.4690115749), z, 11.9400905721), z, 0.607771387771)), y, x);
	} else {
		tmp = fma((3.13060547623 + (-1.0 * ((36.52704169880642 + (-1.0 * ((457.9610022158428 + t) / z))) / z))), y, x);
	}
	return tmp;
}

Julia

function code(x, y, z, t, a, b)
	tmp = 0.0
	if (Float64(x + Float64(Float64(y * Float64(Float64(Float64(Float64(Float64(Float64(Float64(Float64(z * 3.13060547623) + 11.1667541262) * z) + t) * z) + a) * z) + b)) / Float64(Float64(Float64(Float64(Float64(Float64(Float64(z + 15.234687407) * z) + 31.4690115749) * z) + 11.9400905721) * z) + 0.607771387771))) <= Inf)
		tmp = fma(Float64(fma(fma(fma(Float64(3.13060547623 * z), z, t), z, a), z, b) / fma(fma(fma(Float64(z - -15.234687407), z, 31.4690115749), z, 11.9400905721), z, 0.607771387771)), y, x);
	else
		tmp = fma(Float64(3.13060547623 + Float64(-1.0 * Float64(Float64(36.52704169880642 + Float64(-1.0 * Float64(Float64(457.9610022158428 + t) / z))) / z))), y, x);
	end
	return tmp
end

Wolfram

code[x_, y_, z_, t_, a_, b_] := If[LessEqual[N[(x + N[(N[(y * N[(N[(N[(N[(N[(N[(N[(N[(z * 3.13060547623), $MachinePrecision] + 11.1667541262), $MachinePrecision] * z), $MachinePrecision] + t), $MachinePrecision] * z), $MachinePrecision] + a), $MachinePrecision] * z), $MachinePrecision] + b), $MachinePrecision]), $MachinePrecision] / N[(N[(N[(N[(N[(N[(N[(z + 15.234687407), $MachinePrecision] * z), $MachinePrecision] + 31.4690115749), $MachinePrecision] * z), $MachinePrecision] + 11.9400905721), $MachinePrecision] * z), $MachinePrecision] + 0.607771387771), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], Infinity], N[(N[(N[(N[(N[(N[(3.13060547623 * z), $MachinePrecision] * z + t), $MachinePrecision] * z + a), $MachinePrecision] * z + b), $MachinePrecision] / N[(N[(N[(N[(z - -15.234687407), $MachinePrecision] * z + 31.4690115749), $MachinePrecision] * z + 11.9400905721), $MachinePrecision] * z + 0.607771387771), $MachinePrecision]), $MachinePrecision] * y + x), $MachinePrecision], N[(N[(3.13060547623 + N[(-1.0 * N[(N[(36.52704169880642 + N[(-1.0 * N[(N[(457.9610022158428 + t), $MachinePrecision] / z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * y + x), $MachinePrecision]]

Derivation

1. Split input into 2 regimes

2. if (+.f64 x (/.f64 (*.f64 y (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 z #s(literal 313060547623/100000000000 binary64)) #s(literal 55833770631/5000000000 binary64)) z) t) z) a) z) b)) (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 z #s(literal 15234687407/1000000000 binary64)) z) #s(literal 314690115749/10000000000 binary64)) z) #s(literal 119400905721/10000000000 binary64)) z) #s(literal 607771387771/1000000000000 binary64)))) < +inf.0

   1. Initial program 58.7%

      \[x + \frac{y \cdot \left(\left(\left(\left(z \cdot 3.13060547623 + 11.1667541262\right) \cdot z + t\right) \cdot z + a\right) \cdot z + b\right)}{\left(\left(\left(z + 15.234687407\right) \cdot z + 31.4690115749\right) \cdot z + 11.9400905721\right) \cdot z + 0.607771387771} \]

   3. Applied rewrites 60.9%

      \[\leadsto \color{blue}{\mathsf{fma}\left(\frac{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(3.13060547623, z, 11.1667541262\right), z, t\right), z, a\right), z, b\right)}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(z - -15.234687407, z, 31.4690115749\right), z, 11.9400905721\right), z, 0.607771387771\right)}, y, x\right)} \]

   4. Taylor expanded in z around inf

      \[\leadsto \mathsf{fma}\left(\frac{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(\color{blue}{\frac{313060547623}{100000000000} \cdot z}, z, t\right), z, a\right), z, b\right)}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(z - \frac{-15234687407}{1000000000}, z, \frac{314690115749}{10000000000}\right), z, \frac{119400905721}{10000000000}\right), z, \frac{607771387771}{1000000000000}\right)}, y, x\right) \]
   5. Step-by-step derivation

      1. lower-*.f64 60.6

         \[\leadsto \mathsf{fma}\left(\frac{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(3.13060547623 \cdot \color{blue}{z}, z, t\right), z, a\right), z, b\right)}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(z - -15.234687407, z, 31.4690115749\right), z, 11.9400905721\right), z, 0.607771387771\right)}, y, x\right) \]

   6. Applied rewrites 60.6%

      \[\leadsto \mathsf{fma}\left(\frac{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(\color{blue}{3.13060547623 \cdot z}, z, t\right), z, a\right), z, b\right)}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(z - -15.234687407, z, 31.4690115749\right), z, 11.9400905721\right), z, 0.607771387771\right)}, y, x\right) \]

   if +inf.0 < (+.f64 x (/.f64 (*.f64 y (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 z #s(literal 313060547623/100000000000 binary64)) #s(literal 55833770631/5000000000 binary64)) z) t) z) a) z) b)) (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 z #s(literal 15234687407/1000000000 binary64)) z) #s(literal 314690115749/10000000000 binary64)) z) #s(literal 119400905721/10000000000 binary64)) z) #s(literal 607771387771/1000000000000 binary64))))

   1. Initial program 58.7%

      \[x + \frac{y \cdot \left(\left(\left(\left(z \cdot 3.13060547623 + 11.1667541262\right) \cdot z + t\right) \cdot z + a\right) \cdot z + b\right)}{\left(\left(\left(z + 15.234687407\right) \cdot z + 31.4690115749\right) \cdot z + 11.9400905721\right) \cdot z + 0.607771387771} \]

   3. Applied rewrites 60.9%

      \[\leadsto \color{blue}{\mathsf{fma}\left(\frac{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(3.13060547623, z, 11.1667541262\right), z, t\right), z, a\right), z, b\right)}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(z - -15.234687407, z, 31.4690115749\right), z, 11.9400905721\right), z, 0.607771387771\right)}, y, x\right)} \]

   4. Taylor expanded in z around -inf

      \[\leadsto \mathsf{fma}\left(\color{blue}{\frac{313060547623}{100000000000} + -1 \cdot \frac{\frac{3652704169880641883561}{100000000000000000000} + -1 \cdot \frac{\frac{45796100221584283915100827016327}{100000000000000000000000000000} + t}{z}}{z}}, y, x\right) \]
   5. Step-by-step derivation

      1. lower-+.f64 N/A

         \[\leadsto \mathsf{fma}\left(\frac{313060547623}{100000000000} + \color{blue}{-1 \cdot \frac{\frac{3652704169880641883561}{100000000000000000000} + -1 \cdot \frac{\frac{45796100221584283915100827016327}{100000000000000000000000000000} + t}{z}}{z}}, y, x\right) \]

      2. lower-*.f64 N/A

         \[\leadsto \mathsf{fma}\left(\frac{313060547623}{100000000000} + -1 \cdot \color{blue}{\frac{\frac{3652704169880641883561}{100000000000000000000} + -1 \cdot \frac{\frac{45796100221584283915100827016327}{100000000000000000000000000000} + t}{z}}{z}}, y, x\right) \]

      3. lower-/.f64 N/A

         \[\leadsto \mathsf{fma}\left(\frac{313060547623}{100000000000} + -1 \cdot \frac{\frac{3652704169880641883561}{100000000000000000000} + -1 \cdot \frac{\frac{45796100221584283915100827016327}{100000000000000000000000000000} + t}{z}}{\color{blue}{z}}, y, x\right) \]

      4. lower-+.f64 N/A

         \[\leadsto \mathsf{fma}\left(\frac{313060547623}{100000000000} + -1 \cdot \frac{\frac{3652704169880641883561}{100000000000000000000} + -1 \cdot \frac{\frac{45796100221584283915100827016327}{100000000000000000000000000000} + t}{z}}{z}, y, x\right) \]

      5. lower-*.f64 N/A

         \[\leadsto \mathsf{fma}\left(\frac{313060547623}{100000000000} + -1 \cdot \frac{\frac{3652704169880641883561}{100000000000000000000} + -1 \cdot \frac{\frac{45796100221584283915100827016327}{100000000000000000000000000000} + t}{z}}{z}, y, x\right) \]

      6. lower-/.f64 N/A

         \[\leadsto \mathsf{fma}\left(\frac{313060547623}{100000000000} + -1 \cdot \frac{\frac{3652704169880641883561}{100000000000000000000} + -1 \cdot \frac{\frac{45796100221584283915100827016327}{100000000000000000000000000000} + t}{z}}{z}, y, x\right) \]

      7. lower-+.f64 55.8

         \[\leadsto \mathsf{fma}\left(3.13060547623 + -1 \cdot \frac{36.52704169880642 + -1 \cdot \frac{457.9610022158428 + t}{z}}{z}, y, x\right) \]

   6. Applied rewrites 55.8%

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

Alternative 4: 96.9% accurate, 1.2× speedup

Math

\[\begin{array}{l} \\ \begin{array}{l} t_1 := \mathsf{fma}\left(3.13060547623 + -1 \cdot \frac{36.52704169880642 + -1 \cdot \frac{457.9610022158428 + t}{z}}{z}, y, x\right)\\ \mathbf{if}\;z \leq -1.85 \cdot 10^{+36}:\\ \;\;\;\;t\_1\\ \mathbf{elif}\;z \leq 1.9 \cdot 10^{+20}:\\ \;\;\;\;x - \frac{\mathsf{fma}\left(\mathsf{fma}\left(t, z, a\right), z, b\right) \cdot y}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(-15.234687407 - z, z, -31.4690115749\right), z, -11.9400905721\right), z, -0.607771387771\right)}\\ \mathbf{else}:\\ \;\;\;\;t\_1\\ \end{array} \end{array} \]

FPCore

(FPCore (x y z t a b)
 :precision binary64
 (let* ((t_1
         (fma
          (+
           3.13060547623
           (*
            -1.0
            (/
             (+ 36.52704169880642 (* -1.0 (/ (+ 457.9610022158428 t) z)))
             z)))
          y
          x)))
   (if (<= z -1.85e+36)
     t_1
     (if (<= z 1.9e+20)
       (-
        x
        (/
         (* (fma (fma t z a) z b) y)
         (fma
          (fma (fma (- -15.234687407 z) z -31.4690115749) z -11.9400905721)
          z
          -0.607771387771)))
       t_1))))

C

double code(double x, double y, double z, double t, double a, double b) {
	double t_1 = fma((3.13060547623 + (-1.0 * ((36.52704169880642 + (-1.0 * ((457.9610022158428 + t) / z))) / z))), y, x);
	double tmp;
	if (z <= -1.85e+36) {
		tmp = t_1;
	} else if (z <= 1.9e+20) {
		tmp = x - ((fma(fma(t, z, a), z, b) * y) / fma(fma(fma((-15.234687407 - z), z, -31.4690115749), z, -11.9400905721), z, -0.607771387771));
	} else {
		tmp = t_1;
	}
	return tmp;
}

Julia

function code(x, y, z, t, a, b)
	t_1 = fma(Float64(3.13060547623 + Float64(-1.0 * Float64(Float64(36.52704169880642 + Float64(-1.0 * Float64(Float64(457.9610022158428 + t) / z))) / z))), y, x)
	tmp = 0.0
	if (z <= -1.85e+36)
		tmp = t_1;
	elseif (z <= 1.9e+20)
		tmp = Float64(x - Float64(Float64(fma(fma(t, z, a), z, b) * y) / fma(fma(fma(Float64(-15.234687407 - z), z, -31.4690115749), z, -11.9400905721), z, -0.607771387771)));
	else
		tmp = t_1;
	end
	return tmp
end

Wolfram

code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(N[(3.13060547623 + N[(-1.0 * N[(N[(36.52704169880642 + N[(-1.0 * N[(N[(457.9610022158428 + t), $MachinePrecision] / z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * y + x), $MachinePrecision]}, If[LessEqual[z, -1.85e+36], t$95$1, If[LessEqual[z, 1.9e+20], N[(x - N[(N[(N[(N[(t * z + a), $MachinePrecision] * z + b), $MachinePrecision] * y), $MachinePrecision] / N[(N[(N[(N[(-15.234687407 - z), $MachinePrecision] * z + -31.4690115749), $MachinePrecision] * z + -11.9400905721), $MachinePrecision] * z + -0.607771387771), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$1]]]

Derivation

1. Split input into 2 regimes

2. if z < -1.85000000000000014e36 or 1.9e20 < z

   1. Initial program 58.7%

      \[x + \frac{y \cdot \left(\left(\left(\left(z \cdot 3.13060547623 + 11.1667541262\right) \cdot z + t\right) \cdot z + a\right) \cdot z + b\right)}{\left(\left(\left(z + 15.234687407\right) \cdot z + 31.4690115749\right) \cdot z + 11.9400905721\right) \cdot z + 0.607771387771} \]

   3. Applied rewrites 60.9%

      \[\leadsto \color{blue}{\mathsf{fma}\left(\frac{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(3.13060547623, z, 11.1667541262\right), z, t\right), z, a\right), z, b\right)}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(z - -15.234687407, z, 31.4690115749\right), z, 11.9400905721\right), z, 0.607771387771\right)}, y, x\right)} \]

   4. Taylor expanded in z around -inf

      \[\leadsto \mathsf{fma}\left(\color{blue}{\frac{313060547623}{100000000000} + -1 \cdot \frac{\frac{3652704169880641883561}{100000000000000000000} + -1 \cdot \frac{\frac{45796100221584283915100827016327}{100000000000000000000000000000} + t}{z}}{z}}, y, x\right) \]
   5. Step-by-step derivation

      1. lower-+.f64 N/A

         \[\leadsto \mathsf{fma}\left(\frac{313060547623}{100000000000} + \color{blue}{-1 \cdot \frac{\frac{3652704169880641883561}{100000000000000000000} + -1 \cdot \frac{\frac{45796100221584283915100827016327}{100000000000000000000000000000} + t}{z}}{z}}, y, x\right) \]

      2. lower-*.f64 N/A

         \[\leadsto \mathsf{fma}\left(\frac{313060547623}{100000000000} + -1 \cdot \color{blue}{\frac{\frac{3652704169880641883561}{100000000000000000000} + -1 \cdot \frac{\frac{45796100221584283915100827016327}{100000000000000000000000000000} + t}{z}}{z}}, y, x\right) \]

      3. lower-/.f64 N/A

         \[\leadsto \mathsf{fma}\left(\frac{313060547623}{100000000000} + -1 \cdot \frac{\frac{3652704169880641883561}{100000000000000000000} + -1 \cdot \frac{\frac{45796100221584283915100827016327}{100000000000000000000000000000} + t}{z}}{\color{blue}{z}}, y, x\right) \]

      4. lower-+.f64 N/A

         \[\leadsto \mathsf{fma}\left(\frac{313060547623}{100000000000} + -1 \cdot \frac{\frac{3652704169880641883561}{100000000000000000000} + -1 \cdot \frac{\frac{45796100221584283915100827016327}{100000000000000000000000000000} + t}{z}}{z}, y, x\right) \]

      5. lower-*.f64 N/A

         \[\leadsto \mathsf{fma}\left(\frac{313060547623}{100000000000} + -1 \cdot \frac{\frac{3652704169880641883561}{100000000000000000000} + -1 \cdot \frac{\frac{45796100221584283915100827016327}{100000000000000000000000000000} + t}{z}}{z}, y, x\right) \]

      6. lower-/.f64 N/A

         \[\leadsto \mathsf{fma}\left(\frac{313060547623}{100000000000} + -1 \cdot \frac{\frac{3652704169880641883561}{100000000000000000000} + -1 \cdot \frac{\frac{45796100221584283915100827016327}{100000000000000000000000000000} + t}{z}}{z}, y, x\right) \]

      7. lower-+.f64 55.8

         \[\leadsto \mathsf{fma}\left(3.13060547623 + -1 \cdot \frac{36.52704169880642 + -1 \cdot \frac{457.9610022158428 + t}{z}}{z}, y, x\right) \]

   6. Applied rewrites 55.8%

      \[\leadsto \mathsf{fma}\left(\color{blue}{3.13060547623 + -1 \cdot \frac{36.52704169880642 + -1 \cdot \frac{457.9610022158428 + t}{z}}{z}}, y, x\right) \]

   if -1.85000000000000014e36 < z < 1.9e20

   1. Initial program 58.7%

      \[x + \frac{y \cdot \left(\left(\left(\left(z \cdot 3.13060547623 + 11.1667541262\right) \cdot z + t\right) \cdot z + a\right) \cdot z + b\right)}{\left(\left(\left(z + 15.234687407\right) \cdot z + 31.4690115749\right) \cdot z + 11.9400905721\right) \cdot z + 0.607771387771} \]

   2. Taylor expanded in z around 0

      \[\leadsto x + \frac{\color{blue}{b \cdot y + z \cdot \left(a \cdot y + t \cdot \left(y \cdot z\right)\right)}}{\left(\left(\left(z + \frac{15234687407}{1000000000}\right) \cdot z + \frac{314690115749}{10000000000}\right) \cdot z + \frac{119400905721}{10000000000}\right) \cdot z + \frac{607771387771}{1000000000000}} \]
   3. Step-by-step derivation

      1. lower-fma.f64 N/A

         \[\leadsto x + \frac{\mathsf{fma}\left(b, \color{blue}{y}, z \cdot \left(a \cdot y + t \cdot \left(y \cdot z\right)\right)\right)}{\left(\left(\left(z + \frac{15234687407}{1000000000}\right) \cdot z + \frac{314690115749}{10000000000}\right) \cdot z + \frac{119400905721}{10000000000}\right) \cdot z + \frac{607771387771}{1000000000000}} \]

      2. lower-*.f64 N/A

         \[\leadsto x + \frac{\mathsf{fma}\left(b, y, z \cdot \left(a \cdot y + t \cdot \left(y \cdot z\right)\right)\right)}{\left(\left(\left(z + \frac{15234687407}{1000000000}\right) \cdot z + \frac{314690115749}{10000000000}\right) \cdot z + \frac{119400905721}{10000000000}\right) \cdot z + \frac{607771387771}{1000000000000}} \]

      3. lower-fma.f64 N/A

         \[\leadsto x + \frac{\mathsf{fma}\left(b, y, z \cdot \mathsf{fma}\left(a, y, t \cdot \left(y \cdot z\right)\right)\right)}{\left(\left(\left(z + \frac{15234687407}{1000000000}\right) \cdot z + \frac{314690115749}{10000000000}\right) \cdot z + \frac{119400905721}{10000000000}\right) \cdot z + \frac{607771387771}{1000000000000}} \]

      4. lower-*.f64 N/A

         \[\leadsto x + \frac{\mathsf{fma}\left(b, y, z \cdot \mathsf{fma}\left(a, y, t \cdot \left(y \cdot z\right)\right)\right)}{\left(\left(\left(z + \frac{15234687407}{1000000000}\right) \cdot z + \frac{314690115749}{10000000000}\right) \cdot z + \frac{119400905721}{10000000000}\right) \cdot z + \frac{607771387771}{1000000000000}} \]

      5. lower-*.f64 62.7

         \[\leadsto x + \frac{\mathsf{fma}\left(b, y, z \cdot \mathsf{fma}\left(a, y, t \cdot \left(y \cdot z\right)\right)\right)}{\left(\left(\left(z + 15.234687407\right) \cdot z + 31.4690115749\right) \cdot z + 11.9400905721\right) \cdot z + 0.607771387771} \]

   4. Applied rewrites 62.7%

      \[\leadsto x + \frac{\color{blue}{\mathsf{fma}\left(b, y, z \cdot \mathsf{fma}\left(a, y, t \cdot \left(y \cdot z\right)\right)\right)}}{\left(\left(\left(z + 15.234687407\right) \cdot z + 31.4690115749\right) \cdot z + 11.9400905721\right) \cdot z + 0.607771387771} \]

   6. Applied rewrites 61.6%

      \[\leadsto \color{blue}{x - \frac{\mathsf{fma}\left(y \cdot \mathsf{fma}\left(t, z, a\right), z, b \cdot y\right)}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(-15.234687407 - z, z, -31.4690115749\right), z, -11.9400905721\right), z, -0.607771387771\right)}} \]
   7. Step-by-step derivation

      1. lift-fma.f64 N/A

         \[\leadsto x - \frac{\left(y \cdot \mathsf{fma}\left(t, z, a\right)\right) \cdot z + \color{blue}{b \cdot y}}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(\frac{-15234687407}{1000000000} - z, z, \frac{-314690115749}{10000000000}\right), z, \frac{-119400905721}{10000000000}\right), z, \frac{-607771387771}{1000000000000}\right)} \]

      2. lift-*.f64 N/A

         \[\leadsto x - \frac{\left(y \cdot \mathsf{fma}\left(t, z, a\right)\right) \cdot z + b \cdot y}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(\frac{-15234687407}{1000000000} - z, z, \frac{-314690115749}{10000000000}\right), z, \frac{-119400905721}{10000000000}\right), z, \frac{-607771387771}{1000000000000}\right)} \]

      3. associate-*l* N/A

         \[\leadsto x - \frac{y \cdot \left(\mathsf{fma}\left(t, z, a\right) \cdot z\right) + \color{blue}{b} \cdot y}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(\frac{-15234687407}{1000000000} - z, z, \frac{-314690115749}{10000000000}\right), z, \frac{-119400905721}{10000000000}\right), z, \frac{-607771387771}{1000000000000}\right)} \]

      4. lift-*.f64 N/A

         \[\leadsto x - \frac{y \cdot \left(\mathsf{fma}\left(t, z, a\right) \cdot z\right) + b \cdot \color{blue}{y}}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(\frac{-15234687407}{1000000000} - z, z, \frac{-314690115749}{10000000000}\right), z, \frac{-119400905721}{10000000000}\right), z, \frac{-607771387771}{1000000000000}\right)} \]

      5. *-commutative N/A

         \[\leadsto x - \frac{y \cdot \left(\mathsf{fma}\left(t, z, a\right) \cdot z\right) + y \cdot \color{blue}{b}}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(\frac{-15234687407}{1000000000} - z, z, \frac{-314690115749}{10000000000}\right), z, \frac{-119400905721}{10000000000}\right), z, \frac{-607771387771}{1000000000000}\right)} \]

      6. distribute-lft-out N/A

         \[\leadsto x - \frac{y \cdot \color{blue}{\left(\mathsf{fma}\left(t, z, a\right) \cdot z + b\right)}}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(\frac{-15234687407}{1000000000} - z, z, \frac{-314690115749}{10000000000}\right), z, \frac{-119400905721}{10000000000}\right), z, \frac{-607771387771}{1000000000000}\right)} \]

      7. *-commutative N/A

         \[\leadsto x - \frac{\left(\mathsf{fma}\left(t, z, a\right) \cdot z + b\right) \cdot \color{blue}{y}}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(\frac{-15234687407}{1000000000} - z, z, \frac{-314690115749}{10000000000}\right), z, \frac{-119400905721}{10000000000}\right), z, \frac{-607771387771}{1000000000000}\right)} \]

      8. lower-*.f64 N/A

         \[\leadsto x - \frac{\left(\mathsf{fma}\left(t, z, a\right) \cdot z + b\right) \cdot \color{blue}{y}}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(\frac{-15234687407}{1000000000} - z, z, \frac{-314690115749}{10000000000}\right), z, \frac{-119400905721}{10000000000}\right), z, \frac{-607771387771}{1000000000000}\right)} \]

      9. lower-fma.f64 61.8

         \[\leadsto x - \frac{\mathsf{fma}\left(\mathsf{fma}\left(t, z, a\right), z, b\right) \cdot y}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(-15.234687407 - z, z, -31.4690115749\right), z, -11.9400905721\right), z, -0.607771387771\right)} \]

   8. Applied rewrites 61.8%

      \[\leadsto x - \frac{\mathsf{fma}\left(\mathsf{fma}\left(t, z, a\right), z, b\right) \cdot \color{blue}{y}}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(-15.234687407 - z, z, -31.4690115749\right), z, -11.9400905721\right), z, -0.607771387771\right)} \]
3. Recombined 2 regimes into one program.
4. Add Preprocessing

Alternative 5: 96.0% accurate, 1.2× speedup

Math

\[\begin{array}{l} \\ \begin{array}{l} t_1 := \mathsf{fma}\left(3.13060547623 + -1 \cdot \frac{36.52704169880642 + -1 \cdot \frac{457.9610022158428 + t}{z}}{z}, y, x\right)\\ \mathbf{if}\;z \leq -7 \cdot 10^{+32}:\\ \;\;\;\;t\_1\\ \mathbf{elif}\;z \leq -7 \cdot 10^{-15}:\\ \;\;\;\;\mathsf{fma}\left(\frac{\mathsf{fma}\left(a, z, b\right)}{\mathsf{fma}\left(\mathsf{fma}\left(z - -15.234687407, z, 31.4690115749\right), z \cdot z, \mathsf{fma}\left(11.9400905721, z, 0.607771387771\right)\right)}, y, x\right)\\ \mathbf{elif}\;z \leq 5.7 \cdot 10^{+19}:\\ \;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(3.13060547623, z, 11.1667541262\right), z, t\right), z, a\right), z, b\right), \left(-y\right) \cdot -1.6453555072203998, x\right)\\ \mathbf{else}:\\ \;\;\;\;t\_1\\ \end{array} \end{array} \]

FPCore

(FPCore (x y z t a b)
 :precision binary64
 (let* ((t_1
         (fma
          (+
           3.13060547623
           (*
            -1.0
            (/
             (+ 36.52704169880642 (* -1.0 (/ (+ 457.9610022158428 t) z)))
             z)))
          y
          x)))
   (if (<= z -7e+32)
     t_1
     (if (<= z -7e-15)
       (fma
        (/
         (fma a z b)
         (fma
          (fma (- z -15.234687407) z 31.4690115749)
          (* z z)
          (fma 11.9400905721 z 0.607771387771)))
        y
        x)
       (if (<= z 5.7e+19)
         (fma
          (fma (fma (fma (fma 3.13060547623 z 11.1667541262) z t) z a) z b)
          (* (- y) -1.6453555072203998)
          x)
         t_1)))))

C

double code(double x, double y, double z, double t, double a, double b) {
	double t_1 = fma((3.13060547623 + (-1.0 * ((36.52704169880642 + (-1.0 * ((457.9610022158428 + t) / z))) / z))), y, x);
	double tmp;
	if (z <= -7e+32) {
		tmp = t_1;
	} else if (z <= -7e-15) {
		tmp = fma((fma(a, z, b) / fma(fma((z - -15.234687407), z, 31.4690115749), (z * z), fma(11.9400905721, z, 0.607771387771))), y, x);
	} else if (z <= 5.7e+19) {
		tmp = fma(fma(fma(fma(fma(3.13060547623, z, 11.1667541262), z, t), z, a), z, b), (-y * -1.6453555072203998), x);
	} else {
		tmp = t_1;
	}
	return tmp;
}

Julia

function code(x, y, z, t, a, b)
	t_1 = fma(Float64(3.13060547623 + Float64(-1.0 * Float64(Float64(36.52704169880642 + Float64(-1.0 * Float64(Float64(457.9610022158428 + t) / z))) / z))), y, x)
	tmp = 0.0
	if (z <= -7e+32)
		tmp = t_1;
	elseif (z <= -7e-15)
		tmp = fma(Float64(fma(a, z, b) / fma(fma(Float64(z - -15.234687407), z, 31.4690115749), Float64(z * z), fma(11.9400905721, z, 0.607771387771))), y, x);
	elseif (z <= 5.7e+19)
		tmp = fma(fma(fma(fma(fma(3.13060547623, z, 11.1667541262), z, t), z, a), z, b), Float64(Float64(-y) * -1.6453555072203998), x);
	else
		tmp = t_1;
	end
	return tmp
end

Wolfram

code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(N[(3.13060547623 + N[(-1.0 * N[(N[(36.52704169880642 + N[(-1.0 * N[(N[(457.9610022158428 + t), $MachinePrecision] / z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * y + x), $MachinePrecision]}, If[LessEqual[z, -7e+32], t$95$1, If[LessEqual[z, -7e-15], N[(N[(N[(a * z + b), $MachinePrecision] / N[(N[(N[(z - -15.234687407), $MachinePrecision] * z + 31.4690115749), $MachinePrecision] * N[(z * z), $MachinePrecision] + N[(11.9400905721 * z + 0.607771387771), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * y + x), $MachinePrecision], If[LessEqual[z, 5.7e+19], N[(N[(N[(N[(N[(3.13060547623 * z + 11.1667541262), $MachinePrecision] * z + t), $MachinePrecision] * z + a), $MachinePrecision] * z + b), $MachinePrecision] * N[((-y) * -1.6453555072203998), $MachinePrecision] + x), $MachinePrecision], t$95$1]]]]

Derivation

1. Split input into 3 regimes

2. if z < -7.0000000000000002e32 or 5.7e19 < z

   1. Initial program 58.7%

      \[x + \frac{y \cdot \left(\left(\left(\left(z \cdot 3.13060547623 + 11.1667541262\right) \cdot z + t\right) \cdot z + a\right) \cdot z + b\right)}{\left(\left(\left(z + 15.234687407\right) \cdot z + 31.4690115749\right) \cdot z + 11.9400905721\right) \cdot z + 0.607771387771} \]

   3. Applied rewrites 60.9%

      \[\leadsto \color{blue}{\mathsf{fma}\left(\frac{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(3.13060547623, z, 11.1667541262\right), z, t\right), z, a\right), z, b\right)}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(z - -15.234687407, z, 31.4690115749\right), z, 11.9400905721\right), z, 0.607771387771\right)}, y, x\right)} \]

   4. Taylor expanded in z around -inf

      \[\leadsto \mathsf{fma}\left(\color{blue}{\frac{313060547623}{100000000000} + -1 \cdot \frac{\frac{3652704169880641883561}{100000000000000000000} + -1 \cdot \frac{\frac{45796100221584283915100827016327}{100000000000000000000000000000} + t}{z}}{z}}, y, x\right) \]
   5. Step-by-step derivation

      1. lower-+.f64 N/A

         \[\leadsto \mathsf{fma}\left(\frac{313060547623}{100000000000} + \color{blue}{-1 \cdot \frac{\frac{3652704169880641883561}{100000000000000000000} + -1 \cdot \frac{\frac{45796100221584283915100827016327}{100000000000000000000000000000} + t}{z}}{z}}, y, x\right) \]

      2. lower-*.f64 N/A

         \[\leadsto \mathsf{fma}\left(\frac{313060547623}{100000000000} + -1 \cdot \color{blue}{\frac{\frac{3652704169880641883561}{100000000000000000000} + -1 \cdot \frac{\frac{45796100221584283915100827016327}{100000000000000000000000000000} + t}{z}}{z}}, y, x\right) \]

      3. lower-/.f64 N/A

         \[\leadsto \mathsf{fma}\left(\frac{313060547623}{100000000000} + -1 \cdot \frac{\frac{3652704169880641883561}{100000000000000000000} + -1 \cdot \frac{\frac{45796100221584283915100827016327}{100000000000000000000000000000} + t}{z}}{\color{blue}{z}}, y, x\right) \]

      4. lower-+.f64 N/A

         \[\leadsto \mathsf{fma}\left(\frac{313060547623}{100000000000} + -1 \cdot \frac{\frac{3652704169880641883561}{100000000000000000000} + -1 \cdot \frac{\frac{45796100221584283915100827016327}{100000000000000000000000000000} + t}{z}}{z}, y, x\right) \]

      5. lower-*.f64 N/A

         \[\leadsto \mathsf{fma}\left(\frac{313060547623}{100000000000} + -1 \cdot \frac{\frac{3652704169880641883561}{100000000000000000000} + -1 \cdot \frac{\frac{45796100221584283915100827016327}{100000000000000000000000000000} + t}{z}}{z}, y, x\right) \]

      6. lower-/.f64 N/A

         \[\leadsto \mathsf{fma}\left(\frac{313060547623}{100000000000} + -1 \cdot \frac{\frac{3652704169880641883561}{100000000000000000000} + -1 \cdot \frac{\frac{45796100221584283915100827016327}{100000000000000000000000000000} + t}{z}}{z}, y, x\right) \]

      7. lower-+.f64 55.8

         \[\leadsto \mathsf{fma}\left(3.13060547623 + -1 \cdot \frac{36.52704169880642 + -1 \cdot \frac{457.9610022158428 + t}{z}}{z}, y, x\right) \]

   6. Applied rewrites 55.8%

      \[\leadsto \mathsf{fma}\left(\color{blue}{3.13060547623 + -1 \cdot \frac{36.52704169880642 + -1 \cdot \frac{457.9610022158428 + t}{z}}{z}}, y, x\right) \]

   if -7.0000000000000002e32 < z < -7.0000000000000001e-15

   1. Initial program 58.7%

      \[x + \frac{y \cdot \left(\left(\left(\left(z \cdot 3.13060547623 + 11.1667541262\right) \cdot z + t\right) \cdot z + a\right) \cdot z + b\right)}{\left(\left(\left(z + 15.234687407\right) \cdot z + 31.4690115749\right) \cdot z + 11.9400905721\right) \cdot z + 0.607771387771} \]

   3. Applied rewrites 60.9%

      \[\leadsto \color{blue}{\mathsf{fma}\left(\frac{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(3.13060547623, z, 11.1667541262\right), z, t\right), z, a\right), z, b\right)}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(z - -15.234687407, z, 31.4690115749\right), z, 11.9400905721\right), z, 0.607771387771\right)}, y, x\right)} \]

   4. Taylor expanded in z around 0

      \[\leadsto \mathsf{fma}\left(\frac{\mathsf{fma}\left(\color{blue}{a}, z, b\right)}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(z - \frac{-15234687407}{1000000000}, z, \frac{314690115749}{10000000000}\right), z, \frac{119400905721}{10000000000}\right), z, \frac{607771387771}{1000000000000}\right)}, y, x\right) \]
   5. Step-by-step derivation

   6. Applied rewrites 67.3%

      \[\leadsto \mathsf{fma}\left(\frac{\mathsf{fma}\left(\color{blue}{a}, z, b\right)}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(z - -15.234687407, z, 31.4690115749\right), z, 11.9400905721\right), z, 0.607771387771\right)}, y, x\right) \]
   7. Step-by-step derivation

      1. lift-fma.f64 N/A

         \[\leadsto \mathsf{fma}\left(\frac{\mathsf{fma}\left(a, z, b\right)}{\color{blue}{\mathsf{fma}\left(\mathsf{fma}\left(z - \frac{-15234687407}{1000000000}, z, \frac{314690115749}{10000000000}\right), z, \frac{119400905721}{10000000000}\right) \cdot z + \frac{607771387771}{1000000000000}}}, y, x\right) \]

      2. *-commutative N/A

         \[\leadsto \mathsf{fma}\left(\frac{\mathsf{fma}\left(a, z, b\right)}{\color{blue}{z \cdot \mathsf{fma}\left(\mathsf{fma}\left(z - \frac{-15234687407}{1000000000}, z, \frac{314690115749}{10000000000}\right), z, \frac{119400905721}{10000000000}\right)} + \frac{607771387771}{1000000000000}}, y, x\right) \]

      3. lift-fma.f64 N/A

         \[\leadsto \mathsf{fma}\left(\frac{\mathsf{fma}\left(a, z, b\right)}{z \cdot \color{blue}{\left(\mathsf{fma}\left(z - \frac{-15234687407}{1000000000}, z, \frac{314690115749}{10000000000}\right) \cdot z + \frac{119400905721}{10000000000}\right)} + \frac{607771387771}{1000000000000}}, y, x\right) \]

      4. distribute-rgt-in N/A

         \[\leadsto \mathsf{fma}\left(\frac{\mathsf{fma}\left(a, z, b\right)}{\color{blue}{\left(\left(\mathsf{fma}\left(z - \frac{-15234687407}{1000000000}, z, \frac{314690115749}{10000000000}\right) \cdot z\right) \cdot z + \frac{119400905721}{10000000000} \cdot z\right)} + \frac{607771387771}{1000000000000}}, y, x\right) \]

      5. associate-+l+ N/A

         \[\leadsto \mathsf{fma}\left(\frac{\mathsf{fma}\left(a, z, b\right)}{\color{blue}{\left(\mathsf{fma}\left(z - \frac{-15234687407}{1000000000}, z, \frac{314690115749}{10000000000}\right) \cdot z\right) \cdot z + \left(\frac{119400905721}{10000000000} \cdot z + \frac{607771387771}{1000000000000}\right)}}, y, x\right) \]

      6. *-commutative N/A

         \[\leadsto \mathsf{fma}\left(\frac{\mathsf{fma}\left(a, z, b\right)}{\left(\mathsf{fma}\left(z - \frac{-15234687407}{1000000000}, z, \frac{314690115749}{10000000000}\right) \cdot z\right) \cdot z + \left(\color{blue}{z \cdot \frac{119400905721}{10000000000}} + \frac{607771387771}{1000000000000}\right)}, y, x\right) \]

      7. associate-*l* N/A

         \[\leadsto \mathsf{fma}\left(\frac{\mathsf{fma}\left(a, z, b\right)}{\color{blue}{\mathsf{fma}\left(z - \frac{-15234687407}{1000000000}, z, \frac{314690115749}{10000000000}\right) \cdot \left(z \cdot z\right)} + \left(z \cdot \frac{119400905721}{10000000000} + \frac{607771387771}{1000000000000}\right)}, y, x\right) \]

      8. lower-fma.f64 N/A

         \[\leadsto \mathsf{fma}\left(\frac{\mathsf{fma}\left(a, z, b\right)}{\color{blue}{\mathsf{fma}\left(\mathsf{fma}\left(z - \frac{-15234687407}{1000000000}, z, \frac{314690115749}{10000000000}\right), z \cdot z, z \cdot \frac{119400905721}{10000000000} + \frac{607771387771}{1000000000000}\right)}}, y, x\right) \]

      9. lower-*.f64 N/A

         \[\leadsto \mathsf{fma}\left(\frac{\mathsf{fma}\left(a, z, b\right)}{\mathsf{fma}\left(\mathsf{fma}\left(z - \frac{-15234687407}{1000000000}, z, \frac{314690115749}{10000000000}\right), \color{blue}{z \cdot z}, z \cdot \frac{119400905721}{10000000000} + \frac{607771387771}{1000000000000}\right)}, y, x\right) \]

      10. *-commutative N/A

          \[\leadsto \mathsf{fma}\left(\frac{\mathsf{fma}\left(a, z, b\right)}{\mathsf{fma}\left(\mathsf{fma}\left(z - \frac{-15234687407}{1000000000}, z, \frac{314690115749}{10000000000}\right), z \cdot z, \color{blue}{\frac{119400905721}{10000000000} \cdot z} + \frac{607771387771}{1000000000000}\right)}, y, x\right) \]

      11. lower-fma.f64 67.2

          \[\leadsto \mathsf{fma}\left(\frac{\mathsf{fma}\left(a, z, b\right)}{\mathsf{fma}\left(\mathsf{fma}\left(z - -15.234687407, z, 31.4690115749\right), z \cdot z, \color{blue}{\mathsf{fma}\left(11.9400905721, z, 0.607771387771\right)}\right)}, y, x\right) \]

   8. Applied rewrites 67.2%

      \[\leadsto \mathsf{fma}\left(\frac{\mathsf{fma}\left(a, z, b\right)}{\color{blue}{\mathsf{fma}\left(\mathsf{fma}\left(z - -15.234687407, z, 31.4690115749\right), z \cdot z, \mathsf{fma}\left(11.9400905721, z, 0.607771387771\right)\right)}}, y, x\right) \]

   if -7.0000000000000001e-15 < z < 5.7e19

   1. Initial program 58.7%

      \[x + \frac{y \cdot \left(\left(\left(\left(z \cdot 3.13060547623 + 11.1667541262\right) \cdot z + t\right) \cdot z + a\right) \cdot z + b\right)}{\left(\left(\left(z + 15.234687407\right) \cdot z + 31.4690115749\right) \cdot z + 11.9400905721\right) \cdot z + 0.607771387771} \]

   3. Applied rewrites 60.6%

      \[\leadsto \color{blue}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(3.13060547623, z, 11.1667541262\right), z, t\right), z, a\right), z, b\right), \left(-y\right) \cdot \frac{-1}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(z - -15.234687407, z, 31.4690115749\right), z, 11.9400905721\right), z, 0.607771387771\right)}, x\right)} \]

   4. Taylor expanded in z around 0

      \[\leadsto \mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(\frac{313060547623}{100000000000}, z, \frac{55833770631}{5000000000}\right), z, t\right), z, a\right), z, b\right), \left(-y\right) \cdot \color{blue}{\frac{-1000000000000}{607771387771}}, x\right) \]
   5. Step-by-step derivation

   6. Applied rewrites 55.0%

      \[\leadsto \mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(3.13060547623, z, 11.1667541262\right), z, t\right), z, a\right), z, b\right), \left(-y\right) \cdot \color{blue}{-1.6453555072203998}, x\right) \]
3. Recombined 3 regimes into one program.
4. Add Preprocessing

Alternative 6: 95.9% accurate, 1.3× speedup

Math

\[\begin{array}{l} \\ \begin{array}{l} t_1 := \mathsf{fma}\left(3.13060547623 + -1 \cdot \frac{36.52704169880642 + -1 \cdot \frac{457.9610022158428 + t}{z}}{z}, y, x\right)\\ \mathbf{if}\;z \leq -7 \cdot 10^{+32}:\\ \;\;\;\;t\_1\\ \mathbf{elif}\;z \leq -7 \cdot 10^{-15}:\\ \;\;\;\;\mathsf{fma}\left(\frac{\mathsf{fma}\left(a, z, b\right)}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(z - -15.234687407, z, 31.4690115749\right), z, 11.9400905721\right), z, 0.607771387771\right)}, y, x\right)\\ \mathbf{elif}\;z \leq 5.7 \cdot 10^{+19}:\\ \;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(3.13060547623, z, 11.1667541262\right), z, t\right), z, a\right), z, b\right), \left(-y\right) \cdot -1.6453555072203998, x\right)\\ \mathbf{else}:\\ \;\;\;\;t\_1\\ \end{array} \end{array} \]

FPCore

(FPCore (x y z t a b)
 :precision binary64
 (let* ((t_1
         (fma
          (+
           3.13060547623
           (*
            -1.0
            (/
             (+ 36.52704169880642 (* -1.0 (/ (+ 457.9610022158428 t) z)))
             z)))
          y
          x)))
   (if (<= z -7e+32)
     t_1
     (if (<= z -7e-15)
       (fma
        (/
         (fma a z b)
         (fma
          (fma (fma (- z -15.234687407) z 31.4690115749) z 11.9400905721)
          z
          0.607771387771))
        y
        x)
       (if (<= z 5.7e+19)
         (fma
          (fma (fma (fma (fma 3.13060547623 z 11.1667541262) z t) z a) z b)
          (* (- y) -1.6453555072203998)
          x)
         t_1)))))

C

double code(double x, double y, double z, double t, double a, double b) {
	double t_1 = fma((3.13060547623 + (-1.0 * ((36.52704169880642 + (-1.0 * ((457.9610022158428 + t) / z))) / z))), y, x);
	double tmp;
	if (z <= -7e+32) {
		tmp = t_1;
	} else if (z <= -7e-15) {
		tmp = fma((fma(a, z, b) / fma(fma(fma((z - -15.234687407), z, 31.4690115749), z, 11.9400905721), z, 0.607771387771)), y, x);
	} else if (z <= 5.7e+19) {
		tmp = fma(fma(fma(fma(fma(3.13060547623, z, 11.1667541262), z, t), z, a), z, b), (-y * -1.6453555072203998), x);
	} else {
		tmp = t_1;
	}
	return tmp;
}

Julia

function code(x, y, z, t, a, b)
	t_1 = fma(Float64(3.13060547623 + Float64(-1.0 * Float64(Float64(36.52704169880642 + Float64(-1.0 * Float64(Float64(457.9610022158428 + t) / z))) / z))), y, x)
	tmp = 0.0
	if (z <= -7e+32)
		tmp = t_1;
	elseif (z <= -7e-15)
		tmp = fma(Float64(fma(a, z, b) / fma(fma(fma(Float64(z - -15.234687407), z, 31.4690115749), z, 11.9400905721), z, 0.607771387771)), y, x);
	elseif (z <= 5.7e+19)
		tmp = fma(fma(fma(fma(fma(3.13060547623, z, 11.1667541262), z, t), z, a), z, b), Float64(Float64(-y) * -1.6453555072203998), x);
	else
		tmp = t_1;
	end
	return tmp
end

Wolfram

code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(N[(3.13060547623 + N[(-1.0 * N[(N[(36.52704169880642 + N[(-1.0 * N[(N[(457.9610022158428 + t), $MachinePrecision] / z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * y + x), $MachinePrecision]}, If[LessEqual[z, -7e+32], t$95$1, If[LessEqual[z, -7e-15], N[(N[(N[(a * z + b), $MachinePrecision] / N[(N[(N[(N[(z - -15.234687407), $MachinePrecision] * z + 31.4690115749), $MachinePrecision] * z + 11.9400905721), $MachinePrecision] * z + 0.607771387771), $MachinePrecision]), $MachinePrecision] * y + x), $MachinePrecision], If[LessEqual[z, 5.7e+19], N[(N[(N[(N[(N[(3.13060547623 * z + 11.1667541262), $MachinePrecision] * z + t), $MachinePrecision] * z + a), $MachinePrecision] * z + b), $MachinePrecision] * N[((-y) * -1.6453555072203998), $MachinePrecision] + x), $MachinePrecision], t$95$1]]]]

Derivation

1. Split input into 3 regimes

2. if z < -7.0000000000000002e32 or 5.7e19 < z

   1. Initial program 58.7%

      \[x + \frac{y \cdot \left(\left(\left(\left(z \cdot 3.13060547623 + 11.1667541262\right) \cdot z + t\right) \cdot z + a\right) \cdot z + b\right)}{\left(\left(\left(z + 15.234687407\right) \cdot z + 31.4690115749\right) \cdot z + 11.9400905721\right) \cdot z + 0.607771387771} \]

   3. Applied rewrites 60.9%

      \[\leadsto \color{blue}{\mathsf{fma}\left(\frac{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(3.13060547623, z, 11.1667541262\right), z, t\right), z, a\right), z, b\right)}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(z - -15.234687407, z, 31.4690115749\right), z, 11.9400905721\right), z, 0.607771387771\right)}, y, x\right)} \]

   4. Taylor expanded in z around -inf

      \[\leadsto \mathsf{fma}\left(\color{blue}{\frac{313060547623}{100000000000} + -1 \cdot \frac{\frac{3652704169880641883561}{100000000000000000000} + -1 \cdot \frac{\frac{45796100221584283915100827016327}{100000000000000000000000000000} + t}{z}}{z}}, y, x\right) \]
   5. Step-by-step derivation

      1. lower-+.f64 N/A

         \[\leadsto \mathsf{fma}\left(\frac{313060547623}{100000000000} + \color{blue}{-1 \cdot \frac{\frac{3652704169880641883561}{100000000000000000000} + -1 \cdot \frac{\frac{45796100221584283915100827016327}{100000000000000000000000000000} + t}{z}}{z}}, y, x\right) \]

      2. lower-*.f64 N/A

         \[\leadsto \mathsf{fma}\left(\frac{313060547623}{100000000000} + -1 \cdot \color{blue}{\frac{\frac{3652704169880641883561}{100000000000000000000} + -1 \cdot \frac{\frac{45796100221584283915100827016327}{100000000000000000000000000000} + t}{z}}{z}}, y, x\right) \]

      3. lower-/.f64 N/A

         \[\leadsto \mathsf{fma}\left(\frac{313060547623}{100000000000} + -1 \cdot \frac{\frac{3652704169880641883561}{100000000000000000000} + -1 \cdot \frac{\frac{45796100221584283915100827016327}{100000000000000000000000000000} + t}{z}}{\color{blue}{z}}, y, x\right) \]

      4. lower-+.f64 N/A

         \[\leadsto \mathsf{fma}\left(\frac{313060547623}{100000000000} + -1 \cdot \frac{\frac{3652704169880641883561}{100000000000000000000} + -1 \cdot \frac{\frac{45796100221584283915100827016327}{100000000000000000000000000000} + t}{z}}{z}, y, x\right) \]

      5. lower-*.f64 N/A

         \[\leadsto \mathsf{fma}\left(\frac{313060547623}{100000000000} + -1 \cdot \frac{\frac{3652704169880641883561}{100000000000000000000} + -1 \cdot \frac{\frac{45796100221584283915100827016327}{100000000000000000000000000000} + t}{z}}{z}, y, x\right) \]

      6. lower-/.f64 N/A

         \[\leadsto \mathsf{fma}\left(\frac{313060547623}{100000000000} + -1 \cdot \frac{\frac{3652704169880641883561}{100000000000000000000} + -1 \cdot \frac{\frac{45796100221584283915100827016327}{100000000000000000000000000000} + t}{z}}{z}, y, x\right) \]

      7. lower-+.f64 55.8

         \[\leadsto \mathsf{fma}\left(3.13060547623 + -1 \cdot \frac{36.52704169880642 + -1 \cdot \frac{457.9610022158428 + t}{z}}{z}, y, x\right) \]

   6. Applied rewrites 55.8%

      \[\leadsto \mathsf{fma}\left(\color{blue}{3.13060547623 + -1 \cdot \frac{36.52704169880642 + -1 \cdot \frac{457.9610022158428 + t}{z}}{z}}, y, x\right) \]

   if -7.0000000000000002e32 < z < -7.0000000000000001e-15

   1. Initial program 58.7%

      \[x + \frac{y \cdot \left(\left(\left(\left(z \cdot 3.13060547623 + 11.1667541262\right) \cdot z + t\right) \cdot z + a\right) \cdot z + b\right)}{\left(\left(\left(z + 15.234687407\right) \cdot z + 31.4690115749\right) \cdot z + 11.9400905721\right) \cdot z + 0.607771387771} \]

   3. Applied rewrites 60.9%

      \[\leadsto \color{blue}{\mathsf{fma}\left(\frac{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(3.13060547623, z, 11.1667541262\right), z, t\right), z, a\right), z, b\right)}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(z - -15.234687407, z, 31.4690115749\right), z, 11.9400905721\right), z, 0.607771387771\right)}, y, x\right)} \]

   4. Taylor expanded in z around 0

      \[\leadsto \mathsf{fma}\left(\frac{\mathsf{fma}\left(\color{blue}{a}, z, b\right)}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(z - \frac{-15234687407}{1000000000}, z, \frac{314690115749}{10000000000}\right), z, \frac{119400905721}{10000000000}\right), z, \frac{607771387771}{1000000000000}\right)}, y, x\right) \]
   5. Step-by-step derivation

   6. Applied rewrites 67.3%

      \[\leadsto \mathsf{fma}\left(\frac{\mathsf{fma}\left(\color{blue}{a}, z, b\right)}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(z - -15.234687407, z, 31.4690115749\right), z, 11.9400905721\right), z, 0.607771387771\right)}, y, x\right) \]

   if -7.0000000000000001e-15 < z < 5.7e19

   1. Initial program 58.7%

      \[x + \frac{y \cdot \left(\left(\left(\left(z \cdot 3.13060547623 + 11.1667541262\right) \cdot z + t\right) \cdot z + a\right) \cdot z + b\right)}{\left(\left(\left(z + 15.234687407\right) \cdot z + 31.4690115749\right) \cdot z + 11.9400905721\right) \cdot z + 0.607771387771} \]

   3. Applied rewrites 60.6%

      \[\leadsto \color{blue}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(3.13060547623, z, 11.1667541262\right), z, t\right), z, a\right), z, b\right), \left(-y\right) \cdot \frac{-1}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(z - -15.234687407, z, 31.4690115749\right), z, 11.9400905721\right), z, 0.607771387771\right)}, x\right)} \]

   4. Taylor expanded in z around 0

      \[\leadsto \mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(\frac{313060547623}{100000000000}, z, \frac{55833770631}{5000000000}\right), z, t\right), z, a\right), z, b\right), \left(-y\right) \cdot \color{blue}{\frac{-1000000000000}{607771387771}}, x\right) \]
   5. Step-by-step derivation

   6. Applied rewrites 55.0%

      \[\leadsto \mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(3.13060547623, z, 11.1667541262\right), z, t\right), z, a\right), z, b\right), \left(-y\right) \cdot \color{blue}{-1.6453555072203998}, x\right) \]
3. Recombined 3 regimes into one program.
4. Add Preprocessing

Alternative 7: 95.9% accurate, 1.4× speedup

Math

\[\begin{array}{l} \\ \begin{array}{l} t_1 := \mathsf{fma}\left(3.13060547623 + -1 \cdot \frac{36.52704169880642 + -1 \cdot \frac{457.9610022158428 + t}{z}}{z}, y, x\right)\\ \mathbf{if}\;z \leq -0.06:\\ \;\;\;\;t\_1\\ \mathbf{elif}\;z \leq 5.7 \cdot 10^{+19}:\\ \;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(3.13060547623, z, 11.1667541262\right), z, t\right), z, a\right), z, b\right), \left(-y\right) \cdot -1.6453555072203998, x\right)\\ \mathbf{else}:\\ \;\;\;\;t\_1\\ \end{array} \end{array} \]

FPCore

(FPCore (x y z t a b)
 :precision binary64
 (let* ((t_1
         (fma
          (+
           3.13060547623
           (*
            -1.0
            (/
             (+ 36.52704169880642 (* -1.0 (/ (+ 457.9610022158428 t) z)))
             z)))
          y
          x)))
   (if (<= z -0.06)
     t_1
     (if (<= z 5.7e+19)
       (fma
        (fma (fma (fma (fma 3.13060547623 z 11.1667541262) z t) z a) z b)
        (* (- y) -1.6453555072203998)
        x)
       t_1))))

C

double code(double x, double y, double z, double t, double a, double b) {
	double t_1 = fma((3.13060547623 + (-1.0 * ((36.52704169880642 + (-1.0 * ((457.9610022158428 + t) / z))) / z))), y, x);
	double tmp;
	if (z <= -0.06) {
		tmp = t_1;
	} else if (z <= 5.7e+19) {
		tmp = fma(fma(fma(fma(fma(3.13060547623, z, 11.1667541262), z, t), z, a), z, b), (-y * -1.6453555072203998), x);
	} else {
		tmp = t_1;
	}
	return tmp;
}

Julia

function code(x, y, z, t, a, b)
	t_1 = fma(Float64(3.13060547623 + Float64(-1.0 * Float64(Float64(36.52704169880642 + Float64(-1.0 * Float64(Float64(457.9610022158428 + t) / z))) / z))), y, x)
	tmp = 0.0
	if (z <= -0.06)
		tmp = t_1;
	elseif (z <= 5.7e+19)
		tmp = fma(fma(fma(fma(fma(3.13060547623, z, 11.1667541262), z, t), z, a), z, b), Float64(Float64(-y) * -1.6453555072203998), x);
	else
		tmp = t_1;
	end
	return tmp
end

Wolfram

code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(N[(3.13060547623 + N[(-1.0 * N[(N[(36.52704169880642 + N[(-1.0 * N[(N[(457.9610022158428 + t), $MachinePrecision] / z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * y + x), $MachinePrecision]}, If[LessEqual[z, -0.06], t$95$1, If[LessEqual[z, 5.7e+19], N[(N[(N[(N[(N[(3.13060547623 * z + 11.1667541262), $MachinePrecision] * z + t), $MachinePrecision] * z + a), $MachinePrecision] * z + b), $MachinePrecision] * N[((-y) * -1.6453555072203998), $MachinePrecision] + x), $MachinePrecision], t$95$1]]]

Derivation

1. Split input into 2 regimes

2. if z < -0.059999999999999998 or 5.7e19 < z

   1. Initial program 58.7%

      \[x + \frac{y \cdot \left(\left(\left(\left(z \cdot 3.13060547623 + 11.1667541262\right) \cdot z + t\right) \cdot z + a\right) \cdot z + b\right)}{\left(\left(\left(z + 15.234687407\right) \cdot z + 31.4690115749\right) \cdot z + 11.9400905721\right) \cdot z + 0.607771387771} \]

   3. Applied rewrites 60.9%

      \[\leadsto \color{blue}{\mathsf{fma}\left(\frac{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(3.13060547623, z, 11.1667541262\right), z, t\right), z, a\right), z, b\right)}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(z - -15.234687407, z, 31.4690115749\right), z, 11.9400905721\right), z, 0.607771387771\right)}, y, x\right)} \]

   4. Taylor expanded in z around -inf

      \[\leadsto \mathsf{fma}\left(\color{blue}{\frac{313060547623}{100000000000} + -1 \cdot \frac{\frac{3652704169880641883561}{100000000000000000000} + -1 \cdot \frac{\frac{45796100221584283915100827016327}{100000000000000000000000000000} + t}{z}}{z}}, y, x\right) \]
   5. Step-by-step derivation

      1. lower-+.f64 N/A

         \[\leadsto \mathsf{fma}\left(\frac{313060547623}{100000000000} + \color{blue}{-1 \cdot \frac{\frac{3652704169880641883561}{100000000000000000000} + -1 \cdot \frac{\frac{45796100221584283915100827016327}{100000000000000000000000000000} + t}{z}}{z}}, y, x\right) \]

      2. lower-*.f64 N/A

         \[\leadsto \mathsf{fma}\left(\frac{313060547623}{100000000000} + -1 \cdot \color{blue}{\frac{\frac{3652704169880641883561}{100000000000000000000} + -1 \cdot \frac{\frac{45796100221584283915100827016327}{100000000000000000000000000000} + t}{z}}{z}}, y, x\right) \]

      3. lower-/.f64 N/A

         \[\leadsto \mathsf{fma}\left(\frac{313060547623}{100000000000} + -1 \cdot \frac{\frac{3652704169880641883561}{100000000000000000000} + -1 \cdot \frac{\frac{45796100221584283915100827016327}{100000000000000000000000000000} + t}{z}}{\color{blue}{z}}, y, x\right) \]

      4. lower-+.f64 N/A

         \[\leadsto \mathsf{fma}\left(\frac{313060547623}{100000000000} + -1 \cdot \frac{\frac{3652704169880641883561}{100000000000000000000} + -1 \cdot \frac{\frac{45796100221584283915100827016327}{100000000000000000000000000000} + t}{z}}{z}, y, x\right) \]

      5. lower-*.f64 N/A

         \[\leadsto \mathsf{fma}\left(\frac{313060547623}{100000000000} + -1 \cdot \frac{\frac{3652704169880641883561}{100000000000000000000} + -1 \cdot \frac{\frac{45796100221584283915100827016327}{100000000000000000000000000000} + t}{z}}{z}, y, x\right) \]

      6. lower-/.f64 N/A

         \[\leadsto \mathsf{fma}\left(\frac{313060547623}{100000000000} + -1 \cdot \frac{\frac{3652704169880641883561}{100000000000000000000} + -1 \cdot \frac{\frac{45796100221584283915100827016327}{100000000000000000000000000000} + t}{z}}{z}, y, x\right) \]

      7. lower-+.f64 55.8

         \[\leadsto \mathsf{fma}\left(3.13060547623 + -1 \cdot \frac{36.52704169880642 + -1 \cdot \frac{457.9610022158428 + t}{z}}{z}, y, x\right) \]

   6. Applied rewrites 55.8%

      \[\leadsto \mathsf{fma}\left(\color{blue}{3.13060547623 + -1 \cdot \frac{36.52704169880642 + -1 \cdot \frac{457.9610022158428 + t}{z}}{z}}, y, x\right) \]

   if -0.059999999999999998 < z < 5.7e19

   1. Initial program 58.7%

      \[x + \frac{y \cdot \left(\left(\left(\left(z \cdot 3.13060547623 + 11.1667541262\right) \cdot z + t\right) \cdot z + a\right) \cdot z + b\right)}{\left(\left(\left(z + 15.234687407\right) \cdot z + 31.4690115749\right) \cdot z + 11.9400905721\right) \cdot z + 0.607771387771} \]

   3. Applied rewrites 60.6%

      \[\leadsto \color{blue}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(3.13060547623, z, 11.1667541262\right), z, t\right), z, a\right), z, b\right), \left(-y\right) \cdot \frac{-1}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(z - -15.234687407, z, 31.4690115749\right), z, 11.9400905721\right), z, 0.607771387771\right)}, x\right)} \]

   4. Taylor expanded in z around 0

      \[\leadsto \mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(\frac{313060547623}{100000000000}, z, \frac{55833770631}{5000000000}\right), z, t\right), z, a\right), z, b\right), \left(-y\right) \cdot \color{blue}{\frac{-1000000000000}{607771387771}}, x\right) \]
   5. Step-by-step derivation

   6. Applied rewrites 55.0%

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

Alternative 8: 93.8% accurate, 1.4× speedup

Math

\[\begin{array}{l} \\ \begin{array}{l} t_1 := \mathsf{fma}\left(3.13060547623 + -1 \cdot \frac{36.52704169880642 + -1 \cdot \frac{457.9610022158428 + t}{z}}{z}, y, x\right)\\ \mathbf{if}\;z \leq -13:\\ \;\;\;\;t\_1\\ \mathbf{elif}\;z \leq -2.5 \cdot 10^{-127}:\\ \;\;\;\;x - \frac{\mathsf{fma}\left(y \cdot \mathsf{fma}\left(t, z, a\right), z, b \cdot y\right)}{\mathsf{fma}\left(-11.9400905721, z, -0.607771387771\right)}\\ \mathbf{elif}\;z \leq 5.7 \cdot 10^{+19}:\\ \;\;\;\;\mathsf{fma}\left(\frac{\mathsf{fma}\left(a, z, b\right)}{\mathsf{fma}\left(\mathsf{fma}\left(31.4690115749, z, 11.9400905721\right), z, 0.607771387771\right)}, y, x\right)\\ \mathbf{else}:\\ \;\;\;\;t\_1\\ \end{array} \end{array} \]

FPCore

(FPCore (x y z t a b)
 :precision binary64
 (let* ((t_1
         (fma
          (+
           3.13060547623
           (*
            -1.0
            (/
             (+ 36.52704169880642 (* -1.0 (/ (+ 457.9610022158428 t) z)))
             z)))
          y
          x)))
   (if (<= z -13.0)
     t_1
     (if (<= z -2.5e-127)
       (-
        x
        (/
         (fma (* y (fma t z a)) z (* b y))
         (fma -11.9400905721 z -0.607771387771)))
       (if (<= z 5.7e+19)
         (fma
          (/
           (fma a z b)
           (fma (fma 31.4690115749 z 11.9400905721) z 0.607771387771))
          y
          x)
         t_1)))))

C

double code(double x, double y, double z, double t, double a, double b) {
	double t_1 = fma((3.13060547623 + (-1.0 * ((36.52704169880642 + (-1.0 * ((457.9610022158428 + t) / z))) / z))), y, x);
	double tmp;
	if (z <= -13.0) {
		tmp = t_1;
	} else if (z <= -2.5e-127) {
		tmp = x - (fma((y * fma(t, z, a)), z, (b * y)) / fma(-11.9400905721, z, -0.607771387771));
	} else if (z <= 5.7e+19) {
		tmp = fma((fma(a, z, b) / fma(fma(31.4690115749, z, 11.9400905721), z, 0.607771387771)), y, x);
	} else {
		tmp = t_1;
	}
	return tmp;
}

Julia

function code(x, y, z, t, a, b)
	t_1 = fma(Float64(3.13060547623 + Float64(-1.0 * Float64(Float64(36.52704169880642 + Float64(-1.0 * Float64(Float64(457.9610022158428 + t) / z))) / z))), y, x)
	tmp = 0.0
	if (z <= -13.0)
		tmp = t_1;
	elseif (z <= -2.5e-127)
		tmp = Float64(x - Float64(fma(Float64(y * fma(t, z, a)), z, Float64(b * y)) / fma(-11.9400905721, z, -0.607771387771)));
	elseif (z <= 5.7e+19)
		tmp = fma(Float64(fma(a, z, b) / fma(fma(31.4690115749, z, 11.9400905721), z, 0.607771387771)), y, x);
	else
		tmp = t_1;
	end
	return tmp
end

Wolfram

code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(N[(3.13060547623 + N[(-1.0 * N[(N[(36.52704169880642 + N[(-1.0 * N[(N[(457.9610022158428 + t), $MachinePrecision] / z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * y + x), $MachinePrecision]}, If[LessEqual[z, -13.0], t$95$1, If[LessEqual[z, -2.5e-127], N[(x - N[(N[(N[(y * N[(t * z + a), $MachinePrecision]), $MachinePrecision] * z + N[(b * y), $MachinePrecision]), $MachinePrecision] / N[(-11.9400905721 * z + -0.607771387771), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[z, 5.7e+19], N[(N[(N[(a * z + b), $MachinePrecision] / N[(N[(31.4690115749 * z + 11.9400905721), $MachinePrecision] * z + 0.607771387771), $MachinePrecision]), $MachinePrecision] * y + x), $MachinePrecision], t$95$1]]]]

Derivation

1. Split input into 3 regimes

2. if z < -13 or 5.7e19 < z

   1. Initial program 58.7%

      \[x + \frac{y \cdot \left(\left(\left(\left(z \cdot 3.13060547623 + 11.1667541262\right) \cdot z + t\right) \cdot z + a\right) \cdot z + b\right)}{\left(\left(\left(z + 15.234687407\right) \cdot z + 31.4690115749\right) \cdot z + 11.9400905721\right) \cdot z + 0.607771387771} \]

   3. Applied rewrites 60.9%

      \[\leadsto \color{blue}{\mathsf{fma}\left(\frac{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(3.13060547623, z, 11.1667541262\right), z, t\right), z, a\right), z, b\right)}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(z - -15.234687407, z, 31.4690115749\right), z, 11.9400905721\right), z, 0.607771387771\right)}, y, x\right)} \]

   4. Taylor expanded in z around -inf

      \[\leadsto \mathsf{fma}\left(\color{blue}{\frac{313060547623}{100000000000} + -1 \cdot \frac{\frac{3652704169880641883561}{100000000000000000000} + -1 \cdot \frac{\frac{45796100221584283915100827016327}{100000000000000000000000000000} + t}{z}}{z}}, y, x\right) \]
   5. Step-by-step derivation

      1. lower-+.f64 N/A

         \[\leadsto \mathsf{fma}\left(\frac{313060547623}{100000000000} + \color{blue}{-1 \cdot \frac{\frac{3652704169880641883561}{100000000000000000000} + -1 \cdot \frac{\frac{45796100221584283915100827016327}{100000000000000000000000000000} + t}{z}}{z}}, y, x\right) \]

      2. lower-*.f64 N/A

         \[\leadsto \mathsf{fma}\left(\frac{313060547623}{100000000000} + -1 \cdot \color{blue}{\frac{\frac{3652704169880641883561}{100000000000000000000} + -1 \cdot \frac{\frac{45796100221584283915100827016327}{100000000000000000000000000000} + t}{z}}{z}}, y, x\right) \]

      3. lower-/.f64 N/A

         \[\leadsto \mathsf{fma}\left(\frac{313060547623}{100000000000} + -1 \cdot \frac{\frac{3652704169880641883561}{100000000000000000000} + -1 \cdot \frac{\frac{45796100221584283915100827016327}{100000000000000000000000000000} + t}{z}}{\color{blue}{z}}, y, x\right) \]

      4. lower-+.f64 N/A

         \[\leadsto \mathsf{fma}\left(\frac{313060547623}{100000000000} + -1 \cdot \frac{\frac{3652704169880641883561}{100000000000000000000} + -1 \cdot \frac{\frac{45796100221584283915100827016327}{100000000000000000000000000000} + t}{z}}{z}, y, x\right) \]

      5. lower-*.f64 N/A

         \[\leadsto \mathsf{fma}\left(\frac{313060547623}{100000000000} + -1 \cdot \frac{\frac{3652704169880641883561}{100000000000000000000} + -1 \cdot \frac{\frac{45796100221584283915100827016327}{100000000000000000000000000000} + t}{z}}{z}, y, x\right) \]

      6. lower-/.f64 N/A

         \[\leadsto \mathsf{fma}\left(\frac{313060547623}{100000000000} + -1 \cdot \frac{\frac{3652704169880641883561}{100000000000000000000} + -1 \cdot \frac{\frac{45796100221584283915100827016327}{100000000000000000000000000000} + t}{z}}{z}, y, x\right) \]

      7. lower-+.f64 55.8

         \[\leadsto \mathsf{fma}\left(3.13060547623 + -1 \cdot \frac{36.52704169880642 + -1 \cdot \frac{457.9610022158428 + t}{z}}{z}, y, x\right) \]

   6. Applied rewrites 55.8%

      \[\leadsto \mathsf{fma}\left(\color{blue}{3.13060547623 + -1 \cdot \frac{36.52704169880642 + -1 \cdot \frac{457.9610022158428 + t}{z}}{z}}, y, x\right) \]

   if -13 < z < -2.4999999999999999e-127

   1. Initial program 58.7%

      \[x + \frac{y \cdot \left(\left(\left(\left(z \cdot 3.13060547623 + 11.1667541262\right) \cdot z + t\right) \cdot z + a\right) \cdot z + b\right)}{\left(\left(\left(z + 15.234687407\right) \cdot z + 31.4690115749\right) \cdot z + 11.9400905721\right) \cdot z + 0.607771387771} \]

   2. Taylor expanded in z around 0

      \[\leadsto x + \frac{\color{blue}{b \cdot y + z \cdot \left(a \cdot y + t \cdot \left(y \cdot z\right)\right)}}{\left(\left(\left(z + \frac{15234687407}{1000000000}\right) \cdot z + \frac{314690115749}{10000000000}\right) \cdot z + \frac{119400905721}{10000000000}\right) \cdot z + \frac{607771387771}{1000000000000}} \]
   3. Step-by-step derivation

      1. lower-fma.f64 N/A

         \[\leadsto x + \frac{\mathsf{fma}\left(b, \color{blue}{y}, z \cdot \left(a \cdot y + t \cdot \left(y \cdot z\right)\right)\right)}{\left(\left(\left(z + \frac{15234687407}{1000000000}\right) \cdot z + \frac{314690115749}{10000000000}\right) \cdot z + \frac{119400905721}{10000000000}\right) \cdot z + \frac{607771387771}{1000000000000}} \]

      2. lower-*.f64 N/A

         \[\leadsto x + \frac{\mathsf{fma}\left(b, y, z \cdot \left(a \cdot y + t \cdot \left(y \cdot z\right)\right)\right)}{\left(\left(\left(z + \frac{15234687407}{1000000000}\right) \cdot z + \frac{314690115749}{10000000000}\right) \cdot z + \frac{119400905721}{10000000000}\right) \cdot z + \frac{607771387771}{1000000000000}} \]

      3. lower-fma.f64 N/A

         \[\leadsto x + \frac{\mathsf{fma}\left(b, y, z \cdot \mathsf{fma}\left(a, y, t \cdot \left(y \cdot z\right)\right)\right)}{\left(\left(\left(z + \frac{15234687407}{1000000000}\right) \cdot z + \frac{314690115749}{10000000000}\right) \cdot z + \frac{119400905721}{10000000000}\right) \cdot z + \frac{607771387771}{1000000000000}} \]

      4. lower-*.f64 N/A

         \[\leadsto x + \frac{\mathsf{fma}\left(b, y, z \cdot \mathsf{fma}\left(a, y, t \cdot \left(y \cdot z\right)\right)\right)}{\left(\left(\left(z + \frac{15234687407}{1000000000}\right) \cdot z + \frac{314690115749}{10000000000}\right) \cdot z + \frac{119400905721}{10000000000}\right) \cdot z + \frac{607771387771}{1000000000000}} \]

      5. lower-*.f64 62.7

         \[\leadsto x + \frac{\mathsf{fma}\left(b, y, z \cdot \mathsf{fma}\left(a, y, t \cdot \left(y \cdot z\right)\right)\right)}{\left(\left(\left(z + 15.234687407\right) \cdot z + 31.4690115749\right) \cdot z + 11.9400905721\right) \cdot z + 0.607771387771} \]

   4. Applied rewrites 62.7%

      \[\leadsto x + \frac{\color{blue}{\mathsf{fma}\left(b, y, z \cdot \mathsf{fma}\left(a, y, t \cdot \left(y \cdot z\right)\right)\right)}}{\left(\left(\left(z + 15.234687407\right) \cdot z + 31.4690115749\right) \cdot z + 11.9400905721\right) \cdot z + 0.607771387771} \]

   6. Applied rewrites 61.6%

      \[\leadsto \color{blue}{x - \frac{\mathsf{fma}\left(y \cdot \mathsf{fma}\left(t, z, a\right), z, b \cdot y\right)}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(-15.234687407 - z, z, -31.4690115749\right), z, -11.9400905721\right), z, -0.607771387771\right)}} \]

   7. Taylor expanded in z around 0

      \[\leadsto x - \frac{\mathsf{fma}\left(y \cdot \mathsf{fma}\left(t, z, a\right), z, b \cdot y\right)}{\mathsf{fma}\left(\color{blue}{\frac{-119400905721}{10000000000}}, z, \frac{-607771387771}{1000000000000}\right)} \]
   8. Step-by-step derivation

   9. Applied rewrites 56.1%

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

   if -2.4999999999999999e-127 < z < 5.7e19

   1. Initial program 58.7%

      \[x + \frac{y \cdot \left(\left(\left(\left(z \cdot 3.13060547623 + 11.1667541262\right) \cdot z + t\right) \cdot z + a\right) \cdot z + b\right)}{\left(\left(\left(z + 15.234687407\right) \cdot z + 31.4690115749\right) \cdot z + 11.9400905721\right) \cdot z + 0.607771387771} \]

   3. Applied rewrites 60.9%

      \[\leadsto \color{blue}{\mathsf{fma}\left(\frac{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(3.13060547623, z, 11.1667541262\right), z, t\right), z, a\right), z, b\right)}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(z - -15.234687407, z, 31.4690115749\right), z, 11.9400905721\right), z, 0.607771387771\right)}, y, x\right)} \]

   4. Taylor expanded in z around 0

      \[\leadsto \mathsf{fma}\left(\frac{\mathsf{fma}\left(\color{blue}{a}, z, b\right)}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(z - \frac{-15234687407}{1000000000}, z, \frac{314690115749}{10000000000}\right), z, \frac{119400905721}{10000000000}\right), z, \frac{607771387771}{1000000000000}\right)}, y, x\right) \]
   5. Step-by-step derivation

   6. Applied rewrites 67.3%

      \[\leadsto \mathsf{fma}\left(\frac{\mathsf{fma}\left(\color{blue}{a}, z, b\right)}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(z - -15.234687407, z, 31.4690115749\right), z, 11.9400905721\right), z, 0.607771387771\right)}, y, x\right) \]

   7. Taylor expanded in z around 0

      \[\leadsto \mathsf{fma}\left(\frac{\mathsf{fma}\left(a, z, b\right)}{\mathsf{fma}\left(\mathsf{fma}\left(\color{blue}{\frac{314690115749}{10000000000}}, z, \frac{119400905721}{10000000000}\right), z, \frac{607771387771}{1000000000000}\right)}, y, x\right) \]
   8. Step-by-step derivation

   9. Applied rewrites 65.0%

      \[\leadsto \mathsf{fma}\left(\frac{\mathsf{fma}\left(a, z, b\right)}{\mathsf{fma}\left(\mathsf{fma}\left(\color{blue}{31.4690115749}, z, 11.9400905721\right), z, 0.607771387771\right)}, y, x\right) \]
3. Recombined 3 regimes into one program.
4. Add Preprocessing

Alternative 9: 93.5% accurate, 1.4× speedup

Math

\[\begin{array}{l} \\ \begin{array}{l} t_1 := \mathsf{fma}\left(3.13060547623 + -1 \cdot \frac{36.52704169880642 + -1 \cdot \frac{457.9610022158428 + t}{z}}{z}, y, x\right)\\ \mathbf{if}\;z \leq -13:\\ \;\;\;\;t\_1\\ \mathbf{elif}\;z \leq 5.7 \cdot 10^{+19}:\\ \;\;\;\;\mathsf{fma}\left(\frac{\mathsf{fma}\left(a, z, b\right)}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(15.234687407, z, 31.4690115749\right), z, 11.9400905721\right), z, 0.607771387771\right)}, y, x\right)\\ \mathbf{else}:\\ \;\;\;\;t\_1\\ \end{array} \end{array} \]

FPCore

(FPCore (x y z t a b)
 :precision binary64
 (let* ((t_1
         (fma
          (+
           3.13060547623
           (*
            -1.0
            (/
             (+ 36.52704169880642 (* -1.0 (/ (+ 457.9610022158428 t) z)))
             z)))
          y
          x)))
   (if (<= z -13.0)
     t_1
     (if (<= z 5.7e+19)
       (fma
        (/
         (fma a z b)
         (fma
          (fma (fma 15.234687407 z 31.4690115749) z 11.9400905721)
          z
          0.607771387771))
        y
        x)
       t_1))))

C

double code(double x, double y, double z, double t, double a, double b) {
	double t_1 = fma((3.13060547623 + (-1.0 * ((36.52704169880642 + (-1.0 * ((457.9610022158428 + t) / z))) / z))), y, x);
	double tmp;
	if (z <= -13.0) {
		tmp = t_1;
	} else if (z <= 5.7e+19) {
		tmp = fma((fma(a, z, b) / fma(fma(fma(15.234687407, z, 31.4690115749), z, 11.9400905721), z, 0.607771387771)), y, x);
	} else {
		tmp = t_1;
	}
	return tmp;
}

Julia

function code(x, y, z, t, a, b)
	t_1 = fma(Float64(3.13060547623 + Float64(-1.0 * Float64(Float64(36.52704169880642 + Float64(-1.0 * Float64(Float64(457.9610022158428 + t) / z))) / z))), y, x)
	tmp = 0.0
	if (z <= -13.0)
		tmp = t_1;
	elseif (z <= 5.7e+19)
		tmp = fma(Float64(fma(a, z, b) / fma(fma(fma(15.234687407, z, 31.4690115749), z, 11.9400905721), z, 0.607771387771)), y, x);
	else
		tmp = t_1;
	end
	return tmp
end

Wolfram

code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(N[(3.13060547623 + N[(-1.0 * N[(N[(36.52704169880642 + N[(-1.0 * N[(N[(457.9610022158428 + t), $MachinePrecision] / z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * y + x), $MachinePrecision]}, If[LessEqual[z, -13.0], t$95$1, If[LessEqual[z, 5.7e+19], N[(N[(N[(a * z + b), $MachinePrecision] / N[(N[(N[(15.234687407 * z + 31.4690115749), $MachinePrecision] * z + 11.9400905721), $MachinePrecision] * z + 0.607771387771), $MachinePrecision]), $MachinePrecision] * y + x), $MachinePrecision], t$95$1]]]

Derivation

1. Split input into 2 regimes

2. if z < -13 or 5.7e19 < z

   1. Initial program 58.7%

      \[x + \frac{y \cdot \left(\left(\left(\left(z \cdot 3.13060547623 + 11.1667541262\right) \cdot z + t\right) \cdot z + a\right) \cdot z + b\right)}{\left(\left(\left(z + 15.234687407\right) \cdot z + 31.4690115749\right) \cdot z + 11.9400905721\right) \cdot z + 0.607771387771} \]

   3. Applied rewrites 60.9%

      \[\leadsto \color{blue}{\mathsf{fma}\left(\frac{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(3.13060547623, z, 11.1667541262\right), z, t\right), z, a\right), z, b\right)}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(z - -15.234687407, z, 31.4690115749\right), z, 11.9400905721\right), z, 0.607771387771\right)}, y, x\right)} \]

   4. Taylor expanded in z around -inf

      \[\leadsto \mathsf{fma}\left(\color{blue}{\frac{313060547623}{100000000000} + -1 \cdot \frac{\frac{3652704169880641883561}{100000000000000000000} + -1 \cdot \frac{\frac{45796100221584283915100827016327}{100000000000000000000000000000} + t}{z}}{z}}, y, x\right) \]
   5. Step-by-step derivation

      1. lower-+.f64 N/A

         \[\leadsto \mathsf{fma}\left(\frac{313060547623}{100000000000} + \color{blue}{-1 \cdot \frac{\frac{3652704169880641883561}{100000000000000000000} + -1 \cdot \frac{\frac{45796100221584283915100827016327}{100000000000000000000000000000} + t}{z}}{z}}, y, x\right) \]

      2. lower-*.f64 N/A

         \[\leadsto \mathsf{fma}\left(\frac{313060547623}{100000000000} + -1 \cdot \color{blue}{\frac{\frac{3652704169880641883561}{100000000000000000000} + -1 \cdot \frac{\frac{45796100221584283915100827016327}{100000000000000000000000000000} + t}{z}}{z}}, y, x\right) \]

      3. lower-/.f64 N/A

         \[\leadsto \mathsf{fma}\left(\frac{313060547623}{100000000000} + -1 \cdot \frac{\frac{3652704169880641883561}{100000000000000000000} + -1 \cdot \frac{\frac{45796100221584283915100827016327}{100000000000000000000000000000} + t}{z}}{\color{blue}{z}}, y, x\right) \]

      4. lower-+.f64 N/A

         \[\leadsto \mathsf{fma}\left(\frac{313060547623}{100000000000} + -1 \cdot \frac{\frac{3652704169880641883561}{100000000000000000000} + -1 \cdot \frac{\frac{45796100221584283915100827016327}{100000000000000000000000000000} + t}{z}}{z}, y, x\right) \]

      5. lower-*.f64 N/A

         \[\leadsto \mathsf{fma}\left(\frac{313060547623}{100000000000} + -1 \cdot \frac{\frac{3652704169880641883561}{100000000000000000000} + -1 \cdot \frac{\frac{45796100221584283915100827016327}{100000000000000000000000000000} + t}{z}}{z}, y, x\right) \]

      6. lower-/.f64 N/A

         \[\leadsto \mathsf{fma}\left(\frac{313060547623}{100000000000} + -1 \cdot \frac{\frac{3652704169880641883561}{100000000000000000000} + -1 \cdot \frac{\frac{45796100221584283915100827016327}{100000000000000000000000000000} + t}{z}}{z}, y, x\right) \]

      7. lower-+.f64 55.8

         \[\leadsto \mathsf{fma}\left(3.13060547623 + -1 \cdot \frac{36.52704169880642 + -1 \cdot \frac{457.9610022158428 + t}{z}}{z}, y, x\right) \]

   6. Applied rewrites 55.8%

      \[\leadsto \mathsf{fma}\left(\color{blue}{3.13060547623 + -1 \cdot \frac{36.52704169880642 + -1 \cdot \frac{457.9610022158428 + t}{z}}{z}}, y, x\right) \]

   if -13 < z < 5.7e19

   1. Initial program 58.7%

      \[x + \frac{y \cdot \left(\left(\left(\left(z \cdot 3.13060547623 + 11.1667541262\right) \cdot z + t\right) \cdot z + a\right) \cdot z + b\right)}{\left(\left(\left(z + 15.234687407\right) \cdot z + 31.4690115749\right) \cdot z + 11.9400905721\right) \cdot z + 0.607771387771} \]

   3. Applied rewrites 60.9%

      \[\leadsto \color{blue}{\mathsf{fma}\left(\frac{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(3.13060547623, z, 11.1667541262\right), z, t\right), z, a\right), z, b\right)}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(z - -15.234687407, z, 31.4690115749\right), z, 11.9400905721\right), z, 0.607771387771\right)}, y, x\right)} \]

   4. Taylor expanded in z around 0

      \[\leadsto \mathsf{fma}\left(\frac{\mathsf{fma}\left(\color{blue}{a}, z, b\right)}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(z - \frac{-15234687407}{1000000000}, z, \frac{314690115749}{10000000000}\right), z, \frac{119400905721}{10000000000}\right), z, \frac{607771387771}{1000000000000}\right)}, y, x\right) \]
   5. Step-by-step derivation

   6. Applied rewrites 67.3%

      \[\leadsto \mathsf{fma}\left(\frac{\mathsf{fma}\left(\color{blue}{a}, z, b\right)}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(z - -15.234687407, z, 31.4690115749\right), z, 11.9400905721\right), z, 0.607771387771\right)}, y, x\right) \]

   7. Taylor expanded in z around 0

      \[\leadsto \mathsf{fma}\left(\frac{\mathsf{fma}\left(a, z, b\right)}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(\color{blue}{\frac{15234687407}{1000000000}}, z, \frac{314690115749}{10000000000}\right), z, \frac{119400905721}{10000000000}\right), z, \frac{607771387771}{1000000000000}\right)}, y, x\right) \]
   8. Step-by-step derivation

   9. Applied rewrites 64.9%

      \[\leadsto \mathsf{fma}\left(\frac{\mathsf{fma}\left(a, z, b\right)}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(\color{blue}{15.234687407}, z, 31.4690115749\right), z, 11.9400905721\right), z, 0.607771387771\right)}, y, x\right) \]
3. Recombined 2 regimes into one program.
4. Add Preprocessing

Alternative 10: 93.4% accurate, 1.5× speedup

Math

\[\begin{array}{l} \\ \begin{array}{l} t_1 := \mathsf{fma}\left(3.13060547623 + -1 \cdot \frac{36.52704169880642 + -1 \cdot \frac{457.9610022158428 + t}{z}}{z}, y, x\right)\\ \mathbf{if}\;z \leq -13:\\ \;\;\;\;t\_1\\ \mathbf{elif}\;z \leq 5.7 \cdot 10^{+19}:\\ \;\;\;\;\mathsf{fma}\left(\frac{\mathsf{fma}\left(a, z, b\right)}{\mathsf{fma}\left(11.9400905721, z, 0.607771387771\right)}, y, x\right)\\ \mathbf{else}:\\ \;\;\;\;t\_1\\ \end{array} \end{array} \]

FPCore

(FPCore (x y z t a b)
 :precision binary64
 (let* ((t_1
         (fma
          (+
           3.13060547623
           (*
            -1.0
            (/
             (+ 36.52704169880642 (* -1.0 (/ (+ 457.9610022158428 t) z)))
             z)))
          y
          x)))
   (if (<= z -13.0)
     t_1
     (if (<= z 5.7e+19)
       (fma (/ (fma a z b) (fma 11.9400905721 z 0.607771387771)) y x)
       t_1))))

C

double code(double x, double y, double z, double t, double a, double b) {
	double t_1 = fma((3.13060547623 + (-1.0 * ((36.52704169880642 + (-1.0 * ((457.9610022158428 + t) / z))) / z))), y, x);
	double tmp;
	if (z <= -13.0) {
		tmp = t_1;
	} else if (z <= 5.7e+19) {
		tmp = fma((fma(a, z, b) / fma(11.9400905721, z, 0.607771387771)), y, x);
	} else {
		tmp = t_1;
	}
	return tmp;
}

Julia

function code(x, y, z, t, a, b)
	t_1 = fma(Float64(3.13060547623 + Float64(-1.0 * Float64(Float64(36.52704169880642 + Float64(-1.0 * Float64(Float64(457.9610022158428 + t) / z))) / z))), y, x)
	tmp = 0.0
	if (z <= -13.0)
		tmp = t_1;
	elseif (z <= 5.7e+19)
		tmp = fma(Float64(fma(a, z, b) / fma(11.9400905721, z, 0.607771387771)), y, x);
	else
		tmp = t_1;
	end
	return tmp
end

Wolfram

code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(N[(3.13060547623 + N[(-1.0 * N[(N[(36.52704169880642 + N[(-1.0 * N[(N[(457.9610022158428 + t), $MachinePrecision] / z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * y + x), $MachinePrecision]}, If[LessEqual[z, -13.0], t$95$1, If[LessEqual[z, 5.7e+19], N[(N[(N[(a * z + b), $MachinePrecision] / N[(11.9400905721 * z + 0.607771387771), $MachinePrecision]), $MachinePrecision] * y + x), $MachinePrecision], t$95$1]]]

Derivation

1. Split input into 2 regimes

2. if z < -13 or 5.7e19 < z

   1. Initial program 58.7%

      \[x + \frac{y \cdot \left(\left(\left(\left(z \cdot 3.13060547623 + 11.1667541262\right) \cdot z + t\right) \cdot z + a\right) \cdot z + b\right)}{\left(\left(\left(z + 15.234687407\right) \cdot z + 31.4690115749\right) \cdot z + 11.9400905721\right) \cdot z + 0.607771387771} \]

   3. Applied rewrites 60.9%

      \[\leadsto \color{blue}{\mathsf{fma}\left(\frac{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(3.13060547623, z, 11.1667541262\right), z, t\right), z, a\right), z, b\right)}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(z - -15.234687407, z, 31.4690115749\right), z, 11.9400905721\right), z, 0.607771387771\right)}, y, x\right)} \]

   4. Taylor expanded in z around -inf

      \[\leadsto \mathsf{fma}\left(\color{blue}{\frac{313060547623}{100000000000} + -1 \cdot \frac{\frac{3652704169880641883561}{100000000000000000000} + -1 \cdot \frac{\frac{45796100221584283915100827016327}{100000000000000000000000000000} + t}{z}}{z}}, y, x\right) \]
   5. Step-by-step derivation

      1. lower-+.f64 N/A

         \[\leadsto \mathsf{fma}\left(\frac{313060547623}{100000000000} + \color{blue}{-1 \cdot \frac{\frac{3652704169880641883561}{100000000000000000000} + -1 \cdot \frac{\frac{45796100221584283915100827016327}{100000000000000000000000000000} + t}{z}}{z}}, y, x\right) \]

      2. lower-*.f64 N/A

         \[\leadsto \mathsf{fma}\left(\frac{313060547623}{100000000000} + -1 \cdot \color{blue}{\frac{\frac{3652704169880641883561}{100000000000000000000} + -1 \cdot \frac{\frac{45796100221584283915100827016327}{100000000000000000000000000000} + t}{z}}{z}}, y, x\right) \]

      3. lower-/.f64 N/A

         \[\leadsto \mathsf{fma}\left(\frac{313060547623}{100000000000} + -1 \cdot \frac{\frac{3652704169880641883561}{100000000000000000000} + -1 \cdot \frac{\frac{45796100221584283915100827016327}{100000000000000000000000000000} + t}{z}}{\color{blue}{z}}, y, x\right) \]

      4. lower-+.f64 N/A

         \[\leadsto \mathsf{fma}\left(\frac{313060547623}{100000000000} + -1 \cdot \frac{\frac{3652704169880641883561}{100000000000000000000} + -1 \cdot \frac{\frac{45796100221584283915100827016327}{100000000000000000000000000000} + t}{z}}{z}, y, x\right) \]

      5. lower-*.f64 N/A

         \[\leadsto \mathsf{fma}\left(\frac{313060547623}{100000000000} + -1 \cdot \frac{\frac{3652704169880641883561}{100000000000000000000} + -1 \cdot \frac{\frac{45796100221584283915100827016327}{100000000000000000000000000000} + t}{z}}{z}, y, x\right) \]

      6. lower-/.f64 N/A

         \[\leadsto \mathsf{fma}\left(\frac{313060547623}{100000000000} + -1 \cdot \frac{\frac{3652704169880641883561}{100000000000000000000} + -1 \cdot \frac{\frac{45796100221584283915100827016327}{100000000000000000000000000000} + t}{z}}{z}, y, x\right) \]

      7. lower-+.f64 55.8

         \[\leadsto \mathsf{fma}\left(3.13060547623 + -1 \cdot \frac{36.52704169880642 + -1 \cdot \frac{457.9610022158428 + t}{z}}{z}, y, x\right) \]

   6. Applied rewrites 55.8%

      \[\leadsto \mathsf{fma}\left(\color{blue}{3.13060547623 + -1 \cdot \frac{36.52704169880642 + -1 \cdot \frac{457.9610022158428 + t}{z}}{z}}, y, x\right) \]

   if -13 < z < 5.7e19

   1. Initial program 58.7%

      \[x + \frac{y \cdot \left(\left(\left(\left(z \cdot 3.13060547623 + 11.1667541262\right) \cdot z + t\right) \cdot z + a\right) \cdot z + b\right)}{\left(\left(\left(z + 15.234687407\right) \cdot z + 31.4690115749\right) \cdot z + 11.9400905721\right) \cdot z + 0.607771387771} \]

   3. Applied rewrites 60.9%

      \[\leadsto \color{blue}{\mathsf{fma}\left(\frac{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(3.13060547623, z, 11.1667541262\right), z, t\right), z, a\right), z, b\right)}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(z - -15.234687407, z, 31.4690115749\right), z, 11.9400905721\right), z, 0.607771387771\right)}, y, x\right)} \]

   4. Taylor expanded in z around 0

      \[\leadsto \mathsf{fma}\left(\frac{\mathsf{fma}\left(\color{blue}{a}, z, b\right)}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(z - \frac{-15234687407}{1000000000}, z, \frac{314690115749}{10000000000}\right), z, \frac{119400905721}{10000000000}\right), z, \frac{607771387771}{1000000000000}\right)}, y, x\right) \]
   5. Step-by-step derivation

   6. Applied rewrites 67.3%

      \[\leadsto \mathsf{fma}\left(\frac{\mathsf{fma}\left(\color{blue}{a}, z, b\right)}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(z - -15.234687407, z, 31.4690115749\right), z, 11.9400905721\right), z, 0.607771387771\right)}, y, x\right) \]

   7. Taylor expanded in z around 0

      \[\leadsto \mathsf{fma}\left(\frac{\mathsf{fma}\left(a, z, b\right)}{\mathsf{fma}\left(\color{blue}{\frac{119400905721}{10000000000}}, z, \frac{607771387771}{1000000000000}\right)}, y, x\right) \]
   8. Step-by-step derivation

   9. Applied rewrites 63.4%

      \[\leadsto \mathsf{fma}\left(\frac{\mathsf{fma}\left(a, z, b\right)}{\mathsf{fma}\left(\color{blue}{11.9400905721}, z, 0.607771387771\right)}, y, x\right) \]
3. Recombined 2 regimes into one program.
4. Add Preprocessing

Alternative 11: 90.2% accurate, 1.9× speedup

Math

\[\begin{array}{l} \\ \begin{array}{l} \mathbf{if}\;z \leq -0.7:\\ \;\;\;\;x - \mathsf{fma}\left(-3.13060547623, y, -1 \cdot \frac{11.1667541262 \cdot y - 47.69379582500642 \cdot y}{z}\right)\\ \mathbf{elif}\;z \leq 6.2 \cdot 10^{+27}:\\ \;\;\;\;\mathsf{fma}\left(\frac{\mathsf{fma}\left(a, z, b\right)}{\mathsf{fma}\left(11.9400905721, z, 0.607771387771\right)}, y, x\right)\\ \mathbf{else}:\\ \;\;\;\;\mathsf{fma}\left(3.13060547623, y, x\right)\\ \end{array} \end{array} \]

FPCore

(FPCore (x y z t a b)
 :precision binary64
 (if (<= z -0.7)
   (-
    x
    (fma
     -3.13060547623
     y
     (* -1.0 (/ (- (* 11.1667541262 y) (* 47.69379582500642 y)) z))))
   (if (<= z 6.2e+27)
     (fma (/ (fma a z b) (fma 11.9400905721 z 0.607771387771)) y x)
     (fma 3.13060547623 y x))))

C

double code(double x, double y, double z, double t, double a, double b) {
	double tmp;
	if (z <= -0.7) {
		tmp = x - fma(-3.13060547623, y, (-1.0 * (((11.1667541262 * y) - (47.69379582500642 * y)) / z)));
	} else if (z <= 6.2e+27) {
		tmp = fma((fma(a, z, b) / fma(11.9400905721, z, 0.607771387771)), y, x);
	} else {
		tmp = fma(3.13060547623, y, x);
	}
	return tmp;
}

Julia

function code(x, y, z, t, a, b)
	tmp = 0.0
	if (z <= -0.7)
		tmp = Float64(x - fma(-3.13060547623, y, Float64(-1.0 * Float64(Float64(Float64(11.1667541262 * y) - Float64(47.69379582500642 * y)) / z))));
	elseif (z <= 6.2e+27)
		tmp = fma(Float64(fma(a, z, b) / fma(11.9400905721, z, 0.607771387771)), y, x);
	else
		tmp = fma(3.13060547623, y, x);
	end
	return tmp
end

Wolfram

code[x_, y_, z_, t_, a_, b_] := If[LessEqual[z, -0.7], N[(x - N[(-3.13060547623 * y + N[(-1.0 * N[(N[(N[(11.1667541262 * y), $MachinePrecision] - N[(47.69379582500642 * y), $MachinePrecision]), $MachinePrecision] / z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[z, 6.2e+27], N[(N[(N[(a * z + b), $MachinePrecision] / N[(11.9400905721 * z + 0.607771387771), $MachinePrecision]), $MachinePrecision] * y + x), $MachinePrecision], N[(3.13060547623 * y + x), $MachinePrecision]]]

Derivation

1. Split input into 3 regimes

2. if z < -0.69999999999999996

   1. Initial program 58.7%

      \[x + \frac{y \cdot \left(\left(\left(\left(z \cdot 3.13060547623 + 11.1667541262\right) \cdot z + t\right) \cdot z + a\right) \cdot z + b\right)}{\left(\left(\left(z + 15.234687407\right) \cdot z + 31.4690115749\right) \cdot z + 11.9400905721\right) \cdot z + 0.607771387771} \]

   2. Taylor expanded in z around 0

      \[\leadsto x + \frac{\color{blue}{b \cdot y + z \cdot \left(a \cdot y + t \cdot \left(y \cdot z\right)\right)}}{\left(\left(\left(z + \frac{15234687407}{1000000000}\right) \cdot z + \frac{314690115749}{10000000000}\right) \cdot z + \frac{119400905721}{10000000000}\right) \cdot z + \frac{607771387771}{1000000000000}} \]
   3. Step-by-step derivation

      1. lower-fma.f64 N/A

         \[\leadsto x + \frac{\mathsf{fma}\left(b, \color{blue}{y}, z \cdot \left(a \cdot y + t \cdot \left(y \cdot z\right)\right)\right)}{\left(\left(\left(z + \frac{15234687407}{1000000000}\right) \cdot z + \frac{314690115749}{10000000000}\right) \cdot z + \frac{119400905721}{10000000000}\right) \cdot z + \frac{607771387771}{1000000000000}} \]

      2. lower-*.f64 N/A

         \[\leadsto x + \frac{\mathsf{fma}\left(b, y, z \cdot \left(a \cdot y + t \cdot \left(y \cdot z\right)\right)\right)}{\left(\left(\left(z + \frac{15234687407}{1000000000}\right) \cdot z + \frac{314690115749}{10000000000}\right) \cdot z + \frac{119400905721}{10000000000}\right) \cdot z + \frac{607771387771}{1000000000000}} \]

      3. lower-fma.f64 N/A

         \[\leadsto x + \frac{\mathsf{fma}\left(b, y, z \cdot \mathsf{fma}\left(a, y, t \cdot \left(y \cdot z\right)\right)\right)}{\left(\left(\left(z + \frac{15234687407}{1000000000}\right) \cdot z + \frac{314690115749}{10000000000}\right) \cdot z + \frac{119400905721}{10000000000}\right) \cdot z + \frac{607771387771}{1000000000000}} \]

      4. lower-*.f64 N/A

         \[\leadsto x + \frac{\mathsf{fma}\left(b, y, z \cdot \mathsf{fma}\left(a, y, t \cdot \left(y \cdot z\right)\right)\right)}{\left(\left(\left(z + \frac{15234687407}{1000000000}\right) \cdot z + \frac{314690115749}{10000000000}\right) \cdot z + \frac{119400905721}{10000000000}\right) \cdot z + \frac{607771387771}{1000000000000}} \]

      5. lower-*.f64 62.7

         \[\leadsto x + \frac{\mathsf{fma}\left(b, y, z \cdot \mathsf{fma}\left(a, y, t \cdot \left(y \cdot z\right)\right)\right)}{\left(\left(\left(z + 15.234687407\right) \cdot z + 31.4690115749\right) \cdot z + 11.9400905721\right) \cdot z + 0.607771387771} \]

   4. Applied rewrites 62.7%

      \[\leadsto x + \frac{\color{blue}{\mathsf{fma}\left(b, y, z \cdot \mathsf{fma}\left(a, y, t \cdot \left(y \cdot z\right)\right)\right)}}{\left(\left(\left(z + 15.234687407\right) \cdot z + 31.4690115749\right) \cdot z + 11.9400905721\right) \cdot z + 0.607771387771} \]

   6. Applied rewrites 61.6%

      \[\leadsto \color{blue}{x - \frac{\mathsf{fma}\left(y \cdot \mathsf{fma}\left(t, z, a\right), z, b \cdot y\right)}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(-15.234687407 - z, z, -31.4690115749\right), z, -11.9400905721\right), z, -0.607771387771\right)}} \]

   7. Taylor expanded in z around -inf

      \[\leadsto x - \color{blue}{\left(\frac{-313060547623}{100000000000} \cdot y + -1 \cdot \frac{\frac{55833770631}{5000000000} \cdot y - \frac{4769379582500641883561}{100000000000000000000} \cdot y}{z}\right)} \]
   8. Step-by-step derivation

      1. lower-fma.f64 N/A

         \[\leadsto x - \mathsf{fma}\left(\frac{-313060547623}{100000000000}, \color{blue}{y}, -1 \cdot \frac{\frac{55833770631}{5000000000} \cdot y - \frac{4769379582500641883561}{100000000000000000000} \cdot y}{z}\right) \]

      2. lower-*.f64 N/A

         \[\leadsto x - \mathsf{fma}\left(\frac{-313060547623}{100000000000}, y, -1 \cdot \frac{\frac{55833770631}{5000000000} \cdot y - \frac{4769379582500641883561}{100000000000000000000} \cdot y}{z}\right) \]

      3. lower-/.f64 N/A

         \[\leadsto x - \mathsf{fma}\left(\frac{-313060547623}{100000000000}, y, -1 \cdot \frac{\frac{55833770631}{5000000000} \cdot y - \frac{4769379582500641883561}{100000000000000000000} \cdot y}{z}\right) \]

      4. lower--.f64 N/A

         \[\leadsto x - \mathsf{fma}\left(\frac{-313060547623}{100000000000}, y, -1 \cdot \frac{\frac{55833770631}{5000000000} \cdot y - \frac{4769379582500641883561}{100000000000000000000} \cdot y}{z}\right) \]

      5. lower-*.f64 N/A

         \[\leadsto x - \mathsf{fma}\left(\frac{-313060547623}{100000000000}, y, -1 \cdot \frac{\frac{55833770631}{5000000000} \cdot y - \frac{4769379582500641883561}{100000000000000000000} \cdot y}{z}\right) \]

      6. lower-*.f64 59.0

         \[\leadsto x - \mathsf{fma}\left(-3.13060547623, y, -1 \cdot \frac{11.1667541262 \cdot y - 47.69379582500642 \cdot y}{z}\right) \]

   9. Applied rewrites 59.0%

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

   if -0.69999999999999996 < z < 6.19999999999999992e27

   1. Initial program 58.7%

      \[x + \frac{y \cdot \left(\left(\left(\left(z \cdot 3.13060547623 + 11.1667541262\right) \cdot z + t\right) \cdot z + a\right) \cdot z + b\right)}{\left(\left(\left(z + 15.234687407\right) \cdot z + 31.4690115749\right) \cdot z + 11.9400905721\right) \cdot z + 0.607771387771} \]

   3. Applied rewrites 60.9%

      \[\leadsto \color{blue}{\mathsf{fma}\left(\frac{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(3.13060547623, z, 11.1667541262\right), z, t\right), z, a\right), z, b\right)}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(z - -15.234687407, z, 31.4690115749\right), z, 11.9400905721\right), z, 0.607771387771\right)}, y, x\right)} \]

   4. Taylor expanded in z around 0

      \[\leadsto \mathsf{fma}\left(\frac{\mathsf{fma}\left(\color{blue}{a}, z, b\right)}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(z - \frac{-15234687407}{1000000000}, z, \frac{314690115749}{10000000000}\right), z, \frac{119400905721}{10000000000}\right), z, \frac{607771387771}{1000000000000}\right)}, y, x\right) \]
   5. Step-by-step derivation

   6. Applied rewrites 67.3%

      \[\leadsto \mathsf{fma}\left(\frac{\mathsf{fma}\left(\color{blue}{a}, z, b\right)}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(z - -15.234687407, z, 31.4690115749\right), z, 11.9400905721\right), z, 0.607771387771\right)}, y, x\right) \]

   7. Taylor expanded in z around 0

      \[\leadsto \mathsf{fma}\left(\frac{\mathsf{fma}\left(a, z, b\right)}{\mathsf{fma}\left(\color{blue}{\frac{119400905721}{10000000000}}, z, \frac{607771387771}{1000000000000}\right)}, y, x\right) \]
   8. Step-by-step derivation

   9. Applied rewrites 63.4%

      \[\leadsto \mathsf{fma}\left(\frac{\mathsf{fma}\left(a, z, b\right)}{\mathsf{fma}\left(\color{blue}{11.9400905721}, z, 0.607771387771\right)}, y, x\right) \]

   if 6.19999999999999992e27 < z

   1. Initial program 58.7%

      \[x + \frac{y \cdot \left(\left(\left(\left(z \cdot 3.13060547623 + 11.1667541262\right) \cdot z + t\right) \cdot z + a\right) \cdot z + b\right)}{\left(\left(\left(z + 15.234687407\right) \cdot z + 31.4690115749\right) \cdot z + 11.9400905721\right) \cdot z + 0.607771387771} \]

   3. Applied rewrites 60.9%

      \[\leadsto \color{blue}{\mathsf{fma}\left(\frac{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(3.13060547623, z, 11.1667541262\right), z, t\right), z, a\right), z, b\right)}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(z - -15.234687407, z, 31.4690115749\right), z, 11.9400905721\right), z, 0.607771387771\right)}, y, x\right)} \]

   4. Taylor expanded in z around 0

      \[\leadsto \mathsf{fma}\left(\frac{\mathsf{fma}\left(\color{blue}{a}, z, b\right)}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(z - \frac{-15234687407}{1000000000}, z, \frac{314690115749}{10000000000}\right), z, \frac{119400905721}{10000000000}\right), z, \frac{607771387771}{1000000000000}\right)}, y, x\right) \]
   5. Step-by-step derivation

   6. Applied rewrites 67.3%

      \[\leadsto \mathsf{fma}\left(\frac{\mathsf{fma}\left(\color{blue}{a}, z, b\right)}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(z - -15.234687407, z, 31.4690115749\right), z, 11.9400905721\right), z, 0.607771387771\right)}, y, x\right) \]

   7. Taylor expanded in z around inf

      \[\leadsto \mathsf{fma}\left(\color{blue}{\frac{313060547623}{100000000000}}, y, x\right) \]
   8. Step-by-step derivation

   9. Applied rewrites 62.8%

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

Alternative 12: 90.2% accurate, 2.0× speedup

Math

\[\begin{array}{l} \\ \begin{array}{l} \mathbf{if}\;z \leq -0.7:\\ \;\;\;\;\mathsf{fma}\left(3.13060547623 - 36.52704169880642 \cdot \frac{1}{z}, y, x\right)\\ \mathbf{elif}\;z \leq 6.2 \cdot 10^{+27}:\\ \;\;\;\;\mathsf{fma}\left(\frac{\mathsf{fma}\left(a, z, b\right)}{\mathsf{fma}\left(11.9400905721, z, 0.607771387771\right)}, y, x\right)\\ \mathbf{else}:\\ \;\;\;\;\mathsf{fma}\left(3.13060547623, y, x\right)\\ \end{array} \end{array} \]

FPCore

(FPCore (x y z t a b)
 :precision binary64
 (if (<= z -0.7)
   (fma (- 3.13060547623 (* 36.52704169880642 (/ 1.0 z))) y x)
   (if (<= z 6.2e+27)
     (fma (/ (fma a z b) (fma 11.9400905721 z 0.607771387771)) y x)
     (fma 3.13060547623 y x))))

C

double code(double x, double y, double z, double t, double a, double b) {
	double tmp;
	if (z <= -0.7) {
		tmp = fma((3.13060547623 - (36.52704169880642 * (1.0 / z))), y, x);
	} else if (z <= 6.2e+27) {
		tmp = fma((fma(a, z, b) / fma(11.9400905721, z, 0.607771387771)), y, x);
	} else {
		tmp = fma(3.13060547623, y, x);
	}
	return tmp;
}

Julia

function code(x, y, z, t, a, b)
	tmp = 0.0
	if (z <= -0.7)
		tmp = fma(Float64(3.13060547623 - Float64(36.52704169880642 * Float64(1.0 / z))), y, x);
	elseif (z <= 6.2e+27)
		tmp = fma(Float64(fma(a, z, b) / fma(11.9400905721, z, 0.607771387771)), y, x);
	else
		tmp = fma(3.13060547623, y, x);
	end
	return tmp
end

Wolfram

code[x_, y_, z_, t_, a_, b_] := If[LessEqual[z, -0.7], N[(N[(3.13060547623 - N[(36.52704169880642 * N[(1.0 / z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * y + x), $MachinePrecision], If[LessEqual[z, 6.2e+27], N[(N[(N[(a * z + b), $MachinePrecision] / N[(11.9400905721 * z + 0.607771387771), $MachinePrecision]), $MachinePrecision] * y + x), $MachinePrecision], N[(3.13060547623 * y + x), $MachinePrecision]]]

Derivation

1. Split input into 3 regimes

2. if z < -0.69999999999999996

   1. Initial program 58.7%

      \[x + \frac{y \cdot \left(\left(\left(\left(z \cdot 3.13060547623 + 11.1667541262\right) \cdot z + t\right) \cdot z + a\right) \cdot z + b\right)}{\left(\left(\left(z + 15.234687407\right) \cdot z + 31.4690115749\right) \cdot z + 11.9400905721\right) \cdot z + 0.607771387771} \]

   3. Applied rewrites 60.9%

      \[\leadsto \color{blue}{\mathsf{fma}\left(\frac{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(3.13060547623, z, 11.1667541262\right), z, t\right), z, a\right), z, b\right)}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(z - -15.234687407, z, 31.4690115749\right), z, 11.9400905721\right), z, 0.607771387771\right)}, y, x\right)} \]

   4. Taylor expanded in z around 0

      \[\leadsto \mathsf{fma}\left(\frac{\mathsf{fma}\left(\color{blue}{a}, z, b\right)}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(z - \frac{-15234687407}{1000000000}, z, \frac{314690115749}{10000000000}\right), z, \frac{119400905721}{10000000000}\right), z, \frac{607771387771}{1000000000000}\right)}, y, x\right) \]
   5. Step-by-step derivation

   6. Applied rewrites 67.3%

      \[\leadsto \mathsf{fma}\left(\frac{\mathsf{fma}\left(\color{blue}{a}, z, b\right)}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(z - -15.234687407, z, 31.4690115749\right), z, 11.9400905721\right), z, 0.607771387771\right)}, y, x\right) \]

   7. Taylor expanded in z around inf

      \[\leadsto \mathsf{fma}\left(\color{blue}{\frac{313060547623}{100000000000} - \frac{3652704169880641883561}{100000000000000000000} \cdot \frac{1}{z}}, y, x\right) \]
   8. Step-by-step derivation

      1. lower--.f64 N/A

         \[\leadsto \mathsf{fma}\left(\frac{313060547623}{100000000000} - \color{blue}{\frac{3652704169880641883561}{100000000000000000000} \cdot \frac{1}{z}}, y, x\right) \]

      2. lower-*.f64 N/A

         \[\leadsto \mathsf{fma}\left(\frac{313060547623}{100000000000} - \frac{3652704169880641883561}{100000000000000000000} \cdot \color{blue}{\frac{1}{z}}, y, x\right) \]

      3. lower-/.f64 59.0

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

   9. Applied rewrites 59.0%

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

   if -0.69999999999999996 < z < 6.19999999999999992e27

   1. Initial program 58.7%

      \[x + \frac{y \cdot \left(\left(\left(\left(z \cdot 3.13060547623 + 11.1667541262\right) \cdot z + t\right) \cdot z + a\right) \cdot z + b\right)}{\left(\left(\left(z + 15.234687407\right) \cdot z + 31.4690115749\right) \cdot z + 11.9400905721\right) \cdot z + 0.607771387771} \]

   3. Applied rewrites 60.9%

      \[\leadsto \color{blue}{\mathsf{fma}\left(\frac{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(3.13060547623, z, 11.1667541262\right), z, t\right), z, a\right), z, b\right)}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(z - -15.234687407, z, 31.4690115749\right), z, 11.9400905721\right), z, 0.607771387771\right)}, y, x\right)} \]

   4. Taylor expanded in z around 0

      \[\leadsto \mathsf{fma}\left(\frac{\mathsf{fma}\left(\color{blue}{a}, z, b\right)}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(z - \frac{-15234687407}{1000000000}, z, \frac{314690115749}{10000000000}\right), z, \frac{119400905721}{10000000000}\right), z, \frac{607771387771}{1000000000000}\right)}, y, x\right) \]
   5. Step-by-step derivation

   6. Applied rewrites 67.3%

      \[\leadsto \mathsf{fma}\left(\frac{\mathsf{fma}\left(\color{blue}{a}, z, b\right)}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(z - -15.234687407, z, 31.4690115749\right), z, 11.9400905721\right), z, 0.607771387771\right)}, y, x\right) \]

   7. Taylor expanded in z around 0

      \[\leadsto \mathsf{fma}\left(\frac{\mathsf{fma}\left(a, z, b\right)}{\mathsf{fma}\left(\color{blue}{\frac{119400905721}{10000000000}}, z, \frac{607771387771}{1000000000000}\right)}, y, x\right) \]
   8. Step-by-step derivation

   9. Applied rewrites 63.4%

      \[\leadsto \mathsf{fma}\left(\frac{\mathsf{fma}\left(a, z, b\right)}{\mathsf{fma}\left(\color{blue}{11.9400905721}, z, 0.607771387771\right)}, y, x\right) \]

   if 6.19999999999999992e27 < z

   1. Initial program 58.7%

      \[x + \frac{y \cdot \left(\left(\left(\left(z \cdot 3.13060547623 + 11.1667541262\right) \cdot z + t\right) \cdot z + a\right) \cdot z + b\right)}{\left(\left(\left(z + 15.234687407\right) \cdot z + 31.4690115749\right) \cdot z + 11.9400905721\right) \cdot z + 0.607771387771} \]

   3. Applied rewrites 60.9%

      \[\leadsto \color{blue}{\mathsf{fma}\left(\frac{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(3.13060547623, z, 11.1667541262\right), z, t\right), z, a\right), z, b\right)}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(z - -15.234687407, z, 31.4690115749\right), z, 11.9400905721\right), z, 0.607771387771\right)}, y, x\right)} \]

   4. Taylor expanded in z around 0

      \[\leadsto \mathsf{fma}\left(\frac{\mathsf{fma}\left(\color{blue}{a}, z, b\right)}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(z - \frac{-15234687407}{1000000000}, z, \frac{314690115749}{10000000000}\right), z, \frac{119400905721}{10000000000}\right), z, \frac{607771387771}{1000000000000}\right)}, y, x\right) \]
   5. Step-by-step derivation

   6. Applied rewrites 67.3%

      \[\leadsto \mathsf{fma}\left(\frac{\mathsf{fma}\left(\color{blue}{a}, z, b\right)}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(z - -15.234687407, z, 31.4690115749\right), z, 11.9400905721\right), z, 0.607771387771\right)}, y, x\right) \]

   7. Taylor expanded in z around inf

      \[\leadsto \mathsf{fma}\left(\color{blue}{\frac{313060547623}{100000000000}}, y, x\right) \]
   8. Step-by-step derivation

   9. Applied rewrites 62.8%

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

Alternative 13: 90.1% accurate, 2.1× speedup

Math

\[\begin{array}{l} \\ \begin{array}{l} \mathbf{if}\;z \leq -3100:\\ \;\;\;\;\mathsf{fma}\left(3.13060547623 - 36.52704169880642 \cdot \frac{1}{z}, y, x\right)\\ \mathbf{elif}\;z \leq 6.2 \cdot 10^{+27}:\\ \;\;\;\;\mathsf{fma}\left(\frac{1}{0.607771387771}, y \cdot \mathsf{fma}\left(a, z, b\right), x\right)\\ \mathbf{else}:\\ \;\;\;\;\mathsf{fma}\left(3.13060547623, y, x\right)\\ \end{array} \end{array} \]

FPCore

(FPCore (x y z t a b)
 :precision binary64
 (if (<= z -3100.0)
   (fma (- 3.13060547623 (* 36.52704169880642 (/ 1.0 z))) y x)
   (if (<= z 6.2e+27)
     (fma (/ 1.0 0.607771387771) (* y (fma a z b)) x)
     (fma 3.13060547623 y x))))

C

double code(double x, double y, double z, double t, double a, double b) {
	double tmp;
	if (z <= -3100.0) {
		tmp = fma((3.13060547623 - (36.52704169880642 * (1.0 / z))), y, x);
	} else if (z <= 6.2e+27) {
		tmp = fma((1.0 / 0.607771387771), (y * fma(a, z, b)), x);
	} else {
		tmp = fma(3.13060547623, y, x);
	}
	return tmp;
}

Julia

function code(x, y, z, t, a, b)
	tmp = 0.0
	if (z <= -3100.0)
		tmp = fma(Float64(3.13060547623 - Float64(36.52704169880642 * Float64(1.0 / z))), y, x);
	elseif (z <= 6.2e+27)
		tmp = fma(Float64(1.0 / 0.607771387771), Float64(y * fma(a, z, b)), x);
	else
		tmp = fma(3.13060547623, y, x);
	end
	return tmp
end

Wolfram

code[x_, y_, z_, t_, a_, b_] := If[LessEqual[z, -3100.0], N[(N[(3.13060547623 - N[(36.52704169880642 * N[(1.0 / z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * y + x), $MachinePrecision], If[LessEqual[z, 6.2e+27], N[(N[(1.0 / 0.607771387771), $MachinePrecision] * N[(y * N[(a * z + b), $MachinePrecision]), $MachinePrecision] + x), $MachinePrecision], N[(3.13060547623 * y + x), $MachinePrecision]]]

Derivation

1. Split input into 3 regimes

2. if z < -3100

   1. Initial program 58.7%

      \[x + \frac{y \cdot \left(\left(\left(\left(z \cdot 3.13060547623 + 11.1667541262\right) \cdot z + t\right) \cdot z + a\right) \cdot z + b\right)}{\left(\left(\left(z + 15.234687407\right) \cdot z + 31.4690115749\right) \cdot z + 11.9400905721\right) \cdot z + 0.607771387771} \]

   3. Applied rewrites 60.9%

      \[\leadsto \color{blue}{\mathsf{fma}\left(\frac{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(3.13060547623, z, 11.1667541262\right), z, t\right), z, a\right), z, b\right)}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(z - -15.234687407, z, 31.4690115749\right), z, 11.9400905721\right), z, 0.607771387771\right)}, y, x\right)} \]

   4. Taylor expanded in z around 0

      \[\leadsto \mathsf{fma}\left(\frac{\mathsf{fma}\left(\color{blue}{a}, z, b\right)}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(z - \frac{-15234687407}{1000000000}, z, \frac{314690115749}{10000000000}\right), z, \frac{119400905721}{10000000000}\right), z, \frac{607771387771}{1000000000000}\right)}, y, x\right) \]
   5. Step-by-step derivation

   6. Applied rewrites 67.3%

      \[\leadsto \mathsf{fma}\left(\frac{\mathsf{fma}\left(\color{blue}{a}, z, b\right)}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(z - -15.234687407, z, 31.4690115749\right), z, 11.9400905721\right), z, 0.607771387771\right)}, y, x\right) \]

   7. Taylor expanded in z around inf

      \[\leadsto \mathsf{fma}\left(\color{blue}{\frac{313060547623}{100000000000} - \frac{3652704169880641883561}{100000000000000000000} \cdot \frac{1}{z}}, y, x\right) \]
   8. Step-by-step derivation

      1. lower--.f64 N/A

         \[\leadsto \mathsf{fma}\left(\frac{313060547623}{100000000000} - \color{blue}{\frac{3652704169880641883561}{100000000000000000000} \cdot \frac{1}{z}}, y, x\right) \]

      2. lower-*.f64 N/A

         \[\leadsto \mathsf{fma}\left(\frac{313060547623}{100000000000} - \frac{3652704169880641883561}{100000000000000000000} \cdot \color{blue}{\frac{1}{z}}, y, x\right) \]

      3. lower-/.f64 59.0

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

   9. Applied rewrites 59.0%

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

   if -3100 < z < 6.19999999999999992e27

   1. Initial program 58.7%

      \[x + \frac{y \cdot \left(\left(\left(\left(z \cdot 3.13060547623 + 11.1667541262\right) \cdot z + t\right) \cdot z + a\right) \cdot z + b\right)}{\left(\left(\left(z + 15.234687407\right) \cdot z + 31.4690115749\right) \cdot z + 11.9400905721\right) \cdot z + 0.607771387771} \]

   3. Applied rewrites 63.3%

      \[\leadsto \color{blue}{\mathsf{fma}\left(y, z \cdot \frac{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(3.13060547623, z, 11.1667541262\right), z, t\right), z, a\right)}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(z - -15.234687407, z, 31.4690115749\right), z, 11.9400905721\right), z, 0.607771387771\right)}, \mathsf{fma}\left(b, \frac{y}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(z - -15.234687407, z, 31.4690115749\right), z, 11.9400905721\right), z, 0.607771387771\right)}, x\right)\right)} \]

   4. Taylor expanded in z around 0

      \[\leadsto \mathsf{fma}\left(y, z \cdot \frac{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(\frac{313060547623}{100000000000}, z, \frac{55833770631}{5000000000}\right), z, t\right), z, a\right)}{\color{blue}{\frac{607771387771}{1000000000000}}}, \mathsf{fma}\left(b, \frac{y}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(z - \frac{-15234687407}{1000000000}, z, \frac{314690115749}{10000000000}\right), z, \frac{119400905721}{10000000000}\right), z, \frac{607771387771}{1000000000000}\right)}, x\right)\right) \]
   5. Step-by-step derivation

   6. Applied rewrites 55.1%

      \[\leadsto \mathsf{fma}\left(y, z \cdot \frac{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(3.13060547623, z, 11.1667541262\right), z, t\right), z, a\right)}{\color{blue}{0.607771387771}}, \mathsf{fma}\left(b, \frac{y}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(z - -15.234687407, z, 31.4690115749\right), z, 11.9400905721\right), z, 0.607771387771\right)}, x\right)\right) \]

   7. Taylor expanded in z around 0

      \[\leadsto \mathsf{fma}\left(y, z \cdot \frac{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(\frac{313060547623}{100000000000}, z, \frac{55833770631}{5000000000}\right), z, t\right), z, a\right)}{\frac{607771387771}{1000000000000}}, \mathsf{fma}\left(b, \frac{y}{\color{blue}{\frac{607771387771}{1000000000000}}}, x\right)\right) \]
   8. Step-by-step derivation

   9. Applied rewrites 54.3%

      \[\leadsto \mathsf{fma}\left(y, z \cdot \frac{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(3.13060547623, z, 11.1667541262\right), z, t\right), z, a\right)}{0.607771387771}, \mathsf{fma}\left(b, \frac{y}{\color{blue}{0.607771387771}}, x\right)\right) \]

   10. Taylor expanded in z around 0

       \[\leadsto \mathsf{fma}\left(y, z \cdot \frac{\color{blue}{a}}{\frac{607771387771}{1000000000000}}, \mathsf{fma}\left(b, \frac{y}{\frac{607771387771}{1000000000000}}, x\right)\right) \]
   11. Step-by-step derivation

   12. Applied rewrites 59.3%

       \[\leadsto \mathsf{fma}\left(y, z \cdot \frac{\color{blue}{a}}{0.607771387771}, \mathsf{fma}\left(b, \frac{y}{0.607771387771}, x\right)\right) \]
   13. Step-by-step derivation

       1. lift-fma.f64 N/A

          \[\leadsto \color{blue}{y \cdot \left(z \cdot \frac{a}{\frac{607771387771}{1000000000000}}\right) + \mathsf{fma}\left(b, \frac{y}{\frac{607771387771}{1000000000000}}, x\right)} \]

       2. lift-fma.f64 N/A

          \[\leadsto y \cdot \left(z \cdot \frac{a}{\frac{607771387771}{1000000000000}}\right) + \color{blue}{\left(b \cdot \frac{y}{\frac{607771387771}{1000000000000}} + x\right)} \]

       3. associate-+r+ N/A

          \[\leadsto \color{blue}{\left(y \cdot \left(z \cdot \frac{a}{\frac{607771387771}{1000000000000}}\right) + b \cdot \frac{y}{\frac{607771387771}{1000000000000}}\right) + x} \]

   14. Applied rewrites 59.7%

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

   if 6.19999999999999992e27 < z

   1. Initial program 58.7%

      \[x + \frac{y \cdot \left(\left(\left(\left(z \cdot 3.13060547623 + 11.1667541262\right) \cdot z + t\right) \cdot z + a\right) \cdot z + b\right)}{\left(\left(\left(z + 15.234687407\right) \cdot z + 31.4690115749\right) \cdot z + 11.9400905721\right) \cdot z + 0.607771387771} \]

   3. Applied rewrites 60.9%

      \[\leadsto \color{blue}{\mathsf{fma}\left(\frac{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(3.13060547623, z, 11.1667541262\right), z, t\right), z, a\right), z, b\right)}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(z - -15.234687407, z, 31.4690115749\right), z, 11.9400905721\right), z, 0.607771387771\right)}, y, x\right)} \]

   4. Taylor expanded in z around 0

      \[\leadsto \mathsf{fma}\left(\frac{\mathsf{fma}\left(\color{blue}{a}, z, b\right)}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(z - \frac{-15234687407}{1000000000}, z, \frac{314690115749}{10000000000}\right), z, \frac{119400905721}{10000000000}\right), z, \frac{607771387771}{1000000000000}\right)}, y, x\right) \]
   5. Step-by-step derivation

   6. Applied rewrites 67.3%

      \[\leadsto \mathsf{fma}\left(\frac{\mathsf{fma}\left(\color{blue}{a}, z, b\right)}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(z - -15.234687407, z, 31.4690115749\right), z, 11.9400905721\right), z, 0.607771387771\right)}, y, x\right) \]

   7. Taylor expanded in z around inf

      \[\leadsto \mathsf{fma}\left(\color{blue}{\frac{313060547623}{100000000000}}, y, x\right) \]
   8. Step-by-step derivation

   9. Applied rewrites 62.8%

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

Alternative 14: 90.1% accurate, 2.3× speedup

Math

\[\begin{array}{l} \\ \begin{array}{l} \mathbf{if}\;z \leq -3100:\\ \;\;\;\;\mathsf{fma}\left(3.13060547623 - 36.52704169880642 \cdot \frac{1}{z}, y, x\right)\\ \mathbf{elif}\;z \leq 6.2 \cdot 10^{+27}:\\ \;\;\;\;\frac{y \cdot \mathsf{fma}\left(a, z, b\right)}{0.607771387771} + x\\ \mathbf{else}:\\ \;\;\;\;\mathsf{fma}\left(3.13060547623, y, x\right)\\ \end{array} \end{array} \]

FPCore

(FPCore (x y z t a b)
 :precision binary64
 (if (<= z -3100.0)
   (fma (- 3.13060547623 (* 36.52704169880642 (/ 1.0 z))) y x)
   (if (<= z 6.2e+27)
     (+ (/ (* y (fma a z b)) 0.607771387771) x)
     (fma 3.13060547623 y x))))

C

double code(double x, double y, double z, double t, double a, double b) {
	double tmp;
	if (z <= -3100.0) {
		tmp = fma((3.13060547623 - (36.52704169880642 * (1.0 / z))), y, x);
	} else if (z <= 6.2e+27) {
		tmp = ((y * fma(a, z, b)) / 0.607771387771) + x;
	} else {
		tmp = fma(3.13060547623, y, x);
	}
	return tmp;
}

Julia

function code(x, y, z, t, a, b)
	tmp = 0.0
	if (z <= -3100.0)
		tmp = fma(Float64(3.13060547623 - Float64(36.52704169880642 * Float64(1.0 / z))), y, x);
	elseif (z <= 6.2e+27)
		tmp = Float64(Float64(Float64(y * fma(a, z, b)) / 0.607771387771) + x);
	else
		tmp = fma(3.13060547623, y, x);
	end
	return tmp
end

Wolfram

code[x_, y_, z_, t_, a_, b_] := If[LessEqual[z, -3100.0], N[(N[(3.13060547623 - N[(36.52704169880642 * N[(1.0 / z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * y + x), $MachinePrecision], If[LessEqual[z, 6.2e+27], N[(N[(N[(y * N[(a * z + b), $MachinePrecision]), $MachinePrecision] / 0.607771387771), $MachinePrecision] + x), $MachinePrecision], N[(3.13060547623 * y + x), $MachinePrecision]]]

Derivation

1. Split input into 3 regimes

2. if z < -3100

   1. Initial program 58.7%

      \[x + \frac{y \cdot \left(\left(\left(\left(z \cdot 3.13060547623 + 11.1667541262\right) \cdot z + t\right) \cdot z + a\right) \cdot z + b\right)}{\left(\left(\left(z + 15.234687407\right) \cdot z + 31.4690115749\right) \cdot z + 11.9400905721\right) \cdot z + 0.607771387771} \]

   3. Applied rewrites 60.9%

      \[\leadsto \color{blue}{\mathsf{fma}\left(\frac{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(3.13060547623, z, 11.1667541262\right), z, t\right), z, a\right), z, b\right)}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(z - -15.234687407, z, 31.4690115749\right), z, 11.9400905721\right), z, 0.607771387771\right)}, y, x\right)} \]

   4. Taylor expanded in z around 0

      \[\leadsto \mathsf{fma}\left(\frac{\mathsf{fma}\left(\color{blue}{a}, z, b\right)}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(z - \frac{-15234687407}{1000000000}, z, \frac{314690115749}{10000000000}\right), z, \frac{119400905721}{10000000000}\right), z, \frac{607771387771}{1000000000000}\right)}, y, x\right) \]
   5. Step-by-step derivation

   6. Applied rewrites 67.3%

      \[\leadsto \mathsf{fma}\left(\frac{\mathsf{fma}\left(\color{blue}{a}, z, b\right)}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(z - -15.234687407, z, 31.4690115749\right), z, 11.9400905721\right), z, 0.607771387771\right)}, y, x\right) \]

   7. Taylor expanded in z around inf

      \[\leadsto \mathsf{fma}\left(\color{blue}{\frac{313060547623}{100000000000} - \frac{3652704169880641883561}{100000000000000000000} \cdot \frac{1}{z}}, y, x\right) \]
   8. Step-by-step derivation

      1. lower--.f64 N/A

         \[\leadsto \mathsf{fma}\left(\frac{313060547623}{100000000000} - \color{blue}{\frac{3652704169880641883561}{100000000000000000000} \cdot \frac{1}{z}}, y, x\right) \]

      2. lower-*.f64 N/A

         \[\leadsto \mathsf{fma}\left(\frac{313060547623}{100000000000} - \frac{3652704169880641883561}{100000000000000000000} \cdot \color{blue}{\frac{1}{z}}, y, x\right) \]

      3. lower-/.f64 59.0

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

   9. Applied rewrites 59.0%

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

   if -3100 < z < 6.19999999999999992e27

   1. Initial program 58.7%

      \[x + \frac{y \cdot \left(\left(\left(\left(z \cdot 3.13060547623 + 11.1667541262\right) \cdot z + t\right) \cdot z + a\right) \cdot z + b\right)}{\left(\left(\left(z + 15.234687407\right) \cdot z + 31.4690115749\right) \cdot z + 11.9400905721\right) \cdot z + 0.607771387771} \]

   3. Applied rewrites 63.3%

      \[\leadsto \color{blue}{\mathsf{fma}\left(y, z \cdot \frac{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(3.13060547623, z, 11.1667541262\right), z, t\right), z, a\right)}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(z - -15.234687407, z, 31.4690115749\right), z, 11.9400905721\right), z, 0.607771387771\right)}, \mathsf{fma}\left(b, \frac{y}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(z - -15.234687407, z, 31.4690115749\right), z, 11.9400905721\right), z, 0.607771387771\right)}, x\right)\right)} \]

   4. Taylor expanded in z around 0

      \[\leadsto \mathsf{fma}\left(y, z \cdot \frac{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(\frac{313060547623}{100000000000}, z, \frac{55833770631}{5000000000}\right), z, t\right), z, a\right)}{\color{blue}{\frac{607771387771}{1000000000000}}}, \mathsf{fma}\left(b, \frac{y}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(z - \frac{-15234687407}{1000000000}, z, \frac{314690115749}{10000000000}\right), z, \frac{119400905721}{10000000000}\right), z, \frac{607771387771}{1000000000000}\right)}, x\right)\right) \]
   5. Step-by-step derivation

   6. Applied rewrites 55.1%

      \[\leadsto \mathsf{fma}\left(y, z \cdot \frac{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(3.13060547623, z, 11.1667541262\right), z, t\right), z, a\right)}{\color{blue}{0.607771387771}}, \mathsf{fma}\left(b, \frac{y}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(z - -15.234687407, z, 31.4690115749\right), z, 11.9400905721\right), z, 0.607771387771\right)}, x\right)\right) \]

   7. Taylor expanded in z around 0

      \[\leadsto \mathsf{fma}\left(y, z \cdot \frac{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(\frac{313060547623}{100000000000}, z, \frac{55833770631}{5000000000}\right), z, t\right), z, a\right)}{\frac{607771387771}{1000000000000}}, \mathsf{fma}\left(b, \frac{y}{\color{blue}{\frac{607771387771}{1000000000000}}}, x\right)\right) \]
   8. Step-by-step derivation

   9. Applied rewrites 54.3%

      \[\leadsto \mathsf{fma}\left(y, z \cdot \frac{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(3.13060547623, z, 11.1667541262\right), z, t\right), z, a\right)}{0.607771387771}, \mathsf{fma}\left(b, \frac{y}{\color{blue}{0.607771387771}}, x\right)\right) \]

   10. Taylor expanded in z around 0

       \[\leadsto \mathsf{fma}\left(y, z \cdot \frac{\color{blue}{a}}{\frac{607771387771}{1000000000000}}, \mathsf{fma}\left(b, \frac{y}{\frac{607771387771}{1000000000000}}, x\right)\right) \]
   11. Step-by-step derivation

   12. Applied rewrites 59.3%

       \[\leadsto \mathsf{fma}\left(y, z \cdot \frac{\color{blue}{a}}{0.607771387771}, \mathsf{fma}\left(b, \frac{y}{0.607771387771}, x\right)\right) \]
   13. Step-by-step derivation

       1. lift-fma.f64 N/A

          \[\leadsto \color{blue}{y \cdot \left(z \cdot \frac{a}{\frac{607771387771}{1000000000000}}\right) + \mathsf{fma}\left(b, \frac{y}{\frac{607771387771}{1000000000000}}, x\right)} \]

       2. lift-fma.f64 N/A

          \[\leadsto y \cdot \left(z \cdot \frac{a}{\frac{607771387771}{1000000000000}}\right) + \color{blue}{\left(b \cdot \frac{y}{\frac{607771387771}{1000000000000}} + x\right)} \]

       3. associate-+r+ N/A

          \[\leadsto \color{blue}{\left(y \cdot \left(z \cdot \frac{a}{\frac{607771387771}{1000000000000}}\right) + b \cdot \frac{y}{\frac{607771387771}{1000000000000}}\right) + x} \]

       4. lower-+.f64 N/A

          \[\leadsto \color{blue}{\left(y \cdot \left(z \cdot \frac{a}{\frac{607771387771}{1000000000000}}\right) + b \cdot \frac{y}{\frac{607771387771}{1000000000000}}\right) + x} \]

   14. Applied rewrites 59.7%

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

   if 6.19999999999999992e27 < z

   1. Initial program 58.7%

      \[x + \frac{y \cdot \left(\left(\left(\left(z \cdot 3.13060547623 + 11.1667541262\right) \cdot z + t\right) \cdot z + a\right) \cdot z + b\right)}{\left(\left(\left(z + 15.234687407\right) \cdot z + 31.4690115749\right) \cdot z + 11.9400905721\right) \cdot z + 0.607771387771} \]

   3. Applied rewrites 60.9%

      \[\leadsto \color{blue}{\mathsf{fma}\left(\frac{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(3.13060547623, z, 11.1667541262\right), z, t\right), z, a\right), z, b\right)}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(z - -15.234687407, z, 31.4690115749\right), z, 11.9400905721\right), z, 0.607771387771\right)}, y, x\right)} \]

   4. Taylor expanded in z around 0

      \[\leadsto \mathsf{fma}\left(\frac{\mathsf{fma}\left(\color{blue}{a}, z, b\right)}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(z - \frac{-15234687407}{1000000000}, z, \frac{314690115749}{10000000000}\right), z, \frac{119400905721}{10000000000}\right), z, \frac{607771387771}{1000000000000}\right)}, y, x\right) \]
   5. Step-by-step derivation

   6. Applied rewrites 67.3%

      \[\leadsto \mathsf{fma}\left(\frac{\mathsf{fma}\left(\color{blue}{a}, z, b\right)}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(z - -15.234687407, z, 31.4690115749\right), z, 11.9400905721\right), z, 0.607771387771\right)}, y, x\right) \]

   7. Taylor expanded in z around inf

      \[\leadsto \mathsf{fma}\left(\color{blue}{\frac{313060547623}{100000000000}}, y, x\right) \]
   8. Step-by-step derivation

   9. Applied rewrites 62.8%

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

Alternative 15: 83.2% accurate, 2.4× speedup

Math

\[\begin{array}{l} \\ \begin{array}{l} \mathbf{if}\;z \leq -0.55:\\ \;\;\;\;\mathsf{fma}\left(3.13060547623 - 36.52704169880642 \cdot \frac{1}{z}, y, x\right)\\ \mathbf{elif}\;z \leq 6.2 \cdot 10^{+27}:\\ \;\;\;\;\mathsf{fma}\left(\frac{y}{\mathsf{fma}\left(11.9400905721, z, 0.607771387771\right)}, b, x\right)\\ \mathbf{else}:\\ \;\;\;\;\mathsf{fma}\left(3.13060547623, y, x\right)\\ \end{array} \end{array} \]

FPCore

(FPCore (x y z t a b)
 :precision binary64
 (if (<= z -0.55)
   (fma (- 3.13060547623 (* 36.52704169880642 (/ 1.0 z))) y x)
   (if (<= z 6.2e+27)
     (fma (/ y (fma 11.9400905721 z 0.607771387771)) b x)
     (fma 3.13060547623 y x))))

C

double code(double x, double y, double z, double t, double a, double b) {
	double tmp;
	if (z <= -0.55) {
		tmp = fma((3.13060547623 - (36.52704169880642 * (1.0 / z))), y, x);
	} else if (z <= 6.2e+27) {
		tmp = fma((y / fma(11.9400905721, z, 0.607771387771)), b, x);
	} else {
		tmp = fma(3.13060547623, y, x);
	}
	return tmp;
}

Julia

function code(x, y, z, t, a, b)
	tmp = 0.0
	if (z <= -0.55)
		tmp = fma(Float64(3.13060547623 - Float64(36.52704169880642 * Float64(1.0 / z))), y, x);
	elseif (z <= 6.2e+27)
		tmp = fma(Float64(y / fma(11.9400905721, z, 0.607771387771)), b, x);
	else
		tmp = fma(3.13060547623, y, x);
	end
	return tmp
end

Wolfram

code[x_, y_, z_, t_, a_, b_] := If[LessEqual[z, -0.55], N[(N[(3.13060547623 - N[(36.52704169880642 * N[(1.0 / z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * y + x), $MachinePrecision], If[LessEqual[z, 6.2e+27], N[(N[(y / N[(11.9400905721 * z + 0.607771387771), $MachinePrecision]), $MachinePrecision] * b + x), $MachinePrecision], N[(3.13060547623 * y + x), $MachinePrecision]]]

Derivation

1. Split input into 3 regimes

2. if z < -0.55000000000000004

   1. Initial program 58.7%

      \[x + \frac{y \cdot \left(\left(\left(\left(z \cdot 3.13060547623 + 11.1667541262\right) \cdot z + t\right) \cdot z + a\right) \cdot z + b\right)}{\left(\left(\left(z + 15.234687407\right) \cdot z + 31.4690115749\right) \cdot z + 11.9400905721\right) \cdot z + 0.607771387771} \]

   3. Applied rewrites 60.9%

      \[\leadsto \color{blue}{\mathsf{fma}\left(\frac{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(3.13060547623, z, 11.1667541262\right), z, t\right), z, a\right), z, b\right)}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(z - -15.234687407, z, 31.4690115749\right), z, 11.9400905721\right), z, 0.607771387771\right)}, y, x\right)} \]

   4. Taylor expanded in z around 0

      \[\leadsto \mathsf{fma}\left(\frac{\mathsf{fma}\left(\color{blue}{a}, z, b\right)}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(z - \frac{-15234687407}{1000000000}, z, \frac{314690115749}{10000000000}\right), z, \frac{119400905721}{10000000000}\right), z, \frac{607771387771}{1000000000000}\right)}, y, x\right) \]
   5. Step-by-step derivation

   6. Applied rewrites 67.3%

      \[\leadsto \mathsf{fma}\left(\frac{\mathsf{fma}\left(\color{blue}{a}, z, b\right)}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(z - -15.234687407, z, 31.4690115749\right), z, 11.9400905721\right), z, 0.607771387771\right)}, y, x\right) \]

   7. Taylor expanded in z around inf

      \[\leadsto \mathsf{fma}\left(\color{blue}{\frac{313060547623}{100000000000} - \frac{3652704169880641883561}{100000000000000000000} \cdot \frac{1}{z}}, y, x\right) \]
   8. Step-by-step derivation

      1. lower--.f64 N/A

         \[\leadsto \mathsf{fma}\left(\frac{313060547623}{100000000000} - \color{blue}{\frac{3652704169880641883561}{100000000000000000000} \cdot \frac{1}{z}}, y, x\right) \]

      2. lower-*.f64 N/A

         \[\leadsto \mathsf{fma}\left(\frac{313060547623}{100000000000} - \frac{3652704169880641883561}{100000000000000000000} \cdot \color{blue}{\frac{1}{z}}, y, x\right) \]

      3. lower-/.f64 59.0

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

   9. Applied rewrites 59.0%

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

   if -0.55000000000000004 < z < 6.19999999999999992e27

   1. Initial program 58.7%

      \[x + \frac{y \cdot \left(\left(\left(\left(z \cdot 3.13060547623 + 11.1667541262\right) \cdot z + t\right) \cdot z + a\right) \cdot z + b\right)}{\left(\left(\left(z + 15.234687407\right) \cdot z + 31.4690115749\right) \cdot z + 11.9400905721\right) \cdot z + 0.607771387771} \]

   2. Taylor expanded in z around 0

      \[\leadsto x + \frac{y \cdot \color{blue}{b}}{\left(\left(\left(z + \frac{15234687407}{1000000000}\right) \cdot z + \frac{314690115749}{10000000000}\right) \cdot z + \frac{119400905721}{10000000000}\right) \cdot z + \frac{607771387771}{1000000000000}} \]
   3. Step-by-step derivation

   4. Applied rewrites 64.2%

      \[\leadsto x + \frac{y \cdot \color{blue}{b}}{\left(\left(\left(z + 15.234687407\right) \cdot z + 31.4690115749\right) \cdot z + 11.9400905721\right) \cdot z + 0.607771387771} \]

   6. Applied rewrites 65.2%

      \[\leadsto \color{blue}{\mathsf{fma}\left(\frac{y}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(z - -15.234687407, z, 31.4690115749\right), z, 11.9400905721\right), z, 0.607771387771\right)}, b, x\right)} \]

   7. Taylor expanded in z around 0

      \[\leadsto \mathsf{fma}\left(\frac{y}{\mathsf{fma}\left(\color{blue}{\frac{119400905721}{10000000000}}, z, \frac{607771387771}{1000000000000}\right)}, b, x\right) \]
   8. Step-by-step derivation

   9. Applied rewrites 63.3%

      \[\leadsto \mathsf{fma}\left(\frac{y}{\mathsf{fma}\left(\color{blue}{11.9400905721}, z, 0.607771387771\right)}, b, x\right) \]

   if 6.19999999999999992e27 < z

   1. Initial program 58.7%

      \[x + \frac{y \cdot \left(\left(\left(\left(z \cdot 3.13060547623 + 11.1667541262\right) \cdot z + t\right) \cdot z + a\right) \cdot z + b\right)}{\left(\left(\left(z + 15.234687407\right) \cdot z + 31.4690115749\right) \cdot z + 11.9400905721\right) \cdot z + 0.607771387771} \]

   3. Applied rewrites 60.9%

      \[\leadsto \color{blue}{\mathsf{fma}\left(\frac{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(3.13060547623, z, 11.1667541262\right), z, t\right), z, a\right), z, b\right)}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(z - -15.234687407, z, 31.4690115749\right), z, 11.9400905721\right), z, 0.607771387771\right)}, y, x\right)} \]

   4. Taylor expanded in z around 0

      \[\leadsto \mathsf{fma}\left(\frac{\mathsf{fma}\left(\color{blue}{a}, z, b\right)}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(z - \frac{-15234687407}{1000000000}, z, \frac{314690115749}{10000000000}\right), z, \frac{119400905721}{10000000000}\right), z, \frac{607771387771}{1000000000000}\right)}, y, x\right) \]
   5. Step-by-step derivation

   6. Applied rewrites 67.3%

      \[\leadsto \mathsf{fma}\left(\frac{\mathsf{fma}\left(\color{blue}{a}, z, b\right)}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(z - -15.234687407, z, 31.4690115749\right), z, 11.9400905721\right), z, 0.607771387771\right)}, y, x\right) \]

   7. Taylor expanded in z around inf

      \[\leadsto \mathsf{fma}\left(\color{blue}{\frac{313060547623}{100000000000}}, y, x\right) \]
   8. Step-by-step derivation

   9. Applied rewrites 62.8%

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

Alternative 16: 83.2% accurate, 2.4× speedup

Math

\[\begin{array}{l} \\ \begin{array}{l} \mathbf{if}\;z \leq -0.55:\\ \;\;\;\;\mathsf{fma}\left(3.13060547623, y, x\right)\\ \mathbf{elif}\;z \leq 6.2 \cdot 10^{+27}:\\ \;\;\;\;\mathsf{fma}\left(\frac{y}{\mathsf{fma}\left(11.9400905721, z, 0.607771387771\right)}, b, x\right)\\ \mathbf{else}:\\ \;\;\;\;\mathsf{fma}\left(3.13060547623, y, x\right)\\ \end{array} \end{array} \]

FPCore

(FPCore (x y z t a b)
 :precision binary64
 (if (<= z -0.55)
   (fma 3.13060547623 y x)
   (if (<= z 6.2e+27)
     (fma (/ y (fma 11.9400905721 z 0.607771387771)) b x)
     (fma 3.13060547623 y x))))

C

double code(double x, double y, double z, double t, double a, double b) {
	double tmp;
	if (z <= -0.55) {
		tmp = fma(3.13060547623, y, x);
	} else if (z <= 6.2e+27) {
		tmp = fma((y / fma(11.9400905721, z, 0.607771387771)), b, x);
	} else {
		tmp = fma(3.13060547623, y, x);
	}
	return tmp;
}

Julia

function code(x, y, z, t, a, b)
	tmp = 0.0
	if (z <= -0.55)
		tmp = fma(3.13060547623, y, x);
	elseif (z <= 6.2e+27)
		tmp = fma(Float64(y / fma(11.9400905721, z, 0.607771387771)), b, x);
	else
		tmp = fma(3.13060547623, y, x);
	end
	return tmp
end

Wolfram

code[x_, y_, z_, t_, a_, b_] := If[LessEqual[z, -0.55], N[(3.13060547623 * y + x), $MachinePrecision], If[LessEqual[z, 6.2e+27], N[(N[(y / N[(11.9400905721 * z + 0.607771387771), $MachinePrecision]), $MachinePrecision] * b + x), $MachinePrecision], N[(3.13060547623 * y + x), $MachinePrecision]]]

Derivation

1. Split input into 2 regimes

2. if z < -0.55000000000000004 or 6.19999999999999992e27 < z

   1. Initial program 58.7%

      \[x + \frac{y \cdot \left(\left(\left(\left(z \cdot 3.13060547623 + 11.1667541262\right) \cdot z + t\right) \cdot z + a\right) \cdot z + b\right)}{\left(\left(\left(z + 15.234687407\right) \cdot z + 31.4690115749\right) \cdot z + 11.9400905721\right) \cdot z + 0.607771387771} \]

   3. Applied rewrites 60.9%

      \[\leadsto \color{blue}{\mathsf{fma}\left(\frac{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(3.13060547623, z, 11.1667541262\right), z, t\right), z, a\right), z, b\right)}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(z - -15.234687407, z, 31.4690115749\right), z, 11.9400905721\right), z, 0.607771387771\right)}, y, x\right)} \]

   4. Taylor expanded in z around 0

      \[\leadsto \mathsf{fma}\left(\frac{\mathsf{fma}\left(\color{blue}{a}, z, b\right)}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(z - \frac{-15234687407}{1000000000}, z, \frac{314690115749}{10000000000}\right), z, \frac{119400905721}{10000000000}\right), z, \frac{607771387771}{1000000000000}\right)}, y, x\right) \]
   5. Step-by-step derivation

   6. Applied rewrites 67.3%

      \[\leadsto \mathsf{fma}\left(\frac{\mathsf{fma}\left(\color{blue}{a}, z, b\right)}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(z - -15.234687407, z, 31.4690115749\right), z, 11.9400905721\right), z, 0.607771387771\right)}, y, x\right) \]

   7. Taylor expanded in z around inf

      \[\leadsto \mathsf{fma}\left(\color{blue}{\frac{313060547623}{100000000000}}, y, x\right) \]
   8. Step-by-step derivation

   9. Applied rewrites 62.8%

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

   if -0.55000000000000004 < z < 6.19999999999999992e27

   1. Initial program 58.7%

      \[x + \frac{y \cdot \left(\left(\left(\left(z \cdot 3.13060547623 + 11.1667541262\right) \cdot z + t\right) \cdot z + a\right) \cdot z + b\right)}{\left(\left(\left(z + 15.234687407\right) \cdot z + 31.4690115749\right) \cdot z + 11.9400905721\right) \cdot z + 0.607771387771} \]

   2. Taylor expanded in z around 0

      \[\leadsto x + \frac{y \cdot \color{blue}{b}}{\left(\left(\left(z + \frac{15234687407}{1000000000}\right) \cdot z + \frac{314690115749}{10000000000}\right) \cdot z + \frac{119400905721}{10000000000}\right) \cdot z + \frac{607771387771}{1000000000000}} \]
   3. Step-by-step derivation

   4. Applied rewrites 64.2%

      \[\leadsto x + \frac{y \cdot \color{blue}{b}}{\left(\left(\left(z + 15.234687407\right) \cdot z + 31.4690115749\right) \cdot z + 11.9400905721\right) \cdot z + 0.607771387771} \]

   6. Applied rewrites 65.2%

      \[\leadsto \color{blue}{\mathsf{fma}\left(\frac{y}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(z - -15.234687407, z, 31.4690115749\right), z, 11.9400905721\right), z, 0.607771387771\right)}, b, x\right)} \]

   7. Taylor expanded in z around 0

      \[\leadsto \mathsf{fma}\left(\frac{y}{\mathsf{fma}\left(\color{blue}{\frac{119400905721}{10000000000}}, z, \frac{607771387771}{1000000000000}\right)}, b, x\right) \]
   8. Step-by-step derivation

   9. Applied rewrites 63.3%

      \[\leadsto \mathsf{fma}\left(\frac{y}{\mathsf{fma}\left(\color{blue}{11.9400905721}, z, 0.607771387771\right)}, b, x\right) \]
3. Recombined 2 regimes into one program.
4. Add Preprocessing

Alternative 17: 83.2% accurate, 3.2× speedup

Math

\[\begin{array}{l} \\ \begin{array}{l} \mathbf{if}\;z \leq -620:\\ \;\;\;\;\mathsf{fma}\left(3.13060547623, y, x\right)\\ \mathbf{elif}\;z \leq 6.2 \cdot 10^{+27}:\\ \;\;\;\;\mathsf{fma}\left(1.6453555072203998 \cdot b, y, x\right)\\ \mathbf{else}:\\ \;\;\;\;\mathsf{fma}\left(3.13060547623, y, x\right)\\ \end{array} \end{array} \]

FPCore

(FPCore (x y z t a b)
 :precision binary64
 (if (<= z -620.0)
   (fma 3.13060547623 y x)
   (if (<= z 6.2e+27)
     (fma (* 1.6453555072203998 b) y x)
     (fma 3.13060547623 y x))))

C

double code(double x, double y, double z, double t, double a, double b) {
	double tmp;
	if (z <= -620.0) {
		tmp = fma(3.13060547623, y, x);
	} else if (z <= 6.2e+27) {
		tmp = fma((1.6453555072203998 * b), y, x);
	} else {
		tmp = fma(3.13060547623, y, x);
	}
	return tmp;
}

Julia

function code(x, y, z, t, a, b)
	tmp = 0.0
	if (z <= -620.0)
		tmp = fma(3.13060547623, y, x);
	elseif (z <= 6.2e+27)
		tmp = fma(Float64(1.6453555072203998 * b), y, x);
	else
		tmp = fma(3.13060547623, y, x);
	end
	return tmp
end

Wolfram

code[x_, y_, z_, t_, a_, b_] := If[LessEqual[z, -620.0], N[(3.13060547623 * y + x), $MachinePrecision], If[LessEqual[z, 6.2e+27], N[(N[(1.6453555072203998 * b), $MachinePrecision] * y + x), $MachinePrecision], N[(3.13060547623 * y + x), $MachinePrecision]]]

Derivation

1. Split input into 2 regimes

2. if z < -620 or 6.19999999999999992e27 < z

   1. Initial program 58.7%

      \[x + \frac{y \cdot \left(\left(\left(\left(z \cdot 3.13060547623 + 11.1667541262\right) \cdot z + t\right) \cdot z + a\right) \cdot z + b\right)}{\left(\left(\left(z + 15.234687407\right) \cdot z + 31.4690115749\right) \cdot z + 11.9400905721\right) \cdot z + 0.607771387771} \]

   3. Applied rewrites 60.9%

      \[\leadsto \color{blue}{\mathsf{fma}\left(\frac{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(3.13060547623, z, 11.1667541262\right), z, t\right), z, a\right), z, b\right)}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(z - -15.234687407, z, 31.4690115749\right), z, 11.9400905721\right), z, 0.607771387771\right)}, y, x\right)} \]

   4. Taylor expanded in z around 0

      \[\leadsto \mathsf{fma}\left(\frac{\mathsf{fma}\left(\color{blue}{a}, z, b\right)}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(z - \frac{-15234687407}{1000000000}, z, \frac{314690115749}{10000000000}\right), z, \frac{119400905721}{10000000000}\right), z, \frac{607771387771}{1000000000000}\right)}, y, x\right) \]
   5. Step-by-step derivation

   6. Applied rewrites 67.3%

      \[\leadsto \mathsf{fma}\left(\frac{\mathsf{fma}\left(\color{blue}{a}, z, b\right)}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(z - -15.234687407, z, 31.4690115749\right), z, 11.9400905721\right), z, 0.607771387771\right)}, y, x\right) \]

   7. Taylor expanded in z around inf

      \[\leadsto \mathsf{fma}\left(\color{blue}{\frac{313060547623}{100000000000}}, y, x\right) \]
   8. Step-by-step derivation

   9. Applied rewrites 62.8%

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

   if -620 < z < 6.19999999999999992e27

   1. Initial program 58.7%

      \[x + \frac{y \cdot \left(\left(\left(\left(z \cdot 3.13060547623 + 11.1667541262\right) \cdot z + t\right) \cdot z + a\right) \cdot z + b\right)}{\left(\left(\left(z + 15.234687407\right) \cdot z + 31.4690115749\right) \cdot z + 11.9400905721\right) \cdot z + 0.607771387771} \]

   3. Applied rewrites 60.9%

      \[\leadsto \color{blue}{\mathsf{fma}\left(\frac{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(3.13060547623, z, 11.1667541262\right), z, t\right), z, a\right), z, b\right)}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(z - -15.234687407, z, 31.4690115749\right), z, 11.9400905721\right), z, 0.607771387771\right)}, y, x\right)} \]

   4. Taylor expanded in z around 0

      \[\leadsto \mathsf{fma}\left(\color{blue}{\frac{1000000000000}{607771387771} \cdot b}, y, x\right) \]
   5. Step-by-step derivation

      1. lower-*.f64 60.0

         \[\leadsto \mathsf{fma}\left(1.6453555072203998 \cdot \color{blue}{b}, y, x\right) \]

   6. Applied rewrites 60.0%

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

Alternative 18: 83.2% accurate, 3.2× speedup

Math

\[\begin{array}{l} \\ \begin{array}{l} \mathbf{if}\;z \leq -620:\\ \;\;\;\;\mathsf{fma}\left(3.13060547623, y, x\right)\\ \mathbf{elif}\;z \leq 6.2 \cdot 10^{+27}:\\ \;\;\;\;\mathsf{fma}\left(1.6453555072203998 \cdot y, b, x\right)\\ \mathbf{else}:\\ \;\;\;\;\mathsf{fma}\left(3.13060547623, y, x\right)\\ \end{array} \end{array} \]

FPCore

(FPCore (x y z t a b)
 :precision binary64
 (if (<= z -620.0)
   (fma 3.13060547623 y x)
   (if (<= z 6.2e+27)
     (fma (* 1.6453555072203998 y) b x)
     (fma 3.13060547623 y x))))

C

double code(double x, double y, double z, double t, double a, double b) {
	double tmp;
	if (z <= -620.0) {
		tmp = fma(3.13060547623, y, x);
	} else if (z <= 6.2e+27) {
		tmp = fma((1.6453555072203998 * y), b, x);
	} else {
		tmp = fma(3.13060547623, y, x);
	}
	return tmp;
}

Julia

function code(x, y, z, t, a, b)
	tmp = 0.0
	if (z <= -620.0)
		tmp = fma(3.13060547623, y, x);
	elseif (z <= 6.2e+27)
		tmp = fma(Float64(1.6453555072203998 * y), b, x);
	else
		tmp = fma(3.13060547623, y, x);
	end
	return tmp
end

Wolfram

code[x_, y_, z_, t_, a_, b_] := If[LessEqual[z, -620.0], N[(3.13060547623 * y + x), $MachinePrecision], If[LessEqual[z, 6.2e+27], N[(N[(1.6453555072203998 * y), $MachinePrecision] * b + x), $MachinePrecision], N[(3.13060547623 * y + x), $MachinePrecision]]]

Derivation

1. Split input into 2 regimes

2. if z < -620 or 6.19999999999999992e27 < z

   1. Initial program 58.7%

      \[x + \frac{y \cdot \left(\left(\left(\left(z \cdot 3.13060547623 + 11.1667541262\right) \cdot z + t\right) \cdot z + a\right) \cdot z + b\right)}{\left(\left(\left(z + 15.234687407\right) \cdot z + 31.4690115749\right) \cdot z + 11.9400905721\right) \cdot z + 0.607771387771} \]

   3. Applied rewrites 60.9%

      \[\leadsto \color{blue}{\mathsf{fma}\left(\frac{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(3.13060547623, z, 11.1667541262\right), z, t\right), z, a\right), z, b\right)}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(z - -15.234687407, z, 31.4690115749\right), z, 11.9400905721\right), z, 0.607771387771\right)}, y, x\right)} \]

   4. Taylor expanded in z around 0

      \[\leadsto \mathsf{fma}\left(\frac{\mathsf{fma}\left(\color{blue}{a}, z, b\right)}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(z - \frac{-15234687407}{1000000000}, z, \frac{314690115749}{10000000000}\right), z, \frac{119400905721}{10000000000}\right), z, \frac{607771387771}{1000000000000}\right)}, y, x\right) \]
   5. Step-by-step derivation

   6. Applied rewrites 67.3%

      \[\leadsto \mathsf{fma}\left(\frac{\mathsf{fma}\left(\color{blue}{a}, z, b\right)}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(z - -15.234687407, z, 31.4690115749\right), z, 11.9400905721\right), z, 0.607771387771\right)}, y, x\right) \]

   7. Taylor expanded in z around inf

      \[\leadsto \mathsf{fma}\left(\color{blue}{\frac{313060547623}{100000000000}}, y, x\right) \]
   8. Step-by-step derivation

   9. Applied rewrites 62.8%

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

   if -620 < z < 6.19999999999999992e27

   1. Initial program 58.7%

      \[x + \frac{y \cdot \left(\left(\left(\left(z \cdot 3.13060547623 + 11.1667541262\right) \cdot z + t\right) \cdot z + a\right) \cdot z + b\right)}{\left(\left(\left(z + 15.234687407\right) \cdot z + 31.4690115749\right) \cdot z + 11.9400905721\right) \cdot z + 0.607771387771} \]

   2. Taylor expanded in z around 0

      \[\leadsto x + \frac{y \cdot \color{blue}{b}}{\left(\left(\left(z + \frac{15234687407}{1000000000}\right) \cdot z + \frac{314690115749}{10000000000}\right) \cdot z + \frac{119400905721}{10000000000}\right) \cdot z + \frac{607771387771}{1000000000000}} \]
   3. Step-by-step derivation

   4. Applied rewrites 64.2%

      \[\leadsto x + \frac{y \cdot \color{blue}{b}}{\left(\left(\left(z + 15.234687407\right) \cdot z + 31.4690115749\right) \cdot z + 11.9400905721\right) \cdot z + 0.607771387771} \]

   6. Applied rewrites 65.2%

      \[\leadsto \color{blue}{\mathsf{fma}\left(\frac{y}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(z - -15.234687407, z, 31.4690115749\right), z, 11.9400905721\right), z, 0.607771387771\right)}, b, x\right)} \]

   7. Taylor expanded in z around 0

      \[\leadsto \mathsf{fma}\left(\color{blue}{\frac{1000000000000}{607771387771} \cdot y}, b, x\right) \]
   8. Step-by-step derivation

      1. lower-*.f64 60.0

         \[\leadsto \mathsf{fma}\left(1.6453555072203998 \cdot \color{blue}{y}, b, x\right) \]

   9. Applied rewrites 60.0%

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

Alternative 19: 62.8% accurate, 8.8× speedup

Math

\[\begin{array}{l} \\ \mathsf{fma}\left(3.13060547623, y, x\right) \end{array} \]

FPCore

(FPCore (x y z t a b) :precision binary64 (fma 3.13060547623 y x))

C

double code(double x, double y, double z, double t, double a, double b) {
	return fma(3.13060547623, y, x);
}

Julia

function code(x, y, z, t, a, b)
	return fma(3.13060547623, y, x)
end

Wolfram

code[x_, y_, z_, t_, a_, b_] := N[(3.13060547623 * y + x), $MachinePrecision]

Derivation

1. Initial program 58.7%

   \[x + \frac{y \cdot \left(\left(\left(\left(z \cdot 3.13060547623 + 11.1667541262\right) \cdot z + t\right) \cdot z + a\right) \cdot z + b\right)}{\left(\left(\left(z + 15.234687407\right) \cdot z + 31.4690115749\right) \cdot z + 11.9400905721\right) \cdot z + 0.607771387771} \]

3. Applied rewrites 60.9%

   \[\leadsto \color{blue}{\mathsf{fma}\left(\frac{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(3.13060547623, z, 11.1667541262\right), z, t\right), z, a\right), z, b\right)}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(z - -15.234687407, z, 31.4690115749\right), z, 11.9400905721\right), z, 0.607771387771\right)}, y, x\right)} \]

4. Taylor expanded in z around 0

   \[\leadsto \mathsf{fma}\left(\frac{\mathsf{fma}\left(\color{blue}{a}, z, b\right)}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(z - \frac{-15234687407}{1000000000}, z, \frac{314690115749}{10000000000}\right), z, \frac{119400905721}{10000000000}\right), z, \frac{607771387771}{1000000000000}\right)}, y, x\right) \]
5. Step-by-step derivation

6. Applied rewrites 67.3%

   \[\leadsto \mathsf{fma}\left(\frac{\mathsf{fma}\left(\color{blue}{a}, z, b\right)}{\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(z - -15.234687407, z, 31.4690115749\right), z, 11.9400905721\right), z, 0.607771387771\right)}, y, x\right) \]

7. Taylor expanded in z around inf

   \[\leadsto \mathsf{fma}\left(\color{blue}{\frac{313060547623}{100000000000}}, y, x\right) \]
8. Step-by-step derivation

9. Applied rewrites 62.8%

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

Reproduce

herbie shell --seed 2025154 
(FPCore (x y z t a b)
  :name "Numeric.SpecFunctions:logGamma from math-functions-0.1.5.2, D"
  :precision binary64
  (+ x (/ (* y (+ (* (+ (* (+ (* (+ (* z 3.13060547623) 11.1667541262) z) t) z) a) z) b)) (+ (* (+ (* (+ (* (+ z 15.234687407) z) 31.4690115749) z) 11.9400905721) z) 0.607771387771))))