x + \frac{y \cdot \left(z - x\right)}{t}\begin{array}{l}
\mathbf{if}\;t \le -4.024110690636155670064376141757827793703 \cdot 10^{-54}:\\
\;\;\;\;x + y \cdot \left(\frac{z}{t} - \frac{x}{t}\right)\\
\mathbf{elif}\;t \le 4.482541642658298470975327978871003617787 \cdot 10^{-266}:\\
\;\;\;\;x + \frac{1}{\frac{t}{y \cdot \left(z - x\right)}}\\
\mathbf{else}:\\
\;\;\;\;x + \frac{\sqrt[3]{y} \cdot \sqrt[3]{y}}{\sqrt[3]{t}} \cdot \left(\frac{\sqrt[3]{y}}{\sqrt[3]{t}} \cdot \frac{z - x}{\sqrt[3]{t}}\right)\\
\end{array}double f(double x, double y, double z, double t) {
double r266194 = x;
double r266195 = y;
double r266196 = z;
double r266197 = r266196 - r266194;
double r266198 = r266195 * r266197;
double r266199 = t;
double r266200 = r266198 / r266199;
double r266201 = r266194 + r266200;
return r266201;
}
double f(double x, double y, double z, double t) {
double r266202 = t;
double r266203 = -4.0241106906361557e-54;
bool r266204 = r266202 <= r266203;
double r266205 = x;
double r266206 = y;
double r266207 = z;
double r266208 = r266207 / r266202;
double r266209 = r266205 / r266202;
double r266210 = r266208 - r266209;
double r266211 = r266206 * r266210;
double r266212 = r266205 + r266211;
double r266213 = 4.4825416426582985e-266;
bool r266214 = r266202 <= r266213;
double r266215 = 1.0;
double r266216 = r266207 - r266205;
double r266217 = r266206 * r266216;
double r266218 = r266202 / r266217;
double r266219 = r266215 / r266218;
double r266220 = r266205 + r266219;
double r266221 = cbrt(r266206);
double r266222 = r266221 * r266221;
double r266223 = cbrt(r266202);
double r266224 = r266222 / r266223;
double r266225 = r266221 / r266223;
double r266226 = r266216 / r266223;
double r266227 = r266225 * r266226;
double r266228 = r266224 * r266227;
double r266229 = r266205 + r266228;
double r266230 = r266214 ? r266220 : r266229;
double r266231 = r266204 ? r266212 : r266230;
return r266231;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t
Results
| Original | 6.4 |
|---|---|
| Target | 2.2 |
| Herbie | 1.8 |
if t < -4.0241106906361557e-54Initial program 8.6
rmApplied add-cube-cbrt8.9
Applied times-frac0.9
Taylor expanded around 0 8.6
Simplified1.3
if -4.0241106906361557e-54 < t < 4.4825416426582985e-266Initial program 2.3
rmApplied clear-num2.4
if 4.4825416426582985e-266 < t Initial program 6.2
rmApplied add-cube-cbrt6.7
Applied times-frac2.8
rmApplied add-cube-cbrt2.9
Applied times-frac2.9
Applied associate-*l*2.0
Final simplification1.8
herbie shell --seed 2019323
(FPCore (x y z t)
:name "Optimisation.CirclePacking:place from circle-packing-0.1.0.4, D"
:precision binary64
:herbie-target
(- x (+ (* x (/ y t)) (* (- z) (/ y t))))
(+ x (/ (* y (- z x)) t)))