\frac{\log \left(\sqrt{re \cdot re + im \cdot im}\right) \cdot \log base + \tan^{-1}_* \frac{im}{re} \cdot 0.0}{\log base \cdot \log base + 0.0 \cdot 0.0}\begin{array}{l}
\mathbf{if}\;re \le -1.1225739876249368 \cdot 10^{21}:\\
\;\;\;\;\frac{\frac{\log \left(-1 \cdot re\right) \cdot \log base + \tan^{-1}_* \frac{im}{re} \cdot 0.0}{\sqrt{\log base \cdot \log base + 0.0 \cdot 0.0}}}{\sqrt{\log base \cdot \log base + 0.0 \cdot 0.0}}\\
\mathbf{elif}\;re \le -4.27711132562022042 \cdot 10^{-185}:\\
\;\;\;\;\frac{\frac{\log \left(\left|\sqrt[3]{re \cdot re + im \cdot im}\right| \cdot \sqrt{\sqrt[3]{re \cdot re + im \cdot im}}\right) \cdot \log base + \tan^{-1}_* \frac{im}{re} \cdot 0.0}{\sqrt{\log base \cdot \log base + 0.0 \cdot 0.0}}}{\sqrt{\log base \cdot \log base + 0.0 \cdot 0.0}}\\
\mathbf{elif}\;re \le 3.510677278861231 \cdot 10^{-284}:\\
\;\;\;\;\frac{\log im}{\log base}\\
\mathbf{elif}\;re \le 5.9451169291533458 \cdot 10^{55}:\\
\;\;\;\;\frac{\frac{\log \left(\left|\sqrt[3]{re \cdot re + im \cdot im}\right| \cdot \sqrt{\sqrt[3]{re \cdot re + im \cdot im}}\right) \cdot \log base + \tan^{-1}_* \frac{im}{re} \cdot 0.0}{\sqrt{\log base \cdot \log base + 0.0 \cdot 0.0}}}{\sqrt{\log base \cdot \log base + 0.0 \cdot 0.0}}\\
\mathbf{else}:\\
\;\;\;\;\frac{\log re \cdot \log base + \tan^{-1}_* \frac{im}{re} \cdot 0.0}{\log base \cdot \log base + 0.0 \cdot 0.0}\\
\end{array}double code(double re, double im, double base) {
return (((log(sqrt(((re * re) + (im * im)))) * log(base)) + (atan2(im, re) * 0.0)) / ((log(base) * log(base)) + (0.0 * 0.0)));
}
double code(double re, double im, double base) {
double temp;
if ((re <= -1.1225739876249368e+21)) {
temp = ((((log((-1.0 * re)) * log(base)) + (atan2(im, re) * 0.0)) / sqrt(((log(base) * log(base)) + (0.0 * 0.0)))) / sqrt(((log(base) * log(base)) + (0.0 * 0.0))));
} else {
double temp_1;
if ((re <= -4.2771113256202204e-185)) {
temp_1 = ((((log((fabs(cbrt(((re * re) + (im * im)))) * sqrt(cbrt(((re * re) + (im * im)))))) * log(base)) + (atan2(im, re) * 0.0)) / sqrt(((log(base) * log(base)) + (0.0 * 0.0)))) / sqrt(((log(base) * log(base)) + (0.0 * 0.0))));
} else {
double temp_2;
if ((re <= 3.510677278861231e-284)) {
temp_2 = (log(im) / log(base));
} else {
double temp_3;
if ((re <= 5.945116929153346e+55)) {
temp_3 = ((((log((fabs(cbrt(((re * re) + (im * im)))) * sqrt(cbrt(((re * re) + (im * im)))))) * log(base)) + (atan2(im, re) * 0.0)) / sqrt(((log(base) * log(base)) + (0.0 * 0.0)))) / sqrt(((log(base) * log(base)) + (0.0 * 0.0))));
} else {
temp_3 = (((log(re) * log(base)) + (atan2(im, re) * 0.0)) / ((log(base) * log(base)) + (0.0 * 0.0)));
}
temp_2 = temp_3;
}
temp_1 = temp_2;
}
temp = temp_1;
}
return temp;
}



Bits error versus re



Bits error versus im



Bits error versus base
Results
if re < -1.1225739876249368e+21Initial program 43.2
rmApplied add-sqr-sqrt43.2
Applied associate-/r*43.2
Taylor expanded around -inf 12.5
if -1.1225739876249368e+21 < re < -4.2771113256202204e-185 or 3.510677278861231e-284 < re < 5.945116929153346e+55Initial program 20.3
rmApplied add-sqr-sqrt20.3
Applied associate-/r*20.2
rmApplied add-cube-cbrt20.2
Applied sqrt-prod20.2
Simplified20.2
if -4.2771113256202204e-185 < re < 3.510677278861231e-284Initial program 32.3
Taylor expanded around 0 35.1
if 5.945116929153346e+55 < re Initial program 45.5
Taylor expanded around inf 10.8
Final simplification18.3
herbie shell --seed 2020057
(FPCore (re im base)
:name "math.log/2 on complex, real part"
:precision binary64
(/ (+ (* (log (sqrt (+ (* re re) (* im im)))) (log base)) (* (atan2 im re) 0.0)) (+ (* (log base) (log base)) (* 0.0 0.0))))