
(FPCore (re im base) :precision binary64 (/ (+ (* (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) {
return ((log(sqrt(((re * re) + (im * im)))) * log(base)) + (atan2(im, re) * 0.0)) / ((log(base) * log(base)) + (0.0 * 0.0));
}
real(8) function code(re, im, base)
real(8), intent (in) :: re
real(8), intent (in) :: im
real(8), intent (in) :: base
code = ((log(sqrt(((re * re) + (im * im)))) * log(base)) + (atan2(im, re) * 0.0d0)) / ((log(base) * log(base)) + (0.0d0 * 0.0d0))
end function
public static double code(double re, double im, double base) {
return ((Math.log(Math.sqrt(((re * re) + (im * im)))) * Math.log(base)) + (Math.atan2(im, re) * 0.0)) / ((Math.log(base) * Math.log(base)) + (0.0 * 0.0));
}
def code(re, im, base): return ((math.log(math.sqrt(((re * re) + (im * im)))) * math.log(base)) + (math.atan2(im, re) * 0.0)) / ((math.log(base) * math.log(base)) + (0.0 * 0.0))
function code(re, im, base) return Float64(Float64(Float64(log(sqrt(Float64(Float64(re * re) + Float64(im * im)))) * log(base)) + Float64(atan(im, re) * 0.0)) / Float64(Float64(log(base) * log(base)) + Float64(0.0 * 0.0))) end
function tmp = code(re, im, base) tmp = ((log(sqrt(((re * re) + (im * im)))) * log(base)) + (atan2(im, re) * 0.0)) / ((log(base) * log(base)) + (0.0 * 0.0)); end
code[re_, im_, base_] := N[(N[(N[(N[Log[N[Sqrt[N[(N[(re * re), $MachinePrecision] + N[(im * im), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]], $MachinePrecision] * N[Log[base], $MachinePrecision]), $MachinePrecision] + N[(N[ArcTan[im / re], $MachinePrecision] * 0.0), $MachinePrecision]), $MachinePrecision] / N[(N[(N[Log[base], $MachinePrecision] * N[Log[base], $MachinePrecision]), $MachinePrecision] + N[(0.0 * 0.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{\log \left(\sqrt{re \cdot re + im \cdot im}\right) \cdot \log base + \tan^{-1}_* \frac{im}{re} \cdot 0}{\log base \cdot \log base + 0 \cdot 0}
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 6 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (re im base) :precision binary64 (/ (+ (* (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) {
return ((log(sqrt(((re * re) + (im * im)))) * log(base)) + (atan2(im, re) * 0.0)) / ((log(base) * log(base)) + (0.0 * 0.0));
}
real(8) function code(re, im, base)
real(8), intent (in) :: re
real(8), intent (in) :: im
real(8), intent (in) :: base
code = ((log(sqrt(((re * re) + (im * im)))) * log(base)) + (atan2(im, re) * 0.0d0)) / ((log(base) * log(base)) + (0.0d0 * 0.0d0))
end function
public static double code(double re, double im, double base) {
return ((Math.log(Math.sqrt(((re * re) + (im * im)))) * Math.log(base)) + (Math.atan2(im, re) * 0.0)) / ((Math.log(base) * Math.log(base)) + (0.0 * 0.0));
}
def code(re, im, base): return ((math.log(math.sqrt(((re * re) + (im * im)))) * math.log(base)) + (math.atan2(im, re) * 0.0)) / ((math.log(base) * math.log(base)) + (0.0 * 0.0))
function code(re, im, base) return Float64(Float64(Float64(log(sqrt(Float64(Float64(re * re) + Float64(im * im)))) * log(base)) + Float64(atan(im, re) * 0.0)) / Float64(Float64(log(base) * log(base)) + Float64(0.0 * 0.0))) end
function tmp = code(re, im, base) tmp = ((log(sqrt(((re * re) + (im * im)))) * log(base)) + (atan2(im, re) * 0.0)) / ((log(base) * log(base)) + (0.0 * 0.0)); end
code[re_, im_, base_] := N[(N[(N[(N[Log[N[Sqrt[N[(N[(re * re), $MachinePrecision] + N[(im * im), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]], $MachinePrecision] * N[Log[base], $MachinePrecision]), $MachinePrecision] + N[(N[ArcTan[im / re], $MachinePrecision] * 0.0), $MachinePrecision]), $MachinePrecision] / N[(N[(N[Log[base], $MachinePrecision] * N[Log[base], $MachinePrecision]), $MachinePrecision] + N[(0.0 * 0.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{\log \left(\sqrt{re \cdot re + im \cdot im}\right) \cdot \log base + \tan^{-1}_* \frac{im}{re} \cdot 0}{\log base \cdot \log base + 0 \cdot 0}
\end{array}
(FPCore (re im base) :precision binary64 (/ (log (hypot re im)) (log base)))
double code(double re, double im, double base) {
return log(hypot(re, im)) / log(base);
}
public static double code(double re, double im, double base) {
return Math.log(Math.hypot(re, im)) / Math.log(base);
}
def code(re, im, base): return math.log(math.hypot(re, im)) / math.log(base)
function code(re, im, base) return Float64(log(hypot(re, im)) / log(base)) end
function tmp = code(re, im, base) tmp = log(hypot(re, im)) / log(base); end
code[re_, im_, base_] := N[(N[Log[N[Sqrt[re ^ 2 + im ^ 2], $MachinePrecision]], $MachinePrecision] / N[Log[base], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{\log \left(\mathsf{hypot}\left(re, im\right)\right)}{\log base}
\end{array}
Initial program 46.6%
fma-define46.6%
+-commutative46.6%
mul0-rgt46.6%
mul0-rgt46.6%
fma-define46.6%
Simplified99.5%
(FPCore (re im base) :precision binary64 (* 0.5 (/ (+ (* (/ re im) (/ re im)) (* (log im) 2.0)) (log base))))
double code(double re, double im, double base) {
return 0.5 * ((((re / im) * (re / im)) + (log(im) * 2.0)) / log(base));
}
real(8) function code(re, im, base)
real(8), intent (in) :: re
real(8), intent (in) :: im
real(8), intent (in) :: base
code = 0.5d0 * ((((re / im) * (re / im)) + (log(im) * 2.0d0)) / log(base))
end function
public static double code(double re, double im, double base) {
return 0.5 * ((((re / im) * (re / im)) + (Math.log(im) * 2.0)) / Math.log(base));
}
def code(re, im, base): return 0.5 * ((((re / im) * (re / im)) + (math.log(im) * 2.0)) / math.log(base))
function code(re, im, base) return Float64(0.5 * Float64(Float64(Float64(Float64(re / im) * Float64(re / im)) + Float64(log(im) * 2.0)) / log(base))) end
function tmp = code(re, im, base) tmp = 0.5 * ((((re / im) * (re / im)) + (log(im) * 2.0)) / log(base)); end
code[re_, im_, base_] := N[(0.5 * N[(N[(N[(N[(re / im), $MachinePrecision] * N[(re / im), $MachinePrecision]), $MachinePrecision] + N[(N[Log[im], $MachinePrecision] * 2.0), $MachinePrecision]), $MachinePrecision] / N[Log[base], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
0.5 \cdot \frac{\frac{re}{im} \cdot \frac{re}{im} + \log im \cdot 2}{\log base}
\end{array}
Initial program 46.6%
fma-define46.6%
+-commutative46.6%
mul0-rgt46.6%
mul0-rgt46.6%
fma-define46.6%
Simplified99.5%
expm1-log1p-u74.4%
expm1-undefine74.1%
hypot-define38.6%
pow1/238.6%
log-pow38.6%
Applied egg-rr38.6%
sub-neg38.6%
metadata-eval38.6%
+-commutative38.6%
log1p-undefine38.6%
rem-exp-log46.3%
*-commutative46.3%
associate-/l*46.3%
Simplified46.3%
Taylor expanded in im around inf 26.2%
+-commutative26.2%
unpow226.2%
unpow226.2%
times-frac28.7%
*-commutative28.7%
log-rec28.7%
sub0-neg28.7%
Simplified28.7%
Taylor expanded in base around 0 26.4%
+-commutative26.4%
unpow226.4%
unpow226.4%
times-frac28.8%
*-commutative28.8%
Simplified28.8%
(FPCore (re im base) :precision binary64 (/ (log (/ 1.0 im)) (- (log base))))
double code(double re, double im, double base) {
return log((1.0 / im)) / -log(base);
}
real(8) function code(re, im, base)
real(8), intent (in) :: re
real(8), intent (in) :: im
real(8), intent (in) :: base
code = log((1.0d0 / im)) / -log(base)
end function
public static double code(double re, double im, double base) {
return Math.log((1.0 / im)) / -Math.log(base);
}
def code(re, im, base): return math.log((1.0 / im)) / -math.log(base)
function code(re, im, base) return Float64(log(Float64(1.0 / im)) / Float64(-log(base))) end
function tmp = code(re, im, base) tmp = log((1.0 / im)) / -log(base); end
code[re_, im_, base_] := N[(N[Log[N[(1.0 / im), $MachinePrecision]], $MachinePrecision] / (-N[Log[base], $MachinePrecision])), $MachinePrecision]
\begin{array}{l}
\\
\frac{\log \left(\frac{1}{im}\right)}{-\log base}
\end{array}
Initial program 46.6%
fma-define46.6%
+-commutative46.6%
mul0-rgt46.6%
mul0-rgt46.6%
fma-define46.6%
Simplified99.5%
Taylor expanded in im around inf 30.9%
Final simplification30.9%
(FPCore (re im base) :precision binary64 (/ (log im) (log base)))
double code(double re, double im, double base) {
return log(im) / log(base);
}
real(8) function code(re, im, base)
real(8), intent (in) :: re
real(8), intent (in) :: im
real(8), intent (in) :: base
code = log(im) / log(base)
end function
public static double code(double re, double im, double base) {
return Math.log(im) / Math.log(base);
}
def code(re, im, base): return math.log(im) / math.log(base)
function code(re, im, base) return Float64(log(im) / log(base)) end
function tmp = code(re, im, base) tmp = log(im) / log(base); end
code[re_, im_, base_] := N[(N[Log[im], $MachinePrecision] / N[Log[base], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{\log im}{\log base}
\end{array}
Initial program 46.6%
fma-define46.6%
+-commutative46.6%
mul0-rgt46.6%
mul0-rgt46.6%
fma-define46.6%
Simplified99.5%
Taylor expanded in re around 0 30.9%
(FPCore (re im base) :precision binary64 (* (* (/ re im) (/ re im)) (/ 0.5 (log base))))
double code(double re, double im, double base) {
return ((re / im) * (re / im)) * (0.5 / log(base));
}
real(8) function code(re, im, base)
real(8), intent (in) :: re
real(8), intent (in) :: im
real(8), intent (in) :: base
code = ((re / im) * (re / im)) * (0.5d0 / log(base))
end function
public static double code(double re, double im, double base) {
return ((re / im) * (re / im)) * (0.5 / Math.log(base));
}
def code(re, im, base): return ((re / im) * (re / im)) * (0.5 / math.log(base))
function code(re, im, base) return Float64(Float64(Float64(re / im) * Float64(re / im)) * Float64(0.5 / log(base))) end
function tmp = code(re, im, base) tmp = ((re / im) * (re / im)) * (0.5 / log(base)); end
code[re_, im_, base_] := N[(N[(N[(re / im), $MachinePrecision] * N[(re / im), $MachinePrecision]), $MachinePrecision] * N[(0.5 / N[Log[base], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\left(\frac{re}{im} \cdot \frac{re}{im}\right) \cdot \frac{0.5}{\log base}
\end{array}
Initial program 46.6%
fma-define46.6%
+-commutative46.6%
mul0-rgt46.6%
mul0-rgt46.6%
fma-define46.6%
Simplified99.5%
expm1-log1p-u74.4%
expm1-undefine74.1%
hypot-define38.6%
pow1/238.6%
log-pow38.6%
Applied egg-rr38.6%
sub-neg38.6%
metadata-eval38.6%
+-commutative38.6%
log1p-undefine38.6%
rem-exp-log46.3%
*-commutative46.3%
associate-/l*46.3%
Simplified46.3%
Taylor expanded in im around inf 26.2%
+-commutative26.2%
unpow226.2%
unpow226.2%
times-frac28.7%
*-commutative28.7%
log-rec28.7%
sub0-neg28.7%
Simplified28.7%
Taylor expanded in re around inf 2.8%
associate-*r/2.8%
unpow22.8%
*-commutative2.8%
unpow22.8%
Simplified2.8%
times-frac2.8%
frac-times3.3%
Applied egg-rr3.3%
Final simplification3.3%
(FPCore (re im base) :precision binary64 (* re (/ (/ re im) (* im (* (log base) 2.0)))))
double code(double re, double im, double base) {
return re * ((re / im) / (im * (log(base) * 2.0)));
}
real(8) function code(re, im, base)
real(8), intent (in) :: re
real(8), intent (in) :: im
real(8), intent (in) :: base
code = re * ((re / im) / (im * (log(base) * 2.0d0)))
end function
public static double code(double re, double im, double base) {
return re * ((re / im) / (im * (Math.log(base) * 2.0)));
}
def code(re, im, base): return re * ((re / im) / (im * (math.log(base) * 2.0)))
function code(re, im, base) return Float64(re * Float64(Float64(re / im) / Float64(im * Float64(log(base) * 2.0)))) end
function tmp = code(re, im, base) tmp = re * ((re / im) / (im * (log(base) * 2.0))); end
code[re_, im_, base_] := N[(re * N[(N[(re / im), $MachinePrecision] / N[(im * N[(N[Log[base], $MachinePrecision] * 2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
re \cdot \frac{\frac{re}{im}}{im \cdot \left(\log base \cdot 2\right)}
\end{array}
Initial program 46.6%
fma-define46.6%
+-commutative46.6%
mul0-rgt46.6%
mul0-rgt46.6%
fma-define46.6%
Simplified99.5%
expm1-log1p-u74.4%
expm1-undefine74.1%
hypot-define38.6%
pow1/238.6%
log-pow38.6%
Applied egg-rr38.6%
sub-neg38.6%
metadata-eval38.6%
+-commutative38.6%
log1p-undefine38.6%
rem-exp-log46.3%
*-commutative46.3%
associate-/l*46.3%
Simplified46.3%
Taylor expanded in im around inf 26.2%
+-commutative26.2%
unpow226.2%
unpow226.2%
times-frac28.7%
*-commutative28.7%
log-rec28.7%
sub0-neg28.7%
Simplified28.7%
Taylor expanded in re around inf 2.8%
associate-*r/2.8%
unpow22.8%
*-commutative2.8%
unpow22.8%
Simplified2.8%
times-frac2.8%
clear-num2.8%
frac-times2.8%
*-un-lft-identity2.8%
div-inv2.8%
metadata-eval2.8%
Applied egg-rr2.8%
unpow22.8%
times-frac3.0%
unpow23.0%
associate-/r*3.2%
times-frac3.3%
associate-/l*3.2%
*-commutative3.2%
*-commutative3.2%
Simplified3.2%
Final simplification3.2%
herbie shell --seed 2024107
(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))))