\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(\sqrt{re^2 + im^2}^*\right)\right)double f(double re, double im) {
double r825429 = re;
double r825430 = r825429 * r825429;
double r825431 = im;
double r825432 = r825431 * r825431;
double r825433 = r825430 + r825432;
double r825434 = sqrt(r825433);
double r825435 = log(r825434);
double r825436 = 10.0;
double r825437 = log(r825436);
double r825438 = r825435 / r825437;
return r825438;
}
double f(double re, double im) {
double r825439 = 1.0;
double r825440 = 10.0;
double r825441 = log(r825440);
double r825442 = sqrt(r825441);
double r825443 = r825439 / r825442;
double r825444 = re;
double r825445 = im;
double r825446 = hypot(r825444, r825445);
double r825447 = log(r825446);
double r825448 = r825443 * r825447;
double r825449 = r825443 * r825448;
return r825449;
}



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.6
rmApplied div-inv0.4
Applied associate-*r*0.4
Final simplification0.4
herbie shell --seed 2019119 +o rules:numerics
(FPCore (re im)
:name "math.log10 on complex, real part"
(/ (log (sqrt (+ (* re re) (* im im)))) (log 10)))