x + \frac{\left(y - z\right) \cdot \left(t - x\right)}{a - z}\begin{array}{l}
\mathbf{if}\;z \le -6.8292046976100912 \cdot 10^{155} \lor \neg \left(z \le 7.37344260106273469 \cdot 10^{217}\right):\\
\;\;\;\;y \cdot \left(\frac{x}{z} - \frac{t}{z}\right) + t\\
\mathbf{else}:\\
\;\;\;\;x + \frac{\sqrt[3]{y - z} \cdot \sqrt[3]{y - z}}{\sqrt[3]{a - z}} \cdot \left(\frac{\sqrt[3]{y - z}}{\sqrt[3]{a - z}} \cdot \frac{t - x}{\sqrt[3]{a - z}}\right)\\
\end{array}double f(double x, double y, double z, double t, double a) {
double r609000 = x;
double r609001 = y;
double r609002 = z;
double r609003 = r609001 - r609002;
double r609004 = t;
double r609005 = r609004 - r609000;
double r609006 = r609003 * r609005;
double r609007 = a;
double r609008 = r609007 - r609002;
double r609009 = r609006 / r609008;
double r609010 = r609000 + r609009;
return r609010;
}
double f(double x, double y, double z, double t, double a) {
double r609011 = z;
double r609012 = -6.829204697610091e+155;
bool r609013 = r609011 <= r609012;
double r609014 = 7.373442601062735e+217;
bool r609015 = r609011 <= r609014;
double r609016 = !r609015;
bool r609017 = r609013 || r609016;
double r609018 = y;
double r609019 = x;
double r609020 = r609019 / r609011;
double r609021 = t;
double r609022 = r609021 / r609011;
double r609023 = r609020 - r609022;
double r609024 = r609018 * r609023;
double r609025 = r609024 + r609021;
double r609026 = r609018 - r609011;
double r609027 = cbrt(r609026);
double r609028 = r609027 * r609027;
double r609029 = a;
double r609030 = r609029 - r609011;
double r609031 = cbrt(r609030);
double r609032 = r609028 / r609031;
double r609033 = r609027 / r609031;
double r609034 = r609021 - r609019;
double r609035 = r609034 / r609031;
double r609036 = r609033 * r609035;
double r609037 = r609032 * r609036;
double r609038 = r609019 + r609037;
double r609039 = r609017 ? r609025 : r609038;
return r609039;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t




Bits error versus a
Results
| Original | 24.5 |
|---|---|
| Target | 11.5 |
| Herbie | 10.0 |
if z < -6.829204697610091e+155 or 7.373442601062735e+217 < z Initial program 49.2
rmApplied add-cube-cbrt49.5
Applied times-frac24.5
rmApplied add-cube-cbrt24.1
Applied times-frac24.1
Applied associate-*l*24.1
Taylor expanded around inf 24.4
Simplified15.1
if -6.829204697610091e+155 < z < 7.373442601062735e+217Initial program 17.1
rmApplied add-cube-cbrt17.6
Applied times-frac8.9
rmApplied add-cube-cbrt8.8
Applied times-frac8.9
Applied associate-*l*8.4
Final simplification10.0
herbie shell --seed 2020081
(FPCore (x y z t a)
:name "Graphics.Rendering.Chart.Axis.Types:invLinMap from Chart-1.5.3"
:precision binary64
:herbie-target
(if (< z -1.2536131056095036e+188) (- t (* (/ y z) (- t x))) (if (< z 4.446702369113811e+64) (+ x (/ (- y z) (/ (- a z) (- t x)))) (- t (* (/ y z) (- t x)))))
(+ x (/ (* (- y z) (- t x)) (- a z))))