x + \frac{\left(y - x\right) \cdot \left(z - t\right)}{a - t}\begin{array}{l}
\mathbf{if}\;a \le -1.903940587482351271448540463407352524949 \cdot 10^{-94}:\\
\;\;\;\;x + \frac{\frac{1}{\frac{\sqrt[3]{a - t}}{\sqrt[3]{y - x} \cdot \sqrt[3]{z - t}} \cdot \frac{\sqrt[3]{a - t}}{\left(\sqrt[3]{\sqrt[3]{y - x}} \cdot \sqrt[3]{\sqrt[3]{y - x} \cdot \sqrt[3]{y - x}}\right) \cdot \sqrt[3]{z - t}}}}{\frac{\frac{\sqrt[3]{a - t}}{\sqrt[3]{z - t}}}{\sqrt[3]{y - x}}}\\
\mathbf{elif}\;a \le 2.486579097351087244398754404400467109927 \cdot 10^{-154}:\\
\;\;\;\;\left(y + \frac{z \cdot x}{t}\right) - \frac{z \cdot y}{t}\\
\mathbf{else}:\\
\;\;\;\;\frac{1}{\frac{\frac{a - t}{z - t}}{y - x}} + x\\
\end{array}double f(double x, double y, double z, double t, double a) {
double r35484476 = x;
double r35484477 = y;
double r35484478 = r35484477 - r35484476;
double r35484479 = z;
double r35484480 = t;
double r35484481 = r35484479 - r35484480;
double r35484482 = r35484478 * r35484481;
double r35484483 = a;
double r35484484 = r35484483 - r35484480;
double r35484485 = r35484482 / r35484484;
double r35484486 = r35484476 + r35484485;
return r35484486;
}
double f(double x, double y, double z, double t, double a) {
double r35484487 = a;
double r35484488 = -1.9039405874823513e-94;
bool r35484489 = r35484487 <= r35484488;
double r35484490 = x;
double r35484491 = 1.0;
double r35484492 = t;
double r35484493 = r35484487 - r35484492;
double r35484494 = cbrt(r35484493);
double r35484495 = y;
double r35484496 = r35484495 - r35484490;
double r35484497 = cbrt(r35484496);
double r35484498 = z;
double r35484499 = r35484498 - r35484492;
double r35484500 = cbrt(r35484499);
double r35484501 = r35484497 * r35484500;
double r35484502 = r35484494 / r35484501;
double r35484503 = cbrt(r35484497);
double r35484504 = r35484497 * r35484497;
double r35484505 = cbrt(r35484504);
double r35484506 = r35484503 * r35484505;
double r35484507 = r35484506 * r35484500;
double r35484508 = r35484494 / r35484507;
double r35484509 = r35484502 * r35484508;
double r35484510 = r35484491 / r35484509;
double r35484511 = r35484494 / r35484500;
double r35484512 = r35484511 / r35484497;
double r35484513 = r35484510 / r35484512;
double r35484514 = r35484490 + r35484513;
double r35484515 = 2.4865790973510872e-154;
bool r35484516 = r35484487 <= r35484515;
double r35484517 = r35484498 * r35484490;
double r35484518 = r35484517 / r35484492;
double r35484519 = r35484495 + r35484518;
double r35484520 = r35484498 * r35484495;
double r35484521 = r35484520 / r35484492;
double r35484522 = r35484519 - r35484521;
double r35484523 = r35484493 / r35484499;
double r35484524 = r35484523 / r35484496;
double r35484525 = r35484491 / r35484524;
double r35484526 = r35484525 + r35484490;
double r35484527 = r35484516 ? r35484522 : r35484526;
double r35484528 = r35484489 ? r35484514 : r35484527;
return r35484528;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t




Bits error versus a
Results
| Original | 24.3 |
|---|---|
| Target | 9.0 |
| Herbie | 10.3 |
if a < -1.9039405874823513e-94Initial program 22.0
rmApplied associate-/l*7.8
rmApplied clear-num7.9
rmApplied add-cube-cbrt8.4
Applied add-cube-cbrt8.5
Applied add-cube-cbrt8.5
Applied times-frac8.5
Applied times-frac8.1
Applied associate-/r*8.1
Simplified8.1
rmApplied add-cube-cbrt8.2
Applied cbrt-prod8.2
if -1.9039405874823513e-94 < a < 2.4865790973510872e-154Initial program 29.5
Taylor expanded around inf 14.5
if 2.4865790973510872e-154 < a Initial program 22.9
rmApplied associate-/l*9.3
rmApplied clear-num9.4
Final simplification10.3
herbie shell --seed 2019171
(FPCore (x y z t a)
:name "Graphics.Rendering.Chart.Axis.Types:linMap from Chart-1.5.3"
:herbie-target
(if (< a -1.6153062845442575e-142) (+ x (* (/ (- y x) 1.0) (/ (- z t) (- a t)))) (if (< a 3.774403170083174e-182) (- y (* (/ z t) (- y x))) (+ x (* (/ (- y x) 1.0) (/ (- z t) (- a t))))))
(+ x (/ (* (- y x) (- z t)) (- a t))))