\frac{\log \left(\sqrt{re \cdot re + im \cdot im}\right)}{\log 10}\begin{array}{l}
\mathbf{if}\;re \le -1.3496059084460043 \cdot 10^{+95}:\\
\;\;\;\;\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 5.956366146560391 \cdot 10^{+103}:\\
\;\;\;\;\left(\frac{\sqrt[3]{\frac{1}{2}}}{\sqrt{\sqrt[3]{\log 10}}} \cdot \frac{\log \left(re \cdot re + im \cdot im\right)}{\sqrt{\log 10}}\right) \cdot \frac{\sqrt[3]{\frac{1}{2}} \cdot \sqrt[3]{\frac{1}{2}}}{\sqrt{\sqrt[3]{\log 10} \cdot \sqrt[3]{\log 10}}}\\
\mathbf{else}:\\
\;\;\;\;\left(\sqrt{\frac{1}{\log 10}} \cdot \left(\log re \cdot 2\right)\right) \cdot \frac{\frac{1}{2}}{\sqrt{\log 10}}\\
\end{array}double f(double re, double im) {
double r993266 = re;
double r993267 = r993266 * r993266;
double r993268 = im;
double r993269 = r993268 * r993268;
double r993270 = r993267 + r993269;
double r993271 = sqrt(r993270);
double r993272 = log(r993271);
double r993273 = 10.0;
double r993274 = log(r993273);
double r993275 = r993272 / r993274;
return r993275;
}
double f(double re, double im) {
double r993276 = re;
double r993277 = -1.3496059084460043e+95;
bool r993278 = r993276 <= r993277;
double r993279 = 0.5;
double r993280 = 10.0;
double r993281 = log(r993280);
double r993282 = sqrt(r993281);
double r993283 = r993279 / r993282;
double r993284 = 1.0;
double r993285 = r993284 / r993281;
double r993286 = sqrt(r993285);
double r993287 = -1.0;
double r993288 = r993287 / r993276;
double r993289 = log(r993288);
double r993290 = r993286 * r993289;
double r993291 = -2.0;
double r993292 = r993290 * r993291;
double r993293 = r993283 * r993292;
double r993294 = 5.956366146560391e+103;
bool r993295 = r993276 <= r993294;
double r993296 = cbrt(r993279);
double r993297 = cbrt(r993281);
double r993298 = sqrt(r993297);
double r993299 = r993296 / r993298;
double r993300 = r993276 * r993276;
double r993301 = im;
double r993302 = r993301 * r993301;
double r993303 = r993300 + r993302;
double r993304 = log(r993303);
double r993305 = r993304 / r993282;
double r993306 = r993299 * r993305;
double r993307 = r993296 * r993296;
double r993308 = r993297 * r993297;
double r993309 = sqrt(r993308);
double r993310 = r993307 / r993309;
double r993311 = r993306 * r993310;
double r993312 = log(r993276);
double r993313 = 2.0;
double r993314 = r993312 * r993313;
double r993315 = r993286 * r993314;
double r993316 = r993315 * r993283;
double r993317 = r993295 ? r993311 : r993316;
double r993318 = r993278 ? r993293 : r993317;
return r993318;
}



Bits error versus re



Bits error versus im
Results
if re < -1.3496059084460043e+95Initial program 48.9
rmApplied add-sqr-sqrt48.9
Applied pow1/248.9
Applied log-pow48.9
Applied times-frac48.9
Taylor expanded around -inf 8.4
if -1.3496059084460043e+95 < re < 5.956366146560391e+103Initial program 20.8
rmApplied add-sqr-sqrt20.8
Applied pow1/220.8
Applied log-pow20.8
Applied times-frac20.8
rmApplied add-cube-cbrt21.3
Applied sqrt-prod21.3
Applied add-cube-cbrt20.8
Applied times-frac20.8
Applied associate-*l*20.7
if 5.956366146560391e+103 < re Initial program 51.8
rmApplied add-sqr-sqrt51.8
Applied pow1/251.8
Applied log-pow51.8
Applied times-frac51.8
Taylor expanded around inf 8.5
Simplified8.5
Final simplification16.6
herbie shell --seed 2019152
(FPCore (re im)
:name "math.log10 on complex, real part"
(/ (log (sqrt (+ (* re re) (* im im)))) (log 10)))