
(FPCore (lo hi x) :precision binary64 (/ (- x lo) (- hi lo)))
double code(double lo, double hi, double x) {
return (x - lo) / (hi - lo);
}
real(8) function code(lo, hi, x)
real(8), intent (in) :: lo
real(8), intent (in) :: hi
real(8), intent (in) :: x
code = (x - lo) / (hi - lo)
end function
public static double code(double lo, double hi, double x) {
return (x - lo) / (hi - lo);
}
def code(lo, hi, x): return (x - lo) / (hi - lo)
function code(lo, hi, x) return Float64(Float64(x - lo) / Float64(hi - lo)) end
function tmp = code(lo, hi, x) tmp = (x - lo) / (hi - lo); end
code[lo_, hi_, x_] := N[(N[(x - lo), $MachinePrecision] / N[(hi - lo), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{x - lo}{hi - lo}
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 5 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (lo hi x) :precision binary64 (/ (- x lo) (- hi lo)))
double code(double lo, double hi, double x) {
return (x - lo) / (hi - lo);
}
real(8) function code(lo, hi, x)
real(8), intent (in) :: lo
real(8), intent (in) :: hi
real(8), intent (in) :: x
code = (x - lo) / (hi - lo)
end function
public static double code(double lo, double hi, double x) {
return (x - lo) / (hi - lo);
}
def code(lo, hi, x): return (x - lo) / (hi - lo)
function code(lo, hi, x) return Float64(Float64(x - lo) / Float64(hi - lo)) end
function tmp = code(lo, hi, x) tmp = (x - lo) / (hi - lo); end
code[lo_, hi_, x_] := N[(N[(x - lo), $MachinePrecision] / N[(hi - lo), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{x - lo}{hi - lo}
\end{array}
(FPCore (lo hi x) :precision binary64 (- 1.0 (fabs (* hi (/ (- -1.0 (/ hi lo)) lo)))))
double code(double lo, double hi, double x) {
return 1.0 - fabs((hi * ((-1.0 - (hi / lo)) / lo)));
}
real(8) function code(lo, hi, x)
real(8), intent (in) :: lo
real(8), intent (in) :: hi
real(8), intent (in) :: x
code = 1.0d0 - abs((hi * (((-1.0d0) - (hi / lo)) / lo)))
end function
public static double code(double lo, double hi, double x) {
return 1.0 - Math.abs((hi * ((-1.0 - (hi / lo)) / lo)));
}
def code(lo, hi, x): return 1.0 - math.fabs((hi * ((-1.0 - (hi / lo)) / lo)))
function code(lo, hi, x) return Float64(1.0 - abs(Float64(hi * Float64(Float64(-1.0 - Float64(hi / lo)) / lo)))) end
function tmp = code(lo, hi, x) tmp = 1.0 - abs((hi * ((-1.0 - (hi / lo)) / lo))); end
code[lo_, hi_, x_] := N[(1.0 - N[Abs[N[(hi * N[(N[(-1.0 - N[(hi / lo), $MachinePrecision]), $MachinePrecision] / lo), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
1 - \left|hi \cdot \frac{-1 - \frac{hi}{lo}}{lo}\right|
\end{array}
Initial program 3.1%
Taylor expanded in lo around -inf 3.1%
mul-1-neg3.1%
unsub-neg3.1%
+-commutative3.1%
associate-/l*14.7%
fma-define14.7%
Simplified14.7%
add-sqr-sqrt9.5%
sqrt-unprod14.1%
pow214.1%
Applied egg-rr14.1%
unpow214.1%
rem-sqrt-square14.1%
Simplified19.8%
Taylor expanded in x around 0 19.8%
associate-*r/19.8%
mul-1-neg19.8%
distribute-rgt-neg-in19.8%
distribute-neg-in19.8%
metadata-eval19.8%
unsub-neg19.8%
Simplified19.8%
associate-/l*19.8%
Applied egg-rr19.8%
(FPCore (lo hi x) :precision binary64 (+ 1.0 (* hi (/ (- (/ hi lo) -1.0) lo))))
double code(double lo, double hi, double x) {
return 1.0 + (hi * (((hi / lo) - -1.0) / lo));
}
real(8) function code(lo, hi, x)
real(8), intent (in) :: lo
real(8), intent (in) :: hi
real(8), intent (in) :: x
code = 1.0d0 + (hi * (((hi / lo) - (-1.0d0)) / lo))
end function
public static double code(double lo, double hi, double x) {
return 1.0 + (hi * (((hi / lo) - -1.0) / lo));
}
def code(lo, hi, x): return 1.0 + (hi * (((hi / lo) - -1.0) / lo))
function code(lo, hi, x) return Float64(1.0 + Float64(hi * Float64(Float64(Float64(hi / lo) - -1.0) / lo))) end
function tmp = code(lo, hi, x) tmp = 1.0 + (hi * (((hi / lo) - -1.0) / lo)); end
code[lo_, hi_, x_] := N[(1.0 + N[(hi * N[(N[(N[(hi / lo), $MachinePrecision] - -1.0), $MachinePrecision] / lo), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
1 + hi \cdot \frac{\frac{hi}{lo} - -1}{lo}
\end{array}
Initial program 3.1%
Taylor expanded in lo around -inf 3.1%
mul-1-neg3.1%
unsub-neg3.1%
+-commutative3.1%
associate-/l*14.7%
fma-define14.7%
Simplified14.7%
add-sqr-sqrt9.5%
sqrt-unprod14.1%
pow214.1%
Applied egg-rr14.1%
unpow214.1%
rem-sqrt-square14.1%
Simplified19.8%
Taylor expanded in x around 0 19.8%
associate-*r/19.8%
mul-1-neg19.8%
distribute-rgt-neg-in19.8%
distribute-neg-in19.8%
metadata-eval19.8%
unsub-neg19.8%
Simplified19.8%
sub-neg19.8%
add-sqr-sqrt9.5%
fabs-sqr9.5%
add-sqr-sqrt18.9%
associate-/l*18.9%
Applied egg-rr18.9%
Final simplification18.9%
(FPCore (lo hi x) :precision binary64 (* lo (/ (+ -1.0 (/ x hi)) hi)))
double code(double lo, double hi, double x) {
return lo * ((-1.0 + (x / hi)) / hi);
}
real(8) function code(lo, hi, x)
real(8), intent (in) :: lo
real(8), intent (in) :: hi
real(8), intent (in) :: x
code = lo * (((-1.0d0) + (x / hi)) / hi)
end function
public static double code(double lo, double hi, double x) {
return lo * ((-1.0 + (x / hi)) / hi);
}
def code(lo, hi, x): return lo * ((-1.0 + (x / hi)) / hi)
function code(lo, hi, x) return Float64(lo * Float64(Float64(-1.0 + Float64(x / hi)) / hi)) end
function tmp = code(lo, hi, x) tmp = lo * ((-1.0 + (x / hi)) / hi); end
code[lo_, hi_, x_] := N[(lo * N[(N[(-1.0 + N[(x / hi), $MachinePrecision]), $MachinePrecision] / hi), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
lo \cdot \frac{-1 + \frac{x}{hi}}{hi}
\end{array}
Initial program 3.1%
Taylor expanded in lo around 0 18.8%
+-commutative18.8%
mul-1-neg18.8%
unsub-neg18.8%
+-commutative18.8%
mul-1-neg18.8%
unsub-neg18.8%
Simplified18.8%
Taylor expanded in lo around inf 18.8%
rem-square-sqrt9.8%
unpow29.8%
times-frac9.8%
associate-*r/9.8%
*-commutative9.8%
associate-*r/9.8%
rem-square-sqrt18.8%
div-sub18.8%
sub-neg18.8%
metadata-eval18.8%
Simplified18.8%
Final simplification18.8%
(FPCore (lo hi x) :precision binary64 (/ lo (- hi)))
double code(double lo, double hi, double x) {
return lo / -hi;
}
real(8) function code(lo, hi, x)
real(8), intent (in) :: lo
real(8), intent (in) :: hi
real(8), intent (in) :: x
code = lo / -hi
end function
public static double code(double lo, double hi, double x) {
return lo / -hi;
}
def code(lo, hi, x): return lo / -hi
function code(lo, hi, x) return Float64(lo / Float64(-hi)) end
function tmp = code(lo, hi, x) tmp = lo / -hi; end
code[lo_, hi_, x_] := N[(lo / (-hi)), $MachinePrecision]
\begin{array}{l}
\\
\frac{lo}{-hi}
\end{array}
Initial program 3.1%
Taylor expanded in hi around inf 18.7%
Taylor expanded in x around 0 18.8%
neg-mul-118.8%
distribute-neg-frac18.8%
Simplified18.8%
Final simplification18.8%
(FPCore (lo hi x) :precision binary64 1.0)
double code(double lo, double hi, double x) {
return 1.0;
}
real(8) function code(lo, hi, x)
real(8), intent (in) :: lo
real(8), intent (in) :: hi
real(8), intent (in) :: x
code = 1.0d0
end function
public static double code(double lo, double hi, double x) {
return 1.0;
}
def code(lo, hi, x): return 1.0
function code(lo, hi, x) return 1.0 end
function tmp = code(lo, hi, x) tmp = 1.0; end
code[lo_, hi_, x_] := 1.0
\begin{array}{l}
\\
1
\end{array}
Initial program 3.1%
Taylor expanded in lo around inf 18.7%
herbie shell --seed 2024129
(FPCore (lo hi x)
:name "xlohi (overflows)"
:precision binary64
:pre (and (< lo -1e+308) (> hi 1e+308))
(/ (- x lo) (- hi lo)))