\frac{\left(\left(x \cdot 9\right) \cdot y - \left(\left(z \cdot 4\right) \cdot t\right) \cdot a\right) + b}{z \cdot c}\begin{array}{l}
\mathbf{if}\;c \le -1.1394805299673187 \cdot 10^{46} \lor \neg \left(c \le 4.2411716606081791 \cdot 10^{-37}\right):\\
\;\;\;\;\mathsf{fma}\left(-4, a \cdot \frac{t}{c}, \frac{\frac{\mathsf{fma}\left(x, 9 \cdot y, b\right)}{c}}{z}\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(-4, \frac{t \cdot a}{c}, \frac{\frac{\mathsf{fma}\left(x, 9 \cdot y, b\right)}{z}}{c}\right)\\
\end{array}double code(double x, double y, double z, double t, double a, double b, double c) {
return ((double) (((double) (((double) (((double) (((double) (x * 9.0)) * y)) - ((double) (((double) (((double) (z * 4.0)) * t)) * a)))) + b)) / ((double) (z * c))));
}
double code(double x, double y, double z, double t, double a, double b, double c) {
double VAR;
if (((c <= -1.1394805299673187e+46) || !(c <= 4.241171660608179e-37))) {
VAR = ((double) fma(((double) -(4.0)), ((double) (a * ((double) (t / c)))), ((double) (((double) (((double) fma(x, ((double) (9.0 * y)), b)) / c)) / z))));
} else {
VAR = ((double) fma(((double) -(4.0)), ((double) (((double) (t * a)) / c)), ((double) (((double) (((double) fma(x, ((double) (9.0 * y)), b)) / z)) / c))));
}
return VAR;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t




Bits error versus a




Bits error versus b




Bits error versus c
Results
| Original | 20.5 |
|---|---|
| Target | 14.4 |
| Herbie | 6.9 |
if c < -1.1394805299673187e+46 or 4.241171660608179e-37 < c Initial program 23.4
Simplified15.0
rmApplied *-commutative15.0
Applied associate-/l*11.9
rmApplied *-commutative11.9
Applied associate-/r*8.7
rmApplied clear-num8.8
Applied associate-/r/8.6
Simplified8.6
if -1.1394805299673187e+46 < c < 4.241171660608179e-37Initial program 14.7
Simplified5.9
rmApplied associate-/r*3.6
Final simplification6.9
herbie shell --seed 2020114 +o rules:numerics
(FPCore (x y z t a b c)
:name "Diagrams.Solve.Polynomial:cubForm from diagrams-solve-0.1, J"
:precision binary64
:herbie-target
(if (< (/ (+ (- (* (* x 9) y) (* (* (* z 4) t) a)) b) (* z c)) -1.1001567408041051e-171) (/ (+ (- (* (* x 9) y) (* (* z 4) (* t a))) b) (* z c)) (if (< (/ (+ (- (* (* x 9) y) (* (* (* z 4) t) a)) b) (* z c)) -0.0) (/ (/ (+ (- (* (* x 9) y) (* (* (* z 4) t) a)) b) z) c) (if (< (/ (+ (- (* (* x 9) y) (* (* (* z 4) t) a)) b) (* z c)) 1.1708877911747488e-53) (/ (+ (- (* (* x 9) y) (* (* z 4) (* t a))) b) (* z c)) (if (< (/ (+ (- (* (* x 9) y) (* (* (* z 4) t) a)) b) (* z c)) 2.876823679546137e+130) (- (+ (* (* 9 (/ y c)) (/ x z)) (/ b (* c z))) (* 4 (/ (* a t) c))) (if (< (/ (+ (- (* (* x 9) y) (* (* (* z 4) t) a)) b) (* z c)) 1.3838515042456319e+158) (/ (+ (- (* (* x 9) y) (* (* z 4) (* t a))) b) (* z c)) (- (+ (* 9 (* (/ y (* c z)) x)) (/ b (* c z))) (* 4 (/ (* a t) c))))))))
(/ (+ (- (* (* x 9) y) (* (* (* z 4) t) a)) b) (* z c)))