
(FPCore (x) :precision binary64 (- x (/ 1.0 3.0)))
double code(double x) {
return x - (1.0 / 3.0);
}
real(8) function code(x)
real(8), intent (in) :: x
code = x - (1.0d0 / 3.0d0)
end function
public static double code(double x) {
return x - (1.0 / 3.0);
}
def code(x): return x - (1.0 / 3.0)
function code(x) return Float64(x - Float64(1.0 / 3.0)) end
function tmp = code(x) tmp = x - (1.0 / 3.0); end
code[x_] := N[(x - N[(1.0 / 3.0), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
x - \frac{1}{3}
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 3 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (x) :precision binary64 (- x (/ 1.0 3.0)))
double code(double x) {
return x - (1.0 / 3.0);
}
real(8) function code(x)
real(8), intent (in) :: x
code = x - (1.0d0 / 3.0d0)
end function
public static double code(double x) {
return x - (1.0 / 3.0);
}
def code(x): return x - (1.0 / 3.0)
function code(x) return Float64(x - Float64(1.0 / 3.0)) end
function tmp = code(x) tmp = x - (1.0 / 3.0); end
code[x_] := N[(x - N[(1.0 / 3.0), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
x - \frac{1}{3}
\end{array}
(FPCore (x) :precision binary64 (+ x -0.3333333333333333))
double code(double x) {
return x + -0.3333333333333333;
}
real(8) function code(x)
real(8), intent (in) :: x
code = x + (-0.3333333333333333d0)
end function
public static double code(double x) {
return x + -0.3333333333333333;
}
def code(x): return x + -0.3333333333333333
function code(x) return Float64(x + -0.3333333333333333) end
function tmp = code(x) tmp = x + -0.3333333333333333; end
code[x_] := N[(x + -0.3333333333333333), $MachinePrecision]
\begin{array}{l}
\\
x + -0.3333333333333333
\end{array}
Initial program 100.0%
sub-neg100.0%
metadata-eval100.0%
metadata-eval100.0%
Simplified100.0%
(FPCore (x) :precision binary64 (if (<= x -0.33) x (if (<= x 0.33) -0.3333333333333333 x)))
double code(double x) {
double tmp;
if (x <= -0.33) {
tmp = x;
} else if (x <= 0.33) {
tmp = -0.3333333333333333;
} else {
tmp = x;
}
return tmp;
}
real(8) function code(x)
real(8), intent (in) :: x
real(8) :: tmp
if (x <= (-0.33d0)) then
tmp = x
else if (x <= 0.33d0) then
tmp = -0.3333333333333333d0
else
tmp = x
end if
code = tmp
end function
public static double code(double x) {
double tmp;
if (x <= -0.33) {
tmp = x;
} else if (x <= 0.33) {
tmp = -0.3333333333333333;
} else {
tmp = x;
}
return tmp;
}
def code(x): tmp = 0 if x <= -0.33: tmp = x elif x <= 0.33: tmp = -0.3333333333333333 else: tmp = x return tmp
function code(x) tmp = 0.0 if (x <= -0.33) tmp = x; elseif (x <= 0.33) tmp = -0.3333333333333333; else tmp = x; end return tmp end
function tmp_2 = code(x) tmp = 0.0; if (x <= -0.33) tmp = x; elseif (x <= 0.33) tmp = -0.3333333333333333; else tmp = x; end tmp_2 = tmp; end
code[x_] := If[LessEqual[x, -0.33], x, If[LessEqual[x, 0.33], -0.3333333333333333, x]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq -0.33:\\
\;\;\;\;x\\
\mathbf{elif}\;x \leq 0.33:\\
\;\;\;\;-0.3333333333333333\\
\mathbf{else}:\\
\;\;\;\;x\\
\end{array}
\end{array}
if x < -0.330000000000000016 or 0.330000000000000016 < x Initial program 100.0%
sub-neg100.0%
metadata-eval100.0%
metadata-eval100.0%
Simplified100.0%
Taylor expanded in x around inf 97.6%
if -0.330000000000000016 < x < 0.330000000000000016Initial program 100.0%
sub-neg100.0%
metadata-eval100.0%
metadata-eval100.0%
Simplified100.0%
Taylor expanded in x around 0 98.0%
(FPCore (x) :precision binary64 -0.3333333333333333)
double code(double x) {
return -0.3333333333333333;
}
real(8) function code(x)
real(8), intent (in) :: x
code = -0.3333333333333333d0
end function
public static double code(double x) {
return -0.3333333333333333;
}
def code(x): return -0.3333333333333333
function code(x) return -0.3333333333333333 end
function tmp = code(x) tmp = -0.3333333333333333; end
code[x_] := -0.3333333333333333
\begin{array}{l}
\\
-0.3333333333333333
\end{array}
Initial program 100.0%
sub-neg100.0%
metadata-eval100.0%
metadata-eval100.0%
Simplified100.0%
Taylor expanded in x around 0 51.0%
herbie shell --seed 2024106
(FPCore (x)
:name "Data.Colour.RGBSpace.HSL:hsl from colour-2.3.3, G"
:precision binary64
(- x (/ 1.0 3.0)))