
(FPCore (x) :precision binary64 (- (log (- (/ 1.0 x) 1.0))))
double code(double x) {
return -log(((1.0 / x) - 1.0));
}
real(8) function code(x)
real(8), intent (in) :: x
code = -log(((1.0d0 / x) - 1.0d0))
end function
public static double code(double x) {
return -Math.log(((1.0 / x) - 1.0));
}
def code(x): return -math.log(((1.0 / x) - 1.0))
function code(x) return Float64(-log(Float64(Float64(1.0 / x) - 1.0))) end
function tmp = code(x) tmp = -log(((1.0 / x) - 1.0)); end
code[x_] := (-N[Log[N[(N[(1.0 / x), $MachinePrecision] - 1.0), $MachinePrecision]], $MachinePrecision])
\begin{array}{l}
\\
-\log \left(\frac{1}{x} - 1\right)
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 8 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (x) :precision binary64 (- (log (- (/ 1.0 x) 1.0))))
double code(double x) {
return -log(((1.0 / x) - 1.0));
}
real(8) function code(x)
real(8), intent (in) :: x
code = -log(((1.0d0 / x) - 1.0d0))
end function
public static double code(double x) {
return -Math.log(((1.0 / x) - 1.0));
}
def code(x): return -math.log(((1.0 / x) - 1.0))
function code(x) return Float64(-log(Float64(Float64(1.0 / x) - 1.0))) end
function tmp = code(x) tmp = -log(((1.0 / x) - 1.0)); end
code[x_] := (-N[Log[N[(N[(1.0 / x), $MachinePrecision] - 1.0), $MachinePrecision]], $MachinePrecision])
\begin{array}{l}
\\
-\log \left(\frac{1}{x} - 1\right)
\end{array}
(FPCore (x) :precision binary64 (- (log (+ (/ 1.0 x) -1.0))))
double code(double x) {
return -log(((1.0 / x) + -1.0));
}
real(8) function code(x)
real(8), intent (in) :: x
code = -log(((1.0d0 / x) + (-1.0d0)))
end function
public static double code(double x) {
return -Math.log(((1.0 / x) + -1.0));
}
def code(x): return -math.log(((1.0 / x) + -1.0))
function code(x) return Float64(-log(Float64(Float64(1.0 / x) + -1.0))) end
function tmp = code(x) tmp = -log(((1.0 / x) + -1.0)); end
code[x_] := (-N[Log[N[(N[(1.0 / x), $MachinePrecision] + -1.0), $MachinePrecision]], $MachinePrecision])
\begin{array}{l}
\\
-\log \left(\frac{1}{x} + -1\right)
\end{array}
Initial program 100.0%
Final simplification100.0%
(FPCore (x) :precision binary64 (fma (* x x) 0.5 (+ x (log x))))
double code(double x) {
return fma((x * x), 0.5, (x + log(x)));
}
function code(x) return fma(Float64(x * x), 0.5, Float64(x + log(x))) end
code[x_] := N[(N[(x * x), $MachinePrecision] * 0.5 + N[(x + N[Log[x], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\mathsf{fma}\left(x \cdot x, 0.5, x + \log x\right)
\end{array}
Initial program 100.0%
Taylor expanded in x around 0
cancel-sign-sub-invN/A
metadata-evalN/A
*-lft-identityN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower-log.f6499.7
Applied rewrites99.7%
distribute-lft-inN/A
*-rgt-identityN/A
lift-log.f64N/A
associate-+l+N/A
associate-*r*N/A
lift-+.f64N/A
lower-fma.f64N/A
lower-*.f6499.7
Applied rewrites99.7%
(FPCore (x) :precision binary64 (fma x (fma x 0.5 1.0) (log x)))
double code(double x) {
return fma(x, fma(x, 0.5, 1.0), log(x));
}
function code(x) return fma(x, fma(x, 0.5, 1.0), log(x)) end
code[x_] := N[(x * N[(x * 0.5 + 1.0), $MachinePrecision] + N[Log[x], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\mathsf{fma}\left(x, \mathsf{fma}\left(x, 0.5, 1\right), \log x\right)
\end{array}
Initial program 100.0%
Taylor expanded in x around 0
cancel-sign-sub-invN/A
metadata-evalN/A
*-lft-identityN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower-log.f6499.7
Applied rewrites99.7%
(FPCore (x) :precision binary64 (+ x (log x)))
double code(double x) {
return x + log(x);
}
real(8) function code(x)
real(8), intent (in) :: x
code = x + log(x)
end function
public static double code(double x) {
return x + Math.log(x);
}
def code(x): return x + math.log(x)
function code(x) return Float64(x + log(x)) end
function tmp = code(x) tmp = x + log(x); end
code[x_] := N[(x + N[Log[x], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
x + \log x
\end{array}
Initial program 100.0%
Taylor expanded in x around 0
sub-negN/A
mul-1-negN/A
remove-double-negN/A
lower-+.f64N/A
lower-log.f6499.4
Applied rewrites99.4%
(FPCore (x) :precision binary64 (log x))
double code(double x) {
return log(x);
}
real(8) function code(x)
real(8), intent (in) :: x
code = log(x)
end function
public static double code(double x) {
return Math.log(x);
}
def code(x): return math.log(x)
function code(x) return log(x) end
function tmp = code(x) tmp = log(x); end
code[x_] := N[Log[x], $MachinePrecision]
\begin{array}{l}
\\
\log x
\end{array}
Initial program 100.0%
Taylor expanded in x around 0
lower-log.f6498.5
Applied rewrites98.5%
(FPCore (x) :precision binary64 (* (* x x) 0.5))
double code(double x) {
return (x * x) * 0.5;
}
real(8) function code(x)
real(8), intent (in) :: x
code = (x * x) * 0.5d0
end function
public static double code(double x) {
return (x * x) * 0.5;
}
def code(x): return (x * x) * 0.5
function code(x) return Float64(Float64(x * x) * 0.5) end
function tmp = code(x) tmp = (x * x) * 0.5; end
code[x_] := N[(N[(x * x), $MachinePrecision] * 0.5), $MachinePrecision]
\begin{array}{l}
\\
\left(x \cdot x\right) \cdot 0.5
\end{array}
Initial program 100.0%
Taylor expanded in x around 0
cancel-sign-sub-invN/A
metadata-evalN/A
*-lft-identityN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower-log.f6499.7
Applied rewrites99.7%
Taylor expanded in x around inf
lower-*.f64N/A
unpow2N/A
lower-*.f642.7
Applied rewrites2.7%
Final simplification2.7%
(FPCore (x) :precision binary64 (* x 0.5))
double code(double x) {
return x * 0.5;
}
real(8) function code(x)
real(8), intent (in) :: x
code = x * 0.5d0
end function
public static double code(double x) {
return x * 0.5;
}
def code(x): return x * 0.5
function code(x) return Float64(x * 0.5) end
function tmp = code(x) tmp = x * 0.5; end
code[x_] := N[(x * 0.5), $MachinePrecision]
\begin{array}{l}
\\
x \cdot 0.5
\end{array}
Initial program 100.0%
Taylor expanded in x around 0
cancel-sign-sub-invN/A
metadata-evalN/A
*-lft-identityN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower-log.f6499.7
Applied rewrites99.7%
Taylor expanded in x around inf
lower-*.f64N/A
unpow2N/A
lower-*.f642.7
Applied rewrites2.7%
associate-*r*N/A
lower-*.f64N/A
lower-*.f642.7
Applied rewrites2.7%
Applied rewrites2.2%
Final simplification2.2%
(FPCore (x) :precision binary64 0.5)
double code(double x) {
return 0.5;
}
real(8) function code(x)
real(8), intent (in) :: x
code = 0.5d0
end function
public static double code(double x) {
return 0.5;
}
def code(x): return 0.5
function code(x) return 0.5 end
function tmp = code(x) tmp = 0.5; end
code[x_] := 0.5
\begin{array}{l}
\\
0.5
\end{array}
Initial program 100.0%
Taylor expanded in x around 0
cancel-sign-sub-invN/A
metadata-evalN/A
*-lft-identityN/A
lower-fma.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower-log.f6499.7
Applied rewrites99.7%
Taylor expanded in x around inf
lower-*.f64N/A
unpow2N/A
lower-*.f642.7
Applied rewrites2.7%
associate-*r*N/A
lower-*.f64N/A
lower-*.f642.7
Applied rewrites2.7%
Applied rewrites1.6%
herbie shell --seed 2024219
(FPCore (x)
:name "neg log"
:precision binary64
(- (log (- (/ 1.0 x) 1.0))))