
(FPCore (t) :precision binary64 (- (* 1.7e+308 t) 1.7e+308))
double code(double t) {
return (1.7e+308 * t) - 1.7e+308;
}
real(8) function code(t)
real(8), intent (in) :: t
code = (1.7d+308 * t) - 1.7d+308
end function
public static double code(double t) {
return (1.7e+308 * t) - 1.7e+308;
}
def code(t): return (1.7e+308 * t) - 1.7e+308
function code(t) return Float64(Float64(1.7e+308 * t) - 1.7e+308) end
function tmp = code(t) tmp = (1.7e+308 * t) - 1.7e+308; end
code[t_] := N[(N[(1.7e+308 * t), $MachinePrecision] - 1.7e+308), $MachinePrecision]
\begin{array}{l}
\\
1.7 \cdot 10^{+308} \cdot t - 1.7 \cdot 10^{+308}
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 4 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (t) :precision binary64 (- (* 1.7e+308 t) 1.7e+308))
double code(double t) {
return (1.7e+308 * t) - 1.7e+308;
}
real(8) function code(t)
real(8), intent (in) :: t
code = (1.7d+308 * t) - 1.7d+308
end function
public static double code(double t) {
return (1.7e+308 * t) - 1.7e+308;
}
def code(t): return (1.7e+308 * t) - 1.7e+308
function code(t) return Float64(Float64(1.7e+308 * t) - 1.7e+308) end
function tmp = code(t) tmp = (1.7e+308 * t) - 1.7e+308; end
code[t_] := N[(N[(1.7e+308 * t), $MachinePrecision] - 1.7e+308), $MachinePrecision]
\begin{array}{l}
\\
1.7 \cdot 10^{+308} \cdot t - 1.7 \cdot 10^{+308}
\end{array}
(FPCore (t) :precision binary64 (fma 1.7e+308 t -1.7e+308))
double code(double t) {
return fma(1.7e+308, t, -1.7e+308);
}
function code(t) return fma(1.7e+308, t, -1.7e+308) end
code[t_] := N[(1.7e+308 * t + -1.7e+308), $MachinePrecision]
\begin{array}{l}
\\
\mathsf{fma}\left(1.7 \cdot 10^{+308}, t, -1.7 \cdot 10^{+308}\right)
\end{array}
Initial program 36.7%
fma-neg99.6%
metadata-eval99.6%
Simplified99.6%
Final simplification99.6%
(FPCore (t) :precision binary64 (if (<= t 2.03) 1.7e+308 (* 1.7e+308 t)))
double code(double t) {
double tmp;
if (t <= 2.03) {
tmp = 1.7e+308;
} else {
tmp = 1.7e+308 * t;
}
return tmp;
}
real(8) function code(t)
real(8), intent (in) :: t
real(8) :: tmp
if (t <= 2.03d0) then
tmp = 1.7d+308
else
tmp = 1.7d+308 * t
end if
code = tmp
end function
public static double code(double t) {
double tmp;
if (t <= 2.03) {
tmp = 1.7e+308;
} else {
tmp = 1.7e+308 * t;
}
return tmp;
}
def code(t): tmp = 0 if t <= 2.03: tmp = 1.7e+308 else: tmp = 1.7e+308 * t return tmp
function code(t) tmp = 0.0 if (t <= 2.03) tmp = 1.7e+308; else tmp = Float64(1.7e+308 * t); end return tmp end
function tmp_2 = code(t) tmp = 0.0; if (t <= 2.03) tmp = 1.7e+308; else tmp = 1.7e+308 * t; end tmp_2 = tmp; end
code[t_] := If[LessEqual[t, 2.03], 1.7e+308, N[(1.7e+308 * t), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;t \leq 2.03:\\
\;\;\;\;1.7 \cdot 10^{+308}\\
\mathbf{else}:\\
\;\;\;\;1.7 \cdot 10^{+308} \cdot t\\
\end{array}
\end{array}
if t < 2.0299999999999998Initial program 22.7%
fma-neg99.5%
metadata-eval99.5%
add-cube-cbrt97.7%
pow297.7%
Applied egg-rr97.7%
pow-plus97.7%
metadata-eval97.7%
sqr-pow97.6%
pow-prod-down97.6%
unpow297.6%
add-cube-cbrt96.6%
pow396.6%
pow-pow96.6%
unpow296.6%
cbrt-prod96.3%
pow296.3%
metadata-eval96.3%
metadata-eval96.3%
Applied egg-rr96.3%
fma-def22.7%
*-commutative22.7%
fma-udef96.3%
Simplified96.3%
Taylor expanded in t around 0 25.0%
if 2.0299999999999998 < t Initial program 78.6%
Taylor expanded in t around inf 78.6%
Final simplification38.4%
(FPCore (t) :precision binary64 -1.7e+308)
double code(double t) {
return -1.7e+308;
}
real(8) function code(t)
real(8), intent (in) :: t
code = -1.7d+308
end function
public static double code(double t) {
return -1.7e+308;
}
def code(t): return -1.7e+308
function code(t) return -1.7e+308 end
function tmp = code(t) tmp = -1.7e+308; end
code[t_] := -1.7e+308
\begin{array}{l}
\\
-1.7 \cdot 10^{+308}
\end{array}
Initial program 36.7%
Taylor expanded in t around 0 0.0%
Final simplification0.0%
(FPCore (t) :precision binary64 1.7e+308)
double code(double t) {
return 1.7e+308;
}
real(8) function code(t)
real(8), intent (in) :: t
code = 1.7d+308
end function
public static double code(double t) {
return 1.7e+308;
}
def code(t): return 1.7e+308
function code(t) return 1.7e+308 end
function tmp = code(t) tmp = 1.7e+308; end
code[t_] := 1.7e+308
\begin{array}{l}
\\
1.7 \cdot 10^{+308}
\end{array}
Initial program 36.7%
fma-neg99.6%
metadata-eval99.6%
add-cube-cbrt98.1%
pow298.1%
Applied egg-rr98.1%
pow-plus98.1%
metadata-eval98.1%
sqr-pow98.0%
pow-prod-down98.0%
unpow298.0%
add-cube-cbrt97.2%
pow397.2%
pow-pow97.2%
unpow297.2%
cbrt-prod96.9%
pow296.9%
metadata-eval96.9%
metadata-eval96.9%
Applied egg-rr96.9%
fma-def36.7%
*-commutative36.7%
fma-udef96.9%
Simplified96.9%
Taylor expanded in t around 0 24.8%
Final simplification24.8%
(FPCore (t) :precision binary64 (fma 1.7e+308 t (- 1.7e+308)))
double code(double t) {
return fma(1.7e+308, t, -1.7e+308);
}
function code(t) return fma(1.7e+308, t, Float64(-1.7e+308)) end
code[t_] := N[(1.7e+308 * t + (-1.7e+308)), $MachinePrecision]
\begin{array}{l}
\\
\mathsf{fma}\left(1.7 \cdot 10^{+308}, t, -1.7 \cdot 10^{+308}\right)
\end{array}
herbie shell --seed 2024041
(FPCore (t)
:name "fma_test2"
:precision binary64
:pre (and (<= 1.9 t) (<= t 2.1))
:herbie-target
(fma 1.7e+308 t (- 1.7e+308))
(- (* 1.7e+308 t) 1.7e+308))