\frac{\log \left(\sqrt{re \cdot re + im \cdot im}\right)}{\log 10}\begin{array}{l}
\mathbf{if}\;re \le -4.844506481309793477015811086474695118347 \cdot 10^{101}:\\
\;\;\;\;\frac{\frac{1}{2}}{\sqrt{\log 10}} \cdot \left(\left(\sqrt{\frac{1}{\log 10}} \cdot \log \left(\frac{-1}{re}\right)\right) \cdot -2\right)\\
\mathbf{elif}\;re \le 6.3735182711102215089283798883990759293 \cdot 10^{95}:\\
\;\;\;\;\frac{\sqrt[3]{\frac{1}{2}} \cdot \sqrt[3]{\frac{1}{2}}}{\sqrt{\sqrt{\log 10}}} \cdot \left(\frac{\log \left(re \cdot re + im \cdot im\right)}{\sqrt{\log 10}} \cdot \frac{\sqrt[3]{\frac{1}{2}}}{\sqrt{\sqrt{\log 10}}}\right)\\
\mathbf{else}:\\
\;\;\;\;\frac{\sqrt[3]{\frac{1}{2}} \cdot \sqrt[3]{\frac{1}{2}}}{\sqrt{\sqrt[3]{\log 10} \cdot \sqrt[3]{\log 10}}} \cdot \left(\frac{\sqrt[3]{\frac{1}{2}}}{\sqrt{\sqrt[3]{\log 10}}} \cdot \frac{2 \cdot \log re}{\sqrt{\log 10}}\right)\\
\end{array}double f(double re, double im) {
double r999034 = re;
double r999035 = r999034 * r999034;
double r999036 = im;
double r999037 = r999036 * r999036;
double r999038 = r999035 + r999037;
double r999039 = sqrt(r999038);
double r999040 = log(r999039);
double r999041 = 10.0;
double r999042 = log(r999041);
double r999043 = r999040 / r999042;
return r999043;
}
double f(double re, double im) {
double r999044 = re;
double r999045 = -4.8445064813097935e+101;
bool r999046 = r999044 <= r999045;
double r999047 = 0.5;
double r999048 = 10.0;
double r999049 = log(r999048);
double r999050 = sqrt(r999049);
double r999051 = r999047 / r999050;
double r999052 = 1.0;
double r999053 = r999052 / r999049;
double r999054 = sqrt(r999053);
double r999055 = -1.0;
double r999056 = r999055 / r999044;
double r999057 = log(r999056);
double r999058 = r999054 * r999057;
double r999059 = -2.0;
double r999060 = r999058 * r999059;
double r999061 = r999051 * r999060;
double r999062 = 6.373518271110222e+95;
bool r999063 = r999044 <= r999062;
double r999064 = cbrt(r999047);
double r999065 = r999064 * r999064;
double r999066 = sqrt(r999050);
double r999067 = r999065 / r999066;
double r999068 = r999044 * r999044;
double r999069 = im;
double r999070 = r999069 * r999069;
double r999071 = r999068 + r999070;
double r999072 = log(r999071);
double r999073 = r999072 / r999050;
double r999074 = r999064 / r999066;
double r999075 = r999073 * r999074;
double r999076 = r999067 * r999075;
double r999077 = cbrt(r999049);
double r999078 = r999077 * r999077;
double r999079 = sqrt(r999078);
double r999080 = r999065 / r999079;
double r999081 = sqrt(r999077);
double r999082 = r999064 / r999081;
double r999083 = 2.0;
double r999084 = log(r999044);
double r999085 = r999083 * r999084;
double r999086 = r999085 / r999050;
double r999087 = r999082 * r999086;
double r999088 = r999080 * r999087;
double r999089 = r999063 ? r999076 : r999088;
double r999090 = r999046 ? r999061 : r999089;
return r999090;
}



Bits error versus re



Bits error versus im
Results
if re < -4.8445064813097935e+101Initial program 52.5
rmApplied add-sqr-sqrt52.5
Applied pow1/252.5
Applied log-pow52.5
Applied times-frac52.4
Taylor expanded around -inf 8.0
if -4.8445064813097935e+101 < re < 6.373518271110222e+95Initial program 22.3
rmApplied add-sqr-sqrt22.3
Applied pow1/222.3
Applied log-pow22.3
Applied times-frac22.3
rmApplied add-sqr-sqrt22.3
Applied sqrt-prod22.8
Applied add-cube-cbrt22.3
Applied times-frac22.3
Applied associate-*l*22.2
if 6.373518271110222e+95 < re Initial program 51.2
rmApplied add-sqr-sqrt51.2
Applied pow1/251.2
Applied log-pow51.2
Applied times-frac51.2
rmApplied add-cube-cbrt51.4
Applied sqrt-prod51.4
Applied add-cube-cbrt51.2
Applied times-frac51.2
Applied associate-*l*51.2
Taylor expanded around inf 9.0
Simplified9.0
Final simplification17.5
herbie shell --seed 2019170
(FPCore (re im)
:name "math.log10 on complex, real part"
(/ (log (sqrt (+ (* re re) (* im im)))) (log 10.0)))