\frac{\log \left(\sqrt{re \cdot re + im \cdot im}\right)}{\log 10}\begin{array}{l}
\mathbf{if}\;re \le -1.300812438992646141617859246198844532718 \cdot 10^{100}:\\
\;\;\;\;\left(\left(\log \left(\frac{-1}{re}\right) \cdot -2\right) \cdot \frac{1}{\sqrt{\log 10}}\right) \cdot \frac{\frac{1}{2}}{\sqrt{\log 10}}\\
\mathbf{elif}\;re \le 5.714402214507161350041984173167312711037 \cdot 10^{91}:\\
\;\;\;\;\left(\frac{1}{\sqrt{\log 10}} \cdot \log \left(im \cdot im + re \cdot re\right)\right) \cdot \frac{\frac{1}{2}}{\sqrt{\log 10}}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{1}{2}}{\sqrt{\log 10}} \cdot \left(\left(2 \cdot \log re\right) \cdot \frac{1}{\sqrt{\log 10}}\right)\\
\end{array}double f(double re, double im) {
double r1121695 = re;
double r1121696 = r1121695 * r1121695;
double r1121697 = im;
double r1121698 = r1121697 * r1121697;
double r1121699 = r1121696 + r1121698;
double r1121700 = sqrt(r1121699);
double r1121701 = log(r1121700);
double r1121702 = 10.0;
double r1121703 = log(r1121702);
double r1121704 = r1121701 / r1121703;
return r1121704;
}
double f(double re, double im) {
double r1121705 = re;
double r1121706 = -1.3008124389926461e+100;
bool r1121707 = r1121705 <= r1121706;
double r1121708 = -1.0;
double r1121709 = r1121708 / r1121705;
double r1121710 = log(r1121709);
double r1121711 = -2.0;
double r1121712 = r1121710 * r1121711;
double r1121713 = 1.0;
double r1121714 = 10.0;
double r1121715 = log(r1121714);
double r1121716 = sqrt(r1121715);
double r1121717 = r1121713 / r1121716;
double r1121718 = r1121712 * r1121717;
double r1121719 = 0.5;
double r1121720 = r1121719 / r1121716;
double r1121721 = r1121718 * r1121720;
double r1121722 = 5.714402214507161e+91;
bool r1121723 = r1121705 <= r1121722;
double r1121724 = im;
double r1121725 = r1121724 * r1121724;
double r1121726 = r1121705 * r1121705;
double r1121727 = r1121725 + r1121726;
double r1121728 = log(r1121727);
double r1121729 = r1121717 * r1121728;
double r1121730 = r1121729 * r1121720;
double r1121731 = 2.0;
double r1121732 = log(r1121705);
double r1121733 = r1121731 * r1121732;
double r1121734 = r1121733 * r1121717;
double r1121735 = r1121720 * r1121734;
double r1121736 = r1121723 ? r1121730 : r1121735;
double r1121737 = r1121707 ? r1121721 : r1121736;
return r1121737;
}



Bits error versus re



Bits error versus im
Results
if re < -1.3008124389926461e+100Initial program 51.9
rmApplied add-sqr-sqrt51.9
Applied pow1/251.9
Applied log-pow51.9
Applied times-frac51.9
rmApplied div-inv51.9
Taylor expanded around -inf 9.2
Simplified9.2
if -1.3008124389926461e+100 < re < 5.714402214507161e+91Initial program 21.9
rmApplied add-sqr-sqrt21.9
Applied pow1/221.9
Applied log-pow21.9
Applied times-frac21.9
rmApplied div-inv21.8
if 5.714402214507161e+91 < re Initial program 49.9
rmApplied add-sqr-sqrt49.9
Applied pow1/249.9
Applied log-pow49.9
Applied times-frac49.9
rmApplied div-inv49.9
Taylor expanded around inf 8.9
Simplified8.9
Final simplification17.3
herbie shell --seed 2019171
(FPCore (re im)
:name "math.log10 on complex, real part"
(/ (log (sqrt (+ (* re re) (* im im)))) (log 10.0)))