
(FPCore (x y z) :precision binary64 (+ x (* (* (- y x) 6.0) (- (/ 2.0 3.0) z))))
double code(double x, double y, double z) {
return x + (((y - x) * 6.0) * ((2.0 / 3.0) - 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 + (((y - x) * 6.0d0) * ((2.0d0 / 3.0d0) - z))
end function
public static double code(double x, double y, double z) {
return x + (((y - x) * 6.0) * ((2.0 / 3.0) - z));
}
def code(x, y, z): return x + (((y - x) * 6.0) * ((2.0 / 3.0) - z))
function code(x, y, z) return Float64(x + Float64(Float64(Float64(y - x) * 6.0) * Float64(Float64(2.0 / 3.0) - z))) end
function tmp = code(x, y, z) tmp = x + (((y - x) * 6.0) * ((2.0 / 3.0) - z)); end
code[x_, y_, z_] := N[(x + N[(N[(N[(y - x), $MachinePrecision] * 6.0), $MachinePrecision] * N[(N[(2.0 / 3.0), $MachinePrecision] - z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
x + \left(\left(y - x\right) \cdot 6\right) \cdot \left(\frac{2}{3} - z\right)
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 11 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (x y z) :precision binary64 (+ x (* (* (- y x) 6.0) (- (/ 2.0 3.0) z))))
double code(double x, double y, double z) {
return x + (((y - x) * 6.0) * ((2.0 / 3.0) - 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 + (((y - x) * 6.0d0) * ((2.0d0 / 3.0d0) - z))
end function
public static double code(double x, double y, double z) {
return x + (((y - x) * 6.0) * ((2.0 / 3.0) - z));
}
def code(x, y, z): return x + (((y - x) * 6.0) * ((2.0 / 3.0) - z))
function code(x, y, z) return Float64(x + Float64(Float64(Float64(y - x) * 6.0) * Float64(Float64(2.0 / 3.0) - z))) end
function tmp = code(x, y, z) tmp = x + (((y - x) * 6.0) * ((2.0 / 3.0) - z)); end
code[x_, y_, z_] := N[(x + N[(N[(N[(y - x), $MachinePrecision] * 6.0), $MachinePrecision] * N[(N[(2.0 / 3.0), $MachinePrecision] - z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
x + \left(\left(y - x\right) \cdot 6\right) \cdot \left(\frac{2}{3} - z\right)
\end{array}
(FPCore (x y z) :precision binary64 (+ x (* (* (- y x) 6.0) (- (/ 2.0 3.0) z))))
double code(double x, double y, double z) {
return x + (((y - x) * 6.0) * ((2.0 / 3.0) - 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 + (((y - x) * 6.0d0) * ((2.0d0 / 3.0d0) - z))
end function
public static double code(double x, double y, double z) {
return x + (((y - x) * 6.0) * ((2.0 / 3.0) - z));
}
def code(x, y, z): return x + (((y - x) * 6.0) * ((2.0 / 3.0) - z))
function code(x, y, z) return Float64(x + Float64(Float64(Float64(y - x) * 6.0) * Float64(Float64(2.0 / 3.0) - z))) end
function tmp = code(x, y, z) tmp = x + (((y - x) * 6.0) * ((2.0 / 3.0) - z)); end
code[x_, y_, z_] := N[(x + N[(N[(N[(y - x), $MachinePrecision] * 6.0), $MachinePrecision] * N[(N[(2.0 / 3.0), $MachinePrecision] - z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
x + \left(\left(y - x\right) \cdot 6\right) \cdot \left(\frac{2}{3} - z\right)
\end{array}
Initial program 99.6%
(FPCore (x y z) :precision binary64 (let* ((t_0 (- (/ 2.0 3.0) z)) (t_1 (* 6.0 (* z (- x y))))) (if (<= t_0 -10000000.0) t_1 (if (<= t_0 1.0) (fma 4.0 (- y x) x) t_1))))
double code(double x, double y, double z) {
double t_0 = (2.0 / 3.0) - z;
double t_1 = 6.0 * (z * (x - y));
double tmp;
if (t_0 <= -10000000.0) {
tmp = t_1;
} else if (t_0 <= 1.0) {
tmp = fma(4.0, (y - x), x);
} else {
tmp = t_1;
}
return tmp;
}
function code(x, y, z) t_0 = Float64(Float64(2.0 / 3.0) - z) t_1 = Float64(6.0 * Float64(z * Float64(x - y))) tmp = 0.0 if (t_0 <= -10000000.0) tmp = t_1; elseif (t_0 <= 1.0) tmp = fma(4.0, Float64(y - x), x); else tmp = t_1; end return tmp end
code[x_, y_, z_] := Block[{t$95$0 = N[(N[(2.0 / 3.0), $MachinePrecision] - z), $MachinePrecision]}, Block[{t$95$1 = N[(6.0 * N[(z * N[(x - y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$0, -10000000.0], t$95$1, If[LessEqual[t$95$0, 1.0], N[(4.0 * N[(y - x), $MachinePrecision] + x), $MachinePrecision], t$95$1]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \frac{2}{3} - z\\
t_1 := 6 \cdot \left(z \cdot \left(x - y\right)\right)\\
\mathbf{if}\;t\_0 \leq -10000000:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;t\_0 \leq 1:\\
\;\;\;\;\mathsf{fma}\left(4, y - x, x\right)\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if (-.f64 (/.f64 #s(literal 2 binary64) #s(literal 3 binary64)) z) < -1e7 or 1 < (-.f64 (/.f64 #s(literal 2 binary64) #s(literal 3 binary64)) z) Initial program 99.8%
Taylor expanded in z around inf
metadata-evalN/A
distribute-lft-neg-inN/A
distribute-rgt-neg-inN/A
lower-*.f64N/A
distribute-rgt-neg-inN/A
neg-mul-1N/A
lower-*.f64N/A
sub-negN/A
neg-mul-1N/A
+-commutativeN/A
distribute-lft-inN/A
associate-*r*N/A
metadata-evalN/A
*-lft-identityN/A
mul-1-negN/A
sub-negN/A
lower--.f6499.3
Applied rewrites99.3%
if -1e7 < (-.f64 (/.f64 #s(literal 2 binary64) #s(literal 3 binary64)) z) < 1Initial program 99.4%
Taylor expanded in z around 0
+-commutativeN/A
lower-fma.f64N/A
lower--.f6498.0
Applied rewrites98.0%
(FPCore (x y z)
:precision binary64
(let* ((t_0 (- (/ 2.0 3.0) z)) (t_1 (* x (fma 6.0 z -3.0))))
(if (<= t_0 -10000000.0)
t_1
(if (<= t_0 0.66666666667) (fma 4.0 (- y x) x) t_1))))
double code(double x, double y, double z) {
double t_0 = (2.0 / 3.0) - z;
double t_1 = x * fma(6.0, z, -3.0);
double tmp;
if (t_0 <= -10000000.0) {
tmp = t_1;
} else if (t_0 <= 0.66666666667) {
tmp = fma(4.0, (y - x), x);
} else {
tmp = t_1;
}
return tmp;
}
function code(x, y, z) t_0 = Float64(Float64(2.0 / 3.0) - z) t_1 = Float64(x * fma(6.0, z, -3.0)) tmp = 0.0 if (t_0 <= -10000000.0) tmp = t_1; elseif (t_0 <= 0.66666666667) tmp = fma(4.0, Float64(y - x), x); else tmp = t_1; end return tmp end
code[x_, y_, z_] := Block[{t$95$0 = N[(N[(2.0 / 3.0), $MachinePrecision] - z), $MachinePrecision]}, Block[{t$95$1 = N[(x * N[(6.0 * z + -3.0), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$0, -10000000.0], t$95$1, If[LessEqual[t$95$0, 0.66666666667], N[(4.0 * N[(y - x), $MachinePrecision] + x), $MachinePrecision], t$95$1]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \frac{2}{3} - z\\
t_1 := x \cdot \mathsf{fma}\left(6, z, -3\right)\\
\mathbf{if}\;t\_0 \leq -10000000:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;t\_0 \leq 0.66666666667:\\
\;\;\;\;\mathsf{fma}\left(4, y - x, x\right)\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if (-.f64 (/.f64 #s(literal 2 binary64) #s(literal 3 binary64)) z) < -1e7 or 0.666666666669999963 < (-.f64 (/.f64 #s(literal 2 binary64) #s(literal 3 binary64)) z) Initial program 99.8%
Taylor expanded in x around inf
remove-double-negN/A
neg-mul-1N/A
distribute-lft-neg-inN/A
distribute-rgt-neg-inN/A
distribute-neg-inN/A
metadata-evalN/A
distribute-lft-neg-inN/A
metadata-evalN/A
+-commutativeN/A
metadata-evalN/A
sub-negN/A
*-commutativeN/A
associate-*l*N/A
neg-mul-1N/A
lower-*.f64N/A
sub-negN/A
metadata-evalN/A
distribute-neg-inN/A
distribute-lft-neg-inN/A
metadata-evalN/A
Applied rewrites59.0%
if -1e7 < (-.f64 (/.f64 #s(literal 2 binary64) #s(literal 3 binary64)) z) < 0.666666666669999963Initial program 99.4%
Taylor expanded in z around 0
+-commutativeN/A
lower-fma.f64N/A
lower--.f6499.2
Applied rewrites99.2%
(FPCore (x y z)
:precision binary64
(let* ((t_0 (- (/ 2.0 3.0) z)))
(if (<= t_0 -10000000.0)
(* z (* x 6.0))
(if (<= t_0 2e+19) (fma 4.0 (- y x) x) (* x (* 6.0 z))))))
double code(double x, double y, double z) {
double t_0 = (2.0 / 3.0) - z;
double tmp;
if (t_0 <= -10000000.0) {
tmp = z * (x * 6.0);
} else if (t_0 <= 2e+19) {
tmp = fma(4.0, (y - x), x);
} else {
tmp = x * (6.0 * z);
}
return tmp;
}
function code(x, y, z) t_0 = Float64(Float64(2.0 / 3.0) - z) tmp = 0.0 if (t_0 <= -10000000.0) tmp = Float64(z * Float64(x * 6.0)); elseif (t_0 <= 2e+19) tmp = fma(4.0, Float64(y - x), x); else tmp = Float64(x * Float64(6.0 * z)); end return tmp end
code[x_, y_, z_] := Block[{t$95$0 = N[(N[(2.0 / 3.0), $MachinePrecision] - z), $MachinePrecision]}, If[LessEqual[t$95$0, -10000000.0], N[(z * N[(x * 6.0), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$0, 2e+19], N[(4.0 * N[(y - x), $MachinePrecision] + x), $MachinePrecision], N[(x * N[(6.0 * z), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \frac{2}{3} - z\\
\mathbf{if}\;t\_0 \leq -10000000:\\
\;\;\;\;z \cdot \left(x \cdot 6\right)\\
\mathbf{elif}\;t\_0 \leq 2 \cdot 10^{+19}:\\
\;\;\;\;\mathsf{fma}\left(4, y - x, x\right)\\
\mathbf{else}:\\
\;\;\;\;x \cdot \left(6 \cdot z\right)\\
\end{array}
\end{array}
if (-.f64 (/.f64 #s(literal 2 binary64) #s(literal 3 binary64)) z) < -1e7Initial program 99.8%
Taylor expanded in y around 0
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
lower-*.f64N/A
sub-negN/A
mul-1-negN/A
distribute-lft-inN/A
metadata-evalN/A
+-commutativeN/A
associate-*r*N/A
metadata-evalN/A
lower-fma.f6457.5
Applied rewrites57.5%
Taylor expanded in z around inf
*-commutativeN/A
associate-*l*N/A
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f6456.8
Applied rewrites56.8%
*-commutativeN/A
associate-*r*N/A
lower-*.f64N/A
lower-*.f6456.9
Applied rewrites56.9%
if -1e7 < (-.f64 (/.f64 #s(literal 2 binary64) #s(literal 3 binary64)) z) < 2e19Initial program 99.4%
Taylor expanded in z around 0
+-commutativeN/A
lower-fma.f64N/A
lower--.f6497.3
Applied rewrites97.3%
if 2e19 < (-.f64 (/.f64 #s(literal 2 binary64) #s(literal 3 binary64)) z) Initial program 99.7%
Taylor expanded in y around 0
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
lower-*.f64N/A
sub-negN/A
mul-1-negN/A
distribute-lft-inN/A
metadata-evalN/A
+-commutativeN/A
associate-*r*N/A
metadata-evalN/A
lower-fma.f6461.5
Applied rewrites61.5%
Taylor expanded in z around inf
*-commutativeN/A
associate-*l*N/A
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f6461.5
Applied rewrites61.5%
Final simplification78.0%
(FPCore (x y z) :precision binary64 (let* ((t_0 (- (/ 2.0 3.0) z)) (t_1 (* x (* 6.0 z)))) (if (<= t_0 -10000000.0) t_1 (if (<= t_0 2e+19) (fma 4.0 (- y x) x) t_1))))
double code(double x, double y, double z) {
double t_0 = (2.0 / 3.0) - z;
double t_1 = x * (6.0 * z);
double tmp;
if (t_0 <= -10000000.0) {
tmp = t_1;
} else if (t_0 <= 2e+19) {
tmp = fma(4.0, (y - x), x);
} else {
tmp = t_1;
}
return tmp;
}
function code(x, y, z) t_0 = Float64(Float64(2.0 / 3.0) - z) t_1 = Float64(x * Float64(6.0 * z)) tmp = 0.0 if (t_0 <= -10000000.0) tmp = t_1; elseif (t_0 <= 2e+19) tmp = fma(4.0, Float64(y - x), x); else tmp = t_1; end return tmp end
code[x_, y_, z_] := Block[{t$95$0 = N[(N[(2.0 / 3.0), $MachinePrecision] - z), $MachinePrecision]}, Block[{t$95$1 = N[(x * N[(6.0 * z), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$0, -10000000.0], t$95$1, If[LessEqual[t$95$0, 2e+19], N[(4.0 * N[(y - x), $MachinePrecision] + x), $MachinePrecision], t$95$1]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \frac{2}{3} - z\\
t_1 := x \cdot \left(6 \cdot z\right)\\
\mathbf{if}\;t\_0 \leq -10000000:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;t\_0 \leq 2 \cdot 10^{+19}:\\
\;\;\;\;\mathsf{fma}\left(4, y - x, x\right)\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if (-.f64 (/.f64 #s(literal 2 binary64) #s(literal 3 binary64)) z) < -1e7 or 2e19 < (-.f64 (/.f64 #s(literal 2 binary64) #s(literal 3 binary64)) z) Initial program 99.8%
Taylor expanded in y around 0
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
lower-*.f64N/A
sub-negN/A
mul-1-negN/A
distribute-lft-inN/A
metadata-evalN/A
+-commutativeN/A
associate-*r*N/A
metadata-evalN/A
lower-fma.f6459.3
Applied rewrites59.3%
Taylor expanded in z around inf
*-commutativeN/A
associate-*l*N/A
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f6458.9
Applied rewrites58.9%
if -1e7 < (-.f64 (/.f64 #s(literal 2 binary64) #s(literal 3 binary64)) z) < 2e19Initial program 99.4%
Taylor expanded in z around 0
+-commutativeN/A
lower-fma.f64N/A
lower--.f6497.3
Applied rewrites97.3%
Final simplification77.9%
(FPCore (x y z)
:precision binary64
(let* ((t_0 (* y (fma z -6.0 4.0))))
(if (<= y -1.26e+47)
t_0
(if (<= y 3.1e+103) (fma (fma z 6.0 -4.0) x x) t_0))))
double code(double x, double y, double z) {
double t_0 = y * fma(z, -6.0, 4.0);
double tmp;
if (y <= -1.26e+47) {
tmp = t_0;
} else if (y <= 3.1e+103) {
tmp = fma(fma(z, 6.0, -4.0), x, x);
} else {
tmp = t_0;
}
return tmp;
}
function code(x, y, z) t_0 = Float64(y * fma(z, -6.0, 4.0)) tmp = 0.0 if (y <= -1.26e+47) tmp = t_0; elseif (y <= 3.1e+103) tmp = fma(fma(z, 6.0, -4.0), x, x); else tmp = t_0; end return tmp end
code[x_, y_, z_] := Block[{t$95$0 = N[(y * N[(z * -6.0 + 4.0), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -1.26e+47], t$95$0, If[LessEqual[y, 3.1e+103], N[(N[(z * 6.0 + -4.0), $MachinePrecision] * x + x), $MachinePrecision], t$95$0]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := y \cdot \mathsf{fma}\left(z, -6, 4\right)\\
\mathbf{if}\;y \leq -1.26 \cdot 10^{+47}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;y \leq 3.1 \cdot 10^{+103}:\\
\;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(z, 6, -4\right), x, x\right)\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
if y < -1.26e47 or 3.1000000000000002e103 < y Initial program 99.7%
Taylor expanded in x around 0
associate-*r*N/A
*-commutativeN/A
associate-*r*N/A
lower-*.f64N/A
sub-negN/A
mul-1-negN/A
+-commutativeN/A
distribute-lft-inN/A
associate-*r*N/A
metadata-evalN/A
*-commutativeN/A
metadata-evalN/A
lower-fma.f6489.5
Applied rewrites89.5%
if -1.26e47 < y < 3.1000000000000002e103Initial program 99.5%
Taylor expanded in y around 0
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
lower-*.f64N/A
sub-negN/A
mul-1-negN/A
distribute-lft-inN/A
metadata-evalN/A
+-commutativeN/A
associate-*r*N/A
metadata-evalN/A
lower-fma.f6478.7
Applied rewrites78.7%
lift-fma.f64N/A
*-commutativeN/A
distribute-rgt1-inN/A
distribute-lft1-inN/A
lower-fma.f6478.7
lift-fma.f64N/A
*-commutativeN/A
lower-fma.f6478.7
Applied rewrites78.7%
(FPCore (x y z) :precision binary64 (let* ((t_0 (* y (fma z -6.0 4.0)))) (if (<= y -1.26e+47) t_0 (if (<= y 3.1e+103) (* x (fma 6.0 z -3.0)) t_0))))
double code(double x, double y, double z) {
double t_0 = y * fma(z, -6.0, 4.0);
double tmp;
if (y <= -1.26e+47) {
tmp = t_0;
} else if (y <= 3.1e+103) {
tmp = x * fma(6.0, z, -3.0);
} else {
tmp = t_0;
}
return tmp;
}
function code(x, y, z) t_0 = Float64(y * fma(z, -6.0, 4.0)) tmp = 0.0 if (y <= -1.26e+47) tmp = t_0; elseif (y <= 3.1e+103) tmp = Float64(x * fma(6.0, z, -3.0)); else tmp = t_0; end return tmp end
code[x_, y_, z_] := Block[{t$95$0 = N[(y * N[(z * -6.0 + 4.0), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -1.26e+47], t$95$0, If[LessEqual[y, 3.1e+103], N[(x * N[(6.0 * z + -3.0), $MachinePrecision]), $MachinePrecision], t$95$0]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := y \cdot \mathsf{fma}\left(z, -6, 4\right)\\
\mathbf{if}\;y \leq -1.26 \cdot 10^{+47}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;y \leq 3.1 \cdot 10^{+103}:\\
\;\;\;\;x \cdot \mathsf{fma}\left(6, z, -3\right)\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
if y < -1.26e47 or 3.1000000000000002e103 < y Initial program 99.7%
Taylor expanded in x around 0
associate-*r*N/A
*-commutativeN/A
associate-*r*N/A
lower-*.f64N/A
sub-negN/A
mul-1-negN/A
+-commutativeN/A
distribute-lft-inN/A
associate-*r*N/A
metadata-evalN/A
*-commutativeN/A
metadata-evalN/A
lower-fma.f6489.5
Applied rewrites89.5%
if -1.26e47 < y < 3.1000000000000002e103Initial program 99.5%
Taylor expanded in x around inf
remove-double-negN/A
neg-mul-1N/A
distribute-lft-neg-inN/A
distribute-rgt-neg-inN/A
distribute-neg-inN/A
metadata-evalN/A
distribute-lft-neg-inN/A
metadata-evalN/A
+-commutativeN/A
metadata-evalN/A
sub-negN/A
*-commutativeN/A
associate-*l*N/A
neg-mul-1N/A
lower-*.f64N/A
sub-negN/A
metadata-evalN/A
distribute-neg-inN/A
distribute-lft-neg-inN/A
metadata-evalN/A
Applied rewrites78.7%
(FPCore (x y z) :precision binary64 (if (<= y -1e+133) (* y 4.0) (if (<= y 2e+96) (* x -3.0) (* y 4.0))))
double code(double x, double y, double z) {
double tmp;
if (y <= -1e+133) {
tmp = y * 4.0;
} else if (y <= 2e+96) {
tmp = x * -3.0;
} else {
tmp = y * 4.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) :: tmp
if (y <= (-1d+133)) then
tmp = y * 4.0d0
else if (y <= 2d+96) then
tmp = x * (-3.0d0)
else
tmp = y * 4.0d0
end if
code = tmp
end function
public static double code(double x, double y, double z) {
double tmp;
if (y <= -1e+133) {
tmp = y * 4.0;
} else if (y <= 2e+96) {
tmp = x * -3.0;
} else {
tmp = y * 4.0;
}
return tmp;
}
def code(x, y, z): tmp = 0 if y <= -1e+133: tmp = y * 4.0 elif y <= 2e+96: tmp = x * -3.0 else: tmp = y * 4.0 return tmp
function code(x, y, z) tmp = 0.0 if (y <= -1e+133) tmp = Float64(y * 4.0); elseif (y <= 2e+96) tmp = Float64(x * -3.0); else tmp = Float64(y * 4.0); end return tmp end
function tmp_2 = code(x, y, z) tmp = 0.0; if (y <= -1e+133) tmp = y * 4.0; elseif (y <= 2e+96) tmp = x * -3.0; else tmp = y * 4.0; end tmp_2 = tmp; end
code[x_, y_, z_] := If[LessEqual[y, -1e+133], N[(y * 4.0), $MachinePrecision], If[LessEqual[y, 2e+96], N[(x * -3.0), $MachinePrecision], N[(y * 4.0), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq -1 \cdot 10^{+133}:\\
\;\;\;\;y \cdot 4\\
\mathbf{elif}\;y \leq 2 \cdot 10^{+96}:\\
\;\;\;\;x \cdot -3\\
\mathbf{else}:\\
\;\;\;\;y \cdot 4\\
\end{array}
\end{array}
if y < -1e133 or 2.0000000000000001e96 < y Initial program 99.7%
Taylor expanded in x around 0
associate-*r*N/A
*-commutativeN/A
associate-*r*N/A
lower-*.f64N/A
sub-negN/A
mul-1-negN/A
+-commutativeN/A
distribute-lft-inN/A
associate-*r*N/A
metadata-evalN/A
*-commutativeN/A
metadata-evalN/A
lower-fma.f6493.5
Applied rewrites93.5%
Taylor expanded in z around 0
Applied rewrites50.9%
if -1e133 < y < 2.0000000000000001e96Initial program 99.6%
Taylor expanded in z around 0
+-commutativeN/A
lower-fma.f64N/A
lower--.f6448.5
Applied rewrites48.5%
Taylor expanded in y around 0
distribute-rgt1-inN/A
metadata-evalN/A
*-commutativeN/A
lower-*.f6438.3
Applied rewrites38.3%
(FPCore (x y z) :precision binary64 (fma (- 0.6666666666666666 z) (* (- y x) 6.0) x))
double code(double x, double y, double z) {
return fma((0.6666666666666666 - z), ((y - x) * 6.0), x);
}
function code(x, y, z) return fma(Float64(0.6666666666666666 - z), Float64(Float64(y - x) * 6.0), x) end
code[x_, y_, z_] := N[(N[(0.6666666666666666 - z), $MachinePrecision] * N[(N[(y - x), $MachinePrecision] * 6.0), $MachinePrecision] + x), $MachinePrecision]
\begin{array}{l}
\\
\mathsf{fma}\left(0.6666666666666666 - z, \left(y - x\right) \cdot 6, x\right)
\end{array}
Initial program 99.6%
lift--.f64N/A
lift-/.f64N/A
lift--.f64N/A
lift-*.f64N/A
lift-*.f64N/A
+-commutativeN/A
lift-*.f64N/A
*-commutativeN/A
lower-fma.f6499.5
lift-/.f64N/A
metadata-eval99.5
Applied rewrites99.5%
(FPCore (x y z) :precision binary64 (fma 4.0 (- y x) x))
double code(double x, double y, double z) {
return fma(4.0, (y - x), x);
}
function code(x, y, z) return fma(4.0, Float64(y - x), x) end
code[x_, y_, z_] := N[(4.0 * N[(y - x), $MachinePrecision] + x), $MachinePrecision]
\begin{array}{l}
\\
\mathsf{fma}\left(4, y - x, x\right)
\end{array}
Initial program 99.6%
Taylor expanded in z around 0
+-commutativeN/A
lower-fma.f64N/A
lower--.f6449.9
Applied rewrites49.9%
(FPCore (x y z) :precision binary64 (* x -3.0))
double code(double x, double y, double z) {
return x * -3.0;
}
real(8) function code(x, y, z)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
code = x * (-3.0d0)
end function
public static double code(double x, double y, double z) {
return x * -3.0;
}
def code(x, y, z): return x * -3.0
function code(x, y, z) return Float64(x * -3.0) end
function tmp = code(x, y, z) tmp = x * -3.0; end
code[x_, y_, z_] := N[(x * -3.0), $MachinePrecision]
\begin{array}{l}
\\
x \cdot -3
\end{array}
Initial program 99.6%
Taylor expanded in z around 0
+-commutativeN/A
lower-fma.f64N/A
lower--.f6449.9
Applied rewrites49.9%
Taylor expanded in y around 0
distribute-rgt1-inN/A
metadata-evalN/A
*-commutativeN/A
lower-*.f6426.7
Applied rewrites26.7%
herbie shell --seed 2024219
(FPCore (x y z)
:name "Data.Colour.RGBSpace.HSL:hsl from colour-2.3.3, D"
:precision binary64
(+ x (* (* (- y x) 6.0) (- (/ 2.0 3.0) z))))