\frac{6 \cdot \left(x - 1\right)}{\left(x + 1\right) + 4 \cdot \sqrt{x}}\frac{\sqrt{x} + \sqrt{1}}{\sqrt{\mathsf{fma}\left(\sqrt{x}, 4, x + 1\right)}} \cdot \left(\frac{\sqrt{x} - \sqrt{1}}{\sqrt{\mathsf{fma}\left(\sqrt{x}, 4, x + 1\right)}} \cdot 6\right)double f(double x) {
double r1072372 = 6.0;
double r1072373 = x;
double r1072374 = 1.0;
double r1072375 = r1072373 - r1072374;
double r1072376 = r1072372 * r1072375;
double r1072377 = r1072373 + r1072374;
double r1072378 = 4.0;
double r1072379 = sqrt(r1072373);
double r1072380 = r1072378 * r1072379;
double r1072381 = r1072377 + r1072380;
double r1072382 = r1072376 / r1072381;
return r1072382;
}
double f(double x) {
double r1072383 = x;
double r1072384 = sqrt(r1072383);
double r1072385 = 1.0;
double r1072386 = sqrt(r1072385);
double r1072387 = r1072384 + r1072386;
double r1072388 = 4.0;
double r1072389 = r1072383 + r1072385;
double r1072390 = fma(r1072384, r1072388, r1072389);
double r1072391 = sqrt(r1072390);
double r1072392 = r1072387 / r1072391;
double r1072393 = r1072384 - r1072386;
double r1072394 = r1072393 / r1072391;
double r1072395 = 6.0;
double r1072396 = r1072394 * r1072395;
double r1072397 = r1072392 * r1072396;
return r1072397;
}




Bits error versus x
| Original | 0.2 |
|---|---|
| Target | 0.1 |
| Herbie | 0.1 |
Initial program 0.2
Simplified0.0
rmApplied *-un-lft-identity0.0
Applied add-sqr-sqrt0.3
Applied times-frac0.3
Applied add-sqr-sqrt0.3
Applied add-sqr-sqrt0.2
Applied difference-of-squares0.2
Applied times-frac0.1
Simplified0.1
rmApplied associate-/r/0.1
Final simplification0.1
herbie shell --seed 2019354 +o rules:numerics
(FPCore (x)
:name "Data.Approximate.Numerics:blog from approximate-0.2.2.1"
:precision binary64
:herbie-target
(/ 6 (/ (+ (+ x 1) (* 4 (sqrt x))) (- x 1)))
(/ (* 6 (- x 1)) (+ (+ x 1) (* 4 (sqrt x)))))