\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|\sqrt[3]{\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)} \cdot \left(\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)} \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 \left(\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)} \cdot \sqrt[3]{\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)} \cdot \left(\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)} \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)}\right)\right|double f(double a) {
double r894015 = a;
double r894016 = expm1(r894015);
double r894017 = sin(r894016);
double r894018 = expm1(r894017);
double r894019 = atan(r894015);
double r894020 = atan2(r894018, r894019);
double r894021 = fmod(r894020, r894015);
double r894022 = fabs(r894021);
return r894022;
}
double f(double a) {
double r894023 = a;
double r894024 = expm1(r894023);
double r894025 = sin(r894024);
double r894026 = expm1(r894025);
double r894027 = atan(r894023);
double r894028 = atan2(r894026, r894027);
double r894029 = fmod(r894028, r894023);
double r894030 = cbrt(r894029);
double r894031 = r894030 * r894030;
double r894032 = r894030 * r894031;
double r894033 = cbrt(r894032);
double r894034 = r894030 * r894033;
double r894035 = r894033 * r894034;
double r894036 = fabs(r894035);
return r894036;
}



Bits error versus a
Initial program 33.8
rmApplied add-cube-cbrt33.8
rmApplied add-cube-cbrt33.8
rmApplied add-cube-cbrt33.8
Final simplification33.8
herbie shell --seed 2019171 +o rules:numerics
(FPCore (a)
:name "Random Jason Timeout Test 006"
(fabs (fmod (atan2 (expm1 (sin (expm1 a))) (atan a)) a)))