
(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 8 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-sin.f64N/A
lift-+.f64N/A
lift-cos.f64N/A
lift-*.f64N/A
+-commutativeN/A
lift-*.f64N/A
*-commutativeN/A
lower-fma.f6499.9
Applied rewrites99.9%
(FPCore (x y z)
:precision binary64
(let* ((t_0 (+ (* (cos y) z) (+ x (sin y)))))
(if (<= t_0 -2e+28)
(+ z x)
(if (<= t_0 -0.005)
(sin y)
(if (<= t_0 5e-12)
(+ y (+ z x))
(if (<= t_0 2e+18) (sin y) (+ z x)))))))
double code(double x, double y, double z) {
double t_0 = (cos(y) * z) + (x + sin(y));
double tmp;
if (t_0 <= -2e+28) {
tmp = z + x;
} else if (t_0 <= -0.005) {
tmp = sin(y);
} else if (t_0 <= 5e-12) {
tmp = y + (z + x);
} else if (t_0 <= 2e+18) {
tmp = sin(y);
} else {
tmp = z + x;
}
return tmp;
}
real(8) function code(x, y, z)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8) :: t_0
real(8) :: tmp
t_0 = (cos(y) * z) + (x + sin(y))
if (t_0 <= (-2d+28)) then
tmp = z + x
else if (t_0 <= (-0.005d0)) then
tmp = sin(y)
else if (t_0 <= 5d-12) then
tmp = y + (z + x)
else if (t_0 <= 2d+18) then
tmp = sin(y)
else
tmp = z + x
end if
code = tmp
end function
public static double code(double x, double y, double z) {
double t_0 = (Math.cos(y) * z) + (x + Math.sin(y));
double tmp;
if (t_0 <= -2e+28) {
tmp = z + x;
} else if (t_0 <= -0.005) {
tmp = Math.sin(y);
} else if (t_0 <= 5e-12) {
tmp = y + (z + x);
} else if (t_0 <= 2e+18) {
tmp = Math.sin(y);
} else {
tmp = z + x;
}
return tmp;
}
def code(x, y, z): t_0 = (math.cos(y) * z) + (x + math.sin(y)) tmp = 0 if t_0 <= -2e+28: tmp = z + x elif t_0 <= -0.005: tmp = math.sin(y) elif t_0 <= 5e-12: tmp = y + (z + x) elif t_0 <= 2e+18: tmp = math.sin(y) else: tmp = z + x return tmp
function code(x, y, z) t_0 = Float64(Float64(cos(y) * z) + Float64(x + sin(y))) tmp = 0.0 if (t_0 <= -2e+28) tmp = Float64(z + x); elseif (t_0 <= -0.005) tmp = sin(y); elseif (t_0 <= 5e-12) tmp = Float64(y + Float64(z + x)); elseif (t_0 <= 2e+18) tmp = sin(y); else tmp = Float64(z + x); end return tmp end
function tmp_2 = code(x, y, z) t_0 = (cos(y) * z) + (x + sin(y)); tmp = 0.0; if (t_0 <= -2e+28) tmp = z + x; elseif (t_0 <= -0.005) tmp = sin(y); elseif (t_0 <= 5e-12) tmp = y + (z + x); elseif (t_0 <= 2e+18) tmp = sin(y); else tmp = z + x; end tmp_2 = tmp; end
code[x_, y_, z_] := Block[{t$95$0 = N[(N[(N[Cos[y], $MachinePrecision] * z), $MachinePrecision] + N[(x + N[Sin[y], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$0, -2e+28], N[(z + x), $MachinePrecision], If[LessEqual[t$95$0, -0.005], N[Sin[y], $MachinePrecision], If[LessEqual[t$95$0, 5e-12], N[(y + N[(z + x), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$0, 2e+18], N[Sin[y], $MachinePrecision], N[(z + x), $MachinePrecision]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \cos y \cdot z + \left(x + \sin y\right)\\
\mathbf{if}\;t\_0 \leq -2 \cdot 10^{+28}:\\
\;\;\;\;z + x\\
\mathbf{elif}\;t\_0 \leq -0.005:\\
\;\;\;\;\sin y\\
\mathbf{elif}\;t\_0 \leq 5 \cdot 10^{-12}:\\
\;\;\;\;y + \left(z + x\right)\\
\mathbf{elif}\;t\_0 \leq 2 \cdot 10^{+18}:\\
\;\;\;\;\sin y\\
\mathbf{else}:\\
\;\;\;\;z + x\\
\end{array}
\end{array}
if (+.f64 (+.f64 x (sin.f64 y)) (*.f64 z (cos.f64 y))) < -1.99999999999999992e28 or 2e18 < (+.f64 (+.f64 x (sin.f64 y)) (*.f64 z (cos.f64 y))) Initial program 99.9%
Taylor expanded in y around 0
+-commutativeN/A
lower-+.f6480.9
Applied rewrites80.9%
if -1.99999999999999992e28 < (+.f64 (+.f64 x (sin.f64 y)) (*.f64 z (cos.f64 y))) < -0.0050000000000000001 or 4.9999999999999997e-12 < (+.f64 (+.f64 x (sin.f64 y)) (*.f64 z (cos.f64 y))) < 2e18Initial program 99.9%
Taylor expanded in x around 0
+-commutativeN/A
lower-fma.f64N/A
lower-cos.f64N/A
lower-sin.f6494.1
Applied rewrites94.1%
Taylor expanded in z around 0
lower-sin.f6486.2
Applied rewrites86.2%
if -0.0050000000000000001 < (+.f64 (+.f64 x (sin.f64 y)) (*.f64 z (cos.f64 y))) < 4.9999999999999997e-12Initial program 100.0%
Taylor expanded in y around 0
+-commutativeN/A
associate-+l+N/A
+-commutativeN/A
lower-+.f64N/A
+-commutativeN/A
lower-+.f64100.0
Applied rewrites100.0%
Final simplification83.3%
(FPCore (x y z)
:precision binary64
(let* ((t_0 (* x (+ (/ (* (cos y) z) x) 1.0))))
(if (<= x -14500000.0)
t_0
(if (<= x 2.8e-23) (fma z (cos y) (sin y)) t_0))))
double code(double x, double y, double z) {
double t_0 = x * (((cos(y) * z) / x) + 1.0);
double tmp;
if (x <= -14500000.0) {
tmp = t_0;
} else if (x <= 2.8e-23) {
tmp = fma(z, cos(y), sin(y));
} else {
tmp = t_0;
}
return tmp;
}
function code(x, y, z) t_0 = Float64(x * Float64(Float64(Float64(cos(y) * z) / x) + 1.0)) tmp = 0.0 if (x <= -14500000.0) tmp = t_0; elseif (x <= 2.8e-23) tmp = fma(z, cos(y), sin(y)); else tmp = t_0; end return tmp end
code[x_, y_, z_] := Block[{t$95$0 = N[(x * N[(N[(N[(N[Cos[y], $MachinePrecision] * z), $MachinePrecision] / x), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[x, -14500000.0], t$95$0, If[LessEqual[x, 2.8e-23], N[(z * N[Cos[y], $MachinePrecision] + N[Sin[y], $MachinePrecision]), $MachinePrecision], t$95$0]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := x \cdot \left(\frac{\cos y \cdot z}{x} + 1\right)\\
\mathbf{if}\;x \leq -14500000:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;x \leq 2.8 \cdot 10^{-23}:\\
\;\;\;\;\mathsf{fma}\left(z, \cos y, \sin y\right)\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
if x < -1.45e7 or 2.7999999999999997e-23 < x Initial program 99.9%
lift-sin.f64N/A
lift-+.f64N/A
lift-cos.f64N/A
lift-*.f64N/A
+-commutativeN/A
lift-*.f64N/A
*-commutativeN/A
lower-fma.f64100.0
Applied rewrites100.0%
Taylor expanded in x around -inf
mul-1-negN/A
distribute-rgt-neg-inN/A
mul-1-negN/A
lower-*.f64N/A
sub-negN/A
metadata-evalN/A
distribute-lft-inN/A
neg-mul-1N/A
mul-1-negN/A
remove-double-negN/A
metadata-evalN/A
lower-+.f64N/A
Applied rewrites99.9%
Taylor expanded in z around inf
lower-/.f64N/A
lower-*.f64N/A
lower-cos.f6499.1
Applied rewrites99.1%
if -1.45e7 < x < 2.7999999999999997e-23Initial program 99.8%
Taylor expanded in x around 0
+-commutativeN/A
lower-fma.f64N/A
lower-cos.f64N/A
lower-sin.f6495.3
Applied rewrites95.3%
Final simplification97.4%
(FPCore (x y z)
:precision binary64
(let* ((t_0 (* (cos y) z)) (t_1 (* x (+ (/ t_0 x) 1.0))))
(if (<= z -8.5e+67)
t_0
(if (<= z -1.45e-25)
t_1
(if (<= z 3.2e-20) (+ x (sin y)) (if (<= z 6.8e+152) t_1 t_0))))))
double code(double x, double y, double z) {
double t_0 = cos(y) * z;
double t_1 = x * ((t_0 / x) + 1.0);
double tmp;
if (z <= -8.5e+67) {
tmp = t_0;
} else if (z <= -1.45e-25) {
tmp = t_1;
} else if (z <= 3.2e-20) {
tmp = x + sin(y);
} else if (z <= 6.8e+152) {
tmp = t_1;
} 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) :: t_1
real(8) :: tmp
t_0 = cos(y) * z
t_1 = x * ((t_0 / x) + 1.0d0)
if (z <= (-8.5d+67)) then
tmp = t_0
else if (z <= (-1.45d-25)) then
tmp = t_1
else if (z <= 3.2d-20) then
tmp = x + sin(y)
else if (z <= 6.8d+152) then
tmp = t_1
else
tmp = t_0
end if
code = tmp
end function
public static double code(double x, double y, double z) {
double t_0 = Math.cos(y) * z;
double t_1 = x * ((t_0 / x) + 1.0);
double tmp;
if (z <= -8.5e+67) {
tmp = t_0;
} else if (z <= -1.45e-25) {
tmp = t_1;
} else if (z <= 3.2e-20) {
tmp = x + Math.sin(y);
} else if (z <= 6.8e+152) {
tmp = t_1;
} else {
tmp = t_0;
}
return tmp;
}
def code(x, y, z): t_0 = math.cos(y) * z t_1 = x * ((t_0 / x) + 1.0) tmp = 0 if z <= -8.5e+67: tmp = t_0 elif z <= -1.45e-25: tmp = t_1 elif z <= 3.2e-20: tmp = x + math.sin(y) elif z <= 6.8e+152: tmp = t_1 else: tmp = t_0 return tmp
function code(x, y, z) t_0 = Float64(cos(y) * z) t_1 = Float64(x * Float64(Float64(t_0 / x) + 1.0)) tmp = 0.0 if (z <= -8.5e+67) tmp = t_0; elseif (z <= -1.45e-25) tmp = t_1; elseif (z <= 3.2e-20) tmp = Float64(x + sin(y)); elseif (z <= 6.8e+152) tmp = t_1; else tmp = t_0; end return tmp end
function tmp_2 = code(x, y, z) t_0 = cos(y) * z; t_1 = x * ((t_0 / x) + 1.0); tmp = 0.0; if (z <= -8.5e+67) tmp = t_0; elseif (z <= -1.45e-25) tmp = t_1; elseif (z <= 3.2e-20) tmp = x + sin(y); elseif (z <= 6.8e+152) tmp = t_1; else tmp = t_0; end tmp_2 = tmp; end
code[x_, y_, z_] := Block[{t$95$0 = N[(N[Cos[y], $MachinePrecision] * z), $MachinePrecision]}, Block[{t$95$1 = N[(x * N[(N[(t$95$0 / x), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[z, -8.5e+67], t$95$0, If[LessEqual[z, -1.45e-25], t$95$1, If[LessEqual[z, 3.2e-20], N[(x + N[Sin[y], $MachinePrecision]), $MachinePrecision], If[LessEqual[z, 6.8e+152], t$95$1, t$95$0]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \cos y \cdot z\\
t_1 := x \cdot \left(\frac{t\_0}{x} + 1\right)\\
\mathbf{if}\;z \leq -8.5 \cdot 10^{+67}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;z \leq -1.45 \cdot 10^{-25}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;z \leq 3.2 \cdot 10^{-20}:\\
\;\;\;\;x + \sin y\\
\mathbf{elif}\;z \leq 6.8 \cdot 10^{+152}:\\
\;\;\;\;t\_1\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
if z < -8.50000000000000038e67 or 6.80000000000000041e152 < z Initial program 99.8%
Taylor expanded in z around inf
lower-*.f64N/A
lower-cos.f6483.6
Applied rewrites83.6%
if -8.50000000000000038e67 < z < -1.45e-25 or 3.1999999999999997e-20 < z < 6.80000000000000041e152Initial program 99.8%
lift-sin.f64N/A
lift-+.f64N/A
lift-cos.f64N/A
lift-*.f64N/A
+-commutativeN/A
lift-*.f64N/A
*-commutativeN/A
lower-fma.f6499.8
Applied rewrites99.8%
Taylor expanded in x around -inf
mul-1-negN/A
distribute-rgt-neg-inN/A
mul-1-negN/A
lower-*.f64N/A
sub-negN/A
metadata-evalN/A
distribute-lft-inN/A
neg-mul-1N/A
mul-1-negN/A
remove-double-negN/A
metadata-evalN/A
lower-+.f64N/A
Applied rewrites94.3%
Taylor expanded in z around inf
lower-/.f64N/A
lower-*.f64N/A
lower-cos.f6490.8
Applied rewrites90.8%
if -1.45e-25 < z < 3.1999999999999997e-20Initial program 100.0%
Taylor expanded in z around 0
+-commutativeN/A
lower-+.f64N/A
lower-sin.f6493.6
Applied rewrites93.6%
Final simplification89.7%
(FPCore (x y z)
:precision binary64
(let* ((t_0 (* (cos y) z)))
(if (<= z -9e+230)
t_0
(if (<= z -1.6e-25) (+ z x) (if (<= z 3.85e+65) (+ x (sin y)) t_0)))))
double code(double x, double y, double z) {
double t_0 = cos(y) * z;
double tmp;
if (z <= -9e+230) {
tmp = t_0;
} else if (z <= -1.6e-25) {
tmp = z + x;
} else if (z <= 3.85e+65) {
tmp = x + sin(y);
} 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 = cos(y) * z
if (z <= (-9d+230)) then
tmp = t_0
else if (z <= (-1.6d-25)) then
tmp = z + x
else if (z <= 3.85d+65) then
tmp = x + sin(y)
else
tmp = t_0
end if
code = tmp
end function
public static double code(double x, double y, double z) {
double t_0 = Math.cos(y) * z;
double tmp;
if (z <= -9e+230) {
tmp = t_0;
} else if (z <= -1.6e-25) {
tmp = z + x;
} else if (z <= 3.85e+65) {
tmp = x + Math.sin(y);
} else {
tmp = t_0;
}
return tmp;
}
def code(x, y, z): t_0 = math.cos(y) * z tmp = 0 if z <= -9e+230: tmp = t_0 elif z <= -1.6e-25: tmp = z + x elif z <= 3.85e+65: tmp = x + math.sin(y) else: tmp = t_0 return tmp
function code(x, y, z) t_0 = Float64(cos(y) * z) tmp = 0.0 if (z <= -9e+230) tmp = t_0; elseif (z <= -1.6e-25) tmp = Float64(z + x); elseif (z <= 3.85e+65) tmp = Float64(x + sin(y)); else tmp = t_0; end return tmp end
function tmp_2 = code(x, y, z) t_0 = cos(y) * z; tmp = 0.0; if (z <= -9e+230) tmp = t_0; elseif (z <= -1.6e-25) tmp = z + x; elseif (z <= 3.85e+65) tmp = x + sin(y); else tmp = t_0; end tmp_2 = tmp; end
code[x_, y_, z_] := Block[{t$95$0 = N[(N[Cos[y], $MachinePrecision] * z), $MachinePrecision]}, If[LessEqual[z, -9e+230], t$95$0, If[LessEqual[z, -1.6e-25], N[(z + x), $MachinePrecision], If[LessEqual[z, 3.85e+65], N[(x + N[Sin[y], $MachinePrecision]), $MachinePrecision], t$95$0]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \cos y \cdot z\\
\mathbf{if}\;z \leq -9 \cdot 10^{+230}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;z \leq -1.6 \cdot 10^{-25}:\\
\;\;\;\;z + x\\
\mathbf{elif}\;z \leq 3.85 \cdot 10^{+65}:\\
\;\;\;\;x + \sin y\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
if z < -8.9999999999999998e230 or 3.85000000000000019e65 < z Initial program 99.8%
Taylor expanded in z around inf
lower-*.f64N/A
lower-cos.f6487.0
Applied rewrites87.0%
if -8.9999999999999998e230 < z < -1.6000000000000001e-25Initial program 99.9%
Taylor expanded in y around 0
+-commutativeN/A
lower-+.f6477.1
Applied rewrites77.1%
if -1.6000000000000001e-25 < z < 3.85000000000000019e65Initial program 100.0%
Taylor expanded in z around 0
+-commutativeN/A
lower-+.f64N/A
lower-sin.f6490.5
Applied rewrites90.5%
Final simplification86.7%
(FPCore (x y z) :precision binary64 (let* ((t_0 (+ x (sin y)))) (if (<= y -2.25e+50) t_0 (if (<= y 0.0305) (+ y (+ z x)) t_0))))
double code(double x, double y, double z) {
double t_0 = x + sin(y);
double tmp;
if (y <= -2.25e+50) {
tmp = t_0;
} else if (y <= 0.0305) {
tmp = y + (z + x);
} 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 = x + sin(y)
if (y <= (-2.25d+50)) then
tmp = t_0
else if (y <= 0.0305d0) then
tmp = y + (z + x)
else
tmp = t_0
end if
code = tmp
end function
public static double code(double x, double y, double z) {
double t_0 = x + Math.sin(y);
double tmp;
if (y <= -2.25e+50) {
tmp = t_0;
} else if (y <= 0.0305) {
tmp = y + (z + x);
} else {
tmp = t_0;
}
return tmp;
}
def code(x, y, z): t_0 = x + math.sin(y) tmp = 0 if y <= -2.25e+50: tmp = t_0 elif y <= 0.0305: tmp = y + (z + x) else: tmp = t_0 return tmp
function code(x, y, z) t_0 = Float64(x + sin(y)) tmp = 0.0 if (y <= -2.25e+50) tmp = t_0; elseif (y <= 0.0305) tmp = Float64(y + Float64(z + x)); else tmp = t_0; end return tmp end
function tmp_2 = code(x, y, z) t_0 = x + sin(y); tmp = 0.0; if (y <= -2.25e+50) tmp = t_0; elseif (y <= 0.0305) tmp = y + (z + x); else tmp = t_0; end tmp_2 = tmp; end
code[x_, y_, z_] := Block[{t$95$0 = N[(x + N[Sin[y], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -2.25e+50], t$95$0, If[LessEqual[y, 0.0305], N[(y + N[(z + x), $MachinePrecision]), $MachinePrecision], t$95$0]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := x + \sin y\\
\mathbf{if}\;y \leq -2.25 \cdot 10^{+50}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;y \leq 0.0305:\\
\;\;\;\;y + \left(z + x\right)\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
if y < -2.25000000000000007e50 or 0.030499999999999999 < y Initial program 99.8%
Taylor expanded in z around 0
+-commutativeN/A
lower-+.f64N/A
lower-sin.f6464.8
Applied rewrites64.8%
if -2.25000000000000007e50 < y < 0.030499999999999999Initial program 100.0%
Taylor expanded in y around 0
+-commutativeN/A
associate-+l+N/A
+-commutativeN/A
lower-+.f64N/A
+-commutativeN/A
lower-+.f6498.1
Applied rewrites98.1%
Final simplification82.9%
(FPCore (x y z) :precision binary64 (+ z x))
double code(double x, double y, double z) {
return z + x;
}
real(8) function code(x, y, z)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
code = z + x
end function
public static double code(double x, double y, double z) {
return z + x;
}
def code(x, y, z): return z + x
function code(x, y, z) return Float64(z + x) end
function tmp = code(x, y, z) tmp = z + x; end
code[x_, y_, z_] := N[(z + x), $MachinePrecision]
\begin{array}{l}
\\
z + x
\end{array}
Initial program 99.9%
Taylor expanded in y around 0
+-commutativeN/A
lower-+.f6470.6
Applied rewrites70.6%
(FPCore (x y z) :precision binary64 x)
double code(double x, double y, double z) {
return x;
}
real(8) function code(x, y, z)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
code = x
end function
public static double code(double x, double y, double z) {
return x;
}
def code(x, y, z): return x
function code(x, y, z) return x end
function tmp = code(x, y, z) tmp = x; end
code[x_, y_, z_] := x
\begin{array}{l}
\\
x
\end{array}
Initial program 99.9%
lift-sin.f64N/A
lift-+.f64N/A
lift-cos.f64N/A
lift-*.f64N/A
+-commutativeN/A
lift-*.f64N/A
*-commutativeN/A
lower-fma.f6499.9
Applied rewrites99.9%
Taylor expanded in x around -inf
mul-1-negN/A
distribute-rgt-neg-inN/A
mul-1-negN/A
lower-*.f64N/A
sub-negN/A
metadata-evalN/A
distribute-lft-inN/A
neg-mul-1N/A
mul-1-negN/A
remove-double-negN/A
metadata-evalN/A
lower-+.f64N/A
Applied rewrites85.9%
Taylor expanded in x around inf
Applied rewrites43.7%
*-rgt-identity43.7
Applied rewrites43.7%
herbie shell --seed 2024219
(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))))