double f(double re, double im) {
double r1435691 = re;
double r1435692 = r1435691 * r1435691;
double r1435693 = im;
double r1435694 = r1435693 * r1435693;
double r1435695 = r1435692 + r1435694;
double r1435696 = sqrt(r1435695);
double r1435697 = log(r1435696);
double r1435698 = 10.0;
double r1435699 = log(r1435698);
double r1435700 = r1435697 / r1435699;
return r1435700;
}
double f(double re, double im) {
double r1435701 = re;
double r1435702 = im;
double r1435703 = hypot(r1435701, r1435702);
double r1435704 = cbrt(r1435703);
double r1435705 = log(r1435704);
double r1435706 = 1.0;
double r1435707 = 10.0;
double r1435708 = log(r1435707);
double r1435709 = sqrt(r1435708);
double r1435710 = r1435706 / r1435709;
double r1435711 = 3.0;
double r1435712 = r1435710 * r1435711;
double r1435713 = r1435705 * r1435712;
double r1435714 = r1435713 * r1435710;
return r1435714;
}
\frac{\log \left(\sqrt{re \cdot re + im \cdot im}\right)}{\log 10}\left(\log \left(\sqrt[3]{\sqrt{re^2 + im^2}^*}\right) \cdot \left(\frac{1}{\sqrt{\log 10}} \cdot 3\right)\right) \cdot \frac{1}{\sqrt{\log 10}}


Bits error versus re



Bits error versus im
Initial program 30.5
Simplified0.6
rmApplied add-cube-cbrt0.6
rmApplied pow10.6
Applied log-pow0.6
Applied pow10.6
Applied pow10.6
Applied pow-sqr0.6
Applied pow-plus0.6
Applied log-pow0.6
Applied times-frac0.6
Simplified0.6
rmApplied add-sqr-sqrt0.6
Applied *-un-lft-identity0.6
Applied times-frac0.6
Applied associate-*r*0.5
rmApplied div-inv0.4
Applied associate-*r*0.4
Final simplification0.4
herbie shell --seed 2019101 +o rules:numerics
(FPCore (re im)
:name "math.log10 on complex, real part"
(/ (log (sqrt (+ (* re re) (* im im)))) (log 10)))