\frac{\log \left(\sqrt{re \cdot re + im \cdot im}\right)}{\log 10}\frac{1}{\sqrt{\log 10}} \cdot \left(\frac{1}{\sqrt{\log 10}} \cdot \log \left(\mathsf{hypot}\left(re, im\right)\right)\right)double f(double re, double im) {
double r712814 = re;
double r712815 = r712814 * r712814;
double r712816 = im;
double r712817 = r712816 * r712816;
double r712818 = r712815 + r712817;
double r712819 = sqrt(r712818);
double r712820 = log(r712819);
double r712821 = 10.0;
double r712822 = log(r712821);
double r712823 = r712820 / r712822;
return r712823;
}
double f(double re, double im) {
double r712824 = 1.0;
double r712825 = 10.0;
double r712826 = log(r712825);
double r712827 = sqrt(r712826);
double r712828 = r712824 / r712827;
double r712829 = re;
double r712830 = im;
double r712831 = hypot(r712829, r712830);
double r712832 = log(r712831);
double r712833 = r712828 * r712832;
double r712834 = r712828 * r712833;
return r712834;
}



Bits error versus re



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