\frac{\log \left(\sqrt{re \cdot re + im \cdot im}\right)}{\log 10}\frac{1}{\sqrt{\log 10}} \cdot \log \left({\left(\mathsf{hypot}\left(re, im\right)\right)}^{\left(\frac{1}{\sqrt{\log 10}}\right)}\right)double f(double re, double im) {
double r36886 = re;
double r36887 = r36886 * r36886;
double r36888 = im;
double r36889 = r36888 * r36888;
double r36890 = r36887 + r36889;
double r36891 = sqrt(r36890);
double r36892 = log(r36891);
double r36893 = 10.0;
double r36894 = log(r36893);
double r36895 = r36892 / r36894;
return r36895;
}
double f(double re, double im) {
double r36896 = 1.0;
double r36897 = 10.0;
double r36898 = log(r36897);
double r36899 = sqrt(r36898);
double r36900 = r36896 / r36899;
double r36901 = re;
double r36902 = im;
double r36903 = hypot(r36901, r36902);
double r36904 = pow(r36903, r36900);
double r36905 = log(r36904);
double r36906 = r36900 * r36905;
return r36906;
}



Bits error versus re



Bits error versus im
Results
Initial program 32.0
Simplified0.6
rmApplied add-sqr-sqrt0.6
Applied pow10.6
Applied log-pow0.6
Applied times-frac0.6
rmApplied div-inv0.4
rmApplied add-log-exp0.4
Simplified0.3
Final simplification0.3
herbie shell --seed 2019306 +o rules:numerics
(FPCore (re im)
:name "math.log10 on complex, real part"
:precision binary64
(/ (log (sqrt (+ (* re re) (* im im)))) (log 10)))