\tan^{-1} \left(N + 1\right) - \tan^{-1} N\tan^{-1}_* \frac{1}{1 + \left(\sqrt[3]{N + 1} \cdot \sqrt[3]{N + 1}\right) \cdot \left(\left(\sqrt[3]{N \cdot N - 1 \cdot 1} \cdot \frac{1}{\sqrt[3]{N - 1}}\right) \cdot N\right)}double f(double N) {
double r138440 = N;
double r138441 = 1.0;
double r138442 = r138440 + r138441;
double r138443 = atan(r138442);
double r138444 = atan(r138440);
double r138445 = r138443 - r138444;
return r138445;
}
double f(double N) {
double r138446 = 1.0;
double r138447 = 1.0;
double r138448 = N;
double r138449 = r138448 + r138446;
double r138450 = cbrt(r138449);
double r138451 = r138450 * r138450;
double r138452 = r138448 * r138448;
double r138453 = r138446 * r138446;
double r138454 = r138452 - r138453;
double r138455 = cbrt(r138454);
double r138456 = r138448 - r138446;
double r138457 = cbrt(r138456);
double r138458 = r138447 / r138457;
double r138459 = r138455 * r138458;
double r138460 = r138459 * r138448;
double r138461 = r138451 * r138460;
double r138462 = r138447 + r138461;
double r138463 = atan2(r138446, r138462);
return r138463;
}




Bits error versus N
Results
| Original | 14.9 |
|---|---|
| Target | 0.3 |
| Herbie | 0.5 |
Initial program 14.9
rmApplied diff-atan13.8
Simplified0.3
rmApplied add-cube-cbrt0.6
Applied associate-*l*0.6
rmApplied flip-+0.6
Applied cbrt-div0.5
rmApplied div-inv0.5
Final simplification0.5
herbie shell --seed 2020047
(FPCore (N)
:name "2atan (example 3.5)"
:precision binary64
:herbie-target
(atan (/ 1 (+ 1 (* N (+ N 1)))))
(- (atan (+ N 1)) (atan N)))