\frac{a \cdot c + b \cdot d}{c \cdot c + d \cdot d}\frac{\frac{1}{\frac{\sqrt{c \cdot c + d \cdot d}}{d \cdot b + a \cdot c}}}{\sqrt{c \cdot c + d \cdot d}}double f(double a, double b, double c, double d) {
double r93235 = a;
double r93236 = c;
double r93237 = r93235 * r93236;
double r93238 = b;
double r93239 = d;
double r93240 = r93238 * r93239;
double r93241 = r93237 + r93240;
double r93242 = r93236 * r93236;
double r93243 = r93239 * r93239;
double r93244 = r93242 + r93243;
double r93245 = r93241 / r93244;
return r93245;
}
double f(double a, double b, double c, double d) {
double r93246 = 1.0;
double r93247 = c;
double r93248 = r93247 * r93247;
double r93249 = d;
double r93250 = r93249 * r93249;
double r93251 = r93248 + r93250;
double r93252 = sqrt(r93251);
double r93253 = b;
double r93254 = r93249 * r93253;
double r93255 = a;
double r93256 = r93255 * r93247;
double r93257 = r93254 + r93256;
double r93258 = r93252 / r93257;
double r93259 = r93246 / r93258;
double r93260 = r93259 / r93252;
return r93260;
}




Bits error versus a




Bits error versus b




Bits error versus c




Bits error versus d
Results
| Original | 25.9 |
|---|---|
| Target | 0.4 |
| Herbie | 25.9 |
Initial program 25.9
rmApplied add-sqr-sqrt25.9
Applied associate-/r*25.8
Simplified25.8
rmApplied clear-num25.9
Simplified25.9
Final simplification25.9
herbie shell --seed 2019198
(FPCore (a b c d)
:name "Complex division, real part"
:herbie-target
(if (< (fabs d) (fabs c)) (/ (+ a (* b (/ d c))) (+ c (* d (/ d c)))) (/ (+ b (* a (/ c d))) (+ d (* c (/ c d)))))
(/ (+ (* a c) (* b d)) (+ (* c c) (* d d))))