\sqrt{\left(x \cdot x + y \cdot y\right) + z \cdot z}\begin{array}{l}
\mathbf{if}\;x \le -7.97378979562191825 \cdot 10^{148}:\\
\;\;\;\;-1 \cdot x\\
\mathbf{elif}\;x \le 7.0809856866698759 \cdot 10^{84}:\\
\;\;\;\;\sqrt{\left(x \cdot x + y \cdot y\right) + z \cdot z}\\
\mathbf{else}:\\
\;\;\;\;x\\
\end{array}double code(double x, double y, double z) {
return sqrt((((x * x) + (y * y)) + (z * z)));
}
double code(double x, double y, double z) {
double temp;
if ((x <= -7.973789795621918e+148)) {
temp = (-1.0 * x);
} else {
double temp_1;
if ((x <= 7.080985686669876e+84)) {
temp_1 = sqrt((((x * x) + (y * y)) + (z * z)));
} else {
temp_1 = x;
}
temp = temp_1;
}
return temp;
}




Bits error versus x




Bits error versus y




Bits error versus z
Results
| Original | 37.7 |
|---|---|
| Target | 25.1 |
| Herbie | 25.4 |
if x < -7.973789795621918e+148Initial program 63.2
Taylor expanded around -inf 14.1
if -7.973789795621918e+148 < x < 7.080985686669876e+84Initial program 29.2
if 7.080985686669876e+84 < x Initial program 52.9
Taylor expanded around inf 18.8
Final simplification25.4
herbie shell --seed 2020057
(FPCore (x y z)
:name "FRP.Yampa.Vector3:vector3Rho from Yampa-0.10.2"
:precision binary64
:herbie-target
(if (< z -6.396479394109776e+136) (- z) (if (< z 7.320293694404182e+117) (sqrt (+ (+ (* z z) (* x x)) (* y y))) z))
(sqrt (+ (+ (* x x) (* y y)) (* z z))))