
(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.3%
+-commutative99.3%
associate-*l*99.8%
fma-define99.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 (* y z))) (t_1 (* 6.0 (* x z))))
(if (<= z -1.06e+236)
t_0
(if (<= z -7.2e+68)
t_1
(if (<= z -0.072)
t_0
(if (<= z -3e-127)
(* y 4.0)
(if (<= z -1.2e-159)
(* x -3.0)
(if (<= z 1.22e-253)
(* y 4.0)
(if (<= z 9.6e-198)
(* x -3.0)
(if (<= z 0.5) (* y 4.0) t_1))))))))))
double code(double x, double y, double z) {
double t_0 = -6.0 * (y * z);
double t_1 = 6.0 * (x * z);
double tmp;
if (z <= -1.06e+236) {
tmp = t_0;
} else if (z <= -7.2e+68) {
tmp = t_1;
} else if (z <= -0.072) {
tmp = t_0;
} else if (z <= -3e-127) {
tmp = y * 4.0;
} else if (z <= -1.2e-159) {
tmp = x * -3.0;
} else if (z <= 1.22e-253) {
tmp = y * 4.0;
} else if (z <= 9.6e-198) {
tmp = x * -3.0;
} else if (z <= 0.5) {
tmp = y * 4.0;
} else {
tmp = t_1;
}
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) * (y * z)
t_1 = 6.0d0 * (x * z)
if (z <= (-1.06d+236)) then
tmp = t_0
else if (z <= (-7.2d+68)) then
tmp = t_1
else if (z <= (-0.072d0)) then
tmp = t_0
else if (z <= (-3d-127)) then
tmp = y * 4.0d0
else if (z <= (-1.2d-159)) then
tmp = x * (-3.0d0)
else if (z <= 1.22d-253) then
tmp = y * 4.0d0
else if (z <= 9.6d-198) then
tmp = x * (-3.0d0)
else if (z <= 0.5d0) then
tmp = y * 4.0d0
else
tmp = t_1
end if
code = tmp
end function
public static double code(double x, double y, double z) {
double t_0 = -6.0 * (y * z);
double t_1 = 6.0 * (x * z);
double tmp;
if (z <= -1.06e+236) {
tmp = t_0;
} else if (z <= -7.2e+68) {
tmp = t_1;
} else if (z <= -0.072) {
tmp = t_0;
} else if (z <= -3e-127) {
tmp = y * 4.0;
} else if (z <= -1.2e-159) {
tmp = x * -3.0;
} else if (z <= 1.22e-253) {
tmp = y * 4.0;
} else if (z <= 9.6e-198) {
tmp = x * -3.0;
} else if (z <= 0.5) {
tmp = y * 4.0;
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z): t_0 = -6.0 * (y * z) t_1 = 6.0 * (x * z) tmp = 0 if z <= -1.06e+236: tmp = t_0 elif z <= -7.2e+68: tmp = t_1 elif z <= -0.072: tmp = t_0 elif z <= -3e-127: tmp = y * 4.0 elif z <= -1.2e-159: tmp = x * -3.0 elif z <= 1.22e-253: tmp = y * 4.0 elif z <= 9.6e-198: tmp = x * -3.0 elif z <= 0.5: tmp = y * 4.0 else: tmp = t_1 return tmp
function code(x, y, z) t_0 = Float64(-6.0 * Float64(y * z)) t_1 = Float64(6.0 * Float64(x * z)) tmp = 0.0 if (z <= -1.06e+236) tmp = t_0; elseif (z <= -7.2e+68) tmp = t_1; elseif (z <= -0.072) tmp = t_0; elseif (z <= -3e-127) tmp = Float64(y * 4.0); elseif (z <= -1.2e-159) tmp = Float64(x * -3.0); elseif (z <= 1.22e-253) tmp = Float64(y * 4.0); elseif (z <= 9.6e-198) tmp = Float64(x * -3.0); elseif (z <= 0.5) tmp = Float64(y * 4.0); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z) t_0 = -6.0 * (y * z); t_1 = 6.0 * (x * z); tmp = 0.0; if (z <= -1.06e+236) tmp = t_0; elseif (z <= -7.2e+68) tmp = t_1; elseif (z <= -0.072) tmp = t_0; elseif (z <= -3e-127) tmp = y * 4.0; elseif (z <= -1.2e-159) tmp = x * -3.0; elseif (z <= 1.22e-253) tmp = y * 4.0; elseif (z <= 9.6e-198) tmp = x * -3.0; elseif (z <= 0.5) tmp = y * 4.0; else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_] := Block[{t$95$0 = N[(-6.0 * N[(y * z), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(6.0 * N[(x * z), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[z, -1.06e+236], t$95$0, If[LessEqual[z, -7.2e+68], t$95$1, If[LessEqual[z, -0.072], t$95$0, If[LessEqual[z, -3e-127], N[(y * 4.0), $MachinePrecision], If[LessEqual[z, -1.2e-159], N[(x * -3.0), $MachinePrecision], If[LessEqual[z, 1.22e-253], N[(y * 4.0), $MachinePrecision], If[LessEqual[z, 9.6e-198], N[(x * -3.0), $MachinePrecision], If[LessEqual[z, 0.5], N[(y * 4.0), $MachinePrecision], t$95$1]]]]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := -6 \cdot \left(y \cdot z\right)\\
t_1 := 6 \cdot \left(x \cdot z\right)\\
\mathbf{if}\;z \leq -1.06 \cdot 10^{+236}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;z \leq -7.2 \cdot 10^{+68}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;z \leq -0.072:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;z \leq -3 \cdot 10^{-127}:\\
\;\;\;\;y \cdot 4\\
\mathbf{elif}\;z \leq -1.2 \cdot 10^{-159}:\\
\;\;\;\;x \cdot -3\\
\mathbf{elif}\;z \leq 1.22 \cdot 10^{-253}:\\
\;\;\;\;y \cdot 4\\
\mathbf{elif}\;z \leq 9.6 \cdot 10^{-198}:\\
\;\;\;\;x \cdot -3\\
\mathbf{elif}\;z \leq 0.5:\\
\;\;\;\;y \cdot 4\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if z < -1.05999999999999991e236 or -7.1999999999999998e68 < z < -0.0719999999999999946Initial program 99.8%
+-commutative99.8%
associate-*l*99.6%
fma-define99.6%
sub-neg99.6%
distribute-rgt-in99.6%
metadata-eval99.6%
metadata-eval99.6%
distribute-lft-neg-out99.6%
distribute-rgt-neg-in99.6%
metadata-eval99.6%
Simplified99.6%
Taylor expanded in y around inf 69.8%
Taylor expanded in z around inf 65.1%
if -1.05999999999999991e236 < z < -7.1999999999999998e68 or 0.5 < z Initial program 99.8%
metadata-eval99.8%
Simplified99.8%
Taylor expanded in x around inf 58.2%
+-commutative58.2%
sub-neg58.2%
distribute-rgt-in58.2%
metadata-eval58.2%
distribute-lft-neg-in58.2%
metadata-eval58.2%
neg-mul-158.2%
distribute-lft-in58.2%
*-commutative58.2%
+-commutative58.2%
mul-1-neg58.2%
*-commutative58.2%
distribute-neg-in58.2%
metadata-eval58.2%
associate-+r+58.2%
metadata-eval58.2%
distribute-rgt-neg-in58.2%
metadata-eval58.2%
Simplified58.2%
Taylor expanded in z around inf 56.7%
if -0.0719999999999999946 < z < -3.00000000000000009e-127 or -1.19999999999999999e-159 < z < 1.22e-253 or 9.59999999999999946e-198 < z < 0.5Initial program 98.6%
+-commutative98.6%
associate-*l*99.9%
fma-define99.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 63.3%
Taylor expanded in z around 0 61.5%
if -3.00000000000000009e-127 < z < -1.19999999999999999e-159 or 1.22e-253 < z < 9.59999999999999946e-198Initial program 99.5%
metadata-eval99.5%
Simplified99.5%
Taylor expanded in x around inf 85.4%
+-commutative85.4%
sub-neg85.4%
distribute-rgt-in85.4%
metadata-eval85.4%
distribute-lft-neg-in85.4%
metadata-eval85.4%
neg-mul-185.4%
distribute-lft-in85.4%
*-commutative85.4%
+-commutative85.4%
mul-1-neg85.4%
*-commutative85.4%
distribute-neg-in85.4%
metadata-eval85.4%
associate-+r+85.4%
metadata-eval85.4%
distribute-rgt-neg-in85.4%
metadata-eval85.4%
Simplified85.4%
Taylor expanded in z around 0 85.4%
*-commutative85.4%
Simplified85.4%
Final simplification62.0%
(FPCore (x y z)
:precision binary64
(let* ((t_0 (* -6.0 (* y z))) (t_1 (* x (* z 6.0))))
(if (<= z -1.26e+237)
t_0
(if (<= z -9.5e+65)
t_1
(if (<= z -0.072)
t_0
(if (<= z -6.2e-124)
(* y 4.0)
(if (<= z -8.2e-151)
(* x -3.0)
(if (<= z 2.65e-253)
(* y 4.0)
(if (<= z 3.3e-193)
(* x -3.0)
(if (<= z 0.62) (* y 4.0) t_1))))))))))
double code(double x, double y, double z) {
double t_0 = -6.0 * (y * z);
double t_1 = x * (z * 6.0);
double tmp;
if (z <= -1.26e+237) {
tmp = t_0;
} else if (z <= -9.5e+65) {
tmp = t_1;
} else if (z <= -0.072) {
tmp = t_0;
} else if (z <= -6.2e-124) {
tmp = y * 4.0;
} else if (z <= -8.2e-151) {
tmp = x * -3.0;
} else if (z <= 2.65e-253) {
tmp = y * 4.0;
} else if (z <= 3.3e-193) {
tmp = x * -3.0;
} else if (z <= 0.62) {
tmp = y * 4.0;
} else {
tmp = t_1;
}
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) * (y * z)
t_1 = x * (z * 6.0d0)
if (z <= (-1.26d+237)) then
tmp = t_0
else if (z <= (-9.5d+65)) then
tmp = t_1
else if (z <= (-0.072d0)) then
tmp = t_0
else if (z <= (-6.2d-124)) then
tmp = y * 4.0d0
else if (z <= (-8.2d-151)) then
tmp = x * (-3.0d0)
else if (z <= 2.65d-253) then
tmp = y * 4.0d0
else if (z <= 3.3d-193) then
tmp = x * (-3.0d0)
else if (z <= 0.62d0) then
tmp = y * 4.0d0
else
tmp = t_1
end if
code = tmp
end function
public static double code(double x, double y, double z) {
double t_0 = -6.0 * (y * z);
double t_1 = x * (z * 6.0);
double tmp;
if (z <= -1.26e+237) {
tmp = t_0;
} else if (z <= -9.5e+65) {
tmp = t_1;
} else if (z <= -0.072) {
tmp = t_0;
} else if (z <= -6.2e-124) {
tmp = y * 4.0;
} else if (z <= -8.2e-151) {
tmp = x * -3.0;
} else if (z <= 2.65e-253) {
tmp = y * 4.0;
} else if (z <= 3.3e-193) {
tmp = x * -3.0;
} else if (z <= 0.62) {
tmp = y * 4.0;
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z): t_0 = -6.0 * (y * z) t_1 = x * (z * 6.0) tmp = 0 if z <= -1.26e+237: tmp = t_0 elif z <= -9.5e+65: tmp = t_1 elif z <= -0.072: tmp = t_0 elif z <= -6.2e-124: tmp = y * 4.0 elif z <= -8.2e-151: tmp = x * -3.0 elif z <= 2.65e-253: tmp = y * 4.0 elif z <= 3.3e-193: tmp = x * -3.0 elif z <= 0.62: tmp = y * 4.0 else: tmp = t_1 return tmp
function code(x, y, z) t_0 = Float64(-6.0 * Float64(y * z)) t_1 = Float64(x * Float64(z * 6.0)) tmp = 0.0 if (z <= -1.26e+237) tmp = t_0; elseif (z <= -9.5e+65) tmp = t_1; elseif (z <= -0.072) tmp = t_0; elseif (z <= -6.2e-124) tmp = Float64(y * 4.0); elseif (z <= -8.2e-151) tmp = Float64(x * -3.0); elseif (z <= 2.65e-253) tmp = Float64(y * 4.0); elseif (z <= 3.3e-193) tmp = Float64(x * -3.0); elseif (z <= 0.62) tmp = Float64(y * 4.0); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z) t_0 = -6.0 * (y * z); t_1 = x * (z * 6.0); tmp = 0.0; if (z <= -1.26e+237) tmp = t_0; elseif (z <= -9.5e+65) tmp = t_1; elseif (z <= -0.072) tmp = t_0; elseif (z <= -6.2e-124) tmp = y * 4.0; elseif (z <= -8.2e-151) tmp = x * -3.0; elseif (z <= 2.65e-253) tmp = y * 4.0; elseif (z <= 3.3e-193) tmp = x * -3.0; elseif (z <= 0.62) tmp = y * 4.0; else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_] := Block[{t$95$0 = N[(-6.0 * N[(y * z), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(x * N[(z * 6.0), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[z, -1.26e+237], t$95$0, If[LessEqual[z, -9.5e+65], t$95$1, If[LessEqual[z, -0.072], t$95$0, If[LessEqual[z, -6.2e-124], N[(y * 4.0), $MachinePrecision], If[LessEqual[z, -8.2e-151], N[(x * -3.0), $MachinePrecision], If[LessEqual[z, 2.65e-253], N[(y * 4.0), $MachinePrecision], If[LessEqual[z, 3.3e-193], N[(x * -3.0), $MachinePrecision], If[LessEqual[z, 0.62], N[(y * 4.0), $MachinePrecision], t$95$1]]]]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := -6 \cdot \left(y \cdot z\right)\\
t_1 := x \cdot \left(z \cdot 6\right)\\
\mathbf{if}\;z \leq -1.26 \cdot 10^{+237}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;z \leq -9.5 \cdot 10^{+65}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;z \leq -0.072:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;z \leq -6.2 \cdot 10^{-124}:\\
\;\;\;\;y \cdot 4\\
\mathbf{elif}\;z \leq -8.2 \cdot 10^{-151}:\\
\;\;\;\;x \cdot -3\\
\mathbf{elif}\;z \leq 2.65 \cdot 10^{-253}:\\
\;\;\;\;y \cdot 4\\
\mathbf{elif}\;z \leq 3.3 \cdot 10^{-193}:\\
\;\;\;\;x \cdot -3\\
\mathbf{elif}\;z \leq 0.62:\\
\;\;\;\;y \cdot 4\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if z < -1.2599999999999999e237 or -9.5000000000000005e65 < z < -0.0719999999999999946Initial program 99.8%
+-commutative99.8%
associate-*l*99.6%
fma-define99.6%
sub-neg99.6%
distribute-rgt-in99.6%
metadata-eval99.6%
metadata-eval99.6%
distribute-lft-neg-out99.6%
distribute-rgt-neg-in99.6%
metadata-eval99.6%
Simplified99.6%
Taylor expanded in y around inf 71.0%
Taylor expanded in z around inf 66.0%
if -1.2599999999999999e237 < z < -9.5000000000000005e65 or 0.619999999999999996 < z Initial program 99.8%
metadata-eval99.8%
Simplified99.8%
Taylor expanded in x around inf 58.1%
+-commutative58.1%
sub-neg58.1%
distribute-rgt-in58.1%
metadata-eval58.1%
distribute-lft-neg-in58.1%
metadata-eval58.1%
neg-mul-158.1%
distribute-lft-in58.1%
*-commutative58.1%
+-commutative58.1%
mul-1-neg58.1%
*-commutative58.1%
distribute-neg-in58.1%
metadata-eval58.1%
associate-+r+58.1%
metadata-eval58.1%
distribute-rgt-neg-in58.1%
metadata-eval58.1%
Simplified58.1%
Taylor expanded in z around inf 56.6%
*-commutative56.6%
associate-*r*56.7%
Simplified56.7%
if -0.0719999999999999946 < z < -6.1999999999999996e-124 or -8.2000000000000002e-151 < z < 2.6500000000000001e-253 or 3.2999999999999999e-193 < z < 0.619999999999999996Initial program 98.6%
+-commutative98.6%
associate-*l*99.9%
fma-define99.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 63.3%
Taylor expanded in z around 0 61.5%
if -6.1999999999999996e-124 < z < -8.2000000000000002e-151 or 2.6500000000000001e-253 < z < 3.2999999999999999e-193Initial program 99.5%
metadata-eval99.5%
Simplified99.5%
Taylor expanded in x around inf 85.4%
+-commutative85.4%
sub-neg85.4%
distribute-rgt-in85.4%
metadata-eval85.4%
distribute-lft-neg-in85.4%
metadata-eval85.4%
neg-mul-185.4%
distribute-lft-in85.4%
*-commutative85.4%
+-commutative85.4%
mul-1-neg85.4%
*-commutative85.4%
distribute-neg-in85.4%
metadata-eval85.4%
associate-+r+85.4%
metadata-eval85.4%
distribute-rgt-neg-in85.4%
metadata-eval85.4%
Simplified85.4%
Taylor expanded in z around 0 85.4%
*-commutative85.4%
Simplified85.4%
Final simplification62.0%
(FPCore (x y z)
:precision binary64
(let* ((t_0 (* x (* z 6.0))))
(if (<= z -6.5e+228)
(* -6.0 (* y z))
(if (<= z -4.1e+73)
t_0
(if (<= z -0.072)
(* y (* z -6.0))
(if (<= z -3.9e-127)
(* y 4.0)
(if (<= z -2.75e-157)
(* x -3.0)
(if (<= z 2.5e-253)
(* y 4.0)
(if (<= z 1.75e-197)
(* x -3.0)
(if (<= z 0.5) (* y 4.0) t_0))))))))))
double code(double x, double y, double z) {
double t_0 = x * (z * 6.0);
double tmp;
if (z <= -6.5e+228) {
tmp = -6.0 * (y * z);
} else if (z <= -4.1e+73) {
tmp = t_0;
} else if (z <= -0.072) {
tmp = y * (z * -6.0);
} else if (z <= -3.9e-127) {
tmp = y * 4.0;
} else if (z <= -2.75e-157) {
tmp = x * -3.0;
} else if (z <= 2.5e-253) {
tmp = y * 4.0;
} else if (z <= 1.75e-197) {
tmp = x * -3.0;
} else if (z <= 0.5) {
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 = x * (z * 6.0d0)
if (z <= (-6.5d+228)) then
tmp = (-6.0d0) * (y * z)
else if (z <= (-4.1d+73)) then
tmp = t_0
else if (z <= (-0.072d0)) then
tmp = y * (z * (-6.0d0))
else if (z <= (-3.9d-127)) then
tmp = y * 4.0d0
else if (z <= (-2.75d-157)) then
tmp = x * (-3.0d0)
else if (z <= 2.5d-253) then
tmp = y * 4.0d0
else if (z <= 1.75d-197) then
tmp = x * (-3.0d0)
else if (z <= 0.5d0) 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 = x * (z * 6.0);
double tmp;
if (z <= -6.5e+228) {
tmp = -6.0 * (y * z);
} else if (z <= -4.1e+73) {
tmp = t_0;
} else if (z <= -0.072) {
tmp = y * (z * -6.0);
} else if (z <= -3.9e-127) {
tmp = y * 4.0;
} else if (z <= -2.75e-157) {
tmp = x * -3.0;
} else if (z <= 2.5e-253) {
tmp = y * 4.0;
} else if (z <= 1.75e-197) {
tmp = x * -3.0;
} else if (z <= 0.5) {
tmp = y * 4.0;
} else {
tmp = t_0;
}
return tmp;
}
def code(x, y, z): t_0 = x * (z * 6.0) tmp = 0 if z <= -6.5e+228: tmp = -6.0 * (y * z) elif z <= -4.1e+73: tmp = t_0 elif z <= -0.072: tmp = y * (z * -6.0) elif z <= -3.9e-127: tmp = y * 4.0 elif z <= -2.75e-157: tmp = x * -3.0 elif z <= 2.5e-253: tmp = y * 4.0 elif z <= 1.75e-197: tmp = x * -3.0 elif z <= 0.5: tmp = y * 4.0 else: tmp = t_0 return tmp
function code(x, y, z) t_0 = Float64(x * Float64(z * 6.0)) tmp = 0.0 if (z <= -6.5e+228) tmp = Float64(-6.0 * Float64(y * z)); elseif (z <= -4.1e+73) tmp = t_0; elseif (z <= -0.072) tmp = Float64(y * Float64(z * -6.0)); elseif (z <= -3.9e-127) tmp = Float64(y * 4.0); elseif (z <= -2.75e-157) tmp = Float64(x * -3.0); elseif (z <= 2.5e-253) tmp = Float64(y * 4.0); elseif (z <= 1.75e-197) tmp = Float64(x * -3.0); elseif (z <= 0.5) tmp = Float64(y * 4.0); else tmp = t_0; end return tmp end
function tmp_2 = code(x, y, z) t_0 = x * (z * 6.0); tmp = 0.0; if (z <= -6.5e+228) tmp = -6.0 * (y * z); elseif (z <= -4.1e+73) tmp = t_0; elseif (z <= -0.072) tmp = y * (z * -6.0); elseif (z <= -3.9e-127) tmp = y * 4.0; elseif (z <= -2.75e-157) tmp = x * -3.0; elseif (z <= 2.5e-253) tmp = y * 4.0; elseif (z <= 1.75e-197) tmp = x * -3.0; elseif (z <= 0.5) tmp = y * 4.0; else tmp = t_0; end tmp_2 = tmp; end
code[x_, y_, z_] := Block[{t$95$0 = N[(x * N[(z * 6.0), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[z, -6.5e+228], N[(-6.0 * N[(y * z), $MachinePrecision]), $MachinePrecision], If[LessEqual[z, -4.1e+73], t$95$0, If[LessEqual[z, -0.072], N[(y * N[(z * -6.0), $MachinePrecision]), $MachinePrecision], If[LessEqual[z, -3.9e-127], N[(y * 4.0), $MachinePrecision], If[LessEqual[z, -2.75e-157], N[(x * -3.0), $MachinePrecision], If[LessEqual[z, 2.5e-253], N[(y * 4.0), $MachinePrecision], If[LessEqual[z, 1.75e-197], N[(x * -3.0), $MachinePrecision], If[LessEqual[z, 0.5], N[(y * 4.0), $MachinePrecision], t$95$0]]]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := x \cdot \left(z \cdot 6\right)\\
\mathbf{if}\;z \leq -6.5 \cdot 10^{+228}:\\
\;\;\;\;-6 \cdot \left(y \cdot z\right)\\
\mathbf{elif}\;z \leq -4.1 \cdot 10^{+73}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;z \leq -0.072:\\
\;\;\;\;y \cdot \left(z \cdot -6\right)\\
\mathbf{elif}\;z \leq -3.9 \cdot 10^{-127}:\\
\;\;\;\;y \cdot 4\\
\mathbf{elif}\;z \leq -2.75 \cdot 10^{-157}:\\
\;\;\;\;x \cdot -3\\
\mathbf{elif}\;z \leq 2.5 \cdot 10^{-253}:\\
\;\;\;\;y \cdot 4\\
\mathbf{elif}\;z \leq 1.75 \cdot 10^{-197}:\\
\;\;\;\;x \cdot -3\\
\mathbf{elif}\;z \leq 0.5:\\
\;\;\;\;y \cdot 4\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
if z < -6.5e228Initial program 100.0%
+-commutative100.0%
associate-*l*99.8%
fma-define99.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 75.4%
Taylor expanded in z around inf 75.4%
if -6.5e228 < z < -4.0999999999999998e73 or 0.5 < z Initial program 99.8%
metadata-eval99.8%
Simplified99.8%
Taylor expanded in x around inf 58.2%
+-commutative58.2%
sub-neg58.2%
distribute-rgt-in58.2%
metadata-eval58.2%
distribute-lft-neg-in58.2%
metadata-eval58.2%
neg-mul-158.2%
distribute-lft-in58.2%
*-commutative58.2%
+-commutative58.2%
mul-1-neg58.2%
*-commutative58.2%
distribute-neg-in58.2%
metadata-eval58.2%
associate-+r+58.2%
metadata-eval58.2%
distribute-rgt-neg-in58.2%
metadata-eval58.2%
Simplified58.2%
Taylor expanded in z around inf 56.7%
*-commutative56.7%
associate-*r*56.8%
Simplified56.8%
if -4.0999999999999998e73 < z < -0.0719999999999999946Initial program 99.6%
+-commutative99.6%
associate-*l*99.4%
fma-define99.4%
sub-neg99.4%
distribute-rgt-in99.5%
metadata-eval99.5%
metadata-eval99.5%
distribute-lft-neg-out99.5%
distribute-rgt-neg-in99.5%
metadata-eval99.5%
Simplified99.5%
Taylor expanded in y around inf 64.9%
Taylor expanded in z around inf 56.0%
*-commutative56.0%
associate-*r*56.1%
*-commutative56.1%
Simplified56.1%
if -0.0719999999999999946 < z < -3.89999999999999966e-127 or -2.7499999999999999e-157 < z < 2.49999999999999986e-253 or 1.7499999999999999e-197 < z < 0.5Initial program 98.6%
+-commutative98.6%
associate-*l*99.9%
fma-define99.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 63.3%
Taylor expanded in z around 0 61.5%
if -3.89999999999999966e-127 < z < -2.7499999999999999e-157 or 2.49999999999999986e-253 < z < 1.7499999999999999e-197Initial program 99.5%
metadata-eval99.5%
Simplified99.5%
Taylor expanded in x around inf 85.4%
+-commutative85.4%
sub-neg85.4%
distribute-rgt-in85.4%
metadata-eval85.4%
distribute-lft-neg-in85.4%
metadata-eval85.4%
neg-mul-185.4%
distribute-lft-in85.4%
*-commutative85.4%
+-commutative85.4%
mul-1-neg85.4%
*-commutative85.4%
distribute-neg-in85.4%
metadata-eval85.4%
associate-+r+85.4%
metadata-eval85.4%
distribute-rgt-neg-in85.4%
metadata-eval85.4%
Simplified85.4%
Taylor expanded in z around 0 85.4%
*-commutative85.4%
Simplified85.4%
Final simplification62.0%
(FPCore (x y z)
:precision binary64
(if (<= z -1.12e+228)
(* -6.0 (* y z))
(if (<= z -4.9e+71)
(* x (* z 6.0))
(if (<= z -0.072)
(* y (* z -6.0))
(if (<= z -1.7e-123)
(* y 4.0)
(if (<= z -2.95e-151)
(* x -3.0)
(if (<= z 2.7e-253)
(* y 4.0)
(if (<= z 5.6e-196)
(* x -3.0)
(if (<= z 0.65) (* y 4.0) (* z (* x 6.0)))))))))))
double code(double x, double y, double z) {
double tmp;
if (z <= -1.12e+228) {
tmp = -6.0 * (y * z);
} else if (z <= -4.9e+71) {
tmp = x * (z * 6.0);
} else if (z <= -0.072) {
tmp = y * (z * -6.0);
} else if (z <= -1.7e-123) {
tmp = y * 4.0;
} else if (z <= -2.95e-151) {
tmp = x * -3.0;
} else if (z <= 2.7e-253) {
tmp = y * 4.0;
} else if (z <= 5.6e-196) {
tmp = x * -3.0;
} else if (z <= 0.65) {
tmp = y * 4.0;
} else {
tmp = z * (x * 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 (z <= (-1.12d+228)) then
tmp = (-6.0d0) * (y * z)
else if (z <= (-4.9d+71)) then
tmp = x * (z * 6.0d0)
else if (z <= (-0.072d0)) then
tmp = y * (z * (-6.0d0))
else if (z <= (-1.7d-123)) then
tmp = y * 4.0d0
else if (z <= (-2.95d-151)) then
tmp = x * (-3.0d0)
else if (z <= 2.7d-253) then
tmp = y * 4.0d0
else if (z <= 5.6d-196) then
tmp = x * (-3.0d0)
else if (z <= 0.65d0) then
tmp = y * 4.0d0
else
tmp = z * (x * 6.0d0)
end if
code = tmp
end function
public static double code(double x, double y, double z) {
double tmp;
if (z <= -1.12e+228) {
tmp = -6.0 * (y * z);
} else if (z <= -4.9e+71) {
tmp = x * (z * 6.0);
} else if (z <= -0.072) {
tmp = y * (z * -6.0);
} else if (z <= -1.7e-123) {
tmp = y * 4.0;
} else if (z <= -2.95e-151) {
tmp = x * -3.0;
} else if (z <= 2.7e-253) {
tmp = y * 4.0;
} else if (z <= 5.6e-196) {
tmp = x * -3.0;
} else if (z <= 0.65) {
tmp = y * 4.0;
} else {
tmp = z * (x * 6.0);
}
return tmp;
}
def code(x, y, z): tmp = 0 if z <= -1.12e+228: tmp = -6.0 * (y * z) elif z <= -4.9e+71: tmp = x * (z * 6.0) elif z <= -0.072: tmp = y * (z * -6.0) elif z <= -1.7e-123: tmp = y * 4.0 elif z <= -2.95e-151: tmp = x * -3.0 elif z <= 2.7e-253: tmp = y * 4.0 elif z <= 5.6e-196: tmp = x * -3.0 elif z <= 0.65: tmp = y * 4.0 else: tmp = z * (x * 6.0) return tmp
function code(x, y, z) tmp = 0.0 if (z <= -1.12e+228) tmp = Float64(-6.0 * Float64(y * z)); elseif (z <= -4.9e+71) tmp = Float64(x * Float64(z * 6.0)); elseif (z <= -0.072) tmp = Float64(y * Float64(z * -6.0)); elseif (z <= -1.7e-123) tmp = Float64(y * 4.0); elseif (z <= -2.95e-151) tmp = Float64(x * -3.0); elseif (z <= 2.7e-253) tmp = Float64(y * 4.0); elseif (z <= 5.6e-196) tmp = Float64(x * -3.0); elseif (z <= 0.65) tmp = Float64(y * 4.0); else tmp = Float64(z * Float64(x * 6.0)); end return tmp end
function tmp_2 = code(x, y, z) tmp = 0.0; if (z <= -1.12e+228) tmp = -6.0 * (y * z); elseif (z <= -4.9e+71) tmp = x * (z * 6.0); elseif (z <= -0.072) tmp = y * (z * -6.0); elseif (z <= -1.7e-123) tmp = y * 4.0; elseif (z <= -2.95e-151) tmp = x * -3.0; elseif (z <= 2.7e-253) tmp = y * 4.0; elseif (z <= 5.6e-196) tmp = x * -3.0; elseif (z <= 0.65) tmp = y * 4.0; else tmp = z * (x * 6.0); end tmp_2 = tmp; end
code[x_, y_, z_] := If[LessEqual[z, -1.12e+228], N[(-6.0 * N[(y * z), $MachinePrecision]), $MachinePrecision], If[LessEqual[z, -4.9e+71], N[(x * N[(z * 6.0), $MachinePrecision]), $MachinePrecision], If[LessEqual[z, -0.072], N[(y * N[(z * -6.0), $MachinePrecision]), $MachinePrecision], If[LessEqual[z, -1.7e-123], N[(y * 4.0), $MachinePrecision], If[LessEqual[z, -2.95e-151], N[(x * -3.0), $MachinePrecision], If[LessEqual[z, 2.7e-253], N[(y * 4.0), $MachinePrecision], If[LessEqual[z, 5.6e-196], N[(x * -3.0), $MachinePrecision], If[LessEqual[z, 0.65], N[(y * 4.0), $MachinePrecision], N[(z * N[(x * 6.0), $MachinePrecision]), $MachinePrecision]]]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;z \leq -1.12 \cdot 10^{+228}:\\
\;\;\;\;-6 \cdot \left(y \cdot z\right)\\
\mathbf{elif}\;z \leq -4.9 \cdot 10^{+71}:\\
\;\;\;\;x \cdot \left(z \cdot 6\right)\\
\mathbf{elif}\;z \leq -0.072:\\
\;\;\;\;y \cdot \left(z \cdot -6\right)\\
\mathbf{elif}\;z \leq -1.7 \cdot 10^{-123}:\\
\;\;\;\;y \cdot 4\\
\mathbf{elif}\;z \leq -2.95 \cdot 10^{-151}:\\
\;\;\;\;x \cdot -3\\
\mathbf{elif}\;z \leq 2.7 \cdot 10^{-253}:\\
\;\;\;\;y \cdot 4\\
\mathbf{elif}\;z \leq 5.6 \cdot 10^{-196}:\\
\;\;\;\;x \cdot -3\\
\mathbf{elif}\;z \leq 0.65:\\
\;\;\;\;y \cdot 4\\
\mathbf{else}:\\
\;\;\;\;z \cdot \left(x \cdot 6\right)\\
\end{array}
\end{array}
if z < -1.11999999999999994e228Initial program 100.0%
+-commutative100.0%
associate-*l*99.8%
fma-define99.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 75.4%
Taylor expanded in z around inf 75.4%
if -1.11999999999999994e228 < z < -4.8999999999999997e71Initial program 99.8%
metadata-eval99.8%
Simplified99.8%
Taylor expanded in x around inf 63.8%
+-commutative63.8%
sub-neg63.8%
distribute-rgt-in63.8%
metadata-eval63.8%
distribute-lft-neg-in63.8%
metadata-eval63.8%
neg-mul-163.8%
distribute-lft-in63.8%
*-commutative63.8%
+-commutative63.8%
mul-1-neg63.8%
*-commutative63.8%
distribute-neg-in63.8%
metadata-eval63.8%
associate-+r+63.8%
metadata-eval63.8%
distribute-rgt-neg-in63.8%
metadata-eval63.8%
Simplified63.8%
Taylor expanded in z around inf 63.7%
*-commutative63.7%
associate-*r*63.8%
Simplified63.8%
if -4.8999999999999997e71 < z < -0.0719999999999999946Initial program 99.6%
+-commutative99.6%
associate-*l*99.4%
fma-define99.4%
sub-neg99.4%
distribute-rgt-in99.5%
metadata-eval99.5%
metadata-eval99.5%
distribute-lft-neg-out99.5%
distribute-rgt-neg-in99.5%
metadata-eval99.5%
Simplified99.5%
Taylor expanded in y around inf 64.9%
Taylor expanded in z around inf 56.0%
*-commutative56.0%
associate-*r*56.1%
*-commutative56.1%
Simplified56.1%
if -0.0719999999999999946 < z < -1.7e-123 or -2.95e-151 < z < 2.69999999999999999e-253 or 5.5999999999999997e-196 < z < 0.650000000000000022Initial program 98.6%
+-commutative98.6%
associate-*l*99.9%
fma-define99.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 63.3%
Taylor expanded in z around 0 61.5%
if -1.7e-123 < z < -2.95e-151 or 2.69999999999999999e-253 < z < 5.5999999999999997e-196Initial program 99.5%
metadata-eval99.5%
Simplified99.5%
Taylor expanded in x around inf 85.4%
+-commutative85.4%
sub-neg85.4%
distribute-rgt-in85.4%
metadata-eval85.4%
distribute-lft-neg-in85.4%
metadata-eval85.4%
neg-mul-185.4%
distribute-lft-in85.4%
*-commutative85.4%
+-commutative85.4%
mul-1-neg85.4%
*-commutative85.4%
distribute-neg-in85.4%
metadata-eval85.4%
associate-+r+85.4%
metadata-eval85.4%
distribute-rgt-neg-in85.4%
metadata-eval85.4%
Simplified85.4%
Taylor expanded in z around 0 85.4%
*-commutative85.4%
Simplified85.4%
if 0.650000000000000022 < z Initial program 99.7%
metadata-eval99.7%
Simplified99.7%
Taylor expanded in x around inf 56.1%
+-commutative56.1%
sub-neg56.1%
distribute-rgt-in56.1%
metadata-eval56.1%
distribute-lft-neg-in56.1%
metadata-eval56.1%
neg-mul-156.1%
distribute-lft-in56.1%
*-commutative56.1%
+-commutative56.1%
mul-1-neg56.1%
*-commutative56.1%
distribute-neg-in56.1%
metadata-eval56.1%
associate-+r+56.1%
metadata-eval56.1%
distribute-rgt-neg-in56.1%
metadata-eval56.1%
Simplified56.1%
Taylor expanded in z around inf 54.1%
*-commutative54.1%
associate-*r*54.1%
Simplified54.1%
Taylor expanded in x around 0 54.1%
associate-*r*54.1%
*-commutative54.1%
Simplified54.1%
Final simplification62.0%
(FPCore (x y z)
:precision binary64
(let* ((t_0 (* -6.0 (* (- y x) z))))
(if (<= z -0.04)
t_0
(if (<= z -7.5e-128)
(* y 4.0)
(if (<= z -1.52e-161)
(* x -3.0)
(if (<= z 9.8e-254)
(* y 4.0)
(if (<= z 2.95e-197) (* x -3.0) (if (<= z 0.5) (* y 4.0) t_0))))))))
double code(double x, double y, double z) {
double t_0 = -6.0 * ((y - x) * z);
double tmp;
if (z <= -0.04) {
tmp = t_0;
} else if (z <= -7.5e-128) {
tmp = y * 4.0;
} else if (z <= -1.52e-161) {
tmp = x * -3.0;
} else if (z <= 9.8e-254) {
tmp = y * 4.0;
} else if (z <= 2.95e-197) {
tmp = x * -3.0;
} else if (z <= 0.5) {
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 - x) * z)
if (z <= (-0.04d0)) then
tmp = t_0
else if (z <= (-7.5d-128)) then
tmp = y * 4.0d0
else if (z <= (-1.52d-161)) then
tmp = x * (-3.0d0)
else if (z <= 9.8d-254) then
tmp = y * 4.0d0
else if (z <= 2.95d-197) then
tmp = x * (-3.0d0)
else if (z <= 0.5d0) 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 - x) * z);
double tmp;
if (z <= -0.04) {
tmp = t_0;
} else if (z <= -7.5e-128) {
tmp = y * 4.0;
} else if (z <= -1.52e-161) {
tmp = x * -3.0;
} else if (z <= 9.8e-254) {
tmp = y * 4.0;
} else if (z <= 2.95e-197) {
tmp = x * -3.0;
} else if (z <= 0.5) {
tmp = y * 4.0;
} else {
tmp = t_0;
}
return tmp;
}
def code(x, y, z): t_0 = -6.0 * ((y - x) * z) tmp = 0 if z <= -0.04: tmp = t_0 elif z <= -7.5e-128: tmp = y * 4.0 elif z <= -1.52e-161: tmp = x * -3.0 elif z <= 9.8e-254: tmp = y * 4.0 elif z <= 2.95e-197: tmp = x * -3.0 elif z <= 0.5: tmp = y * 4.0 else: tmp = t_0 return tmp
function code(x, y, z) t_0 = Float64(-6.0 * Float64(Float64(y - x) * z)) tmp = 0.0 if (z <= -0.04) tmp = t_0; elseif (z <= -7.5e-128) tmp = Float64(y * 4.0); elseif (z <= -1.52e-161) tmp = Float64(x * -3.0); elseif (z <= 9.8e-254) tmp = Float64(y * 4.0); elseif (z <= 2.95e-197) tmp = Float64(x * -3.0); elseif (z <= 0.5) 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 - x) * z); tmp = 0.0; if (z <= -0.04) tmp = t_0; elseif (z <= -7.5e-128) tmp = y * 4.0; elseif (z <= -1.52e-161) tmp = x * -3.0; elseif (z <= 9.8e-254) tmp = y * 4.0; elseif (z <= 2.95e-197) tmp = x * -3.0; elseif (z <= 0.5) 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[(N[(y - x), $MachinePrecision] * z), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[z, -0.04], t$95$0, If[LessEqual[z, -7.5e-128], N[(y * 4.0), $MachinePrecision], If[LessEqual[z, -1.52e-161], N[(x * -3.0), $MachinePrecision], If[LessEqual[z, 9.8e-254], N[(y * 4.0), $MachinePrecision], If[LessEqual[z, 2.95e-197], N[(x * -3.0), $MachinePrecision], If[LessEqual[z, 0.5], N[(y * 4.0), $MachinePrecision], t$95$0]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := -6 \cdot \left(\left(y - x\right) \cdot z\right)\\
\mathbf{if}\;z \leq -0.04:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;z \leq -7.5 \cdot 10^{-128}:\\
\;\;\;\;y \cdot 4\\
\mathbf{elif}\;z \leq -1.52 \cdot 10^{-161}:\\
\;\;\;\;x \cdot -3\\
\mathbf{elif}\;z \leq 9.8 \cdot 10^{-254}:\\
\;\;\;\;y \cdot 4\\
\mathbf{elif}\;z \leq 2.95 \cdot 10^{-197}:\\
\;\;\;\;x \cdot -3\\
\mathbf{elif}\;z \leq 0.5:\\
\;\;\;\;y \cdot 4\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
if z < -0.0400000000000000008 or 0.5 < z Initial program 99.8%
metadata-eval99.8%
Simplified99.8%
Taylor expanded in z around inf 96.3%
Taylor expanded in z around inf 96.4%
if -0.0400000000000000008 < z < -7.50000000000000021e-128 or -1.52000000000000002e-161 < z < 9.79999999999999959e-254 or 2.95000000000000024e-197 < z < 0.5Initial program 98.6%
+-commutative98.6%
associate-*l*99.9%
fma-define99.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 63.3%
Taylor expanded in z around 0 61.5%
if -7.50000000000000021e-128 < z < -1.52000000000000002e-161 or 9.79999999999999959e-254 < z < 2.95000000000000024e-197Initial program 99.5%
metadata-eval99.5%
Simplified99.5%
Taylor expanded in x around inf 85.4%
+-commutative85.4%
sub-neg85.4%
distribute-rgt-in85.4%
metadata-eval85.4%
distribute-lft-neg-in85.4%
metadata-eval85.4%
neg-mul-185.4%
distribute-lft-in85.4%
*-commutative85.4%
+-commutative85.4%
mul-1-neg85.4%
*-commutative85.4%
distribute-neg-in85.4%
metadata-eval85.4%
associate-+r+85.4%
metadata-eval85.4%
distribute-rgt-neg-in85.4%
metadata-eval85.4%
Simplified85.4%
Taylor expanded in z around 0 85.4%
*-commutative85.4%
Simplified85.4%
Final simplification81.5%
(FPCore (x y z)
:precision binary64
(let* ((t_0 (* -6.0 (* y z))))
(if (<= z -0.072)
t_0
(if (<= z -2e-126)
(* y 4.0)
(if (<= z -8.8e-154)
(* x -3.0)
(if (<= z 1.75e-254)
(* y 4.0)
(if (<= z 2.5e-197)
(* x -3.0)
(if (<= z 42000.0) (* 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.072) {
tmp = t_0;
} else if (z <= -2e-126) {
tmp = y * 4.0;
} else if (z <= -8.8e-154) {
tmp = x * -3.0;
} else if (z <= 1.75e-254) {
tmp = y * 4.0;
} else if (z <= 2.5e-197) {
tmp = x * -3.0;
} else if (z <= 42000.0) {
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.072d0)) then
tmp = t_0
else if (z <= (-2d-126)) then
tmp = y * 4.0d0
else if (z <= (-8.8d-154)) then
tmp = x * (-3.0d0)
else if (z <= 1.75d-254) then
tmp = y * 4.0d0
else if (z <= 2.5d-197) then
tmp = x * (-3.0d0)
else if (z <= 42000.0d0) 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.072) {
tmp = t_0;
} else if (z <= -2e-126) {
tmp = y * 4.0;
} else if (z <= -8.8e-154) {
tmp = x * -3.0;
} else if (z <= 1.75e-254) {
tmp = y * 4.0;
} else if (z <= 2.5e-197) {
tmp = x * -3.0;
} else if (z <= 42000.0) {
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.072: tmp = t_0 elif z <= -2e-126: tmp = y * 4.0 elif z <= -8.8e-154: tmp = x * -3.0 elif z <= 1.75e-254: tmp = y * 4.0 elif z <= 2.5e-197: tmp = x * -3.0 elif z <= 42000.0: 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.072) tmp = t_0; elseif (z <= -2e-126) tmp = Float64(y * 4.0); elseif (z <= -8.8e-154) tmp = Float64(x * -3.0); elseif (z <= 1.75e-254) tmp = Float64(y * 4.0); elseif (z <= 2.5e-197) tmp = Float64(x * -3.0); elseif (z <= 42000.0) 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.072) tmp = t_0; elseif (z <= -2e-126) tmp = y * 4.0; elseif (z <= -8.8e-154) tmp = x * -3.0; elseif (z <= 1.75e-254) tmp = y * 4.0; elseif (z <= 2.5e-197) tmp = x * -3.0; elseif (z <= 42000.0) 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.072], t$95$0, If[LessEqual[z, -2e-126], N[(y * 4.0), $MachinePrecision], If[LessEqual[z, -8.8e-154], N[(x * -3.0), $MachinePrecision], If[LessEqual[z, 1.75e-254], N[(y * 4.0), $MachinePrecision], If[LessEqual[z, 2.5e-197], N[(x * -3.0), $MachinePrecision], If[LessEqual[z, 42000.0], 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.072:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;z \leq -2 \cdot 10^{-126}:\\
\;\;\;\;y \cdot 4\\
\mathbf{elif}\;z \leq -8.8 \cdot 10^{-154}:\\
\;\;\;\;x \cdot -3\\
\mathbf{elif}\;z \leq 1.75 \cdot 10^{-254}:\\
\;\;\;\;y \cdot 4\\
\mathbf{elif}\;z \leq 2.5 \cdot 10^{-197}:\\
\;\;\;\;x \cdot -3\\
\mathbf{elif}\;z \leq 42000:\\
\;\;\;\;y \cdot 4\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
if z < -0.0719999999999999946 or 42000 < z Initial program 99.8%
+-commutative99.8%
associate-*l*99.7%
fma-define99.7%
sub-neg99.7%
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.7%
Taylor expanded in z around inf 51.1%
if -0.0719999999999999946 < z < -1.9999999999999999e-126 or -8.80000000000000029e-154 < z < 1.75000000000000004e-254 or 2.5000000000000001e-197 < z < 42000Initial program 98.6%
+-commutative98.6%
associate-*l*99.9%
fma-define99.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 62.7%
Taylor expanded in z around 0 60.9%
if -1.9999999999999999e-126 < z < -8.80000000000000029e-154 or 1.75000000000000004e-254 < z < 2.5000000000000001e-197Initial program 99.5%
metadata-eval99.5%
Simplified99.5%
Taylor expanded in x around inf 85.4%
+-commutative85.4%
sub-neg85.4%
distribute-rgt-in85.4%
metadata-eval85.4%
distribute-lft-neg-in85.4%
metadata-eval85.4%
neg-mul-185.4%
distribute-lft-in85.4%
*-commutative85.4%
+-commutative85.4%
mul-1-neg85.4%
*-commutative85.4%
distribute-neg-in85.4%
metadata-eval85.4%
associate-+r+85.4%
metadata-eval85.4%
distribute-rgt-neg-in85.4%
metadata-eval85.4%
Simplified85.4%
Taylor expanded in z around 0 85.4%
*-commutative85.4%
Simplified85.4%
Final simplification57.8%
(FPCore (x y z) :precision binary64 (if (or (<= x -1040000.0) (not (<= x 1.6e+53))) (* x (+ -3.0 (* z 6.0))) (* y (+ 4.0 (* z -6.0)))))
double code(double x, double y, double z) {
double tmp;
if ((x <= -1040000.0) || !(x <= 1.6e+53)) {
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 <= (-1040000.0d0)) .or. (.not. (x <= 1.6d+53))) 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 <= -1040000.0) || !(x <= 1.6e+53)) {
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 <= -1040000.0) or not (x <= 1.6e+53): 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 <= -1040000.0) || !(x <= 1.6e+53)) 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 <= -1040000.0) || ~((x <= 1.6e+53))) 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, -1040000.0], N[Not[LessEqual[x, 1.6e+53]], $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 -1040000 \lor \neg \left(x \leq 1.6 \cdot 10^{+53}\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 < -1.04e6 or 1.6e53 < x Initial program 98.9%
metadata-eval98.9%
Simplified98.9%
Taylor expanded in x around inf 84.0%
+-commutative84.0%
sub-neg84.0%
distribute-rgt-in84.0%
metadata-eval84.0%
distribute-lft-neg-in84.0%
metadata-eval84.0%
neg-mul-184.0%
distribute-lft-in84.0%
*-commutative84.0%
+-commutative84.0%
mul-1-neg84.0%
*-commutative84.0%
distribute-neg-in84.0%
metadata-eval84.0%
associate-+r+84.0%
metadata-eval84.0%
distribute-rgt-neg-in84.0%
metadata-eval84.0%
Simplified84.0%
if -1.04e6 < x < 1.6e53Initial program 99.6%
+-commutative99.6%
associate-*l*99.8%
fma-define99.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 81.3%
Final simplification82.5%
(FPCore (x y z) :precision binary64 (if (<= z -0.6) (* (- y x) (* z -6.0)) (if (<= z 0.5) (+ x (* (- y x) 4.0)) (* -6.0 (* (- y x) z)))))
double code(double x, double y, double z) {
double tmp;
if (z <= -0.6) {
tmp = (y - x) * (z * -6.0);
} else if (z <= 0.5) {
tmp = x + ((y - x) * 4.0);
} else {
tmp = -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.6d0)) then
tmp = (y - x) * (z * (-6.0d0))
else if (z <= 0.5d0) then
tmp = x + ((y - x) * 4.0d0)
else
tmp = (-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.6) {
tmp = (y - x) * (z * -6.0);
} else if (z <= 0.5) {
tmp = x + ((y - x) * 4.0);
} else {
tmp = -6.0 * ((y - x) * z);
}
return tmp;
}
def code(x, y, z): tmp = 0 if z <= -0.6: tmp = (y - x) * (z * -6.0) elif z <= 0.5: tmp = x + ((y - x) * 4.0) else: tmp = -6.0 * ((y - x) * z) return tmp
function code(x, y, z) tmp = 0.0 if (z <= -0.6) tmp = Float64(Float64(y - x) * Float64(z * -6.0)); elseif (z <= 0.5) tmp = Float64(x + Float64(Float64(y - x) * 4.0)); else tmp = 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.6) tmp = (y - x) * (z * -6.0); elseif (z <= 0.5) tmp = x + ((y - x) * 4.0); else tmp = -6.0 * ((y - x) * z); end tmp_2 = tmp; end
code[x_, y_, z_] := If[LessEqual[z, -0.6], N[(N[(y - x), $MachinePrecision] * N[(z * -6.0), $MachinePrecision]), $MachinePrecision], If[LessEqual[z, 0.5], N[(x + N[(N[(y - x), $MachinePrecision] * 4.0), $MachinePrecision]), $MachinePrecision], N[(-6.0 * N[(N[(y - x), $MachinePrecision] * z), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;z \leq -0.6:\\
\;\;\;\;\left(y - x\right) \cdot \left(z \cdot -6\right)\\
\mathbf{elif}\;z \leq 0.5:\\
\;\;\;\;x + \left(y - x\right) \cdot 4\\
\mathbf{else}:\\
\;\;\;\;-6 \cdot \left(\left(y - x\right) \cdot z\right)\\
\end{array}
\end{array}
if z < -0.599999999999999978Initial program 99.8%
metadata-eval99.8%
Simplified99.8%
Taylor expanded in z around inf 96.9%
Taylor expanded in z around inf 96.9%
associate-*r*97.1%
Simplified97.1%
if -0.599999999999999978 < z < 0.5Initial program 98.8%
metadata-eval98.8%
Simplified98.8%
Taylor expanded in z around 0 97.1%
if 0.5 < z Initial program 99.7%
metadata-eval99.7%
Simplified99.7%
Taylor expanded in z around inf 97.0%
Taylor expanded in z around inf 97.0%
Final simplification97.1%
(FPCore (x y z) :precision binary64 (if (or (<= x -7.8e+130) (not (<= x 6.4e+94))) (* x -3.0) (* y 4.0)))
double code(double x, double y, double z) {
double tmp;
if ((x <= -7.8e+130) || !(x <= 6.4e+94)) {
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 <= (-7.8d+130)) .or. (.not. (x <= 6.4d+94))) 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 <= -7.8e+130) || !(x <= 6.4e+94)) {
tmp = x * -3.0;
} else {
tmp = y * 4.0;
}
return tmp;
}
def code(x, y, z): tmp = 0 if (x <= -7.8e+130) or not (x <= 6.4e+94): tmp = x * -3.0 else: tmp = y * 4.0 return tmp
function code(x, y, z) tmp = 0.0 if ((x <= -7.8e+130) || !(x <= 6.4e+94)) 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 <= -7.8e+130) || ~((x <= 6.4e+94))) tmp = x * -3.0; else tmp = y * 4.0; end tmp_2 = tmp; end
code[x_, y_, z_] := If[Or[LessEqual[x, -7.8e+130], N[Not[LessEqual[x, 6.4e+94]], $MachinePrecision]], N[(x * -3.0), $MachinePrecision], N[(y * 4.0), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq -7.8 \cdot 10^{+130} \lor \neg \left(x \leq 6.4 \cdot 10^{+94}\right):\\
\;\;\;\;x \cdot -3\\
\mathbf{else}:\\
\;\;\;\;y \cdot 4\\
\end{array}
\end{array}
if x < -7.8000000000000004e130 or 6.40000000000000028e94 < x Initial program 98.6%
metadata-eval98.6%
Simplified98.6%
Taylor expanded in x around inf 92.3%
+-commutative92.3%
sub-neg92.3%
distribute-rgt-in92.3%
metadata-eval92.3%
distribute-lft-neg-in92.3%
metadata-eval92.3%
neg-mul-192.3%
distribute-lft-in92.3%
*-commutative92.3%
+-commutative92.3%
mul-1-neg92.3%
*-commutative92.3%
distribute-neg-in92.3%
metadata-eval92.3%
associate-+r+92.3%
metadata-eval92.3%
distribute-rgt-neg-in92.3%
metadata-eval92.3%
Simplified92.3%
Taylor expanded in z around 0 49.8%
*-commutative49.8%
Simplified49.8%
if -7.8000000000000004e130 < x < 6.40000000000000028e94Initial program 99.6%
+-commutative99.6%
associate-*l*99.8%
fma-define99.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 72.7%
Taylor expanded in z around 0 37.1%
Final simplification41.1%
(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.3%
metadata-eval99.3%
Simplified99.3%
Final simplification99.3%
(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.3%
metadata-eval99.3%
Simplified99.3%
Taylor expanded in x around inf 48.1%
+-commutative48.1%
sub-neg48.1%
distribute-rgt-in48.1%
metadata-eval48.1%
distribute-lft-neg-in48.1%
metadata-eval48.1%
neg-mul-148.1%
distribute-lft-in48.1%
*-commutative48.1%
+-commutative48.1%
mul-1-neg48.1%
*-commutative48.1%
distribute-neg-in48.1%
metadata-eval48.1%
associate-+r+48.1%
metadata-eval48.1%
distribute-rgt-neg-in48.1%
metadata-eval48.1%
Simplified48.1%
Taylor expanded in z around 0 22.1%
*-commutative22.1%
Simplified22.1%
Final simplification22.1%
(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.3%
metadata-eval99.3%
Simplified99.3%
Taylor expanded in z around inf 51.3%
Taylor expanded in z around 0 2.8%
Final simplification2.8%
herbie shell --seed 2024044
(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))))