\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(\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) \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 \left(\left(\sqrt[3]{\sqrt[3]{\mathsf{expm1}\left(\sin \left(\mathsf{expm1}\left(a\right)\right)\right)}} \cdot \sqrt[3]{\sqrt[3]{\mathsf{expm1}\left(\sin \left(\mathsf{expm1}\left(a\right)\right)\right)}}\right) \cdot \sqrt[3]{\sqrt[3]{\mathsf{expm1}\left(\sin \left(\mathsf{expm1}\left(a\right)\right)\right)}}\right)}{\tan^{-1} a}\right) \bmod a\right)}\right|double f(double a) {
double r16016 = a;
double r16017 = expm1(r16016);
double r16018 = sin(r16017);
double r16019 = expm1(r16018);
double r16020 = atan(r16016);
double r16021 = atan2(r16019, r16020);
double r16022 = fmod(r16021, r16016);
double r16023 = fabs(r16022);
return r16023;
}
double f(double a) {
double r16024 = a;
double r16025 = expm1(r16024);
double r16026 = sin(r16025);
double r16027 = expm1(r16026);
double r16028 = cbrt(r16027);
double r16029 = r16028 * r16028;
double r16030 = r16029 * r16028;
double r16031 = atan(r16024);
double r16032 = atan2(r16030, r16031);
double r16033 = fmod(r16032, r16024);
double r16034 = cbrt(r16033);
double r16035 = r16034 * r16034;
double r16036 = r16035 * r16034;
double r16037 = cbrt(r16036);
double r16038 = r16034 * r16037;
double r16039 = cbrt(r16028);
double r16040 = r16039 * r16039;
double r16041 = r16040 * r16039;
double r16042 = r16029 * r16041;
double r16043 = atan2(r16042, r16031);
double r16044 = fmod(r16043, r16024);
double r16045 = cbrt(r16044);
double r16046 = r16038 * r16045;
double r16047 = fabs(r16046);
return r16047;
}



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