\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(\sqrt[3]{{\left(\tan^{-1}_* \frac{\mathsf{expm1}\left(\sin \left(\mathsf{expm1}\left(a\right)\right)\right)}{\tan^{-1} a}\right)}^{3}}\right) \bmod a\right)} \cdot \sqrt[3]{\left(\left(\sqrt[3]{{\left(\tan^{-1}_* \frac{\mathsf{expm1}\left(\sin \left(\mathsf{expm1}\left(a\right)\right)\right)}{\tan^{-1} a}\right)}^{3}}\right) \bmod a\right)}\right) \cdot \sqrt[3]{\left(\left(\sqrt[3]{\tan^{-1}_* \frac{\mathsf{expm1}\left(\sin \left(\mathsf{expm1}\left(a\right)\right)\right)}{\tan^{-1} a} \cdot {\left(\tan^{-1}_* \frac{\mathsf{expm1}\left(\sin \left(\mathsf{expm1}\left(a\right)\right)\right)}{\tan^{-1} a}\right)}^{2}}\right) \bmod a\right)}\right|double f(double a) {
double r14485 = a;
double r14486 = expm1(r14485);
double r14487 = sin(r14486);
double r14488 = expm1(r14487);
double r14489 = atan(r14485);
double r14490 = atan2(r14488, r14489);
double r14491 = fmod(r14490, r14485);
double r14492 = fabs(r14491);
return r14492;
}
double f(double a) {
double r14493 = a;
double r14494 = expm1(r14493);
double r14495 = sin(r14494);
double r14496 = expm1(r14495);
double r14497 = atan(r14493);
double r14498 = atan2(r14496, r14497);
double r14499 = 3.0;
double r14500 = pow(r14498, r14499);
double r14501 = cbrt(r14500);
double r14502 = fmod(r14501, r14493);
double r14503 = cbrt(r14502);
double r14504 = r14503 * r14503;
double r14505 = 2.0;
double r14506 = pow(r14498, r14505);
double r14507 = r14498 * r14506;
double r14508 = cbrt(r14507);
double r14509 = fmod(r14508, r14493);
double r14510 = cbrt(r14509);
double r14511 = r14504 * r14510;
double r14512 = fabs(r14511);
return r14512;
}



Bits error versus a
Initial program 33.7
rmApplied add-cbrt-cube33.7
Simplified33.7
rmApplied add-cube-cbrt33.7
rmApplied cube-mult33.7
Simplified33.7
Final simplification33.7
herbie shell --seed 2019294
(FPCore (a)
:name "Random Jason Timeout Test 006"
:precision binary64
(fabs (fmod (atan2 (expm1 (sin (expm1 a))) (atan a)) a)))