\left|\left(\left(\tan^{-1}_* \frac{\mathsf{expm1}\left(\sin \left(\mathsf{expm1}\left(a\right)\right)\right)}{\tan^{-1} a}\right) \bmod a\right)\right|\left|\left(\sqrt[3]{\left(\left(\left(\sqrt[3]{\tan^{-1}_* \frac{\mathsf{expm1}\left(\sin \left(\mathsf{expm1}\left(a\right)\right)\right)}{\tan^{-1} a}} \cdot \sqrt[3]{\tan^{-1}_* \frac{\mathsf{expm1}\left(\sin \left(\mathsf{expm1}\left(a\right)\right)\right)}{\tan^{-1} a}}\right) \cdot \sqrt[3]{\tan^{-1}_* \frac{\mathsf{expm1}\left(\sin \left(\mathsf{expm1}\left(a\right)\right)\right)}{\tan^{-1} a}}\right) \bmod a\right)} \cdot \sqrt[3]{\left(\left(\tan^{-1}_* \frac{\mathsf{expm1}\left(\sin \left(\mathsf{expm1}\left(a\right)\right)\right)}{\tan^{-1} a}\right) \bmod a\right)}\right) \cdot \sqrt[3]{\left(\left(\left(\sqrt[3]{\tan^{-1}_* \frac{\mathsf{expm1}\left(\sin \left(\mathsf{expm1}\left(a\right)\right)\right)}{\tan^{-1} a}} \cdot \sqrt[3]{\tan^{-1}_* \frac{\mathsf{expm1}\left(\sin \left(\mathsf{expm1}\left(a\right)\right)\right)}{\tan^{-1} a}}\right) \cdot \left(\left(\sqrt[3]{\sqrt[3]{\tan^{-1}_* \frac{\mathsf{expm1}\left(\sin \left(\mathsf{expm1}\left(a\right)\right)\right)}{\tan^{-1} a}}} \cdot \sqrt[3]{\sqrt[3]{\tan^{-1}_* \frac{\mathsf{expm1}\left(\sin \left(\mathsf{expm1}\left(a\right)\right)\right)}{\tan^{-1} a}}}\right) \cdot \sqrt[3]{\sqrt[3]{\tan^{-1}_* \frac{\mathsf{expm1}\left(\sin \left(\mathsf{expm1}\left(a\right)\right)\right)}{\tan^{-1} a}}}\right)\right) \bmod a\right)}\right|double f(double a) {
double r12448 = a;
double r12449 = expm1(r12448);
double r12450 = sin(r12449);
double r12451 = expm1(r12450);
double r12452 = atan(r12448);
double r12453 = atan2(r12451, r12452);
double r12454 = fmod(r12453, r12448);
double r12455 = fabs(r12454);
return r12455;
}
double f(double a) {
double r12456 = a;
double r12457 = expm1(r12456);
double r12458 = sin(r12457);
double r12459 = expm1(r12458);
double r12460 = atan(r12456);
double r12461 = atan2(r12459, r12460);
double r12462 = cbrt(r12461);
double r12463 = r12462 * r12462;
double r12464 = r12463 * r12462;
double r12465 = fmod(r12464, r12456);
double r12466 = cbrt(r12465);
double r12467 = fmod(r12461, r12456);
double r12468 = cbrt(r12467);
double r12469 = r12466 * r12468;
double r12470 = cbrt(r12462);
double r12471 = r12470 * r12470;
double r12472 = r12471 * r12470;
double r12473 = r12463 * r12472;
double r12474 = fmod(r12473, r12456);
double r12475 = cbrt(r12474);
double r12476 = r12469 * r12475;
double r12477 = fabs(r12476);
return r12477;
}



Bits error versus a
Initial program 33.3
rmApplied add-cube-cbrt33.4
rmApplied add-cube-cbrt33.4
Taylor expanded around 0 33.4
rmApplied add-cube-cbrt33.4
Final simplification33.4
herbie shell --seed 2020065 +o rules:numerics
(FPCore (a)
:name "Random Jason Timeout Test 006"
:precision binary64
(fabs (fmod (atan2 (expm1 (sin (expm1 a))) (atan a)) a)))