double code(double re, double im, double base) {
return ((double) (((double) (((double) (((double) log(((double) sqrt(((double) (((double) (re * re)) + ((double) (im * im)))))))) * ((double) log(base)))) + ((double) (((double) atan2(im, re)) * 0.0)))) / ((double) (((double) (((double) log(base)) * ((double) log(base)))) + ((double) (0.0 * 0.0))))));
}
double code(double re, double im, double base) {
double VAR;
if ((re <= -3.2792927917367154e+72)) {
VAR = ((double) (((double) (((double) log(1.0)) - ((double) log(((double) (-1.0 / re)))))) / ((double) log(base))));
} else {
double VAR_1;
if ((re <= 1.3644360600461893e-269)) {
VAR_1 = ((double) (((double) (((double) (((double) log(base)) * ((double) log(((double) sqrt(((double) (((double) (re * re)) + ((double) (im * im)))))))))) + ((double) (((double) atan2(im, re)) * 0.0)))) / ((double) (((double) cbrt(((double) pow(((double) log(base)), 6.0)))) + ((double) (0.0 * 0.0))))));
} else {
double VAR_2;
if ((re <= 3.4160082491612644e-169)) {
VAR_2 = ((double) (((double) (((double) log(1.0)) + ((double) log(im)))) / ((double) log(base))));
} else {
double VAR_3;
if ((re <= 7.175104160147713e+17)) {
VAR_3 = ((double) cbrt(((double) pow(((double) (((double) (((double) (((double) log(base)) * ((double) log(((double) sqrt(((double) (((double) (re * re)) + ((double) (im * im)))))))))) + ((double) (((double) atan2(im, re)) * 0.0)))) / ((double) (((double) (0.0 * 0.0)) + ((double) pow(((double) log(base)), 2.0)))))), 3.0))));
} else {
VAR_3 = ((double) (((double) log(re)) / ((double) log(base))));
}
VAR_2 = VAR_3;
}
VAR_1 = VAR_2;
}
VAR = VAR_1;
}
return VAR;
}



Bits error versus re



Bits error versus im



Bits error versus base
Results
if re < -3.2792927917367154e72Initial program 47.5
Taylor expanded around -inf 64.0
Simplified10.2
if -3.2792927917367154e72 < re < 1.36443606004618934e-269Initial program 23.6
rmApplied add-cbrt-cube23.8
Applied add-cbrt-cube23.9
Applied cbrt-unprod23.8
Simplified23.7
if 1.36443606004618934e-269 < re < 3.41600824916126442e-169Initial program 34.7
Taylor expanded around 0 35.7
Simplified35.7
if 3.41600824916126442e-169 < re < 717510416014771328Initial program 16.7
rmApplied add-cbrt-cube16.9
Applied add-cbrt-cube17.0
Applied cbrt-undiv16.9
Simplified16.9
if 717510416014771328 < re Initial program 42.3
Taylor expanded around inf 12.7
Simplified12.7
Final simplification18.5
herbie shell --seed 2020191
(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))))