\frac{\log \left(\sqrt{re \cdot re + im \cdot im}\right)}{\log 10}\left(\left(\sqrt{\frac{1}{\sqrt{\log 10}}} \cdot \left(\frac{1}{\sqrt{\log 10}} \cdot \log \left(\mathsf{hypot}\left(re, im\right)\right)\right)\right) \cdot \sqrt{\sqrt{\frac{1}{\sqrt{\log 10}}}}\right) \cdot \sqrt{\sqrt{\frac{1}{\sqrt{\log 10}}}}double f(double re, double im) {
double r757042 = re;
double r757043 = r757042 * r757042;
double r757044 = im;
double r757045 = r757044 * r757044;
double r757046 = r757043 + r757045;
double r757047 = sqrt(r757046);
double r757048 = log(r757047);
double r757049 = 10.0;
double r757050 = log(r757049);
double r757051 = r757048 / r757050;
return r757051;
}
double f(double re, double im) {
double r757052 = 1.0;
double r757053 = 10.0;
double r757054 = log(r757053);
double r757055 = sqrt(r757054);
double r757056 = r757052 / r757055;
double r757057 = sqrt(r757056);
double r757058 = re;
double r757059 = im;
double r757060 = hypot(r757058, r757059);
double r757061 = log(r757060);
double r757062 = r757056 * r757061;
double r757063 = r757057 * r757062;
double r757064 = sqrt(r757057);
double r757065 = r757063 * r757064;
double r757066 = r757065 * r757064;
return r757066;
}



Bits error versus re



Bits error versus im
Results
Initial program 31.1
Simplified0.6
rmApplied add-sqr-sqrt0.6
Applied *-un-lft-identity0.6
Applied times-frac0.5
rmApplied div-inv0.4
rmApplied add-sqr-sqrt0.4
Applied associate-*l*0.5
rmApplied add-sqr-sqrt0.5
Applied sqrt-prod0.6
Applied associate-*l*0.5
Final simplification0.5
herbie shell --seed 2019144 +o rules:numerics
(FPCore (re im)
:name "math.log10 on complex, real part"
(/ (log (sqrt (+ (* re re) (* im im)))) (log 10)))