
(FPCore (x y z t) :precision binary64 (+ x (* (- y x) (/ z t))))
double code(double x, double y, double z, double t) {
return x + ((y - x) * (z / t));
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
code = x + ((y - x) * (z / t))
end function
public static double code(double x, double y, double z, double t) {
return x + ((y - x) * (z / t));
}
def code(x, y, z, t): return x + ((y - x) * (z / t))
function code(x, y, z, t) return Float64(x + Float64(Float64(y - x) * Float64(z / t))) end
function tmp = code(x, y, z, t) tmp = x + ((y - x) * (z / t)); end
code[x_, y_, z_, t_] := N[(x + N[(N[(y - x), $MachinePrecision] * N[(z / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
x + \left(y - x\right) \cdot \frac{z}{t}
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 9 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (x y z t) :precision binary64 (+ x (* (- y x) (/ z t))))
double code(double x, double y, double z, double t) {
return x + ((y - x) * (z / t));
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
code = x + ((y - x) * (z / t))
end function
public static double code(double x, double y, double z, double t) {
return x + ((y - x) * (z / t));
}
def code(x, y, z, t): return x + ((y - x) * (z / t))
function code(x, y, z, t) return Float64(x + Float64(Float64(y - x) * Float64(z / t))) end
function tmp = code(x, y, z, t) tmp = x + ((y - x) * (z / t)); end
code[x_, y_, z_, t_] := N[(x + N[(N[(y - x), $MachinePrecision] * N[(z / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
x + \left(y - x\right) \cdot \frac{z}{t}
\end{array}
(FPCore (x y z t) :precision binary64 (fma (/ z t) (- y x) x))
double code(double x, double y, double z, double t) {
return fma((z / t), (y - x), x);
}
function code(x, y, z, t) return fma(Float64(z / t), Float64(y - x), x) end
code[x_, y_, z_, t_] := N[(N[(z / t), $MachinePrecision] * N[(y - x), $MachinePrecision] + x), $MachinePrecision]
\begin{array}{l}
\\
\mathsf{fma}\left(\frac{z}{t}, y - x, x\right)
\end{array}
Initial program 98.3%
lift-+.f64N/A
+-commutativeN/A
lift-*.f64N/A
*-commutativeN/A
lower-fma.f6498.3
Applied rewrites98.3%
(FPCore (x y z t) :precision binary64 (if (or (<= (/ z t) -4e-33) (not (<= (/ z t) 50000.0))) (/ (* (- y x) z) t) (fma (/ y t) z x)))
double code(double x, double y, double z, double t) {
double tmp;
if (((z / t) <= -4e-33) || !((z / t) <= 50000.0)) {
tmp = ((y - x) * z) / t;
} else {
tmp = fma((y / t), z, x);
}
return tmp;
}
function code(x, y, z, t) tmp = 0.0 if ((Float64(z / t) <= -4e-33) || !(Float64(z / t) <= 50000.0)) tmp = Float64(Float64(Float64(y - x) * z) / t); else tmp = fma(Float64(y / t), z, x); end return tmp end
code[x_, y_, z_, t_] := If[Or[LessEqual[N[(z / t), $MachinePrecision], -4e-33], N[Not[LessEqual[N[(z / t), $MachinePrecision], 50000.0]], $MachinePrecision]], N[(N[(N[(y - x), $MachinePrecision] * z), $MachinePrecision] / t), $MachinePrecision], N[(N[(y / t), $MachinePrecision] * z + x), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\frac{z}{t} \leq -4 \cdot 10^{-33} \lor \neg \left(\frac{z}{t} \leq 50000\right):\\
\;\;\;\;\frac{\left(y - x\right) \cdot z}{t}\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(\frac{y}{t}, z, x\right)\\
\end{array}
\end{array}
if (/.f64 z t) < -4.0000000000000002e-33 or 5e4 < (/.f64 z t) Initial program 96.9%
Taylor expanded in z around inf
div-subN/A
*-commutativeN/A
lower-*.f64N/A
lower-/.f64N/A
lower--.f6489.4
Applied rewrites89.4%
Applied rewrites92.4%
if -4.0000000000000002e-33 < (/.f64 z t) < 5e4Initial program 99.8%
lift-+.f64N/A
+-commutativeN/A
lift-*.f64N/A
*-commutativeN/A
lower-fma.f6499.9
Applied rewrites99.9%
lift-fma.f64N/A
*-commutativeN/A
lift-/.f64N/A
associate-*r/N/A
associate-*l/N/A
lift-/.f64N/A
lower-fma.f6494.4
Applied rewrites94.4%
Taylor expanded in x around 0
lower-/.f6495.5
Applied rewrites95.5%
Final simplification93.9%
(FPCore (x y z t) :precision binary64 (if (or (<= (/ z t) -1000000.0) (not (<= (/ z t) 0.0005))) (* (/ (- y x) t) z) (fma (/ y t) z x)))
double code(double x, double y, double z, double t) {
double tmp;
if (((z / t) <= -1000000.0) || !((z / t) <= 0.0005)) {
tmp = ((y - x) / t) * z;
} else {
tmp = fma((y / t), z, x);
}
return tmp;
}
function code(x, y, z, t) tmp = 0.0 if ((Float64(z / t) <= -1000000.0) || !(Float64(z / t) <= 0.0005)) tmp = Float64(Float64(Float64(y - x) / t) * z); else tmp = fma(Float64(y / t), z, x); end return tmp end
code[x_, y_, z_, t_] := If[Or[LessEqual[N[(z / t), $MachinePrecision], -1000000.0], N[Not[LessEqual[N[(z / t), $MachinePrecision], 0.0005]], $MachinePrecision]], N[(N[(N[(y - x), $MachinePrecision] / t), $MachinePrecision] * z), $MachinePrecision], N[(N[(y / t), $MachinePrecision] * z + x), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\frac{z}{t} \leq -1000000 \lor \neg \left(\frac{z}{t} \leq 0.0005\right):\\
\;\;\;\;\frac{y - x}{t} \cdot z\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(\frac{y}{t}, z, x\right)\\
\end{array}
\end{array}
if (/.f64 z t) < -1e6 or 5.0000000000000001e-4 < (/.f64 z t) Initial program 96.8%
Taylor expanded in z around inf
div-subN/A
*-commutativeN/A
lower-*.f64N/A
lower-/.f64N/A
lower--.f6492.6
Applied rewrites92.6%
if -1e6 < (/.f64 z t) < 5.0000000000000001e-4Initial program 99.8%
lift-+.f64N/A
+-commutativeN/A
lift-*.f64N/A
*-commutativeN/A
lower-fma.f6499.9
Applied rewrites99.9%
lift-fma.f64N/A
*-commutativeN/A
lift-/.f64N/A
associate-*r/N/A
associate-*l/N/A
lift-/.f64N/A
lower-fma.f6492.4
Applied rewrites92.4%
Taylor expanded in x around 0
lower-/.f6492.9
Applied rewrites92.9%
Final simplification92.7%
(FPCore (x y z t) :precision binary64 (if (<= (/ z t) 200000000.0) (fma (/ y t) z x) (if (<= (/ z t) 1e+236) (/ (* (- x) z) t) (* y (/ z t)))))
double code(double x, double y, double z, double t) {
double tmp;
if ((z / t) <= 200000000.0) {
tmp = fma((y / t), z, x);
} else if ((z / t) <= 1e+236) {
tmp = (-x * z) / t;
} else {
tmp = y * (z / t);
}
return tmp;
}
function code(x, y, z, t) tmp = 0.0 if (Float64(z / t) <= 200000000.0) tmp = fma(Float64(y / t), z, x); elseif (Float64(z / t) <= 1e+236) tmp = Float64(Float64(Float64(-x) * z) / t); else tmp = Float64(y * Float64(z / t)); end return tmp end
code[x_, y_, z_, t_] := If[LessEqual[N[(z / t), $MachinePrecision], 200000000.0], N[(N[(y / t), $MachinePrecision] * z + x), $MachinePrecision], If[LessEqual[N[(z / t), $MachinePrecision], 1e+236], N[(N[((-x) * z), $MachinePrecision] / t), $MachinePrecision], N[(y * N[(z / t), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\frac{z}{t} \leq 200000000:\\
\;\;\;\;\mathsf{fma}\left(\frac{y}{t}, z, x\right)\\
\mathbf{elif}\;\frac{z}{t} \leq 10^{+236}:\\
\;\;\;\;\frac{\left(-x\right) \cdot z}{t}\\
\mathbf{else}:\\
\;\;\;\;y \cdot \frac{z}{t}\\
\end{array}
\end{array}
if (/.f64 z t) < 2e8Initial program 98.3%
lift-+.f64N/A
+-commutativeN/A
lift-*.f64N/A
*-commutativeN/A
lower-fma.f6498.3
Applied rewrites98.3%
lift-fma.f64N/A
*-commutativeN/A
lift-/.f64N/A
associate-*r/N/A
associate-*l/N/A
lift-/.f64N/A
lower-fma.f6492.7
Applied rewrites92.7%
Taylor expanded in x around 0
lower-/.f6482.9
Applied rewrites82.9%
if 2e8 < (/.f64 z t) < 1.00000000000000005e236Initial program 99.6%
Taylor expanded in z around inf
div-subN/A
*-commutativeN/A
lower-*.f64N/A
lower-/.f64N/A
lower--.f6486.7
Applied rewrites86.7%
Applied rewrites91.8%
Taylor expanded in x around inf
Applied rewrites65.0%
if 1.00000000000000005e236 < (/.f64 z t) Initial program 97.0%
Taylor expanded in x around 0
associate-*l/N/A
lower-*.f64N/A
lower-/.f6470.5
Applied rewrites70.5%
Applied rewrites78.9%
(FPCore (x y z t) :precision binary64 (if (<= (/ z t) 200000000.0) (fma (/ y t) z x) (if (<= (/ z t) 1e+236) (* (/ (- x) t) z) (* y (/ z t)))))
double code(double x, double y, double z, double t) {
double tmp;
if ((z / t) <= 200000000.0) {
tmp = fma((y / t), z, x);
} else if ((z / t) <= 1e+236) {
tmp = (-x / t) * z;
} else {
tmp = y * (z / t);
}
return tmp;
}
function code(x, y, z, t) tmp = 0.0 if (Float64(z / t) <= 200000000.0) tmp = fma(Float64(y / t), z, x); elseif (Float64(z / t) <= 1e+236) tmp = Float64(Float64(Float64(-x) / t) * z); else tmp = Float64(y * Float64(z / t)); end return tmp end
code[x_, y_, z_, t_] := If[LessEqual[N[(z / t), $MachinePrecision], 200000000.0], N[(N[(y / t), $MachinePrecision] * z + x), $MachinePrecision], If[LessEqual[N[(z / t), $MachinePrecision], 1e+236], N[(N[((-x) / t), $MachinePrecision] * z), $MachinePrecision], N[(y * N[(z / t), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\frac{z}{t} \leq 200000000:\\
\;\;\;\;\mathsf{fma}\left(\frac{y}{t}, z, x\right)\\
\mathbf{elif}\;\frac{z}{t} \leq 10^{+236}:\\
\;\;\;\;\frac{-x}{t} \cdot z\\
\mathbf{else}:\\
\;\;\;\;y \cdot \frac{z}{t}\\
\end{array}
\end{array}
if (/.f64 z t) < 2e8Initial program 98.3%
lift-+.f64N/A
+-commutativeN/A
lift-*.f64N/A
*-commutativeN/A
lower-fma.f6498.3
Applied rewrites98.3%
lift-fma.f64N/A
*-commutativeN/A
lift-/.f64N/A
associate-*r/N/A
associate-*l/N/A
lift-/.f64N/A
lower-fma.f6492.7
Applied rewrites92.7%
Taylor expanded in x around 0
lower-/.f6482.9
Applied rewrites82.9%
if 2e8 < (/.f64 z t) < 1.00000000000000005e236Initial program 99.6%
Taylor expanded in z around inf
div-subN/A
*-commutativeN/A
lower-*.f64N/A
lower-/.f64N/A
lower--.f6486.7
Applied rewrites86.7%
Taylor expanded in x around inf
Applied rewrites57.4%
if 1.00000000000000005e236 < (/.f64 z t) Initial program 97.0%
Taylor expanded in x around 0
associate-*l/N/A
lower-*.f64N/A
lower-/.f6470.5
Applied rewrites70.5%
Applied rewrites78.9%
(FPCore (x y z t) :precision binary64 (if (<= (/ z t) 50000.0) (fma (/ y t) z x) (* y (/ z t))))
double code(double x, double y, double z, double t) {
double tmp;
if ((z / t) <= 50000.0) {
tmp = fma((y / t), z, x);
} else {
tmp = y * (z / t);
}
return tmp;
}
function code(x, y, z, t) tmp = 0.0 if (Float64(z / t) <= 50000.0) tmp = fma(Float64(y / t), z, x); else tmp = Float64(y * Float64(z / t)); end return tmp end
code[x_, y_, z_, t_] := If[LessEqual[N[(z / t), $MachinePrecision], 50000.0], N[(N[(y / t), $MachinePrecision] * z + x), $MachinePrecision], N[(y * N[(z / t), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;\frac{z}{t} \leq 50000:\\
\;\;\;\;\mathsf{fma}\left(\frac{y}{t}, z, x\right)\\
\mathbf{else}:\\
\;\;\;\;y \cdot \frac{z}{t}\\
\end{array}
\end{array}
if (/.f64 z t) < 5e4Initial program 98.3%
lift-+.f64N/A
+-commutativeN/A
lift-*.f64N/A
*-commutativeN/A
lower-fma.f6498.3
Applied rewrites98.3%
lift-fma.f64N/A
*-commutativeN/A
lift-/.f64N/A
associate-*r/N/A
associate-*l/N/A
lift-/.f64N/A
lower-fma.f6492.7
Applied rewrites92.7%
Taylor expanded in x around 0
lower-/.f6483.2
Applied rewrites83.2%
if 5e4 < (/.f64 z t) Initial program 98.3%
Taylor expanded in x around 0
associate-*l/N/A
lower-*.f64N/A
lower-/.f6450.2
Applied rewrites50.2%
Applied rewrites55.6%
(FPCore (x y z t) :precision binary64 (if (or (<= x -9.5e-60) (not (<= x 6.8e+34))) (* (- 1.0 (/ z t)) x) (fma (/ y t) z x)))
double code(double x, double y, double z, double t) {
double tmp;
if ((x <= -9.5e-60) || !(x <= 6.8e+34)) {
tmp = (1.0 - (z / t)) * x;
} else {
tmp = fma((y / t), z, x);
}
return tmp;
}
function code(x, y, z, t) tmp = 0.0 if ((x <= -9.5e-60) || !(x <= 6.8e+34)) tmp = Float64(Float64(1.0 - Float64(z / t)) * x); else tmp = fma(Float64(y / t), z, x); end return tmp end
code[x_, y_, z_, t_] := If[Or[LessEqual[x, -9.5e-60], N[Not[LessEqual[x, 6.8e+34]], $MachinePrecision]], N[(N[(1.0 - N[(z / t), $MachinePrecision]), $MachinePrecision] * x), $MachinePrecision], N[(N[(y / t), $MachinePrecision] * z + x), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq -9.5 \cdot 10^{-60} \lor \neg \left(x \leq 6.8 \cdot 10^{+34}\right):\\
\;\;\;\;\left(1 - \frac{z}{t}\right) \cdot x\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(\frac{y}{t}, z, x\right)\\
\end{array}
\end{array}
if x < -9.49999999999999958e-60 or 6.7999999999999999e34 < x Initial program 99.2%
Taylor expanded in x around inf
*-commutativeN/A
lower-*.f64N/A
fp-cancel-sign-sub-invN/A
metadata-evalN/A
*-lft-identityN/A
lower--.f64N/A
lower-/.f6487.0
Applied rewrites87.0%
if -9.49999999999999958e-60 < x < 6.7999999999999999e34Initial program 97.4%
lift-+.f64N/A
+-commutativeN/A
lift-*.f64N/A
*-commutativeN/A
lower-fma.f6497.4
Applied rewrites97.4%
lift-fma.f64N/A
*-commutativeN/A
lift-/.f64N/A
associate-*r/N/A
associate-*l/N/A
lift-/.f64N/A
lower-fma.f6494.2
Applied rewrites94.2%
Taylor expanded in x around 0
lower-/.f6489.4
Applied rewrites89.4%
Final simplification88.1%
(FPCore (x y z t) :precision binary64 (fma z (/ (- y x) t) x))
double code(double x, double y, double z, double t) {
return fma(z, ((y - x) / t), x);
}
function code(x, y, z, t) return fma(z, Float64(Float64(y - x) / t), x) end
code[x_, y_, z_, t_] := N[(z * N[(N[(y - x), $MachinePrecision] / t), $MachinePrecision] + x), $MachinePrecision]
\begin{array}{l}
\\
\mathsf{fma}\left(z, \frac{y - x}{t}, x\right)
\end{array}
Initial program 98.3%
lift-+.f64N/A
+-commutativeN/A
lift-*.f64N/A
*-commutativeN/A
lift-/.f64N/A
associate-*l/N/A
associate-/l*N/A
lower-fma.f64N/A
lower-/.f6492.8
Applied rewrites92.8%
(FPCore (x y z t) :precision binary64 (* y (/ z t)))
double code(double x, double y, double z, double t) {
return y * (z / t);
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
code = y * (z / t)
end function
public static double code(double x, double y, double z, double t) {
return y * (z / t);
}
def code(x, y, z, t): return y * (z / t)
function code(x, y, z, t) return Float64(y * Float64(z / t)) end
function tmp = code(x, y, z, t) tmp = y * (z / t); end
code[x_, y_, z_, t_] := N[(y * N[(z / t), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
y \cdot \frac{z}{t}
\end{array}
Initial program 98.3%
Taylor expanded in x around 0
associate-*l/N/A
lower-*.f64N/A
lower-/.f6442.2
Applied rewrites42.2%
Applied rewrites45.8%
(FPCore (x y z t)
:precision binary64
(let* ((t_1 (* (- y x) (/ z t))) (t_2 (+ x (/ (- y x) (/ t z)))))
(if (< t_1 -1013646692435.8867)
t_2
(if (< t_1 0.0) (+ x (/ (* (- y x) z) t)) t_2))))
double code(double x, double y, double z, double t) {
double t_1 = (y - x) * (z / t);
double t_2 = x + ((y - x) / (t / z));
double tmp;
if (t_1 < -1013646692435.8867) {
tmp = t_2;
} else if (t_1 < 0.0) {
tmp = x + (((y - x) * z) / t);
} else {
tmp = t_2;
}
return tmp;
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8) :: t_1
real(8) :: t_2
real(8) :: tmp
t_1 = (y - x) * (z / t)
t_2 = x + ((y - x) / (t / z))
if (t_1 < (-1013646692435.8867d0)) then
tmp = t_2
else if (t_1 < 0.0d0) then
tmp = x + (((y - x) * z) / t)
else
tmp = t_2
end if
code = tmp
end function
public static double code(double x, double y, double z, double t) {
double t_1 = (y - x) * (z / t);
double t_2 = x + ((y - x) / (t / z));
double tmp;
if (t_1 < -1013646692435.8867) {
tmp = t_2;
} else if (t_1 < 0.0) {
tmp = x + (((y - x) * z) / t);
} else {
tmp = t_2;
}
return tmp;
}
def code(x, y, z, t): t_1 = (y - x) * (z / t) t_2 = x + ((y - x) / (t / z)) tmp = 0 if t_1 < -1013646692435.8867: tmp = t_2 elif t_1 < 0.0: tmp = x + (((y - x) * z) / t) else: tmp = t_2 return tmp
function code(x, y, z, t) t_1 = Float64(Float64(y - x) * Float64(z / t)) t_2 = Float64(x + Float64(Float64(y - x) / Float64(t / z))) tmp = 0.0 if (t_1 < -1013646692435.8867) tmp = t_2; elseif (t_1 < 0.0) tmp = Float64(x + Float64(Float64(Float64(y - x) * z) / t)); else tmp = t_2; end return tmp end
function tmp_2 = code(x, y, z, t) t_1 = (y - x) * (z / t); t_2 = x + ((y - x) / (t / z)); tmp = 0.0; if (t_1 < -1013646692435.8867) tmp = t_2; elseif (t_1 < 0.0) tmp = x + (((y - x) * z) / t); else tmp = t_2; end tmp_2 = tmp; end
code[x_, y_, z_, t_] := Block[{t$95$1 = N[(N[(y - x), $MachinePrecision] * N[(z / t), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(x + N[(N[(y - x), $MachinePrecision] / N[(t / z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[Less[t$95$1, -1013646692435.8867], t$95$2, If[Less[t$95$1, 0.0], N[(x + N[(N[(N[(y - x), $MachinePrecision] * z), $MachinePrecision] / t), $MachinePrecision]), $MachinePrecision], t$95$2]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \left(y - x\right) \cdot \frac{z}{t}\\
t_2 := x + \frac{y - x}{\frac{t}{z}}\\
\mathbf{if}\;t\_1 < -1013646692435.8867:\\
\;\;\;\;t\_2\\
\mathbf{elif}\;t\_1 < 0:\\
\;\;\;\;x + \frac{\left(y - x\right) \cdot z}{t}\\
\mathbf{else}:\\
\;\;\;\;t\_2\\
\end{array}
\end{array}
herbie shell --seed 2024329
(FPCore (x y z t)
:name "Graphics.Rendering.Plot.Render.Plot.Axis:tickPosition from plot-0.2.3.4"
:precision binary64
:alt
(! :herbie-platform default (if (< (* (- y x) (/ z t)) -10136466924358867/10000) (+ x (/ (- y x) (/ t z))) (if (< (* (- y x) (/ z t)) 0) (+ x (/ (* (- y x) z) t)) (+ x (/ (- y x) (/ t z))))))
(+ x (* (- y x) (/ z t))))