\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(\tan^{-1}_* \frac{\left(\sqrt[3]{\mathsf{expm1}\left(\sin \left(\mathsf{expm1}\left(a\right)\right)\right)} \cdot \sqrt[3]{\mathsf{expm1}\left(\sin \left(\mathsf{expm1}\left(a\right)\right)\right)}\right) \cdot \sqrt[3]{\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{\left(\sqrt[3]{\mathsf{expm1}\left(\sin \left(\mathsf{expm1}\left(a\right)\right)\right)} \cdot \sqrt[3]{\mathsf{expm1}\left(\sin \left(\mathsf{expm1}\left(a\right)\right)\right)}\right) \cdot \sqrt[3]{\mathsf{expm1}\left(\sin \left(\mathsf{expm1}\left(a\right)\right)\right)}}{\tan^{-1} a}\right) \bmod a\right)}\right) \cdot \sqrt[3]{\left(\sqrt[3]{\left(\left(\tan^{-1}_* \frac{\left(\sqrt[3]{\mathsf{expm1}\left(\sin \left(\mathsf{expm1}\left(a\right)\right)\right)} \cdot \sqrt[3]{\mathsf{expm1}\left(\sin \left(\mathsf{expm1}\left(a\right)\right)\right)}\right) \cdot \sqrt[3]{\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{\left(\sqrt[3]{\mathsf{expm1}\left(\sin \left(\mathsf{expm1}\left(a\right)\right)\right)} \cdot \sqrt[3]{\mathsf{expm1}\left(\sin \left(\mathsf{expm1}\left(a\right)\right)\right)}\right) \cdot \sqrt[3]{\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(\tan^{-1}_* \frac{\left(\sqrt[3]{\mathsf{expm1}\left(\sin \left(\mathsf{expm1}\left(a\right)\right)\right)} \cdot \sqrt[3]{\mathsf{expm1}\left(\sin \left(\mathsf{expm1}\left(a\right)\right)\right)}\right) \cdot \sqrt[3]{\mathsf{expm1}\left(\sin \left(\mathsf{expm1}\left(a\right)\right)\right)}}{\tan^{-1} a}\right) \bmod a\right)}}\right|double f(double a) {
double r19970 = a;
double r19971 = expm1(r19970);
double r19972 = sin(r19971);
double r19973 = expm1(r19972);
double r19974 = atan(r19970);
double r19975 = atan2(r19973, r19974);
double r19976 = fmod(r19975, r19970);
double r19977 = fabs(r19976);
return r19977;
}
double f(double a) {
double r19978 = a;
double r19979 = expm1(r19978);
double r19980 = sin(r19979);
double r19981 = expm1(r19980);
double r19982 = cbrt(r19981);
double r19983 = r19982 * r19982;
double r19984 = r19983 * r19982;
double r19985 = atan(r19978);
double r19986 = atan2(r19984, r19985);
double r19987 = fmod(r19986, r19978);
double r19988 = cbrt(r19987);
double r19989 = r19988 * r19988;
double r19990 = r19989 * r19988;
double r19991 = cbrt(r19990);
double r19992 = r19989 * r19991;
double r19993 = fabs(r19992);
return r19993;
}



Bits error versus a
Initial program 33.5
rmApplied add-cube-cbrt33.5
rmApplied add-cube-cbrt33.5
rmApplied add-cube-cbrt33.5
Final simplification33.5
herbie shell --seed 2019297
(FPCore (a)
:name "Random Jason Timeout Test 006"
:precision binary64
(fabs (fmod (atan2 (expm1 (sin (expm1 a))) (atan a)) a)))