
(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.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%
Final simplification99.8%
(FPCore (x y z)
:precision binary64
(let* ((t_0 (* -6.0 (* y z))) (t_1 (+ x (* y 4.0))))
(if (<= z -0.3)
t_0
(if (<= z -9e-203)
t_1
(if (<= z -2.5e-299)
(* x -3.0)
(if (<= z 7e-300)
t_1
(if (<= z 3.9e-200)
(* x -3.0)
(if (<= z 1.3e-22)
t_1
(if (<= z 0.5)
(* x -3.0)
(if (<= z 3.3e+164) (* x (* z 6.0)) t_0))))))))))
double code(double x, double y, double z) {
double t_0 = -6.0 * (y * z);
double t_1 = x + (y * 4.0);
double tmp;
if (z <= -0.3) {
tmp = t_0;
} else if (z <= -9e-203) {
tmp = t_1;
} else if (z <= -2.5e-299) {
tmp = x * -3.0;
} else if (z <= 7e-300) {
tmp = t_1;
} else if (z <= 3.9e-200) {
tmp = x * -3.0;
} else if (z <= 1.3e-22) {
tmp = t_1;
} else if (z <= 0.5) {
tmp = x * -3.0;
} else if (z <= 3.3e+164) {
tmp = x * (z * 6.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) :: t_1
real(8) :: tmp
t_0 = (-6.0d0) * (y * z)
t_1 = x + (y * 4.0d0)
if (z <= (-0.3d0)) then
tmp = t_0
else if (z <= (-9d-203)) then
tmp = t_1
else if (z <= (-2.5d-299)) then
tmp = x * (-3.0d0)
else if (z <= 7d-300) then
tmp = t_1
else if (z <= 3.9d-200) then
tmp = x * (-3.0d0)
else if (z <= 1.3d-22) then
tmp = t_1
else if (z <= 0.5d0) then
tmp = x * (-3.0d0)
else if (z <= 3.3d+164) then
tmp = x * (z * 6.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 t_1 = x + (y * 4.0);
double tmp;
if (z <= -0.3) {
tmp = t_0;
} else if (z <= -9e-203) {
tmp = t_1;
} else if (z <= -2.5e-299) {
tmp = x * -3.0;
} else if (z <= 7e-300) {
tmp = t_1;
} else if (z <= 3.9e-200) {
tmp = x * -3.0;
} else if (z <= 1.3e-22) {
tmp = t_1;
} else if (z <= 0.5) {
tmp = x * -3.0;
} else if (z <= 3.3e+164) {
tmp = x * (z * 6.0);
} else {
tmp = t_0;
}
return tmp;
}
def code(x, y, z): t_0 = -6.0 * (y * z) t_1 = x + (y * 4.0) tmp = 0 if z <= -0.3: tmp = t_0 elif z <= -9e-203: tmp = t_1 elif z <= -2.5e-299: tmp = x * -3.0 elif z <= 7e-300: tmp = t_1 elif z <= 3.9e-200: tmp = x * -3.0 elif z <= 1.3e-22: tmp = t_1 elif z <= 0.5: tmp = x * -3.0 elif z <= 3.3e+164: tmp = x * (z * 6.0) else: tmp = t_0 return tmp
function code(x, y, z) t_0 = Float64(-6.0 * Float64(y * z)) t_1 = Float64(x + Float64(y * 4.0)) tmp = 0.0 if (z <= -0.3) tmp = t_0; elseif (z <= -9e-203) tmp = t_1; elseif (z <= -2.5e-299) tmp = Float64(x * -3.0); elseif (z <= 7e-300) tmp = t_1; elseif (z <= 3.9e-200) tmp = Float64(x * -3.0); elseif (z <= 1.3e-22) tmp = t_1; elseif (z <= 0.5) tmp = Float64(x * -3.0); elseif (z <= 3.3e+164) tmp = Float64(x * Float64(z * 6.0)); else tmp = t_0; end return tmp end
function tmp_2 = code(x, y, z) t_0 = -6.0 * (y * z); t_1 = x + (y * 4.0); tmp = 0.0; if (z <= -0.3) tmp = t_0; elseif (z <= -9e-203) tmp = t_1; elseif (z <= -2.5e-299) tmp = x * -3.0; elseif (z <= 7e-300) tmp = t_1; elseif (z <= 3.9e-200) tmp = x * -3.0; elseif (z <= 1.3e-22) tmp = t_1; elseif (z <= 0.5) tmp = x * -3.0; elseif (z <= 3.3e+164) tmp = x * (z * 6.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]}, Block[{t$95$1 = N[(x + N[(y * 4.0), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[z, -0.3], t$95$0, If[LessEqual[z, -9e-203], t$95$1, If[LessEqual[z, -2.5e-299], N[(x * -3.0), $MachinePrecision], If[LessEqual[z, 7e-300], t$95$1, If[LessEqual[z, 3.9e-200], N[(x * -3.0), $MachinePrecision], If[LessEqual[z, 1.3e-22], t$95$1, If[LessEqual[z, 0.5], N[(x * -3.0), $MachinePrecision], If[LessEqual[z, 3.3e+164], N[(x * N[(z * 6.0), $MachinePrecision]), $MachinePrecision], t$95$0]]]]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := -6 \cdot \left(y \cdot z\right)\\
t_1 := x + y \cdot 4\\
\mathbf{if}\;z \leq -0.3:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;z \leq -9 \cdot 10^{-203}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;z \leq -2.5 \cdot 10^{-299}:\\
\;\;\;\;x \cdot -3\\
\mathbf{elif}\;z \leq 7 \cdot 10^{-300}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;z \leq 3.9 \cdot 10^{-200}:\\
\;\;\;\;x \cdot -3\\
\mathbf{elif}\;z \leq 1.3 \cdot 10^{-22}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;z \leq 0.5:\\
\;\;\;\;x \cdot -3\\
\mathbf{elif}\;z \leq 3.3 \cdot 10^{+164}:\\
\;\;\;\;x \cdot \left(z \cdot 6\right)\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
if z < -0.299999999999999989 or 3.29999999999999995e164 < z Initial program 99.8%
metadata-eval99.8%
Simplified99.8%
Taylor expanded in z around 0 99.8%
Taylor expanded in z around inf 98.0%
Taylor expanded in y around inf 60.3%
*-commutative60.3%
Simplified60.3%
if -0.299999999999999989 < z < -9.0000000000000003e-203 or -2.49999999999999978e-299 < z < 7.0000000000000003e-300 or 3.89999999999999999e-200 < z < 1.3e-22Initial program 99.4%
metadata-eval99.4%
Simplified99.4%
Taylor expanded in y around inf 65.8%
associate-*r*65.8%
Simplified65.8%
Taylor expanded in z around 0 66.0%
*-commutative66.0%
Simplified66.0%
if -9.0000000000000003e-203 < z < -2.49999999999999978e-299 or 7.0000000000000003e-300 < z < 3.89999999999999999e-200 or 1.3e-22 < z < 0.5Initial program 99.3%
metadata-eval99.3%
Simplified99.3%
Taylor expanded in x around inf 69.2%
sub-neg69.2%
distribute-lft-in69.2%
metadata-eval69.2%
metadata-eval69.2%
neg-mul-169.2%
*-commutative69.2%
associate-*l*69.2%
distribute-rgt-in69.2%
distribute-lft-in69.2%
metadata-eval69.2%
metadata-eval69.2%
associate-+r+69.2%
metadata-eval69.2%
metadata-eval69.2%
associate-*r*69.2%
metadata-eval69.2%
*-commutative69.2%
Simplified69.2%
Taylor expanded in z around 0 67.1%
*-commutative67.1%
Simplified67.1%
if 0.5 < z < 3.29999999999999995e164Initial program 99.5%
metadata-eval99.5%
Simplified99.5%
Taylor expanded in x around inf 67.1%
sub-neg67.1%
distribute-lft-in67.1%
metadata-eval67.1%
metadata-eval67.1%
neg-mul-167.1%
*-commutative67.1%
associate-*l*67.1%
distribute-rgt-in67.1%
distribute-lft-in67.1%
metadata-eval67.1%
metadata-eval67.1%
associate-+r+67.1%
metadata-eval67.1%
metadata-eval67.1%
associate-*r*67.1%
metadata-eval67.1%
*-commutative67.1%
Simplified67.1%
Taylor expanded in z around inf 62.8%
Final simplification63.6%
(FPCore (x y z)
:precision binary64
(let* ((t_0 (+ x (* y 4.0))) (t_1 (* -6.0 (* (- y x) z))))
(if (<= z -0.0295)
t_1
(if (<= z -2.25e-203)
t_0
(if (<= z -5.5e-299)
(* x -3.0)
(if (<= z 1.05e-294)
t_0
(if (<= z 3.5e-200)
(* x -3.0)
(if (<= z 1.05e-22) t_0 (if (<= z 0.5) (* x -3.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 <= -0.0295) {
tmp = t_1;
} else if (z <= -2.25e-203) {
tmp = t_0;
} else if (z <= -5.5e-299) {
tmp = x * -3.0;
} else if (z <= 1.05e-294) {
tmp = t_0;
} else if (z <= 3.5e-200) {
tmp = x * -3.0;
} else if (z <= 1.05e-22) {
tmp = t_0;
} else if (z <= 0.5) {
tmp = x * -3.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 <= (-0.0295d0)) then
tmp = t_1
else if (z <= (-2.25d-203)) then
tmp = t_0
else if (z <= (-5.5d-299)) then
tmp = x * (-3.0d0)
else if (z <= 1.05d-294) then
tmp = t_0
else if (z <= 3.5d-200) then
tmp = x * (-3.0d0)
else if (z <= 1.05d-22) then
tmp = t_0
else if (z <= 0.5d0) then
tmp = x * (-3.0d0)
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 <= -0.0295) {
tmp = t_1;
} else if (z <= -2.25e-203) {
tmp = t_0;
} else if (z <= -5.5e-299) {
tmp = x * -3.0;
} else if (z <= 1.05e-294) {
tmp = t_0;
} else if (z <= 3.5e-200) {
tmp = x * -3.0;
} else if (z <= 1.05e-22) {
tmp = t_0;
} else if (z <= 0.5) {
tmp = x * -3.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 <= -0.0295: tmp = t_1 elif z <= -2.25e-203: tmp = t_0 elif z <= -5.5e-299: tmp = x * -3.0 elif z <= 1.05e-294: tmp = t_0 elif z <= 3.5e-200: tmp = x * -3.0 elif z <= 1.05e-22: tmp = t_0 elif z <= 0.5: tmp = x * -3.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 <= -0.0295) tmp = t_1; elseif (z <= -2.25e-203) tmp = t_0; elseif (z <= -5.5e-299) tmp = Float64(x * -3.0); elseif (z <= 1.05e-294) tmp = t_0; elseif (z <= 3.5e-200) tmp = Float64(x * -3.0); elseif (z <= 1.05e-22) tmp = t_0; elseif (z <= 0.5) tmp = Float64(x * -3.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 <= -0.0295) tmp = t_1; elseif (z <= -2.25e-203) tmp = t_0; elseif (z <= -5.5e-299) tmp = x * -3.0; elseif (z <= 1.05e-294) tmp = t_0; elseif (z <= 3.5e-200) tmp = x * -3.0; elseif (z <= 1.05e-22) tmp = t_0; elseif (z <= 0.5) tmp = x * -3.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, -0.0295], t$95$1, If[LessEqual[z, -2.25e-203], t$95$0, If[LessEqual[z, -5.5e-299], N[(x * -3.0), $MachinePrecision], If[LessEqual[z, 1.05e-294], t$95$0, If[LessEqual[z, 3.5e-200], N[(x * -3.0), $MachinePrecision], If[LessEqual[z, 1.05e-22], t$95$0, If[LessEqual[z, 0.5], N[(x * -3.0), $MachinePrecision], 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 -0.0295:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;z \leq -2.25 \cdot 10^{-203}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;z \leq -5.5 \cdot 10^{-299}:\\
\;\;\;\;x \cdot -3\\
\mathbf{elif}\;z \leq 1.05 \cdot 10^{-294}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;z \leq 3.5 \cdot 10^{-200}:\\
\;\;\;\;x \cdot -3\\
\mathbf{elif}\;z \leq 1.05 \cdot 10^{-22}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;z \leq 0.5:\\
\;\;\;\;x \cdot -3\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if z < -0.029499999999999998 or 0.5 < z Initial program 99.7%
metadata-eval99.7%
Simplified99.7%
Taylor expanded in z around 0 99.7%
Taylor expanded in z around inf 97.2%
if -0.029499999999999998 < z < -2.2500000000000001e-203 or -5.5e-299 < z < 1.04999999999999992e-294 or 3.50000000000000023e-200 < z < 1.05000000000000004e-22Initial program 99.4%
metadata-eval99.4%
Simplified99.4%
Taylor expanded in y around inf 65.8%
associate-*r*65.8%
Simplified65.8%
Taylor expanded in z around 0 66.0%
*-commutative66.0%
Simplified66.0%
if -2.2500000000000001e-203 < z < -5.5e-299 or 1.04999999999999992e-294 < z < 3.50000000000000023e-200 or 1.05000000000000004e-22 < z < 0.5Initial program 99.3%
metadata-eval99.3%
Simplified99.3%
Taylor expanded in x around inf 69.2%
sub-neg69.2%
distribute-lft-in69.2%
metadata-eval69.2%
metadata-eval69.2%
neg-mul-169.2%
*-commutative69.2%
associate-*l*69.2%
distribute-rgt-in69.2%
distribute-lft-in69.2%
metadata-eval69.2%
metadata-eval69.2%
associate-+r+69.2%
metadata-eval69.2%
metadata-eval69.2%
associate-*r*69.2%
metadata-eval69.2%
*-commutative69.2%
Simplified69.2%
Taylor expanded in z around 0 67.1%
*-commutative67.1%
Simplified67.1%
Final simplification82.0%
(FPCore (x y z)
:precision binary64
(let* ((t_0 (+ x (* y 4.0))) (t_1 (* -6.0 (* (- y x) z))))
(if (<= z -0.008)
t_1
(if (<= z -8e-203)
t_0
(if (<= z -1.2e-299)
(* x -3.0)
(if (<= z 2e-294)
t_0
(if (<= z 4.9e-198)
(* x -3.0)
(if (<= z 1.3e-22)
t_0
(if (<= z 1100000.0) (* x (+ -3.0 (* z 6.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 <= -0.008) {
tmp = t_1;
} else if (z <= -8e-203) {
tmp = t_0;
} else if (z <= -1.2e-299) {
tmp = x * -3.0;
} else if (z <= 2e-294) {
tmp = t_0;
} else if (z <= 4.9e-198) {
tmp = x * -3.0;
} else if (z <= 1.3e-22) {
tmp = t_0;
} else if (z <= 1100000.0) {
tmp = x * (-3.0 + (z * 6.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 <= (-0.008d0)) then
tmp = t_1
else if (z <= (-8d-203)) then
tmp = t_0
else if (z <= (-1.2d-299)) then
tmp = x * (-3.0d0)
else if (z <= 2d-294) then
tmp = t_0
else if (z <= 4.9d-198) then
tmp = x * (-3.0d0)
else if (z <= 1.3d-22) then
tmp = t_0
else if (z <= 1100000.0d0) then
tmp = x * ((-3.0d0) + (z * 6.0d0))
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 <= -0.008) {
tmp = t_1;
} else if (z <= -8e-203) {
tmp = t_0;
} else if (z <= -1.2e-299) {
tmp = x * -3.0;
} else if (z <= 2e-294) {
tmp = t_0;
} else if (z <= 4.9e-198) {
tmp = x * -3.0;
} else if (z <= 1.3e-22) {
tmp = t_0;
} else if (z <= 1100000.0) {
tmp = x * (-3.0 + (z * 6.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 <= -0.008: tmp = t_1 elif z <= -8e-203: tmp = t_0 elif z <= -1.2e-299: tmp = x * -3.0 elif z <= 2e-294: tmp = t_0 elif z <= 4.9e-198: tmp = x * -3.0 elif z <= 1.3e-22: tmp = t_0 elif z <= 1100000.0: tmp = x * (-3.0 + (z * 6.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 <= -0.008) tmp = t_1; elseif (z <= -8e-203) tmp = t_0; elseif (z <= -1.2e-299) tmp = Float64(x * -3.0); elseif (z <= 2e-294) tmp = t_0; elseif (z <= 4.9e-198) tmp = Float64(x * -3.0); elseif (z <= 1.3e-22) tmp = t_0; elseif (z <= 1100000.0) tmp = Float64(x * Float64(-3.0 + Float64(z * 6.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 <= -0.008) tmp = t_1; elseif (z <= -8e-203) tmp = t_0; elseif (z <= -1.2e-299) tmp = x * -3.0; elseif (z <= 2e-294) tmp = t_0; elseif (z <= 4.9e-198) tmp = x * -3.0; elseif (z <= 1.3e-22) tmp = t_0; elseif (z <= 1100000.0) tmp = x * (-3.0 + (z * 6.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, -0.008], t$95$1, If[LessEqual[z, -8e-203], t$95$0, If[LessEqual[z, -1.2e-299], N[(x * -3.0), $MachinePrecision], If[LessEqual[z, 2e-294], t$95$0, If[LessEqual[z, 4.9e-198], N[(x * -3.0), $MachinePrecision], If[LessEqual[z, 1.3e-22], t$95$0, If[LessEqual[z, 1100000.0], N[(x * N[(-3.0 + N[(z * 6.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 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 -0.008:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;z \leq -8 \cdot 10^{-203}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;z \leq -1.2 \cdot 10^{-299}:\\
\;\;\;\;x \cdot -3\\
\mathbf{elif}\;z \leq 2 \cdot 10^{-294}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;z \leq 4.9 \cdot 10^{-198}:\\
\;\;\;\;x \cdot -3\\
\mathbf{elif}\;z \leq 1.3 \cdot 10^{-22}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;z \leq 1100000:\\
\;\;\;\;x \cdot \left(-3 + z \cdot 6\right)\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if z < -0.0080000000000000002 or 1.1e6 < 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%
if -0.0080000000000000002 < z < -8.0000000000000003e-203 or -1.2000000000000001e-299 < z < 2.00000000000000003e-294 or 4.9000000000000002e-198 < z < 1.3e-22Initial program 99.4%
metadata-eval99.4%
Simplified99.4%
Taylor expanded in y around inf 65.8%
associate-*r*65.8%
Simplified65.8%
Taylor expanded in z around 0 66.0%
*-commutative66.0%
Simplified66.0%
if -8.0000000000000003e-203 < z < -1.2000000000000001e-299 or 2.00000000000000003e-294 < z < 4.9000000000000002e-198Initial program 99.3%
metadata-eval99.3%
Simplified99.3%
Taylor expanded in x around inf 68.9%
sub-neg68.9%
distribute-lft-in68.9%
metadata-eval68.9%
metadata-eval68.9%
neg-mul-168.9%
*-commutative68.9%
associate-*l*68.9%
distribute-rgt-in68.9%
distribute-lft-in68.9%
metadata-eval68.9%
metadata-eval68.9%
associate-+r+68.9%
metadata-eval68.9%
metadata-eval68.9%
associate-*r*68.9%
metadata-eval68.9%
*-commutative68.9%
Simplified68.9%
Taylor expanded in z around 0 68.9%
*-commutative68.9%
Simplified68.9%
if 1.3e-22 < z < 1.1e6Initial program 99.6%
metadata-eval99.6%
Simplified99.6%
Taylor expanded in x around inf 79.6%
sub-neg79.6%
distribute-lft-in79.6%
metadata-eval79.6%
metadata-eval79.6%
neg-mul-179.6%
*-commutative79.6%
associate-*l*79.6%
distribute-rgt-in79.6%
distribute-lft-in79.6%
metadata-eval79.6%
metadata-eval79.6%
associate-+r+79.6%
metadata-eval79.6%
metadata-eval79.6%
associate-*r*79.6%
metadata-eval79.6%
*-commutative79.6%
Simplified79.6%
Final simplification82.9%
(FPCore (x y z)
:precision binary64
(let* ((t_0 (* y (+ 4.0 (* z -6.0)))))
(if (<= y -1.45e+69)
t_0
(if (<= y -3.5e-49)
(* -6.0 (* (- y x) z))
(if (or (<= y -8.5e-53) (not (<= y 2.25e+34)))
t_0
(* x (+ -3.0 (* z 6.0))))))))
double code(double x, double y, double z) {
double t_0 = y * (4.0 + (z * -6.0));
double tmp;
if (y <= -1.45e+69) {
tmp = t_0;
} else if (y <= -3.5e-49) {
tmp = -6.0 * ((y - x) * z);
} else if ((y <= -8.5e-53) || !(y <= 2.25e+34)) {
tmp = t_0;
} else {
tmp = x * (-3.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) :: t_0
real(8) :: tmp
t_0 = y * (4.0d0 + (z * (-6.0d0)))
if (y <= (-1.45d+69)) then
tmp = t_0
else if (y <= (-3.5d-49)) then
tmp = (-6.0d0) * ((y - x) * z)
else if ((y <= (-8.5d-53)) .or. (.not. (y <= 2.25d+34))) then
tmp = t_0
else
tmp = x * ((-3.0d0) + (z * 6.0d0))
end if
code = tmp
end function
public static double code(double x, double y, double z) {
double t_0 = y * (4.0 + (z * -6.0));
double tmp;
if (y <= -1.45e+69) {
tmp = t_0;
} else if (y <= -3.5e-49) {
tmp = -6.0 * ((y - x) * z);
} else if ((y <= -8.5e-53) || !(y <= 2.25e+34)) {
tmp = t_0;
} else {
tmp = x * (-3.0 + (z * 6.0));
}
return tmp;
}
def code(x, y, z): t_0 = y * (4.0 + (z * -6.0)) tmp = 0 if y <= -1.45e+69: tmp = t_0 elif y <= -3.5e-49: tmp = -6.0 * ((y - x) * z) elif (y <= -8.5e-53) or not (y <= 2.25e+34): tmp = t_0 else: tmp = x * (-3.0 + (z * 6.0)) return tmp
function code(x, y, z) t_0 = Float64(y * Float64(4.0 + Float64(z * -6.0))) tmp = 0.0 if (y <= -1.45e+69) tmp = t_0; elseif (y <= -3.5e-49) tmp = Float64(-6.0 * Float64(Float64(y - x) * z)); elseif ((y <= -8.5e-53) || !(y <= 2.25e+34)) tmp = t_0; else tmp = Float64(x * Float64(-3.0 + Float64(z * 6.0))); end return tmp end
function tmp_2 = code(x, y, z) t_0 = y * (4.0 + (z * -6.0)); tmp = 0.0; if (y <= -1.45e+69) tmp = t_0; elseif (y <= -3.5e-49) tmp = -6.0 * ((y - x) * z); elseif ((y <= -8.5e-53) || ~((y <= 2.25e+34))) tmp = t_0; else tmp = x * (-3.0 + (z * 6.0)); end tmp_2 = tmp; end
code[x_, y_, z_] := Block[{t$95$0 = N[(y * N[(4.0 + N[(z * -6.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -1.45e+69], t$95$0, If[LessEqual[y, -3.5e-49], N[(-6.0 * N[(N[(y - x), $MachinePrecision] * z), $MachinePrecision]), $MachinePrecision], If[Or[LessEqual[y, -8.5e-53], N[Not[LessEqual[y, 2.25e+34]], $MachinePrecision]], t$95$0, N[(x * N[(-3.0 + N[(z * 6.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := y \cdot \left(4 + z \cdot -6\right)\\
\mathbf{if}\;y \leq -1.45 \cdot 10^{+69}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;y \leq -3.5 \cdot 10^{-49}:\\
\;\;\;\;-6 \cdot \left(\left(y - x\right) \cdot z\right)\\
\mathbf{elif}\;y \leq -8.5 \cdot 10^{-53} \lor \neg \left(y \leq 2.25 \cdot 10^{+34}\right):\\
\;\;\;\;t\_0\\
\mathbf{else}:\\
\;\;\;\;x \cdot \left(-3 + z \cdot 6\right)\\
\end{array}
\end{array}
if y < -1.4499999999999999e69 or -3.50000000000000006e-49 < y < -8.50000000000000044e-53 or 2.25e34 < y Initial program 99.6%
+-commutative99.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 84.5%
if -1.4499999999999999e69 < y < -3.50000000000000006e-49Initial program 99.4%
metadata-eval99.4%
Simplified99.4%
Taylor expanded in z around 0 99.7%
Taylor expanded in z around inf 75.0%
if -8.50000000000000044e-53 < y < 2.25e34Initial program 99.5%
metadata-eval99.5%
Simplified99.5%
Taylor expanded in x around inf 76.6%
sub-neg76.6%
distribute-lft-in76.6%
metadata-eval76.6%
metadata-eval76.6%
neg-mul-176.6%
*-commutative76.6%
associate-*l*76.6%
distribute-rgt-in76.6%
distribute-lft-in76.6%
metadata-eval76.6%
metadata-eval76.6%
associate-+r+76.6%
metadata-eval76.6%
metadata-eval76.6%
associate-*r*76.6%
metadata-eval76.6%
*-commutative76.6%
Simplified76.6%
Final simplification80.2%
(FPCore (x y z)
:precision binary64
(let* ((t_0 (* -6.0 (* y z))))
(if (<= z -1360000000.0)
t_0
(if (<= z 0.5) (* x -3.0) (if (<= z 7.6e+164) (* 6.0 (* x z)) t_0)))))
double code(double x, double y, double z) {
double t_0 = -6.0 * (y * z);
double tmp;
if (z <= -1360000000.0) {
tmp = t_0;
} else if (z <= 0.5) {
tmp = x * -3.0;
} else if (z <= 7.6e+164) {
tmp = 6.0 * (x * z);
} 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 <= (-1360000000.0d0)) then
tmp = t_0
else if (z <= 0.5d0) then
tmp = x * (-3.0d0)
else if (z <= 7.6d+164) then
tmp = 6.0d0 * (x * z)
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 <= -1360000000.0) {
tmp = t_0;
} else if (z <= 0.5) {
tmp = x * -3.0;
} else if (z <= 7.6e+164) {
tmp = 6.0 * (x * z);
} else {
tmp = t_0;
}
return tmp;
}
def code(x, y, z): t_0 = -6.0 * (y * z) tmp = 0 if z <= -1360000000.0: tmp = t_0 elif z <= 0.5: tmp = x * -3.0 elif z <= 7.6e+164: tmp = 6.0 * (x * z) else: tmp = t_0 return tmp
function code(x, y, z) t_0 = Float64(-6.0 * Float64(y * z)) tmp = 0.0 if (z <= -1360000000.0) tmp = t_0; elseif (z <= 0.5) tmp = Float64(x * -3.0); elseif (z <= 7.6e+164) tmp = Float64(6.0 * Float64(x * z)); 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 <= -1360000000.0) tmp = t_0; elseif (z <= 0.5) tmp = x * -3.0; elseif (z <= 7.6e+164) tmp = 6.0 * (x * z); 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, -1360000000.0], t$95$0, If[LessEqual[z, 0.5], N[(x * -3.0), $MachinePrecision], If[LessEqual[z, 7.6e+164], N[(6.0 * N[(x * z), $MachinePrecision]), $MachinePrecision], t$95$0]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := -6 \cdot \left(y \cdot z\right)\\
\mathbf{if}\;z \leq -1360000000:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;z \leq 0.5:\\
\;\;\;\;x \cdot -3\\
\mathbf{elif}\;z \leq 7.6 \cdot 10^{+164}:\\
\;\;\;\;6 \cdot \left(x \cdot z\right)\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
if z < -1.36e9 or 7.60000000000000042e164 < z Initial program 99.8%
metadata-eval99.8%
Simplified99.8%
Taylor expanded in z around 0 99.8%
Taylor expanded in z around inf 99.7%
Taylor expanded in y around inf 62.0%
*-commutative62.0%
Simplified62.0%
if -1.36e9 < z < 0.5Initial program 99.3%
metadata-eval99.3%
Simplified99.3%
Taylor expanded in x around inf 48.8%
sub-neg48.8%
distribute-lft-in48.8%
metadata-eval48.8%
metadata-eval48.8%
neg-mul-148.8%
*-commutative48.8%
associate-*l*48.8%
distribute-rgt-in48.8%
distribute-lft-in48.8%
metadata-eval48.8%
metadata-eval48.8%
associate-+r+48.8%
metadata-eval48.8%
metadata-eval48.8%
associate-*r*48.8%
metadata-eval48.8%
*-commutative48.8%
Simplified48.8%
Taylor expanded in z around 0 46.0%
*-commutative46.0%
Simplified46.0%
if 0.5 < z < 7.60000000000000042e164Initial program 99.5%
metadata-eval99.5%
Simplified99.5%
Taylor expanded in z around 0 99.4%
Taylor expanded in z around inf 94.4%
Taylor expanded in y around 0 62.7%
Final simplification54.0%
(FPCore (x y z)
:precision binary64
(let* ((t_0 (* -6.0 (* y z))))
(if (<= z -1360000000.0)
t_0
(if (<= z 0.5) (* x -3.0) (if (<= z 6.2e+166) (* x (* z 6.0)) t_0)))))
double code(double x, double y, double z) {
double t_0 = -6.0 * (y * z);
double tmp;
if (z <= -1360000000.0) {
tmp = t_0;
} else if (z <= 0.5) {
tmp = x * -3.0;
} else if (z <= 6.2e+166) {
tmp = x * (z * 6.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 <= (-1360000000.0d0)) then
tmp = t_0
else if (z <= 0.5d0) then
tmp = x * (-3.0d0)
else if (z <= 6.2d+166) then
tmp = x * (z * 6.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 <= -1360000000.0) {
tmp = t_0;
} else if (z <= 0.5) {
tmp = x * -3.0;
} else if (z <= 6.2e+166) {
tmp = x * (z * 6.0);
} else {
tmp = t_0;
}
return tmp;
}
def code(x, y, z): t_0 = -6.0 * (y * z) tmp = 0 if z <= -1360000000.0: tmp = t_0 elif z <= 0.5: tmp = x * -3.0 elif z <= 6.2e+166: tmp = x * (z * 6.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 <= -1360000000.0) tmp = t_0; elseif (z <= 0.5) tmp = Float64(x * -3.0); elseif (z <= 6.2e+166) tmp = Float64(x * Float64(z * 6.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 <= -1360000000.0) tmp = t_0; elseif (z <= 0.5) tmp = x * -3.0; elseif (z <= 6.2e+166) tmp = x * (z * 6.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, -1360000000.0], t$95$0, If[LessEqual[z, 0.5], N[(x * -3.0), $MachinePrecision], If[LessEqual[z, 6.2e+166], N[(x * N[(z * 6.0), $MachinePrecision]), $MachinePrecision], t$95$0]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := -6 \cdot \left(y \cdot z\right)\\
\mathbf{if}\;z \leq -1360000000:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;z \leq 0.5:\\
\;\;\;\;x \cdot -3\\
\mathbf{elif}\;z \leq 6.2 \cdot 10^{+166}:\\
\;\;\;\;x \cdot \left(z \cdot 6\right)\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
if z < -1.36e9 or 6.19999999999999966e166 < z Initial program 99.8%
metadata-eval99.8%
Simplified99.8%
Taylor expanded in z around 0 99.8%
Taylor expanded in z around inf 99.7%
Taylor expanded in y around inf 62.0%
*-commutative62.0%
Simplified62.0%
if -1.36e9 < z < 0.5Initial program 99.3%
metadata-eval99.3%
Simplified99.3%
Taylor expanded in x around inf 48.8%
sub-neg48.8%
distribute-lft-in48.8%
metadata-eval48.8%
metadata-eval48.8%
neg-mul-148.8%
*-commutative48.8%
associate-*l*48.8%
distribute-rgt-in48.8%
distribute-lft-in48.8%
metadata-eval48.8%
metadata-eval48.8%
associate-+r+48.8%
metadata-eval48.8%
metadata-eval48.8%
associate-*r*48.8%
metadata-eval48.8%
*-commutative48.8%
Simplified48.8%
Taylor expanded in z around 0 46.0%
*-commutative46.0%
Simplified46.0%
if 0.5 < z < 6.19999999999999966e166Initial program 99.5%
metadata-eval99.5%
Simplified99.5%
Taylor expanded in x around inf 67.1%
sub-neg67.1%
distribute-lft-in67.1%
metadata-eval67.1%
metadata-eval67.1%
neg-mul-167.1%
*-commutative67.1%
associate-*l*67.1%
distribute-rgt-in67.1%
distribute-lft-in67.1%
metadata-eval67.1%
metadata-eval67.1%
associate-+r+67.1%
metadata-eval67.1%
metadata-eval67.1%
associate-*r*67.1%
metadata-eval67.1%
*-commutative67.1%
Simplified67.1%
Taylor expanded in z around inf 62.8%
Final simplification54.0%
(FPCore (x y z) :precision binary64 (if (or (<= z -0.59) (not (<= z 0.66))) (* -6.0 (* (- y x) z)) (- x (* 4.0 (- x y)))))
double code(double x, double y, double z) {
double tmp;
if ((z <= -0.59) || !(z <= 0.66)) {
tmp = -6.0 * ((y - x) * z);
} else {
tmp = x - (4.0 * (x - y));
}
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.59d0)) .or. (.not. (z <= 0.66d0))) then
tmp = (-6.0d0) * ((y - x) * z)
else
tmp = x - (4.0d0 * (x - y))
end if
code = tmp
end function
public static double code(double x, double y, double z) {
double tmp;
if ((z <= -0.59) || !(z <= 0.66)) {
tmp = -6.0 * ((y - x) * z);
} else {
tmp = x - (4.0 * (x - y));
}
return tmp;
}
def code(x, y, z): tmp = 0 if (z <= -0.59) or not (z <= 0.66): tmp = -6.0 * ((y - x) * z) else: tmp = x - (4.0 * (x - y)) return tmp
function code(x, y, z) tmp = 0.0 if ((z <= -0.59) || !(z <= 0.66)) tmp = Float64(-6.0 * Float64(Float64(y - x) * z)); else tmp = Float64(x - Float64(4.0 * Float64(x - y))); end return tmp end
function tmp_2 = code(x, y, z) tmp = 0.0; if ((z <= -0.59) || ~((z <= 0.66))) tmp = -6.0 * ((y - x) * z); else tmp = x - (4.0 * (x - y)); end tmp_2 = tmp; end
code[x_, y_, z_] := If[Or[LessEqual[z, -0.59], N[Not[LessEqual[z, 0.66]], $MachinePrecision]], N[(-6.0 * N[(N[(y - x), $MachinePrecision] * z), $MachinePrecision]), $MachinePrecision], N[(x - N[(4.0 * N[(x - y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;z \leq -0.59 \lor \neg \left(z \leq 0.66\right):\\
\;\;\;\;-6 \cdot \left(\left(y - x\right) \cdot z\right)\\
\mathbf{else}:\\
\;\;\;\;x - 4 \cdot \left(x - y\right)\\
\end{array}
\end{array}
if z < -0.589999999999999969 or 0.660000000000000031 < z Initial program 99.7%
metadata-eval99.7%
Simplified99.7%
Taylor expanded in z around 0 99.7%
Taylor expanded in z around inf 97.8%
if -0.589999999999999969 < z < 0.660000000000000031Initial program 99.3%
metadata-eval99.3%
Simplified99.3%
Taylor expanded in z around 0 97.7%
Final simplification97.7%
(FPCore (x y z) :precision binary64 (if (or (<= z -1360000000.0) (not (<= z 0.65))) (* -6.0 (* y z)) (* x -3.0)))
double code(double x, double y, double z) {
double tmp;
if ((z <= -1360000000.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 <= (-1360000000.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 <= -1360000000.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 <= -1360000000.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 <= -1360000000.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 <= -1360000000.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, -1360000000.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 -1360000000 \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 < -1.36e9 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.5%
Taylor expanded in y around inf 56.3%
*-commutative56.3%
Simplified56.3%
if -1.36e9 < z < 0.650000000000000022Initial program 99.3%
metadata-eval99.3%
Simplified99.3%
Taylor expanded in x around inf 48.8%
sub-neg48.8%
distribute-lft-in48.8%
metadata-eval48.8%
metadata-eval48.8%
neg-mul-148.8%
*-commutative48.8%
associate-*l*48.8%
distribute-rgt-in48.8%
distribute-lft-in48.8%
metadata-eval48.8%
metadata-eval48.8%
associate-+r+48.8%
metadata-eval48.8%
metadata-eval48.8%
associate-*r*48.8%
metadata-eval48.8%
*-commutative48.8%
Simplified48.8%
Taylor expanded in z around 0 46.0%
*-commutative46.0%
Simplified46.0%
Final simplification51.1%
(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%
Final simplification99.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 49.8%
sub-neg49.8%
distribute-lft-in49.8%
metadata-eval49.8%
metadata-eval49.8%
neg-mul-149.8%
*-commutative49.8%
associate-*l*49.8%
distribute-rgt-in49.8%
distribute-lft-in49.8%
metadata-eval49.8%
metadata-eval49.8%
associate-+r+49.8%
metadata-eval49.8%
metadata-eval49.8%
associate-*r*49.8%
metadata-eval49.8%
*-commutative49.8%
Simplified49.8%
Taylor expanded in z around 0 24.7%
*-commutative24.7%
Simplified24.7%
Final simplification24.7%
(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.5%
metadata-eval99.5%
Simplified99.5%
Taylor expanded in y around inf 54.3%
associate-*r*54.3%
Simplified54.3%
Taylor expanded in x around inf 2.6%
Final simplification2.6%
herbie shell --seed 2024046
(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))))