
(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 13 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 (fma (- y x) (+ 4.0 (* z -6.0)) x))
double code(double x, double y, double z) {
return fma((y - x), (4.0 + (z * -6.0)), x);
}
function code(x, y, z) return fma(Float64(y - x), Float64(4.0 + Float64(z * -6.0)), x) end
code[x_, y_, z_] := N[(N[(y - x), $MachinePrecision] * N[(4.0 + N[(z * -6.0), $MachinePrecision]), $MachinePrecision] + x), $MachinePrecision]
\begin{array}{l}
\\
\mathsf{fma}\left(y - x, 4 + z \cdot -6, x\right)
\end{array}
Initial program 99.6%
+-commutative99.6%
associate-*l*99.8%
fma-def99.8%
sub-neg99.8%
distribute-rgt-in99.8%
metadata-eval99.8%
metadata-eval99.8%
distribute-lft-neg-out99.8%
distribute-rgt-neg-in99.8%
metadata-eval99.8%
Simplified99.8%
Final simplification99.8%
(FPCore (x y z)
:precision binary64
(let* ((t_0 (* 6.0 (* x z))) (t_1 (* -6.0 (* y z))))
(if (<= z -1e+65)
t_0
(if (<= z -0.66)
t_1
(if (<= z -7e-223)
(* y 4.0)
(if (<= z -8.2e-281)
(* x -3.0)
(if (<= z 1.95e-268)
(* y 4.0)
(if (<= z 3.7e-160)
(* x -3.0)
(if (<= z 0.086)
(* y 4.0)
(if (or (<= z 5.2e+197) (not (<= z 2.8e+270)))
t_1
t_0))))))))))
double code(double x, double y, double z) {
double t_0 = 6.0 * (x * z);
double t_1 = -6.0 * (y * z);
double tmp;
if (z <= -1e+65) {
tmp = t_0;
} else if (z <= -0.66) {
tmp = t_1;
} else if (z <= -7e-223) {
tmp = y * 4.0;
} else if (z <= -8.2e-281) {
tmp = x * -3.0;
} else if (z <= 1.95e-268) {
tmp = y * 4.0;
} else if (z <= 3.7e-160) {
tmp = x * -3.0;
} else if (z <= 0.086) {
tmp = y * 4.0;
} else if ((z <= 5.2e+197) || !(z <= 2.8e+270)) {
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 = 6.0d0 * (x * z)
t_1 = (-6.0d0) * (y * z)
if (z <= (-1d+65)) then
tmp = t_0
else if (z <= (-0.66d0)) then
tmp = t_1
else if (z <= (-7d-223)) then
tmp = y * 4.0d0
else if (z <= (-8.2d-281)) then
tmp = x * (-3.0d0)
else if (z <= 1.95d-268) then
tmp = y * 4.0d0
else if (z <= 3.7d-160) then
tmp = x * (-3.0d0)
else if (z <= 0.086d0) then
tmp = y * 4.0d0
else if ((z <= 5.2d+197) .or. (.not. (z <= 2.8d+270))) 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 = 6.0 * (x * z);
double t_1 = -6.0 * (y * z);
double tmp;
if (z <= -1e+65) {
tmp = t_0;
} else if (z <= -0.66) {
tmp = t_1;
} else if (z <= -7e-223) {
tmp = y * 4.0;
} else if (z <= -8.2e-281) {
tmp = x * -3.0;
} else if (z <= 1.95e-268) {
tmp = y * 4.0;
} else if (z <= 3.7e-160) {
tmp = x * -3.0;
} else if (z <= 0.086) {
tmp = y * 4.0;
} else if ((z <= 5.2e+197) || !(z <= 2.8e+270)) {
tmp = t_1;
} else {
tmp = t_0;
}
return tmp;
}
def code(x, y, z): t_0 = 6.0 * (x * z) t_1 = -6.0 * (y * z) tmp = 0 if z <= -1e+65: tmp = t_0 elif z <= -0.66: tmp = t_1 elif z <= -7e-223: tmp = y * 4.0 elif z <= -8.2e-281: tmp = x * -3.0 elif z <= 1.95e-268: tmp = y * 4.0 elif z <= 3.7e-160: tmp = x * -3.0 elif z <= 0.086: tmp = y * 4.0 elif (z <= 5.2e+197) or not (z <= 2.8e+270): tmp = t_1 else: tmp = t_0 return tmp
function code(x, y, z) t_0 = Float64(6.0 * Float64(x * z)) t_1 = Float64(-6.0 * Float64(y * z)) tmp = 0.0 if (z <= -1e+65) tmp = t_0; elseif (z <= -0.66) tmp = t_1; elseif (z <= -7e-223) tmp = Float64(y * 4.0); elseif (z <= -8.2e-281) tmp = Float64(x * -3.0); elseif (z <= 1.95e-268) tmp = Float64(y * 4.0); elseif (z <= 3.7e-160) tmp = Float64(x * -3.0); elseif (z <= 0.086) tmp = Float64(y * 4.0); elseif ((z <= 5.2e+197) || !(z <= 2.8e+270)) tmp = t_1; else tmp = t_0; end return tmp end
function tmp_2 = code(x, y, z) t_0 = 6.0 * (x * z); t_1 = -6.0 * (y * z); tmp = 0.0; if (z <= -1e+65) tmp = t_0; elseif (z <= -0.66) tmp = t_1; elseif (z <= -7e-223) tmp = y * 4.0; elseif (z <= -8.2e-281) tmp = x * -3.0; elseif (z <= 1.95e-268) tmp = y * 4.0; elseif (z <= 3.7e-160) tmp = x * -3.0; elseif (z <= 0.086) tmp = y * 4.0; elseif ((z <= 5.2e+197) || ~((z <= 2.8e+270))) tmp = t_1; else tmp = t_0; end tmp_2 = tmp; end
code[x_, y_, z_] := Block[{t$95$0 = N[(6.0 * N[(x * z), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(-6.0 * N[(y * z), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[z, -1e+65], t$95$0, If[LessEqual[z, -0.66], t$95$1, If[LessEqual[z, -7e-223], N[(y * 4.0), $MachinePrecision], If[LessEqual[z, -8.2e-281], N[(x * -3.0), $MachinePrecision], If[LessEqual[z, 1.95e-268], N[(y * 4.0), $MachinePrecision], If[LessEqual[z, 3.7e-160], N[(x * -3.0), $MachinePrecision], If[LessEqual[z, 0.086], N[(y * 4.0), $MachinePrecision], If[Or[LessEqual[z, 5.2e+197], N[Not[LessEqual[z, 2.8e+270]], $MachinePrecision]], t$95$1, t$95$0]]]]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := 6 \cdot \left(x \cdot z\right)\\
t_1 := -6 \cdot \left(y \cdot z\right)\\
\mathbf{if}\;z \leq -1 \cdot 10^{+65}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;z \leq -0.66:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;z \leq -7 \cdot 10^{-223}:\\
\;\;\;\;y \cdot 4\\
\mathbf{elif}\;z \leq -8.2 \cdot 10^{-281}:\\
\;\;\;\;x \cdot -3\\
\mathbf{elif}\;z \leq 1.95 \cdot 10^{-268}:\\
\;\;\;\;y \cdot 4\\
\mathbf{elif}\;z \leq 3.7 \cdot 10^{-160}:\\
\;\;\;\;x \cdot -3\\
\mathbf{elif}\;z \leq 0.086:\\
\;\;\;\;y \cdot 4\\
\mathbf{elif}\;z \leq 5.2 \cdot 10^{+197} \lor \neg \left(z \leq 2.8 \cdot 10^{+270}\right):\\
\;\;\;\;t\_1\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
if z < -9.9999999999999999e64 or 5.19999999999999975e197 < z < 2.8000000000000001e270Initial program 99.8%
metadata-eval99.8%
Simplified99.8%
Taylor expanded in x around inf 71.2%
sub-neg71.2%
distribute-rgt-in71.2%
metadata-eval71.2%
metadata-eval71.2%
distribute-lft-neg-in71.2%
associate-+r+71.2%
metadata-eval71.2%
metadata-eval71.2%
distribute-rgt-neg-in71.2%
metadata-eval71.2%
Simplified71.2%
Taylor expanded in z around inf 71.2%
if -9.9999999999999999e64 < z < -0.660000000000000031 or 0.085999999999999993 < z < 5.19999999999999975e197 or 2.8000000000000001e270 < z Initial program 99.7%
+-commutative99.7%
associate-*l*99.8%
fma-def99.8%
sub-neg99.8%
distribute-rgt-in99.8%
metadata-eval99.8%
metadata-eval99.8%
distribute-lft-neg-out99.8%
distribute-rgt-neg-in99.8%
metadata-eval99.8%
Simplified99.8%
Taylor expanded in y around inf 69.2%
Taylor expanded in z around inf 68.2%
*-commutative68.2%
Simplified68.2%
if -0.660000000000000031 < z < -7.00000000000000018e-223 or -8.1999999999999998e-281 < z < 1.9499999999999999e-268 or 3.69999999999999977e-160 < z < 0.085999999999999993Initial program 99.4%
+-commutative99.4%
associate-*l*99.8%
fma-def99.9%
sub-neg99.9%
distribute-rgt-in99.9%
metadata-eval99.9%
metadata-eval99.9%
distribute-lft-neg-out99.9%
distribute-rgt-neg-in99.9%
metadata-eval99.9%
Simplified99.9%
Taylor expanded in y around inf 61.0%
Taylor expanded in z around 0 59.0%
*-commutative59.0%
Simplified59.0%
if -7.00000000000000018e-223 < z < -8.1999999999999998e-281 or 1.9499999999999999e-268 < z < 3.69999999999999977e-160Initial program 99.3%
metadata-eval99.3%
Simplified99.3%
Taylor expanded in x around inf 70.0%
sub-neg70.0%
distribute-rgt-in70.0%
metadata-eval70.0%
metadata-eval70.0%
distribute-lft-neg-in70.0%
associate-+r+70.0%
metadata-eval70.0%
metadata-eval70.0%
distribute-rgt-neg-in70.0%
metadata-eval70.0%
Simplified70.0%
Taylor expanded in z around 0 70.0%
*-commutative70.0%
Simplified70.0%
Final simplification65.5%
(FPCore (x y z)
:precision binary64
(let* ((t_0 (* 6.0 (* x z))) (t_1 (* y (* z -6.0))))
(if (<= z -8.5e+64)
t_0
(if (<= z -0.66)
t_1
(if (<= z -4.3e-223)
(* y 4.0)
(if (<= z -9.8e-281)
(* x -3.0)
(if (<= z 1.1e-268)
(* y 4.0)
(if (<= z 4.3e-159)
(* x -3.0)
(if (<= z 0.086)
(* y 4.0)
(if (<= z 6e+194)
t_1
(if (<= z 3.6e+264) t_0 (* -6.0 (* y z)))))))))))))
double code(double x, double y, double z) {
double t_0 = 6.0 * (x * z);
double t_1 = y * (z * -6.0);
double tmp;
if (z <= -8.5e+64) {
tmp = t_0;
} else if (z <= -0.66) {
tmp = t_1;
} else if (z <= -4.3e-223) {
tmp = y * 4.0;
} else if (z <= -9.8e-281) {
tmp = x * -3.0;
} else if (z <= 1.1e-268) {
tmp = y * 4.0;
} else if (z <= 4.3e-159) {
tmp = x * -3.0;
} else if (z <= 0.086) {
tmp = y * 4.0;
} else if (z <= 6e+194) {
tmp = t_1;
} else if (z <= 3.6e+264) {
tmp = t_0;
} else {
tmp = -6.0 * (y * 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) :: t_0
real(8) :: t_1
real(8) :: tmp
t_0 = 6.0d0 * (x * z)
t_1 = y * (z * (-6.0d0))
if (z <= (-8.5d+64)) then
tmp = t_0
else if (z <= (-0.66d0)) then
tmp = t_1
else if (z <= (-4.3d-223)) then
tmp = y * 4.0d0
else if (z <= (-9.8d-281)) then
tmp = x * (-3.0d0)
else if (z <= 1.1d-268) then
tmp = y * 4.0d0
else if (z <= 4.3d-159) then
tmp = x * (-3.0d0)
else if (z <= 0.086d0) then
tmp = y * 4.0d0
else if (z <= 6d+194) then
tmp = t_1
else if (z <= 3.6d+264) then
tmp = t_0
else
tmp = (-6.0d0) * (y * z)
end if
code = tmp
end function
public static double code(double x, double y, double z) {
double t_0 = 6.0 * (x * z);
double t_1 = y * (z * -6.0);
double tmp;
if (z <= -8.5e+64) {
tmp = t_0;
} else if (z <= -0.66) {
tmp = t_1;
} else if (z <= -4.3e-223) {
tmp = y * 4.0;
} else if (z <= -9.8e-281) {
tmp = x * -3.0;
} else if (z <= 1.1e-268) {
tmp = y * 4.0;
} else if (z <= 4.3e-159) {
tmp = x * -3.0;
} else if (z <= 0.086) {
tmp = y * 4.0;
} else if (z <= 6e+194) {
tmp = t_1;
} else if (z <= 3.6e+264) {
tmp = t_0;
} else {
tmp = -6.0 * (y * z);
}
return tmp;
}
def code(x, y, z): t_0 = 6.0 * (x * z) t_1 = y * (z * -6.0) tmp = 0 if z <= -8.5e+64: tmp = t_0 elif z <= -0.66: tmp = t_1 elif z <= -4.3e-223: tmp = y * 4.0 elif z <= -9.8e-281: tmp = x * -3.0 elif z <= 1.1e-268: tmp = y * 4.0 elif z <= 4.3e-159: tmp = x * -3.0 elif z <= 0.086: tmp = y * 4.0 elif z <= 6e+194: tmp = t_1 elif z <= 3.6e+264: tmp = t_0 else: tmp = -6.0 * (y * z) return tmp
function code(x, y, z) t_0 = Float64(6.0 * Float64(x * z)) t_1 = Float64(y * Float64(z * -6.0)) tmp = 0.0 if (z <= -8.5e+64) tmp = t_0; elseif (z <= -0.66) tmp = t_1; elseif (z <= -4.3e-223) tmp = Float64(y * 4.0); elseif (z <= -9.8e-281) tmp = Float64(x * -3.0); elseif (z <= 1.1e-268) tmp = Float64(y * 4.0); elseif (z <= 4.3e-159) tmp = Float64(x * -3.0); elseif (z <= 0.086) tmp = Float64(y * 4.0); elseif (z <= 6e+194) tmp = t_1; elseif (z <= 3.6e+264) tmp = t_0; else tmp = Float64(-6.0 * Float64(y * z)); end return tmp end
function tmp_2 = code(x, y, z) t_0 = 6.0 * (x * z); t_1 = y * (z * -6.0); tmp = 0.0; if (z <= -8.5e+64) tmp = t_0; elseif (z <= -0.66) tmp = t_1; elseif (z <= -4.3e-223) tmp = y * 4.0; elseif (z <= -9.8e-281) tmp = x * -3.0; elseif (z <= 1.1e-268) tmp = y * 4.0; elseif (z <= 4.3e-159) tmp = x * -3.0; elseif (z <= 0.086) tmp = y * 4.0; elseif (z <= 6e+194) tmp = t_1; elseif (z <= 3.6e+264) tmp = t_0; else tmp = -6.0 * (y * z); end tmp_2 = tmp; end
code[x_, y_, z_] := Block[{t$95$0 = N[(6.0 * N[(x * z), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(y * N[(z * -6.0), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[z, -8.5e+64], t$95$0, If[LessEqual[z, -0.66], t$95$1, If[LessEqual[z, -4.3e-223], N[(y * 4.0), $MachinePrecision], If[LessEqual[z, -9.8e-281], N[(x * -3.0), $MachinePrecision], If[LessEqual[z, 1.1e-268], N[(y * 4.0), $MachinePrecision], If[LessEqual[z, 4.3e-159], N[(x * -3.0), $MachinePrecision], If[LessEqual[z, 0.086], N[(y * 4.0), $MachinePrecision], If[LessEqual[z, 6e+194], t$95$1, If[LessEqual[z, 3.6e+264], t$95$0, N[(-6.0 * N[(y * z), $MachinePrecision]), $MachinePrecision]]]]]]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := 6 \cdot \left(x \cdot z\right)\\
t_1 := y \cdot \left(z \cdot -6\right)\\
\mathbf{if}\;z \leq -8.5 \cdot 10^{+64}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;z \leq -0.66:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;z \leq -4.3 \cdot 10^{-223}:\\
\;\;\;\;y \cdot 4\\
\mathbf{elif}\;z \leq -9.8 \cdot 10^{-281}:\\
\;\;\;\;x \cdot -3\\
\mathbf{elif}\;z \leq 1.1 \cdot 10^{-268}:\\
\;\;\;\;y \cdot 4\\
\mathbf{elif}\;z \leq 4.3 \cdot 10^{-159}:\\
\;\;\;\;x \cdot -3\\
\mathbf{elif}\;z \leq 0.086:\\
\;\;\;\;y \cdot 4\\
\mathbf{elif}\;z \leq 6 \cdot 10^{+194}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;z \leq 3.6 \cdot 10^{+264}:\\
\;\;\;\;t\_0\\
\mathbf{else}:\\
\;\;\;\;-6 \cdot \left(y \cdot z\right)\\
\end{array}
\end{array}
if z < -8.4999999999999998e64 or 6.0000000000000006e194 < z < 3.60000000000000011e264Initial program 99.8%
metadata-eval99.8%
Simplified99.8%
Taylor expanded in x around inf 71.2%
sub-neg71.2%
distribute-rgt-in71.2%
metadata-eval71.2%
metadata-eval71.2%
distribute-lft-neg-in71.2%
associate-+r+71.2%
metadata-eval71.2%
metadata-eval71.2%
distribute-rgt-neg-in71.2%
metadata-eval71.2%
Simplified71.2%
Taylor expanded in z around inf 71.2%
if -8.4999999999999998e64 < z < -0.660000000000000031 or 0.085999999999999993 < z < 6.0000000000000006e194Initial program 99.7%
+-commutative99.7%
associate-*l*99.8%
fma-def99.8%
sub-neg99.8%
distribute-rgt-in99.7%
metadata-eval99.7%
metadata-eval99.7%
distribute-lft-neg-out99.7%
distribute-rgt-neg-in99.7%
metadata-eval99.7%
Simplified99.7%
Taylor expanded in y around inf 67.7%
Taylor expanded in z around inf 66.7%
if -0.660000000000000031 < z < -4.2999999999999999e-223 or -9.7999999999999999e-281 < z < 1.10000000000000002e-268 or 4.3e-159 < z < 0.085999999999999993Initial program 99.4%
+-commutative99.4%
associate-*l*99.8%
fma-def99.9%
sub-neg99.9%
distribute-rgt-in99.9%
metadata-eval99.9%
metadata-eval99.9%
distribute-lft-neg-out99.9%
distribute-rgt-neg-in99.9%
metadata-eval99.9%
Simplified99.9%
Taylor expanded in y around inf 61.0%
Taylor expanded in z around 0 59.0%
*-commutative59.0%
Simplified59.0%
if -4.2999999999999999e-223 < z < -9.7999999999999999e-281 or 1.10000000000000002e-268 < z < 4.3e-159Initial program 99.3%
metadata-eval99.3%
Simplified99.3%
Taylor expanded in x around inf 70.0%
sub-neg70.0%
distribute-rgt-in70.0%
metadata-eval70.0%
metadata-eval70.0%
distribute-lft-neg-in70.0%
associate-+r+70.0%
metadata-eval70.0%
metadata-eval70.0%
distribute-rgt-neg-in70.0%
metadata-eval70.0%
Simplified70.0%
Taylor expanded in z around 0 70.0%
*-commutative70.0%
Simplified70.0%
if 3.60000000000000011e264 < z Initial program 99.7%
+-commutative99.7%
associate-*l*100.0%
fma-def100.0%
sub-neg100.0%
distribute-rgt-in100.0%
metadata-eval100.0%
metadata-eval100.0%
distribute-lft-neg-out100.0%
distribute-rgt-neg-in100.0%
metadata-eval100.0%
Simplified100.0%
Taylor expanded in y around inf 83.3%
Taylor expanded in z around inf 83.3%
*-commutative83.3%
Simplified83.3%
Final simplification65.5%
(FPCore (x y z)
:precision binary64
(let* ((t_0 (* y (* z -6.0))))
(if (<= z -8.3e+64)
(* z (* x 6.0))
(if (<= z -0.66)
t_0
(if (<= z -7.2e-222)
(* y 4.0)
(if (<= z -2.2e-279)
(* x -3.0)
(if (<= z 1.4e-269)
(* y 4.0)
(if (<= z 4.3e-159)
(* x -3.0)
(if (<= z 0.086)
(* y 4.0)
(if (<= z 2.7e+202)
t_0
(if (<= z 2.6e+266)
(* 6.0 (* x z))
(* -6.0 (* y z)))))))))))))
double code(double x, double y, double z) {
double t_0 = y * (z * -6.0);
double tmp;
if (z <= -8.3e+64) {
tmp = z * (x * 6.0);
} else if (z <= -0.66) {
tmp = t_0;
} else if (z <= -7.2e-222) {
tmp = y * 4.0;
} else if (z <= -2.2e-279) {
tmp = x * -3.0;
} else if (z <= 1.4e-269) {
tmp = y * 4.0;
} else if (z <= 4.3e-159) {
tmp = x * -3.0;
} else if (z <= 0.086) {
tmp = y * 4.0;
} else if (z <= 2.7e+202) {
tmp = t_0;
} else if (z <= 2.6e+266) {
tmp = 6.0 * (x * z);
} else {
tmp = -6.0 * (y * 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) :: t_0
real(8) :: tmp
t_0 = y * (z * (-6.0d0))
if (z <= (-8.3d+64)) then
tmp = z * (x * 6.0d0)
else if (z <= (-0.66d0)) then
tmp = t_0
else if (z <= (-7.2d-222)) then
tmp = y * 4.0d0
else if (z <= (-2.2d-279)) then
tmp = x * (-3.0d0)
else if (z <= 1.4d-269) then
tmp = y * 4.0d0
else if (z <= 4.3d-159) then
tmp = x * (-3.0d0)
else if (z <= 0.086d0) then
tmp = y * 4.0d0
else if (z <= 2.7d+202) then
tmp = t_0
else if (z <= 2.6d+266) then
tmp = 6.0d0 * (x * z)
else
tmp = (-6.0d0) * (y * z)
end if
code = tmp
end function
public static double code(double x, double y, double z) {
double t_0 = y * (z * -6.0);
double tmp;
if (z <= -8.3e+64) {
tmp = z * (x * 6.0);
} else if (z <= -0.66) {
tmp = t_0;
} else if (z <= -7.2e-222) {
tmp = y * 4.0;
} else if (z <= -2.2e-279) {
tmp = x * -3.0;
} else if (z <= 1.4e-269) {
tmp = y * 4.0;
} else if (z <= 4.3e-159) {
tmp = x * -3.0;
} else if (z <= 0.086) {
tmp = y * 4.0;
} else if (z <= 2.7e+202) {
tmp = t_0;
} else if (z <= 2.6e+266) {
tmp = 6.0 * (x * z);
} else {
tmp = -6.0 * (y * z);
}
return tmp;
}
def code(x, y, z): t_0 = y * (z * -6.0) tmp = 0 if z <= -8.3e+64: tmp = z * (x * 6.0) elif z <= -0.66: tmp = t_0 elif z <= -7.2e-222: tmp = y * 4.0 elif z <= -2.2e-279: tmp = x * -3.0 elif z <= 1.4e-269: tmp = y * 4.0 elif z <= 4.3e-159: tmp = x * -3.0 elif z <= 0.086: tmp = y * 4.0 elif z <= 2.7e+202: tmp = t_0 elif z <= 2.6e+266: tmp = 6.0 * (x * z) else: tmp = -6.0 * (y * z) return tmp
function code(x, y, z) t_0 = Float64(y * Float64(z * -6.0)) tmp = 0.0 if (z <= -8.3e+64) tmp = Float64(z * Float64(x * 6.0)); elseif (z <= -0.66) tmp = t_0; elseif (z <= -7.2e-222) tmp = Float64(y * 4.0); elseif (z <= -2.2e-279) tmp = Float64(x * -3.0); elseif (z <= 1.4e-269) tmp = Float64(y * 4.0); elseif (z <= 4.3e-159) tmp = Float64(x * -3.0); elseif (z <= 0.086) tmp = Float64(y * 4.0); elseif (z <= 2.7e+202) tmp = t_0; elseif (z <= 2.6e+266) tmp = Float64(6.0 * Float64(x * z)); else tmp = Float64(-6.0 * Float64(y * z)); end return tmp end
function tmp_2 = code(x, y, z) t_0 = y * (z * -6.0); tmp = 0.0; if (z <= -8.3e+64) tmp = z * (x * 6.0); elseif (z <= -0.66) tmp = t_0; elseif (z <= -7.2e-222) tmp = y * 4.0; elseif (z <= -2.2e-279) tmp = x * -3.0; elseif (z <= 1.4e-269) tmp = y * 4.0; elseif (z <= 4.3e-159) tmp = x * -3.0; elseif (z <= 0.086) tmp = y * 4.0; elseif (z <= 2.7e+202) tmp = t_0; elseif (z <= 2.6e+266) tmp = 6.0 * (x * z); else tmp = -6.0 * (y * z); end tmp_2 = tmp; end
code[x_, y_, z_] := Block[{t$95$0 = N[(y * N[(z * -6.0), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[z, -8.3e+64], N[(z * N[(x * 6.0), $MachinePrecision]), $MachinePrecision], If[LessEqual[z, -0.66], t$95$0, If[LessEqual[z, -7.2e-222], N[(y * 4.0), $MachinePrecision], If[LessEqual[z, -2.2e-279], N[(x * -3.0), $MachinePrecision], If[LessEqual[z, 1.4e-269], N[(y * 4.0), $MachinePrecision], If[LessEqual[z, 4.3e-159], N[(x * -3.0), $MachinePrecision], If[LessEqual[z, 0.086], N[(y * 4.0), $MachinePrecision], If[LessEqual[z, 2.7e+202], t$95$0, If[LessEqual[z, 2.6e+266], N[(6.0 * N[(x * z), $MachinePrecision]), $MachinePrecision], N[(-6.0 * N[(y * z), $MachinePrecision]), $MachinePrecision]]]]]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := y \cdot \left(z \cdot -6\right)\\
\mathbf{if}\;z \leq -8.3 \cdot 10^{+64}:\\
\;\;\;\;z \cdot \left(x \cdot 6\right)\\
\mathbf{elif}\;z \leq -0.66:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;z \leq -7.2 \cdot 10^{-222}:\\
\;\;\;\;y \cdot 4\\
\mathbf{elif}\;z \leq -2.2 \cdot 10^{-279}:\\
\;\;\;\;x \cdot -3\\
\mathbf{elif}\;z \leq 1.4 \cdot 10^{-269}:\\
\;\;\;\;y \cdot 4\\
\mathbf{elif}\;z \leq 4.3 \cdot 10^{-159}:\\
\;\;\;\;x \cdot -3\\
\mathbf{elif}\;z \leq 0.086:\\
\;\;\;\;y \cdot 4\\
\mathbf{elif}\;z \leq 2.7 \cdot 10^{+202}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;z \leq 2.6 \cdot 10^{+266}:\\
\;\;\;\;6 \cdot \left(x \cdot z\right)\\
\mathbf{else}:\\
\;\;\;\;-6 \cdot \left(y \cdot z\right)\\
\end{array}
\end{array}
if z < -8.3000000000000004e64Initial program 99.8%
metadata-eval99.8%
Simplified99.8%
Taylor expanded in x around inf 67.7%
sub-neg67.7%
distribute-rgt-in67.7%
metadata-eval67.7%
metadata-eval67.7%
distribute-lft-neg-in67.7%
associate-+r+67.7%
metadata-eval67.7%
metadata-eval67.7%
distribute-rgt-neg-in67.7%
metadata-eval67.7%
Simplified67.7%
Taylor expanded in z around inf 67.7%
*-commutative67.7%
associate-*l*67.7%
*-commutative67.7%
associate-*r*67.7%
Simplified67.7%
if -8.3000000000000004e64 < z < -0.660000000000000031 or 0.085999999999999993 < z < 2.69999999999999995e202Initial program 99.7%
+-commutative99.7%
associate-*l*99.8%
fma-def99.8%
sub-neg99.8%
distribute-rgt-in99.7%
metadata-eval99.7%
metadata-eval99.7%
distribute-lft-neg-out99.7%
distribute-rgt-neg-in99.7%
metadata-eval99.7%
Simplified99.7%
Taylor expanded in y around inf 67.7%
Taylor expanded in z around inf 66.7%
if -0.660000000000000031 < z < -7.19999999999999948e-222 or -2.2e-279 < z < 1.39999999999999997e-269 or 4.3e-159 < z < 0.085999999999999993Initial program 99.4%
+-commutative99.4%
associate-*l*99.8%
fma-def99.9%
sub-neg99.9%
distribute-rgt-in99.9%
metadata-eval99.9%
metadata-eval99.9%
distribute-lft-neg-out99.9%
distribute-rgt-neg-in99.9%
metadata-eval99.9%
Simplified99.9%
Taylor expanded in y around inf 61.0%
Taylor expanded in z around 0 59.0%
*-commutative59.0%
Simplified59.0%
if -7.19999999999999948e-222 < z < -2.2e-279 or 1.39999999999999997e-269 < z < 4.3e-159Initial program 99.3%
metadata-eval99.3%
Simplified99.3%
Taylor expanded in x around inf 70.0%
sub-neg70.0%
distribute-rgt-in70.0%
metadata-eval70.0%
metadata-eval70.0%
distribute-lft-neg-in70.0%
associate-+r+70.0%
metadata-eval70.0%
metadata-eval70.0%
distribute-rgt-neg-in70.0%
metadata-eval70.0%
Simplified70.0%
Taylor expanded in z around 0 70.0%
*-commutative70.0%
Simplified70.0%
if 2.69999999999999995e202 < z < 2.60000000000000014e266Initial program 100.0%
metadata-eval100.0%
Simplified100.0%
Taylor expanded in x around inf 89.5%
sub-neg89.5%
distribute-rgt-in89.5%
metadata-eval89.5%
metadata-eval89.5%
distribute-lft-neg-in89.5%
associate-+r+89.5%
metadata-eval89.5%
metadata-eval89.5%
distribute-rgt-neg-in89.5%
metadata-eval89.5%
Simplified89.5%
Taylor expanded in z around inf 89.6%
if 2.60000000000000014e266 < z Initial program 99.7%
+-commutative99.7%
associate-*l*100.0%
fma-def100.0%
sub-neg100.0%
distribute-rgt-in100.0%
metadata-eval100.0%
metadata-eval100.0%
distribute-lft-neg-out100.0%
distribute-rgt-neg-in100.0%
metadata-eval100.0%
Simplified100.0%
Taylor expanded in y around inf 83.3%
Taylor expanded in z around inf 83.3%
*-commutative83.3%
Simplified83.3%
Final simplification65.5%
(FPCore (x y z)
:precision binary64
(let* ((t_0 (* -6.0 (* y z))))
(if (<= z -0.66)
t_0
(if (<= z -5.3e-222)
(* y 4.0)
(if (<= z -2.5e-277)
(* x -3.0)
(if (<= z 1.45e-272)
(* y 4.0)
(if (<= z 3.5e-156)
(* x -3.0)
(if (<= z 0.086) (* y 4.0) t_0))))))))
double code(double x, double y, double z) {
double t_0 = -6.0 * (y * z);
double tmp;
if (z <= -0.66) {
tmp = t_0;
} else if (z <= -5.3e-222) {
tmp = y * 4.0;
} else if (z <= -2.5e-277) {
tmp = x * -3.0;
} else if (z <= 1.45e-272) {
tmp = y * 4.0;
} else if (z <= 3.5e-156) {
tmp = x * -3.0;
} else if (z <= 0.086) {
tmp = y * 4.0;
} 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 = (-6.0d0) * (y * z)
if (z <= (-0.66d0)) then
tmp = t_0
else if (z <= (-5.3d-222)) then
tmp = y * 4.0d0
else if (z <= (-2.5d-277)) then
tmp = x * (-3.0d0)
else if (z <= 1.45d-272) then
tmp = y * 4.0d0
else if (z <= 3.5d-156) then
tmp = x * (-3.0d0)
else if (z <= 0.086d0) then
tmp = y * 4.0d0
else
tmp = t_0
end if
code = tmp
end function
public static double code(double x, double y, double z) {
double t_0 = -6.0 * (y * z);
double tmp;
if (z <= -0.66) {
tmp = t_0;
} else if (z <= -5.3e-222) {
tmp = y * 4.0;
} else if (z <= -2.5e-277) {
tmp = x * -3.0;
} else if (z <= 1.45e-272) {
tmp = y * 4.0;
} else if (z <= 3.5e-156) {
tmp = x * -3.0;
} else if (z <= 0.086) {
tmp = y * 4.0;
} else {
tmp = t_0;
}
return tmp;
}
def code(x, y, z): t_0 = -6.0 * (y * z) tmp = 0 if z <= -0.66: tmp = t_0 elif z <= -5.3e-222: tmp = y * 4.0 elif z <= -2.5e-277: tmp = x * -3.0 elif z <= 1.45e-272: tmp = y * 4.0 elif z <= 3.5e-156: tmp = x * -3.0 elif z <= 0.086: tmp = y * 4.0 else: tmp = t_0 return tmp
function code(x, y, z) t_0 = Float64(-6.0 * Float64(y * z)) tmp = 0.0 if (z <= -0.66) tmp = t_0; elseif (z <= -5.3e-222) tmp = Float64(y * 4.0); elseif (z <= -2.5e-277) tmp = Float64(x * -3.0); elseif (z <= 1.45e-272) tmp = Float64(y * 4.0); elseif (z <= 3.5e-156) tmp = Float64(x * -3.0); elseif (z <= 0.086) tmp = Float64(y * 4.0); else tmp = t_0; end return tmp end
function tmp_2 = code(x, y, z) t_0 = -6.0 * (y * z); tmp = 0.0; if (z <= -0.66) tmp = t_0; elseif (z <= -5.3e-222) tmp = y * 4.0; elseif (z <= -2.5e-277) tmp = x * -3.0; elseif (z <= 1.45e-272) tmp = y * 4.0; elseif (z <= 3.5e-156) tmp = x * -3.0; elseif (z <= 0.086) tmp = y * 4.0; else tmp = t_0; end tmp_2 = tmp; end
code[x_, y_, z_] := Block[{t$95$0 = N[(-6.0 * N[(y * z), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[z, -0.66], t$95$0, If[LessEqual[z, -5.3e-222], N[(y * 4.0), $MachinePrecision], If[LessEqual[z, -2.5e-277], N[(x * -3.0), $MachinePrecision], If[LessEqual[z, 1.45e-272], N[(y * 4.0), $MachinePrecision], If[LessEqual[z, 3.5e-156], N[(x * -3.0), $MachinePrecision], If[LessEqual[z, 0.086], N[(y * 4.0), $MachinePrecision], t$95$0]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := -6 \cdot \left(y \cdot z\right)\\
\mathbf{if}\;z \leq -0.66:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;z \leq -5.3 \cdot 10^{-222}:\\
\;\;\;\;y \cdot 4\\
\mathbf{elif}\;z \leq -2.5 \cdot 10^{-277}:\\
\;\;\;\;x \cdot -3\\
\mathbf{elif}\;z \leq 1.45 \cdot 10^{-272}:\\
\;\;\;\;y \cdot 4\\
\mathbf{elif}\;z \leq 3.5 \cdot 10^{-156}:\\
\;\;\;\;x \cdot -3\\
\mathbf{elif}\;z \leq 0.086:\\
\;\;\;\;y \cdot 4\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
if z < -0.660000000000000031 or 0.085999999999999993 < z Initial program 99.8%
+-commutative99.8%
associate-*l*99.8%
fma-def99.8%
sub-neg99.8%
distribute-rgt-in99.7%
metadata-eval99.7%
metadata-eval99.7%
distribute-lft-neg-out99.7%
distribute-rgt-neg-in99.7%
metadata-eval99.7%
Simplified99.7%
Taylor expanded in y around inf 52.3%
Taylor expanded in z around inf 51.8%
*-commutative51.8%
Simplified51.8%
if -0.660000000000000031 < z < -5.29999999999999981e-222 or -2.5e-277 < z < 1.44999999999999997e-272 or 3.4999999999999999e-156 < z < 0.085999999999999993Initial program 99.4%
+-commutative99.4%
associate-*l*99.8%
fma-def99.9%
sub-neg99.9%
distribute-rgt-in99.9%
metadata-eval99.9%
metadata-eval99.9%
distribute-lft-neg-out99.9%
distribute-rgt-neg-in99.9%
metadata-eval99.9%
Simplified99.9%
Taylor expanded in y around inf 61.0%
Taylor expanded in z around 0 59.0%
*-commutative59.0%
Simplified59.0%
if -5.29999999999999981e-222 < z < -2.5e-277 or 1.44999999999999997e-272 < z < 3.4999999999999999e-156Initial program 99.3%
metadata-eval99.3%
Simplified99.3%
Taylor expanded in x around inf 70.0%
sub-neg70.0%
distribute-rgt-in70.0%
metadata-eval70.0%
metadata-eval70.0%
distribute-lft-neg-in70.0%
associate-+r+70.0%
metadata-eval70.0%
metadata-eval70.0%
distribute-rgt-neg-in70.0%
metadata-eval70.0%
Simplified70.0%
Taylor expanded in z around 0 70.0%
*-commutative70.0%
Simplified70.0%
Final simplification57.2%
(FPCore (x y z)
:precision binary64
(let* ((t_0 (* z (* x 6.0))))
(if (<= x -1.25e+229)
(* x -3.0)
(if (<= x -2.7e+93)
t_0
(if (<= x 1.55e+66)
(* 6.0 (* y (- 0.6666666666666666 z)))
(if (<= x 7.5e+213) t_0 (* x -3.0)))))))
double code(double x, double y, double z) {
double t_0 = z * (x * 6.0);
double tmp;
if (x <= -1.25e+229) {
tmp = x * -3.0;
} else if (x <= -2.7e+93) {
tmp = t_0;
} else if (x <= 1.55e+66) {
tmp = 6.0 * (y * (0.6666666666666666 - z));
} else if (x <= 7.5e+213) {
tmp = t_0;
} else {
tmp = x * -3.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 * (x * 6.0d0)
if (x <= (-1.25d+229)) then
tmp = x * (-3.0d0)
else if (x <= (-2.7d+93)) then
tmp = t_0
else if (x <= 1.55d+66) then
tmp = 6.0d0 * (y * (0.6666666666666666d0 - z))
else if (x <= 7.5d+213) then
tmp = t_0
else
tmp = x * (-3.0d0)
end if
code = tmp
end function
public static double code(double x, double y, double z) {
double t_0 = z * (x * 6.0);
double tmp;
if (x <= -1.25e+229) {
tmp = x * -3.0;
} else if (x <= -2.7e+93) {
tmp = t_0;
} else if (x <= 1.55e+66) {
tmp = 6.0 * (y * (0.6666666666666666 - z));
} else if (x <= 7.5e+213) {
tmp = t_0;
} else {
tmp = x * -3.0;
}
return tmp;
}
def code(x, y, z): t_0 = z * (x * 6.0) tmp = 0 if x <= -1.25e+229: tmp = x * -3.0 elif x <= -2.7e+93: tmp = t_0 elif x <= 1.55e+66: tmp = 6.0 * (y * (0.6666666666666666 - z)) elif x <= 7.5e+213: tmp = t_0 else: tmp = x * -3.0 return tmp
function code(x, y, z) t_0 = Float64(z * Float64(x * 6.0)) tmp = 0.0 if (x <= -1.25e+229) tmp = Float64(x * -3.0); elseif (x <= -2.7e+93) tmp = t_0; elseif (x <= 1.55e+66) tmp = Float64(6.0 * Float64(y * Float64(0.6666666666666666 - z))); elseif (x <= 7.5e+213) tmp = t_0; else tmp = Float64(x * -3.0); end return tmp end
function tmp_2 = code(x, y, z) t_0 = z * (x * 6.0); tmp = 0.0; if (x <= -1.25e+229) tmp = x * -3.0; elseif (x <= -2.7e+93) tmp = t_0; elseif (x <= 1.55e+66) tmp = 6.0 * (y * (0.6666666666666666 - z)); elseif (x <= 7.5e+213) tmp = t_0; else tmp = x * -3.0; end tmp_2 = tmp; end
code[x_, y_, z_] := Block[{t$95$0 = N[(z * N[(x * 6.0), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[x, -1.25e+229], N[(x * -3.0), $MachinePrecision], If[LessEqual[x, -2.7e+93], t$95$0, If[LessEqual[x, 1.55e+66], N[(6.0 * N[(y * N[(0.6666666666666666 - z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[x, 7.5e+213], t$95$0, N[(x * -3.0), $MachinePrecision]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := z \cdot \left(x \cdot 6\right)\\
\mathbf{if}\;x \leq -1.25 \cdot 10^{+229}:\\
\;\;\;\;x \cdot -3\\
\mathbf{elif}\;x \leq -2.7 \cdot 10^{+93}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;x \leq 1.55 \cdot 10^{+66}:\\
\;\;\;\;6 \cdot \left(y \cdot \left(0.6666666666666666 - z\right)\right)\\
\mathbf{elif}\;x \leq 7.5 \cdot 10^{+213}:\\
\;\;\;\;t\_0\\
\mathbf{else}:\\
\;\;\;\;x \cdot -3\\
\end{array}
\end{array}
if x < -1.25000000000000012e229 or 7.5000000000000003e213 < x Initial program 99.7%
metadata-eval99.7%
Simplified99.7%
Taylor expanded in x around inf 91.6%
sub-neg91.6%
distribute-rgt-in91.6%
metadata-eval91.6%
metadata-eval91.6%
distribute-lft-neg-in91.6%
associate-+r+91.6%
metadata-eval91.6%
metadata-eval91.6%
distribute-rgt-neg-in91.6%
metadata-eval91.6%
Simplified91.6%
Taylor expanded in z around 0 63.0%
*-commutative63.0%
Simplified63.0%
if -1.25000000000000012e229 < x < -2.6999999999999999e93 or 1.55000000000000009e66 < x < 7.5000000000000003e213Initial program 99.6%
metadata-eval99.6%
Simplified99.6%
Taylor expanded in x around inf 83.9%
sub-neg83.9%
distribute-rgt-in83.9%
metadata-eval83.9%
metadata-eval83.9%
distribute-lft-neg-in83.9%
associate-+r+83.9%
metadata-eval83.9%
metadata-eval83.9%
distribute-rgt-neg-in83.9%
metadata-eval83.9%
Simplified83.9%
Taylor expanded in z around inf 52.9%
*-commutative52.9%
associate-*l*52.9%
*-commutative52.9%
associate-*r*52.9%
Simplified52.9%
if -2.6999999999999999e93 < x < 1.55000000000000009e66Initial program 99.5%
metadata-eval99.5%
Simplified99.5%
Taylor expanded in x around 0 99.5%
Taylor expanded in y around inf 75.4%
*-commutative75.4%
Simplified75.4%
Final simplification68.0%
(FPCore (x y z) :precision binary64 (if (or (<= z -0.58) (not (<= z 0.53))) (* z (+ (* y -6.0) (* x 6.0))) (+ x (* (- y x) 4.0))))
double code(double x, double y, double z) {
double tmp;
if ((z <= -0.58) || !(z <= 0.53)) {
tmp = z * ((y * -6.0) + (x * 6.0));
} else {
tmp = x + ((y - x) * 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 ((z <= (-0.58d0)) .or. (.not. (z <= 0.53d0))) then
tmp = z * ((y * (-6.0d0)) + (x * 6.0d0))
else
tmp = x + ((y - x) * 4.0d0)
end if
code = tmp
end function
public static double code(double x, double y, double z) {
double tmp;
if ((z <= -0.58) || !(z <= 0.53)) {
tmp = z * ((y * -6.0) + (x * 6.0));
} else {
tmp = x + ((y - x) * 4.0);
}
return tmp;
}
def code(x, y, z): tmp = 0 if (z <= -0.58) or not (z <= 0.53): tmp = z * ((y * -6.0) + (x * 6.0)) else: tmp = x + ((y - x) * 4.0) return tmp
function code(x, y, z) tmp = 0.0 if ((z <= -0.58) || !(z <= 0.53)) tmp = Float64(z * Float64(Float64(y * -6.0) + Float64(x * 6.0))); else tmp = Float64(x + Float64(Float64(y - x) * 4.0)); end return tmp end
function tmp_2 = code(x, y, z) tmp = 0.0; if ((z <= -0.58) || ~((z <= 0.53))) tmp = z * ((y * -6.0) + (x * 6.0)); else tmp = x + ((y - x) * 4.0); end tmp_2 = tmp; end
code[x_, y_, z_] := If[Or[LessEqual[z, -0.58], N[Not[LessEqual[z, 0.53]], $MachinePrecision]], N[(z * N[(N[(y * -6.0), $MachinePrecision] + N[(x * 6.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x + N[(N[(y - x), $MachinePrecision] * 4.0), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;z \leq -0.58 \lor \neg \left(z \leq 0.53\right):\\
\;\;\;\;z \cdot \left(y \cdot -6 + x \cdot 6\right)\\
\mathbf{else}:\\
\;\;\;\;x + \left(y - x\right) \cdot 4\\
\end{array}
\end{array}
if z < -0.57999999999999996 or 0.53000000000000003 < z Initial program 99.8%
metadata-eval99.8%
Simplified99.8%
Taylor expanded in x around 0 96.3%
Taylor expanded in z around inf 99.2%
if -0.57999999999999996 < z < 0.53000000000000003Initial program 99.4%
metadata-eval99.4%
Simplified99.4%
Taylor expanded in z around 0 96.4%
Final simplification97.7%
(FPCore (x y z) :precision binary64 (if (<= z -0.58) (* z (+ (* y -6.0) (* x 6.0))) (if (<= z 0.55) (+ x (* (- y x) 4.0)) (+ x (* -6.0 (* (- y x) z))))))
double code(double x, double y, double z) {
double tmp;
if (z <= -0.58) {
tmp = z * ((y * -6.0) + (x * 6.0));
} else if (z <= 0.55) {
tmp = x + ((y - x) * 4.0);
} else {
tmp = x + (-6.0 * ((y - 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 (z <= (-0.58d0)) then
tmp = z * ((y * (-6.0d0)) + (x * 6.0d0))
else if (z <= 0.55d0) then
tmp = x + ((y - x) * 4.0d0)
else
tmp = x + ((-6.0d0) * ((y - x) * z))
end if
code = tmp
end function
public static double code(double x, double y, double z) {
double tmp;
if (z <= -0.58) {
tmp = z * ((y * -6.0) + (x * 6.0));
} else if (z <= 0.55) {
tmp = x + ((y - x) * 4.0);
} else {
tmp = x + (-6.0 * ((y - x) * z));
}
return tmp;
}
def code(x, y, z): tmp = 0 if z <= -0.58: tmp = z * ((y * -6.0) + (x * 6.0)) elif z <= 0.55: tmp = x + ((y - x) * 4.0) else: tmp = x + (-6.0 * ((y - x) * z)) return tmp
function code(x, y, z) tmp = 0.0 if (z <= -0.58) tmp = Float64(z * Float64(Float64(y * -6.0) + Float64(x * 6.0))); elseif (z <= 0.55) tmp = Float64(x + Float64(Float64(y - x) * 4.0)); else tmp = Float64(x + Float64(-6.0 * Float64(Float64(y - x) * z))); end return tmp end
function tmp_2 = code(x, y, z) tmp = 0.0; if (z <= -0.58) tmp = z * ((y * -6.0) + (x * 6.0)); elseif (z <= 0.55) tmp = x + ((y - x) * 4.0); else tmp = x + (-6.0 * ((y - x) * z)); end tmp_2 = tmp; end
code[x_, y_, z_] := If[LessEqual[z, -0.58], N[(z * N[(N[(y * -6.0), $MachinePrecision] + N[(x * 6.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[z, 0.55], N[(x + N[(N[(y - x), $MachinePrecision] * 4.0), $MachinePrecision]), $MachinePrecision], N[(x + N[(-6.0 * N[(N[(y - x), $MachinePrecision] * z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;z \leq -0.58:\\
\;\;\;\;z \cdot \left(y \cdot -6 + x \cdot 6\right)\\
\mathbf{elif}\;z \leq 0.55:\\
\;\;\;\;x + \left(y - x\right) \cdot 4\\
\mathbf{else}:\\
\;\;\;\;x + -6 \cdot \left(\left(y - x\right) \cdot z\right)\\
\end{array}
\end{array}
if z < -0.57999999999999996Initial program 99.8%
metadata-eval99.8%
Simplified99.8%
Taylor expanded in x around 0 96.6%
Taylor expanded in z around inf 98.7%
if -0.57999999999999996 < z < 0.55000000000000004Initial program 99.4%
metadata-eval99.4%
Simplified99.4%
Taylor expanded in z around 0 96.4%
if 0.55000000000000004 < z Initial program 99.8%
metadata-eval99.8%
Simplified99.8%
Taylor expanded in z around inf 99.9%
Final simplification97.7%
(FPCore (x y z) :precision binary64 (if (or (<= x -410000.0) (not (<= x 2.7e+42))) (* x (+ -3.0 (* z 6.0))) (* 6.0 (* y (- 0.6666666666666666 z)))))
double code(double x, double y, double z) {
double tmp;
if ((x <= -410000.0) || !(x <= 2.7e+42)) {
tmp = x * (-3.0 + (z * 6.0));
} else {
tmp = 6.0 * (y * (0.6666666666666666 - 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 ((x <= (-410000.0d0)) .or. (.not. (x <= 2.7d+42))) then
tmp = x * ((-3.0d0) + (z * 6.0d0))
else
tmp = 6.0d0 * (y * (0.6666666666666666d0 - z))
end if
code = tmp
end function
public static double code(double x, double y, double z) {
double tmp;
if ((x <= -410000.0) || !(x <= 2.7e+42)) {
tmp = x * (-3.0 + (z * 6.0));
} else {
tmp = 6.0 * (y * (0.6666666666666666 - z));
}
return tmp;
}
def code(x, y, z): tmp = 0 if (x <= -410000.0) or not (x <= 2.7e+42): tmp = x * (-3.0 + (z * 6.0)) else: tmp = 6.0 * (y * (0.6666666666666666 - z)) return tmp
function code(x, y, z) tmp = 0.0 if ((x <= -410000.0) || !(x <= 2.7e+42)) tmp = Float64(x * Float64(-3.0 + Float64(z * 6.0))); else tmp = Float64(6.0 * Float64(y * Float64(0.6666666666666666 - z))); end return tmp end
function tmp_2 = code(x, y, z) tmp = 0.0; if ((x <= -410000.0) || ~((x <= 2.7e+42))) tmp = x * (-3.0 + (z * 6.0)); else tmp = 6.0 * (y * (0.6666666666666666 - z)); end tmp_2 = tmp; end
code[x_, y_, z_] := If[Or[LessEqual[x, -410000.0], N[Not[LessEqual[x, 2.7e+42]], $MachinePrecision]], N[(x * N[(-3.0 + N[(z * 6.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(6.0 * N[(y * N[(0.6666666666666666 - z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq -410000 \lor \neg \left(x \leq 2.7 \cdot 10^{+42}\right):\\
\;\;\;\;x \cdot \left(-3 + z \cdot 6\right)\\
\mathbf{else}:\\
\;\;\;\;6 \cdot \left(y \cdot \left(0.6666666666666666 - z\right)\right)\\
\end{array}
\end{array}
if x < -4.1e5 or 2.7000000000000001e42 < x Initial program 99.6%
metadata-eval99.6%
Simplified99.6%
Taylor expanded in x around inf 81.7%
sub-neg81.7%
distribute-rgt-in81.7%
metadata-eval81.7%
metadata-eval81.7%
distribute-lft-neg-in81.7%
associate-+r+81.7%
metadata-eval81.7%
metadata-eval81.7%
distribute-rgt-neg-in81.7%
metadata-eval81.7%
Simplified81.7%
if -4.1e5 < x < 2.7000000000000001e42Initial program 99.6%
metadata-eval99.6%
Simplified99.6%
Taylor expanded in x around 0 99.5%
Taylor expanded in y around inf 81.9%
*-commutative81.9%
Simplified81.9%
Final simplification81.8%
(FPCore (x y z) :precision binary64 (if (or (<= x -59000.0) (not (<= x 3.5e+42))) (* x (+ -3.0 (* z 6.0))) (* y (+ 4.0 (* z -6.0)))))
double code(double x, double y, double z) {
double tmp;
if ((x <= -59000.0) || !(x <= 3.5e+42)) {
tmp = x * (-3.0 + (z * 6.0));
} else {
tmp = y * (4.0 + (z * -6.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 ((x <= (-59000.0d0)) .or. (.not. (x <= 3.5d+42))) then
tmp = x * ((-3.0d0) + (z * 6.0d0))
else
tmp = y * (4.0d0 + (z * (-6.0d0)))
end if
code = tmp
end function
public static double code(double x, double y, double z) {
double tmp;
if ((x <= -59000.0) || !(x <= 3.5e+42)) {
tmp = x * (-3.0 + (z * 6.0));
} else {
tmp = y * (4.0 + (z * -6.0));
}
return tmp;
}
def code(x, y, z): tmp = 0 if (x <= -59000.0) or not (x <= 3.5e+42): tmp = x * (-3.0 + (z * 6.0)) else: tmp = y * (4.0 + (z * -6.0)) return tmp
function code(x, y, z) tmp = 0.0 if ((x <= -59000.0) || !(x <= 3.5e+42)) tmp = Float64(x * Float64(-3.0 + Float64(z * 6.0))); else tmp = Float64(y * Float64(4.0 + Float64(z * -6.0))); end return tmp end
function tmp_2 = code(x, y, z) tmp = 0.0; if ((x <= -59000.0) || ~((x <= 3.5e+42))) tmp = x * (-3.0 + (z * 6.0)); else tmp = y * (4.0 + (z * -6.0)); end tmp_2 = tmp; end
code[x_, y_, z_] := If[Or[LessEqual[x, -59000.0], N[Not[LessEqual[x, 3.5e+42]], $MachinePrecision]], N[(x * N[(-3.0 + N[(z * 6.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(y * N[(4.0 + N[(z * -6.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq -59000 \lor \neg \left(x \leq 3.5 \cdot 10^{+42}\right):\\
\;\;\;\;x \cdot \left(-3 + z \cdot 6\right)\\
\mathbf{else}:\\
\;\;\;\;y \cdot \left(4 + z \cdot -6\right)\\
\end{array}
\end{array}
if x < -59000 or 3.50000000000000023e42 < x Initial program 99.6%
metadata-eval99.6%
Simplified99.6%
Taylor expanded in x around inf 81.7%
sub-neg81.7%
distribute-rgt-in81.7%
metadata-eval81.7%
metadata-eval81.7%
distribute-lft-neg-in81.7%
associate-+r+81.7%
metadata-eval81.7%
metadata-eval81.7%
distribute-rgt-neg-in81.7%
metadata-eval81.7%
Simplified81.7%
if -59000 < x < 3.50000000000000023e42Initial program 99.6%
+-commutative99.6%
associate-*l*99.8%
fma-def99.8%
sub-neg99.8%
distribute-rgt-in99.8%
metadata-eval99.8%
metadata-eval99.8%
distribute-lft-neg-out99.8%
distribute-rgt-neg-in99.8%
metadata-eval99.8%
Simplified99.8%
Taylor expanded in y around inf 82.2%
Final simplification82.0%
(FPCore (x y z) :precision binary64 (if (or (<= x -0.23) (not (<= x 1e-18))) (* x -3.0) (* y 4.0)))
double code(double x, double y, double z) {
double tmp;
if ((x <= -0.23) || !(x <= 1e-18)) {
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 ((x <= (-0.23d0)) .or. (.not. (x <= 1d-18))) 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 ((x <= -0.23) || !(x <= 1e-18)) {
tmp = x * -3.0;
} else {
tmp = y * 4.0;
}
return tmp;
}
def code(x, y, z): tmp = 0 if (x <= -0.23) or not (x <= 1e-18): tmp = x * -3.0 else: tmp = y * 4.0 return tmp
function code(x, y, z) tmp = 0.0 if ((x <= -0.23) || !(x <= 1e-18)) 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 ((x <= -0.23) || ~((x <= 1e-18))) tmp = x * -3.0; else tmp = y * 4.0; end tmp_2 = tmp; end
code[x_, y_, z_] := If[Or[LessEqual[x, -0.23], N[Not[LessEqual[x, 1e-18]], $MachinePrecision]], N[(x * -3.0), $MachinePrecision], N[(y * 4.0), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq -0.23 \lor \neg \left(x \leq 10^{-18}\right):\\
\;\;\;\;x \cdot -3\\
\mathbf{else}:\\
\;\;\;\;y \cdot 4\\
\end{array}
\end{array}
if x < -0.23000000000000001 or 1.0000000000000001e-18 < x Initial program 99.6%
metadata-eval99.6%
Simplified99.6%
Taylor expanded in x around inf 77.5%
sub-neg77.5%
distribute-rgt-in77.5%
metadata-eval77.5%
metadata-eval77.5%
distribute-lft-neg-in77.5%
associate-+r+77.5%
metadata-eval77.5%
metadata-eval77.5%
distribute-rgt-neg-in77.5%
metadata-eval77.5%
Simplified77.5%
Taylor expanded in z around 0 40.7%
*-commutative40.7%
Simplified40.7%
if -0.23000000000000001 < x < 1.0000000000000001e-18Initial program 99.6%
+-commutative99.6%
associate-*l*99.8%
fma-def99.8%
sub-neg99.8%
distribute-rgt-in99.8%
metadata-eval99.8%
metadata-eval99.8%
distribute-lft-neg-out99.8%
distribute-rgt-neg-in99.8%
metadata-eval99.8%
Simplified99.8%
Taylor expanded in y around inf 82.8%
Taylor expanded in z around 0 47.0%
*-commutative47.0%
Simplified47.0%
Final simplification43.6%
(FPCore (x y z) :precision binary64 (+ x (* (* (- y x) 6.0) (- 0.6666666666666666 z))))
double code(double x, double y, double z) {
return x + (((y - x) * 6.0) * (0.6666666666666666 - 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) * (0.6666666666666666d0 - z))
end function
public static double code(double x, double y, double z) {
return x + (((y - x) * 6.0) * (0.6666666666666666 - z));
}
def code(x, y, z): return x + (((y - x) * 6.0) * (0.6666666666666666 - z))
function code(x, y, z) return Float64(x + Float64(Float64(Float64(y - x) * 6.0) * Float64(0.6666666666666666 - z))) end
function tmp = code(x, y, z) tmp = x + (((y - x) * 6.0) * (0.6666666666666666 - z)); end
code[x_, y_, z_] := N[(x + N[(N[(N[(y - x), $MachinePrecision] * 6.0), $MachinePrecision] * N[(0.6666666666666666 - z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
x + \left(\left(y - x\right) \cdot 6\right) \cdot \left(0.6666666666666666 - z\right)
\end{array}
Initial program 99.6%
metadata-eval99.6%
Simplified99.6%
Final simplification99.6%
(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%
metadata-eval99.6%
Simplified99.6%
Taylor expanded in x around inf 50.1%
sub-neg50.1%
distribute-rgt-in50.1%
metadata-eval50.1%
metadata-eval50.1%
distribute-lft-neg-in50.1%
associate-+r+50.1%
metadata-eval50.1%
metadata-eval50.1%
distribute-rgt-neg-in50.1%
metadata-eval50.1%
Simplified50.1%
Taylor expanded in z around 0 26.6%
*-commutative26.6%
Simplified26.6%
Final simplification26.6%
herbie shell --seed 2024041
(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))))