\left(x - \frac{y}{z \cdot 3}\right) + \frac{t}{\left(z \cdot 3\right) \cdot y}\begin{array}{l}
\mathbf{if}\;t \le -3.1361300744902736 \cdot 10^{22}:\\
\;\;\;\;x - \left(\frac{\frac{y}{z}}{3} - t \cdot \frac{\frac{1}{z \cdot 3}}{y}\right)\\
\mathbf{elif}\;t \le 5.46097822804766331 \cdot 10^{49}:\\
\;\;\;\;x - \left(\frac{\frac{y}{z}}{3} - \frac{1}{z} \cdot \frac{\frac{t}{3}}{y}\right)\\
\mathbf{else}:\\
\;\;\;\;\left(x - \frac{y}{z \cdot 3}\right) + \frac{t}{z \cdot \left(3 \cdot y\right)}\\
\end{array}double f(double x, double y, double z, double t) {
double r1308321 = x;
double r1308322 = y;
double r1308323 = z;
double r1308324 = 3.0;
double r1308325 = r1308323 * r1308324;
double r1308326 = r1308322 / r1308325;
double r1308327 = r1308321 - r1308326;
double r1308328 = t;
double r1308329 = r1308325 * r1308322;
double r1308330 = r1308328 / r1308329;
double r1308331 = r1308327 + r1308330;
return r1308331;
}
double f(double x, double y, double z, double t) {
double r1308332 = t;
double r1308333 = -3.1361300744902736e+22;
bool r1308334 = r1308332 <= r1308333;
double r1308335 = x;
double r1308336 = y;
double r1308337 = z;
double r1308338 = r1308336 / r1308337;
double r1308339 = 3.0;
double r1308340 = r1308338 / r1308339;
double r1308341 = 1.0;
double r1308342 = r1308337 * r1308339;
double r1308343 = r1308341 / r1308342;
double r1308344 = r1308343 / r1308336;
double r1308345 = r1308332 * r1308344;
double r1308346 = r1308340 - r1308345;
double r1308347 = r1308335 - r1308346;
double r1308348 = 5.460978228047663e+49;
bool r1308349 = r1308332 <= r1308348;
double r1308350 = r1308341 / r1308337;
double r1308351 = r1308332 / r1308339;
double r1308352 = r1308351 / r1308336;
double r1308353 = r1308350 * r1308352;
double r1308354 = r1308340 - r1308353;
double r1308355 = r1308335 - r1308354;
double r1308356 = r1308336 / r1308342;
double r1308357 = r1308335 - r1308356;
double r1308358 = r1308339 * r1308336;
double r1308359 = r1308337 * r1308358;
double r1308360 = r1308332 / r1308359;
double r1308361 = r1308357 + r1308360;
double r1308362 = r1308349 ? r1308355 : r1308361;
double r1308363 = r1308334 ? r1308347 : r1308362;
return r1308363;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t
Results
| Original | 3.9 |
|---|---|
| Target | 1.7 |
| Herbie | 0.4 |
if t < -3.1361300744902736e+22Initial program 0.5
rmApplied associate-/r*2.7
rmApplied associate-/r*2.7
rmApplied associate-+l-2.7
rmApplied *-un-lft-identity2.7
Applied div-inv2.8
Applied times-frac0.4
Simplified0.4
if -3.1361300744902736e+22 < t < 5.460978228047663e+49Initial program 5.8
rmApplied associate-/r*1.1
rmApplied associate-/r*1.1
rmApplied associate-+l-1.1
rmApplied *-un-lft-identity1.1
Applied *-un-lft-identity1.1
Applied times-frac1.1
Applied times-frac0.3
Simplified0.3
if 5.460978228047663e+49 < t Initial program 0.6
rmApplied associate-*l*0.6
Final simplification0.4
herbie shell --seed 2020047
(FPCore (x y z t)
:name "Diagrams.Solve.Polynomial:cubForm from diagrams-solve-0.1, H"
:precision binary64
:herbie-target
(+ (- x (/ y (* z 3))) (/ (/ t (* z 3)) y))
(+ (- x (/ y (* z 3))) (/ t (* (* z 3) y))))