\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 r627177 = a;
double r627178 = expm1(r627177);
double r627179 = sin(r627178);
double r627180 = expm1(r627179);
double r627181 = atan(r627177);
double r627182 = atan2(r627180, r627181);
double r627183 = fmod(r627182, r627177);
double r627184 = fabs(r627183);
return r627184;
}
double f(double a) {
double r627185 = a;
double r627186 = expm1(r627185);
double r627187 = sin(r627186);
double r627188 = expm1(r627187);
double r627189 = atan(r627185);
double r627190 = atan2(r627188, r627189);
double r627191 = fmod(r627190, r627185);
double r627192 = cbrt(r627191);
double r627193 = r627192 * r627192;
double r627194 = r627192 * r627193;
double r627195 = cbrt(r627194);
double r627196 = r627192 * r627195;
double r627197 = r627195 * r627196;
double r627198 = fabs(r627197);
return r627198;
}



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)))