
(FPCore (x y z t a) :precision binary64 (+ (/ (* 60.0 (- x y)) (- z t)) (* a 120.0)))
double code(double x, double y, double z, double t, double a) {
return ((60.0 * (x - y)) / (z - t)) + (a * 120.0);
}
real(8) function code(x, y, z, t, a)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
code = ((60.0d0 * (x - y)) / (z - t)) + (a * 120.0d0)
end function
public static double code(double x, double y, double z, double t, double a) {
return ((60.0 * (x - y)) / (z - t)) + (a * 120.0);
}
def code(x, y, z, t, a): return ((60.0 * (x - y)) / (z - t)) + (a * 120.0)
function code(x, y, z, t, a) return Float64(Float64(Float64(60.0 * Float64(x - y)) / Float64(z - t)) + Float64(a * 120.0)) end
function tmp = code(x, y, z, t, a) tmp = ((60.0 * (x - y)) / (z - t)) + (a * 120.0); end
code[x_, y_, z_, t_, a_] := N[(N[(N[(60.0 * N[(x - y), $MachinePrecision]), $MachinePrecision] / N[(z - t), $MachinePrecision]), $MachinePrecision] + N[(a * 120.0), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{60 \cdot \left(x - y\right)}{z - t} + a \cdot 120
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 20 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (x y z t a) :precision binary64 (+ (/ (* 60.0 (- x y)) (- z t)) (* a 120.0)))
double code(double x, double y, double z, double t, double a) {
return ((60.0 * (x - y)) / (z - t)) + (a * 120.0);
}
real(8) function code(x, y, z, t, a)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
code = ((60.0d0 * (x - y)) / (z - t)) + (a * 120.0d0)
end function
public static double code(double x, double y, double z, double t, double a) {
return ((60.0 * (x - y)) / (z - t)) + (a * 120.0);
}
def code(x, y, z, t, a): return ((60.0 * (x - y)) / (z - t)) + (a * 120.0)
function code(x, y, z, t, a) return Float64(Float64(Float64(60.0 * Float64(x - y)) / Float64(z - t)) + Float64(a * 120.0)) end
function tmp = code(x, y, z, t, a) tmp = ((60.0 * (x - y)) / (z - t)) + (a * 120.0); end
code[x_, y_, z_, t_, a_] := N[(N[(N[(60.0 * N[(x - y), $MachinePrecision]), $MachinePrecision] / N[(z - t), $MachinePrecision]), $MachinePrecision] + N[(a * 120.0), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{60 \cdot \left(x - y\right)}{z - t} + a \cdot 120
\end{array}
(FPCore (x y z t a) :precision binary64 (+ (/ 60.0 (/ (- z t) (- x y))) (* a 120.0)))
double code(double x, double y, double z, double t, double a) {
return (60.0 / ((z - t) / (x - y))) + (a * 120.0);
}
real(8) function code(x, y, z, t, a)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
code = (60.0d0 / ((z - t) / (x - y))) + (a * 120.0d0)
end function
public static double code(double x, double y, double z, double t, double a) {
return (60.0 / ((z - t) / (x - y))) + (a * 120.0);
}
def code(x, y, z, t, a): return (60.0 / ((z - t) / (x - y))) + (a * 120.0)
function code(x, y, z, t, a) return Float64(Float64(60.0 / Float64(Float64(z - t) / Float64(x - y))) + Float64(a * 120.0)) end
function tmp = code(x, y, z, t, a) tmp = (60.0 / ((z - t) / (x - y))) + (a * 120.0); end
code[x_, y_, z_, t_, a_] := N[(N[(60.0 / N[(N[(z - t), $MachinePrecision] / N[(x - y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(a * 120.0), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{60}{\frac{z - t}{x - y}} + a \cdot 120
\end{array}
Initial program 99.1%
lift--.f64N/A
lift-*.f64N/A
remove-double-negN/A
lift--.f64N/A
remove-double-negN/A
lift-*.f64N/A
associate-/l*N/A
clear-numN/A
un-div-invN/A
lower-/.f64N/A
lower-/.f6499.8
Applied rewrites99.8%
(FPCore (x y z t a)
:precision binary64
(let* ((t_1 (/ (* 60.0 (- x y)) (- z t))))
(if (<= t_1 -2e+146)
(/ (* y -60.0) z)
(if (<= t_1 1e+295) (* a 120.0) (/ (* 60.0 y) t)))))
double code(double x, double y, double z, double t, double a) {
double t_1 = (60.0 * (x - y)) / (z - t);
double tmp;
if (t_1 <= -2e+146) {
tmp = (y * -60.0) / z;
} else if (t_1 <= 1e+295) {
tmp = a * 120.0;
} else {
tmp = (60.0 * y) / t;
}
return tmp;
}
real(8) function code(x, y, z, t, a)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8) :: t_1
real(8) :: tmp
t_1 = (60.0d0 * (x - y)) / (z - t)
if (t_1 <= (-2d+146)) then
tmp = (y * (-60.0d0)) / z
else if (t_1 <= 1d+295) then
tmp = a * 120.0d0
else
tmp = (60.0d0 * y) / t
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double t_1 = (60.0 * (x - y)) / (z - t);
double tmp;
if (t_1 <= -2e+146) {
tmp = (y * -60.0) / z;
} else if (t_1 <= 1e+295) {
tmp = a * 120.0;
} else {
tmp = (60.0 * y) / t;
}
return tmp;
}
def code(x, y, z, t, a): t_1 = (60.0 * (x - y)) / (z - t) tmp = 0 if t_1 <= -2e+146: tmp = (y * -60.0) / z elif t_1 <= 1e+295: tmp = a * 120.0 else: tmp = (60.0 * y) / t return tmp
function code(x, y, z, t, a) t_1 = Float64(Float64(60.0 * Float64(x - y)) / Float64(z - t)) tmp = 0.0 if (t_1 <= -2e+146) tmp = Float64(Float64(y * -60.0) / z); elseif (t_1 <= 1e+295) tmp = Float64(a * 120.0); else tmp = Float64(Float64(60.0 * y) / t); end return tmp end
function tmp_2 = code(x, y, z, t, a) t_1 = (60.0 * (x - y)) / (z - t); tmp = 0.0; if (t_1 <= -2e+146) tmp = (y * -60.0) / z; elseif (t_1 <= 1e+295) tmp = a * 120.0; else tmp = (60.0 * y) / t; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := Block[{t$95$1 = N[(N[(60.0 * N[(x - y), $MachinePrecision]), $MachinePrecision] / N[(z - t), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$1, -2e+146], N[(N[(y * -60.0), $MachinePrecision] / z), $MachinePrecision], If[LessEqual[t$95$1, 1e+295], N[(a * 120.0), $MachinePrecision], N[(N[(60.0 * y), $MachinePrecision] / t), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{60 \cdot \left(x - y\right)}{z - t}\\
\mathbf{if}\;t\_1 \leq -2 \cdot 10^{+146}:\\
\;\;\;\;\frac{y \cdot -60}{z}\\
\mathbf{elif}\;t\_1 \leq 10^{+295}:\\
\;\;\;\;a \cdot 120\\
\mathbf{else}:\\
\;\;\;\;\frac{60 \cdot y}{t}\\
\end{array}
\end{array}
if (/.f64 (*.f64 #s(literal 60 binary64) (-.f64 x y)) (-.f64 z t)) < -1.99999999999999987e146Initial program 93.6%
lift--.f64N/A
lift-*.f64N/A
remove-double-negN/A
lift--.f64N/A
remove-double-negN/A
clear-numN/A
lower-/.f64N/A
lower-/.f6493.5
Applied rewrites93.5%
Taylor expanded in y around inf
associate-*r/N/A
lower-/.f64N/A
lower-*.f64N/A
lower--.f6452.6
Applied rewrites52.6%
Taylor expanded in z around inf
associate-*r/N/A
lower-/.f64N/A
lower-*.f6442.9
Applied rewrites42.9%
if -1.99999999999999987e146 < (/.f64 (*.f64 #s(literal 60 binary64) (-.f64 x y)) (-.f64 z t)) < 9.9999999999999998e294Initial program 99.8%
Taylor expanded in z around inf
lower-*.f6455.8
Applied rewrites55.8%
if 9.9999999999999998e294 < (/.f64 (*.f64 #s(literal 60 binary64) (-.f64 x y)) (-.f64 z t)) Initial program 100.0%
Taylor expanded in z around 0
lower-fma.f64N/A
lower-/.f64N/A
lower--.f64N/A
lower-*.f64100.0
Applied rewrites100.0%
Taylor expanded in y around inf
associate-*r/N/A
lower-/.f64N/A
*-commutativeN/A
lower-*.f6483.4
Applied rewrites83.4%
Final simplification55.5%
(FPCore (x y z t a)
:precision binary64
(let* ((t_1 (/ (* 60.0 (- x y)) (- z t))))
(if (<= t_1 -2e+146)
(/ (* y -60.0) z)
(if (<= t_1 1e+295) (* a 120.0) (* 60.0 (/ y t))))))
double code(double x, double y, double z, double t, double a) {
double t_1 = (60.0 * (x - y)) / (z - t);
double tmp;
if (t_1 <= -2e+146) {
tmp = (y * -60.0) / z;
} else if (t_1 <= 1e+295) {
tmp = a * 120.0;
} else {
tmp = 60.0 * (y / t);
}
return tmp;
}
real(8) function code(x, y, z, t, a)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8) :: t_1
real(8) :: tmp
t_1 = (60.0d0 * (x - y)) / (z - t)
if (t_1 <= (-2d+146)) then
tmp = (y * (-60.0d0)) / z
else if (t_1 <= 1d+295) then
tmp = a * 120.0d0
else
tmp = 60.0d0 * (y / t)
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double t_1 = (60.0 * (x - y)) / (z - t);
double tmp;
if (t_1 <= -2e+146) {
tmp = (y * -60.0) / z;
} else if (t_1 <= 1e+295) {
tmp = a * 120.0;
} else {
tmp = 60.0 * (y / t);
}
return tmp;
}
def code(x, y, z, t, a): t_1 = (60.0 * (x - y)) / (z - t) tmp = 0 if t_1 <= -2e+146: tmp = (y * -60.0) / z elif t_1 <= 1e+295: tmp = a * 120.0 else: tmp = 60.0 * (y / t) return tmp
function code(x, y, z, t, a) t_1 = Float64(Float64(60.0 * Float64(x - y)) / Float64(z - t)) tmp = 0.0 if (t_1 <= -2e+146) tmp = Float64(Float64(y * -60.0) / z); elseif (t_1 <= 1e+295) tmp = Float64(a * 120.0); else tmp = Float64(60.0 * Float64(y / t)); end return tmp end
function tmp_2 = code(x, y, z, t, a) t_1 = (60.0 * (x - y)) / (z - t); tmp = 0.0; if (t_1 <= -2e+146) tmp = (y * -60.0) / z; elseif (t_1 <= 1e+295) tmp = a * 120.0; else tmp = 60.0 * (y / t); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := Block[{t$95$1 = N[(N[(60.0 * N[(x - y), $MachinePrecision]), $MachinePrecision] / N[(z - t), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$1, -2e+146], N[(N[(y * -60.0), $MachinePrecision] / z), $MachinePrecision], If[LessEqual[t$95$1, 1e+295], N[(a * 120.0), $MachinePrecision], N[(60.0 * N[(y / t), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{60 \cdot \left(x - y\right)}{z - t}\\
\mathbf{if}\;t\_1 \leq -2 \cdot 10^{+146}:\\
\;\;\;\;\frac{y \cdot -60}{z}\\
\mathbf{elif}\;t\_1 \leq 10^{+295}:\\
\;\;\;\;a \cdot 120\\
\mathbf{else}:\\
\;\;\;\;60 \cdot \frac{y}{t}\\
\end{array}
\end{array}
if (/.f64 (*.f64 #s(literal 60 binary64) (-.f64 x y)) (-.f64 z t)) < -1.99999999999999987e146Initial program 93.6%
lift--.f64N/A
lift-*.f64N/A
remove-double-negN/A
lift--.f64N/A
remove-double-negN/A
clear-numN/A
lower-/.f64N/A
lower-/.f6493.5
Applied rewrites93.5%
Taylor expanded in y around inf
associate-*r/N/A
lower-/.f64N/A
lower-*.f64N/A
lower--.f6452.6
Applied rewrites52.6%
Taylor expanded in z around inf
associate-*r/N/A
lower-/.f64N/A
lower-*.f6442.9
Applied rewrites42.9%
if -1.99999999999999987e146 < (/.f64 (*.f64 #s(literal 60 binary64) (-.f64 x y)) (-.f64 z t)) < 9.9999999999999998e294Initial program 99.8%
Taylor expanded in z around inf
lower-*.f6455.8
Applied rewrites55.8%
if 9.9999999999999998e294 < (/.f64 (*.f64 #s(literal 60 binary64) (-.f64 x y)) (-.f64 z t)) Initial program 100.0%
Taylor expanded in z around 0
lower-fma.f64N/A
lower-/.f64N/A
lower--.f64N/A
lower-*.f64100.0
Applied rewrites100.0%
Taylor expanded in y around inf
associate-*r/N/A
lower-/.f64N/A
*-commutativeN/A
lower-*.f6483.4
Applied rewrites83.4%
associate-*l/N/A
lift-/.f64N/A
lower-*.f6483.4
Applied rewrites83.4%
Final simplification55.5%
(FPCore (x y z t a)
:precision binary64
(let* ((t_1 (/ (* 60.0 (- x y)) (- z t))))
(if (<= t_1 -1.5e+126)
(* x (/ 60.0 z))
(if (<= t_1 1e+295) (* a 120.0) (* 60.0 (/ y t))))))
double code(double x, double y, double z, double t, double a) {
double t_1 = (60.0 * (x - y)) / (z - t);
double tmp;
if (t_1 <= -1.5e+126) {
tmp = x * (60.0 / z);
} else if (t_1 <= 1e+295) {
tmp = a * 120.0;
} else {
tmp = 60.0 * (y / t);
}
return tmp;
}
real(8) function code(x, y, z, t, a)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8) :: t_1
real(8) :: tmp
t_1 = (60.0d0 * (x - y)) / (z - t)
if (t_1 <= (-1.5d+126)) then
tmp = x * (60.0d0 / z)
else if (t_1 <= 1d+295) then
tmp = a * 120.0d0
else
tmp = 60.0d0 * (y / t)
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double t_1 = (60.0 * (x - y)) / (z - t);
double tmp;
if (t_1 <= -1.5e+126) {
tmp = x * (60.0 / z);
} else if (t_1 <= 1e+295) {
tmp = a * 120.0;
} else {
tmp = 60.0 * (y / t);
}
return tmp;
}
def code(x, y, z, t, a): t_1 = (60.0 * (x - y)) / (z - t) tmp = 0 if t_1 <= -1.5e+126: tmp = x * (60.0 / z) elif t_1 <= 1e+295: tmp = a * 120.0 else: tmp = 60.0 * (y / t) return tmp
function code(x, y, z, t, a) t_1 = Float64(Float64(60.0 * Float64(x - y)) / Float64(z - t)) tmp = 0.0 if (t_1 <= -1.5e+126) tmp = Float64(x * Float64(60.0 / z)); elseif (t_1 <= 1e+295) tmp = Float64(a * 120.0); else tmp = Float64(60.0 * Float64(y / t)); end return tmp end
function tmp_2 = code(x, y, z, t, a) t_1 = (60.0 * (x - y)) / (z - t); tmp = 0.0; if (t_1 <= -1.5e+126) tmp = x * (60.0 / z); elseif (t_1 <= 1e+295) tmp = a * 120.0; else tmp = 60.0 * (y / t); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := Block[{t$95$1 = N[(N[(60.0 * N[(x - y), $MachinePrecision]), $MachinePrecision] / N[(z - t), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$1, -1.5e+126], N[(x * N[(60.0 / z), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$1, 1e+295], N[(a * 120.0), $MachinePrecision], N[(60.0 * N[(y / t), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{60 \cdot \left(x - y\right)}{z - t}\\
\mathbf{if}\;t\_1 \leq -1.5 \cdot 10^{+126}:\\
\;\;\;\;x \cdot \frac{60}{z}\\
\mathbf{elif}\;t\_1 \leq 10^{+295}:\\
\;\;\;\;a \cdot 120\\
\mathbf{else}:\\
\;\;\;\;60 \cdot \frac{y}{t}\\
\end{array}
\end{array}
if (/.f64 (*.f64 #s(literal 60 binary64) (-.f64 x y)) (-.f64 z t)) < -1.5000000000000001e126Initial program 94.8%
Taylor expanded in x around inf
associate-*r/N/A
lower-/.f64N/A
lower-*.f64N/A
lower--.f6446.2
Applied rewrites46.2%
*-commutativeN/A
lift--.f64N/A
associate-/l*N/A
lower-*.f64N/A
lower-/.f6446.1
Applied rewrites46.1%
Taylor expanded in z around inf
associate-*r/N/A
*-commutativeN/A
associate-/l*N/A
metadata-evalN/A
associate-*r/N/A
lower-*.f64N/A
associate-*r/N/A
metadata-evalN/A
lower-/.f6433.3
Applied rewrites33.3%
if -1.5000000000000001e126 < (/.f64 (*.f64 #s(literal 60 binary64) (-.f64 x y)) (-.f64 z t)) < 9.9999999999999998e294Initial program 99.8%
Taylor expanded in z around inf
lower-*.f6457.1
Applied rewrites57.1%
if 9.9999999999999998e294 < (/.f64 (*.f64 #s(literal 60 binary64) (-.f64 x y)) (-.f64 z t)) Initial program 100.0%
Taylor expanded in z around 0
lower-fma.f64N/A
lower-/.f64N/A
lower--.f64N/A
lower-*.f64100.0
Applied rewrites100.0%
Taylor expanded in y around inf
associate-*r/N/A
lower-/.f64N/A
*-commutativeN/A
lower-*.f6483.4
Applied rewrites83.4%
associate-*l/N/A
lift-/.f64N/A
lower-*.f6483.4
Applied rewrites83.4%
Final simplification54.8%
(FPCore (x y z t a)
:precision binary64
(if (<= (* a 120.0) -5e-92)
(* a 120.0)
(if (<= (* a 120.0) 2e-268)
(/ (* (- x y) -60.0) t)
(if (<= (* a 120.0) 400000000.0) (/ (* 60.0 (- x y)) z) (* a 120.0)))))
double code(double x, double y, double z, double t, double a) {
double tmp;
if ((a * 120.0) <= -5e-92) {
tmp = a * 120.0;
} else if ((a * 120.0) <= 2e-268) {
tmp = ((x - y) * -60.0) / t;
} else if ((a * 120.0) <= 400000000.0) {
tmp = (60.0 * (x - y)) / z;
} else {
tmp = a * 120.0;
}
return tmp;
}
real(8) function code(x, y, z, t, a)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8) :: tmp
if ((a * 120.0d0) <= (-5d-92)) then
tmp = a * 120.0d0
else if ((a * 120.0d0) <= 2d-268) then
tmp = ((x - y) * (-60.0d0)) / t
else if ((a * 120.0d0) <= 400000000.0d0) then
tmp = (60.0d0 * (x - y)) / z
else
tmp = a * 120.0d0
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double tmp;
if ((a * 120.0) <= -5e-92) {
tmp = a * 120.0;
} else if ((a * 120.0) <= 2e-268) {
tmp = ((x - y) * -60.0) / t;
} else if ((a * 120.0) <= 400000000.0) {
tmp = (60.0 * (x - y)) / z;
} else {
tmp = a * 120.0;
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if (a * 120.0) <= -5e-92: tmp = a * 120.0 elif (a * 120.0) <= 2e-268: tmp = ((x - y) * -60.0) / t elif (a * 120.0) <= 400000000.0: tmp = (60.0 * (x - y)) / z else: tmp = a * 120.0 return tmp
function code(x, y, z, t, a) tmp = 0.0 if (Float64(a * 120.0) <= -5e-92) tmp = Float64(a * 120.0); elseif (Float64(a * 120.0) <= 2e-268) tmp = Float64(Float64(Float64(x - y) * -60.0) / t); elseif (Float64(a * 120.0) <= 400000000.0) tmp = Float64(Float64(60.0 * Float64(x - y)) / z); else tmp = Float64(a * 120.0); end return tmp end
function tmp_2 = code(x, y, z, t, a) tmp = 0.0; if ((a * 120.0) <= -5e-92) tmp = a * 120.0; elseif ((a * 120.0) <= 2e-268) tmp = ((x - y) * -60.0) / t; elseif ((a * 120.0) <= 400000000.0) tmp = (60.0 * (x - y)) / z; else tmp = a * 120.0; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[LessEqual[N[(a * 120.0), $MachinePrecision], -5e-92], N[(a * 120.0), $MachinePrecision], If[LessEqual[N[(a * 120.0), $MachinePrecision], 2e-268], N[(N[(N[(x - y), $MachinePrecision] * -60.0), $MachinePrecision] / t), $MachinePrecision], If[LessEqual[N[(a * 120.0), $MachinePrecision], 400000000.0], N[(N[(60.0 * N[(x - y), $MachinePrecision]), $MachinePrecision] / z), $MachinePrecision], N[(a * 120.0), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;a \cdot 120 \leq -5 \cdot 10^{-92}:\\
\;\;\;\;a \cdot 120\\
\mathbf{elif}\;a \cdot 120 \leq 2 \cdot 10^{-268}:\\
\;\;\;\;\frac{\left(x - y\right) \cdot -60}{t}\\
\mathbf{elif}\;a \cdot 120 \leq 400000000:\\
\;\;\;\;\frac{60 \cdot \left(x - y\right)}{z}\\
\mathbf{else}:\\
\;\;\;\;a \cdot 120\\
\end{array}
\end{array}
if (*.f64 a #s(literal 120 binary64)) < -5.00000000000000011e-92 or 4e8 < (*.f64 a #s(literal 120 binary64)) Initial program 98.5%
Taylor expanded in z around inf
lower-*.f6474.8
Applied rewrites74.8%
if -5.00000000000000011e-92 < (*.f64 a #s(literal 120 binary64)) < 1.99999999999999992e-268Initial program 99.7%
Taylor expanded in z around 0
lower-fma.f64N/A
lower-/.f64N/A
lower--.f64N/A
lower-*.f6465.1
Applied rewrites65.1%
Taylor expanded in t around 0
associate-*r/N/A
lower-/.f64N/A
lower-*.f64N/A
lower--.f6454.3
Applied rewrites54.3%
if 1.99999999999999992e-268 < (*.f64 a #s(literal 120 binary64)) < 4e8Initial program 99.8%
lift--.f64N/A
lift-*.f64N/A
remove-double-negN/A
lift--.f64N/A
remove-double-negN/A
lift-*.f64N/A
associate-/l*N/A
clear-numN/A
un-div-invN/A
lower-/.f64N/A
lower-/.f6499.7
Applied rewrites99.7%
Taylor expanded in a around 0
associate-*r/N/A
lower-/.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower--.f64N/A
lower--.f6484.1
Applied rewrites84.1%
Taylor expanded in z around inf
associate-*r/N/A
lower-/.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower--.f6453.2
Applied rewrites53.2%
Final simplification64.8%
(FPCore (x y z t a)
:precision binary64
(if (<= (* a 120.0) -5e-80)
(* a 120.0)
(if (<= (* a 120.0) 1e-205)
(/ (* y -60.0) (- z t))
(if (<= (* a 120.0) 5000.0)
(/ x (* (- z t) 0.016666666666666666))
(* a 120.0)))))
double code(double x, double y, double z, double t, double a) {
double tmp;
if ((a * 120.0) <= -5e-80) {
tmp = a * 120.0;
} else if ((a * 120.0) <= 1e-205) {
tmp = (y * -60.0) / (z - t);
} else if ((a * 120.0) <= 5000.0) {
tmp = x / ((z - t) * 0.016666666666666666);
} else {
tmp = a * 120.0;
}
return tmp;
}
real(8) function code(x, y, z, t, a)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8) :: tmp
if ((a * 120.0d0) <= (-5d-80)) then
tmp = a * 120.0d0
else if ((a * 120.0d0) <= 1d-205) then
tmp = (y * (-60.0d0)) / (z - t)
else if ((a * 120.0d0) <= 5000.0d0) then
tmp = x / ((z - t) * 0.016666666666666666d0)
else
tmp = a * 120.0d0
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double tmp;
if ((a * 120.0) <= -5e-80) {
tmp = a * 120.0;
} else if ((a * 120.0) <= 1e-205) {
tmp = (y * -60.0) / (z - t);
} else if ((a * 120.0) <= 5000.0) {
tmp = x / ((z - t) * 0.016666666666666666);
} else {
tmp = a * 120.0;
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if (a * 120.0) <= -5e-80: tmp = a * 120.0 elif (a * 120.0) <= 1e-205: tmp = (y * -60.0) / (z - t) elif (a * 120.0) <= 5000.0: tmp = x / ((z - t) * 0.016666666666666666) else: tmp = a * 120.0 return tmp
function code(x, y, z, t, a) tmp = 0.0 if (Float64(a * 120.0) <= -5e-80) tmp = Float64(a * 120.0); elseif (Float64(a * 120.0) <= 1e-205) tmp = Float64(Float64(y * -60.0) / Float64(z - t)); elseif (Float64(a * 120.0) <= 5000.0) tmp = Float64(x / Float64(Float64(z - t) * 0.016666666666666666)); else tmp = Float64(a * 120.0); end return tmp end
function tmp_2 = code(x, y, z, t, a) tmp = 0.0; if ((a * 120.0) <= -5e-80) tmp = a * 120.0; elseif ((a * 120.0) <= 1e-205) tmp = (y * -60.0) / (z - t); elseif ((a * 120.0) <= 5000.0) tmp = x / ((z - t) * 0.016666666666666666); else tmp = a * 120.0; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[LessEqual[N[(a * 120.0), $MachinePrecision], -5e-80], N[(a * 120.0), $MachinePrecision], If[LessEqual[N[(a * 120.0), $MachinePrecision], 1e-205], N[(N[(y * -60.0), $MachinePrecision] / N[(z - t), $MachinePrecision]), $MachinePrecision], If[LessEqual[N[(a * 120.0), $MachinePrecision], 5000.0], N[(x / N[(N[(z - t), $MachinePrecision] * 0.016666666666666666), $MachinePrecision]), $MachinePrecision], N[(a * 120.0), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;a \cdot 120 \leq -5 \cdot 10^{-80}:\\
\;\;\;\;a \cdot 120\\
\mathbf{elif}\;a \cdot 120 \leq 10^{-205}:\\
\;\;\;\;\frac{y \cdot -60}{z - t}\\
\mathbf{elif}\;a \cdot 120 \leq 5000:\\
\;\;\;\;\frac{x}{\left(z - t\right) \cdot 0.016666666666666666}\\
\mathbf{else}:\\
\;\;\;\;a \cdot 120\\
\end{array}
\end{array}
if (*.f64 a #s(literal 120 binary64)) < -5e-80 or 5e3 < (*.f64 a #s(literal 120 binary64)) Initial program 98.5%
Taylor expanded in z around inf
lower-*.f6475.4
Applied rewrites75.4%
if -5e-80 < (*.f64 a #s(literal 120 binary64)) < 1e-205Initial program 99.7%
lift--.f64N/A
lift-*.f64N/A
remove-double-negN/A
lift--.f64N/A
remove-double-negN/A
clear-numN/A
lower-/.f64N/A
lower-/.f6499.5
Applied rewrites99.5%
Taylor expanded in y around inf
associate-*r/N/A
lower-/.f64N/A
lower-*.f64N/A
lower--.f6454.1
Applied rewrites54.1%
if 1e-205 < (*.f64 a #s(literal 120 binary64)) < 5e3Initial program 99.8%
Taylor expanded in x around inf
associate-*r/N/A
lower-/.f64N/A
lower-*.f64N/A
lower--.f6449.1
Applied rewrites49.1%
*-commutativeN/A
lift--.f64N/A
associate-/l*N/A
clear-numN/A
un-div-invN/A
lower-/.f64N/A
div-invN/A
lower-*.f64N/A
metadata-eval49.1
Applied rewrites49.1%
Final simplification64.2%
(FPCore (x y z t a)
:precision binary64
(if (<= (* a 120.0) -5e-80)
(* a 120.0)
(if (<= (* a 120.0) 1e-205)
(* -60.0 (/ y (- z t)))
(if (<= (* a 120.0) 5000.0)
(/ x (* (- z t) 0.016666666666666666))
(* a 120.0)))))
double code(double x, double y, double z, double t, double a) {
double tmp;
if ((a * 120.0) <= -5e-80) {
tmp = a * 120.0;
} else if ((a * 120.0) <= 1e-205) {
tmp = -60.0 * (y / (z - t));
} else if ((a * 120.0) <= 5000.0) {
tmp = x / ((z - t) * 0.016666666666666666);
} else {
tmp = a * 120.0;
}
return tmp;
}
real(8) function code(x, y, z, t, a)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8) :: tmp
if ((a * 120.0d0) <= (-5d-80)) then
tmp = a * 120.0d0
else if ((a * 120.0d0) <= 1d-205) then
tmp = (-60.0d0) * (y / (z - t))
else if ((a * 120.0d0) <= 5000.0d0) then
tmp = x / ((z - t) * 0.016666666666666666d0)
else
tmp = a * 120.0d0
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double tmp;
if ((a * 120.0) <= -5e-80) {
tmp = a * 120.0;
} else if ((a * 120.0) <= 1e-205) {
tmp = -60.0 * (y / (z - t));
} else if ((a * 120.0) <= 5000.0) {
tmp = x / ((z - t) * 0.016666666666666666);
} else {
tmp = a * 120.0;
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if (a * 120.0) <= -5e-80: tmp = a * 120.0 elif (a * 120.0) <= 1e-205: tmp = -60.0 * (y / (z - t)) elif (a * 120.0) <= 5000.0: tmp = x / ((z - t) * 0.016666666666666666) else: tmp = a * 120.0 return tmp
function code(x, y, z, t, a) tmp = 0.0 if (Float64(a * 120.0) <= -5e-80) tmp = Float64(a * 120.0); elseif (Float64(a * 120.0) <= 1e-205) tmp = Float64(-60.0 * Float64(y / Float64(z - t))); elseif (Float64(a * 120.0) <= 5000.0) tmp = Float64(x / Float64(Float64(z - t) * 0.016666666666666666)); else tmp = Float64(a * 120.0); end return tmp end
function tmp_2 = code(x, y, z, t, a) tmp = 0.0; if ((a * 120.0) <= -5e-80) tmp = a * 120.0; elseif ((a * 120.0) <= 1e-205) tmp = -60.0 * (y / (z - t)); elseif ((a * 120.0) <= 5000.0) tmp = x / ((z - t) * 0.016666666666666666); else tmp = a * 120.0; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[LessEqual[N[(a * 120.0), $MachinePrecision], -5e-80], N[(a * 120.0), $MachinePrecision], If[LessEqual[N[(a * 120.0), $MachinePrecision], 1e-205], N[(-60.0 * N[(y / N[(z - t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[N[(a * 120.0), $MachinePrecision], 5000.0], N[(x / N[(N[(z - t), $MachinePrecision] * 0.016666666666666666), $MachinePrecision]), $MachinePrecision], N[(a * 120.0), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;a \cdot 120 \leq -5 \cdot 10^{-80}:\\
\;\;\;\;a \cdot 120\\
\mathbf{elif}\;a \cdot 120 \leq 10^{-205}:\\
\;\;\;\;-60 \cdot \frac{y}{z - t}\\
\mathbf{elif}\;a \cdot 120 \leq 5000:\\
\;\;\;\;\frac{x}{\left(z - t\right) \cdot 0.016666666666666666}\\
\mathbf{else}:\\
\;\;\;\;a \cdot 120\\
\end{array}
\end{array}
if (*.f64 a #s(literal 120 binary64)) < -5e-80 or 5e3 < (*.f64 a #s(literal 120 binary64)) Initial program 98.5%
Taylor expanded in z around inf
lower-*.f6475.4
Applied rewrites75.4%
if -5e-80 < (*.f64 a #s(literal 120 binary64)) < 1e-205Initial program 99.7%
Taylor expanded in y around inf
lower-*.f64N/A
lower-/.f64N/A
lower--.f6454.1
Applied rewrites54.1%
if 1e-205 < (*.f64 a #s(literal 120 binary64)) < 5e3Initial program 99.8%
Taylor expanded in x around inf
associate-*r/N/A
lower-/.f64N/A
lower-*.f64N/A
lower--.f6449.1
Applied rewrites49.1%
*-commutativeN/A
lift--.f64N/A
associate-/l*N/A
clear-numN/A
un-div-invN/A
lower-/.f64N/A
div-invN/A
lower-*.f64N/A
metadata-eval49.1
Applied rewrites49.1%
Final simplification64.2%
(FPCore (x y z t a)
:precision binary64
(if (<= (* a 120.0) -5e-80)
(* a 120.0)
(if (<= (* a 120.0) -2e-301)
(* -60.0 (/ y (- z t)))
(if (<= (* a 120.0) 5000.0) (* 60.0 (/ x (- z t))) (* a 120.0)))))
double code(double x, double y, double z, double t, double a) {
double tmp;
if ((a * 120.0) <= -5e-80) {
tmp = a * 120.0;
} else if ((a * 120.0) <= -2e-301) {
tmp = -60.0 * (y / (z - t));
} else if ((a * 120.0) <= 5000.0) {
tmp = 60.0 * (x / (z - t));
} else {
tmp = a * 120.0;
}
return tmp;
}
real(8) function code(x, y, z, t, a)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8) :: tmp
if ((a * 120.0d0) <= (-5d-80)) then
tmp = a * 120.0d0
else if ((a * 120.0d0) <= (-2d-301)) then
tmp = (-60.0d0) * (y / (z - t))
else if ((a * 120.0d0) <= 5000.0d0) then
tmp = 60.0d0 * (x / (z - t))
else
tmp = a * 120.0d0
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double tmp;
if ((a * 120.0) <= -5e-80) {
tmp = a * 120.0;
} else if ((a * 120.0) <= -2e-301) {
tmp = -60.0 * (y / (z - t));
} else if ((a * 120.0) <= 5000.0) {
tmp = 60.0 * (x / (z - t));
} else {
tmp = a * 120.0;
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if (a * 120.0) <= -5e-80: tmp = a * 120.0 elif (a * 120.0) <= -2e-301: tmp = -60.0 * (y / (z - t)) elif (a * 120.0) <= 5000.0: tmp = 60.0 * (x / (z - t)) else: tmp = a * 120.0 return tmp
function code(x, y, z, t, a) tmp = 0.0 if (Float64(a * 120.0) <= -5e-80) tmp = Float64(a * 120.0); elseif (Float64(a * 120.0) <= -2e-301) tmp = Float64(-60.0 * Float64(y / Float64(z - t))); elseif (Float64(a * 120.0) <= 5000.0) tmp = Float64(60.0 * Float64(x / Float64(z - t))); else tmp = Float64(a * 120.0); end return tmp end
function tmp_2 = code(x, y, z, t, a) tmp = 0.0; if ((a * 120.0) <= -5e-80) tmp = a * 120.0; elseif ((a * 120.0) <= -2e-301) tmp = -60.0 * (y / (z - t)); elseif ((a * 120.0) <= 5000.0) tmp = 60.0 * (x / (z - t)); else tmp = a * 120.0; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[LessEqual[N[(a * 120.0), $MachinePrecision], -5e-80], N[(a * 120.0), $MachinePrecision], If[LessEqual[N[(a * 120.0), $MachinePrecision], -2e-301], N[(-60.0 * N[(y / N[(z - t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[N[(a * 120.0), $MachinePrecision], 5000.0], N[(60.0 * N[(x / N[(z - t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(a * 120.0), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;a \cdot 120 \leq -5 \cdot 10^{-80}:\\
\;\;\;\;a \cdot 120\\
\mathbf{elif}\;a \cdot 120 \leq -2 \cdot 10^{-301}:\\
\;\;\;\;-60 \cdot \frac{y}{z - t}\\
\mathbf{elif}\;a \cdot 120 \leq 5000:\\
\;\;\;\;60 \cdot \frac{x}{z - t}\\
\mathbf{else}:\\
\;\;\;\;a \cdot 120\\
\end{array}
\end{array}
if (*.f64 a #s(literal 120 binary64)) < -5e-80 or 5e3 < (*.f64 a #s(literal 120 binary64)) Initial program 98.5%
Taylor expanded in z around inf
lower-*.f6475.4
Applied rewrites75.4%
if -5e-80 < (*.f64 a #s(literal 120 binary64)) < -2.00000000000000013e-301Initial program 99.7%
Taylor expanded in y around inf
lower-*.f64N/A
lower-/.f64N/A
lower--.f6457.3
Applied rewrites57.3%
if -2.00000000000000013e-301 < (*.f64 a #s(literal 120 binary64)) < 5e3Initial program 99.8%
Taylor expanded in x around inf
associate-*r/N/A
lower-/.f64N/A
lower-*.f64N/A
lower--.f6448.5
Applied rewrites48.5%
lift--.f64N/A
associate-/l*N/A
*-commutativeN/A
lower-*.f64N/A
lower-/.f6448.4
Applied rewrites48.4%
Final simplification64.2%
(FPCore (x y z t a)
:precision binary64
(if (<= (* a 120.0) -5e-80)
(* a 120.0)
(if (<= (* a 120.0) 1e-205)
(* -60.0 (/ y (- z t)))
(if (<= (* a 120.0) 5000.0) (* x (/ 60.0 (- z t))) (* a 120.0)))))
double code(double x, double y, double z, double t, double a) {
double tmp;
if ((a * 120.0) <= -5e-80) {
tmp = a * 120.0;
} else if ((a * 120.0) <= 1e-205) {
tmp = -60.0 * (y / (z - t));
} else if ((a * 120.0) <= 5000.0) {
tmp = x * (60.0 / (z - t));
} else {
tmp = a * 120.0;
}
return tmp;
}
real(8) function code(x, y, z, t, a)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8) :: tmp
if ((a * 120.0d0) <= (-5d-80)) then
tmp = a * 120.0d0
else if ((a * 120.0d0) <= 1d-205) then
tmp = (-60.0d0) * (y / (z - t))
else if ((a * 120.0d0) <= 5000.0d0) then
tmp = x * (60.0d0 / (z - t))
else
tmp = a * 120.0d0
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double tmp;
if ((a * 120.0) <= -5e-80) {
tmp = a * 120.0;
} else if ((a * 120.0) <= 1e-205) {
tmp = -60.0 * (y / (z - t));
} else if ((a * 120.0) <= 5000.0) {
tmp = x * (60.0 / (z - t));
} else {
tmp = a * 120.0;
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if (a * 120.0) <= -5e-80: tmp = a * 120.0 elif (a * 120.0) <= 1e-205: tmp = -60.0 * (y / (z - t)) elif (a * 120.0) <= 5000.0: tmp = x * (60.0 / (z - t)) else: tmp = a * 120.0 return tmp
function code(x, y, z, t, a) tmp = 0.0 if (Float64(a * 120.0) <= -5e-80) tmp = Float64(a * 120.0); elseif (Float64(a * 120.0) <= 1e-205) tmp = Float64(-60.0 * Float64(y / Float64(z - t))); elseif (Float64(a * 120.0) <= 5000.0) tmp = Float64(x * Float64(60.0 / Float64(z - t))); else tmp = Float64(a * 120.0); end return tmp end
function tmp_2 = code(x, y, z, t, a) tmp = 0.0; if ((a * 120.0) <= -5e-80) tmp = a * 120.0; elseif ((a * 120.0) <= 1e-205) tmp = -60.0 * (y / (z - t)); elseif ((a * 120.0) <= 5000.0) tmp = x * (60.0 / (z - t)); else tmp = a * 120.0; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[LessEqual[N[(a * 120.0), $MachinePrecision], -5e-80], N[(a * 120.0), $MachinePrecision], If[LessEqual[N[(a * 120.0), $MachinePrecision], 1e-205], N[(-60.0 * N[(y / N[(z - t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[N[(a * 120.0), $MachinePrecision], 5000.0], N[(x * N[(60.0 / N[(z - t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(a * 120.0), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;a \cdot 120 \leq -5 \cdot 10^{-80}:\\
\;\;\;\;a \cdot 120\\
\mathbf{elif}\;a \cdot 120 \leq 10^{-205}:\\
\;\;\;\;-60 \cdot \frac{y}{z - t}\\
\mathbf{elif}\;a \cdot 120 \leq 5000:\\
\;\;\;\;x \cdot \frac{60}{z - t}\\
\mathbf{else}:\\
\;\;\;\;a \cdot 120\\
\end{array}
\end{array}
if (*.f64 a #s(literal 120 binary64)) < -5e-80 or 5e3 < (*.f64 a #s(literal 120 binary64)) Initial program 98.5%
Taylor expanded in z around inf
lower-*.f6475.4
Applied rewrites75.4%
if -5e-80 < (*.f64 a #s(literal 120 binary64)) < 1e-205Initial program 99.7%
Taylor expanded in y around inf
lower-*.f64N/A
lower-/.f64N/A
lower--.f6454.1
Applied rewrites54.1%
if 1e-205 < (*.f64 a #s(literal 120 binary64)) < 5e3Initial program 99.8%
Taylor expanded in x around inf
associate-*r/N/A
lower-/.f64N/A
lower-*.f64N/A
lower--.f6449.1
Applied rewrites49.1%
*-commutativeN/A
lift--.f64N/A
associate-/l*N/A
lower-*.f64N/A
lower-/.f6449.0
Applied rewrites49.0%
Final simplification64.2%
(FPCore (x y z t a)
:precision binary64
(let* ((t_1 (fma -60.0 (/ (- x y) t) (* a 120.0))))
(if (<= t -4.6e-26)
t_1
(if (<= t -5e-113)
(/ (* 60.0 (- x y)) (- z t))
(if (<= t 1.05e-100)
(fma 60.0 (/ (- x y) z) (* a 120.0))
(if (<= t 7.5e+57) (+ (* a 120.0) (/ (* 60.0 x) (- z t))) t_1))))))
double code(double x, double y, double z, double t, double a) {
double t_1 = fma(-60.0, ((x - y) / t), (a * 120.0));
double tmp;
if (t <= -4.6e-26) {
tmp = t_1;
} else if (t <= -5e-113) {
tmp = (60.0 * (x - y)) / (z - t);
} else if (t <= 1.05e-100) {
tmp = fma(60.0, ((x - y) / z), (a * 120.0));
} else if (t <= 7.5e+57) {
tmp = (a * 120.0) + ((60.0 * x) / (z - t));
} else {
tmp = t_1;
}
return tmp;
}
function code(x, y, z, t, a) t_1 = fma(-60.0, Float64(Float64(x - y) / t), Float64(a * 120.0)) tmp = 0.0 if (t <= -4.6e-26) tmp = t_1; elseif (t <= -5e-113) tmp = Float64(Float64(60.0 * Float64(x - y)) / Float64(z - t)); elseif (t <= 1.05e-100) tmp = fma(60.0, Float64(Float64(x - y) / z), Float64(a * 120.0)); elseif (t <= 7.5e+57) tmp = Float64(Float64(a * 120.0) + Float64(Float64(60.0 * x) / Float64(z - t))); else tmp = t_1; end return tmp end
code[x_, y_, z_, t_, a_] := Block[{t$95$1 = N[(-60.0 * N[(N[(x - y), $MachinePrecision] / t), $MachinePrecision] + N[(a * 120.0), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t, -4.6e-26], t$95$1, If[LessEqual[t, -5e-113], N[(N[(60.0 * N[(x - y), $MachinePrecision]), $MachinePrecision] / N[(z - t), $MachinePrecision]), $MachinePrecision], If[LessEqual[t, 1.05e-100], N[(60.0 * N[(N[(x - y), $MachinePrecision] / z), $MachinePrecision] + N[(a * 120.0), $MachinePrecision]), $MachinePrecision], If[LessEqual[t, 7.5e+57], N[(N[(a * 120.0), $MachinePrecision] + N[(N[(60.0 * x), $MachinePrecision] / N[(z - t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$1]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \mathsf{fma}\left(-60, \frac{x - y}{t}, a \cdot 120\right)\\
\mathbf{if}\;t \leq -4.6 \cdot 10^{-26}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;t \leq -5 \cdot 10^{-113}:\\
\;\;\;\;\frac{60 \cdot \left(x - y\right)}{z - t}\\
\mathbf{elif}\;t \leq 1.05 \cdot 10^{-100}:\\
\;\;\;\;\mathsf{fma}\left(60, \frac{x - y}{z}, a \cdot 120\right)\\
\mathbf{elif}\;t \leq 7.5 \cdot 10^{+57}:\\
\;\;\;\;a \cdot 120 + \frac{60 \cdot x}{z - t}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if t < -4.60000000000000018e-26 or 7.5000000000000006e57 < t Initial program 99.8%
Taylor expanded in z around 0
lower-fma.f64N/A
lower-/.f64N/A
lower--.f64N/A
lower-*.f6489.6
Applied rewrites89.6%
if -4.60000000000000018e-26 < t < -4.9999999999999997e-113Initial program 95.5%
Taylor expanded in a around 0
associate-*r/N/A
lower-/.f64N/A
lower-*.f64N/A
lower--.f64N/A
lower--.f6481.7
Applied rewrites81.7%
if -4.9999999999999997e-113 < t < 1.05000000000000005e-100Initial program 98.9%
Taylor expanded in z around inf
lower-fma.f64N/A
lower-/.f64N/A
lower--.f64N/A
lower-*.f6490.8
Applied rewrites90.8%
if 1.05000000000000005e-100 < t < 7.5000000000000006e57Initial program 99.9%
Taylor expanded in x around inf
associate-*r/N/A
lower-/.f64N/A
lower-*.f64N/A
lower--.f6489.5
Applied rewrites89.5%
Final simplification89.4%
(FPCore (x y z t a)
:precision binary64
(let* ((t_1 (fma -60.0 (/ (- x y) t) (* a 120.0))))
(if (<= t -4.6e-26)
t_1
(if (<= t -5e-113)
(/ (* 60.0 (- x y)) (- z t))
(if (<= t 1.05e-100)
(fma 60.0 (/ (- x y) z) (* a 120.0))
(if (<= t 7.5e+57)
(fma a 120.0 (/ x (* (- z t) 0.016666666666666666)))
t_1))))))
double code(double x, double y, double z, double t, double a) {
double t_1 = fma(-60.0, ((x - y) / t), (a * 120.0));
double tmp;
if (t <= -4.6e-26) {
tmp = t_1;
} else if (t <= -5e-113) {
tmp = (60.0 * (x - y)) / (z - t);
} else if (t <= 1.05e-100) {
tmp = fma(60.0, ((x - y) / z), (a * 120.0));
} else if (t <= 7.5e+57) {
tmp = fma(a, 120.0, (x / ((z - t) * 0.016666666666666666)));
} else {
tmp = t_1;
}
return tmp;
}
function code(x, y, z, t, a) t_1 = fma(-60.0, Float64(Float64(x - y) / t), Float64(a * 120.0)) tmp = 0.0 if (t <= -4.6e-26) tmp = t_1; elseif (t <= -5e-113) tmp = Float64(Float64(60.0 * Float64(x - y)) / Float64(z - t)); elseif (t <= 1.05e-100) tmp = fma(60.0, Float64(Float64(x - y) / z), Float64(a * 120.0)); elseif (t <= 7.5e+57) tmp = fma(a, 120.0, Float64(x / Float64(Float64(z - t) * 0.016666666666666666))); else tmp = t_1; end return tmp end
code[x_, y_, z_, t_, a_] := Block[{t$95$1 = N[(-60.0 * N[(N[(x - y), $MachinePrecision] / t), $MachinePrecision] + N[(a * 120.0), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t, -4.6e-26], t$95$1, If[LessEqual[t, -5e-113], N[(N[(60.0 * N[(x - y), $MachinePrecision]), $MachinePrecision] / N[(z - t), $MachinePrecision]), $MachinePrecision], If[LessEqual[t, 1.05e-100], N[(60.0 * N[(N[(x - y), $MachinePrecision] / z), $MachinePrecision] + N[(a * 120.0), $MachinePrecision]), $MachinePrecision], If[LessEqual[t, 7.5e+57], N[(a * 120.0 + N[(x / N[(N[(z - t), $MachinePrecision] * 0.016666666666666666), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$1]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \mathsf{fma}\left(-60, \frac{x - y}{t}, a \cdot 120\right)\\
\mathbf{if}\;t \leq -4.6 \cdot 10^{-26}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;t \leq -5 \cdot 10^{-113}:\\
\;\;\;\;\frac{60 \cdot \left(x - y\right)}{z - t}\\
\mathbf{elif}\;t \leq 1.05 \cdot 10^{-100}:\\
\;\;\;\;\mathsf{fma}\left(60, \frac{x - y}{z}, a \cdot 120\right)\\
\mathbf{elif}\;t \leq 7.5 \cdot 10^{+57}:\\
\;\;\;\;\mathsf{fma}\left(a, 120, \frac{x}{\left(z - t\right) \cdot 0.016666666666666666}\right)\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if t < -4.60000000000000018e-26 or 7.5000000000000006e57 < t Initial program 99.8%
Taylor expanded in z around 0
lower-fma.f64N/A
lower-/.f64N/A
lower--.f64N/A
lower-*.f6489.6
Applied rewrites89.6%
if -4.60000000000000018e-26 < t < -4.9999999999999997e-113Initial program 95.5%
Taylor expanded in a around 0
associate-*r/N/A
lower-/.f64N/A
lower-*.f64N/A
lower--.f64N/A
lower--.f6481.7
Applied rewrites81.7%
if -4.9999999999999997e-113 < t < 1.05000000000000005e-100Initial program 98.9%
Taylor expanded in z around inf
lower-fma.f64N/A
lower-/.f64N/A
lower--.f64N/A
lower-*.f6490.8
Applied rewrites90.8%
if 1.05000000000000005e-100 < t < 7.5000000000000006e57Initial program 99.9%
Taylor expanded in x around inf
associate-*r/N/A
lower-/.f64N/A
lower-*.f64N/A
lower--.f6489.5
Applied rewrites89.5%
lift--.f64N/A
lift-*.f64N/A
lift-/.f64N/A
lift-*.f64N/A
+-commutativeN/A
lift-*.f64N/A
lower-fma.f6489.5
lift-/.f64N/A
lift-*.f64N/A
*-commutativeN/A
associate-/l*N/A
clear-numN/A
un-div-invN/A
lower-/.f64N/A
div-invN/A
lower-*.f64N/A
metadata-eval89.5
Applied rewrites89.5%
Final simplification89.4%
(FPCore (x y z t a)
:precision binary64
(if (<= (* a 120.0) -1e-57)
(fma a 120.0 (/ (* 60.0 x) z))
(if (<= (* a 120.0) 1e+35)
(/ (- x y) (* (- z t) 0.016666666666666666))
(* a 120.0))))
double code(double x, double y, double z, double t, double a) {
double tmp;
if ((a * 120.0) <= -1e-57) {
tmp = fma(a, 120.0, ((60.0 * x) / z));
} else if ((a * 120.0) <= 1e+35) {
tmp = (x - y) / ((z - t) * 0.016666666666666666);
} else {
tmp = a * 120.0;
}
return tmp;
}
function code(x, y, z, t, a) tmp = 0.0 if (Float64(a * 120.0) <= -1e-57) tmp = fma(a, 120.0, Float64(Float64(60.0 * x) / z)); elseif (Float64(a * 120.0) <= 1e+35) tmp = Float64(Float64(x - y) / Float64(Float64(z - t) * 0.016666666666666666)); else tmp = Float64(a * 120.0); end return tmp end
code[x_, y_, z_, t_, a_] := If[LessEqual[N[(a * 120.0), $MachinePrecision], -1e-57], N[(a * 120.0 + N[(N[(60.0 * x), $MachinePrecision] / z), $MachinePrecision]), $MachinePrecision], If[LessEqual[N[(a * 120.0), $MachinePrecision], 1e+35], N[(N[(x - y), $MachinePrecision] / N[(N[(z - t), $MachinePrecision] * 0.016666666666666666), $MachinePrecision]), $MachinePrecision], N[(a * 120.0), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;a \cdot 120 \leq -1 \cdot 10^{-57}:\\
\;\;\;\;\mathsf{fma}\left(a, 120, \frac{60 \cdot x}{z}\right)\\
\mathbf{elif}\;a \cdot 120 \leq 10^{+35}:\\
\;\;\;\;\frac{x - y}{\left(z - t\right) \cdot 0.016666666666666666}\\
\mathbf{else}:\\
\;\;\;\;a \cdot 120\\
\end{array}
\end{array}
if (*.f64 a #s(literal 120 binary64)) < -9.99999999999999955e-58Initial program 98.6%
Taylor expanded in x around inf
associate-*r/N/A
lower-/.f64N/A
lower-*.f64N/A
lower--.f6481.1
Applied rewrites81.1%
lift--.f64N/A
lift-*.f64N/A
lift-/.f64N/A
lift-*.f64N/A
+-commutativeN/A
lift-*.f64N/A
lower-fma.f6481.1
lift-/.f64N/A
lift-*.f64N/A
*-commutativeN/A
associate-/l*N/A
clear-numN/A
un-div-invN/A
lower-/.f64N/A
div-invN/A
lower-*.f64N/A
metadata-eval81.1
Applied rewrites81.1%
Taylor expanded in z around inf
associate-*r/N/A
lower-/.f64N/A
lower-*.f6472.3
Applied rewrites72.3%
if -9.99999999999999955e-58 < (*.f64 a #s(literal 120 binary64)) < 9.9999999999999997e34Initial program 99.8%
lift--.f64N/A
lift-*.f64N/A
remove-double-negN/A
lift--.f64N/A
remove-double-negN/A
lift-*.f64N/A
associate-/l*N/A
clear-numN/A
un-div-invN/A
lower-/.f64N/A
lower-/.f6499.7
Applied rewrites99.7%
Taylor expanded in a around 0
associate-*r/N/A
lower-/.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower--.f64N/A
lower--.f6482.9
Applied rewrites82.9%
lift--.f64N/A
lift--.f64N/A
associate-/l*N/A
clear-numN/A
un-div-invN/A
lower-/.f64N/A
div-invN/A
lower-*.f64N/A
metadata-eval82.9
Applied rewrites82.9%
if 9.9999999999999997e34 < (*.f64 a #s(literal 120 binary64)) Initial program 97.9%
Taylor expanded in z around inf
lower-*.f6490.2
Applied rewrites90.2%
Final simplification81.1%
(FPCore (x y z t a) :precision binary64 (if (<= (* a 120.0) -5e-5) (* a 120.0) (if (<= (* a 120.0) 1e+35) (/ (* 60.0 (- x y)) (- z t)) (* a 120.0))))
double code(double x, double y, double z, double t, double a) {
double tmp;
if ((a * 120.0) <= -5e-5) {
tmp = a * 120.0;
} else if ((a * 120.0) <= 1e+35) {
tmp = (60.0 * (x - y)) / (z - t);
} else {
tmp = a * 120.0;
}
return tmp;
}
real(8) function code(x, y, z, t, a)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8) :: tmp
if ((a * 120.0d0) <= (-5d-5)) then
tmp = a * 120.0d0
else if ((a * 120.0d0) <= 1d+35) then
tmp = (60.0d0 * (x - y)) / (z - t)
else
tmp = a * 120.0d0
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double tmp;
if ((a * 120.0) <= -5e-5) {
tmp = a * 120.0;
} else if ((a * 120.0) <= 1e+35) {
tmp = (60.0 * (x - y)) / (z - t);
} else {
tmp = a * 120.0;
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if (a * 120.0) <= -5e-5: tmp = a * 120.0 elif (a * 120.0) <= 1e+35: tmp = (60.0 * (x - y)) / (z - t) else: tmp = a * 120.0 return tmp
function code(x, y, z, t, a) tmp = 0.0 if (Float64(a * 120.0) <= -5e-5) tmp = Float64(a * 120.0); elseif (Float64(a * 120.0) <= 1e+35) tmp = Float64(Float64(60.0 * Float64(x - y)) / Float64(z - t)); else tmp = Float64(a * 120.0); end return tmp end
function tmp_2 = code(x, y, z, t, a) tmp = 0.0; if ((a * 120.0) <= -5e-5) tmp = a * 120.0; elseif ((a * 120.0) <= 1e+35) tmp = (60.0 * (x - y)) / (z - t); else tmp = a * 120.0; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[LessEqual[N[(a * 120.0), $MachinePrecision], -5e-5], N[(a * 120.0), $MachinePrecision], If[LessEqual[N[(a * 120.0), $MachinePrecision], 1e+35], N[(N[(60.0 * N[(x - y), $MachinePrecision]), $MachinePrecision] / N[(z - t), $MachinePrecision]), $MachinePrecision], N[(a * 120.0), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;a \cdot 120 \leq -5 \cdot 10^{-5}:\\
\;\;\;\;a \cdot 120\\
\mathbf{elif}\;a \cdot 120 \leq 10^{+35}:\\
\;\;\;\;\frac{60 \cdot \left(x - y\right)}{z - t}\\
\mathbf{else}:\\
\;\;\;\;a \cdot 120\\
\end{array}
\end{array}
if (*.f64 a #s(literal 120 binary64)) < -5.00000000000000024e-5 or 9.9999999999999997e34 < (*.f64 a #s(literal 120 binary64)) Initial program 98.3%
Taylor expanded in z around inf
lower-*.f6479.4
Applied rewrites79.4%
if -5.00000000000000024e-5 < (*.f64 a #s(literal 120 binary64)) < 9.9999999999999997e34Initial program 99.8%
Taylor expanded in a around 0
associate-*r/N/A
lower-/.f64N/A
lower-*.f64N/A
lower--.f64N/A
lower--.f6481.1
Applied rewrites81.1%
Final simplification80.3%
(FPCore (x y z t a)
:precision binary64
(let* ((t_1 (fma -60.0 (/ (- x y) t) (* a 120.0))))
(if (<= t -4.6e-26)
t_1
(if (<= t -5e-113)
(/ (* 60.0 (- x y)) (- z t))
(if (<= t 4.6e-19) (fma 60.0 (/ (- x y) z) (* a 120.0)) t_1)))))
double code(double x, double y, double z, double t, double a) {
double t_1 = fma(-60.0, ((x - y) / t), (a * 120.0));
double tmp;
if (t <= -4.6e-26) {
tmp = t_1;
} else if (t <= -5e-113) {
tmp = (60.0 * (x - y)) / (z - t);
} else if (t <= 4.6e-19) {
tmp = fma(60.0, ((x - y) / z), (a * 120.0));
} else {
tmp = t_1;
}
return tmp;
}
function code(x, y, z, t, a) t_1 = fma(-60.0, Float64(Float64(x - y) / t), Float64(a * 120.0)) tmp = 0.0 if (t <= -4.6e-26) tmp = t_1; elseif (t <= -5e-113) tmp = Float64(Float64(60.0 * Float64(x - y)) / Float64(z - t)); elseif (t <= 4.6e-19) tmp = fma(60.0, Float64(Float64(x - y) / z), Float64(a * 120.0)); else tmp = t_1; end return tmp end
code[x_, y_, z_, t_, a_] := Block[{t$95$1 = N[(-60.0 * N[(N[(x - y), $MachinePrecision] / t), $MachinePrecision] + N[(a * 120.0), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t, -4.6e-26], t$95$1, If[LessEqual[t, -5e-113], N[(N[(60.0 * N[(x - y), $MachinePrecision]), $MachinePrecision] / N[(z - t), $MachinePrecision]), $MachinePrecision], If[LessEqual[t, 4.6e-19], N[(60.0 * N[(N[(x - y), $MachinePrecision] / z), $MachinePrecision] + N[(a * 120.0), $MachinePrecision]), $MachinePrecision], t$95$1]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \mathsf{fma}\left(-60, \frac{x - y}{t}, a \cdot 120\right)\\
\mathbf{if}\;t \leq -4.6 \cdot 10^{-26}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;t \leq -5 \cdot 10^{-113}:\\
\;\;\;\;\frac{60 \cdot \left(x - y\right)}{z - t}\\
\mathbf{elif}\;t \leq 4.6 \cdot 10^{-19}:\\
\;\;\;\;\mathsf{fma}\left(60, \frac{x - y}{z}, a \cdot 120\right)\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if t < -4.60000000000000018e-26 or 4.5999999999999996e-19 < t Initial program 99.8%
Taylor expanded in z around 0
lower-fma.f64N/A
lower-/.f64N/A
lower--.f64N/A
lower-*.f6487.4
Applied rewrites87.4%
if -4.60000000000000018e-26 < t < -4.9999999999999997e-113Initial program 95.5%
Taylor expanded in a around 0
associate-*r/N/A
lower-/.f64N/A
lower-*.f64N/A
lower--.f64N/A
lower--.f6481.7
Applied rewrites81.7%
if -4.9999999999999997e-113 < t < 4.5999999999999996e-19Initial program 99.0%
Taylor expanded in z around inf
lower-fma.f64N/A
lower-/.f64N/A
lower--.f64N/A
lower-*.f6488.5
Applied rewrites88.5%
Final simplification87.4%
(FPCore (x y z t a) :precision binary64 (if (<= (* a 120.0) -5e-80) (* a 120.0) (if (<= (* a 120.0) 4e+14) (* -60.0 (/ y (- z t))) (* a 120.0))))
double code(double x, double y, double z, double t, double a) {
double tmp;
if ((a * 120.0) <= -5e-80) {
tmp = a * 120.0;
} else if ((a * 120.0) <= 4e+14) {
tmp = -60.0 * (y / (z - t));
} else {
tmp = a * 120.0;
}
return tmp;
}
real(8) function code(x, y, z, t, a)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8) :: tmp
if ((a * 120.0d0) <= (-5d-80)) then
tmp = a * 120.0d0
else if ((a * 120.0d0) <= 4d+14) then
tmp = (-60.0d0) * (y / (z - t))
else
tmp = a * 120.0d0
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double tmp;
if ((a * 120.0) <= -5e-80) {
tmp = a * 120.0;
} else if ((a * 120.0) <= 4e+14) {
tmp = -60.0 * (y / (z - t));
} else {
tmp = a * 120.0;
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if (a * 120.0) <= -5e-80: tmp = a * 120.0 elif (a * 120.0) <= 4e+14: tmp = -60.0 * (y / (z - t)) else: tmp = a * 120.0 return tmp
function code(x, y, z, t, a) tmp = 0.0 if (Float64(a * 120.0) <= -5e-80) tmp = Float64(a * 120.0); elseif (Float64(a * 120.0) <= 4e+14) tmp = Float64(-60.0 * Float64(y / Float64(z - t))); else tmp = Float64(a * 120.0); end return tmp end
function tmp_2 = code(x, y, z, t, a) tmp = 0.0; if ((a * 120.0) <= -5e-80) tmp = a * 120.0; elseif ((a * 120.0) <= 4e+14) tmp = -60.0 * (y / (z - t)); else tmp = a * 120.0; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[LessEqual[N[(a * 120.0), $MachinePrecision], -5e-80], N[(a * 120.0), $MachinePrecision], If[LessEqual[N[(a * 120.0), $MachinePrecision], 4e+14], N[(-60.0 * N[(y / N[(z - t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(a * 120.0), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;a \cdot 120 \leq -5 \cdot 10^{-80}:\\
\;\;\;\;a \cdot 120\\
\mathbf{elif}\;a \cdot 120 \leq 4 \cdot 10^{+14}:\\
\;\;\;\;-60 \cdot \frac{y}{z - t}\\
\mathbf{else}:\\
\;\;\;\;a \cdot 120\\
\end{array}
\end{array}
if (*.f64 a #s(literal 120 binary64)) < -5e-80 or 4e14 < (*.f64 a #s(literal 120 binary64)) Initial program 98.4%
Taylor expanded in z around inf
lower-*.f6476.3
Applied rewrites76.3%
if -5e-80 < (*.f64 a #s(literal 120 binary64)) < 4e14Initial program 99.8%
Taylor expanded in y around inf
lower-*.f64N/A
lower-/.f64N/A
lower--.f6447.8
Applied rewrites47.8%
Final simplification62.4%
(FPCore (x y z t a) :precision binary64 (if (<= x -1.55e+164) (* x (/ -60.0 t)) (if (<= x 6.8e+251) (* a 120.0) (* x (/ 60.0 z)))))
double code(double x, double y, double z, double t, double a) {
double tmp;
if (x <= -1.55e+164) {
tmp = x * (-60.0 / t);
} else if (x <= 6.8e+251) {
tmp = a * 120.0;
} else {
tmp = x * (60.0 / z);
}
return tmp;
}
real(8) function code(x, y, z, t, a)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8) :: tmp
if (x <= (-1.55d+164)) then
tmp = x * ((-60.0d0) / t)
else if (x <= 6.8d+251) then
tmp = a * 120.0d0
else
tmp = x * (60.0d0 / z)
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double tmp;
if (x <= -1.55e+164) {
tmp = x * (-60.0 / t);
} else if (x <= 6.8e+251) {
tmp = a * 120.0;
} else {
tmp = x * (60.0 / z);
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if x <= -1.55e+164: tmp = x * (-60.0 / t) elif x <= 6.8e+251: tmp = a * 120.0 else: tmp = x * (60.0 / z) return tmp
function code(x, y, z, t, a) tmp = 0.0 if (x <= -1.55e+164) tmp = Float64(x * Float64(-60.0 / t)); elseif (x <= 6.8e+251) tmp = Float64(a * 120.0); else tmp = Float64(x * Float64(60.0 / z)); end return tmp end
function tmp_2 = code(x, y, z, t, a) tmp = 0.0; if (x <= -1.55e+164) tmp = x * (-60.0 / t); elseif (x <= 6.8e+251) tmp = a * 120.0; else tmp = x * (60.0 / z); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[LessEqual[x, -1.55e+164], N[(x * N[(-60.0 / t), $MachinePrecision]), $MachinePrecision], If[LessEqual[x, 6.8e+251], N[(a * 120.0), $MachinePrecision], N[(x * N[(60.0 / z), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq -1.55 \cdot 10^{+164}:\\
\;\;\;\;x \cdot \frac{-60}{t}\\
\mathbf{elif}\;x \leq 6.8 \cdot 10^{+251}:\\
\;\;\;\;a \cdot 120\\
\mathbf{else}:\\
\;\;\;\;x \cdot \frac{60}{z}\\
\end{array}
\end{array}
if x < -1.5500000000000001e164Initial program 99.6%
Taylor expanded in x around inf
associate-*r/N/A
lower-/.f64N/A
lower-*.f64N/A
lower--.f6475.6
Applied rewrites75.6%
Taylor expanded in z around 0
associate-*r/N/A
lower-/.f64N/A
lower-*.f6452.8
Applied rewrites52.8%
*-commutativeN/A
associate-/l*N/A
lower-*.f64N/A
lower-/.f6452.8
Applied rewrites52.8%
if -1.5500000000000001e164 < x < 6.80000000000000023e251Initial program 99.4%
Taylor expanded in z around inf
lower-*.f6453.6
Applied rewrites53.6%
if 6.80000000000000023e251 < x Initial program 88.9%
Taylor expanded in x around inf
associate-*r/N/A
lower-/.f64N/A
lower-*.f64N/A
lower--.f6485.5
Applied rewrites85.5%
*-commutativeN/A
lift--.f64N/A
associate-/l*N/A
lower-*.f64N/A
lower-/.f6485.4
Applied rewrites85.4%
Taylor expanded in z around inf
associate-*r/N/A
*-commutativeN/A
associate-/l*N/A
metadata-evalN/A
associate-*r/N/A
lower-*.f64N/A
associate-*r/N/A
metadata-evalN/A
lower-/.f6474.8
Applied rewrites74.8%
Final simplification54.2%
(FPCore (x y z t a) :precision binary64 (fma (/ 60.0 (- z t)) (- x y) (* a 120.0)))
double code(double x, double y, double z, double t, double a) {
return fma((60.0 / (z - t)), (x - y), (a * 120.0));
}
function code(x, y, z, t, a) return fma(Float64(60.0 / Float64(z - t)), Float64(x - y), Float64(a * 120.0)) end
code[x_, y_, z_, t_, a_] := N[(N[(60.0 / N[(z - t), $MachinePrecision]), $MachinePrecision] * N[(x - y), $MachinePrecision] + N[(a * 120.0), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\mathsf{fma}\left(\frac{60}{z - t}, x - y, a \cdot 120\right)
\end{array}
Initial program 99.1%
lift--.f64N/A
lift-*.f64N/A
remove-double-negN/A
lift--.f64N/A
remove-double-negN/A
lift-*.f64N/A
*-commutativeN/A
associate-/l*N/A
*-commutativeN/A
lift-*.f64N/A
lower-fma.f64N/A
lower-/.f6499.8
Applied rewrites99.8%
(FPCore (x y z t a) :precision binary64 (fma a 120.0 (/ (* (- x y) -60.0) (- t z))))
double code(double x, double y, double z, double t, double a) {
return fma(a, 120.0, (((x - y) * -60.0) / (t - z)));
}
function code(x, y, z, t, a) return fma(a, 120.0, Float64(Float64(Float64(x - y) * -60.0) / Float64(t - z))) end
code[x_, y_, z_, t_, a_] := N[(a * 120.0 + N[(N[(N[(x - y), $MachinePrecision] * -60.0), $MachinePrecision] / N[(t - z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\mathsf{fma}\left(a, 120, \frac{\left(x - y\right) \cdot -60}{t - z}\right)
\end{array}
Initial program 99.1%
lift--.f64N/A
lift-*.f64N/A
remove-double-negN/A
lift--.f64N/A
remove-double-negN/A
lift-/.f64N/A
lift-*.f64N/A
+-commutativeN/A
lift-*.f64N/A
lower-fma.f6499.1
lift-/.f64N/A
frac-2negN/A
lower-/.f64N/A
lift-*.f64N/A
*-commutativeN/A
distribute-rgt-neg-inN/A
lower-*.f64N/A
metadata-evalN/A
neg-sub0N/A
lift--.f64N/A
sub-negN/A
+-commutativeN/A
associate--r+N/A
neg-sub0N/A
remove-double-negN/A
lower--.f6499.1
Applied rewrites99.1%
(FPCore (x y z t a) :precision binary64 (if (<= x -1.55e+164) (* x (/ -60.0 t)) (* a 120.0)))
double code(double x, double y, double z, double t, double a) {
double tmp;
if (x <= -1.55e+164) {
tmp = x * (-60.0 / t);
} else {
tmp = a * 120.0;
}
return tmp;
}
real(8) function code(x, y, z, t, a)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8) :: tmp
if (x <= (-1.55d+164)) then
tmp = x * ((-60.0d0) / t)
else
tmp = a * 120.0d0
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a) {
double tmp;
if (x <= -1.55e+164) {
tmp = x * (-60.0 / t);
} else {
tmp = a * 120.0;
}
return tmp;
}
def code(x, y, z, t, a): tmp = 0 if x <= -1.55e+164: tmp = x * (-60.0 / t) else: tmp = a * 120.0 return tmp
function code(x, y, z, t, a) tmp = 0.0 if (x <= -1.55e+164) tmp = Float64(x * Float64(-60.0 / t)); else tmp = Float64(a * 120.0); end return tmp end
function tmp_2 = code(x, y, z, t, a) tmp = 0.0; if (x <= -1.55e+164) tmp = x * (-60.0 / t); else tmp = a * 120.0; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_] := If[LessEqual[x, -1.55e+164], N[(x * N[(-60.0 / t), $MachinePrecision]), $MachinePrecision], N[(a * 120.0), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq -1.55 \cdot 10^{+164}:\\
\;\;\;\;x \cdot \frac{-60}{t}\\
\mathbf{else}:\\
\;\;\;\;a \cdot 120\\
\end{array}
\end{array}
if x < -1.5500000000000001e164Initial program 99.6%
Taylor expanded in x around inf
associate-*r/N/A
lower-/.f64N/A
lower-*.f64N/A
lower--.f6475.6
Applied rewrites75.6%
Taylor expanded in z around 0
associate-*r/N/A
lower-/.f64N/A
lower-*.f6452.8
Applied rewrites52.8%
*-commutativeN/A
associate-/l*N/A
lower-*.f64N/A
lower-/.f6452.8
Applied rewrites52.8%
if -1.5500000000000001e164 < x Initial program 99.0%
Taylor expanded in z around inf
lower-*.f6452.0
Applied rewrites52.0%
Final simplification52.1%
(FPCore (x y z t a) :precision binary64 (* a 120.0))
double code(double x, double y, double z, double t, double a) {
return a * 120.0;
}
real(8) function code(x, y, z, t, a)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
code = a * 120.0d0
end function
public static double code(double x, double y, double z, double t, double a) {
return a * 120.0;
}
def code(x, y, z, t, a): return a * 120.0
function code(x, y, z, t, a) return Float64(a * 120.0) end
function tmp = code(x, y, z, t, a) tmp = a * 120.0; end
code[x_, y_, z_, t_, a_] := N[(a * 120.0), $MachinePrecision]
\begin{array}{l}
\\
a \cdot 120
\end{array}
Initial program 99.1%
Taylor expanded in z around inf
lower-*.f6447.9
Applied rewrites47.9%
Final simplification47.9%
(FPCore (x y z t a) :precision binary64 (+ (/ 60.0 (/ (- z t) (- x y))) (* a 120.0)))
double code(double x, double y, double z, double t, double a) {
return (60.0 / ((z - t) / (x - y))) + (a * 120.0);
}
real(8) function code(x, y, z, t, a)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
code = (60.0d0 / ((z - t) / (x - y))) + (a * 120.0d0)
end function
public static double code(double x, double y, double z, double t, double a) {
return (60.0 / ((z - t) / (x - y))) + (a * 120.0);
}
def code(x, y, z, t, a): return (60.0 / ((z - t) / (x - y))) + (a * 120.0)
function code(x, y, z, t, a) return Float64(Float64(60.0 / Float64(Float64(z - t) / Float64(x - y))) + Float64(a * 120.0)) end
function tmp = code(x, y, z, t, a) tmp = (60.0 / ((z - t) / (x - y))) + (a * 120.0); end
code[x_, y_, z_, t_, a_] := N[(N[(60.0 / N[(N[(z - t), $MachinePrecision] / N[(x - y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(a * 120.0), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{60}{\frac{z - t}{x - y}} + a \cdot 120
\end{array}
herbie shell --seed 2024219
(FPCore (x y z t a)
:name "Data.Colour.RGB:hslsv from colour-2.3.3, B"
:precision binary64
:alt
(! :herbie-platform default (+ (/ 60 (/ (- z t) (- x y))) (* a 120)))
(+ (/ (* 60.0 (- x y)) (- z t)) (* a 120.0)))