\frac{\left(x \cdot x + y \cdot y\right) - z \cdot z}{y \cdot 2}\frac{y + \left(x + z\right) \cdot \frac{x - z}{y}}{2}double f(double x, double y, double z) {
double r550152 = x;
double r550153 = r550152 * r550152;
double r550154 = y;
double r550155 = r550154 * r550154;
double r550156 = r550153 + r550155;
double r550157 = z;
double r550158 = r550157 * r550157;
double r550159 = r550156 - r550158;
double r550160 = 2.0;
double r550161 = r550154 * r550160;
double r550162 = r550159 / r550161;
return r550162;
}
double f(double x, double y, double z) {
double r550163 = y;
double r550164 = x;
double r550165 = z;
double r550166 = r550164 + r550165;
double r550167 = r550164 - r550165;
double r550168 = r550167 / r550163;
double r550169 = r550166 * r550168;
double r550170 = r550163 + r550169;
double r550171 = 2.0;
double r550172 = r550170 / r550171;
return r550172;
}




Bits error versus x




Bits error versus y




Bits error versus z
Results
| Original | 28.2 |
|---|---|
| Target | 0.2 |
| Herbie | 0.1 |
Initial program 28.2
Simplified12.2
rmApplied *-un-lft-identity12.2
Applied difference-of-squares12.2
Applied times-frac0.1
Simplified0.1
Final simplification0.1
herbie shell --seed 2019212
(FPCore (x y z)
:name "Diagrams.TwoD.Apollonian:initialConfig from diagrams-contrib-1.3.0.5, A"
:precision binary64
:herbie-target
(- (* y 0.5) (* (* (/ 0.5 y) (+ z x)) (- z x)))
(/ (- (+ (* x x) (* y y)) (* z z)) (* y 2)))