
(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 3 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 56.0%
mul0-rgt56.0%
+-rgt-identity56.0%
metadata-eval56.0%
+-rgt-identity56.0%
times-frac56.1%
*-inverses56.1%
*-rgt-identity56.1%
hypot-def99.4%
Simplified99.4%
Final simplification99.4%
(FPCore (re im base) :precision binary64 (if (or (<= re -1.05e-36) (and (not (<= re -3.9e-84)) (<= re -3.9e-178))) (/ (log (- re)) (log base)) (/ (log im) (log base))))
double code(double re, double im, double base) {
double tmp;
if ((re <= -1.05e-36) || (!(re <= -3.9e-84) && (re <= -3.9e-178))) {
tmp = log(-re) / log(base);
} else {
tmp = log(im) / log(base);
}
return tmp;
}
real(8) function code(re, im, base)
real(8), intent (in) :: re
real(8), intent (in) :: im
real(8), intent (in) :: base
real(8) :: tmp
if ((re <= (-1.05d-36)) .or. (.not. (re <= (-3.9d-84))) .and. (re <= (-3.9d-178))) then
tmp = log(-re) / log(base)
else
tmp = log(im) / log(base)
end if
code = tmp
end function
public static double code(double re, double im, double base) {
double tmp;
if ((re <= -1.05e-36) || (!(re <= -3.9e-84) && (re <= -3.9e-178))) {
tmp = Math.log(-re) / Math.log(base);
} else {
tmp = Math.log(im) / Math.log(base);
}
return tmp;
}
def code(re, im, base): tmp = 0 if (re <= -1.05e-36) or (not (re <= -3.9e-84) and (re <= -3.9e-178)): tmp = math.log(-re) / math.log(base) else: tmp = math.log(im) / math.log(base) return tmp
function code(re, im, base) tmp = 0.0 if ((re <= -1.05e-36) || (!(re <= -3.9e-84) && (re <= -3.9e-178))) tmp = Float64(log(Float64(-re)) / log(base)); else tmp = Float64(log(im) / log(base)); end return tmp end
function tmp_2 = code(re, im, base) tmp = 0.0; if ((re <= -1.05e-36) || (~((re <= -3.9e-84)) && (re <= -3.9e-178))) tmp = log(-re) / log(base); else tmp = log(im) / log(base); end tmp_2 = tmp; end
code[re_, im_, base_] := If[Or[LessEqual[re, -1.05e-36], And[N[Not[LessEqual[re, -3.9e-84]], $MachinePrecision], LessEqual[re, -3.9e-178]]], N[(N[Log[(-re)], $MachinePrecision] / N[Log[base], $MachinePrecision]), $MachinePrecision], N[(N[Log[im], $MachinePrecision] / N[Log[base], $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;re \leq -1.05 \cdot 10^{-36} \lor \neg \left(re \leq -3.9 \cdot 10^{-84}\right) \land re \leq -3.9 \cdot 10^{-178}:\\
\;\;\;\;\frac{\log \left(-re\right)}{\log base}\\
\mathbf{else}:\\
\;\;\;\;\frac{\log im}{\log base}\\
\end{array}
\end{array}
if re < -1.04999999999999995e-36 or -3.90000000000000023e-84 < re < -3.90000000000000025e-178Initial program 51.0%
mul0-rgt51.0%
+-rgt-identity51.0%
metadata-eval51.0%
+-rgt-identity51.0%
times-frac51.0%
*-inverses51.0%
*-rgt-identity51.0%
hypot-def99.5%
Simplified99.5%
pow199.5%
metadata-eval99.5%
pow-div99.5%
pow299.5%
pow199.5%
associate-/r/99.2%
pow299.2%
Applied egg-rr99.2%
Taylor expanded in re around -inf 69.3%
associate-*r/69.3%
mul-1-neg69.3%
Simplified69.3%
expm1-log1p-u69.3%
expm1-udef62.6%
div-inv62.6%
neg-log62.6%
frac-2neg62.6%
metadata-eval62.6%
remove-double-div62.6%
pow-flip62.6%
metadata-eval62.6%
Applied egg-rr62.6%
expm1-def69.3%
expm1-log1p69.4%
*-commutative69.4%
Simplified69.4%
Taylor expanded in base around 0 69.5%
if -1.04999999999999995e-36 < re < -3.90000000000000023e-84 or -3.90000000000000025e-178 < re Initial program 57.8%
mul0-rgt57.8%
+-rgt-identity57.8%
metadata-eval57.8%
+-rgt-identity57.8%
times-frac57.9%
*-inverses57.9%
*-rgt-identity57.9%
hypot-def99.3%
Simplified99.3%
Taylor expanded in re around 0 38.0%
Final simplification46.3%
(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 56.0%
mul0-rgt56.0%
+-rgt-identity56.0%
metadata-eval56.0%
+-rgt-identity56.0%
times-frac56.1%
*-inverses56.1%
*-rgt-identity56.1%
hypot-def99.4%
Simplified99.4%
Taylor expanded in re around 0 32.0%
Final simplification32.0%
herbie shell --seed 2023178
(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))))