
(FPCore (x y z) :precision binary64 (+ (+ x (sin y)) (* z (cos y))))
double code(double x, double y, double z) {
return (x + sin(y)) + (z * cos(y));
}
real(8) function code(x, y, z)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
code = (x + sin(y)) + (z * cos(y))
end function
public static double code(double x, double y, double z) {
return (x + Math.sin(y)) + (z * Math.cos(y));
}
def code(x, y, z): return (x + math.sin(y)) + (z * math.cos(y))
function code(x, y, z) return Float64(Float64(x + sin(y)) + Float64(z * cos(y))) end
function tmp = code(x, y, z) tmp = (x + sin(y)) + (z * cos(y)); end
code[x_, y_, z_] := N[(N[(x + N[Sin[y], $MachinePrecision]), $MachinePrecision] + N[(z * N[Cos[y], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\left(x + \sin y\right) + z \cdot \cos y
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 12 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (x y z) :precision binary64 (+ (+ x (sin y)) (* z (cos y))))
double code(double x, double y, double z) {
return (x + sin(y)) + (z * cos(y));
}
real(8) function code(x, y, z)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
code = (x + sin(y)) + (z * cos(y))
end function
public static double code(double x, double y, double z) {
return (x + Math.sin(y)) + (z * Math.cos(y));
}
def code(x, y, z): return (x + math.sin(y)) + (z * math.cos(y))
function code(x, y, z) return Float64(Float64(x + sin(y)) + Float64(z * cos(y))) end
function tmp = code(x, y, z) tmp = (x + sin(y)) + (z * cos(y)); end
code[x_, y_, z_] := N[(N[(x + N[Sin[y], $MachinePrecision]), $MachinePrecision] + N[(z * N[Cos[y], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\left(x + \sin y\right) + z \cdot \cos y
\end{array}
(FPCore (x y z) :precision binary64 (fma (cos y) z (+ x (sin y))))
double code(double x, double y, double z) {
return fma(cos(y), z, (x + sin(y)));
}
function code(x, y, z) return fma(cos(y), z, Float64(x + sin(y))) end
code[x_, y_, z_] := N[(N[Cos[y], $MachinePrecision] * z + N[(x + N[Sin[y], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\mathsf{fma}\left(\cos y, z, x + \sin y\right)
\end{array}
Initial program 99.9%
lift-+.f64N/A
+-commutativeN/A
lift-*.f64N/A
*-commutativeN/A
lower-fma.f6499.9
lift-+.f64N/A
+-commutativeN/A
lower-+.f6499.9
Applied rewrites99.9%
Final simplification99.9%
(FPCore (x y z) :precision binary64 (if (<= x -1.3e-67) (fma 1.0 z (+ x (sin y))) (if (<= x 3e+40) (fma (cos y) z (sin y)) (fma (/ z x) x x))))
double code(double x, double y, double z) {
double tmp;
if (x <= -1.3e-67) {
tmp = fma(1.0, z, (x + sin(y)));
} else if (x <= 3e+40) {
tmp = fma(cos(y), z, sin(y));
} else {
tmp = fma((z / x), x, x);
}
return tmp;
}
function code(x, y, z) tmp = 0.0 if (x <= -1.3e-67) tmp = fma(1.0, z, Float64(x + sin(y))); elseif (x <= 3e+40) tmp = fma(cos(y), z, sin(y)); else tmp = fma(Float64(z / x), x, x); end return tmp end
code[x_, y_, z_] := If[LessEqual[x, -1.3e-67], N[(1.0 * z + N[(x + N[Sin[y], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[x, 3e+40], N[(N[Cos[y], $MachinePrecision] * z + N[Sin[y], $MachinePrecision]), $MachinePrecision], N[(N[(z / x), $MachinePrecision] * x + x), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq -1.3 \cdot 10^{-67}:\\
\;\;\;\;\mathsf{fma}\left(1, z, x + \sin y\right)\\
\mathbf{elif}\;x \leq 3 \cdot 10^{+40}:\\
\;\;\;\;\mathsf{fma}\left(\cos y, z, \sin y\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(\frac{z}{x}, x, x\right)\\
\end{array}
\end{array}
if x < -1.2999999999999999e-67Initial program 99.9%
lift-+.f64N/A
+-commutativeN/A
lift-*.f64N/A
*-commutativeN/A
lower-fma.f6499.9
lift-+.f64N/A
+-commutativeN/A
lower-+.f6499.9
Applied rewrites99.9%
Taylor expanded in y around 0
Applied rewrites87.7%
if -1.2999999999999999e-67 < x < 3.0000000000000002e40Initial program 99.9%
Taylor expanded in x around 0
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower-cos.f64N/A
lower-sin.f6494.0
Applied rewrites94.0%
if 3.0000000000000002e40 < x Initial program 100.0%
Taylor expanded in y around 0
+-commutativeN/A
lower-+.f6491.7
Applied rewrites91.7%
Taylor expanded in x around inf
Applied rewrites91.7%
Final simplification91.6%
(FPCore (x y z) :precision binary64 (if (<= z -3.8e+121) (/ 1.0 (/ 1.0 (* z (cos y)))) (if (<= z 4.4e+88) (fma 1.0 z (+ x (sin y))) (fma (cos y) z (+ x y)))))
double code(double x, double y, double z) {
double tmp;
if (z <= -3.8e+121) {
tmp = 1.0 / (1.0 / (z * cos(y)));
} else if (z <= 4.4e+88) {
tmp = fma(1.0, z, (x + sin(y)));
} else {
tmp = fma(cos(y), z, (x + y));
}
return tmp;
}
function code(x, y, z) tmp = 0.0 if (z <= -3.8e+121) tmp = Float64(1.0 / Float64(1.0 / Float64(z * cos(y)))); elseif (z <= 4.4e+88) tmp = fma(1.0, z, Float64(x + sin(y))); else tmp = fma(cos(y), z, Float64(x + y)); end return tmp end
code[x_, y_, z_] := If[LessEqual[z, -3.8e+121], N[(1.0 / N[(1.0 / N[(z * N[Cos[y], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[z, 4.4e+88], N[(1.0 * z + N[(x + N[Sin[y], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[Cos[y], $MachinePrecision] * z + N[(x + y), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;z \leq -3.8 \cdot 10^{+121}:\\
\;\;\;\;\frac{1}{\frac{1}{z \cdot \cos y}}\\
\mathbf{elif}\;z \leq 4.4 \cdot 10^{+88}:\\
\;\;\;\;\mathsf{fma}\left(1, z, x + \sin y\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(\cos y, z, x + y\right)\\
\end{array}
\end{array}
if z < -3.8e121Initial program 99.9%
lift-+.f64N/A
flip3-+N/A
clear-numN/A
lower-/.f64N/A
clear-numN/A
flip3-+N/A
lift-+.f64N/A
inv-powN/A
lower-pow.f6499.9
Applied rewrites99.9%
Taylor expanded in z around inf
associate-/l/N/A
lower-/.f64N/A
lower-/.f64N/A
lower-cos.f6485.4
Applied rewrites85.4%
Applied rewrites85.5%
if -3.8e121 < z < 4.40000000000000017e88Initial program 100.0%
lift-+.f64N/A
+-commutativeN/A
lift-*.f64N/A
*-commutativeN/A
lower-fma.f64100.0
lift-+.f64N/A
+-commutativeN/A
lower-+.f64100.0
Applied rewrites100.0%
Taylor expanded in y around 0
Applied rewrites93.6%
if 4.40000000000000017e88 < z Initial program 99.7%
lift-+.f64N/A
+-commutativeN/A
lift-*.f64N/A
*-commutativeN/A
lower-fma.f6499.8
lift-+.f64N/A
+-commutativeN/A
lower-+.f6499.8
Applied rewrites99.8%
Taylor expanded in y around 0
+-commutativeN/A
lower-+.f6487.1
Applied rewrites87.1%
Final simplification90.9%
(FPCore (x y z)
:precision binary64
(let* ((t_0 (* z (cos y))))
(if (<= z -3.8e+121)
t_0
(if (<= z -3e-20)
(+ x z)
(if (<= z 8e-12) (+ x (sin y)) (if (<= z 1.3e+107) (+ x z) t_0))))))
double code(double x, double y, double z) {
double t_0 = z * cos(y);
double tmp;
if (z <= -3.8e+121) {
tmp = t_0;
} else if (z <= -3e-20) {
tmp = x + z;
} else if (z <= 8e-12) {
tmp = x + sin(y);
} else if (z <= 1.3e+107) {
tmp = x + z;
} else {
tmp = t_0;
}
return tmp;
}
real(8) function code(x, y, z)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8) :: t_0
real(8) :: tmp
t_0 = z * cos(y)
if (z <= (-3.8d+121)) then
tmp = t_0
else if (z <= (-3d-20)) then
tmp = x + z
else if (z <= 8d-12) then
tmp = x + sin(y)
else if (z <= 1.3d+107) then
tmp = x + z
else
tmp = t_0
end if
code = tmp
end function
public static double code(double x, double y, double z) {
double t_0 = z * Math.cos(y);
double tmp;
if (z <= -3.8e+121) {
tmp = t_0;
} else if (z <= -3e-20) {
tmp = x + z;
} else if (z <= 8e-12) {
tmp = x + Math.sin(y);
} else if (z <= 1.3e+107) {
tmp = x + z;
} else {
tmp = t_0;
}
return tmp;
}
def code(x, y, z): t_0 = z * math.cos(y) tmp = 0 if z <= -3.8e+121: tmp = t_0 elif z <= -3e-20: tmp = x + z elif z <= 8e-12: tmp = x + math.sin(y) elif z <= 1.3e+107: tmp = x + z else: tmp = t_0 return tmp
function code(x, y, z) t_0 = Float64(z * cos(y)) tmp = 0.0 if (z <= -3.8e+121) tmp = t_0; elseif (z <= -3e-20) tmp = Float64(x + z); elseif (z <= 8e-12) tmp = Float64(x + sin(y)); elseif (z <= 1.3e+107) tmp = Float64(x + z); else tmp = t_0; end return tmp end
function tmp_2 = code(x, y, z) t_0 = z * cos(y); tmp = 0.0; if (z <= -3.8e+121) tmp = t_0; elseif (z <= -3e-20) tmp = x + z; elseif (z <= 8e-12) tmp = x + sin(y); elseif (z <= 1.3e+107) tmp = x + z; else tmp = t_0; end tmp_2 = tmp; end
code[x_, y_, z_] := Block[{t$95$0 = N[(z * N[Cos[y], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[z, -3.8e+121], t$95$0, If[LessEqual[z, -3e-20], N[(x + z), $MachinePrecision], If[LessEqual[z, 8e-12], N[(x + N[Sin[y], $MachinePrecision]), $MachinePrecision], If[LessEqual[z, 1.3e+107], N[(x + z), $MachinePrecision], t$95$0]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := z \cdot \cos y\\
\mathbf{if}\;z \leq -3.8 \cdot 10^{+121}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;z \leq -3 \cdot 10^{-20}:\\
\;\;\;\;x + z\\
\mathbf{elif}\;z \leq 8 \cdot 10^{-12}:\\
\;\;\;\;x + \sin y\\
\mathbf{elif}\;z \leq 1.3 \cdot 10^{+107}:\\
\;\;\;\;x + z\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
if z < -3.8e121 or 1.3000000000000001e107 < z Initial program 99.8%
Taylor expanded in z around inf
*-commutativeN/A
lower-*.f64N/A
lower-cos.f6485.7
Applied rewrites85.7%
if -3.8e121 < z < -3.00000000000000029e-20 or 7.99999999999999984e-12 < z < 1.3000000000000001e107Initial program 99.9%
Taylor expanded in y around 0
+-commutativeN/A
lower-+.f6482.6
Applied rewrites82.6%
if -3.00000000000000029e-20 < z < 7.99999999999999984e-12Initial program 100.0%
Taylor expanded in z around 0
+-commutativeN/A
lower-+.f64N/A
lower-sin.f6495.5
Applied rewrites95.5%
Final simplification88.9%
(FPCore (x y z) :precision binary64 (if (<= z -3.8e+121) (* z (cos y)) (if (<= z 4.4e+88) (fma 1.0 z (+ x (sin y))) (fma (cos y) z (+ x y)))))
double code(double x, double y, double z) {
double tmp;
if (z <= -3.8e+121) {
tmp = z * cos(y);
} else if (z <= 4.4e+88) {
tmp = fma(1.0, z, (x + sin(y)));
} else {
tmp = fma(cos(y), z, (x + y));
}
return tmp;
}
function code(x, y, z) tmp = 0.0 if (z <= -3.8e+121) tmp = Float64(z * cos(y)); elseif (z <= 4.4e+88) tmp = fma(1.0, z, Float64(x + sin(y))); else tmp = fma(cos(y), z, Float64(x + y)); end return tmp end
code[x_, y_, z_] := If[LessEqual[z, -3.8e+121], N[(z * N[Cos[y], $MachinePrecision]), $MachinePrecision], If[LessEqual[z, 4.4e+88], N[(1.0 * z + N[(x + N[Sin[y], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[Cos[y], $MachinePrecision] * z + N[(x + y), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;z \leq -3.8 \cdot 10^{+121}:\\
\;\;\;\;z \cdot \cos y\\
\mathbf{elif}\;z \leq 4.4 \cdot 10^{+88}:\\
\;\;\;\;\mathsf{fma}\left(1, z, x + \sin y\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(\cos y, z, x + y\right)\\
\end{array}
\end{array}
if z < -3.8e121Initial program 99.9%
Taylor expanded in z around inf
*-commutativeN/A
lower-*.f64N/A
lower-cos.f6485.5
Applied rewrites85.5%
if -3.8e121 < z < 4.40000000000000017e88Initial program 100.0%
lift-+.f64N/A
+-commutativeN/A
lift-*.f64N/A
*-commutativeN/A
lower-fma.f64100.0
lift-+.f64N/A
+-commutativeN/A
lower-+.f64100.0
Applied rewrites100.0%
Taylor expanded in y around 0
Applied rewrites93.6%
if 4.40000000000000017e88 < z Initial program 99.7%
lift-+.f64N/A
+-commutativeN/A
lift-*.f64N/A
*-commutativeN/A
lower-fma.f6499.8
lift-+.f64N/A
+-commutativeN/A
lower-+.f6499.8
Applied rewrites99.8%
Taylor expanded in y around 0
+-commutativeN/A
lower-+.f6487.1
Applied rewrites87.1%
Final simplification90.9%
(FPCore (x y z)
:precision binary64
(let* ((t_0 (* z (cos y))))
(if (<= z -3.8e+121)
t_0
(if (<= z 1.3e+107) (fma 1.0 z (+ x (sin y))) t_0))))
double code(double x, double y, double z) {
double t_0 = z * cos(y);
double tmp;
if (z <= -3.8e+121) {
tmp = t_0;
} else if (z <= 1.3e+107) {
tmp = fma(1.0, z, (x + sin(y)));
} else {
tmp = t_0;
}
return tmp;
}
function code(x, y, z) t_0 = Float64(z * cos(y)) tmp = 0.0 if (z <= -3.8e+121) tmp = t_0; elseif (z <= 1.3e+107) tmp = fma(1.0, z, Float64(x + sin(y))); else tmp = t_0; end return tmp end
code[x_, y_, z_] := Block[{t$95$0 = N[(z * N[Cos[y], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[z, -3.8e+121], t$95$0, If[LessEqual[z, 1.3e+107], N[(1.0 * z + N[(x + N[Sin[y], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$0]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := z \cdot \cos y\\
\mathbf{if}\;z \leq -3.8 \cdot 10^{+121}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;z \leq 1.3 \cdot 10^{+107}:\\
\;\;\;\;\mathsf{fma}\left(1, z, x + \sin y\right)\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
if z < -3.8e121 or 1.3000000000000001e107 < z Initial program 99.8%
Taylor expanded in z around inf
*-commutativeN/A
lower-*.f64N/A
lower-cos.f6485.7
Applied rewrites85.7%
if -3.8e121 < z < 1.3000000000000001e107Initial program 100.0%
lift-+.f64N/A
+-commutativeN/A
lift-*.f64N/A
*-commutativeN/A
lower-fma.f64100.0
lift-+.f64N/A
+-commutativeN/A
lower-+.f64100.0
Applied rewrites100.0%
Taylor expanded in y around 0
Applied rewrites93.7%
Final simplification90.7%
(FPCore (x y z)
:precision binary64
(let* ((t_0 (+ x (sin y))))
(if (<= y -950.0)
t_0
(if (<= y 0.27) (fma (fma (* -0.5 y) z 1.0) y (+ x z)) t_0))))
double code(double x, double y, double z) {
double t_0 = x + sin(y);
double tmp;
if (y <= -950.0) {
tmp = t_0;
} else if (y <= 0.27) {
tmp = fma(fma((-0.5 * y), z, 1.0), y, (x + z));
} else {
tmp = t_0;
}
return tmp;
}
function code(x, y, z) t_0 = Float64(x + sin(y)) tmp = 0.0 if (y <= -950.0) tmp = t_0; elseif (y <= 0.27) tmp = fma(fma(Float64(-0.5 * y), z, 1.0), y, Float64(x + z)); else tmp = t_0; end return tmp end
code[x_, y_, z_] := Block[{t$95$0 = N[(x + N[Sin[y], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -950.0], t$95$0, If[LessEqual[y, 0.27], N[(N[(N[(-0.5 * y), $MachinePrecision] * z + 1.0), $MachinePrecision] * y + N[(x + z), $MachinePrecision]), $MachinePrecision], t$95$0]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := x + \sin y\\
\mathbf{if}\;y \leq -950:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;y \leq 0.27:\\
\;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(-0.5 \cdot y, z, 1\right), y, x + z\right)\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
if y < -950 or 0.27000000000000002 < y Initial program 99.8%
Taylor expanded in z around 0
+-commutativeN/A
lower-+.f64N/A
lower-sin.f6461.7
Applied rewrites61.7%
if -950 < y < 0.27000000000000002Initial program 100.0%
Taylor expanded in y around 0
associate-+r+N/A
+-commutativeN/A
*-commutativeN/A
*-commutativeN/A
associate-*r*N/A
lower-fma.f64N/A
+-commutativeN/A
associate-*r*N/A
*-commutativeN/A
associate-*r*N/A
lower-fma.f64N/A
lower-*.f64N/A
+-commutativeN/A
lower-+.f6499.2
Applied rewrites99.2%
Final simplification80.5%
(FPCore (x y z)
:precision binary64
(if (<= y -9.6e+23)
(+ x z)
(if (<= y 6000.0)
(fma (fma (fma -0.16666666666666666 y (* -0.5 z)) y 1.0) y (+ x z))
(+ x z))))
double code(double x, double y, double z) {
double tmp;
if (y <= -9.6e+23) {
tmp = x + z;
} else if (y <= 6000.0) {
tmp = fma(fma(fma(-0.16666666666666666, y, (-0.5 * z)), y, 1.0), y, (x + z));
} else {
tmp = x + z;
}
return tmp;
}
function code(x, y, z) tmp = 0.0 if (y <= -9.6e+23) tmp = Float64(x + z); elseif (y <= 6000.0) tmp = fma(fma(fma(-0.16666666666666666, y, Float64(-0.5 * z)), y, 1.0), y, Float64(x + z)); else tmp = Float64(x + z); end return tmp end
code[x_, y_, z_] := If[LessEqual[y, -9.6e+23], N[(x + z), $MachinePrecision], If[LessEqual[y, 6000.0], N[(N[(N[(-0.16666666666666666 * y + N[(-0.5 * z), $MachinePrecision]), $MachinePrecision] * y + 1.0), $MachinePrecision] * y + N[(x + z), $MachinePrecision]), $MachinePrecision], N[(x + z), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -9.6 \cdot 10^{+23}:\\
\;\;\;\;x + z\\
\mathbf{elif}\;y \leq 6000:\\
\;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(\mathsf{fma}\left(-0.16666666666666666, y, -0.5 \cdot z\right), y, 1\right), y, x + z\right)\\
\mathbf{else}:\\
\;\;\;\;x + z\\
\end{array}
\end{array}
if y < -9.6e23 or 6e3 < y Initial program 99.8%
Taylor expanded in y around 0
+-commutativeN/A
lower-+.f6441.4
Applied rewrites41.4%
if -9.6e23 < y < 6e3Initial program 100.0%
Taylor expanded in y around 0
associate-+r+N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
lower-*.f64N/A
+-commutativeN/A
lower-+.f6496.5
Applied rewrites96.5%
Final simplification70.3%
(FPCore (x y z) :precision binary64 (if (<= y -1.2e+37) (+ x z) (if (<= y 9000000000000.0) (fma (fma (* -0.5 y) z 1.0) y (+ x z)) (+ x z))))
double code(double x, double y, double z) {
double tmp;
if (y <= -1.2e+37) {
tmp = x + z;
} else if (y <= 9000000000000.0) {
tmp = fma(fma((-0.5 * y), z, 1.0), y, (x + z));
} else {
tmp = x + z;
}
return tmp;
}
function code(x, y, z) tmp = 0.0 if (y <= -1.2e+37) tmp = Float64(x + z); elseif (y <= 9000000000000.0) tmp = fma(fma(Float64(-0.5 * y), z, 1.0), y, Float64(x + z)); else tmp = Float64(x + z); end return tmp end
code[x_, y_, z_] := If[LessEqual[y, -1.2e+37], N[(x + z), $MachinePrecision], If[LessEqual[y, 9000000000000.0], N[(N[(N[(-0.5 * y), $MachinePrecision] * z + 1.0), $MachinePrecision] * y + N[(x + z), $MachinePrecision]), $MachinePrecision], N[(x + z), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -1.2 \cdot 10^{+37}:\\
\;\;\;\;x + z\\
\mathbf{elif}\;y \leq 9000000000000:\\
\;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(-0.5 \cdot y, z, 1\right), y, x + z\right)\\
\mathbf{else}:\\
\;\;\;\;x + z\\
\end{array}
\end{array}
if y < -1.2e37 or 9e12 < y Initial program 99.8%
Taylor expanded in y around 0
+-commutativeN/A
lower-+.f6441.2
Applied rewrites41.2%
if -1.2e37 < y < 9e12Initial program 100.0%
Taylor expanded in y around 0
associate-+r+N/A
+-commutativeN/A
*-commutativeN/A
*-commutativeN/A
associate-*r*N/A
lower-fma.f64N/A
+-commutativeN/A
associate-*r*N/A
*-commutativeN/A
associate-*r*N/A
lower-fma.f64N/A
lower-*.f64N/A
+-commutativeN/A
lower-+.f6495.8
Applied rewrites95.8%
Final simplification70.2%
(FPCore (x y z)
:precision binary64
(if (<= y -9.6e+23)
(+ x z)
(if (<= y 6200.0)
(fma (fma (* -0.16666666666666666 y) y 1.0) y (+ x z))
(+ x z))))
double code(double x, double y, double z) {
double tmp;
if (y <= -9.6e+23) {
tmp = x + z;
} else if (y <= 6200.0) {
tmp = fma(fma((-0.16666666666666666 * y), y, 1.0), y, (x + z));
} else {
tmp = x + z;
}
return tmp;
}
function code(x, y, z) tmp = 0.0 if (y <= -9.6e+23) tmp = Float64(x + z); elseif (y <= 6200.0) tmp = fma(fma(Float64(-0.16666666666666666 * y), y, 1.0), y, Float64(x + z)); else tmp = Float64(x + z); end return tmp end
code[x_, y_, z_] := If[LessEqual[y, -9.6e+23], N[(x + z), $MachinePrecision], If[LessEqual[y, 6200.0], N[(N[(N[(-0.16666666666666666 * y), $MachinePrecision] * y + 1.0), $MachinePrecision] * y + N[(x + z), $MachinePrecision]), $MachinePrecision], N[(x + z), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -9.6 \cdot 10^{+23}:\\
\;\;\;\;x + z\\
\mathbf{elif}\;y \leq 6200:\\
\;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(-0.16666666666666666 \cdot y, y, 1\right), y, x + z\right)\\
\mathbf{else}:\\
\;\;\;\;x + z\\
\end{array}
\end{array}
if y < -9.6e23 or 6200 < y Initial program 99.8%
Taylor expanded in y around 0
+-commutativeN/A
lower-+.f6441.4
Applied rewrites41.4%
if -9.6e23 < y < 6200Initial program 100.0%
Taylor expanded in y around 0
associate-+r+N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
+-commutativeN/A
lower-fma.f64N/A
lower-*.f64N/A
+-commutativeN/A
lower-+.f6496.5
Applied rewrites96.5%
Taylor expanded in z around 0
Applied rewrites96.1%
Final simplification70.0%
(FPCore (x y z) :precision binary64 (if (<= y -2200.0) (+ x z) (if (<= y 9.5e+72) (+ (+ x y) z) (+ x z))))
double code(double x, double y, double z) {
double tmp;
if (y <= -2200.0) {
tmp = x + z;
} else if (y <= 9.5e+72) {
tmp = (x + y) + z;
} else {
tmp = x + z;
}
return tmp;
}
real(8) function code(x, y, z)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8) :: tmp
if (y <= (-2200.0d0)) then
tmp = x + z
else if (y <= 9.5d+72) then
tmp = (x + y) + z
else
tmp = x + z
end if
code = tmp
end function
public static double code(double x, double y, double z) {
double tmp;
if (y <= -2200.0) {
tmp = x + z;
} else if (y <= 9.5e+72) {
tmp = (x + y) + z;
} else {
tmp = x + z;
}
return tmp;
}
def code(x, y, z): tmp = 0 if y <= -2200.0: tmp = x + z elif y <= 9.5e+72: tmp = (x + y) + z else: tmp = x + z return tmp
function code(x, y, z) tmp = 0.0 if (y <= -2200.0) tmp = Float64(x + z); elseif (y <= 9.5e+72) tmp = Float64(Float64(x + y) + z); else tmp = Float64(x + z); end return tmp end
function tmp_2 = code(x, y, z) tmp = 0.0; if (y <= -2200.0) tmp = x + z; elseif (y <= 9.5e+72) tmp = (x + y) + z; else tmp = x + z; end tmp_2 = tmp; end
code[x_, y_, z_] := If[LessEqual[y, -2200.0], N[(x + z), $MachinePrecision], If[LessEqual[y, 9.5e+72], N[(N[(x + y), $MachinePrecision] + z), $MachinePrecision], N[(x + z), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -2200:\\
\;\;\;\;x + z\\
\mathbf{elif}\;y \leq 9.5 \cdot 10^{+72}:\\
\;\;\;\;\left(x + y\right) + z\\
\mathbf{else}:\\
\;\;\;\;x + z\\
\end{array}
\end{array}
if y < -2200 or 9.50000000000000054e72 < y Initial program 99.8%
Taylor expanded in y around 0
+-commutativeN/A
lower-+.f6438.2
Applied rewrites38.2%
if -2200 < y < 9.50000000000000054e72Initial program 100.0%
Taylor expanded in y around 0
associate-+r+N/A
lower-+.f64N/A
+-commutativeN/A
lower-+.f6495.1
Applied rewrites95.1%
Final simplification70.0%
(FPCore (x y z) :precision binary64 (+ x z))
double code(double x, double y, double z) {
return x + z;
}
real(8) function code(x, y, z)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
code = x + z
end function
public static double code(double x, double y, double z) {
return x + z;
}
def code(x, y, z): return x + z
function code(x, y, z) return Float64(x + z) end
function tmp = code(x, y, z) tmp = x + z; end
code[x_, y_, z_] := N[(x + z), $MachinePrecision]
\begin{array}{l}
\\
x + z
\end{array}
Initial program 99.9%
Taylor expanded in y around 0
+-commutativeN/A
lower-+.f6466.0
Applied rewrites66.0%
Final simplification66.0%
herbie shell --seed 2024276
(FPCore (x y z)
:name "Graphics.Rasterific.Svg.PathConverter:segmentToBezier from rasterific-svg-0.2.3.1, C"
:precision binary64
(+ (+ x (sin y)) (* z (cos y))))