
(FPCore (x y z t a b c i j k) :precision binary64 (- (- (+ (- (* (* (* (* x 18.0) y) z) t) (* (* a 4.0) t)) (* b c)) (* (* x 4.0) i)) (* (* j 27.0) k)))
double code(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) {
return (((((((x * 18.0) * y) * z) * t) - ((a * 4.0) * t)) + (b * c)) - ((x * 4.0) * i)) - ((j * 27.0) * k);
}
module fmin_fmax_functions
implicit none
private
public fmax
public fmin
interface fmax
module procedure fmax88
module procedure fmax44
module procedure fmax84
module procedure fmax48
end interface
interface fmin
module procedure fmin88
module procedure fmin44
module procedure fmin84
module procedure fmin48
end interface
contains
real(8) function fmax88(x, y) result (res)
real(8), intent (in) :: x
real(8), intent (in) :: y
res = merge(y, merge(x, max(x, y), y /= y), x /= x)
end function
real(4) function fmax44(x, y) result (res)
real(4), intent (in) :: x
real(4), intent (in) :: y
res = merge(y, merge(x, max(x, y), y /= y), x /= x)
end function
real(8) function fmax84(x, y) result(res)
real(8), intent (in) :: x
real(4), intent (in) :: y
res = merge(dble(y), merge(x, max(x, dble(y)), y /= y), x /= x)
end function
real(8) function fmax48(x, y) result(res)
real(4), intent (in) :: x
real(8), intent (in) :: y
res = merge(y, merge(dble(x), max(dble(x), y), y /= y), x /= x)
end function
real(8) function fmin88(x, y) result (res)
real(8), intent (in) :: x
real(8), intent (in) :: y
res = merge(y, merge(x, min(x, y), y /= y), x /= x)
end function
real(4) function fmin44(x, y) result (res)
real(4), intent (in) :: x
real(4), intent (in) :: y
res = merge(y, merge(x, min(x, y), y /= y), x /= x)
end function
real(8) function fmin84(x, y) result(res)
real(8), intent (in) :: x
real(4), intent (in) :: y
res = merge(dble(y), merge(x, min(x, dble(y)), y /= y), x /= x)
end function
real(8) function fmin48(x, y) result(res)
real(4), intent (in) :: x
real(8), intent (in) :: y
res = merge(y, merge(dble(x), min(dble(x), y), y /= y), x /= x)
end function
end module
real(8) function code(x, y, z, t, a, b, c, i, j, k)
use fmin_fmax_functions
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), intent (in) :: b
real(8), intent (in) :: c
real(8), intent (in) :: i
real(8), intent (in) :: j
real(8), intent (in) :: k
code = (((((((x * 18.0d0) * y) * z) * t) - ((a * 4.0d0) * t)) + (b * c)) - ((x * 4.0d0) * i)) - ((j * 27.0d0) * k)
end function
public static double code(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) {
return (((((((x * 18.0) * y) * z) * t) - ((a * 4.0) * t)) + (b * c)) - ((x * 4.0) * i)) - ((j * 27.0) * k);
}
def code(x, y, z, t, a, b, c, i, j, k): return (((((((x * 18.0) * y) * z) * t) - ((a * 4.0) * t)) + (b * c)) - ((x * 4.0) * i)) - ((j * 27.0) * k)
function code(x, y, z, t, a, b, c, i, j, k) return Float64(Float64(Float64(Float64(Float64(Float64(Float64(Float64(x * 18.0) * y) * z) * t) - Float64(Float64(a * 4.0) * t)) + Float64(b * c)) - Float64(Float64(x * 4.0) * i)) - Float64(Float64(j * 27.0) * k)) end
function tmp = code(x, y, z, t, a, b, c, i, j, k) tmp = (((((((x * 18.0) * y) * z) * t) - ((a * 4.0) * t)) + (b * c)) - ((x * 4.0) * i)) - ((j * 27.0) * k); end
code[x_, y_, z_, t_, a_, b_, c_, i_, j_, k_] := N[(N[(N[(N[(N[(N[(N[(N[(x * 18.0), $MachinePrecision] * y), $MachinePrecision] * z), $MachinePrecision] * t), $MachinePrecision] - N[(N[(a * 4.0), $MachinePrecision] * t), $MachinePrecision]), $MachinePrecision] + N[(b * c), $MachinePrecision]), $MachinePrecision] - N[(N[(x * 4.0), $MachinePrecision] * i), $MachinePrecision]), $MachinePrecision] - N[(N[(j * 27.0), $MachinePrecision] * k), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\left(\left(\left(\left(\left(\left(x \cdot 18\right) \cdot y\right) \cdot z\right) \cdot t - \left(a \cdot 4\right) \cdot t\right) + b \cdot c\right) - \left(x \cdot 4\right) \cdot i\right) - \left(j \cdot 27\right) \cdot k
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 18 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (x y z t a b c i j k) :precision binary64 (- (- (+ (- (* (* (* (* x 18.0) y) z) t) (* (* a 4.0) t)) (* b c)) (* (* x 4.0) i)) (* (* j 27.0) k)))
double code(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) {
return (((((((x * 18.0) * y) * z) * t) - ((a * 4.0) * t)) + (b * c)) - ((x * 4.0) * i)) - ((j * 27.0) * k);
}
module fmin_fmax_functions
implicit none
private
public fmax
public fmin
interface fmax
module procedure fmax88
module procedure fmax44
module procedure fmax84
module procedure fmax48
end interface
interface fmin
module procedure fmin88
module procedure fmin44
module procedure fmin84
module procedure fmin48
end interface
contains
real(8) function fmax88(x, y) result (res)
real(8), intent (in) :: x
real(8), intent (in) :: y
res = merge(y, merge(x, max(x, y), y /= y), x /= x)
end function
real(4) function fmax44(x, y) result (res)
real(4), intent (in) :: x
real(4), intent (in) :: y
res = merge(y, merge(x, max(x, y), y /= y), x /= x)
end function
real(8) function fmax84(x, y) result(res)
real(8), intent (in) :: x
real(4), intent (in) :: y
res = merge(dble(y), merge(x, max(x, dble(y)), y /= y), x /= x)
end function
real(8) function fmax48(x, y) result(res)
real(4), intent (in) :: x
real(8), intent (in) :: y
res = merge(y, merge(dble(x), max(dble(x), y), y /= y), x /= x)
end function
real(8) function fmin88(x, y) result (res)
real(8), intent (in) :: x
real(8), intent (in) :: y
res = merge(y, merge(x, min(x, y), y /= y), x /= x)
end function
real(4) function fmin44(x, y) result (res)
real(4), intent (in) :: x
real(4), intent (in) :: y
res = merge(y, merge(x, min(x, y), y /= y), x /= x)
end function
real(8) function fmin84(x, y) result(res)
real(8), intent (in) :: x
real(4), intent (in) :: y
res = merge(dble(y), merge(x, min(x, dble(y)), y /= y), x /= x)
end function
real(8) function fmin48(x, y) result(res)
real(4), intent (in) :: x
real(8), intent (in) :: y
res = merge(y, merge(dble(x), min(dble(x), y), y /= y), x /= x)
end function
end module
real(8) function code(x, y, z, t, a, b, c, i, j, k)
use fmin_fmax_functions
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), intent (in) :: b
real(8), intent (in) :: c
real(8), intent (in) :: i
real(8), intent (in) :: j
real(8), intent (in) :: k
code = (((((((x * 18.0d0) * y) * z) * t) - ((a * 4.0d0) * t)) + (b * c)) - ((x * 4.0d0) * i)) - ((j * 27.0d0) * k)
end function
public static double code(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) {
return (((((((x * 18.0) * y) * z) * t) - ((a * 4.0) * t)) + (b * c)) - ((x * 4.0) * i)) - ((j * 27.0) * k);
}
def code(x, y, z, t, a, b, c, i, j, k): return (((((((x * 18.0) * y) * z) * t) - ((a * 4.0) * t)) + (b * c)) - ((x * 4.0) * i)) - ((j * 27.0) * k)
function code(x, y, z, t, a, b, c, i, j, k) return Float64(Float64(Float64(Float64(Float64(Float64(Float64(Float64(x * 18.0) * y) * z) * t) - Float64(Float64(a * 4.0) * t)) + Float64(b * c)) - Float64(Float64(x * 4.0) * i)) - Float64(Float64(j * 27.0) * k)) end
function tmp = code(x, y, z, t, a, b, c, i, j, k) tmp = (((((((x * 18.0) * y) * z) * t) - ((a * 4.0) * t)) + (b * c)) - ((x * 4.0) * i)) - ((j * 27.0) * k); end
code[x_, y_, z_, t_, a_, b_, c_, i_, j_, k_] := N[(N[(N[(N[(N[(N[(N[(N[(x * 18.0), $MachinePrecision] * y), $MachinePrecision] * z), $MachinePrecision] * t), $MachinePrecision] - N[(N[(a * 4.0), $MachinePrecision] * t), $MachinePrecision]), $MachinePrecision] + N[(b * c), $MachinePrecision]), $MachinePrecision] - N[(N[(x * 4.0), $MachinePrecision] * i), $MachinePrecision]), $MachinePrecision] - N[(N[(j * 27.0), $MachinePrecision] * k), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\left(\left(\left(\left(\left(\left(x \cdot 18\right) \cdot y\right) \cdot z\right) \cdot t - \left(a \cdot 4\right) \cdot t\right) + b \cdot c\right) - \left(x \cdot 4\right) \cdot i\right) - \left(j \cdot 27\right) \cdot k
\end{array}
NOTE: x, y, z, t, a, b, c, i, j, and k should be sorted in increasing order before calling this function.
(FPCore (x y z t a b c i j k)
:precision binary64
(if (<=
(-
(-
(+ (- (* (* (* (* x 18.0) y) z) t) (* (* a 4.0) t)) (* b c))
(* (* x 4.0) i))
(* (* j 27.0) k))
INFINITY)
(-
(fma (fma z (* y (* 18.0 x)) (* -4.0 a)) t (fma c b (* (* -4.0 x) i)))
(* (* k 27.0) j))
(* (fma y (* z (* t 18.0)) (* i -4.0)) x)))assert(x < y && y < z && z < t && t < a && a < b && b < c && c < i && i < j && j < k);
double code(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) {
double tmp;
if (((((((((x * 18.0) * y) * z) * t) - ((a * 4.0) * t)) + (b * c)) - ((x * 4.0) * i)) - ((j * 27.0) * k)) <= ((double) INFINITY)) {
tmp = fma(fma(z, (y * (18.0 * x)), (-4.0 * a)), t, fma(c, b, ((-4.0 * x) * i))) - ((k * 27.0) * j);
} else {
tmp = fma(y, (z * (t * 18.0)), (i * -4.0)) * x;
}
return tmp;
}
x, y, z, t, a, b, c, i, j, k = sort([x, y, z, t, a, b, c, i, j, k]) function code(x, y, z, t, a, b, c, i, j, k) tmp = 0.0 if (Float64(Float64(Float64(Float64(Float64(Float64(Float64(Float64(x * 18.0) * y) * z) * t) - Float64(Float64(a * 4.0) * t)) + Float64(b * c)) - Float64(Float64(x * 4.0) * i)) - Float64(Float64(j * 27.0) * k)) <= Inf) tmp = Float64(fma(fma(z, Float64(y * Float64(18.0 * x)), Float64(-4.0 * a)), t, fma(c, b, Float64(Float64(-4.0 * x) * i))) - Float64(Float64(k * 27.0) * j)); else tmp = Float64(fma(y, Float64(z * Float64(t * 18.0)), Float64(i * -4.0)) * x); end return tmp end
NOTE: x, y, z, t, a, b, c, i, j, and k should be sorted in increasing order before calling this function. code[x_, y_, z_, t_, a_, b_, c_, i_, j_, k_] := If[LessEqual[N[(N[(N[(N[(N[(N[(N[(N[(x * 18.0), $MachinePrecision] * y), $MachinePrecision] * z), $MachinePrecision] * t), $MachinePrecision] - N[(N[(a * 4.0), $MachinePrecision] * t), $MachinePrecision]), $MachinePrecision] + N[(b * c), $MachinePrecision]), $MachinePrecision] - N[(N[(x * 4.0), $MachinePrecision] * i), $MachinePrecision]), $MachinePrecision] - N[(N[(j * 27.0), $MachinePrecision] * k), $MachinePrecision]), $MachinePrecision], Infinity], N[(N[(N[(z * N[(y * N[(18.0 * x), $MachinePrecision]), $MachinePrecision] + N[(-4.0 * a), $MachinePrecision]), $MachinePrecision] * t + N[(c * b + N[(N[(-4.0 * x), $MachinePrecision] * i), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - N[(N[(k * 27.0), $MachinePrecision] * j), $MachinePrecision]), $MachinePrecision], N[(N[(y * N[(z * N[(t * 18.0), $MachinePrecision]), $MachinePrecision] + N[(i * -4.0), $MachinePrecision]), $MachinePrecision] * x), $MachinePrecision]]
\begin{array}{l}
[x, y, z, t, a, b, c, i, j, k] = \mathsf{sort}([x, y, z, t, a, b, c, i, j, k])\\
\\
\begin{array}{l}
\mathbf{if}\;\left(\left(\left(\left(\left(\left(x \cdot 18\right) \cdot y\right) \cdot z\right) \cdot t - \left(a \cdot 4\right) \cdot t\right) + b \cdot c\right) - \left(x \cdot 4\right) \cdot i\right) - \left(j \cdot 27\right) \cdot k \leq \infty:\\
\;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(z, y \cdot \left(18 \cdot x\right), -4 \cdot a\right), t, \mathsf{fma}\left(c, b, \left(-4 \cdot x\right) \cdot i\right)\right) - \left(k \cdot 27\right) \cdot j\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(y, z \cdot \left(t \cdot 18\right), i \cdot -4\right) \cdot x\\
\end{array}
\end{array}
if (-.f64 (-.f64 (+.f64 (-.f64 (*.f64 (*.f64 (*.f64 (*.f64 x #s(literal 18 binary64)) y) z) t) (*.f64 (*.f64 a #s(literal 4 binary64)) t)) (*.f64 b c)) (*.f64 (*.f64 x #s(literal 4 binary64)) i)) (*.f64 (*.f64 j #s(literal 27 binary64)) k)) < +inf.0Initial program 95.6%
Applied rewrites95.7%
if +inf.0 < (-.f64 (-.f64 (+.f64 (-.f64 (*.f64 (*.f64 (*.f64 (*.f64 x #s(literal 18 binary64)) y) z) t) (*.f64 (*.f64 a #s(literal 4 binary64)) t)) (*.f64 b c)) (*.f64 (*.f64 x #s(literal 4 binary64)) i)) (*.f64 (*.f64 j #s(literal 27 binary64)) k)) Initial program 0.0%
Taylor expanded in x around inf
*-commutativeN/A
lower-*.f64N/A
fp-cancel-sub-sign-invN/A
metadata-evalN/A
+-commutativeN/A
lower-fma.f64N/A
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f6474.5
Applied rewrites74.5%
Applied rewrites74.5%
Final simplification93.1%
NOTE: x, y, z, t, a, b, c, i, j, and k should be sorted in increasing order before calling this function.
(FPCore (x y z t a b c i j k)
:precision binary64
(let* ((t_1 (* (* j 27.0) k)))
(if (<=
(-
(-
(+ (- (* (* (* (* x 18.0) y) z) t) (* (* a 4.0) t)) (* b c))
(* (* x 4.0) i))
t_1)
INFINITY)
(-
(fma (* -4.0 x) i (fma (fma z (* y (* 18.0 x)) (* -4.0 a)) t (* c b)))
t_1)
(* (fma y (* z (* t 18.0)) (* i -4.0)) x))))assert(x < y && y < z && z < t && t < a && a < b && b < c && c < i && i < j && j < k);
double code(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) {
double t_1 = (j * 27.0) * k;
double tmp;
if (((((((((x * 18.0) * y) * z) * t) - ((a * 4.0) * t)) + (b * c)) - ((x * 4.0) * i)) - t_1) <= ((double) INFINITY)) {
tmp = fma((-4.0 * x), i, fma(fma(z, (y * (18.0 * x)), (-4.0 * a)), t, (c * b))) - t_1;
} else {
tmp = fma(y, (z * (t * 18.0)), (i * -4.0)) * x;
}
return tmp;
}
x, y, z, t, a, b, c, i, j, k = sort([x, y, z, t, a, b, c, i, j, k]) function code(x, y, z, t, a, b, c, i, j, k) t_1 = Float64(Float64(j * 27.0) * k) tmp = 0.0 if (Float64(Float64(Float64(Float64(Float64(Float64(Float64(Float64(x * 18.0) * y) * z) * t) - Float64(Float64(a * 4.0) * t)) + Float64(b * c)) - Float64(Float64(x * 4.0) * i)) - t_1) <= Inf) tmp = Float64(fma(Float64(-4.0 * x), i, fma(fma(z, Float64(y * Float64(18.0 * x)), Float64(-4.0 * a)), t, Float64(c * b))) - t_1); else tmp = Float64(fma(y, Float64(z * Float64(t * 18.0)), Float64(i * -4.0)) * x); end return tmp end
NOTE: x, y, z, t, a, b, c, i, j, and k should be sorted in increasing order before calling this function.
code[x_, y_, z_, t_, a_, b_, c_, i_, j_, k_] := Block[{t$95$1 = N[(N[(j * 27.0), $MachinePrecision] * k), $MachinePrecision]}, If[LessEqual[N[(N[(N[(N[(N[(N[(N[(N[(x * 18.0), $MachinePrecision] * y), $MachinePrecision] * z), $MachinePrecision] * t), $MachinePrecision] - N[(N[(a * 4.0), $MachinePrecision] * t), $MachinePrecision]), $MachinePrecision] + N[(b * c), $MachinePrecision]), $MachinePrecision] - N[(N[(x * 4.0), $MachinePrecision] * i), $MachinePrecision]), $MachinePrecision] - t$95$1), $MachinePrecision], Infinity], N[(N[(N[(-4.0 * x), $MachinePrecision] * i + N[(N[(z * N[(y * N[(18.0 * x), $MachinePrecision]), $MachinePrecision] + N[(-4.0 * a), $MachinePrecision]), $MachinePrecision] * t + N[(c * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - t$95$1), $MachinePrecision], N[(N[(y * N[(z * N[(t * 18.0), $MachinePrecision]), $MachinePrecision] + N[(i * -4.0), $MachinePrecision]), $MachinePrecision] * x), $MachinePrecision]]]
\begin{array}{l}
[x, y, z, t, a, b, c, i, j, k] = \mathsf{sort}([x, y, z, t, a, b, c, i, j, k])\\
\\
\begin{array}{l}
t_1 := \left(j \cdot 27\right) \cdot k\\
\mathbf{if}\;\left(\left(\left(\left(\left(\left(x \cdot 18\right) \cdot y\right) \cdot z\right) \cdot t - \left(a \cdot 4\right) \cdot t\right) + b \cdot c\right) - \left(x \cdot 4\right) \cdot i\right) - t\_1 \leq \infty:\\
\;\;\;\;\mathsf{fma}\left(-4 \cdot x, i, \mathsf{fma}\left(\mathsf{fma}\left(z, y \cdot \left(18 \cdot x\right), -4 \cdot a\right), t, c \cdot b\right)\right) - t\_1\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(y, z \cdot \left(t \cdot 18\right), i \cdot -4\right) \cdot x\\
\end{array}
\end{array}
if (-.f64 (-.f64 (+.f64 (-.f64 (*.f64 (*.f64 (*.f64 (*.f64 x #s(literal 18 binary64)) y) z) t) (*.f64 (*.f64 a #s(literal 4 binary64)) t)) (*.f64 b c)) (*.f64 (*.f64 x #s(literal 4 binary64)) i)) (*.f64 (*.f64 j #s(literal 27 binary64)) k)) < +inf.0Initial program 95.6%
lift--.f64N/A
lift-*.f64N/A
fp-cancel-sub-sign-invN/A
+-commutativeN/A
lower-fma.f64N/A
lift-*.f64N/A
*-commutativeN/A
distribute-lft-neg-inN/A
lower-*.f64N/A
metadata-eval95.6
lift-+.f64N/A
lift--.f64N/A
lift-*.f64N/A
lift-*.f64N/A
distribute-rgt-out--N/A
*-commutativeN/A
Applied rewrites95.6%
if +inf.0 < (-.f64 (-.f64 (+.f64 (-.f64 (*.f64 (*.f64 (*.f64 (*.f64 x #s(literal 18 binary64)) y) z) t) (*.f64 (*.f64 a #s(literal 4 binary64)) t)) (*.f64 b c)) (*.f64 (*.f64 x #s(literal 4 binary64)) i)) (*.f64 (*.f64 j #s(literal 27 binary64)) k)) Initial program 0.0%
Taylor expanded in x around inf
*-commutativeN/A
lower-*.f64N/A
fp-cancel-sub-sign-invN/A
metadata-evalN/A
+-commutativeN/A
lower-fma.f64N/A
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f6474.5
Applied rewrites74.5%
Applied rewrites74.5%
NOTE: x, y, z, t, a, b, c, i, j, and k should be sorted in increasing order before calling this function.
(FPCore (x y z t a b c i j k)
:precision binary64
(if (<=
(-
(-
(+ (- (* (* (* (* x 18.0) y) z) t) (* (* a 4.0) t)) (* b c))
(* (* x 4.0) i))
(* (* j 27.0) k))
INFINITY)
(fma
(* -27.0 j)
k
(fma (fma z (* y (* 18.0 x)) (* -4.0 a)) t (fma c b (* (* -4.0 x) i))))
(* (fma y (* z (* t 18.0)) (* i -4.0)) x)))assert(x < y && y < z && z < t && t < a && a < b && b < c && c < i && i < j && j < k);
double code(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) {
double tmp;
if (((((((((x * 18.0) * y) * z) * t) - ((a * 4.0) * t)) + (b * c)) - ((x * 4.0) * i)) - ((j * 27.0) * k)) <= ((double) INFINITY)) {
tmp = fma((-27.0 * j), k, fma(fma(z, (y * (18.0 * x)), (-4.0 * a)), t, fma(c, b, ((-4.0 * x) * i))));
} else {
tmp = fma(y, (z * (t * 18.0)), (i * -4.0)) * x;
}
return tmp;
}
x, y, z, t, a, b, c, i, j, k = sort([x, y, z, t, a, b, c, i, j, k]) function code(x, y, z, t, a, b, c, i, j, k) tmp = 0.0 if (Float64(Float64(Float64(Float64(Float64(Float64(Float64(Float64(x * 18.0) * y) * z) * t) - Float64(Float64(a * 4.0) * t)) + Float64(b * c)) - Float64(Float64(x * 4.0) * i)) - Float64(Float64(j * 27.0) * k)) <= Inf) tmp = fma(Float64(-27.0 * j), k, fma(fma(z, Float64(y * Float64(18.0 * x)), Float64(-4.0 * a)), t, fma(c, b, Float64(Float64(-4.0 * x) * i)))); else tmp = Float64(fma(y, Float64(z * Float64(t * 18.0)), Float64(i * -4.0)) * x); end return tmp end
NOTE: x, y, z, t, a, b, c, i, j, and k should be sorted in increasing order before calling this function. code[x_, y_, z_, t_, a_, b_, c_, i_, j_, k_] := If[LessEqual[N[(N[(N[(N[(N[(N[(N[(N[(x * 18.0), $MachinePrecision] * y), $MachinePrecision] * z), $MachinePrecision] * t), $MachinePrecision] - N[(N[(a * 4.0), $MachinePrecision] * t), $MachinePrecision]), $MachinePrecision] + N[(b * c), $MachinePrecision]), $MachinePrecision] - N[(N[(x * 4.0), $MachinePrecision] * i), $MachinePrecision]), $MachinePrecision] - N[(N[(j * 27.0), $MachinePrecision] * k), $MachinePrecision]), $MachinePrecision], Infinity], N[(N[(-27.0 * j), $MachinePrecision] * k + N[(N[(z * N[(y * N[(18.0 * x), $MachinePrecision]), $MachinePrecision] + N[(-4.0 * a), $MachinePrecision]), $MachinePrecision] * t + N[(c * b + N[(N[(-4.0 * x), $MachinePrecision] * i), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(y * N[(z * N[(t * 18.0), $MachinePrecision]), $MachinePrecision] + N[(i * -4.0), $MachinePrecision]), $MachinePrecision] * x), $MachinePrecision]]
\begin{array}{l}
[x, y, z, t, a, b, c, i, j, k] = \mathsf{sort}([x, y, z, t, a, b, c, i, j, k])\\
\\
\begin{array}{l}
\mathbf{if}\;\left(\left(\left(\left(\left(\left(x \cdot 18\right) \cdot y\right) \cdot z\right) \cdot t - \left(a \cdot 4\right) \cdot t\right) + b \cdot c\right) - \left(x \cdot 4\right) \cdot i\right) - \left(j \cdot 27\right) \cdot k \leq \infty:\\
\;\;\;\;\mathsf{fma}\left(-27 \cdot j, k, \mathsf{fma}\left(\mathsf{fma}\left(z, y \cdot \left(18 \cdot x\right), -4 \cdot a\right), t, \mathsf{fma}\left(c, b, \left(-4 \cdot x\right) \cdot i\right)\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(y, z \cdot \left(t \cdot 18\right), i \cdot -4\right) \cdot x\\
\end{array}
\end{array}
if (-.f64 (-.f64 (+.f64 (-.f64 (*.f64 (*.f64 (*.f64 (*.f64 x #s(literal 18 binary64)) y) z) t) (*.f64 (*.f64 a #s(literal 4 binary64)) t)) (*.f64 b c)) (*.f64 (*.f64 x #s(literal 4 binary64)) i)) (*.f64 (*.f64 j #s(literal 27 binary64)) k)) < +inf.0Initial program 95.6%
Applied rewrites95.6%
if +inf.0 < (-.f64 (-.f64 (+.f64 (-.f64 (*.f64 (*.f64 (*.f64 (*.f64 x #s(literal 18 binary64)) y) z) t) (*.f64 (*.f64 a #s(literal 4 binary64)) t)) (*.f64 b c)) (*.f64 (*.f64 x #s(literal 4 binary64)) i)) (*.f64 (*.f64 j #s(literal 27 binary64)) k)) Initial program 0.0%
Taylor expanded in x around inf
*-commutativeN/A
lower-*.f64N/A
fp-cancel-sub-sign-invN/A
metadata-evalN/A
+-commutativeN/A
lower-fma.f64N/A
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f6474.5
Applied rewrites74.5%
Applied rewrites74.5%
NOTE: x, y, z, t, a, b, c, i, j, and k should be sorted in increasing order before calling this function.
(FPCore (x y z t a b c i j k)
:precision binary64
(if (<=
(-
(-
(+ (- (* (* (* (* x 18.0) y) z) t) (* (* a 4.0) t)) (* b c))
(* (* x 4.0) i))
(* (* j 27.0) k))
INFINITY)
(fma -4.0 (fma t a (* i x)) (fma -27.0 (* k j) (* b c)))
(* (fma y (* z (* t 18.0)) (* i -4.0)) x)))assert(x < y && y < z && z < t && t < a && a < b && b < c && c < i && i < j && j < k);
double code(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) {
double tmp;
if (((((((((x * 18.0) * y) * z) * t) - ((a * 4.0) * t)) + (b * c)) - ((x * 4.0) * i)) - ((j * 27.0) * k)) <= ((double) INFINITY)) {
tmp = fma(-4.0, fma(t, a, (i * x)), fma(-27.0, (k * j), (b * c)));
} else {
tmp = fma(y, (z * (t * 18.0)), (i * -4.0)) * x;
}
return tmp;
}
x, y, z, t, a, b, c, i, j, k = sort([x, y, z, t, a, b, c, i, j, k]) function code(x, y, z, t, a, b, c, i, j, k) tmp = 0.0 if (Float64(Float64(Float64(Float64(Float64(Float64(Float64(Float64(x * 18.0) * y) * z) * t) - Float64(Float64(a * 4.0) * t)) + Float64(b * c)) - Float64(Float64(x * 4.0) * i)) - Float64(Float64(j * 27.0) * k)) <= Inf) tmp = fma(-4.0, fma(t, a, Float64(i * x)), fma(-27.0, Float64(k * j), Float64(b * c))); else tmp = Float64(fma(y, Float64(z * Float64(t * 18.0)), Float64(i * -4.0)) * x); end return tmp end
NOTE: x, y, z, t, a, b, c, i, j, and k should be sorted in increasing order before calling this function. code[x_, y_, z_, t_, a_, b_, c_, i_, j_, k_] := If[LessEqual[N[(N[(N[(N[(N[(N[(N[(N[(x * 18.0), $MachinePrecision] * y), $MachinePrecision] * z), $MachinePrecision] * t), $MachinePrecision] - N[(N[(a * 4.0), $MachinePrecision] * t), $MachinePrecision]), $MachinePrecision] + N[(b * c), $MachinePrecision]), $MachinePrecision] - N[(N[(x * 4.0), $MachinePrecision] * i), $MachinePrecision]), $MachinePrecision] - N[(N[(j * 27.0), $MachinePrecision] * k), $MachinePrecision]), $MachinePrecision], Infinity], N[(-4.0 * N[(t * a + N[(i * x), $MachinePrecision]), $MachinePrecision] + N[(-27.0 * N[(k * j), $MachinePrecision] + N[(b * c), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(y * N[(z * N[(t * 18.0), $MachinePrecision]), $MachinePrecision] + N[(i * -4.0), $MachinePrecision]), $MachinePrecision] * x), $MachinePrecision]]
\begin{array}{l}
[x, y, z, t, a, b, c, i, j, k] = \mathsf{sort}([x, y, z, t, a, b, c, i, j, k])\\
\\
\begin{array}{l}
\mathbf{if}\;\left(\left(\left(\left(\left(\left(x \cdot 18\right) \cdot y\right) \cdot z\right) \cdot t - \left(a \cdot 4\right) \cdot t\right) + b \cdot c\right) - \left(x \cdot 4\right) \cdot i\right) - \left(j \cdot 27\right) \cdot k \leq \infty:\\
\;\;\;\;\mathsf{fma}\left(-4, \mathsf{fma}\left(t, a, i \cdot x\right), \mathsf{fma}\left(-27, k \cdot j, b \cdot c\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(y, z \cdot \left(t \cdot 18\right), i \cdot -4\right) \cdot x\\
\end{array}
\end{array}
if (-.f64 (-.f64 (+.f64 (-.f64 (*.f64 (*.f64 (*.f64 (*.f64 x #s(literal 18 binary64)) y) z) t) (*.f64 (*.f64 a #s(literal 4 binary64)) t)) (*.f64 b c)) (*.f64 (*.f64 x #s(literal 4 binary64)) i)) (*.f64 (*.f64 j #s(literal 27 binary64)) k)) < +inf.0Initial program 95.6%
Taylor expanded in j around inf
associate-*r*N/A
lower-*.f64N/A
lower-*.f6420.3
Applied rewrites20.3%
Taylor expanded in y around 0
associate-+r+N/A
associate--r+N/A
fp-cancel-sign-sub-invN/A
metadata-evalN/A
associate--r-N/A
fp-cancel-sub-sign-invN/A
metadata-evalN/A
+-commutativeN/A
associate-+r+N/A
+-commutativeN/A
associate-+r+N/A
Applied rewrites89.4%
if +inf.0 < (-.f64 (-.f64 (+.f64 (-.f64 (*.f64 (*.f64 (*.f64 (*.f64 x #s(literal 18 binary64)) y) z) t) (*.f64 (*.f64 a #s(literal 4 binary64)) t)) (*.f64 b c)) (*.f64 (*.f64 x #s(literal 4 binary64)) i)) (*.f64 (*.f64 j #s(literal 27 binary64)) k)) Initial program 0.0%
Taylor expanded in x around inf
*-commutativeN/A
lower-*.f64N/A
fp-cancel-sub-sign-invN/A
metadata-evalN/A
+-commutativeN/A
lower-fma.f64N/A
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f6474.5
Applied rewrites74.5%
Applied rewrites74.5%
NOTE: x, y, z, t, a, b, c, i, j, and k should be sorted in increasing order before calling this function.
(FPCore (x y z t a b c i j k)
:precision binary64
(let* ((t_1 (* (* -4.0 a) t)))
(if (<= (* b c) -5e+60)
(* c b)
(if (<= (* b c) -5e-302)
(* (* k j) -27.0)
(if (<= (* b c) 1e-297)
t_1
(if (<= (* b c) 1e-130)
(* (* -4.0 x) i)
(if (<= (* b c) 2e+116) t_1 (* c b))))))))assert(x < y && y < z && z < t && t < a && a < b && b < c && c < i && i < j && j < k);
double code(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) {
double t_1 = (-4.0 * a) * t;
double tmp;
if ((b * c) <= -5e+60) {
tmp = c * b;
} else if ((b * c) <= -5e-302) {
tmp = (k * j) * -27.0;
} else if ((b * c) <= 1e-297) {
tmp = t_1;
} else if ((b * c) <= 1e-130) {
tmp = (-4.0 * x) * i;
} else if ((b * c) <= 2e+116) {
tmp = t_1;
} else {
tmp = c * b;
}
return tmp;
}
NOTE: x, y, z, t, a, b, c, i, j, and k should be sorted in increasing order before calling this function.
module fmin_fmax_functions
implicit none
private
public fmax
public fmin
interface fmax
module procedure fmax88
module procedure fmax44
module procedure fmax84
module procedure fmax48
end interface
interface fmin
module procedure fmin88
module procedure fmin44
module procedure fmin84
module procedure fmin48
end interface
contains
real(8) function fmax88(x, y) result (res)
real(8), intent (in) :: x
real(8), intent (in) :: y
res = merge(y, merge(x, max(x, y), y /= y), x /= x)
end function
real(4) function fmax44(x, y) result (res)
real(4), intent (in) :: x
real(4), intent (in) :: y
res = merge(y, merge(x, max(x, y), y /= y), x /= x)
end function
real(8) function fmax84(x, y) result(res)
real(8), intent (in) :: x
real(4), intent (in) :: y
res = merge(dble(y), merge(x, max(x, dble(y)), y /= y), x /= x)
end function
real(8) function fmax48(x, y) result(res)
real(4), intent (in) :: x
real(8), intent (in) :: y
res = merge(y, merge(dble(x), max(dble(x), y), y /= y), x /= x)
end function
real(8) function fmin88(x, y) result (res)
real(8), intent (in) :: x
real(8), intent (in) :: y
res = merge(y, merge(x, min(x, y), y /= y), x /= x)
end function
real(4) function fmin44(x, y) result (res)
real(4), intent (in) :: x
real(4), intent (in) :: y
res = merge(y, merge(x, min(x, y), y /= y), x /= x)
end function
real(8) function fmin84(x, y) result(res)
real(8), intent (in) :: x
real(4), intent (in) :: y
res = merge(dble(y), merge(x, min(x, dble(y)), y /= y), x /= x)
end function
real(8) function fmin48(x, y) result(res)
real(4), intent (in) :: x
real(8), intent (in) :: y
res = merge(y, merge(dble(x), min(dble(x), y), y /= y), x /= x)
end function
end module
real(8) function code(x, y, z, t, a, b, c, i, j, k)
use fmin_fmax_functions
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), intent (in) :: b
real(8), intent (in) :: c
real(8), intent (in) :: i
real(8), intent (in) :: j
real(8), intent (in) :: k
real(8) :: t_1
real(8) :: tmp
t_1 = ((-4.0d0) * a) * t
if ((b * c) <= (-5d+60)) then
tmp = c * b
else if ((b * c) <= (-5d-302)) then
tmp = (k * j) * (-27.0d0)
else if ((b * c) <= 1d-297) then
tmp = t_1
else if ((b * c) <= 1d-130) then
tmp = ((-4.0d0) * x) * i
else if ((b * c) <= 2d+116) then
tmp = t_1
else
tmp = c * b
end if
code = tmp
end function
assert x < y && y < z && z < t && t < a && a < b && b < c && c < i && i < j && j < k;
public static double code(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) {
double t_1 = (-4.0 * a) * t;
double tmp;
if ((b * c) <= -5e+60) {
tmp = c * b;
} else if ((b * c) <= -5e-302) {
tmp = (k * j) * -27.0;
} else if ((b * c) <= 1e-297) {
tmp = t_1;
} else if ((b * c) <= 1e-130) {
tmp = (-4.0 * x) * i;
} else if ((b * c) <= 2e+116) {
tmp = t_1;
} else {
tmp = c * b;
}
return tmp;
}
[x, y, z, t, a, b, c, i, j, k] = sort([x, y, z, t, a, b, c, i, j, k]) def code(x, y, z, t, a, b, c, i, j, k): t_1 = (-4.0 * a) * t tmp = 0 if (b * c) <= -5e+60: tmp = c * b elif (b * c) <= -5e-302: tmp = (k * j) * -27.0 elif (b * c) <= 1e-297: tmp = t_1 elif (b * c) <= 1e-130: tmp = (-4.0 * x) * i elif (b * c) <= 2e+116: tmp = t_1 else: tmp = c * b return tmp
x, y, z, t, a, b, c, i, j, k = sort([x, y, z, t, a, b, c, i, j, k]) function code(x, y, z, t, a, b, c, i, j, k) t_1 = Float64(Float64(-4.0 * a) * t) tmp = 0.0 if (Float64(b * c) <= -5e+60) tmp = Float64(c * b); elseif (Float64(b * c) <= -5e-302) tmp = Float64(Float64(k * j) * -27.0); elseif (Float64(b * c) <= 1e-297) tmp = t_1; elseif (Float64(b * c) <= 1e-130) tmp = Float64(Float64(-4.0 * x) * i); elseif (Float64(b * c) <= 2e+116) tmp = t_1; else tmp = Float64(c * b); end return tmp end
x, y, z, t, a, b, c, i, j, k = num2cell(sort([x, y, z, t, a, b, c, i, j, k])){:}
function tmp_2 = code(x, y, z, t, a, b, c, i, j, k)
t_1 = (-4.0 * a) * t;
tmp = 0.0;
if ((b * c) <= -5e+60)
tmp = c * b;
elseif ((b * c) <= -5e-302)
tmp = (k * j) * -27.0;
elseif ((b * c) <= 1e-297)
tmp = t_1;
elseif ((b * c) <= 1e-130)
tmp = (-4.0 * x) * i;
elseif ((b * c) <= 2e+116)
tmp = t_1;
else
tmp = c * b;
end
tmp_2 = tmp;
end
NOTE: x, y, z, t, a, b, c, i, j, and k should be sorted in increasing order before calling this function.
code[x_, y_, z_, t_, a_, b_, c_, i_, j_, k_] := Block[{t$95$1 = N[(N[(-4.0 * a), $MachinePrecision] * t), $MachinePrecision]}, If[LessEqual[N[(b * c), $MachinePrecision], -5e+60], N[(c * b), $MachinePrecision], If[LessEqual[N[(b * c), $MachinePrecision], -5e-302], N[(N[(k * j), $MachinePrecision] * -27.0), $MachinePrecision], If[LessEqual[N[(b * c), $MachinePrecision], 1e-297], t$95$1, If[LessEqual[N[(b * c), $MachinePrecision], 1e-130], N[(N[(-4.0 * x), $MachinePrecision] * i), $MachinePrecision], If[LessEqual[N[(b * c), $MachinePrecision], 2e+116], t$95$1, N[(c * b), $MachinePrecision]]]]]]]
\begin{array}{l}
[x, y, z, t, a, b, c, i, j, k] = \mathsf{sort}([x, y, z, t, a, b, c, i, j, k])\\
\\
\begin{array}{l}
t_1 := \left(-4 \cdot a\right) \cdot t\\
\mathbf{if}\;b \cdot c \leq -5 \cdot 10^{+60}:\\
\;\;\;\;c \cdot b\\
\mathbf{elif}\;b \cdot c \leq -5 \cdot 10^{-302}:\\
\;\;\;\;\left(k \cdot j\right) \cdot -27\\
\mathbf{elif}\;b \cdot c \leq 10^{-297}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;b \cdot c \leq 10^{-130}:\\
\;\;\;\;\left(-4 \cdot x\right) \cdot i\\
\mathbf{elif}\;b \cdot c \leq 2 \cdot 10^{+116}:\\
\;\;\;\;t\_1\\
\mathbf{else}:\\
\;\;\;\;c \cdot b\\
\end{array}
\end{array}
if (*.f64 b c) < -4.99999999999999975e60 or 2.00000000000000003e116 < (*.f64 b c) Initial program 84.0%
Taylor expanded in y around 0
distribute-lft-outN/A
fp-cancel-sub-sign-invN/A
*-commutativeN/A
lower-fma.f64N/A
metadata-evalN/A
lower-*.f64N/A
+-commutativeN/A
lower-fma.f64N/A
lower-*.f6489.1
Applied rewrites89.1%
Taylor expanded in a around 0
associate--r+N/A
fp-cancel-sub-sign-invN/A
metadata-evalN/A
+-commutativeN/A
associate-+r+N/A
associate--l+N/A
associate-*r*N/A
associate-*r*N/A
*-commutativeN/A
associate-*r*N/A
associate-*r*N/A
distribute-rgt-inN/A
*-commutativeN/A
Applied rewrites77.8%
Taylor expanded in b around inf
*-commutativeN/A
lower-*.f6460.3
Applied rewrites60.3%
if -4.99999999999999975e60 < (*.f64 b c) < -5.00000000000000033e-302Initial program 84.4%
Taylor expanded in j around inf
associate-*r*N/A
lower-*.f64N/A
lower-*.f6437.4
Applied rewrites37.4%
Applied rewrites37.5%
if -5.00000000000000033e-302 < (*.f64 b c) < 1.00000000000000004e-297 or 1.0000000000000001e-130 < (*.f64 b c) < 2.00000000000000003e116Initial program 84.0%
lift--.f64N/A
lift-*.f64N/A
fp-cancel-sub-sign-invN/A
+-commutativeN/A
lower-fma.f64N/A
lift-*.f64N/A
*-commutativeN/A
distribute-lft-neg-inN/A
lower-*.f64N/A
metadata-eval87.5
lift-+.f64N/A
lift--.f64N/A
lift-*.f64N/A
lift-*.f64N/A
distribute-rgt-out--N/A
*-commutativeN/A
Applied rewrites88.7%
Taylor expanded in a around inf
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f6437.6
Applied rewrites37.6%
Applied rewrites37.6%
if 1.00000000000000004e-297 < (*.f64 b c) < 1.0000000000000001e-130Initial program 83.3%
Taylor expanded in i around inf
*-commutativeN/A
associate-*r*N/A
lower-*.f64N/A
lower-*.f6469.9
Applied rewrites69.9%
NOTE: x, y, z, t, a, b, c, i, j, and k should be sorted in increasing order before calling this function.
(FPCore (x y z t a b c i j k)
:precision binary64
(let* ((t_1 (* (* -4.0 a) t)))
(if (<= (* b c) -5e+60)
(* c b)
(if (<= (* b c) -5e-302)
(* (* k -27.0) j)
(if (<= (* b c) 1e-297)
t_1
(if (<= (* b c) 1e-130)
(* (* -4.0 x) i)
(if (<= (* b c) 2e+116) t_1 (* c b))))))))assert(x < y && y < z && z < t && t < a && a < b && b < c && c < i && i < j && j < k);
double code(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) {
double t_1 = (-4.0 * a) * t;
double tmp;
if ((b * c) <= -5e+60) {
tmp = c * b;
} else if ((b * c) <= -5e-302) {
tmp = (k * -27.0) * j;
} else if ((b * c) <= 1e-297) {
tmp = t_1;
} else if ((b * c) <= 1e-130) {
tmp = (-4.0 * x) * i;
} else if ((b * c) <= 2e+116) {
tmp = t_1;
} else {
tmp = c * b;
}
return tmp;
}
NOTE: x, y, z, t, a, b, c, i, j, and k should be sorted in increasing order before calling this function.
module fmin_fmax_functions
implicit none
private
public fmax
public fmin
interface fmax
module procedure fmax88
module procedure fmax44
module procedure fmax84
module procedure fmax48
end interface
interface fmin
module procedure fmin88
module procedure fmin44
module procedure fmin84
module procedure fmin48
end interface
contains
real(8) function fmax88(x, y) result (res)
real(8), intent (in) :: x
real(8), intent (in) :: y
res = merge(y, merge(x, max(x, y), y /= y), x /= x)
end function
real(4) function fmax44(x, y) result (res)
real(4), intent (in) :: x
real(4), intent (in) :: y
res = merge(y, merge(x, max(x, y), y /= y), x /= x)
end function
real(8) function fmax84(x, y) result(res)
real(8), intent (in) :: x
real(4), intent (in) :: y
res = merge(dble(y), merge(x, max(x, dble(y)), y /= y), x /= x)
end function
real(8) function fmax48(x, y) result(res)
real(4), intent (in) :: x
real(8), intent (in) :: y
res = merge(y, merge(dble(x), max(dble(x), y), y /= y), x /= x)
end function
real(8) function fmin88(x, y) result (res)
real(8), intent (in) :: x
real(8), intent (in) :: y
res = merge(y, merge(x, min(x, y), y /= y), x /= x)
end function
real(4) function fmin44(x, y) result (res)
real(4), intent (in) :: x
real(4), intent (in) :: y
res = merge(y, merge(x, min(x, y), y /= y), x /= x)
end function
real(8) function fmin84(x, y) result(res)
real(8), intent (in) :: x
real(4), intent (in) :: y
res = merge(dble(y), merge(x, min(x, dble(y)), y /= y), x /= x)
end function
real(8) function fmin48(x, y) result(res)
real(4), intent (in) :: x
real(8), intent (in) :: y
res = merge(y, merge(dble(x), min(dble(x), y), y /= y), x /= x)
end function
end module
real(8) function code(x, y, z, t, a, b, c, i, j, k)
use fmin_fmax_functions
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), intent (in) :: b
real(8), intent (in) :: c
real(8), intent (in) :: i
real(8), intent (in) :: j
real(8), intent (in) :: k
real(8) :: t_1
real(8) :: tmp
t_1 = ((-4.0d0) * a) * t
if ((b * c) <= (-5d+60)) then
tmp = c * b
else if ((b * c) <= (-5d-302)) then
tmp = (k * (-27.0d0)) * j
else if ((b * c) <= 1d-297) then
tmp = t_1
else if ((b * c) <= 1d-130) then
tmp = ((-4.0d0) * x) * i
else if ((b * c) <= 2d+116) then
tmp = t_1
else
tmp = c * b
end if
code = tmp
end function
assert x < y && y < z && z < t && t < a && a < b && b < c && c < i && i < j && j < k;
public static double code(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) {
double t_1 = (-4.0 * a) * t;
double tmp;
if ((b * c) <= -5e+60) {
tmp = c * b;
} else if ((b * c) <= -5e-302) {
tmp = (k * -27.0) * j;
} else if ((b * c) <= 1e-297) {
tmp = t_1;
} else if ((b * c) <= 1e-130) {
tmp = (-4.0 * x) * i;
} else if ((b * c) <= 2e+116) {
tmp = t_1;
} else {
tmp = c * b;
}
return tmp;
}
[x, y, z, t, a, b, c, i, j, k] = sort([x, y, z, t, a, b, c, i, j, k]) def code(x, y, z, t, a, b, c, i, j, k): t_1 = (-4.0 * a) * t tmp = 0 if (b * c) <= -5e+60: tmp = c * b elif (b * c) <= -5e-302: tmp = (k * -27.0) * j elif (b * c) <= 1e-297: tmp = t_1 elif (b * c) <= 1e-130: tmp = (-4.0 * x) * i elif (b * c) <= 2e+116: tmp = t_1 else: tmp = c * b return tmp
x, y, z, t, a, b, c, i, j, k = sort([x, y, z, t, a, b, c, i, j, k]) function code(x, y, z, t, a, b, c, i, j, k) t_1 = Float64(Float64(-4.0 * a) * t) tmp = 0.0 if (Float64(b * c) <= -5e+60) tmp = Float64(c * b); elseif (Float64(b * c) <= -5e-302) tmp = Float64(Float64(k * -27.0) * j); elseif (Float64(b * c) <= 1e-297) tmp = t_1; elseif (Float64(b * c) <= 1e-130) tmp = Float64(Float64(-4.0 * x) * i); elseif (Float64(b * c) <= 2e+116) tmp = t_1; else tmp = Float64(c * b); end return tmp end
x, y, z, t, a, b, c, i, j, k = num2cell(sort([x, y, z, t, a, b, c, i, j, k])){:}
function tmp_2 = code(x, y, z, t, a, b, c, i, j, k)
t_1 = (-4.0 * a) * t;
tmp = 0.0;
if ((b * c) <= -5e+60)
tmp = c * b;
elseif ((b * c) <= -5e-302)
tmp = (k * -27.0) * j;
elseif ((b * c) <= 1e-297)
tmp = t_1;
elseif ((b * c) <= 1e-130)
tmp = (-4.0 * x) * i;
elseif ((b * c) <= 2e+116)
tmp = t_1;
else
tmp = c * b;
end
tmp_2 = tmp;
end
NOTE: x, y, z, t, a, b, c, i, j, and k should be sorted in increasing order before calling this function.
code[x_, y_, z_, t_, a_, b_, c_, i_, j_, k_] := Block[{t$95$1 = N[(N[(-4.0 * a), $MachinePrecision] * t), $MachinePrecision]}, If[LessEqual[N[(b * c), $MachinePrecision], -5e+60], N[(c * b), $MachinePrecision], If[LessEqual[N[(b * c), $MachinePrecision], -5e-302], N[(N[(k * -27.0), $MachinePrecision] * j), $MachinePrecision], If[LessEqual[N[(b * c), $MachinePrecision], 1e-297], t$95$1, If[LessEqual[N[(b * c), $MachinePrecision], 1e-130], N[(N[(-4.0 * x), $MachinePrecision] * i), $MachinePrecision], If[LessEqual[N[(b * c), $MachinePrecision], 2e+116], t$95$1, N[(c * b), $MachinePrecision]]]]]]]
\begin{array}{l}
[x, y, z, t, a, b, c, i, j, k] = \mathsf{sort}([x, y, z, t, a, b, c, i, j, k])\\
\\
\begin{array}{l}
t_1 := \left(-4 \cdot a\right) \cdot t\\
\mathbf{if}\;b \cdot c \leq -5 \cdot 10^{+60}:\\
\;\;\;\;c \cdot b\\
\mathbf{elif}\;b \cdot c \leq -5 \cdot 10^{-302}:\\
\;\;\;\;\left(k \cdot -27\right) \cdot j\\
\mathbf{elif}\;b \cdot c \leq 10^{-297}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;b \cdot c \leq 10^{-130}:\\
\;\;\;\;\left(-4 \cdot x\right) \cdot i\\
\mathbf{elif}\;b \cdot c \leq 2 \cdot 10^{+116}:\\
\;\;\;\;t\_1\\
\mathbf{else}:\\
\;\;\;\;c \cdot b\\
\end{array}
\end{array}
if (*.f64 b c) < -4.99999999999999975e60 or 2.00000000000000003e116 < (*.f64 b c) Initial program 84.0%
Taylor expanded in y around 0
distribute-lft-outN/A
fp-cancel-sub-sign-invN/A
*-commutativeN/A
lower-fma.f64N/A
metadata-evalN/A
lower-*.f64N/A
+-commutativeN/A
lower-fma.f64N/A
lower-*.f6489.1
Applied rewrites89.1%
Taylor expanded in a around 0
associate--r+N/A
fp-cancel-sub-sign-invN/A
metadata-evalN/A
+-commutativeN/A
associate-+r+N/A
associate--l+N/A
associate-*r*N/A
associate-*r*N/A
*-commutativeN/A
associate-*r*N/A
associate-*r*N/A
distribute-rgt-inN/A
*-commutativeN/A
Applied rewrites77.8%
Taylor expanded in b around inf
*-commutativeN/A
lower-*.f6460.3
Applied rewrites60.3%
if -4.99999999999999975e60 < (*.f64 b c) < -5.00000000000000033e-302Initial program 84.4%
Taylor expanded in j around inf
associate-*r*N/A
lower-*.f64N/A
lower-*.f6437.4
Applied rewrites37.4%
Applied rewrites37.4%
if -5.00000000000000033e-302 < (*.f64 b c) < 1.00000000000000004e-297 or 1.0000000000000001e-130 < (*.f64 b c) < 2.00000000000000003e116Initial program 84.0%
lift--.f64N/A
lift-*.f64N/A
fp-cancel-sub-sign-invN/A
+-commutativeN/A
lower-fma.f64N/A
lift-*.f64N/A
*-commutativeN/A
distribute-lft-neg-inN/A
lower-*.f64N/A
metadata-eval87.5
lift-+.f64N/A
lift--.f64N/A
lift-*.f64N/A
lift-*.f64N/A
distribute-rgt-out--N/A
*-commutativeN/A
Applied rewrites88.7%
Taylor expanded in a around inf
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f6437.6
Applied rewrites37.6%
Applied rewrites37.6%
if 1.00000000000000004e-297 < (*.f64 b c) < 1.0000000000000001e-130Initial program 83.3%
Taylor expanded in i around inf
*-commutativeN/A
associate-*r*N/A
lower-*.f64N/A
lower-*.f6469.9
Applied rewrites69.9%
NOTE: x, y, z, t, a, b, c, i, j, and k should be sorted in increasing order before calling this function.
(FPCore (x y z t a b c i j k)
:precision binary64
(let* ((t_1 (* (* -4.0 a) t)))
(if (<= (* b c) -5e+60)
(* c b)
(if (<= (* b c) -5e-302)
(* (* -27.0 j) k)
(if (<= (* b c) 1e-297)
t_1
(if (<= (* b c) 1e-130)
(* (* -4.0 x) i)
(if (<= (* b c) 2e+116) t_1 (* c b))))))))assert(x < y && y < z && z < t && t < a && a < b && b < c && c < i && i < j && j < k);
double code(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) {
double t_1 = (-4.0 * a) * t;
double tmp;
if ((b * c) <= -5e+60) {
tmp = c * b;
} else if ((b * c) <= -5e-302) {
tmp = (-27.0 * j) * k;
} else if ((b * c) <= 1e-297) {
tmp = t_1;
} else if ((b * c) <= 1e-130) {
tmp = (-4.0 * x) * i;
} else if ((b * c) <= 2e+116) {
tmp = t_1;
} else {
tmp = c * b;
}
return tmp;
}
NOTE: x, y, z, t, a, b, c, i, j, and k should be sorted in increasing order before calling this function.
module fmin_fmax_functions
implicit none
private
public fmax
public fmin
interface fmax
module procedure fmax88
module procedure fmax44
module procedure fmax84
module procedure fmax48
end interface
interface fmin
module procedure fmin88
module procedure fmin44
module procedure fmin84
module procedure fmin48
end interface
contains
real(8) function fmax88(x, y) result (res)
real(8), intent (in) :: x
real(8), intent (in) :: y
res = merge(y, merge(x, max(x, y), y /= y), x /= x)
end function
real(4) function fmax44(x, y) result (res)
real(4), intent (in) :: x
real(4), intent (in) :: y
res = merge(y, merge(x, max(x, y), y /= y), x /= x)
end function
real(8) function fmax84(x, y) result(res)
real(8), intent (in) :: x
real(4), intent (in) :: y
res = merge(dble(y), merge(x, max(x, dble(y)), y /= y), x /= x)
end function
real(8) function fmax48(x, y) result(res)
real(4), intent (in) :: x
real(8), intent (in) :: y
res = merge(y, merge(dble(x), max(dble(x), y), y /= y), x /= x)
end function
real(8) function fmin88(x, y) result (res)
real(8), intent (in) :: x
real(8), intent (in) :: y
res = merge(y, merge(x, min(x, y), y /= y), x /= x)
end function
real(4) function fmin44(x, y) result (res)
real(4), intent (in) :: x
real(4), intent (in) :: y
res = merge(y, merge(x, min(x, y), y /= y), x /= x)
end function
real(8) function fmin84(x, y) result(res)
real(8), intent (in) :: x
real(4), intent (in) :: y
res = merge(dble(y), merge(x, min(x, dble(y)), y /= y), x /= x)
end function
real(8) function fmin48(x, y) result(res)
real(4), intent (in) :: x
real(8), intent (in) :: y
res = merge(y, merge(dble(x), min(dble(x), y), y /= y), x /= x)
end function
end module
real(8) function code(x, y, z, t, a, b, c, i, j, k)
use fmin_fmax_functions
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), intent (in) :: b
real(8), intent (in) :: c
real(8), intent (in) :: i
real(8), intent (in) :: j
real(8), intent (in) :: k
real(8) :: t_1
real(8) :: tmp
t_1 = ((-4.0d0) * a) * t
if ((b * c) <= (-5d+60)) then
tmp = c * b
else if ((b * c) <= (-5d-302)) then
tmp = ((-27.0d0) * j) * k
else if ((b * c) <= 1d-297) then
tmp = t_1
else if ((b * c) <= 1d-130) then
tmp = ((-4.0d0) * x) * i
else if ((b * c) <= 2d+116) then
tmp = t_1
else
tmp = c * b
end if
code = tmp
end function
assert x < y && y < z && z < t && t < a && a < b && b < c && c < i && i < j && j < k;
public static double code(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) {
double t_1 = (-4.0 * a) * t;
double tmp;
if ((b * c) <= -5e+60) {
tmp = c * b;
} else if ((b * c) <= -5e-302) {
tmp = (-27.0 * j) * k;
} else if ((b * c) <= 1e-297) {
tmp = t_1;
} else if ((b * c) <= 1e-130) {
tmp = (-4.0 * x) * i;
} else if ((b * c) <= 2e+116) {
tmp = t_1;
} else {
tmp = c * b;
}
return tmp;
}
[x, y, z, t, a, b, c, i, j, k] = sort([x, y, z, t, a, b, c, i, j, k]) def code(x, y, z, t, a, b, c, i, j, k): t_1 = (-4.0 * a) * t tmp = 0 if (b * c) <= -5e+60: tmp = c * b elif (b * c) <= -5e-302: tmp = (-27.0 * j) * k elif (b * c) <= 1e-297: tmp = t_1 elif (b * c) <= 1e-130: tmp = (-4.0 * x) * i elif (b * c) <= 2e+116: tmp = t_1 else: tmp = c * b return tmp
x, y, z, t, a, b, c, i, j, k = sort([x, y, z, t, a, b, c, i, j, k]) function code(x, y, z, t, a, b, c, i, j, k) t_1 = Float64(Float64(-4.0 * a) * t) tmp = 0.0 if (Float64(b * c) <= -5e+60) tmp = Float64(c * b); elseif (Float64(b * c) <= -5e-302) tmp = Float64(Float64(-27.0 * j) * k); elseif (Float64(b * c) <= 1e-297) tmp = t_1; elseif (Float64(b * c) <= 1e-130) tmp = Float64(Float64(-4.0 * x) * i); elseif (Float64(b * c) <= 2e+116) tmp = t_1; else tmp = Float64(c * b); end return tmp end
x, y, z, t, a, b, c, i, j, k = num2cell(sort([x, y, z, t, a, b, c, i, j, k])){:}
function tmp_2 = code(x, y, z, t, a, b, c, i, j, k)
t_1 = (-4.0 * a) * t;
tmp = 0.0;
if ((b * c) <= -5e+60)
tmp = c * b;
elseif ((b * c) <= -5e-302)
tmp = (-27.0 * j) * k;
elseif ((b * c) <= 1e-297)
tmp = t_1;
elseif ((b * c) <= 1e-130)
tmp = (-4.0 * x) * i;
elseif ((b * c) <= 2e+116)
tmp = t_1;
else
tmp = c * b;
end
tmp_2 = tmp;
end
NOTE: x, y, z, t, a, b, c, i, j, and k should be sorted in increasing order before calling this function.
code[x_, y_, z_, t_, a_, b_, c_, i_, j_, k_] := Block[{t$95$1 = N[(N[(-4.0 * a), $MachinePrecision] * t), $MachinePrecision]}, If[LessEqual[N[(b * c), $MachinePrecision], -5e+60], N[(c * b), $MachinePrecision], If[LessEqual[N[(b * c), $MachinePrecision], -5e-302], N[(N[(-27.0 * j), $MachinePrecision] * k), $MachinePrecision], If[LessEqual[N[(b * c), $MachinePrecision], 1e-297], t$95$1, If[LessEqual[N[(b * c), $MachinePrecision], 1e-130], N[(N[(-4.0 * x), $MachinePrecision] * i), $MachinePrecision], If[LessEqual[N[(b * c), $MachinePrecision], 2e+116], t$95$1, N[(c * b), $MachinePrecision]]]]]]]
\begin{array}{l}
[x, y, z, t, a, b, c, i, j, k] = \mathsf{sort}([x, y, z, t, a, b, c, i, j, k])\\
\\
\begin{array}{l}
t_1 := \left(-4 \cdot a\right) \cdot t\\
\mathbf{if}\;b \cdot c \leq -5 \cdot 10^{+60}:\\
\;\;\;\;c \cdot b\\
\mathbf{elif}\;b \cdot c \leq -5 \cdot 10^{-302}:\\
\;\;\;\;\left(-27 \cdot j\right) \cdot k\\
\mathbf{elif}\;b \cdot c \leq 10^{-297}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;b \cdot c \leq 10^{-130}:\\
\;\;\;\;\left(-4 \cdot x\right) \cdot i\\
\mathbf{elif}\;b \cdot c \leq 2 \cdot 10^{+116}:\\
\;\;\;\;t\_1\\
\mathbf{else}:\\
\;\;\;\;c \cdot b\\
\end{array}
\end{array}
if (*.f64 b c) < -4.99999999999999975e60 or 2.00000000000000003e116 < (*.f64 b c) Initial program 84.0%
Taylor expanded in y around 0
distribute-lft-outN/A
fp-cancel-sub-sign-invN/A
*-commutativeN/A
lower-fma.f64N/A
metadata-evalN/A
lower-*.f64N/A
+-commutativeN/A
lower-fma.f64N/A
lower-*.f6489.1
Applied rewrites89.1%
Taylor expanded in a around 0
associate--r+N/A
fp-cancel-sub-sign-invN/A
metadata-evalN/A
+-commutativeN/A
associate-+r+N/A
associate--l+N/A
associate-*r*N/A
associate-*r*N/A
*-commutativeN/A
associate-*r*N/A
associate-*r*N/A
distribute-rgt-inN/A
*-commutativeN/A
Applied rewrites77.8%
Taylor expanded in b around inf
*-commutativeN/A
lower-*.f6460.3
Applied rewrites60.3%
if -4.99999999999999975e60 < (*.f64 b c) < -5.00000000000000033e-302Initial program 84.4%
Taylor expanded in j around inf
associate-*r*N/A
lower-*.f64N/A
lower-*.f6437.4
Applied rewrites37.4%
if -5.00000000000000033e-302 < (*.f64 b c) < 1.00000000000000004e-297 or 1.0000000000000001e-130 < (*.f64 b c) < 2.00000000000000003e116Initial program 84.0%
lift--.f64N/A
lift-*.f64N/A
fp-cancel-sub-sign-invN/A
+-commutativeN/A
lower-fma.f64N/A
lift-*.f64N/A
*-commutativeN/A
distribute-lft-neg-inN/A
lower-*.f64N/A
metadata-eval87.5
lift-+.f64N/A
lift--.f64N/A
lift-*.f64N/A
lift-*.f64N/A
distribute-rgt-out--N/A
*-commutativeN/A
Applied rewrites88.7%
Taylor expanded in a around inf
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f6437.6
Applied rewrites37.6%
Applied rewrites37.6%
if 1.00000000000000004e-297 < (*.f64 b c) < 1.0000000000000001e-130Initial program 83.3%
Taylor expanded in i around inf
*-commutativeN/A
associate-*r*N/A
lower-*.f64N/A
lower-*.f6469.9
Applied rewrites69.9%
NOTE: x, y, z, t, a, b, c, i, j, and k should be sorted in increasing order before calling this function.
(FPCore (x y z t a b c i j k)
:precision binary64
(let* ((t_1 (* (* j 27.0) k)) (t_2 (fma -4.0 i (* (* (* z y) t) 18.0))))
(if (<= x -5.4e+250)
(* t_2 x)
(if (<= x 1.5e+70)
(- (fma c b (* -4.0 (fma i x (* a t)))) t_1)
(- (fma t_2 x (* c b)) t_1)))))assert(x < y && y < z && z < t && t < a && a < b && b < c && c < i && i < j && j < k);
double code(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) {
double t_1 = (j * 27.0) * k;
double t_2 = fma(-4.0, i, (((z * y) * t) * 18.0));
double tmp;
if (x <= -5.4e+250) {
tmp = t_2 * x;
} else if (x <= 1.5e+70) {
tmp = fma(c, b, (-4.0 * fma(i, x, (a * t)))) - t_1;
} else {
tmp = fma(t_2, x, (c * b)) - t_1;
}
return tmp;
}
x, y, z, t, a, b, c, i, j, k = sort([x, y, z, t, a, b, c, i, j, k]) function code(x, y, z, t, a, b, c, i, j, k) t_1 = Float64(Float64(j * 27.0) * k) t_2 = fma(-4.0, i, Float64(Float64(Float64(z * y) * t) * 18.0)) tmp = 0.0 if (x <= -5.4e+250) tmp = Float64(t_2 * x); elseif (x <= 1.5e+70) tmp = Float64(fma(c, b, Float64(-4.0 * fma(i, x, Float64(a * t)))) - t_1); else tmp = Float64(fma(t_2, x, Float64(c * b)) - t_1); end return tmp end
NOTE: x, y, z, t, a, b, c, i, j, and k should be sorted in increasing order before calling this function.
code[x_, y_, z_, t_, a_, b_, c_, i_, j_, k_] := Block[{t$95$1 = N[(N[(j * 27.0), $MachinePrecision] * k), $MachinePrecision]}, Block[{t$95$2 = N[(-4.0 * i + N[(N[(N[(z * y), $MachinePrecision] * t), $MachinePrecision] * 18.0), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[x, -5.4e+250], N[(t$95$2 * x), $MachinePrecision], If[LessEqual[x, 1.5e+70], N[(N[(c * b + N[(-4.0 * N[(i * x + N[(a * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - t$95$1), $MachinePrecision], N[(N[(t$95$2 * x + N[(c * b), $MachinePrecision]), $MachinePrecision] - t$95$1), $MachinePrecision]]]]]
\begin{array}{l}
[x, y, z, t, a, b, c, i, j, k] = \mathsf{sort}([x, y, z, t, a, b, c, i, j, k])\\
\\
\begin{array}{l}
t_1 := \left(j \cdot 27\right) \cdot k\\
t_2 := \mathsf{fma}\left(-4, i, \left(\left(z \cdot y\right) \cdot t\right) \cdot 18\right)\\
\mathbf{if}\;x \leq -5.4 \cdot 10^{+250}:\\
\;\;\;\;t\_2 \cdot x\\
\mathbf{elif}\;x \leq 1.5 \cdot 10^{+70}:\\
\;\;\;\;\mathsf{fma}\left(c, b, -4 \cdot \mathsf{fma}\left(i, x, a \cdot t\right)\right) - t\_1\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(t\_2, x, c \cdot b\right) - t\_1\\
\end{array}
\end{array}
if x < -5.4e250Initial program 70.6%
Taylor expanded in x around inf
*-commutativeN/A
lower-*.f64N/A
fp-cancel-sub-sign-invN/A
metadata-evalN/A
+-commutativeN/A
lower-fma.f64N/A
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f6488.1
Applied rewrites88.1%
if -5.4e250 < x < 1.49999999999999988e70Initial program 90.4%
Taylor expanded in y around 0
distribute-lft-outN/A
fp-cancel-sub-sign-invN/A
*-commutativeN/A
lower-fma.f64N/A
metadata-evalN/A
lower-*.f64N/A
+-commutativeN/A
lower-fma.f64N/A
lower-*.f6490.4
Applied rewrites90.4%
if 1.49999999999999988e70 < x Initial program 65.1%
Taylor expanded in a around 0
fp-cancel-sub-sign-invN/A
metadata-evalN/A
+-commutativeN/A
associate-+r+N/A
+-commutativeN/A
*-commutativeN/A
associate-*r*N/A
associate-*l*N/A
associate-*r*N/A
distribute-rgt-inN/A
metadata-evalN/A
fp-cancel-sub-sign-invN/A
*-commutativeN/A
lower-fma.f64N/A
Applied rewrites92.3%
NOTE: x, y, z, t, a, b, c, i, j, and k should be sorted in increasing order before calling this function.
(FPCore (x y z t a b c i j k)
:precision binary64
(if (<= x -5.4e+250)
(* (fma -4.0 i (* (* (* z y) t) 18.0)) x)
(if (<= x 1.5e+70)
(- (fma c b (* -4.0 (fma i x (* a t)))) (* (* j 27.0) k))
(fma (fma (* (* y z) t) 18.0 (* -4.0 i)) x (fma -27.0 (* k j) (* b c))))))assert(x < y && y < z && z < t && t < a && a < b && b < c && c < i && i < j && j < k);
double code(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) {
double tmp;
if (x <= -5.4e+250) {
tmp = fma(-4.0, i, (((z * y) * t) * 18.0)) * x;
} else if (x <= 1.5e+70) {
tmp = fma(c, b, (-4.0 * fma(i, x, (a * t)))) - ((j * 27.0) * k);
} else {
tmp = fma(fma(((y * z) * t), 18.0, (-4.0 * i)), x, fma(-27.0, (k * j), (b * c)));
}
return tmp;
}
x, y, z, t, a, b, c, i, j, k = sort([x, y, z, t, a, b, c, i, j, k]) function code(x, y, z, t, a, b, c, i, j, k) tmp = 0.0 if (x <= -5.4e+250) tmp = Float64(fma(-4.0, i, Float64(Float64(Float64(z * y) * t) * 18.0)) * x); elseif (x <= 1.5e+70) tmp = Float64(fma(c, b, Float64(-4.0 * fma(i, x, Float64(a * t)))) - Float64(Float64(j * 27.0) * k)); else tmp = fma(fma(Float64(Float64(y * z) * t), 18.0, Float64(-4.0 * i)), x, fma(-27.0, Float64(k * j), Float64(b * c))); end return tmp end
NOTE: x, y, z, t, a, b, c, i, j, and k should be sorted in increasing order before calling this function. code[x_, y_, z_, t_, a_, b_, c_, i_, j_, k_] := If[LessEqual[x, -5.4e+250], N[(N[(-4.0 * i + N[(N[(N[(z * y), $MachinePrecision] * t), $MachinePrecision] * 18.0), $MachinePrecision]), $MachinePrecision] * x), $MachinePrecision], If[LessEqual[x, 1.5e+70], N[(N[(c * b + N[(-4.0 * N[(i * x + N[(a * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - N[(N[(j * 27.0), $MachinePrecision] * k), $MachinePrecision]), $MachinePrecision], N[(N[(N[(N[(y * z), $MachinePrecision] * t), $MachinePrecision] * 18.0 + N[(-4.0 * i), $MachinePrecision]), $MachinePrecision] * x + N[(-27.0 * N[(k * j), $MachinePrecision] + N[(b * c), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
[x, y, z, t, a, b, c, i, j, k] = \mathsf{sort}([x, y, z, t, a, b, c, i, j, k])\\
\\
\begin{array}{l}
\mathbf{if}\;x \leq -5.4 \cdot 10^{+250}:\\
\;\;\;\;\mathsf{fma}\left(-4, i, \left(\left(z \cdot y\right) \cdot t\right) \cdot 18\right) \cdot x\\
\mathbf{elif}\;x \leq 1.5 \cdot 10^{+70}:\\
\;\;\;\;\mathsf{fma}\left(c, b, -4 \cdot \mathsf{fma}\left(i, x, a \cdot t\right)\right) - \left(j \cdot 27\right) \cdot k\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(\left(y \cdot z\right) \cdot t, 18, -4 \cdot i\right), x, \mathsf{fma}\left(-27, k \cdot j, b \cdot c\right)\right)\\
\end{array}
\end{array}
if x < -5.4e250Initial program 70.6%
Taylor expanded in x around inf
*-commutativeN/A
lower-*.f64N/A
fp-cancel-sub-sign-invN/A
metadata-evalN/A
+-commutativeN/A
lower-fma.f64N/A
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f6488.1
Applied rewrites88.1%
if -5.4e250 < x < 1.49999999999999988e70Initial program 90.4%
Taylor expanded in y around 0
distribute-lft-outN/A
fp-cancel-sub-sign-invN/A
*-commutativeN/A
lower-fma.f64N/A
metadata-evalN/A
lower-*.f64N/A
+-commutativeN/A
lower-fma.f64N/A
lower-*.f6490.4
Applied rewrites90.4%
if 1.49999999999999988e70 < x Initial program 65.1%
Taylor expanded in j around inf
associate-*r*N/A
lower-*.f64N/A
lower-*.f6416.6
Applied rewrites16.6%
Taylor expanded in a around 0
associate--r+N/A
fp-cancel-sub-sign-invN/A
metadata-evalN/A
+-commutativeN/A
associate-+r+N/A
associate--l+N/A
Applied rewrites92.2%
NOTE: x, y, z, t, a, b, c, i, j, and k should be sorted in increasing order before calling this function.
(FPCore (x y z t a b c i j k)
:precision binary64
(let* ((t_1 (* (* j 27.0) k)))
(if (<= t_1 -5e+196)
(fma -27.0 (* k j) (* b c))
(if (<= t_1 5e+263)
(fma (fma i x (* a t)) -4.0 (* c b))
(* (* -27.0 j) k)))))assert(x < y && y < z && z < t && t < a && a < b && b < c && c < i && i < j && j < k);
double code(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) {
double t_1 = (j * 27.0) * k;
double tmp;
if (t_1 <= -5e+196) {
tmp = fma(-27.0, (k * j), (b * c));
} else if (t_1 <= 5e+263) {
tmp = fma(fma(i, x, (a * t)), -4.0, (c * b));
} else {
tmp = (-27.0 * j) * k;
}
return tmp;
}
x, y, z, t, a, b, c, i, j, k = sort([x, y, z, t, a, b, c, i, j, k]) function code(x, y, z, t, a, b, c, i, j, k) t_1 = Float64(Float64(j * 27.0) * k) tmp = 0.0 if (t_1 <= -5e+196) tmp = fma(-27.0, Float64(k * j), Float64(b * c)); elseif (t_1 <= 5e+263) tmp = fma(fma(i, x, Float64(a * t)), -4.0, Float64(c * b)); else tmp = Float64(Float64(-27.0 * j) * k); end return tmp end
NOTE: x, y, z, t, a, b, c, i, j, and k should be sorted in increasing order before calling this function.
code[x_, y_, z_, t_, a_, b_, c_, i_, j_, k_] := Block[{t$95$1 = N[(N[(j * 27.0), $MachinePrecision] * k), $MachinePrecision]}, If[LessEqual[t$95$1, -5e+196], N[(-27.0 * N[(k * j), $MachinePrecision] + N[(b * c), $MachinePrecision]), $MachinePrecision], If[LessEqual[t$95$1, 5e+263], N[(N[(i * x + N[(a * t), $MachinePrecision]), $MachinePrecision] * -4.0 + N[(c * b), $MachinePrecision]), $MachinePrecision], N[(N[(-27.0 * j), $MachinePrecision] * k), $MachinePrecision]]]]
\begin{array}{l}
[x, y, z, t, a, b, c, i, j, k] = \mathsf{sort}([x, y, z, t, a, b, c, i, j, k])\\
\\
\begin{array}{l}
t_1 := \left(j \cdot 27\right) \cdot k\\
\mathbf{if}\;t\_1 \leq -5 \cdot 10^{+196}:\\
\;\;\;\;\mathsf{fma}\left(-27, k \cdot j, b \cdot c\right)\\
\mathbf{elif}\;t\_1 \leq 5 \cdot 10^{+263}:\\
\;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(i, x, a \cdot t\right), -4, c \cdot b\right)\\
\mathbf{else}:\\
\;\;\;\;\left(-27 \cdot j\right) \cdot k\\
\end{array}
\end{array}
if (*.f64 (*.f64 j #s(literal 27 binary64)) k) < -4.9999999999999998e196Initial program 75.7%
Taylor expanded in t around 0
associate--r+N/A
lower--.f64N/A
fp-cancel-sub-sign-invN/A
metadata-evalN/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f6479.2
Applied rewrites79.2%
Taylor expanded in x around 0
Applied rewrites79.5%
if -4.9999999999999998e196 < (*.f64 (*.f64 j #s(literal 27 binary64)) k) < 5.00000000000000022e263Initial program 86.1%
Taylor expanded in j around inf
associate-*r*N/A
lower-*.f64N/A
lower-*.f649.4
Applied rewrites9.4%
Taylor expanded in y around 0
associate-+r+N/A
associate--r+N/A
fp-cancel-sign-sub-invN/A
metadata-evalN/A
associate--r-N/A
fp-cancel-sub-sign-invN/A
metadata-evalN/A
+-commutativeN/A
associate-+r+N/A
+-commutativeN/A
associate-+r+N/A
Applied rewrites83.3%
Taylor expanded in b around 0
Applied rewrites59.0%
Taylor expanded in j around 0
Applied rewrites74.6%
if 5.00000000000000022e263 < (*.f64 (*.f64 j #s(literal 27 binary64)) k) Initial program 69.2%
Taylor expanded in j around inf
associate-*r*N/A
lower-*.f64N/A
lower-*.f64100.0
Applied rewrites100.0%
NOTE: x, y, z, t, a, b, c, i, j, and k should be sorted in increasing order before calling this function.
(FPCore (x y z t a b c i j k)
:precision binary64
(let* ((t_1 (fma i x (* a t))))
(if (or (<= (* b c) -9e+15) (not (<= (* b c) 5e+26)))
(fma t_1 -4.0 (* c b))
(fma t_1 -4.0 (* -27.0 (* k j))))))assert(x < y && y < z && z < t && t < a && a < b && b < c && c < i && i < j && j < k);
double code(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) {
double t_1 = fma(i, x, (a * t));
double tmp;
if (((b * c) <= -9e+15) || !((b * c) <= 5e+26)) {
tmp = fma(t_1, -4.0, (c * b));
} else {
tmp = fma(t_1, -4.0, (-27.0 * (k * j)));
}
return tmp;
}
x, y, z, t, a, b, c, i, j, k = sort([x, y, z, t, a, b, c, i, j, k]) function code(x, y, z, t, a, b, c, i, j, k) t_1 = fma(i, x, Float64(a * t)) tmp = 0.0 if ((Float64(b * c) <= -9e+15) || !(Float64(b * c) <= 5e+26)) tmp = fma(t_1, -4.0, Float64(c * b)); else tmp = fma(t_1, -4.0, Float64(-27.0 * Float64(k * j))); end return tmp end
NOTE: x, y, z, t, a, b, c, i, j, and k should be sorted in increasing order before calling this function.
code[x_, y_, z_, t_, a_, b_, c_, i_, j_, k_] := Block[{t$95$1 = N[(i * x + N[(a * t), $MachinePrecision]), $MachinePrecision]}, If[Or[LessEqual[N[(b * c), $MachinePrecision], -9e+15], N[Not[LessEqual[N[(b * c), $MachinePrecision], 5e+26]], $MachinePrecision]], N[(t$95$1 * -4.0 + N[(c * b), $MachinePrecision]), $MachinePrecision], N[(t$95$1 * -4.0 + N[(-27.0 * N[(k * j), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
[x, y, z, t, a, b, c, i, j, k] = \mathsf{sort}([x, y, z, t, a, b, c, i, j, k])\\
\\
\begin{array}{l}
t_1 := \mathsf{fma}\left(i, x, a \cdot t\right)\\
\mathbf{if}\;b \cdot c \leq -9 \cdot 10^{+15} \lor \neg \left(b \cdot c \leq 5 \cdot 10^{+26}\right):\\
\;\;\;\;\mathsf{fma}\left(t\_1, -4, c \cdot b\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(t\_1, -4, -27 \cdot \left(k \cdot j\right)\right)\\
\end{array}
\end{array}
if (*.f64 b c) < -9e15 or 5.0000000000000001e26 < (*.f64 b c) Initial program 84.0%
Taylor expanded in j around inf
associate-*r*N/A
lower-*.f64N/A
lower-*.f6412.8
Applied rewrites12.8%
Taylor expanded in y around 0
associate-+r+N/A
associate--r+N/A
fp-cancel-sign-sub-invN/A
metadata-evalN/A
associate--r-N/A
fp-cancel-sub-sign-invN/A
metadata-evalN/A
+-commutativeN/A
associate-+r+N/A
+-commutativeN/A
associate-+r+N/A
Applied rewrites87.2%
Taylor expanded in b around 0
Applied rewrites42.5%
Taylor expanded in j around 0
Applied rewrites80.4%
if -9e15 < (*.f64 b c) < 5.0000000000000001e26Initial program 84.0%
Taylor expanded in j around inf
associate-*r*N/A
lower-*.f64N/A
lower-*.f6427.2
Applied rewrites27.2%
Taylor expanded in y around 0
associate-+r+N/A
associate--r+N/A
fp-cancel-sign-sub-invN/A
metadata-evalN/A
associate--r-N/A
fp-cancel-sub-sign-invN/A
metadata-evalN/A
+-commutativeN/A
associate-+r+N/A
+-commutativeN/A
associate-+r+N/A
Applied rewrites79.7%
Taylor expanded in b around 0
Applied rewrites76.4%
Final simplification78.1%
NOTE: x, y, z, t, a, b, c, i, j, and k should be sorted in increasing order before calling this function.
(FPCore (x y z t a b c i j k)
:precision binary64
(if (<= x -5.4e+250)
(* (fma -4.0 i (* (* (* z y) t) 18.0)) x)
(if (<= x 1.7e+88)
(- (fma c b (* -4.0 (fma i x (* a t)))) (* (* j 27.0) k))
(* (fma y (* z (* t 18.0)) (* i -4.0)) x))))assert(x < y && y < z && z < t && t < a && a < b && b < c && c < i && i < j && j < k);
double code(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) {
double tmp;
if (x <= -5.4e+250) {
tmp = fma(-4.0, i, (((z * y) * t) * 18.0)) * x;
} else if (x <= 1.7e+88) {
tmp = fma(c, b, (-4.0 * fma(i, x, (a * t)))) - ((j * 27.0) * k);
} else {
tmp = fma(y, (z * (t * 18.0)), (i * -4.0)) * x;
}
return tmp;
}
x, y, z, t, a, b, c, i, j, k = sort([x, y, z, t, a, b, c, i, j, k]) function code(x, y, z, t, a, b, c, i, j, k) tmp = 0.0 if (x <= -5.4e+250) tmp = Float64(fma(-4.0, i, Float64(Float64(Float64(z * y) * t) * 18.0)) * x); elseif (x <= 1.7e+88) tmp = Float64(fma(c, b, Float64(-4.0 * fma(i, x, Float64(a * t)))) - Float64(Float64(j * 27.0) * k)); else tmp = Float64(fma(y, Float64(z * Float64(t * 18.0)), Float64(i * -4.0)) * x); end return tmp end
NOTE: x, y, z, t, a, b, c, i, j, and k should be sorted in increasing order before calling this function. code[x_, y_, z_, t_, a_, b_, c_, i_, j_, k_] := If[LessEqual[x, -5.4e+250], N[(N[(-4.0 * i + N[(N[(N[(z * y), $MachinePrecision] * t), $MachinePrecision] * 18.0), $MachinePrecision]), $MachinePrecision] * x), $MachinePrecision], If[LessEqual[x, 1.7e+88], N[(N[(c * b + N[(-4.0 * N[(i * x + N[(a * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - N[(N[(j * 27.0), $MachinePrecision] * k), $MachinePrecision]), $MachinePrecision], N[(N[(y * N[(z * N[(t * 18.0), $MachinePrecision]), $MachinePrecision] + N[(i * -4.0), $MachinePrecision]), $MachinePrecision] * x), $MachinePrecision]]]
\begin{array}{l}
[x, y, z, t, a, b, c, i, j, k] = \mathsf{sort}([x, y, z, t, a, b, c, i, j, k])\\
\\
\begin{array}{l}
\mathbf{if}\;x \leq -5.4 \cdot 10^{+250}:\\
\;\;\;\;\mathsf{fma}\left(-4, i, \left(\left(z \cdot y\right) \cdot t\right) \cdot 18\right) \cdot x\\
\mathbf{elif}\;x \leq 1.7 \cdot 10^{+88}:\\
\;\;\;\;\mathsf{fma}\left(c, b, -4 \cdot \mathsf{fma}\left(i, x, a \cdot t\right)\right) - \left(j \cdot 27\right) \cdot k\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(y, z \cdot \left(t \cdot 18\right), i \cdot -4\right) \cdot x\\
\end{array}
\end{array}
if x < -5.4e250Initial program 70.6%
Taylor expanded in x around inf
*-commutativeN/A
lower-*.f64N/A
fp-cancel-sub-sign-invN/A
metadata-evalN/A
+-commutativeN/A
lower-fma.f64N/A
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f6488.1
Applied rewrites88.1%
if -5.4e250 < x < 1.70000000000000002e88Initial program 90.1%
Taylor expanded in y around 0
distribute-lft-outN/A
fp-cancel-sub-sign-invN/A
*-commutativeN/A
lower-fma.f64N/A
metadata-evalN/A
lower-*.f64N/A
+-commutativeN/A
lower-fma.f64N/A
lower-*.f6490.1
Applied rewrites90.1%
if 1.70000000000000002e88 < x Initial program 64.3%
Taylor expanded in x around inf
*-commutativeN/A
lower-*.f64N/A
fp-cancel-sub-sign-invN/A
metadata-evalN/A
+-commutativeN/A
lower-fma.f64N/A
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f6483.7
Applied rewrites83.7%
Applied rewrites83.7%
NOTE: x, y, z, t, a, b, c, i, j, and k should be sorted in increasing order before calling this function.
(FPCore (x y z t a b c i j k)
:precision binary64
(if (<= x -4.4e+250)
(* (fma -4.0 i (* (* (* z y) t) 18.0)) x)
(if (<= x -1.7e-106)
(fma (fma i x (* a t)) -4.0 (* c b))
(if (<= x 5e+87)
(fma (* t -4.0) a (fma (* k j) -27.0 (* c b)))
(* (fma y (* z (* t 18.0)) (* i -4.0)) x)))))assert(x < y && y < z && z < t && t < a && a < b && b < c && c < i && i < j && j < k);
double code(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) {
double tmp;
if (x <= -4.4e+250) {
tmp = fma(-4.0, i, (((z * y) * t) * 18.0)) * x;
} else if (x <= -1.7e-106) {
tmp = fma(fma(i, x, (a * t)), -4.0, (c * b));
} else if (x <= 5e+87) {
tmp = fma((t * -4.0), a, fma((k * j), -27.0, (c * b)));
} else {
tmp = fma(y, (z * (t * 18.0)), (i * -4.0)) * x;
}
return tmp;
}
x, y, z, t, a, b, c, i, j, k = sort([x, y, z, t, a, b, c, i, j, k]) function code(x, y, z, t, a, b, c, i, j, k) tmp = 0.0 if (x <= -4.4e+250) tmp = Float64(fma(-4.0, i, Float64(Float64(Float64(z * y) * t) * 18.0)) * x); elseif (x <= -1.7e-106) tmp = fma(fma(i, x, Float64(a * t)), -4.0, Float64(c * b)); elseif (x <= 5e+87) tmp = fma(Float64(t * -4.0), a, fma(Float64(k * j), -27.0, Float64(c * b))); else tmp = Float64(fma(y, Float64(z * Float64(t * 18.0)), Float64(i * -4.0)) * x); end return tmp end
NOTE: x, y, z, t, a, b, c, i, j, and k should be sorted in increasing order before calling this function. code[x_, y_, z_, t_, a_, b_, c_, i_, j_, k_] := If[LessEqual[x, -4.4e+250], N[(N[(-4.0 * i + N[(N[(N[(z * y), $MachinePrecision] * t), $MachinePrecision] * 18.0), $MachinePrecision]), $MachinePrecision] * x), $MachinePrecision], If[LessEqual[x, -1.7e-106], N[(N[(i * x + N[(a * t), $MachinePrecision]), $MachinePrecision] * -4.0 + N[(c * b), $MachinePrecision]), $MachinePrecision], If[LessEqual[x, 5e+87], N[(N[(t * -4.0), $MachinePrecision] * a + N[(N[(k * j), $MachinePrecision] * -27.0 + N[(c * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(y * N[(z * N[(t * 18.0), $MachinePrecision]), $MachinePrecision] + N[(i * -4.0), $MachinePrecision]), $MachinePrecision] * x), $MachinePrecision]]]]
\begin{array}{l}
[x, y, z, t, a, b, c, i, j, k] = \mathsf{sort}([x, y, z, t, a, b, c, i, j, k])\\
\\
\begin{array}{l}
\mathbf{if}\;x \leq -4.4 \cdot 10^{+250}:\\
\;\;\;\;\mathsf{fma}\left(-4, i, \left(\left(z \cdot y\right) \cdot t\right) \cdot 18\right) \cdot x\\
\mathbf{elif}\;x \leq -1.7 \cdot 10^{-106}:\\
\;\;\;\;\mathsf{fma}\left(\mathsf{fma}\left(i, x, a \cdot t\right), -4, c \cdot b\right)\\
\mathbf{elif}\;x \leq 5 \cdot 10^{+87}:\\
\;\;\;\;\mathsf{fma}\left(t \cdot -4, a, \mathsf{fma}\left(k \cdot j, -27, c \cdot b\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(y, z \cdot \left(t \cdot 18\right), i \cdot -4\right) \cdot x\\
\end{array}
\end{array}
if x < -4.40000000000000029e250Initial program 70.6%
Taylor expanded in x around inf
*-commutativeN/A
lower-*.f64N/A
fp-cancel-sub-sign-invN/A
metadata-evalN/A
+-commutativeN/A
lower-fma.f64N/A
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f6488.1
Applied rewrites88.1%
if -4.40000000000000029e250 < x < -1.69999999999999991e-106Initial program 88.3%
Taylor expanded in j around inf
associate-*r*N/A
lower-*.f64N/A
lower-*.f6413.9
Applied rewrites13.9%
Taylor expanded in y around 0
associate-+r+N/A
associate--r+N/A
fp-cancel-sign-sub-invN/A
metadata-evalN/A
associate--r-N/A
fp-cancel-sub-sign-invN/A
metadata-evalN/A
+-commutativeN/A
associate-+r+N/A
+-commutativeN/A
associate-+r+N/A
Applied rewrites87.8%
Taylor expanded in b around 0
Applied rewrites67.0%
Taylor expanded in j around 0
Applied rewrites81.2%
if -1.69999999999999991e-106 < x < 4.9999999999999998e87Initial program 90.8%
Taylor expanded in j around inf
associate-*r*N/A
lower-*.f64N/A
lower-*.f6428.9
Applied rewrites28.9%
Taylor expanded in x around 0
associate--r+N/A
fp-cancel-sub-sign-invN/A
metadata-evalN/A
+-commutativeN/A
associate--l+N/A
*-commutativeN/A
lower-fma.f64N/A
*-commutativeN/A
lower-*.f64N/A
fp-cancel-sub-sign-invN/A
metadata-evalN/A
+-commutativeN/A
lower-fma.f64N/A
*-commutativeN/A
lower-*.f64N/A
lower-*.f6481.8
Applied rewrites81.8%
Applied rewrites83.3%
if 4.9999999999999998e87 < x Initial program 64.3%
Taylor expanded in x around inf
*-commutativeN/A
lower-*.f64N/A
fp-cancel-sub-sign-invN/A
metadata-evalN/A
+-commutativeN/A
lower-fma.f64N/A
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f6483.7
Applied rewrites83.7%
Applied rewrites83.7%
NOTE: x, y, z, t, a, b, c, i, j, and k should be sorted in increasing order before calling this function.
(FPCore (x y z t a b c i j k)
:precision binary64
(if (<= (* b c) -5e+60)
(* c b)
(if (<= (* b c) -5e-302)
(* (* -27.0 j) k)
(if (<= (* b c) 2e+116) (* (* -4.0 a) t) (* c b)))))assert(x < y && y < z && z < t && t < a && a < b && b < c && c < i && i < j && j < k);
double code(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) {
double tmp;
if ((b * c) <= -5e+60) {
tmp = c * b;
} else if ((b * c) <= -5e-302) {
tmp = (-27.0 * j) * k;
} else if ((b * c) <= 2e+116) {
tmp = (-4.0 * a) * t;
} else {
tmp = c * b;
}
return tmp;
}
NOTE: x, y, z, t, a, b, c, i, j, and k should be sorted in increasing order before calling this function.
module fmin_fmax_functions
implicit none
private
public fmax
public fmin
interface fmax
module procedure fmax88
module procedure fmax44
module procedure fmax84
module procedure fmax48
end interface
interface fmin
module procedure fmin88
module procedure fmin44
module procedure fmin84
module procedure fmin48
end interface
contains
real(8) function fmax88(x, y) result (res)
real(8), intent (in) :: x
real(8), intent (in) :: y
res = merge(y, merge(x, max(x, y), y /= y), x /= x)
end function
real(4) function fmax44(x, y) result (res)
real(4), intent (in) :: x
real(4), intent (in) :: y
res = merge(y, merge(x, max(x, y), y /= y), x /= x)
end function
real(8) function fmax84(x, y) result(res)
real(8), intent (in) :: x
real(4), intent (in) :: y
res = merge(dble(y), merge(x, max(x, dble(y)), y /= y), x /= x)
end function
real(8) function fmax48(x, y) result(res)
real(4), intent (in) :: x
real(8), intent (in) :: y
res = merge(y, merge(dble(x), max(dble(x), y), y /= y), x /= x)
end function
real(8) function fmin88(x, y) result (res)
real(8), intent (in) :: x
real(8), intent (in) :: y
res = merge(y, merge(x, min(x, y), y /= y), x /= x)
end function
real(4) function fmin44(x, y) result (res)
real(4), intent (in) :: x
real(4), intent (in) :: y
res = merge(y, merge(x, min(x, y), y /= y), x /= x)
end function
real(8) function fmin84(x, y) result(res)
real(8), intent (in) :: x
real(4), intent (in) :: y
res = merge(dble(y), merge(x, min(x, dble(y)), y /= y), x /= x)
end function
real(8) function fmin48(x, y) result(res)
real(4), intent (in) :: x
real(8), intent (in) :: y
res = merge(y, merge(dble(x), min(dble(x), y), y /= y), x /= x)
end function
end module
real(8) function code(x, y, z, t, a, b, c, i, j, k)
use fmin_fmax_functions
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), intent (in) :: b
real(8), intent (in) :: c
real(8), intent (in) :: i
real(8), intent (in) :: j
real(8), intent (in) :: k
real(8) :: tmp
if ((b * c) <= (-5d+60)) then
tmp = c * b
else if ((b * c) <= (-5d-302)) then
tmp = ((-27.0d0) * j) * k
else if ((b * c) <= 2d+116) then
tmp = ((-4.0d0) * a) * t
else
tmp = c * b
end if
code = tmp
end function
assert x < y && y < z && z < t && t < a && a < b && b < c && c < i && i < j && j < k;
public static double code(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) {
double tmp;
if ((b * c) <= -5e+60) {
tmp = c * b;
} else if ((b * c) <= -5e-302) {
tmp = (-27.0 * j) * k;
} else if ((b * c) <= 2e+116) {
tmp = (-4.0 * a) * t;
} else {
tmp = c * b;
}
return tmp;
}
[x, y, z, t, a, b, c, i, j, k] = sort([x, y, z, t, a, b, c, i, j, k]) def code(x, y, z, t, a, b, c, i, j, k): tmp = 0 if (b * c) <= -5e+60: tmp = c * b elif (b * c) <= -5e-302: tmp = (-27.0 * j) * k elif (b * c) <= 2e+116: tmp = (-4.0 * a) * t else: tmp = c * b return tmp
x, y, z, t, a, b, c, i, j, k = sort([x, y, z, t, a, b, c, i, j, k]) function code(x, y, z, t, a, b, c, i, j, k) tmp = 0.0 if (Float64(b * c) <= -5e+60) tmp = Float64(c * b); elseif (Float64(b * c) <= -5e-302) tmp = Float64(Float64(-27.0 * j) * k); elseif (Float64(b * c) <= 2e+116) tmp = Float64(Float64(-4.0 * a) * t); else tmp = Float64(c * b); end return tmp end
x, y, z, t, a, b, c, i, j, k = num2cell(sort([x, y, z, t, a, b, c, i, j, k])){:}
function tmp_2 = code(x, y, z, t, a, b, c, i, j, k)
tmp = 0.0;
if ((b * c) <= -5e+60)
tmp = c * b;
elseif ((b * c) <= -5e-302)
tmp = (-27.0 * j) * k;
elseif ((b * c) <= 2e+116)
tmp = (-4.0 * a) * t;
else
tmp = c * b;
end
tmp_2 = tmp;
end
NOTE: x, y, z, t, a, b, c, i, j, and k should be sorted in increasing order before calling this function. code[x_, y_, z_, t_, a_, b_, c_, i_, j_, k_] := If[LessEqual[N[(b * c), $MachinePrecision], -5e+60], N[(c * b), $MachinePrecision], If[LessEqual[N[(b * c), $MachinePrecision], -5e-302], N[(N[(-27.0 * j), $MachinePrecision] * k), $MachinePrecision], If[LessEqual[N[(b * c), $MachinePrecision], 2e+116], N[(N[(-4.0 * a), $MachinePrecision] * t), $MachinePrecision], N[(c * b), $MachinePrecision]]]]
\begin{array}{l}
[x, y, z, t, a, b, c, i, j, k] = \mathsf{sort}([x, y, z, t, a, b, c, i, j, k])\\
\\
\begin{array}{l}
\mathbf{if}\;b \cdot c \leq -5 \cdot 10^{+60}:\\
\;\;\;\;c \cdot b\\
\mathbf{elif}\;b \cdot c \leq -5 \cdot 10^{-302}:\\
\;\;\;\;\left(-27 \cdot j\right) \cdot k\\
\mathbf{elif}\;b \cdot c \leq 2 \cdot 10^{+116}:\\
\;\;\;\;\left(-4 \cdot a\right) \cdot t\\
\mathbf{else}:\\
\;\;\;\;c \cdot b\\
\end{array}
\end{array}
if (*.f64 b c) < -4.99999999999999975e60 or 2.00000000000000003e116 < (*.f64 b c) Initial program 84.0%
Taylor expanded in y around 0
distribute-lft-outN/A
fp-cancel-sub-sign-invN/A
*-commutativeN/A
lower-fma.f64N/A
metadata-evalN/A
lower-*.f64N/A
+-commutativeN/A
lower-fma.f64N/A
lower-*.f6489.1
Applied rewrites89.1%
Taylor expanded in a around 0
associate--r+N/A
fp-cancel-sub-sign-invN/A
metadata-evalN/A
+-commutativeN/A
associate-+r+N/A
associate--l+N/A
associate-*r*N/A
associate-*r*N/A
*-commutativeN/A
associate-*r*N/A
associate-*r*N/A
distribute-rgt-inN/A
*-commutativeN/A
Applied rewrites77.8%
Taylor expanded in b around inf
*-commutativeN/A
lower-*.f6460.3
Applied rewrites60.3%
if -4.99999999999999975e60 < (*.f64 b c) < -5.00000000000000033e-302Initial program 84.4%
Taylor expanded in j around inf
associate-*r*N/A
lower-*.f64N/A
lower-*.f6437.4
Applied rewrites37.4%
if -5.00000000000000033e-302 < (*.f64 b c) < 2.00000000000000003e116Initial program 83.9%
lift--.f64N/A
lift-*.f64N/A
fp-cancel-sub-sign-invN/A
+-commutativeN/A
lower-fma.f64N/A
lift-*.f64N/A
*-commutativeN/A
distribute-lft-neg-inN/A
lower-*.f64N/A
metadata-eval86.6
lift-+.f64N/A
lift--.f64N/A
lift-*.f64N/A
lift-*.f64N/A
distribute-rgt-out--N/A
*-commutativeN/A
Applied rewrites87.5%
Taylor expanded in a around inf
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f6434.5
Applied rewrites34.5%
Applied rewrites34.5%
NOTE: x, y, z, t, a, b, c, i, j, and k should be sorted in increasing order before calling this function. (FPCore (x y z t a b c i j k) :precision binary64 (if (or (<= (* b c) -5e-9) (not (<= (* b c) 2000000000000.0))) (fma -27.0 (* k j) (* b c)) (* (fma i x (* a t)) -4.0)))
assert(x < y && y < z && z < t && t < a && a < b && b < c && c < i && i < j && j < k);
double code(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) {
double tmp;
if (((b * c) <= -5e-9) || !((b * c) <= 2000000000000.0)) {
tmp = fma(-27.0, (k * j), (b * c));
} else {
tmp = fma(i, x, (a * t)) * -4.0;
}
return tmp;
}
x, y, z, t, a, b, c, i, j, k = sort([x, y, z, t, a, b, c, i, j, k]) function code(x, y, z, t, a, b, c, i, j, k) tmp = 0.0 if ((Float64(b * c) <= -5e-9) || !(Float64(b * c) <= 2000000000000.0)) tmp = fma(-27.0, Float64(k * j), Float64(b * c)); else tmp = Float64(fma(i, x, Float64(a * t)) * -4.0); end return tmp end
NOTE: x, y, z, t, a, b, c, i, j, and k should be sorted in increasing order before calling this function. code[x_, y_, z_, t_, a_, b_, c_, i_, j_, k_] := If[Or[LessEqual[N[(b * c), $MachinePrecision], -5e-9], N[Not[LessEqual[N[(b * c), $MachinePrecision], 2000000000000.0]], $MachinePrecision]], N[(-27.0 * N[(k * j), $MachinePrecision] + N[(b * c), $MachinePrecision]), $MachinePrecision], N[(N[(i * x + N[(a * t), $MachinePrecision]), $MachinePrecision] * -4.0), $MachinePrecision]]
\begin{array}{l}
[x, y, z, t, a, b, c, i, j, k] = \mathsf{sort}([x, y, z, t, a, b, c, i, j, k])\\
\\
\begin{array}{l}
\mathbf{if}\;b \cdot c \leq -5 \cdot 10^{-9} \lor \neg \left(b \cdot c \leq 2000000000000\right):\\
\;\;\;\;\mathsf{fma}\left(-27, k \cdot j, b \cdot c\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(i, x, a \cdot t\right) \cdot -4\\
\end{array}
\end{array}
if (*.f64 b c) < -5.0000000000000001e-9 or 2e12 < (*.f64 b c) Initial program 83.8%
Taylor expanded in t around 0
associate--r+N/A
lower--.f64N/A
fp-cancel-sub-sign-invN/A
metadata-evalN/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f6476.4
Applied rewrites76.4%
Taylor expanded in x around 0
Applied rewrites64.1%
if -5.0000000000000001e-9 < (*.f64 b c) < 2e12Initial program 84.3%
Taylor expanded in j around inf
associate-*r*N/A
lower-*.f64N/A
lower-*.f6423.5
Applied rewrites23.5%
Taylor expanded in y around 0
associate-+r+N/A
associate--r+N/A
fp-cancel-sign-sub-invN/A
metadata-evalN/A
associate--r-N/A
fp-cancel-sub-sign-invN/A
metadata-evalN/A
+-commutativeN/A
associate-+r+N/A
+-commutativeN/A
associate-+r+N/A
Applied rewrites79.6%
Taylor expanded in b around 0
Applied rewrites76.7%
Taylor expanded in j around 0
Applied rewrites58.3%
Final simplification61.0%
NOTE: x, y, z, t, a, b, c, i, j, and k should be sorted in increasing order before calling this function. (FPCore (x y z t a b c i j k) :precision binary64 (if (or (<= (* b c) -1e+153) (not (<= (* b c) 4e+139))) (* c b) (* (fma i x (* a t)) -4.0)))
assert(x < y && y < z && z < t && t < a && a < b && b < c && c < i && i < j && j < k);
double code(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) {
double tmp;
if (((b * c) <= -1e+153) || !((b * c) <= 4e+139)) {
tmp = c * b;
} else {
tmp = fma(i, x, (a * t)) * -4.0;
}
return tmp;
}
x, y, z, t, a, b, c, i, j, k = sort([x, y, z, t, a, b, c, i, j, k]) function code(x, y, z, t, a, b, c, i, j, k) tmp = 0.0 if ((Float64(b * c) <= -1e+153) || !(Float64(b * c) <= 4e+139)) tmp = Float64(c * b); else tmp = Float64(fma(i, x, Float64(a * t)) * -4.0); end return tmp end
NOTE: x, y, z, t, a, b, c, i, j, and k should be sorted in increasing order before calling this function. code[x_, y_, z_, t_, a_, b_, c_, i_, j_, k_] := If[Or[LessEqual[N[(b * c), $MachinePrecision], -1e+153], N[Not[LessEqual[N[(b * c), $MachinePrecision], 4e+139]], $MachinePrecision]], N[(c * b), $MachinePrecision], N[(N[(i * x + N[(a * t), $MachinePrecision]), $MachinePrecision] * -4.0), $MachinePrecision]]
\begin{array}{l}
[x, y, z, t, a, b, c, i, j, k] = \mathsf{sort}([x, y, z, t, a, b, c, i, j, k])\\
\\
\begin{array}{l}
\mathbf{if}\;b \cdot c \leq -1 \cdot 10^{+153} \lor \neg \left(b \cdot c \leq 4 \cdot 10^{+139}\right):\\
\;\;\;\;c \cdot b\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(i, x, a \cdot t\right) \cdot -4\\
\end{array}
\end{array}
if (*.f64 b c) < -1e153 or 4.00000000000000013e139 < (*.f64 b c) Initial program 81.8%
Taylor expanded in y around 0
distribute-lft-outN/A
fp-cancel-sub-sign-invN/A
*-commutativeN/A
lower-fma.f64N/A
metadata-evalN/A
lower-*.f64N/A
+-commutativeN/A
lower-fma.f64N/A
lower-*.f6488.0
Applied rewrites88.0%
Taylor expanded in a around 0
associate--r+N/A
fp-cancel-sub-sign-invN/A
metadata-evalN/A
+-commutativeN/A
associate-+r+N/A
associate--l+N/A
associate-*r*N/A
associate-*r*N/A
*-commutativeN/A
associate-*r*N/A
associate-*r*N/A
distribute-rgt-inN/A
*-commutativeN/A
Applied rewrites77.9%
Taylor expanded in b around inf
*-commutativeN/A
lower-*.f6470.9
Applied rewrites70.9%
if -1e153 < (*.f64 b c) < 4.00000000000000013e139Initial program 84.8%
Taylor expanded in j around inf
associate-*r*N/A
lower-*.f64N/A
lower-*.f6425.6
Applied rewrites25.6%
Taylor expanded in y around 0
associate-+r+N/A
associate--r+N/A
fp-cancel-sign-sub-invN/A
metadata-evalN/A
associate--r-N/A
fp-cancel-sub-sign-invN/A
metadata-evalN/A
+-commutativeN/A
associate-+r+N/A
+-commutativeN/A
associate-+r+N/A
Applied rewrites82.1%
Taylor expanded in b around 0
Applied rewrites74.5%
Taylor expanded in j around 0
Applied rewrites54.4%
Final simplification58.7%
NOTE: x, y, z, t, a, b, c, i, j, and k should be sorted in increasing order before calling this function. (FPCore (x y z t a b c i j k) :precision binary64 (if (or (<= (* b c) -5e+131) (not (<= (* b c) 2e+116))) (* c b) (* (* -4.0 a) t)))
assert(x < y && y < z && z < t && t < a && a < b && b < c && c < i && i < j && j < k);
double code(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) {
double tmp;
if (((b * c) <= -5e+131) || !((b * c) <= 2e+116)) {
tmp = c * b;
} else {
tmp = (-4.0 * a) * t;
}
return tmp;
}
NOTE: x, y, z, t, a, b, c, i, j, and k should be sorted in increasing order before calling this function.
module fmin_fmax_functions
implicit none
private
public fmax
public fmin
interface fmax
module procedure fmax88
module procedure fmax44
module procedure fmax84
module procedure fmax48
end interface
interface fmin
module procedure fmin88
module procedure fmin44
module procedure fmin84
module procedure fmin48
end interface
contains
real(8) function fmax88(x, y) result (res)
real(8), intent (in) :: x
real(8), intent (in) :: y
res = merge(y, merge(x, max(x, y), y /= y), x /= x)
end function
real(4) function fmax44(x, y) result (res)
real(4), intent (in) :: x
real(4), intent (in) :: y
res = merge(y, merge(x, max(x, y), y /= y), x /= x)
end function
real(8) function fmax84(x, y) result(res)
real(8), intent (in) :: x
real(4), intent (in) :: y
res = merge(dble(y), merge(x, max(x, dble(y)), y /= y), x /= x)
end function
real(8) function fmax48(x, y) result(res)
real(4), intent (in) :: x
real(8), intent (in) :: y
res = merge(y, merge(dble(x), max(dble(x), y), y /= y), x /= x)
end function
real(8) function fmin88(x, y) result (res)
real(8), intent (in) :: x
real(8), intent (in) :: y
res = merge(y, merge(x, min(x, y), y /= y), x /= x)
end function
real(4) function fmin44(x, y) result (res)
real(4), intent (in) :: x
real(4), intent (in) :: y
res = merge(y, merge(x, min(x, y), y /= y), x /= x)
end function
real(8) function fmin84(x, y) result(res)
real(8), intent (in) :: x
real(4), intent (in) :: y
res = merge(dble(y), merge(x, min(x, dble(y)), y /= y), x /= x)
end function
real(8) function fmin48(x, y) result(res)
real(4), intent (in) :: x
real(8), intent (in) :: y
res = merge(y, merge(dble(x), min(dble(x), y), y /= y), x /= x)
end function
end module
real(8) function code(x, y, z, t, a, b, c, i, j, k)
use fmin_fmax_functions
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), intent (in) :: b
real(8), intent (in) :: c
real(8), intent (in) :: i
real(8), intent (in) :: j
real(8), intent (in) :: k
real(8) :: tmp
if (((b * c) <= (-5d+131)) .or. (.not. ((b * c) <= 2d+116))) then
tmp = c * b
else
tmp = ((-4.0d0) * a) * t
end if
code = tmp
end function
assert x < y && y < z && z < t && t < a && a < b && b < c && c < i && i < j && j < k;
public static double code(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) {
double tmp;
if (((b * c) <= -5e+131) || !((b * c) <= 2e+116)) {
tmp = c * b;
} else {
tmp = (-4.0 * a) * t;
}
return tmp;
}
[x, y, z, t, a, b, c, i, j, k] = sort([x, y, z, t, a, b, c, i, j, k]) def code(x, y, z, t, a, b, c, i, j, k): tmp = 0 if ((b * c) <= -5e+131) or not ((b * c) <= 2e+116): tmp = c * b else: tmp = (-4.0 * a) * t return tmp
x, y, z, t, a, b, c, i, j, k = sort([x, y, z, t, a, b, c, i, j, k]) function code(x, y, z, t, a, b, c, i, j, k) tmp = 0.0 if ((Float64(b * c) <= -5e+131) || !(Float64(b * c) <= 2e+116)) tmp = Float64(c * b); else tmp = Float64(Float64(-4.0 * a) * t); end return tmp end
x, y, z, t, a, b, c, i, j, k = num2cell(sort([x, y, z, t, a, b, c, i, j, k])){:}
function tmp_2 = code(x, y, z, t, a, b, c, i, j, k)
tmp = 0.0;
if (((b * c) <= -5e+131) || ~(((b * c) <= 2e+116)))
tmp = c * b;
else
tmp = (-4.0 * a) * t;
end
tmp_2 = tmp;
end
NOTE: x, y, z, t, a, b, c, i, j, and k should be sorted in increasing order before calling this function. code[x_, y_, z_, t_, a_, b_, c_, i_, j_, k_] := If[Or[LessEqual[N[(b * c), $MachinePrecision], -5e+131], N[Not[LessEqual[N[(b * c), $MachinePrecision], 2e+116]], $MachinePrecision]], N[(c * b), $MachinePrecision], N[(N[(-4.0 * a), $MachinePrecision] * t), $MachinePrecision]]
\begin{array}{l}
[x, y, z, t, a, b, c, i, j, k] = \mathsf{sort}([x, y, z, t, a, b, c, i, j, k])\\
\\
\begin{array}{l}
\mathbf{if}\;b \cdot c \leq -5 \cdot 10^{+131} \lor \neg \left(b \cdot c \leq 2 \cdot 10^{+116}\right):\\
\;\;\;\;c \cdot b\\
\mathbf{else}:\\
\;\;\;\;\left(-4 \cdot a\right) \cdot t\\
\end{array}
\end{array}
if (*.f64 b c) < -4.99999999999999995e131 or 2.00000000000000003e116 < (*.f64 b c) Initial program 81.9%
Taylor expanded in y around 0
distribute-lft-outN/A
fp-cancel-sub-sign-invN/A
*-commutativeN/A
lower-fma.f64N/A
metadata-evalN/A
lower-*.f64N/A
+-commutativeN/A
lower-fma.f64N/A
lower-*.f6487.6
Applied rewrites87.6%
Taylor expanded in a around 0
associate--r+N/A
fp-cancel-sub-sign-invN/A
metadata-evalN/A
+-commutativeN/A
associate-+r+N/A
associate--l+N/A
associate-*r*N/A
associate-*r*N/A
*-commutativeN/A
associate-*r*N/A
associate-*r*N/A
distribute-rgt-inN/A
*-commutativeN/A
Applied rewrites79.7%
Taylor expanded in b around inf
*-commutativeN/A
lower-*.f6466.5
Applied rewrites66.5%
if -4.99999999999999995e131 < (*.f64 b c) < 2.00000000000000003e116Initial program 84.9%
lift--.f64N/A
lift-*.f64N/A
fp-cancel-sub-sign-invN/A
+-commutativeN/A
lower-fma.f64N/A
lift-*.f64N/A
*-commutativeN/A
distribute-lft-neg-inN/A
lower-*.f64N/A
metadata-eval87.6
lift-+.f64N/A
lift--.f64N/A
lift-*.f64N/A
lift-*.f64N/A
distribute-rgt-out--N/A
*-commutativeN/A
Applied rewrites89.2%
Taylor expanded in a around inf
*-commutativeN/A
lower-*.f64N/A
*-commutativeN/A
lower-*.f6430.7
Applied rewrites30.7%
Applied rewrites30.7%
Final simplification40.8%
NOTE: x, y, z, t, a, b, c, i, j, and k should be sorted in increasing order before calling this function. (FPCore (x y z t a b c i j k) :precision binary64 (* c b))
assert(x < y && y < z && z < t && t < a && a < b && b < c && c < i && i < j && j < k);
double code(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) {
return c * b;
}
NOTE: x, y, z, t, a, b, c, i, j, and k should be sorted in increasing order before calling this function.
module fmin_fmax_functions
implicit none
private
public fmax
public fmin
interface fmax
module procedure fmax88
module procedure fmax44
module procedure fmax84
module procedure fmax48
end interface
interface fmin
module procedure fmin88
module procedure fmin44
module procedure fmin84
module procedure fmin48
end interface
contains
real(8) function fmax88(x, y) result (res)
real(8), intent (in) :: x
real(8), intent (in) :: y
res = merge(y, merge(x, max(x, y), y /= y), x /= x)
end function
real(4) function fmax44(x, y) result (res)
real(4), intent (in) :: x
real(4), intent (in) :: y
res = merge(y, merge(x, max(x, y), y /= y), x /= x)
end function
real(8) function fmax84(x, y) result(res)
real(8), intent (in) :: x
real(4), intent (in) :: y
res = merge(dble(y), merge(x, max(x, dble(y)), y /= y), x /= x)
end function
real(8) function fmax48(x, y) result(res)
real(4), intent (in) :: x
real(8), intent (in) :: y
res = merge(y, merge(dble(x), max(dble(x), y), y /= y), x /= x)
end function
real(8) function fmin88(x, y) result (res)
real(8), intent (in) :: x
real(8), intent (in) :: y
res = merge(y, merge(x, min(x, y), y /= y), x /= x)
end function
real(4) function fmin44(x, y) result (res)
real(4), intent (in) :: x
real(4), intent (in) :: y
res = merge(y, merge(x, min(x, y), y /= y), x /= x)
end function
real(8) function fmin84(x, y) result(res)
real(8), intent (in) :: x
real(4), intent (in) :: y
res = merge(dble(y), merge(x, min(x, dble(y)), y /= y), x /= x)
end function
real(8) function fmin48(x, y) result(res)
real(4), intent (in) :: x
real(8), intent (in) :: y
res = merge(y, merge(dble(x), min(dble(x), y), y /= y), x /= x)
end function
end module
real(8) function code(x, y, z, t, a, b, c, i, j, k)
use fmin_fmax_functions
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), intent (in) :: b
real(8), intent (in) :: c
real(8), intent (in) :: i
real(8), intent (in) :: j
real(8), intent (in) :: k
code = c * b
end function
assert x < y && y < z && z < t && t < a && a < b && b < c && c < i && i < j && j < k;
public static double code(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) {
return c * b;
}
[x, y, z, t, a, b, c, i, j, k] = sort([x, y, z, t, a, b, c, i, j, k]) def code(x, y, z, t, a, b, c, i, j, k): return c * b
x, y, z, t, a, b, c, i, j, k = sort([x, y, z, t, a, b, c, i, j, k]) function code(x, y, z, t, a, b, c, i, j, k) return Float64(c * b) end
x, y, z, t, a, b, c, i, j, k = num2cell(sort([x, y, z, t, a, b, c, i, j, k])){:}
function tmp = code(x, y, z, t, a, b, c, i, j, k)
tmp = c * b;
end
NOTE: x, y, z, t, a, b, c, i, j, and k should be sorted in increasing order before calling this function. code[x_, y_, z_, t_, a_, b_, c_, i_, j_, k_] := N[(c * b), $MachinePrecision]
\begin{array}{l}
[x, y, z, t, a, b, c, i, j, k] = \mathsf{sort}([x, y, z, t, a, b, c, i, j, k])\\
\\
c \cdot b
\end{array}
Initial program 84.0%
Taylor expanded in y around 0
distribute-lft-outN/A
fp-cancel-sub-sign-invN/A
*-commutativeN/A
lower-fma.f64N/A
metadata-evalN/A
lower-*.f64N/A
+-commutativeN/A
lower-fma.f64N/A
lower-*.f6483.2
Applied rewrites83.2%
Taylor expanded in a around 0
associate--r+N/A
fp-cancel-sub-sign-invN/A
metadata-evalN/A
+-commutativeN/A
associate-+r+N/A
associate--l+N/A
associate-*r*N/A
associate-*r*N/A
*-commutativeN/A
associate-*r*N/A
associate-*r*N/A
distribute-rgt-inN/A
*-commutativeN/A
Applied rewrites72.3%
Taylor expanded in b around inf
*-commutativeN/A
lower-*.f6424.9
Applied rewrites24.9%
(FPCore (x y z t a b c i j k)
:precision binary64
(let* ((t_1 (* (+ (* a t) (* i x)) 4.0))
(t_2
(-
(- (* (* 18.0 t) (* (* x y) z)) t_1)
(- (* (* k j) 27.0) (* c b)))))
(if (< t -1.6210815397541398e-69)
t_2
(if (< t 165.68027943805222)
(+ (- (* (* 18.0 y) (* x (* z t))) t_1) (- (* c b) (* 27.0 (* k j))))
t_2))))
double code(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) {
double t_1 = ((a * t) + (i * x)) * 4.0;
double t_2 = (((18.0 * t) * ((x * y) * z)) - t_1) - (((k * j) * 27.0) - (c * b));
double tmp;
if (t < -1.6210815397541398e-69) {
tmp = t_2;
} else if (t < 165.68027943805222) {
tmp = (((18.0 * y) * (x * (z * t))) - t_1) + ((c * b) - (27.0 * (k * j)));
} else {
tmp = t_2;
}
return tmp;
}
module fmin_fmax_functions
implicit none
private
public fmax
public fmin
interface fmax
module procedure fmax88
module procedure fmax44
module procedure fmax84
module procedure fmax48
end interface
interface fmin
module procedure fmin88
module procedure fmin44
module procedure fmin84
module procedure fmin48
end interface
contains
real(8) function fmax88(x, y) result (res)
real(8), intent (in) :: x
real(8), intent (in) :: y
res = merge(y, merge(x, max(x, y), y /= y), x /= x)
end function
real(4) function fmax44(x, y) result (res)
real(4), intent (in) :: x
real(4), intent (in) :: y
res = merge(y, merge(x, max(x, y), y /= y), x /= x)
end function
real(8) function fmax84(x, y) result(res)
real(8), intent (in) :: x
real(4), intent (in) :: y
res = merge(dble(y), merge(x, max(x, dble(y)), y /= y), x /= x)
end function
real(8) function fmax48(x, y) result(res)
real(4), intent (in) :: x
real(8), intent (in) :: y
res = merge(y, merge(dble(x), max(dble(x), y), y /= y), x /= x)
end function
real(8) function fmin88(x, y) result (res)
real(8), intent (in) :: x
real(8), intent (in) :: y
res = merge(y, merge(x, min(x, y), y /= y), x /= x)
end function
real(4) function fmin44(x, y) result (res)
real(4), intent (in) :: x
real(4), intent (in) :: y
res = merge(y, merge(x, min(x, y), y /= y), x /= x)
end function
real(8) function fmin84(x, y) result(res)
real(8), intent (in) :: x
real(4), intent (in) :: y
res = merge(dble(y), merge(x, min(x, dble(y)), y /= y), x /= x)
end function
real(8) function fmin48(x, y) result(res)
real(4), intent (in) :: x
real(8), intent (in) :: y
res = merge(y, merge(dble(x), min(dble(x), y), y /= y), x /= x)
end function
end module
real(8) function code(x, y, z, t, a, b, c, i, j, k)
use fmin_fmax_functions
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), intent (in) :: b
real(8), intent (in) :: c
real(8), intent (in) :: i
real(8), intent (in) :: j
real(8), intent (in) :: k
real(8) :: t_1
real(8) :: t_2
real(8) :: tmp
t_1 = ((a * t) + (i * x)) * 4.0d0
t_2 = (((18.0d0 * t) * ((x * y) * z)) - t_1) - (((k * j) * 27.0d0) - (c * b))
if (t < (-1.6210815397541398d-69)) then
tmp = t_2
else if (t < 165.68027943805222d0) then
tmp = (((18.0d0 * y) * (x * (z * t))) - t_1) + ((c * b) - (27.0d0 * (k * j)))
else
tmp = t_2
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b, double c, double i, double j, double k) {
double t_1 = ((a * t) + (i * x)) * 4.0;
double t_2 = (((18.0 * t) * ((x * y) * z)) - t_1) - (((k * j) * 27.0) - (c * b));
double tmp;
if (t < -1.6210815397541398e-69) {
tmp = t_2;
} else if (t < 165.68027943805222) {
tmp = (((18.0 * y) * (x * (z * t))) - t_1) + ((c * b) - (27.0 * (k * j)));
} else {
tmp = t_2;
}
return tmp;
}
def code(x, y, z, t, a, b, c, i, j, k): t_1 = ((a * t) + (i * x)) * 4.0 t_2 = (((18.0 * t) * ((x * y) * z)) - t_1) - (((k * j) * 27.0) - (c * b)) tmp = 0 if t < -1.6210815397541398e-69: tmp = t_2 elif t < 165.68027943805222: tmp = (((18.0 * y) * (x * (z * t))) - t_1) + ((c * b) - (27.0 * (k * j))) else: tmp = t_2 return tmp
function code(x, y, z, t, a, b, c, i, j, k) t_1 = Float64(Float64(Float64(a * t) + Float64(i * x)) * 4.0) t_2 = Float64(Float64(Float64(Float64(18.0 * t) * Float64(Float64(x * y) * z)) - t_1) - Float64(Float64(Float64(k * j) * 27.0) - Float64(c * b))) tmp = 0.0 if (t < -1.6210815397541398e-69) tmp = t_2; elseif (t < 165.68027943805222) tmp = Float64(Float64(Float64(Float64(18.0 * y) * Float64(x * Float64(z * t))) - t_1) + Float64(Float64(c * b) - Float64(27.0 * Float64(k * j)))); else tmp = t_2; end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c, i, j, k) t_1 = ((a * t) + (i * x)) * 4.0; t_2 = (((18.0 * t) * ((x * y) * z)) - t_1) - (((k * j) * 27.0) - (c * b)); tmp = 0.0; if (t < -1.6210815397541398e-69) tmp = t_2; elseif (t < 165.68027943805222) tmp = (((18.0 * y) * (x * (z * t))) - t_1) + ((c * b) - (27.0 * (k * j))); else tmp = t_2; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_, i_, j_, k_] := Block[{t$95$1 = N[(N[(N[(a * t), $MachinePrecision] + N[(i * x), $MachinePrecision]), $MachinePrecision] * 4.0), $MachinePrecision]}, Block[{t$95$2 = N[(N[(N[(N[(18.0 * t), $MachinePrecision] * N[(N[(x * y), $MachinePrecision] * z), $MachinePrecision]), $MachinePrecision] - t$95$1), $MachinePrecision] - N[(N[(N[(k * j), $MachinePrecision] * 27.0), $MachinePrecision] - N[(c * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[Less[t, -1.6210815397541398e-69], t$95$2, If[Less[t, 165.68027943805222], N[(N[(N[(N[(18.0 * y), $MachinePrecision] * N[(x * N[(z * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - t$95$1), $MachinePrecision] + N[(N[(c * b), $MachinePrecision] - N[(27.0 * N[(k * j), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$2]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \left(a \cdot t + i \cdot x\right) \cdot 4\\
t_2 := \left(\left(18 \cdot t\right) \cdot \left(\left(x \cdot y\right) \cdot z\right) - t\_1\right) - \left(\left(k \cdot j\right) \cdot 27 - c \cdot b\right)\\
\mathbf{if}\;t < -1.6210815397541398 \cdot 10^{-69}:\\
\;\;\;\;t\_2\\
\mathbf{elif}\;t < 165.68027943805222:\\
\;\;\;\;\left(\left(18 \cdot y\right) \cdot \left(x \cdot \left(z \cdot t\right)\right) - t\_1\right) + \left(c \cdot b - 27 \cdot \left(k \cdot j\right)\right)\\
\mathbf{else}:\\
\;\;\;\;t\_2\\
\end{array}
\end{array}
herbie shell --seed 2024363
(FPCore (x y z t a b c i j k)
:name "Diagrams.Solve.Polynomial:cubForm from diagrams-solve-0.1, E"
:precision binary64
:alt
(! :herbie-platform default (if (< t -8105407698770699/5000000000000000000000000000000000000000000000000000000000000000000000000000000000000) (- (- (* (* 18 t) (* (* x y) z)) (* (+ (* a t) (* i x)) 4)) (- (* (* k j) 27) (* c b))) (if (< t 8284013971902611/50000000000000) (+ (- (* (* 18 y) (* x (* z t))) (* (+ (* a t) (* i x)) 4)) (- (* c b) (* 27 (* k j)))) (- (- (* (* 18 t) (* (* x y) z)) (* (+ (* a t) (* i x)) 4)) (- (* (* k j) 27) (* c b))))))
(- (- (+ (- (* (* (* (* x 18.0) y) z) t) (* (* a 4.0) t)) (* b c)) (* (* x 4.0) i)) (* (* j 27.0) k)))