\frac{b \cdot c - a \cdot d}{c \cdot c + d \cdot d}(FPCore (a b c d) :precision binary64 (/ (- (* b c) (* a d)) (+ (* c c) (* d d))))
double code(double a, double b, double c, double d) {
return ((b * c) - (a * d)) / ((c * c) + (d * d));
}
Please include this information when filing a bug report:
herbie shell --seed 2021060
(FPCore (a b c d)
:name "Complex division, imag part"
:precision binary64
:herbie-target
(if (< (fabs d) (fabs c)) (/ (- b (* a (/ d c))) (+ c (* d (/ d c)))) (/ (+ (- a) (* b (/ c d))) (+ d (* c (/ c d)))))
(/ (- (* b c) (* a d)) (+ (* c c) (* d d))))
| Failed to produce proof for (/.f64 (*.f64 (cbrt.f64 d) (cbrt.f64 d)) (/.f64 (*.f64 (cbrt.f64 (+.f64 (*.f64 c c) (*.f64 d d))) (cbrt.f64 (+.f64 (*.f64 c c) (*.f64 d d)))) (sqrt.f64 a))) to (/.f64 (*.f64 (cbrt.f64 d) (sqrt.f64 a)) (/.f64 (cbrt.f64 (+.f64 (*.f64 c c) (*.f64 d d))) (/.f64 (cbrt.f64 d) (cbrt.f64 (+.f64 (*.f64 c c) (*.f64 d d)))))) | L | C | |
|---|---|---|---|
| for-loop | /data/pavpan/nightlies/herbie/egg-generate-proof/src/core/simplify.rkt | 167 | 8 |
| (unnamed) | /data/pavpan/nightlies/herbie/egg-generate-proof/src/core/simplify.rkt | 149 | 6 |
| egraph-add-exprs | /data/pavpan/nightlies/herbie/egg-generate-proof/egg-herbie/main.rkt | 154 | 0 |
| with-egraph | /data/pavpan/nightlies/herbie/egg-generate-proof/egg-herbie/main.rkt | 74 | 0 |
| (unnamed) | /data/pavpan/nightlies/herbie/egg-generate-proof/src/core/simplify.rkt | 79 | 0 |
| simplify! | /data/pavpan/nightlies/herbie/egg-generate-proof/src/mainloop.rkt | 236 | 0 |
| run-iter! | /data/pavpan/nightlies/herbie/egg-generate-proof/src/mainloop.rkt | 360 | 0 |
| run-improve | /data/pavpan/nightlies/herbie/egg-generate-proof/src/mainloop.rkt | 377 | 0 |
| (unnamed) | /data/pavpan/nightlies/herbie/egg-generate-proof/src/sandbox.rkt | 54 | 6 |
| run | /opt/racket-7.9/share/pkgs/profile-lib/main.rkt | 39 | 2 |
| profile-thunk | /opt/racket-7.9/share/pkgs/profile-lib/main.rkt | 9 | 0 |
| (unnamed) | /opt/racket-7.9/collects/racket/engine.rkt | 42 | 24 |