
(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 14 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) (fma z -6.0 4.0) x))
double code(double x, double y, double z) {
return fma((y - x), fma(z, -6.0, 4.0), x);
}
function code(x, y, z) return fma(Float64(y - x), fma(z, -6.0, 4.0), x) end
code[x_, y_, z_] := N[(N[(y - x), $MachinePrecision] * N[(z * -6.0 + 4.0), $MachinePrecision] + x), $MachinePrecision]
\begin{array}{l}
\\
\mathsf{fma}\left(y - x, \mathsf{fma}\left(z, -6, 4\right), x\right)
\end{array}
Initial program 99.5%
+-commutative99.5%
associate-*l*99.8%
fma-define99.8%
sub-neg99.8%
+-commutative99.8%
distribute-lft-in99.8%
distribute-rgt-neg-out99.8%
*-commutative99.8%
distribute-rgt-neg-in99.8%
fma-define99.8%
metadata-eval99.8%
metadata-eval99.8%
metadata-eval99.8%
Simplified99.8%
(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.5%
+-commutative99.5%
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%
(FPCore (x y z)
:precision binary64
(let* ((t_0 (+ x (* y 4.0))) (t_1 (* y (* z -6.0))))
(if (<= z -5.5e+182)
(* z (* x 6.0))
(if (<= z -0.3)
t_1
(if (<= z -4.5e-167)
t_0
(if (<= z -1.25e-193)
(* x -3.0)
(if (<= z 8.2e-259)
t_0
(if (<= z 5.7e-126)
(* x -3.0)
(if (<= z 59000000000000.0)
t_0
(if (<= z 7.6e+59)
t_1
(if (<= z 3.4e+92)
(* 6.0 (* x z))
(if (<= z 3.1e+177)
t_1
(if (<= z 2.5e+221)
(* x (* z 6.0))
(* -6.0 (* y z)))))))))))))))
double code(double x, double y, double z) {
double t_0 = x + (y * 4.0);
double t_1 = y * (z * -6.0);
double tmp;
if (z <= -5.5e+182) {
tmp = z * (x * 6.0);
} else if (z <= -0.3) {
tmp = t_1;
} else if (z <= -4.5e-167) {
tmp = t_0;
} else if (z <= -1.25e-193) {
tmp = x * -3.0;
} else if (z <= 8.2e-259) {
tmp = t_0;
} else if (z <= 5.7e-126) {
tmp = x * -3.0;
} else if (z <= 59000000000000.0) {
tmp = t_0;
} else if (z <= 7.6e+59) {
tmp = t_1;
} else if (z <= 3.4e+92) {
tmp = 6.0 * (x * z);
} else if (z <= 3.1e+177) {
tmp = t_1;
} else if (z <= 2.5e+221) {
tmp = x * (z * 6.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 = x + (y * 4.0d0)
t_1 = y * (z * (-6.0d0))
if (z <= (-5.5d+182)) then
tmp = z * (x * 6.0d0)
else if (z <= (-0.3d0)) then
tmp = t_1
else if (z <= (-4.5d-167)) then
tmp = t_0
else if (z <= (-1.25d-193)) then
tmp = x * (-3.0d0)
else if (z <= 8.2d-259) then
tmp = t_0
else if (z <= 5.7d-126) then
tmp = x * (-3.0d0)
else if (z <= 59000000000000.0d0) then
tmp = t_0
else if (z <= 7.6d+59) then
tmp = t_1
else if (z <= 3.4d+92) then
tmp = 6.0d0 * (x * z)
else if (z <= 3.1d+177) then
tmp = t_1
else if (z <= 2.5d+221) then
tmp = x * (z * 6.0d0)
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 = x + (y * 4.0);
double t_1 = y * (z * -6.0);
double tmp;
if (z <= -5.5e+182) {
tmp = z * (x * 6.0);
} else if (z <= -0.3) {
tmp = t_1;
} else if (z <= -4.5e-167) {
tmp = t_0;
} else if (z <= -1.25e-193) {
tmp = x * -3.0;
} else if (z <= 8.2e-259) {
tmp = t_0;
} else if (z <= 5.7e-126) {
tmp = x * -3.0;
} else if (z <= 59000000000000.0) {
tmp = t_0;
} else if (z <= 7.6e+59) {
tmp = t_1;
} else if (z <= 3.4e+92) {
tmp = 6.0 * (x * z);
} else if (z <= 3.1e+177) {
tmp = t_1;
} else if (z <= 2.5e+221) {
tmp = x * (z * 6.0);
} else {
tmp = -6.0 * (y * z);
}
return tmp;
}
def code(x, y, z): t_0 = x + (y * 4.0) t_1 = y * (z * -6.0) tmp = 0 if z <= -5.5e+182: tmp = z * (x * 6.0) elif z <= -0.3: tmp = t_1 elif z <= -4.5e-167: tmp = t_0 elif z <= -1.25e-193: tmp = x * -3.0 elif z <= 8.2e-259: tmp = t_0 elif z <= 5.7e-126: tmp = x * -3.0 elif z <= 59000000000000.0: tmp = t_0 elif z <= 7.6e+59: tmp = t_1 elif z <= 3.4e+92: tmp = 6.0 * (x * z) elif z <= 3.1e+177: tmp = t_1 elif z <= 2.5e+221: tmp = x * (z * 6.0) else: tmp = -6.0 * (y * z) return tmp
function code(x, y, z) t_0 = Float64(x + Float64(y * 4.0)) t_1 = Float64(y * Float64(z * -6.0)) tmp = 0.0 if (z <= -5.5e+182) tmp = Float64(z * Float64(x * 6.0)); elseif (z <= -0.3) tmp = t_1; elseif (z <= -4.5e-167) tmp = t_0; elseif (z <= -1.25e-193) tmp = Float64(x * -3.0); elseif (z <= 8.2e-259) tmp = t_0; elseif (z <= 5.7e-126) tmp = Float64(x * -3.0); elseif (z <= 59000000000000.0) tmp = t_0; elseif (z <= 7.6e+59) tmp = t_1; elseif (z <= 3.4e+92) tmp = Float64(6.0 * Float64(x * z)); elseif (z <= 3.1e+177) tmp = t_1; elseif (z <= 2.5e+221) tmp = Float64(x * Float64(z * 6.0)); else tmp = Float64(-6.0 * Float64(y * z)); end return tmp end
function tmp_2 = code(x, y, z) t_0 = x + (y * 4.0); t_1 = y * (z * -6.0); tmp = 0.0; if (z <= -5.5e+182) tmp = z * (x * 6.0); elseif (z <= -0.3) tmp = t_1; elseif (z <= -4.5e-167) tmp = t_0; elseif (z <= -1.25e-193) tmp = x * -3.0; elseif (z <= 8.2e-259) tmp = t_0; elseif (z <= 5.7e-126) tmp = x * -3.0; elseif (z <= 59000000000000.0) tmp = t_0; elseif (z <= 7.6e+59) tmp = t_1; elseif (z <= 3.4e+92) tmp = 6.0 * (x * z); elseif (z <= 3.1e+177) tmp = t_1; elseif (z <= 2.5e+221) tmp = x * (z * 6.0); else tmp = -6.0 * (y * z); end tmp_2 = tmp; end
code[x_, y_, z_] := Block[{t$95$0 = N[(x + N[(y * 4.0), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(y * N[(z * -6.0), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[z, -5.5e+182], N[(z * N[(x * 6.0), $MachinePrecision]), $MachinePrecision], If[LessEqual[z, -0.3], t$95$1, If[LessEqual[z, -4.5e-167], t$95$0, If[LessEqual[z, -1.25e-193], N[(x * -3.0), $MachinePrecision], If[LessEqual[z, 8.2e-259], t$95$0, If[LessEqual[z, 5.7e-126], N[(x * -3.0), $MachinePrecision], If[LessEqual[z, 59000000000000.0], t$95$0, If[LessEqual[z, 7.6e+59], t$95$1, If[LessEqual[z, 3.4e+92], N[(6.0 * N[(x * z), $MachinePrecision]), $MachinePrecision], If[LessEqual[z, 3.1e+177], t$95$1, If[LessEqual[z, 2.5e+221], N[(x * N[(z * 6.0), $MachinePrecision]), $MachinePrecision], N[(-6.0 * N[(y * z), $MachinePrecision]), $MachinePrecision]]]]]]]]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := x + y \cdot 4\\
t_1 := y \cdot \left(z \cdot -6\right)\\
\mathbf{if}\;z \leq -5.5 \cdot 10^{+182}:\\
\;\;\;\;z \cdot \left(x \cdot 6\right)\\
\mathbf{elif}\;z \leq -0.3:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;z \leq -4.5 \cdot 10^{-167}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;z \leq -1.25 \cdot 10^{-193}:\\
\;\;\;\;x \cdot -3\\
\mathbf{elif}\;z \leq 8.2 \cdot 10^{-259}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;z \leq 5.7 \cdot 10^{-126}:\\
\;\;\;\;x \cdot -3\\
\mathbf{elif}\;z \leq 59000000000000:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;z \leq 7.6 \cdot 10^{+59}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;z \leq 3.4 \cdot 10^{+92}:\\
\;\;\;\;6 \cdot \left(x \cdot z\right)\\
\mathbf{elif}\;z \leq 3.1 \cdot 10^{+177}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;z \leq 2.5 \cdot 10^{+221}:\\
\;\;\;\;x \cdot \left(z \cdot 6\right)\\
\mathbf{else}:\\
\;\;\;\;-6 \cdot \left(y \cdot z\right)\\
\end{array}
\end{array}
if z < -5.49999999999999977e182Initial program 99.8%
metadata-eval99.8%
Simplified99.8%
Taylor expanded in z around 0 99.6%
Taylor expanded in z around inf 99.6%
associate-*r*99.8%
metadata-eval99.8%
distribute-lft-neg-in99.8%
distribute-lft-neg-in99.8%
*-commutative99.8%
associate-*r*99.8%
rem-cube-cbrt98.8%
rem-cube-cbrt99.8%
*-commutative99.8%
distribute-rgt-neg-in99.8%
*-commutative99.8%
rem-cube-cbrt98.8%
mul-1-neg98.8%
rem-cube-cbrt99.8%
associate-*r*99.8%
metadata-eval99.8%
Simplified99.8%
Taylor expanded in y around 0 66.9%
*-commutative66.9%
*-commutative66.9%
associate-*r*67.0%
Simplified67.0%
if -5.49999999999999977e182 < z < -0.299999999999999989 or 5.9e13 < z < 7.6000000000000002e59 or 3.3999999999999998e92 < z < 3.0999999999999999e177Initial program 99.6%
metadata-eval99.6%
Simplified99.6%
Taylor expanded in z around 0 99.6%
Taylor expanded in z around inf 96.7%
associate-*r*96.9%
metadata-eval96.9%
distribute-lft-neg-in96.9%
distribute-lft-neg-in96.9%
*-commutative96.9%
associate-*r*96.8%
rem-cube-cbrt95.5%
rem-cube-cbrt96.8%
*-commutative96.8%
distribute-rgt-neg-in96.8%
*-commutative96.8%
rem-cube-cbrt95.5%
mul-1-neg95.5%
rem-cube-cbrt96.8%
associate-*r*96.8%
metadata-eval96.8%
Simplified96.8%
Taylor expanded in y around inf 63.1%
*-commutative63.1%
associate-*r*63.3%
*-commutative63.3%
Simplified63.3%
if -0.299999999999999989 < z < -4.5000000000000001e-167 or -1.2500000000000001e-193 < z < 8.1999999999999996e-259 or 5.69999999999999979e-126 < z < 5.9e13Initial program 99.4%
metadata-eval99.4%
Simplified99.4%
Taylor expanded in y around inf 65.4%
Taylor expanded in z around 0 64.4%
if -4.5000000000000001e-167 < z < -1.2500000000000001e-193 or 8.1999999999999996e-259 < z < 5.69999999999999979e-126Initial program 99.4%
metadata-eval99.4%
Simplified99.4%
Taylor expanded in x around inf 70.7%
sub-neg70.7%
distribute-rgt-in70.7%
metadata-eval70.7%
distribute-lft-neg-in70.7%
associate-+r+70.7%
metadata-eval70.7%
distribute-rgt-neg-in70.7%
metadata-eval70.7%
Simplified70.7%
Taylor expanded in z around 0 70.7%
*-commutative70.7%
Simplified70.7%
if 7.6000000000000002e59 < z < 3.3999999999999998e92Initial program 99.4%
metadata-eval99.4%
Simplified99.4%
Taylor expanded in z around 0 99.7%
Taylor expanded in z around inf 99.7%
associate-*r*99.4%
metadata-eval99.4%
distribute-lft-neg-in99.4%
distribute-lft-neg-in99.4%
*-commutative99.4%
associate-*r*99.4%
rem-cube-cbrt98.0%
rem-cube-cbrt99.4%
*-commutative99.4%
distribute-rgt-neg-in99.4%
*-commutative99.4%
rem-cube-cbrt98.0%
mul-1-neg98.0%
rem-cube-cbrt99.4%
associate-*r*99.4%
metadata-eval99.4%
Simplified99.4%
Taylor expanded in y around 0 99.7%
if 3.0999999999999999e177 < z < 2.5000000000000001e221Initial program 99.8%
metadata-eval99.8%
Simplified99.8%
Taylor expanded in x around inf 81.0%
sub-neg81.0%
distribute-rgt-in81.0%
metadata-eval81.0%
distribute-lft-neg-in81.0%
associate-+r+81.0%
metadata-eval81.0%
distribute-rgt-neg-in81.0%
metadata-eval81.0%
Simplified81.0%
Taylor expanded in z around inf 81.0%
if 2.5000000000000001e221 < z Initial program 99.9%
metadata-eval99.9%
Simplified99.9%
Taylor expanded in z around 0 99.9%
Taylor expanded in z around inf 99.9%
associate-*r*99.8%
metadata-eval99.8%
distribute-lft-neg-in99.8%
distribute-lft-neg-in99.8%
*-commutative99.8%
associate-*r*99.9%
rem-cube-cbrt99.3%
rem-cube-cbrt99.9%
*-commutative99.9%
distribute-rgt-neg-in99.9%
*-commutative99.9%
rem-cube-cbrt99.3%
mul-1-neg99.3%
rem-cube-cbrt99.9%
associate-*r*99.9%
metadata-eval99.9%
Simplified99.9%
Taylor expanded in y around inf 64.4%
Final simplification67.0%
(FPCore (x y z)
:precision binary64
(let* ((t_0 (* y (* z -6.0))))
(if (<= z -3.3e+180)
(* z (* x 6.0))
(if (<= z -8100000000.0)
t_0
(if (<= z 0.66)
(* x -3.0)
(if (<= z 7.5e+59)
t_0
(if (<= z 1.4e+92)
(* 6.0 (* x z))
(if (<= z 3.2e+177)
t_0
(if (<= z 1.9e+221) (* x (* z 6.0)) (* -6.0 (* y z)))))))))))
double code(double x, double y, double z) {
double t_0 = y * (z * -6.0);
double tmp;
if (z <= -3.3e+180) {
tmp = z * (x * 6.0);
} else if (z <= -8100000000.0) {
tmp = t_0;
} else if (z <= 0.66) {
tmp = x * -3.0;
} else if (z <= 7.5e+59) {
tmp = t_0;
} else if (z <= 1.4e+92) {
tmp = 6.0 * (x * z);
} else if (z <= 3.2e+177) {
tmp = t_0;
} else if (z <= 1.9e+221) {
tmp = x * (z * 6.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) :: tmp
t_0 = y * (z * (-6.0d0))
if (z <= (-3.3d+180)) then
tmp = z * (x * 6.0d0)
else if (z <= (-8100000000.0d0)) then
tmp = t_0
else if (z <= 0.66d0) then
tmp = x * (-3.0d0)
else if (z <= 7.5d+59) then
tmp = t_0
else if (z <= 1.4d+92) then
tmp = 6.0d0 * (x * z)
else if (z <= 3.2d+177) then
tmp = t_0
else if (z <= 1.9d+221) then
tmp = x * (z * 6.0d0)
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 <= -3.3e+180) {
tmp = z * (x * 6.0);
} else if (z <= -8100000000.0) {
tmp = t_0;
} else if (z <= 0.66) {
tmp = x * -3.0;
} else if (z <= 7.5e+59) {
tmp = t_0;
} else if (z <= 1.4e+92) {
tmp = 6.0 * (x * z);
} else if (z <= 3.2e+177) {
tmp = t_0;
} else if (z <= 1.9e+221) {
tmp = x * (z * 6.0);
} else {
tmp = -6.0 * (y * z);
}
return tmp;
}
def code(x, y, z): t_0 = y * (z * -6.0) tmp = 0 if z <= -3.3e+180: tmp = z * (x * 6.0) elif z <= -8100000000.0: tmp = t_0 elif z <= 0.66: tmp = x * -3.0 elif z <= 7.5e+59: tmp = t_0 elif z <= 1.4e+92: tmp = 6.0 * (x * z) elif z <= 3.2e+177: tmp = t_0 elif z <= 1.9e+221: tmp = x * (z * 6.0) 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 <= -3.3e+180) tmp = Float64(z * Float64(x * 6.0)); elseif (z <= -8100000000.0) tmp = t_0; elseif (z <= 0.66) tmp = Float64(x * -3.0); elseif (z <= 7.5e+59) tmp = t_0; elseif (z <= 1.4e+92) tmp = Float64(6.0 * Float64(x * z)); elseif (z <= 3.2e+177) tmp = t_0; elseif (z <= 1.9e+221) tmp = Float64(x * Float64(z * 6.0)); 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 <= -3.3e+180) tmp = z * (x * 6.0); elseif (z <= -8100000000.0) tmp = t_0; elseif (z <= 0.66) tmp = x * -3.0; elseif (z <= 7.5e+59) tmp = t_0; elseif (z <= 1.4e+92) tmp = 6.0 * (x * z); elseif (z <= 3.2e+177) tmp = t_0; elseif (z <= 1.9e+221) tmp = x * (z * 6.0); 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, -3.3e+180], N[(z * N[(x * 6.0), $MachinePrecision]), $MachinePrecision], If[LessEqual[z, -8100000000.0], t$95$0, If[LessEqual[z, 0.66], N[(x * -3.0), $MachinePrecision], If[LessEqual[z, 7.5e+59], t$95$0, If[LessEqual[z, 1.4e+92], N[(6.0 * N[(x * z), $MachinePrecision]), $MachinePrecision], If[LessEqual[z, 3.2e+177], t$95$0, If[LessEqual[z, 1.9e+221], N[(x * N[(z * 6.0), $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 -3.3 \cdot 10^{+180}:\\
\;\;\;\;z \cdot \left(x \cdot 6\right)\\
\mathbf{elif}\;z \leq -8100000000:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;z \leq 0.66:\\
\;\;\;\;x \cdot -3\\
\mathbf{elif}\;z \leq 7.5 \cdot 10^{+59}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;z \leq 1.4 \cdot 10^{+92}:\\
\;\;\;\;6 \cdot \left(x \cdot z\right)\\
\mathbf{elif}\;z \leq 3.2 \cdot 10^{+177}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;z \leq 1.9 \cdot 10^{+221}:\\
\;\;\;\;x \cdot \left(z \cdot 6\right)\\
\mathbf{else}:\\
\;\;\;\;-6 \cdot \left(y \cdot z\right)\\
\end{array}
\end{array}
if z < -3.29999999999999989e180Initial program 99.8%
metadata-eval99.8%
Simplified99.8%
Taylor expanded in z around 0 99.6%
Taylor expanded in z around inf 99.6%
associate-*r*99.8%
metadata-eval99.8%
distribute-lft-neg-in99.8%
distribute-lft-neg-in99.8%
*-commutative99.8%
associate-*r*99.8%
rem-cube-cbrt98.8%
rem-cube-cbrt99.8%
*-commutative99.8%
distribute-rgt-neg-in99.8%
*-commutative99.8%
rem-cube-cbrt98.8%
mul-1-neg98.8%
rem-cube-cbrt99.8%
associate-*r*99.8%
metadata-eval99.8%
Simplified99.8%
Taylor expanded in y around 0 66.9%
*-commutative66.9%
*-commutative66.9%
associate-*r*67.0%
Simplified67.0%
if -3.29999999999999989e180 < z < -8.1e9 or 0.660000000000000031 < z < 7.4999999999999996e59 or 1.4e92 < z < 3.2e177Initial program 99.6%
metadata-eval99.6%
Simplified99.6%
Taylor expanded in z around 0 99.6%
Taylor expanded in z around inf 96.3%
associate-*r*96.6%
metadata-eval96.6%
distribute-lft-neg-in96.6%
distribute-lft-neg-in96.6%
*-commutative96.6%
associate-*r*96.4%
rem-cube-cbrt95.1%
rem-cube-cbrt96.4%
*-commutative96.4%
distribute-rgt-neg-in96.4%
*-commutative96.4%
rem-cube-cbrt95.1%
mul-1-neg95.1%
rem-cube-cbrt96.4%
associate-*r*96.4%
metadata-eval96.4%
Simplified96.4%
Taylor expanded in y around inf 62.4%
*-commutative62.4%
associate-*r*62.6%
*-commutative62.6%
Simplified62.6%
if -8.1e9 < z < 0.660000000000000031Initial program 99.4%
metadata-eval99.4%
Simplified99.4%
Taylor expanded in x around inf 43.7%
sub-neg43.7%
distribute-rgt-in43.7%
metadata-eval43.7%
distribute-lft-neg-in43.7%
associate-+r+43.7%
metadata-eval43.7%
distribute-rgt-neg-in43.7%
metadata-eval43.7%
Simplified43.7%
Taylor expanded in z around 0 41.6%
*-commutative41.6%
Simplified41.6%
if 7.4999999999999996e59 < z < 1.4e92Initial program 99.4%
metadata-eval99.4%
Simplified99.4%
Taylor expanded in z around 0 99.7%
Taylor expanded in z around inf 99.7%
associate-*r*99.4%
metadata-eval99.4%
distribute-lft-neg-in99.4%
distribute-lft-neg-in99.4%
*-commutative99.4%
associate-*r*99.4%
rem-cube-cbrt98.0%
rem-cube-cbrt99.4%
*-commutative99.4%
distribute-rgt-neg-in99.4%
*-commutative99.4%
rem-cube-cbrt98.0%
mul-1-neg98.0%
rem-cube-cbrt99.4%
associate-*r*99.4%
metadata-eval99.4%
Simplified99.4%
Taylor expanded in y around 0 99.7%
if 3.2e177 < z < 1.90000000000000017e221Initial program 99.8%
metadata-eval99.8%
Simplified99.8%
Taylor expanded in x around inf 81.0%
sub-neg81.0%
distribute-rgt-in81.0%
metadata-eval81.0%
distribute-lft-neg-in81.0%
associate-+r+81.0%
metadata-eval81.0%
distribute-rgt-neg-in81.0%
metadata-eval81.0%
Simplified81.0%
Taylor expanded in z around inf 81.0%
if 1.90000000000000017e221 < z Initial program 99.9%
metadata-eval99.9%
Simplified99.9%
Taylor expanded in z around 0 99.9%
Taylor expanded in z around inf 99.9%
associate-*r*99.8%
metadata-eval99.8%
distribute-lft-neg-in99.8%
distribute-lft-neg-in99.8%
*-commutative99.8%
associate-*r*99.9%
rem-cube-cbrt99.3%
rem-cube-cbrt99.9%
*-commutative99.9%
distribute-rgt-neg-in99.9%
*-commutative99.9%
rem-cube-cbrt99.3%
mul-1-neg99.3%
rem-cube-cbrt99.9%
associate-*r*99.9%
metadata-eval99.9%
Simplified99.9%
Taylor expanded in y around inf 64.4%
Final simplification55.0%
(FPCore (x y z)
:precision binary64
(let* ((t_0 (* x (* z 6.0))) (t_1 (* y (* z -6.0))))
(if (<= z -2e+183)
t_0
(if (<= z -8100000000.0)
t_1
(if (<= z 0.65)
(* x -3.0)
(if (<= z 1e+60)
t_1
(if (<= z 1.75e+92)
(* 6.0 (* x z))
(if (<= z 2e+178)
t_1
(if (<= z 1.15e+224) t_0 (* -6.0 (* y z)))))))))))
double code(double x, double y, double z) {
double t_0 = x * (z * 6.0);
double t_1 = y * (z * -6.0);
double tmp;
if (z <= -2e+183) {
tmp = t_0;
} else if (z <= -8100000000.0) {
tmp = t_1;
} else if (z <= 0.65) {
tmp = x * -3.0;
} else if (z <= 1e+60) {
tmp = t_1;
} else if (z <= 1.75e+92) {
tmp = 6.0 * (x * z);
} else if (z <= 2e+178) {
tmp = t_1;
} else if (z <= 1.15e+224) {
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 = x * (z * 6.0d0)
t_1 = y * (z * (-6.0d0))
if (z <= (-2d+183)) then
tmp = t_0
else if (z <= (-8100000000.0d0)) then
tmp = t_1
else if (z <= 0.65d0) then
tmp = x * (-3.0d0)
else if (z <= 1d+60) then
tmp = t_1
else if (z <= 1.75d+92) then
tmp = 6.0d0 * (x * z)
else if (z <= 2d+178) then
tmp = t_1
else if (z <= 1.15d+224) 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 = x * (z * 6.0);
double t_1 = y * (z * -6.0);
double tmp;
if (z <= -2e+183) {
tmp = t_0;
} else if (z <= -8100000000.0) {
tmp = t_1;
} else if (z <= 0.65) {
tmp = x * -3.0;
} else if (z <= 1e+60) {
tmp = t_1;
} else if (z <= 1.75e+92) {
tmp = 6.0 * (x * z);
} else if (z <= 2e+178) {
tmp = t_1;
} else if (z <= 1.15e+224) {
tmp = t_0;
} else {
tmp = -6.0 * (y * z);
}
return tmp;
}
def code(x, y, z): t_0 = x * (z * 6.0) t_1 = y * (z * -6.0) tmp = 0 if z <= -2e+183: tmp = t_0 elif z <= -8100000000.0: tmp = t_1 elif z <= 0.65: tmp = x * -3.0 elif z <= 1e+60: tmp = t_1 elif z <= 1.75e+92: tmp = 6.0 * (x * z) elif z <= 2e+178: tmp = t_1 elif z <= 1.15e+224: tmp = t_0 else: tmp = -6.0 * (y * z) return tmp
function code(x, y, z) t_0 = Float64(x * Float64(z * 6.0)) t_1 = Float64(y * Float64(z * -6.0)) tmp = 0.0 if (z <= -2e+183) tmp = t_0; elseif (z <= -8100000000.0) tmp = t_1; elseif (z <= 0.65) tmp = Float64(x * -3.0); elseif (z <= 1e+60) tmp = t_1; elseif (z <= 1.75e+92) tmp = Float64(6.0 * Float64(x * z)); elseif (z <= 2e+178) tmp = t_1; elseif (z <= 1.15e+224) tmp = t_0; else tmp = Float64(-6.0 * Float64(y * z)); end return tmp end
function tmp_2 = code(x, y, z) t_0 = x * (z * 6.0); t_1 = y * (z * -6.0); tmp = 0.0; if (z <= -2e+183) tmp = t_0; elseif (z <= -8100000000.0) tmp = t_1; elseif (z <= 0.65) tmp = x * -3.0; elseif (z <= 1e+60) tmp = t_1; elseif (z <= 1.75e+92) tmp = 6.0 * (x * z); elseif (z <= 2e+178) tmp = t_1; elseif (z <= 1.15e+224) tmp = t_0; else tmp = -6.0 * (y * z); end tmp_2 = tmp; end
code[x_, y_, z_] := Block[{t$95$0 = N[(x * N[(z * 6.0), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(y * N[(z * -6.0), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[z, -2e+183], t$95$0, If[LessEqual[z, -8100000000.0], t$95$1, If[LessEqual[z, 0.65], N[(x * -3.0), $MachinePrecision], If[LessEqual[z, 1e+60], t$95$1, If[LessEqual[z, 1.75e+92], N[(6.0 * N[(x * z), $MachinePrecision]), $MachinePrecision], If[LessEqual[z, 2e+178], t$95$1, If[LessEqual[z, 1.15e+224], t$95$0, N[(-6.0 * N[(y * z), $MachinePrecision]), $MachinePrecision]]]]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := x \cdot \left(z \cdot 6\right)\\
t_1 := y \cdot \left(z \cdot -6\right)\\
\mathbf{if}\;z \leq -2 \cdot 10^{+183}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;z \leq -8100000000:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;z \leq 0.65:\\
\;\;\;\;x \cdot -3\\
\mathbf{elif}\;z \leq 10^{+60}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;z \leq 1.75 \cdot 10^{+92}:\\
\;\;\;\;6 \cdot \left(x \cdot z\right)\\
\mathbf{elif}\;z \leq 2 \cdot 10^{+178}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;z \leq 1.15 \cdot 10^{+224}:\\
\;\;\;\;t\_0\\
\mathbf{else}:\\
\;\;\;\;-6 \cdot \left(y \cdot z\right)\\
\end{array}
\end{array}
if z < -1.99999999999999989e183 or 2.0000000000000001e178 < z < 1.1500000000000001e224Initial program 99.8%
metadata-eval99.8%
Simplified99.8%
Taylor expanded in x around inf 72.4%
sub-neg72.4%
distribute-rgt-in72.4%
metadata-eval72.4%
distribute-lft-neg-in72.4%
associate-+r+72.4%
metadata-eval72.4%
distribute-rgt-neg-in72.4%
metadata-eval72.4%
Simplified72.4%
Taylor expanded in z around inf 72.4%
if -1.99999999999999989e183 < z < -8.1e9 or 0.650000000000000022 < z < 9.9999999999999995e59 or 1.74999999999999993e92 < z < 2.0000000000000001e178Initial program 99.6%
metadata-eval99.6%
Simplified99.6%
Taylor expanded in z around 0 99.6%
Taylor expanded in z around inf 96.3%
associate-*r*96.6%
metadata-eval96.6%
distribute-lft-neg-in96.6%
distribute-lft-neg-in96.6%
*-commutative96.6%
associate-*r*96.4%
rem-cube-cbrt95.1%
rem-cube-cbrt96.4%
*-commutative96.4%
distribute-rgt-neg-in96.4%
*-commutative96.4%
rem-cube-cbrt95.1%
mul-1-neg95.1%
rem-cube-cbrt96.4%
associate-*r*96.4%
metadata-eval96.4%
Simplified96.4%
Taylor expanded in y around inf 62.4%
*-commutative62.4%
associate-*r*62.6%
*-commutative62.6%
Simplified62.6%
if -8.1e9 < z < 0.650000000000000022Initial program 99.4%
metadata-eval99.4%
Simplified99.4%
Taylor expanded in x around inf 43.7%
sub-neg43.7%
distribute-rgt-in43.7%
metadata-eval43.7%
distribute-lft-neg-in43.7%
associate-+r+43.7%
metadata-eval43.7%
distribute-rgt-neg-in43.7%
metadata-eval43.7%
Simplified43.7%
Taylor expanded in z around 0 41.6%
*-commutative41.6%
Simplified41.6%
if 9.9999999999999995e59 < z < 1.74999999999999993e92Initial program 99.4%
metadata-eval99.4%
Simplified99.4%
Taylor expanded in z around 0 99.7%
Taylor expanded in z around inf 99.7%
associate-*r*99.4%
metadata-eval99.4%
distribute-lft-neg-in99.4%
distribute-lft-neg-in99.4%
*-commutative99.4%
associate-*r*99.4%
rem-cube-cbrt98.0%
rem-cube-cbrt99.4%
*-commutative99.4%
distribute-rgt-neg-in99.4%
*-commutative99.4%
rem-cube-cbrt98.0%
mul-1-neg98.0%
rem-cube-cbrt99.4%
associate-*r*99.4%
metadata-eval99.4%
Simplified99.4%
Taylor expanded in y around 0 99.7%
if 1.1500000000000001e224 < z Initial program 99.9%
metadata-eval99.9%
Simplified99.9%
Taylor expanded in z around 0 99.9%
Taylor expanded in z around inf 99.9%
associate-*r*99.8%
metadata-eval99.8%
distribute-lft-neg-in99.8%
distribute-lft-neg-in99.8%
*-commutative99.8%
associate-*r*99.9%
rem-cube-cbrt99.3%
rem-cube-cbrt99.9%
*-commutative99.9%
distribute-rgt-neg-in99.9%
*-commutative99.9%
rem-cube-cbrt99.3%
mul-1-neg99.3%
rem-cube-cbrt99.9%
associate-*r*99.9%
metadata-eval99.9%
Simplified99.9%
Taylor expanded in y around inf 64.4%
Final simplification55.0%
(FPCore (x y z)
:precision binary64
(let* ((t_0 (+ x (* y 4.0))) (t_1 (* -6.0 (* (- y x) z))))
(if (<= z -1.65e-5)
t_1
(if (<= z 1.2e-256)
t_0
(if (<= z 1.7e-128) (* x -3.0) (if (<= z 5.8e-16) t_0 t_1))))))
double code(double x, double y, double z) {
double t_0 = x + (y * 4.0);
double t_1 = -6.0 * ((y - x) * z);
double tmp;
if (z <= -1.65e-5) {
tmp = t_1;
} else if (z <= 1.2e-256) {
tmp = t_0;
} else if (z <= 1.7e-128) {
tmp = x * -3.0;
} else if (z <= 5.8e-16) {
tmp = t_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 = x + (y * 4.0d0)
t_1 = (-6.0d0) * ((y - x) * z)
if (z <= (-1.65d-5)) then
tmp = t_1
else if (z <= 1.2d-256) then
tmp = t_0
else if (z <= 1.7d-128) then
tmp = x * (-3.0d0)
else if (z <= 5.8d-16) then
tmp = t_0
else
tmp = t_1
end if
code = tmp
end function
public static double code(double x, double y, double z) {
double t_0 = x + (y * 4.0);
double t_1 = -6.0 * ((y - x) * z);
double tmp;
if (z <= -1.65e-5) {
tmp = t_1;
} else if (z <= 1.2e-256) {
tmp = t_0;
} else if (z <= 1.7e-128) {
tmp = x * -3.0;
} else if (z <= 5.8e-16) {
tmp = t_0;
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z): t_0 = x + (y * 4.0) t_1 = -6.0 * ((y - x) * z) tmp = 0 if z <= -1.65e-5: tmp = t_1 elif z <= 1.2e-256: tmp = t_0 elif z <= 1.7e-128: tmp = x * -3.0 elif z <= 5.8e-16: tmp = t_0 else: tmp = t_1 return tmp
function code(x, y, z) t_0 = Float64(x + Float64(y * 4.0)) t_1 = Float64(-6.0 * Float64(Float64(y - x) * z)) tmp = 0.0 if (z <= -1.65e-5) tmp = t_1; elseif (z <= 1.2e-256) tmp = t_0; elseif (z <= 1.7e-128) tmp = Float64(x * -3.0); elseif (z <= 5.8e-16) tmp = t_0; else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z) t_0 = x + (y * 4.0); t_1 = -6.0 * ((y - x) * z); tmp = 0.0; if (z <= -1.65e-5) tmp = t_1; elseif (z <= 1.2e-256) tmp = t_0; elseif (z <= 1.7e-128) tmp = x * -3.0; elseif (z <= 5.8e-16) tmp = t_0; else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_] := Block[{t$95$0 = N[(x + N[(y * 4.0), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(-6.0 * N[(N[(y - x), $MachinePrecision] * z), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[z, -1.65e-5], t$95$1, If[LessEqual[z, 1.2e-256], t$95$0, If[LessEqual[z, 1.7e-128], N[(x * -3.0), $MachinePrecision], If[LessEqual[z, 5.8e-16], t$95$0, t$95$1]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := x + y \cdot 4\\
t_1 := -6 \cdot \left(\left(y - x\right) \cdot z\right)\\
\mathbf{if}\;z \leq -1.65 \cdot 10^{-5}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;z \leq 1.2 \cdot 10^{-256}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;z \leq 1.7 \cdot 10^{-128}:\\
\;\;\;\;x \cdot -3\\
\mathbf{elif}\;z \leq 5.8 \cdot 10^{-16}:\\
\;\;\;\;t\_0\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if z < -1.6500000000000001e-5 or 5.7999999999999996e-16 < z Initial program 99.7%
metadata-eval99.7%
Simplified99.7%
Taylor expanded in z around 0 99.7%
Taylor expanded in z around inf 96.4%
if -1.6500000000000001e-5 < z < 1.2e-256 or 1.69999999999999987e-128 < z < 5.7999999999999996e-16Initial program 99.4%
metadata-eval99.4%
Simplified99.4%
Taylor expanded in y around inf 64.0%
Taylor expanded in z around 0 64.0%
if 1.2e-256 < z < 1.69999999999999987e-128Initial program 99.3%
metadata-eval99.3%
Simplified99.3%
Taylor expanded in x around inf 67.9%
sub-neg67.9%
distribute-rgt-in67.9%
metadata-eval67.9%
distribute-lft-neg-in67.9%
associate-+r+67.9%
metadata-eval67.9%
distribute-rgt-neg-in67.9%
metadata-eval67.9%
Simplified67.9%
Taylor expanded in z around 0 67.9%
*-commutative67.9%
Simplified67.9%
Final simplification81.6%
(FPCore (x y z)
:precision binary64
(let* ((t_0 (* x (* z 6.0))) (t_1 (* -6.0 (* y z))))
(if (<= z -1.18e+173)
t_0
(if (<= z -8100000000.0)
t_1
(if (<= z 12.0) (* x -3.0) (if (<= z 2.3e+222) t_0 t_1))))))
double code(double x, double y, double z) {
double t_0 = x * (z * 6.0);
double t_1 = -6.0 * (y * z);
double tmp;
if (z <= -1.18e+173) {
tmp = t_0;
} else if (z <= -8100000000.0) {
tmp = t_1;
} else if (z <= 12.0) {
tmp = x * -3.0;
} else if (z <= 2.3e+222) {
tmp = t_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 = x * (z * 6.0d0)
t_1 = (-6.0d0) * (y * z)
if (z <= (-1.18d+173)) then
tmp = t_0
else if (z <= (-8100000000.0d0)) then
tmp = t_1
else if (z <= 12.0d0) then
tmp = x * (-3.0d0)
else if (z <= 2.3d+222) then
tmp = t_0
else
tmp = t_1
end if
code = tmp
end function
public static double code(double x, double y, double z) {
double t_0 = x * (z * 6.0);
double t_1 = -6.0 * (y * z);
double tmp;
if (z <= -1.18e+173) {
tmp = t_0;
} else if (z <= -8100000000.0) {
tmp = t_1;
} else if (z <= 12.0) {
tmp = x * -3.0;
} else if (z <= 2.3e+222) {
tmp = t_0;
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z): t_0 = x * (z * 6.0) t_1 = -6.0 * (y * z) tmp = 0 if z <= -1.18e+173: tmp = t_0 elif z <= -8100000000.0: tmp = t_1 elif z <= 12.0: tmp = x * -3.0 elif z <= 2.3e+222: tmp = t_0 else: tmp = t_1 return tmp
function code(x, y, z) t_0 = Float64(x * Float64(z * 6.0)) t_1 = Float64(-6.0 * Float64(y * z)) tmp = 0.0 if (z <= -1.18e+173) tmp = t_0; elseif (z <= -8100000000.0) tmp = t_1; elseif (z <= 12.0) tmp = Float64(x * -3.0); elseif (z <= 2.3e+222) tmp = t_0; else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z) t_0 = x * (z * 6.0); t_1 = -6.0 * (y * z); tmp = 0.0; if (z <= -1.18e+173) tmp = t_0; elseif (z <= -8100000000.0) tmp = t_1; elseif (z <= 12.0) tmp = x * -3.0; elseif (z <= 2.3e+222) tmp = t_0; else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_] := Block[{t$95$0 = N[(x * N[(z * 6.0), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(-6.0 * N[(y * z), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[z, -1.18e+173], t$95$0, If[LessEqual[z, -8100000000.0], t$95$1, If[LessEqual[z, 12.0], N[(x * -3.0), $MachinePrecision], If[LessEqual[z, 2.3e+222], t$95$0, t$95$1]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := x \cdot \left(z \cdot 6\right)\\
t_1 := -6 \cdot \left(y \cdot z\right)\\
\mathbf{if}\;z \leq -1.18 \cdot 10^{+173}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;z \leq -8100000000:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;z \leq 12:\\
\;\;\;\;x \cdot -3\\
\mathbf{elif}\;z \leq 2.3 \cdot 10^{+222}:\\
\;\;\;\;t\_0\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if z < -1.18e173 or 12 < z < 2.30000000000000011e222Initial program 99.7%
metadata-eval99.7%
Simplified99.7%
Taylor expanded in x around inf 61.1%
sub-neg61.1%
distribute-rgt-in61.1%
metadata-eval61.1%
distribute-lft-neg-in61.1%
associate-+r+61.1%
metadata-eval61.1%
distribute-rgt-neg-in61.1%
metadata-eval61.1%
Simplified61.1%
Taylor expanded in z around inf 59.8%
if -1.18e173 < z < -8.1e9 or 2.30000000000000011e222 < z Initial program 99.7%
metadata-eval99.7%
Simplified99.7%
Taylor expanded in z around 0 99.8%
Taylor expanded in z around inf 99.4%
associate-*r*99.4%
metadata-eval99.4%
distribute-lft-neg-in99.4%
distribute-lft-neg-in99.4%
*-commutative99.4%
associate-*r*99.4%
rem-cube-cbrt98.5%
rem-cube-cbrt99.4%
*-commutative99.4%
distribute-rgt-neg-in99.4%
*-commutative99.4%
rem-cube-cbrt98.5%
mul-1-neg98.5%
rem-cube-cbrt99.4%
associate-*r*99.4%
metadata-eval99.4%
Simplified99.4%
Taylor expanded in y around inf 64.9%
if -8.1e9 < z < 12Initial program 99.4%
metadata-eval99.4%
Simplified99.4%
Taylor expanded in x around inf 43.4%
sub-neg43.4%
distribute-rgt-in43.4%
metadata-eval43.4%
distribute-lft-neg-in43.4%
associate-+r+43.4%
metadata-eval43.4%
distribute-rgt-neg-in43.4%
metadata-eval43.4%
Simplified43.4%
Taylor expanded in z around 0 41.3%
*-commutative41.3%
Simplified41.3%
Final simplification51.8%
(FPCore (x y z)
:precision binary64
(let* ((t_0 (* 6.0 (* x z))) (t_1 (* -6.0 (* y z))))
(if (<= z -1.08e+172)
t_0
(if (<= z -8100000000.0)
t_1
(if (<= z 12.0) (* x -3.0) (if (<= z 1.6e+224) t_0 t_1))))))
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 <= -1.08e+172) {
tmp = t_0;
} else if (z <= -8100000000.0) {
tmp = t_1;
} else if (z <= 12.0) {
tmp = x * -3.0;
} else if (z <= 1.6e+224) {
tmp = t_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 * (x * z)
t_1 = (-6.0d0) * (y * z)
if (z <= (-1.08d+172)) then
tmp = t_0
else if (z <= (-8100000000.0d0)) then
tmp = t_1
else if (z <= 12.0d0) then
tmp = x * (-3.0d0)
else if (z <= 1.6d+224) then
tmp = t_0
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 * (x * z);
double t_1 = -6.0 * (y * z);
double tmp;
if (z <= -1.08e+172) {
tmp = t_0;
} else if (z <= -8100000000.0) {
tmp = t_1;
} else if (z <= 12.0) {
tmp = x * -3.0;
} else if (z <= 1.6e+224) {
tmp = t_0;
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z): t_0 = 6.0 * (x * z) t_1 = -6.0 * (y * z) tmp = 0 if z <= -1.08e+172: tmp = t_0 elif z <= -8100000000.0: tmp = t_1 elif z <= 12.0: tmp = x * -3.0 elif z <= 1.6e+224: tmp = t_0 else: tmp = t_1 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 <= -1.08e+172) tmp = t_0; elseif (z <= -8100000000.0) tmp = t_1; elseif (z <= 12.0) tmp = Float64(x * -3.0); elseif (z <= 1.6e+224) tmp = t_0; else tmp = t_1; 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 <= -1.08e+172) tmp = t_0; elseif (z <= -8100000000.0) tmp = t_1; elseif (z <= 12.0) tmp = x * -3.0; elseif (z <= 1.6e+224) tmp = t_0; else tmp = t_1; 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, -1.08e+172], t$95$0, If[LessEqual[z, -8100000000.0], t$95$1, If[LessEqual[z, 12.0], N[(x * -3.0), $MachinePrecision], If[LessEqual[z, 1.6e+224], t$95$0, t$95$1]]]]]]
\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.08 \cdot 10^{+172}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;z \leq -8100000000:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;z \leq 12:\\
\;\;\;\;x \cdot -3\\
\mathbf{elif}\;z \leq 1.6 \cdot 10^{+224}:\\
\;\;\;\;t\_0\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if z < -1.0799999999999999e172 or 12 < z < 1.60000000000000007e224Initial program 99.7%
metadata-eval99.7%
Simplified99.7%
Taylor expanded in z around 0 99.6%
Taylor expanded in z around inf 98.3%
associate-*r*98.6%
metadata-eval98.6%
distribute-lft-neg-in98.6%
distribute-lft-neg-in98.6%
*-commutative98.6%
associate-*r*98.5%
rem-cube-cbrt97.3%
rem-cube-cbrt98.5%
*-commutative98.5%
distribute-rgt-neg-in98.5%
*-commutative98.5%
rem-cube-cbrt97.3%
mul-1-neg97.3%
rem-cube-cbrt98.5%
associate-*r*98.5%
metadata-eval98.5%
Simplified98.5%
Taylor expanded in y around 0 59.7%
if -1.0799999999999999e172 < z < -8.1e9 or 1.60000000000000007e224 < z Initial program 99.7%
metadata-eval99.7%
Simplified99.7%
Taylor expanded in z around 0 99.8%
Taylor expanded in z around inf 99.4%
associate-*r*99.4%
metadata-eval99.4%
distribute-lft-neg-in99.4%
distribute-lft-neg-in99.4%
*-commutative99.4%
associate-*r*99.4%
rem-cube-cbrt98.5%
rem-cube-cbrt99.4%
*-commutative99.4%
distribute-rgt-neg-in99.4%
*-commutative99.4%
rem-cube-cbrt98.5%
mul-1-neg98.5%
rem-cube-cbrt99.4%
associate-*r*99.4%
metadata-eval99.4%
Simplified99.4%
Taylor expanded in y around inf 64.9%
if -8.1e9 < z < 12Initial program 99.4%
metadata-eval99.4%
Simplified99.4%
Taylor expanded in x around inf 43.4%
sub-neg43.4%
distribute-rgt-in43.4%
metadata-eval43.4%
distribute-lft-neg-in43.4%
associate-+r+43.4%
metadata-eval43.4%
distribute-rgt-neg-in43.4%
metadata-eval43.4%
Simplified43.4%
Taylor expanded in z around 0 41.3%
*-commutative41.3%
Simplified41.3%
(FPCore (x y z) :precision binary64 (if (or (<= z -0.58) (not (<= z 0.5))) (* (- y x) (* z -6.0)) (+ x (* (- y x) 4.0))))
double code(double x, double y, double z) {
double tmp;
if ((z <= -0.58) || !(z <= 0.5)) {
tmp = (y - x) * (z * -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.5d0))) then
tmp = (y - x) * (z * (-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.5)) {
tmp = (y - x) * (z * -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.5): tmp = (y - x) * (z * -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.5)) tmp = Float64(Float64(y - x) * Float64(z * -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.5))) tmp = (y - x) * (z * -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.5]], $MachinePrecision]], N[(N[(y - x), $MachinePrecision] * N[(z * -6.0), $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.5\right):\\
\;\;\;\;\left(y - x\right) \cdot \left(z \cdot -6\right)\\
\mathbf{else}:\\
\;\;\;\;x + \left(y - x\right) \cdot 4\\
\end{array}
\end{array}
if z < -0.57999999999999996 or 0.5 < z Initial program 99.7%
+-commutative99.7%
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%
flip-+75.5%
metadata-eval75.5%
swap-sqr75.4%
pow275.4%
metadata-eval75.4%
Applied egg-rr75.4%
Taylor expanded in z around inf 97.1%
associate-*r*97.2%
Simplified97.2%
if -0.57999999999999996 < z < 0.5Initial program 99.4%
metadata-eval99.4%
Simplified99.4%
Taylor expanded in z around 0 99.6%
Final simplification98.3%
(FPCore (x y z) :precision binary64 (if (or (<= x -1.12e-40) (not (<= x 6.5e+31))) (* x (+ -3.0 (* z 6.0))) (* y (+ 4.0 (* z -6.0)))))
double code(double x, double y, double z) {
double tmp;
if ((x <= -1.12e-40) || !(x <= 6.5e+31)) {
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 <= (-1.12d-40)) .or. (.not. (x <= 6.5d+31))) 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 <= -1.12e-40) || !(x <= 6.5e+31)) {
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 <= -1.12e-40) or not (x <= 6.5e+31): 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 <= -1.12e-40) || !(x <= 6.5e+31)) 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 <= -1.12e-40) || ~((x <= 6.5e+31))) 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, -1.12e-40], N[Not[LessEqual[x, 6.5e+31]], $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 -1.12 \cdot 10^{-40} \lor \neg \left(x \leq 6.5 \cdot 10^{+31}\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.1200000000000001e-40 or 6.5000000000000004e31 < x Initial program 99.6%
metadata-eval99.6%
Simplified99.6%
Taylor expanded in x around inf 79.8%
sub-neg79.8%
distribute-rgt-in79.8%
metadata-eval79.8%
distribute-lft-neg-in79.8%
associate-+r+79.8%
metadata-eval79.8%
distribute-rgt-neg-in79.8%
metadata-eval79.8%
Simplified79.8%
if -1.1200000000000001e-40 < x < 6.5000000000000004e31Initial program 99.5%
+-commutative99.5%
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 84.5%
Final simplification82.1%
(FPCore (x y z) :precision binary64 (if (or (<= z -8100000000.0) (not (<= z 0.65))) (* -6.0 (* y z)) (* x -3.0)))
double code(double x, double y, double z) {
double tmp;
if ((z <= -8100000000.0) || !(z <= 0.65)) {
tmp = -6.0 * (y * z);
} 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) :: tmp
if ((z <= (-8100000000.0d0)) .or. (.not. (z <= 0.65d0))) then
tmp = (-6.0d0) * (y * z)
else
tmp = x * (-3.0d0)
end if
code = tmp
end function
public static double code(double x, double y, double z) {
double tmp;
if ((z <= -8100000000.0) || !(z <= 0.65)) {
tmp = -6.0 * (y * z);
} else {
tmp = x * -3.0;
}
return tmp;
}
def code(x, y, z): tmp = 0 if (z <= -8100000000.0) or not (z <= 0.65): tmp = -6.0 * (y * z) else: tmp = x * -3.0 return tmp
function code(x, y, z) tmp = 0.0 if ((z <= -8100000000.0) || !(z <= 0.65)) tmp = Float64(-6.0 * Float64(y * z)); else tmp = Float64(x * -3.0); end return tmp end
function tmp_2 = code(x, y, z) tmp = 0.0; if ((z <= -8100000000.0) || ~((z <= 0.65))) tmp = -6.0 * (y * z); else tmp = x * -3.0; end tmp_2 = tmp; end
code[x_, y_, z_] := If[Or[LessEqual[z, -8100000000.0], N[Not[LessEqual[z, 0.65]], $MachinePrecision]], N[(-6.0 * N[(y * z), $MachinePrecision]), $MachinePrecision], N[(x * -3.0), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;z \leq -8100000000 \lor \neg \left(z \leq 0.65\right):\\
\;\;\;\;-6 \cdot \left(y \cdot z\right)\\
\mathbf{else}:\\
\;\;\;\;x \cdot -3\\
\end{array}
\end{array}
if z < -8.1e9 or 0.650000000000000022 < z Initial program 99.7%
metadata-eval99.7%
Simplified99.7%
Taylor expanded in z around 0 99.7%
Taylor expanded in z around inf 98.2%
associate-*r*98.3%
metadata-eval98.3%
distribute-lft-neg-in98.3%
distribute-lft-neg-in98.3%
*-commutative98.3%
associate-*r*98.2%
rem-cube-cbrt97.2%
rem-cube-cbrt98.2%
*-commutative98.2%
distribute-rgt-neg-in98.2%
*-commutative98.2%
rem-cube-cbrt97.2%
mul-1-neg97.2%
rem-cube-cbrt98.2%
associate-*r*98.2%
metadata-eval98.2%
Simplified98.2%
Taylor expanded in y around inf 50.5%
if -8.1e9 < z < 0.650000000000000022Initial program 99.4%
metadata-eval99.4%
Simplified99.4%
Taylor expanded in x around inf 43.7%
sub-neg43.7%
distribute-rgt-in43.7%
metadata-eval43.7%
distribute-lft-neg-in43.7%
associate-+r+43.7%
metadata-eval43.7%
distribute-rgt-neg-in43.7%
metadata-eval43.7%
Simplified43.7%
Taylor expanded in z around 0 41.6%
*-commutative41.6%
Simplified41.6%
Final simplification46.2%
(FPCore (x y z) :precision binary64 (+ x (+ (* -6.0 (* (- y x) z)) (* (- y x) 4.0))))
double code(double x, double y, double z) {
return x + ((-6.0 * ((y - x) * z)) + ((y - x) * 4.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 + (((-6.0d0) * ((y - x) * z)) + ((y - x) * 4.0d0))
end function
public static double code(double x, double y, double z) {
return x + ((-6.0 * ((y - x) * z)) + ((y - x) * 4.0));
}
def code(x, y, z): return x + ((-6.0 * ((y - x) * z)) + ((y - x) * 4.0))
function code(x, y, z) return Float64(x + Float64(Float64(-6.0 * Float64(Float64(y - x) * z)) + Float64(Float64(y - x) * 4.0))) end
function tmp = code(x, y, z) tmp = x + ((-6.0 * ((y - x) * z)) + ((y - x) * 4.0)); end
code[x_, y_, z_] := N[(x + N[(N[(-6.0 * N[(N[(y - x), $MachinePrecision] * z), $MachinePrecision]), $MachinePrecision] + N[(N[(y - x), $MachinePrecision] * 4.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
x + \left(-6 \cdot \left(\left(y - x\right) \cdot z\right) + \left(y - x\right) \cdot 4\right)
\end{array}
Initial program 99.5%
metadata-eval99.5%
Simplified99.5%
Taylor expanded in z around 0 99.8%
Final simplification99.8%
(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.5%
metadata-eval99.5%
Simplified99.5%
(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.5%
metadata-eval99.5%
Simplified99.5%
Taylor expanded in x around inf 48.6%
sub-neg48.6%
distribute-rgt-in48.6%
metadata-eval48.6%
distribute-lft-neg-in48.6%
associate-+r+48.6%
metadata-eval48.6%
distribute-rgt-neg-in48.6%
metadata-eval48.6%
Simplified48.6%
Taylor expanded in z around 0 21.5%
*-commutative21.5%
Simplified21.5%
herbie shell --seed 2024086
(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))))