
(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 11 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 (+ x (/ (- y x) (/ t z))))
double code(double x, double y, double z, double t) {
return x + ((y - x) / (t / z));
}
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) / (t / z))
end function
public static double code(double x, double y, double z, double t) {
return x + ((y - x) / (t / z));
}
def code(x, y, z, t): return x + ((y - x) / (t / z))
function code(x, y, z, t) return Float64(x + Float64(Float64(y - x) / Float64(t / z))) end
function tmp = code(x, y, z, t) tmp = x + ((y - x) / (t / z)); end
code[x_, y_, z_, t_] := N[(x + N[(N[(y - x), $MachinePrecision] / N[(t / z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
x + \frac{y - x}{\frac{t}{z}}
\end{array}
Initial program 97.9%
lift-*.f64N/A
lift-/.f64N/A
clear-numN/A
un-div-invN/A
lower-/.f64N/A
lower-/.f6497.9
Applied rewrites97.9%
(FPCore (x y z t)
:precision binary64
(let* ((t_1 (* (- y x) (/ z t))))
(if (<= (/ z t) -40000000.0)
t_1
(if (<= (/ z t) 5e-15) (+ x (/ (* y z) t)) t_1))))
double code(double x, double y, double z, double t) {
double t_1 = (y - x) * (z / t);
double tmp;
if ((z / t) <= -40000000.0) {
tmp = t_1;
} else if ((z / t) <= 5e-15) {
tmp = x + ((y * z) / t);
} else {
tmp = t_1;
}
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) :: tmp
t_1 = (y - x) * (z / t)
if ((z / t) <= (-40000000.0d0)) then
tmp = t_1
else if ((z / t) <= 5d-15) then
tmp = x + ((y * z) / t)
else
tmp = t_1
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 tmp;
if ((z / t) <= -40000000.0) {
tmp = t_1;
} else if ((z / t) <= 5e-15) {
tmp = x + ((y * z) / t);
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t): t_1 = (y - x) * (z / t) tmp = 0 if (z / t) <= -40000000.0: tmp = t_1 elif (z / t) <= 5e-15: tmp = x + ((y * z) / t) else: tmp = t_1 return tmp
function code(x, y, z, t) t_1 = Float64(Float64(y - x) * Float64(z / t)) tmp = 0.0 if (Float64(z / t) <= -40000000.0) tmp = t_1; elseif (Float64(z / t) <= 5e-15) tmp = Float64(x + Float64(Float64(y * z) / t)); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t) t_1 = (y - x) * (z / t); tmp = 0.0; if ((z / t) <= -40000000.0) tmp = t_1; elseif ((z / t) <= 5e-15) tmp = x + ((y * z) / t); else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_] := Block[{t$95$1 = N[(N[(y - x), $MachinePrecision] * N[(z / t), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[N[(z / t), $MachinePrecision], -40000000.0], t$95$1, If[LessEqual[N[(z / t), $MachinePrecision], 5e-15], N[(x + N[(N[(y * z), $MachinePrecision] / t), $MachinePrecision]), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \left(y - x\right) \cdot \frac{z}{t}\\
\mathbf{if}\;\frac{z}{t} \leq -40000000:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;\frac{z}{t} \leq 5 \cdot 10^{-15}:\\
\;\;\;\;x + \frac{y \cdot z}{t}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if (/.f64 z t) < -4e7 or 4.99999999999999999e-15 < (/.f64 z t) Initial program 97.8%
Taylor expanded in z around inf
div-subN/A
lower-*.f64N/A
lower-/.f64N/A
lower--.f6494.1
Applied rewrites94.1%
Applied rewrites96.6%
if -4e7 < (/.f64 z t) < 4.99999999999999999e-15Initial program 98.1%
Taylor expanded in y around inf
lower-/.f64N/A
lower-*.f6497.9
Applied rewrites97.9%
Final simplification97.2%
(FPCore (x y z t)
:precision binary64
(let* ((t_1 (* (- y x) (/ z t))))
(if (<= (/ z t) -4.5e-13)
t_1
(if (<= (/ z t) 5e-15) (fma x (/ z (- t)) x) t_1))))
double code(double x, double y, double z, double t) {
double t_1 = (y - x) * (z / t);
double tmp;
if ((z / t) <= -4.5e-13) {
tmp = t_1;
} else if ((z / t) <= 5e-15) {
tmp = fma(x, (z / -t), x);
} else {
tmp = t_1;
}
return tmp;
}
function code(x, y, z, t) t_1 = Float64(Float64(y - x) * Float64(z / t)) tmp = 0.0 if (Float64(z / t) <= -4.5e-13) tmp = t_1; elseif (Float64(z / t) <= 5e-15) tmp = fma(x, Float64(z / Float64(-t)), x); else tmp = t_1; end return tmp end
code[x_, y_, z_, t_] := Block[{t$95$1 = N[(N[(y - x), $MachinePrecision] * N[(z / t), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[N[(z / t), $MachinePrecision], -4.5e-13], t$95$1, If[LessEqual[N[(z / t), $MachinePrecision], 5e-15], N[(x * N[(z / (-t)), $MachinePrecision] + x), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \left(y - x\right) \cdot \frac{z}{t}\\
\mathbf{if}\;\frac{z}{t} \leq -4.5 \cdot 10^{-13}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;\frac{z}{t} \leq 5 \cdot 10^{-15}:\\
\;\;\;\;\mathsf{fma}\left(x, \frac{z}{-t}, x\right)\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if (/.f64 z t) < -4.5e-13 or 4.99999999999999999e-15 < (/.f64 z t) Initial program 97.8%
Taylor expanded in z around inf
div-subN/A
lower-*.f64N/A
lower-/.f64N/A
lower--.f6492.3
Applied rewrites92.3%
Applied rewrites96.1%
if -4.5e-13 < (/.f64 z t) < 4.99999999999999999e-15Initial program 98.0%
lift-*.f64N/A
lift-/.f64N/A
clear-numN/A
un-div-invN/A
lower-/.f64N/A
lower-/.f6497.5
Applied rewrites97.5%
Taylor expanded in x around inf
+-commutativeN/A
distribute-lft-inN/A
*-rgt-identityN/A
lower-fma.f64N/A
mul-1-negN/A
distribute-neg-frac2N/A
lower-/.f64N/A
lower-neg.f6481.6
Applied rewrites81.6%
Final simplification89.6%
(FPCore (x y z t)
:precision binary64
(let* ((t_1 (* (- y x) (/ z t))))
(if (<= (/ z t) -4.5e-13)
t_1
(if (<= (/ z t) 5e-15) (- x (/ (* x z) t)) t_1))))
double code(double x, double y, double z, double t) {
double t_1 = (y - x) * (z / t);
double tmp;
if ((z / t) <= -4.5e-13) {
tmp = t_1;
} else if ((z / t) <= 5e-15) {
tmp = x - ((x * z) / t);
} else {
tmp = t_1;
}
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) :: tmp
t_1 = (y - x) * (z / t)
if ((z / t) <= (-4.5d-13)) then
tmp = t_1
else if ((z / t) <= 5d-15) then
tmp = x - ((x * z) / t)
else
tmp = t_1
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 tmp;
if ((z / t) <= -4.5e-13) {
tmp = t_1;
} else if ((z / t) <= 5e-15) {
tmp = x - ((x * z) / t);
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t): t_1 = (y - x) * (z / t) tmp = 0 if (z / t) <= -4.5e-13: tmp = t_1 elif (z / t) <= 5e-15: tmp = x - ((x * z) / t) else: tmp = t_1 return tmp
function code(x, y, z, t) t_1 = Float64(Float64(y - x) * Float64(z / t)) tmp = 0.0 if (Float64(z / t) <= -4.5e-13) tmp = t_1; elseif (Float64(z / t) <= 5e-15) tmp = Float64(x - Float64(Float64(x * z) / t)); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t) t_1 = (y - x) * (z / t); tmp = 0.0; if ((z / t) <= -4.5e-13) tmp = t_1; elseif ((z / t) <= 5e-15) tmp = x - ((x * z) / t); else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_] := Block[{t$95$1 = N[(N[(y - x), $MachinePrecision] * N[(z / t), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[N[(z / t), $MachinePrecision], -4.5e-13], t$95$1, If[LessEqual[N[(z / t), $MachinePrecision], 5e-15], N[(x - N[(N[(x * z), $MachinePrecision] / t), $MachinePrecision]), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \left(y - x\right) \cdot \frac{z}{t}\\
\mathbf{if}\;\frac{z}{t} \leq -4.5 \cdot 10^{-13}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;\frac{z}{t} \leq 5 \cdot 10^{-15}:\\
\;\;\;\;x - \frac{x \cdot z}{t}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if (/.f64 z t) < -4.5e-13 or 4.99999999999999999e-15 < (/.f64 z t) Initial program 97.8%
Taylor expanded in z around inf
div-subN/A
lower-*.f64N/A
lower-/.f64N/A
lower--.f6492.3
Applied rewrites92.3%
Applied rewrites96.1%
if -4.5e-13 < (/.f64 z t) < 4.99999999999999999e-15Initial program 98.0%
Taylor expanded in x around inf
mul-1-negN/A
unsub-negN/A
distribute-lft-out--N/A
*-rgt-identityN/A
associate-/l*N/A
lower--.f64N/A
lower-/.f64N/A
*-commutativeN/A
lower-*.f6479.9
Applied rewrites79.9%
Final simplification88.8%
(FPCore (x y z t) :precision binary64 (let* ((t_1 (* y (/ z t)))) (if (<= y -1.45e-97) t_1 (if (<= y 3.8e-44) (* z (/ (- x) t)) t_1))))
double code(double x, double y, double z, double t) {
double t_1 = y * (z / t);
double tmp;
if (y <= -1.45e-97) {
tmp = t_1;
} else if (y <= 3.8e-44) {
tmp = z * (-x / t);
} else {
tmp = t_1;
}
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) :: tmp
t_1 = y * (z / t)
if (y <= (-1.45d-97)) then
tmp = t_1
else if (y <= 3.8d-44) then
tmp = z * (-x / t)
else
tmp = t_1
end if
code = tmp
end function
public static double code(double x, double y, double z, double t) {
double t_1 = y * (z / t);
double tmp;
if (y <= -1.45e-97) {
tmp = t_1;
} else if (y <= 3.8e-44) {
tmp = z * (-x / t);
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t): t_1 = y * (z / t) tmp = 0 if y <= -1.45e-97: tmp = t_1 elif y <= 3.8e-44: tmp = z * (-x / t) else: tmp = t_1 return tmp
function code(x, y, z, t) t_1 = Float64(y * Float64(z / t)) tmp = 0.0 if (y <= -1.45e-97) tmp = t_1; elseif (y <= 3.8e-44) tmp = Float64(z * Float64(Float64(-x) / t)); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t) t_1 = y * (z / t); tmp = 0.0; if (y <= -1.45e-97) tmp = t_1; elseif (y <= 3.8e-44) tmp = z * (-x / t); else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_] := Block[{t$95$1 = N[(y * N[(z / t), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -1.45e-97], t$95$1, If[LessEqual[y, 3.8e-44], N[(z * N[((-x) / t), $MachinePrecision]), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := y \cdot \frac{z}{t}\\
\mathbf{if}\;y \leq -1.45 \cdot 10^{-97}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y \leq 3.8 \cdot 10^{-44}:\\
\;\;\;\;z \cdot \frac{-x}{t}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if y < -1.45e-97 or 3.8000000000000001e-44 < y Initial program 97.9%
Taylor expanded in x around 0
lower-/.f64N/A
lower-*.f6460.0
Applied rewrites60.0%
Applied rewrites61.1%
if -1.45e-97 < y < 3.8000000000000001e-44Initial program 98.0%
Taylor expanded in z around inf
div-subN/A
lower-*.f64N/A
lower-/.f64N/A
lower--.f6452.7
Applied rewrites52.7%
Taylor expanded in y around 0
Applied rewrites48.9%
Final simplification56.4%
(FPCore (x y z t) :precision binary64 (let* ((t_1 (* y (/ z t)))) (if (<= y -1.35e-97) t_1 (if (<= y 3.8e-44) (/ (* x z) (- t)) t_1))))
double code(double x, double y, double z, double t) {
double t_1 = y * (z / t);
double tmp;
if (y <= -1.35e-97) {
tmp = t_1;
} else if (y <= 3.8e-44) {
tmp = (x * z) / -t;
} else {
tmp = t_1;
}
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) :: tmp
t_1 = y * (z / t)
if (y <= (-1.35d-97)) then
tmp = t_1
else if (y <= 3.8d-44) then
tmp = (x * z) / -t
else
tmp = t_1
end if
code = tmp
end function
public static double code(double x, double y, double z, double t) {
double t_1 = y * (z / t);
double tmp;
if (y <= -1.35e-97) {
tmp = t_1;
} else if (y <= 3.8e-44) {
tmp = (x * z) / -t;
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t): t_1 = y * (z / t) tmp = 0 if y <= -1.35e-97: tmp = t_1 elif y <= 3.8e-44: tmp = (x * z) / -t else: tmp = t_1 return tmp
function code(x, y, z, t) t_1 = Float64(y * Float64(z / t)) tmp = 0.0 if (y <= -1.35e-97) tmp = t_1; elseif (y <= 3.8e-44) tmp = Float64(Float64(x * z) / Float64(-t)); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t) t_1 = y * (z / t); tmp = 0.0; if (y <= -1.35e-97) tmp = t_1; elseif (y <= 3.8e-44) tmp = (x * z) / -t; else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_] := Block[{t$95$1 = N[(y * N[(z / t), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -1.35e-97], t$95$1, If[LessEqual[y, 3.8e-44], N[(N[(x * z), $MachinePrecision] / (-t)), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := y \cdot \frac{z}{t}\\
\mathbf{if}\;y \leq -1.35 \cdot 10^{-97}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y \leq 3.8 \cdot 10^{-44}:\\
\;\;\;\;\frac{x \cdot z}{-t}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if y < -1.34999999999999993e-97 or 3.8000000000000001e-44 < y Initial program 97.9%
Taylor expanded in x around 0
lower-/.f64N/A
lower-*.f6460.0
Applied rewrites60.0%
Applied rewrites61.1%
if -1.34999999999999993e-97 < y < 3.8000000000000001e-44Initial program 98.0%
Taylor expanded in x around inf
mul-1-negN/A
unsub-negN/A
distribute-lft-out--N/A
*-rgt-identityN/A
associate-/l*N/A
lower--.f64N/A
lower-/.f64N/A
*-commutativeN/A
lower-*.f6490.1
Applied rewrites90.1%
Taylor expanded in z around inf
Applied rewrites48.3%
Final simplification56.2%
(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 97.9%
lift-+.f64N/A
+-commutativeN/A
lift-*.f64N/A
*-commutativeN/A
lower-fma.f6497.9
Applied rewrites97.9%
(FPCore (x y z t) :precision binary64 (* (- y x) (/ z t)))
double code(double x, double y, double z, double t) {
return (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 = (y - x) * (z / t)
end function
public static double code(double x, double y, double z, double t) {
return (y - x) * (z / t);
}
def code(x, y, z, t): return (y - x) * (z / t)
function code(x, y, z, t) return Float64(Float64(y - x) * Float64(z / t)) end
function tmp = code(x, y, z, t) tmp = (y - x) * (z / t); end
code[x_, y_, z_, t_] := N[(N[(y - x), $MachinePrecision] * N[(z / t), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\left(y - x\right) \cdot \frac{z}{t}
\end{array}
Initial program 97.9%
Taylor expanded in z around inf
div-subN/A
lower-*.f64N/A
lower-/.f64N/A
lower--.f6459.3
Applied rewrites59.3%
Applied rewrites62.0%
Final simplification62.0%
(FPCore (x y z t) :precision binary64 (* z (/ (- y x) t)))
double code(double x, double y, double z, double t) {
return z * ((y - x) / 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 = z * ((y - x) / t)
end function
public static double code(double x, double y, double z, double t) {
return z * ((y - x) / t);
}
def code(x, y, z, t): return z * ((y - x) / t)
function code(x, y, z, t) return Float64(z * Float64(Float64(y - x) / t)) end
function tmp = code(x, y, z, t) tmp = z * ((y - x) / t); end
code[x_, y_, z_, t_] := N[(z * N[(N[(y - x), $MachinePrecision] / t), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
z \cdot \frac{y - x}{t}
\end{array}
Initial program 97.9%
Taylor expanded in z around inf
div-subN/A
lower-*.f64N/A
lower-/.f64N/A
lower--.f6459.3
Applied rewrites59.3%
(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 97.9%
Taylor expanded in x around 0
lower-/.f64N/A
lower-*.f6440.0
Applied rewrites40.0%
Applied rewrites41.4%
Final simplification41.4%
(FPCore (x y z t) :precision binary64 (* z (/ y t)))
double code(double x, double y, double z, double t) {
return z * (y / 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 = z * (y / t)
end function
public static double code(double x, double y, double z, double t) {
return z * (y / t);
}
def code(x, y, z, t): return z * (y / t)
function code(x, y, z, t) return Float64(z * Float64(y / t)) end
function tmp = code(x, y, z, t) tmp = z * (y / t); end
code[x_, y_, z_, t_] := N[(z * N[(y / t), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
z \cdot \frac{y}{t}
\end{array}
Initial program 97.9%
Taylor expanded in x around 0
lower-/.f64N/A
lower-*.f6440.0
Applied rewrites40.0%
Applied rewrites37.9%
Final simplification37.9%
(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 2024238
(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))))