
(FPCore (x y z t a b) :precision binary64 (+ (- (- x (* (- y 1.0) z)) (* (- t 1.0) a)) (* (- (+ y t) 2.0) b)))
double code(double x, double y, double z, double t, double a, double b) {
return ((x - ((y - 1.0) * z)) - ((t - 1.0) * a)) + (((y + t) - 2.0) * b);
}
real(8) function code(x, y, z, t, a, b)
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
code = ((x - ((y - 1.0d0) * z)) - ((t - 1.0d0) * a)) + (((y + t) - 2.0d0) * b)
end function
public static double code(double x, double y, double z, double t, double a, double b) {
return ((x - ((y - 1.0) * z)) - ((t - 1.0) * a)) + (((y + t) - 2.0) * b);
}
def code(x, y, z, t, a, b): return ((x - ((y - 1.0) * z)) - ((t - 1.0) * a)) + (((y + t) - 2.0) * b)
function code(x, y, z, t, a, b) return Float64(Float64(Float64(x - Float64(Float64(y - 1.0) * z)) - Float64(Float64(t - 1.0) * a)) + Float64(Float64(Float64(y + t) - 2.0) * b)) end
function tmp = code(x, y, z, t, a, b) tmp = ((x - ((y - 1.0) * z)) - ((t - 1.0) * a)) + (((y + t) - 2.0) * b); end
code[x_, y_, z_, t_, a_, b_] := N[(N[(N[(x - N[(N[(y - 1.0), $MachinePrecision] * z), $MachinePrecision]), $MachinePrecision] - N[(N[(t - 1.0), $MachinePrecision] * a), $MachinePrecision]), $MachinePrecision] + N[(N[(N[(y + t), $MachinePrecision] - 2.0), $MachinePrecision] * b), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\left(\left(x - \left(y - 1\right) \cdot z\right) - \left(t - 1\right) \cdot a\right) + \left(\left(y + t\right) - 2\right) \cdot b
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 29 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (x y z t a b) :precision binary64 (+ (- (- x (* (- y 1.0) z)) (* (- t 1.0) a)) (* (- (+ y t) 2.0) b)))
double code(double x, double y, double z, double t, double a, double b) {
return ((x - ((y - 1.0) * z)) - ((t - 1.0) * a)) + (((y + t) - 2.0) * b);
}
real(8) function code(x, y, z, t, a, b)
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
code = ((x - ((y - 1.0d0) * z)) - ((t - 1.0d0) * a)) + (((y + t) - 2.0d0) * b)
end function
public static double code(double x, double y, double z, double t, double a, double b) {
return ((x - ((y - 1.0) * z)) - ((t - 1.0) * a)) + (((y + t) - 2.0) * b);
}
def code(x, y, z, t, a, b): return ((x - ((y - 1.0) * z)) - ((t - 1.0) * a)) + (((y + t) - 2.0) * b)
function code(x, y, z, t, a, b) return Float64(Float64(Float64(x - Float64(Float64(y - 1.0) * z)) - Float64(Float64(t - 1.0) * a)) + Float64(Float64(Float64(y + t) - 2.0) * b)) end
function tmp = code(x, y, z, t, a, b) tmp = ((x - ((y - 1.0) * z)) - ((t - 1.0) * a)) + (((y + t) - 2.0) * b); end
code[x_, y_, z_, t_, a_, b_] := N[(N[(N[(x - N[(N[(y - 1.0), $MachinePrecision] * z), $MachinePrecision]), $MachinePrecision] - N[(N[(t - 1.0), $MachinePrecision] * a), $MachinePrecision]), $MachinePrecision] + N[(N[(N[(y + t), $MachinePrecision] - 2.0), $MachinePrecision] * b), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\left(\left(x - \left(y - 1\right) \cdot z\right) - \left(t - 1\right) \cdot a\right) + \left(\left(y + t\right) - 2\right) \cdot b
\end{array}
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (* a (- t 1.0))) (t_2 (* (- (+ t y) 2.0) b)))
(if (<= (+ t_2 (- (- x (* z (- y 1.0))) t_1)) INFINITY)
(- x (fma z (- y 1.0) (- t_1 t_2)))
(- x (fma z (- y 1.0) (* (- a b) t))))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = a * (t - 1.0);
double t_2 = ((t + y) - 2.0) * b;
double tmp;
if ((t_2 + ((x - (z * (y - 1.0))) - t_1)) <= ((double) INFINITY)) {
tmp = x - fma(z, (y - 1.0), (t_1 - t_2));
} else {
tmp = x - fma(z, (y - 1.0), ((a - b) * t));
}
return tmp;
}
function code(x, y, z, t, a, b) t_1 = Float64(a * Float64(t - 1.0)) t_2 = Float64(Float64(Float64(t + y) - 2.0) * b) tmp = 0.0 if (Float64(t_2 + Float64(Float64(x - Float64(z * Float64(y - 1.0))) - t_1)) <= Inf) tmp = Float64(x - fma(z, Float64(y - 1.0), Float64(t_1 - t_2))); else tmp = Float64(x - fma(z, Float64(y - 1.0), Float64(Float64(a - b) * t))); end return tmp end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(a * N[(t - 1.0), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(N[(N[(t + y), $MachinePrecision] - 2.0), $MachinePrecision] * b), $MachinePrecision]}, If[LessEqual[N[(t$95$2 + N[(N[(x - N[(z * N[(y - 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - t$95$1), $MachinePrecision]), $MachinePrecision], Infinity], N[(x - N[(z * N[(y - 1.0), $MachinePrecision] + N[(t$95$1 - t$95$2), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x - N[(z * N[(y - 1.0), $MachinePrecision] + N[(N[(a - b), $MachinePrecision] * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := a \cdot \left(t - 1\right)\\
t_2 := \left(\left(t + y\right) - 2\right) \cdot b\\
\mathbf{if}\;t\_2 + \left(\left(x - z \cdot \left(y - 1\right)\right) - t\_1\right) \leq \infty:\\
\;\;\;\;x - \mathsf{fma}\left(z, y - 1, t\_1 - t\_2\right)\\
\mathbf{else}:\\
\;\;\;\;x - \mathsf{fma}\left(z, y - 1, \left(a - b\right) \cdot t\right)\\
\end{array}
\end{array}
if (+.f64 (-.f64 (-.f64 x (*.f64 (-.f64 y #s(literal 1 binary64)) z)) (*.f64 (-.f64 t #s(literal 1 binary64)) a)) (*.f64 (-.f64 (+.f64 y t) #s(literal 2 binary64)) b)) < +inf.0Initial program 100.0%
lift-+.f64N/A
lift--.f64N/A
associate-+l-N/A
lift--.f64N/A
associate--l-N/A
lower--.f64N/A
lift-*.f64N/A
*-commutativeN/A
lower-fma.f64N/A
lower--.f64100.0
lift-*.f64N/A
*-commutativeN/A
lower-*.f64100.0
lift-*.f64N/A
*-commutativeN/A
lower-*.f64100.0
lift-+.f64N/A
+-commutativeN/A
lower-+.f64100.0
Applied rewrites100.0%
if +inf.0 < (+.f64 (-.f64 (-.f64 x (*.f64 (-.f64 y #s(literal 1 binary64)) z)) (*.f64 (-.f64 t #s(literal 1 binary64)) a)) (*.f64 (-.f64 (+.f64 y t) #s(literal 2 binary64)) b)) Initial program 0.0%
lift-+.f64N/A
lift--.f64N/A
associate-+l-N/A
lift--.f64N/A
associate--l-N/A
lower--.f64N/A
lift-*.f64N/A
*-commutativeN/A
lower-fma.f64N/A
lower--.f648.3
lift-*.f64N/A
*-commutativeN/A
lower-*.f648.3
lift-*.f64N/A
*-commutativeN/A
lower-*.f648.3
lift-+.f64N/A
+-commutativeN/A
lower-+.f648.3
Applied rewrites8.3%
Taylor expanded in t around inf
*-commutativeN/A
sub-negN/A
mul-1-negN/A
lower-*.f64N/A
mul-1-negN/A
sub-negN/A
lower--.f6483.3
Applied rewrites83.3%
Final simplification99.2%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (* (- b a) t)) (t_2 (* (- 1.0 y) z)))
(if (<= t -2.3e+18)
t_1
(if (<= t -1.12e-203)
(* (- b z) y)
(if (<= t 5.25e-263)
(+ (* b y) a)
(if (<= t 5.4e-157)
t_2
(if (<= t 1.32e-89)
(* (- y 2.0) b)
(if (<= t 6.2e+26) t_2 t_1))))))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = (b - a) * t;
double t_2 = (1.0 - y) * z;
double tmp;
if (t <= -2.3e+18) {
tmp = t_1;
} else if (t <= -1.12e-203) {
tmp = (b - z) * y;
} else if (t <= 5.25e-263) {
tmp = (b * y) + a;
} else if (t <= 5.4e-157) {
tmp = t_2;
} else if (t <= 1.32e-89) {
tmp = (y - 2.0) * b;
} else if (t <= 6.2e+26) {
tmp = t_2;
} else {
tmp = t_1;
}
return tmp;
}
real(8) function code(x, y, z, t, a, b)
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) :: t_1
real(8) :: t_2
real(8) :: tmp
t_1 = (b - a) * t
t_2 = (1.0d0 - y) * z
if (t <= (-2.3d+18)) then
tmp = t_1
else if (t <= (-1.12d-203)) then
tmp = (b - z) * y
else if (t <= 5.25d-263) then
tmp = (b * y) + a
else if (t <= 5.4d-157) then
tmp = t_2
else if (t <= 1.32d-89) then
tmp = (y - 2.0d0) * b
else if (t <= 6.2d+26) then
tmp = t_2
else
tmp = t_1
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b) {
double t_1 = (b - a) * t;
double t_2 = (1.0 - y) * z;
double tmp;
if (t <= -2.3e+18) {
tmp = t_1;
} else if (t <= -1.12e-203) {
tmp = (b - z) * y;
} else if (t <= 5.25e-263) {
tmp = (b * y) + a;
} else if (t <= 5.4e-157) {
tmp = t_2;
} else if (t <= 1.32e-89) {
tmp = (y - 2.0) * b;
} else if (t <= 6.2e+26) {
tmp = t_2;
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = (b - a) * t t_2 = (1.0 - y) * z tmp = 0 if t <= -2.3e+18: tmp = t_1 elif t <= -1.12e-203: tmp = (b - z) * y elif t <= 5.25e-263: tmp = (b * y) + a elif t <= 5.4e-157: tmp = t_2 elif t <= 1.32e-89: tmp = (y - 2.0) * b elif t <= 6.2e+26: tmp = t_2 else: tmp = t_1 return tmp
function code(x, y, z, t, a, b) t_1 = Float64(Float64(b - a) * t) t_2 = Float64(Float64(1.0 - y) * z) tmp = 0.0 if (t <= -2.3e+18) tmp = t_1; elseif (t <= -1.12e-203) tmp = Float64(Float64(b - z) * y); elseif (t <= 5.25e-263) tmp = Float64(Float64(b * y) + a); elseif (t <= 5.4e-157) tmp = t_2; elseif (t <= 1.32e-89) tmp = Float64(Float64(y - 2.0) * b); elseif (t <= 6.2e+26) tmp = t_2; else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a, b) t_1 = (b - a) * t; t_2 = (1.0 - y) * z; tmp = 0.0; if (t <= -2.3e+18) tmp = t_1; elseif (t <= -1.12e-203) tmp = (b - z) * y; elseif (t <= 5.25e-263) tmp = (b * y) + a; elseif (t <= 5.4e-157) tmp = t_2; elseif (t <= 1.32e-89) tmp = (y - 2.0) * b; elseif (t <= 6.2e+26) tmp = t_2; else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(N[(b - a), $MachinePrecision] * t), $MachinePrecision]}, Block[{t$95$2 = N[(N[(1.0 - y), $MachinePrecision] * z), $MachinePrecision]}, If[LessEqual[t, -2.3e+18], t$95$1, If[LessEqual[t, -1.12e-203], N[(N[(b - z), $MachinePrecision] * y), $MachinePrecision], If[LessEqual[t, 5.25e-263], N[(N[(b * y), $MachinePrecision] + a), $MachinePrecision], If[LessEqual[t, 5.4e-157], t$95$2, If[LessEqual[t, 1.32e-89], N[(N[(y - 2.0), $MachinePrecision] * b), $MachinePrecision], If[LessEqual[t, 6.2e+26], t$95$2, t$95$1]]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \left(b - a\right) \cdot t\\
t_2 := \left(1 - y\right) \cdot z\\
\mathbf{if}\;t \leq -2.3 \cdot 10^{+18}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;t \leq -1.12 \cdot 10^{-203}:\\
\;\;\;\;\left(b - z\right) \cdot y\\
\mathbf{elif}\;t \leq 5.25 \cdot 10^{-263}:\\
\;\;\;\;b \cdot y + a\\
\mathbf{elif}\;t \leq 5.4 \cdot 10^{-157}:\\
\;\;\;\;t\_2\\
\mathbf{elif}\;t \leq 1.32 \cdot 10^{-89}:\\
\;\;\;\;\left(y - 2\right) \cdot b\\
\mathbf{elif}\;t \leq 6.2 \cdot 10^{+26}:\\
\;\;\;\;t\_2\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if t < -2.3e18 or 6.1999999999999999e26 < t Initial program 91.9%
Taylor expanded in t around inf
*-commutativeN/A
lower-*.f64N/A
lower--.f6470.1
Applied rewrites70.1%
if -2.3e18 < t < -1.12e-203Initial program 97.5%
Taylor expanded in y around inf
*-commutativeN/A
lower-*.f64N/A
lower--.f6449.6
Applied rewrites49.6%
if -1.12e-203 < t < 5.2499999999999999e-263Initial program 100.0%
Taylor expanded in z around 0
sub-negN/A
+-commutativeN/A
*-commutativeN/A
distribute-lft-neg-inN/A
mul-1-negN/A
lower-fma.f64N/A
mul-1-negN/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
distribute-neg-inN/A
metadata-evalN/A
sub-negN/A
lower--.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower--.f64N/A
+-commutativeN/A
lower-+.f6480.4
Applied rewrites80.4%
Taylor expanded in t around 0
Applied rewrites80.4%
Taylor expanded in y around inf
Applied rewrites59.6%
if 5.2499999999999999e-263 < t < 5.4e-157 or 1.32e-89 < t < 6.1999999999999999e26Initial program 98.0%
Taylor expanded in z around inf
*-commutativeN/A
sub-negN/A
metadata-evalN/A
distribute-neg-inN/A
+-commutativeN/A
metadata-evalN/A
sub-negN/A
mul-1-negN/A
lower-*.f64N/A
sub-negN/A
metadata-evalN/A
distribute-lft-inN/A
metadata-evalN/A
+-commutativeN/A
neg-mul-1N/A
sub-negN/A
lower--.f6450.5
Applied rewrites50.5%
if 5.4e-157 < t < 1.32e-89Initial program 100.0%
Taylor expanded in b around inf
*-commutativeN/A
lower-*.f64N/A
lower--.f64N/A
+-commutativeN/A
lower-+.f6483.0
Applied rewrites83.0%
Taylor expanded in t around 0
Applied rewrites83.0%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (+ (fma (- b a) t (fma (- y 2.0) b x)) a)))
(if (<= x -3.9e+118)
t_1
(if (<= x 1.4e+134)
(fma (- b z) y (fma (- t 2.0) b (fma (- 1.0 t) a z)))
t_1))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = fma((b - a), t, fma((y - 2.0), b, x)) + a;
double tmp;
if (x <= -3.9e+118) {
tmp = t_1;
} else if (x <= 1.4e+134) {
tmp = fma((b - z), y, fma((t - 2.0), b, fma((1.0 - t), a, z)));
} else {
tmp = t_1;
}
return tmp;
}
function code(x, y, z, t, a, b) t_1 = Float64(fma(Float64(b - a), t, fma(Float64(y - 2.0), b, x)) + a) tmp = 0.0 if (x <= -3.9e+118) tmp = t_1; elseif (x <= 1.4e+134) tmp = fma(Float64(b - z), y, fma(Float64(t - 2.0), b, fma(Float64(1.0 - t), a, z))); else tmp = t_1; end return tmp end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(N[(N[(b - a), $MachinePrecision] * t + N[(N[(y - 2.0), $MachinePrecision] * b + x), $MachinePrecision]), $MachinePrecision] + a), $MachinePrecision]}, If[LessEqual[x, -3.9e+118], t$95$1, If[LessEqual[x, 1.4e+134], N[(N[(b - z), $MachinePrecision] * y + N[(N[(t - 2.0), $MachinePrecision] * b + N[(N[(1.0 - t), $MachinePrecision] * a + z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \mathsf{fma}\left(b - a, t, \mathsf{fma}\left(y - 2, b, x\right)\right) + a\\
\mathbf{if}\;x \leq -3.9 \cdot 10^{+118}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;x \leq 1.4 \cdot 10^{+134}:\\
\;\;\;\;\mathsf{fma}\left(b - z, y, \mathsf{fma}\left(t - 2, b, \mathsf{fma}\left(1 - t, a, z\right)\right)\right)\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if x < -3.9e118 or 1.3999999999999999e134 < x Initial program 94.9%
Taylor expanded in z around 0
sub-negN/A
+-commutativeN/A
*-commutativeN/A
distribute-lft-neg-inN/A
mul-1-negN/A
lower-fma.f64N/A
mul-1-negN/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
distribute-neg-inN/A
metadata-evalN/A
sub-negN/A
lower--.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower--.f64N/A
+-commutativeN/A
lower-+.f6484.6
Applied rewrites84.6%
Taylor expanded in t around 0
Applied rewrites85.9%
if -3.9e118 < x < 1.3999999999999999e134Initial program 95.5%
Taylor expanded in x around 0
Applied rewrites97.1%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (* (- y 2.0) b)) (t_2 (* (- 1.0 y) z)) (t_3 (* (- b a) t)))
(if (<= t -75000000.0)
t_3
(if (<= t -1.02e-294)
t_1
(if (<= t 5.4e-157)
t_2
(if (<= t 1.32e-89) t_1 (if (<= t 6.2e+26) t_2 t_3)))))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = (y - 2.0) * b;
double t_2 = (1.0 - y) * z;
double t_3 = (b - a) * t;
double tmp;
if (t <= -75000000.0) {
tmp = t_3;
} else if (t <= -1.02e-294) {
tmp = t_1;
} else if (t <= 5.4e-157) {
tmp = t_2;
} else if (t <= 1.32e-89) {
tmp = t_1;
} else if (t <= 6.2e+26) {
tmp = t_2;
} else {
tmp = t_3;
}
return tmp;
}
real(8) function code(x, y, z, t, a, b)
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) :: t_1
real(8) :: t_2
real(8) :: t_3
real(8) :: tmp
t_1 = (y - 2.0d0) * b
t_2 = (1.0d0 - y) * z
t_3 = (b - a) * t
if (t <= (-75000000.0d0)) then
tmp = t_3
else if (t <= (-1.02d-294)) then
tmp = t_1
else if (t <= 5.4d-157) then
tmp = t_2
else if (t <= 1.32d-89) then
tmp = t_1
else if (t <= 6.2d+26) then
tmp = t_2
else
tmp = t_3
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b) {
double t_1 = (y - 2.0) * b;
double t_2 = (1.0 - y) * z;
double t_3 = (b - a) * t;
double tmp;
if (t <= -75000000.0) {
tmp = t_3;
} else if (t <= -1.02e-294) {
tmp = t_1;
} else if (t <= 5.4e-157) {
tmp = t_2;
} else if (t <= 1.32e-89) {
tmp = t_1;
} else if (t <= 6.2e+26) {
tmp = t_2;
} else {
tmp = t_3;
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = (y - 2.0) * b t_2 = (1.0 - y) * z t_3 = (b - a) * t tmp = 0 if t <= -75000000.0: tmp = t_3 elif t <= -1.02e-294: tmp = t_1 elif t <= 5.4e-157: tmp = t_2 elif t <= 1.32e-89: tmp = t_1 elif t <= 6.2e+26: tmp = t_2 else: tmp = t_3 return tmp
function code(x, y, z, t, a, b) t_1 = Float64(Float64(y - 2.0) * b) t_2 = Float64(Float64(1.0 - y) * z) t_3 = Float64(Float64(b - a) * t) tmp = 0.0 if (t <= -75000000.0) tmp = t_3; elseif (t <= -1.02e-294) tmp = t_1; elseif (t <= 5.4e-157) tmp = t_2; elseif (t <= 1.32e-89) tmp = t_1; elseif (t <= 6.2e+26) tmp = t_2; else tmp = t_3; end return tmp end
function tmp_2 = code(x, y, z, t, a, b) t_1 = (y - 2.0) * b; t_2 = (1.0 - y) * z; t_3 = (b - a) * t; tmp = 0.0; if (t <= -75000000.0) tmp = t_3; elseif (t <= -1.02e-294) tmp = t_1; elseif (t <= 5.4e-157) tmp = t_2; elseif (t <= 1.32e-89) tmp = t_1; elseif (t <= 6.2e+26) tmp = t_2; else tmp = t_3; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(N[(y - 2.0), $MachinePrecision] * b), $MachinePrecision]}, Block[{t$95$2 = N[(N[(1.0 - y), $MachinePrecision] * z), $MachinePrecision]}, Block[{t$95$3 = N[(N[(b - a), $MachinePrecision] * t), $MachinePrecision]}, If[LessEqual[t, -75000000.0], t$95$3, If[LessEqual[t, -1.02e-294], t$95$1, If[LessEqual[t, 5.4e-157], t$95$2, If[LessEqual[t, 1.32e-89], t$95$1, If[LessEqual[t, 6.2e+26], t$95$2, t$95$3]]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \left(y - 2\right) \cdot b\\
t_2 := \left(1 - y\right) \cdot z\\
t_3 := \left(b - a\right) \cdot t\\
\mathbf{if}\;t \leq -75000000:\\
\;\;\;\;t\_3\\
\mathbf{elif}\;t \leq -1.02 \cdot 10^{-294}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;t \leq 5.4 \cdot 10^{-157}:\\
\;\;\;\;t\_2\\
\mathbf{elif}\;t \leq 1.32 \cdot 10^{-89}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;t \leq 6.2 \cdot 10^{+26}:\\
\;\;\;\;t\_2\\
\mathbf{else}:\\
\;\;\;\;t\_3\\
\end{array}
\end{array}
if t < -7.5e7 or 6.1999999999999999e26 < t Initial program 92.0%
Taylor expanded in t around inf
*-commutativeN/A
lower-*.f64N/A
lower--.f6469.6
Applied rewrites69.6%
if -7.5e7 < t < -1.01999999999999998e-294 or 5.4e-157 < t < 1.32e-89Initial program 98.5%
Taylor expanded in b around inf
*-commutativeN/A
lower-*.f64N/A
lower--.f64N/A
+-commutativeN/A
lower-+.f6455.5
Applied rewrites55.5%
Taylor expanded in t around 0
Applied rewrites54.8%
if -1.01999999999999998e-294 < t < 5.4e-157 or 1.32e-89 < t < 6.1999999999999999e26Initial program 98.4%
Taylor expanded in z around inf
*-commutativeN/A
sub-negN/A
metadata-evalN/A
distribute-neg-inN/A
+-commutativeN/A
metadata-evalN/A
sub-negN/A
mul-1-negN/A
lower-*.f64N/A
sub-negN/A
metadata-evalN/A
distribute-lft-inN/A
metadata-evalN/A
+-commutativeN/A
neg-mul-1N/A
sub-negN/A
lower--.f6448.3
Applied rewrites48.3%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (fma (- 1.0 t) a x)) (t_2 (* (- (+ t y) 2.0) b)))
(if (<= b -9600000.0)
t_2
(if (<= b -3.4e-161)
t_1
(if (<= b 1.5e-274) (* (- 1.0 y) z) (if (<= b 410000000.0) t_1 t_2))))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = fma((1.0 - t), a, x);
double t_2 = ((t + y) - 2.0) * b;
double tmp;
if (b <= -9600000.0) {
tmp = t_2;
} else if (b <= -3.4e-161) {
tmp = t_1;
} else if (b <= 1.5e-274) {
tmp = (1.0 - y) * z;
} else if (b <= 410000000.0) {
tmp = t_1;
} else {
tmp = t_2;
}
return tmp;
}
function code(x, y, z, t, a, b) t_1 = fma(Float64(1.0 - t), a, x) t_2 = Float64(Float64(Float64(t + y) - 2.0) * b) tmp = 0.0 if (b <= -9600000.0) tmp = t_2; elseif (b <= -3.4e-161) tmp = t_1; elseif (b <= 1.5e-274) tmp = Float64(Float64(1.0 - y) * z); elseif (b <= 410000000.0) tmp = t_1; else tmp = t_2; end return tmp end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(N[(1.0 - t), $MachinePrecision] * a + x), $MachinePrecision]}, Block[{t$95$2 = N[(N[(N[(t + y), $MachinePrecision] - 2.0), $MachinePrecision] * b), $MachinePrecision]}, If[LessEqual[b, -9600000.0], t$95$2, If[LessEqual[b, -3.4e-161], t$95$1, If[LessEqual[b, 1.5e-274], N[(N[(1.0 - y), $MachinePrecision] * z), $MachinePrecision], If[LessEqual[b, 410000000.0], t$95$1, t$95$2]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \mathsf{fma}\left(1 - t, a, x\right)\\
t_2 := \left(\left(t + y\right) - 2\right) \cdot b\\
\mathbf{if}\;b \leq -9600000:\\
\;\;\;\;t\_2\\
\mathbf{elif}\;b \leq -3.4 \cdot 10^{-161}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;b \leq 1.5 \cdot 10^{-274}:\\
\;\;\;\;\left(1 - y\right) \cdot z\\
\mathbf{elif}\;b \leq 410000000:\\
\;\;\;\;t\_1\\
\mathbf{else}:\\
\;\;\;\;t\_2\\
\end{array}
\end{array}
if b < -9.6e6 or 4.1e8 < b Initial program 91.5%
Taylor expanded in b around inf
*-commutativeN/A
lower-*.f64N/A
lower--.f64N/A
+-commutativeN/A
lower-+.f6472.1
Applied rewrites72.1%
if -9.6e6 < b < -3.39999999999999982e-161 or 1.49999999999999989e-274 < b < 4.1e8Initial program 98.9%
Taylor expanded in z around 0
sub-negN/A
+-commutativeN/A
*-commutativeN/A
distribute-lft-neg-inN/A
mul-1-negN/A
lower-fma.f64N/A
mul-1-negN/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
distribute-neg-inN/A
metadata-evalN/A
sub-negN/A
lower--.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower--.f64N/A
+-commutativeN/A
lower-+.f6465.1
Applied rewrites65.1%
Taylor expanded in b around 0
Applied rewrites54.4%
if -3.39999999999999982e-161 < b < 1.49999999999999989e-274Initial program 100.0%
Taylor expanded in z around inf
*-commutativeN/A
sub-negN/A
metadata-evalN/A
distribute-neg-inN/A
+-commutativeN/A
metadata-evalN/A
sub-negN/A
mul-1-negN/A
lower-*.f64N/A
sub-negN/A
metadata-evalN/A
distribute-lft-inN/A
metadata-evalN/A
+-commutativeN/A
neg-mul-1N/A
sub-negN/A
lower--.f6471.4
Applied rewrites71.4%
Final simplification65.4%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (fma (- (+ t y) 2.0) b x)))
(if (<= b -33000000000.0)
t_1
(if (<= b 2.8e-7)
(fma (- 1.0 y) z (fma (- 1.0 t) a x))
(if (<= b 6.5e+203) (+ (fma (- b a) t (* b y)) a) t_1)))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = fma(((t + y) - 2.0), b, x);
double tmp;
if (b <= -33000000000.0) {
tmp = t_1;
} else if (b <= 2.8e-7) {
tmp = fma((1.0 - y), z, fma((1.0 - t), a, x));
} else if (b <= 6.5e+203) {
tmp = fma((b - a), t, (b * y)) + a;
} else {
tmp = t_1;
}
return tmp;
}
function code(x, y, z, t, a, b) t_1 = fma(Float64(Float64(t + y) - 2.0), b, x) tmp = 0.0 if (b <= -33000000000.0) tmp = t_1; elseif (b <= 2.8e-7) tmp = fma(Float64(1.0 - y), z, fma(Float64(1.0 - t), a, x)); elseif (b <= 6.5e+203) tmp = Float64(fma(Float64(b - a), t, Float64(b * y)) + a); else tmp = t_1; end return tmp end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(N[(N[(t + y), $MachinePrecision] - 2.0), $MachinePrecision] * b + x), $MachinePrecision]}, If[LessEqual[b, -33000000000.0], t$95$1, If[LessEqual[b, 2.8e-7], N[(N[(1.0 - y), $MachinePrecision] * z + N[(N[(1.0 - t), $MachinePrecision] * a + x), $MachinePrecision]), $MachinePrecision], If[LessEqual[b, 6.5e+203], N[(N[(N[(b - a), $MachinePrecision] * t + N[(b * y), $MachinePrecision]), $MachinePrecision] + a), $MachinePrecision], t$95$1]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \mathsf{fma}\left(\left(t + y\right) - 2, b, x\right)\\
\mathbf{if}\;b \leq -33000000000:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;b \leq 2.8 \cdot 10^{-7}:\\
\;\;\;\;\mathsf{fma}\left(1 - y, z, \mathsf{fma}\left(1 - t, a, x\right)\right)\\
\mathbf{elif}\;b \leq 6.5 \cdot 10^{+203}:\\
\;\;\;\;\mathsf{fma}\left(b - a, t, b \cdot y\right) + a\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if b < -3.3e10 or 6.5000000000000003e203 < b Initial program 90.9%
Taylor expanded in z around 0
sub-negN/A
+-commutativeN/A
*-commutativeN/A
distribute-lft-neg-inN/A
mul-1-negN/A
lower-fma.f64N/A
mul-1-negN/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
distribute-neg-inN/A
metadata-evalN/A
sub-negN/A
lower--.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower--.f64N/A
+-commutativeN/A
lower-+.f6491.4
Applied rewrites91.4%
Taylor expanded in t around 0
Applied rewrites89.1%
Taylor expanded in a around 0
Applied rewrites88.9%
if -3.3e10 < b < 2.80000000000000019e-7Initial program 99.2%
Taylor expanded in b around 0
+-commutativeN/A
associate--r+N/A
sub-negN/A
+-commutativeN/A
associate-+r-N/A
*-commutativeN/A
distribute-lft-neg-inN/A
mul-1-negN/A
lower-fma.f64N/A
sub-negN/A
metadata-evalN/A
distribute-lft-inN/A
metadata-evalN/A
+-commutativeN/A
neg-mul-1N/A
sub-negN/A
lower--.f64N/A
sub-negN/A
+-commutativeN/A
Applied rewrites91.5%
if 2.80000000000000019e-7 < b < 6.5000000000000003e203Initial program 93.4%
Taylor expanded in z around 0
sub-negN/A
+-commutativeN/A
*-commutativeN/A
distribute-lft-neg-inN/A
mul-1-negN/A
lower-fma.f64N/A
mul-1-negN/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
distribute-neg-inN/A
metadata-evalN/A
sub-negN/A
lower--.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower--.f64N/A
+-commutativeN/A
lower-+.f6485.3
Applied rewrites85.3%
Taylor expanded in t around 0
Applied rewrites89.6%
Taylor expanded in y around inf
Applied rewrites80.8%
Final simplification88.7%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (- x (* (- z b) y))))
(if (<= y -1.6e+103)
t_1
(if (<= y -245000000000.0)
(+ (fma (- b a) t (* b y)) a)
(if (<= y 2e+114) (- x (fma z -1.0 (* (- a b) t))) t_1)))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = x - ((z - b) * y);
double tmp;
if (y <= -1.6e+103) {
tmp = t_1;
} else if (y <= -245000000000.0) {
tmp = fma((b - a), t, (b * y)) + a;
} else if (y <= 2e+114) {
tmp = x - fma(z, -1.0, ((a - b) * t));
} else {
tmp = t_1;
}
return tmp;
}
function code(x, y, z, t, a, b) t_1 = Float64(x - Float64(Float64(z - b) * y)) tmp = 0.0 if (y <= -1.6e+103) tmp = t_1; elseif (y <= -245000000000.0) tmp = Float64(fma(Float64(b - a), t, Float64(b * y)) + a); elseif (y <= 2e+114) tmp = Float64(x - fma(z, -1.0, Float64(Float64(a - b) * t))); else tmp = t_1; end return tmp end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(x - N[(N[(z - b), $MachinePrecision] * y), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -1.6e+103], t$95$1, If[LessEqual[y, -245000000000.0], N[(N[(N[(b - a), $MachinePrecision] * t + N[(b * y), $MachinePrecision]), $MachinePrecision] + a), $MachinePrecision], If[LessEqual[y, 2e+114], N[(x - N[(z * -1.0 + N[(N[(a - b), $MachinePrecision] * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$1]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x - \left(z - b\right) \cdot y\\
\mathbf{if}\;y \leq -1.6 \cdot 10^{+103}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y \leq -245000000000:\\
\;\;\;\;\mathsf{fma}\left(b - a, t, b \cdot y\right) + a\\
\mathbf{elif}\;y \leq 2 \cdot 10^{+114}:\\
\;\;\;\;x - \mathsf{fma}\left(z, -1, \left(a - b\right) \cdot t\right)\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if y < -1.59999999999999996e103 or 2e114 < y Initial program 91.4%
lift-+.f64N/A
lift--.f64N/A
associate-+l-N/A
lift--.f64N/A
associate--l-N/A
lower--.f64N/A
lift-*.f64N/A
*-commutativeN/A
lower-fma.f64N/A
lower--.f6492.7
lift-*.f64N/A
*-commutativeN/A
lower-*.f6492.7
lift-*.f64N/A
*-commutativeN/A
lower-*.f6492.7
lift-+.f64N/A
+-commutativeN/A
lower-+.f6492.7
Applied rewrites92.7%
Taylor expanded in y around inf
*-commutativeN/A
sub-negN/A
mul-1-negN/A
+-commutativeN/A
*-lft-identityN/A
metadata-evalN/A
cancel-sign-sub-invN/A
lower-*.f64N/A
cancel-sign-sub-invN/A
metadata-evalN/A
*-lft-identityN/A
+-commutativeN/A
mul-1-negN/A
sub-negN/A
lower--.f6487.5
Applied rewrites87.5%
if -1.59999999999999996e103 < y < -2.45e11Initial program 95.8%
Taylor expanded in z around 0
sub-negN/A
+-commutativeN/A
*-commutativeN/A
distribute-lft-neg-inN/A
mul-1-negN/A
lower-fma.f64N/A
mul-1-negN/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
distribute-neg-inN/A
metadata-evalN/A
sub-negN/A
lower--.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower--.f64N/A
+-commutativeN/A
lower-+.f6474.9
Applied rewrites74.9%
Taylor expanded in t around 0
Applied rewrites74.9%
Taylor expanded in y around inf
Applied rewrites74.9%
if -2.45e11 < y < 2e114Initial program 97.3%
lift-+.f64N/A
lift--.f64N/A
associate-+l-N/A
lift--.f64N/A
associate--l-N/A
lower--.f64N/A
lift-*.f64N/A
*-commutativeN/A
lower-fma.f64N/A
lower--.f6497.3
lift-*.f64N/A
*-commutativeN/A
lower-*.f6497.3
lift-*.f64N/A
*-commutativeN/A
lower-*.f6497.3
lift-+.f64N/A
+-commutativeN/A
lower-+.f6497.3
Applied rewrites97.3%
Taylor expanded in t around inf
*-commutativeN/A
sub-negN/A
mul-1-negN/A
lower-*.f64N/A
mul-1-negN/A
sub-negN/A
lower--.f6477.5
Applied rewrites77.5%
Taylor expanded in y around 0
Applied rewrites71.4%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (fma (- (+ t y) 2.0) b x)))
(if (<= b -680000000.0)
t_1
(if (<= b 5.2e-99)
(- x (* z (- y 1.0)))
(if (<= b 6.5e+203) (+ (fma (- b a) t (* b y)) a) t_1)))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = fma(((t + y) - 2.0), b, x);
double tmp;
if (b <= -680000000.0) {
tmp = t_1;
} else if (b <= 5.2e-99) {
tmp = x - (z * (y - 1.0));
} else if (b <= 6.5e+203) {
tmp = fma((b - a), t, (b * y)) + a;
} else {
tmp = t_1;
}
return tmp;
}
function code(x, y, z, t, a, b) t_1 = fma(Float64(Float64(t + y) - 2.0), b, x) tmp = 0.0 if (b <= -680000000.0) tmp = t_1; elseif (b <= 5.2e-99) tmp = Float64(x - Float64(z * Float64(y - 1.0))); elseif (b <= 6.5e+203) tmp = Float64(fma(Float64(b - a), t, Float64(b * y)) + a); else tmp = t_1; end return tmp end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(N[(N[(t + y), $MachinePrecision] - 2.0), $MachinePrecision] * b + x), $MachinePrecision]}, If[LessEqual[b, -680000000.0], t$95$1, If[LessEqual[b, 5.2e-99], N[(x - N[(z * N[(y - 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[b, 6.5e+203], N[(N[(N[(b - a), $MachinePrecision] * t + N[(b * y), $MachinePrecision]), $MachinePrecision] + a), $MachinePrecision], t$95$1]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \mathsf{fma}\left(\left(t + y\right) - 2, b, x\right)\\
\mathbf{if}\;b \leq -680000000:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;b \leq 5.2 \cdot 10^{-99}:\\
\;\;\;\;x - z \cdot \left(y - 1\right)\\
\mathbf{elif}\;b \leq 6.5 \cdot 10^{+203}:\\
\;\;\;\;\mathsf{fma}\left(b - a, t, b \cdot y\right) + a\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if b < -6.8e8 or 6.5000000000000003e203 < b Initial program 90.9%
Taylor expanded in z around 0
sub-negN/A
+-commutativeN/A
*-commutativeN/A
distribute-lft-neg-inN/A
mul-1-negN/A
lower-fma.f64N/A
mul-1-negN/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
distribute-neg-inN/A
metadata-evalN/A
sub-negN/A
lower--.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower--.f64N/A
+-commutativeN/A
lower-+.f6491.4
Applied rewrites91.4%
Taylor expanded in t around 0
Applied rewrites89.1%
Taylor expanded in a around 0
Applied rewrites88.9%
if -6.8e8 < b < 5.2000000000000001e-99Initial program 100.0%
lift-+.f64N/A
lift--.f64N/A
associate-+l-N/A
lift--.f64N/A
associate--l-N/A
lower--.f64N/A
lift-*.f64N/A
*-commutativeN/A
lower-fma.f64N/A
lower--.f64100.0
lift-*.f64N/A
*-commutativeN/A
lower-*.f64100.0
lift-*.f64N/A
*-commutativeN/A
lower-*.f64100.0
lift-+.f64N/A
+-commutativeN/A
lower-+.f64100.0
Applied rewrites100.0%
Taylor expanded in z around inf
*-commutativeN/A
lower-*.f64N/A
lower--.f6466.4
Applied rewrites66.4%
if 5.2000000000000001e-99 < b < 6.5000000000000003e203Initial program 94.2%
Taylor expanded in z around 0
sub-negN/A
+-commutativeN/A
*-commutativeN/A
distribute-lft-neg-inN/A
mul-1-negN/A
lower-fma.f64N/A
mul-1-negN/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
distribute-neg-inN/A
metadata-evalN/A
sub-negN/A
lower--.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower--.f64N/A
+-commutativeN/A
lower-+.f6478.3
Applied rewrites78.3%
Taylor expanded in t around 0
Applied rewrites81.2%
Taylor expanded in y around inf
Applied rewrites71.1%
Final simplification75.4%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (+ (fma (- b a) t (fma (- y 2.0) b x)) a)))
(if (<= a -6000.0)
t_1
(if (<= a 2e+158) (fma (- b z) y (fma (- t 2.0) b (+ z x))) t_1))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = fma((b - a), t, fma((y - 2.0), b, x)) + a;
double tmp;
if (a <= -6000.0) {
tmp = t_1;
} else if (a <= 2e+158) {
tmp = fma((b - z), y, fma((t - 2.0), b, (z + x)));
} else {
tmp = t_1;
}
return tmp;
}
function code(x, y, z, t, a, b) t_1 = Float64(fma(Float64(b - a), t, fma(Float64(y - 2.0), b, x)) + a) tmp = 0.0 if (a <= -6000.0) tmp = t_1; elseif (a <= 2e+158) tmp = fma(Float64(b - z), y, fma(Float64(t - 2.0), b, Float64(z + x))); else tmp = t_1; end return tmp end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(N[(N[(b - a), $MachinePrecision] * t + N[(N[(y - 2.0), $MachinePrecision] * b + x), $MachinePrecision]), $MachinePrecision] + a), $MachinePrecision]}, If[LessEqual[a, -6000.0], t$95$1, If[LessEqual[a, 2e+158], N[(N[(b - z), $MachinePrecision] * y + N[(N[(t - 2.0), $MachinePrecision] * b + N[(z + x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \mathsf{fma}\left(b - a, t, \mathsf{fma}\left(y - 2, b, x\right)\right) + a\\
\mathbf{if}\;a \leq -6000:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;a \leq 2 \cdot 10^{+158}:\\
\;\;\;\;\mathsf{fma}\left(b - z, y, \mathsf{fma}\left(t - 2, b, z + x\right)\right)\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if a < -6e3 or 1.99999999999999991e158 < a Initial program 93.0%
Taylor expanded in z around 0
sub-negN/A
+-commutativeN/A
*-commutativeN/A
distribute-lft-neg-inN/A
mul-1-negN/A
lower-fma.f64N/A
mul-1-negN/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
distribute-neg-inN/A
metadata-evalN/A
sub-negN/A
lower--.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower--.f64N/A
+-commutativeN/A
lower-+.f6486.7
Applied rewrites86.7%
Taylor expanded in t around 0
Applied rewrites87.7%
if -6e3 < a < 1.99999999999999991e158Initial program 96.8%
Taylor expanded in a around 0
Applied rewrites95.0%
(FPCore (x y z t a b)
:precision binary64
(if (<= b -680000000.0)
(fma (- 1.0 t) a (fma (- (+ t y) 2.0) b x))
(if (<= b 2.2e-7)
(fma (- 1.0 y) z (fma (- 1.0 t) a x))
(+ (fma (- b a) t (fma (- y 2.0) b x)) a))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (b <= -680000000.0) {
tmp = fma((1.0 - t), a, fma(((t + y) - 2.0), b, x));
} else if (b <= 2.2e-7) {
tmp = fma((1.0 - y), z, fma((1.0 - t), a, x));
} else {
tmp = fma((b - a), t, fma((y - 2.0), b, x)) + a;
}
return tmp;
}
function code(x, y, z, t, a, b) tmp = 0.0 if (b <= -680000000.0) tmp = fma(Float64(1.0 - t), a, fma(Float64(Float64(t + y) - 2.0), b, x)); elseif (b <= 2.2e-7) tmp = fma(Float64(1.0 - y), z, fma(Float64(1.0 - t), a, x)); else tmp = Float64(fma(Float64(b - a), t, fma(Float64(y - 2.0), b, x)) + a); end return tmp end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[b, -680000000.0], N[(N[(1.0 - t), $MachinePrecision] * a + N[(N[(N[(t + y), $MachinePrecision] - 2.0), $MachinePrecision] * b + x), $MachinePrecision]), $MachinePrecision], If[LessEqual[b, 2.2e-7], N[(N[(1.0 - y), $MachinePrecision] * z + N[(N[(1.0 - t), $MachinePrecision] * a + x), $MachinePrecision]), $MachinePrecision], N[(N[(N[(b - a), $MachinePrecision] * t + N[(N[(y - 2.0), $MachinePrecision] * b + x), $MachinePrecision]), $MachinePrecision] + a), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;b \leq -680000000:\\
\;\;\;\;\mathsf{fma}\left(1 - t, a, \mathsf{fma}\left(\left(t + y\right) - 2, b, x\right)\right)\\
\mathbf{elif}\;b \leq 2.2 \cdot 10^{-7}:\\
\;\;\;\;\mathsf{fma}\left(1 - y, z, \mathsf{fma}\left(1 - t, a, x\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left(b - a, t, \mathsf{fma}\left(y - 2, b, x\right)\right) + a\\
\end{array}
\end{array}
if b < -6.8e8Initial program 93.2%
Taylor expanded in z around 0
sub-negN/A
+-commutativeN/A
*-commutativeN/A
distribute-lft-neg-inN/A
mul-1-negN/A
lower-fma.f64N/A
mul-1-negN/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
distribute-neg-inN/A
metadata-evalN/A
sub-negN/A
lower--.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower--.f64N/A
+-commutativeN/A
lower-+.f6493.8
Applied rewrites93.8%
if -6.8e8 < b < 2.2000000000000001e-7Initial program 99.2%
Taylor expanded in b around 0
+-commutativeN/A
associate--r+N/A
sub-negN/A
+-commutativeN/A
associate-+r-N/A
*-commutativeN/A
distribute-lft-neg-inN/A
mul-1-negN/A
lower-fma.f64N/A
sub-negN/A
metadata-evalN/A
distribute-lft-inN/A
metadata-evalN/A
+-commutativeN/A
neg-mul-1N/A
sub-negN/A
lower--.f64N/A
sub-negN/A
+-commutativeN/A
Applied rewrites91.5%
if 2.2000000000000001e-7 < b Initial program 90.6%
Taylor expanded in z around 0
sub-negN/A
+-commutativeN/A
*-commutativeN/A
distribute-lft-neg-inN/A
mul-1-negN/A
lower-fma.f64N/A
mul-1-negN/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
distribute-neg-inN/A
metadata-evalN/A
sub-negN/A
lower--.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower--.f64N/A
+-commutativeN/A
lower-+.f6485.8
Applied rewrites85.8%
Taylor expanded in t around 0
Applied rewrites87.1%
Final simplification90.7%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (+ (fma (- b a) t (fma (- y 2.0) b x)) a)))
(if (<= b -680000000.0)
t_1
(if (<= b 2.2e-7) (fma (- 1.0 y) z (fma (- 1.0 t) a x)) t_1))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = fma((b - a), t, fma((y - 2.0), b, x)) + a;
double tmp;
if (b <= -680000000.0) {
tmp = t_1;
} else if (b <= 2.2e-7) {
tmp = fma((1.0 - y), z, fma((1.0 - t), a, x));
} else {
tmp = t_1;
}
return tmp;
}
function code(x, y, z, t, a, b) t_1 = Float64(fma(Float64(b - a), t, fma(Float64(y - 2.0), b, x)) + a) tmp = 0.0 if (b <= -680000000.0) tmp = t_1; elseif (b <= 2.2e-7) tmp = fma(Float64(1.0 - y), z, fma(Float64(1.0 - t), a, x)); else tmp = t_1; end return tmp end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(N[(N[(b - a), $MachinePrecision] * t + N[(N[(y - 2.0), $MachinePrecision] * b + x), $MachinePrecision]), $MachinePrecision] + a), $MachinePrecision]}, If[LessEqual[b, -680000000.0], t$95$1, If[LessEqual[b, 2.2e-7], N[(N[(1.0 - y), $MachinePrecision] * z + N[(N[(1.0 - t), $MachinePrecision] * a + x), $MachinePrecision]), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \mathsf{fma}\left(b - a, t, \mathsf{fma}\left(y - 2, b, x\right)\right) + a\\
\mathbf{if}\;b \leq -680000000:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;b \leq 2.2 \cdot 10^{-7}:\\
\;\;\;\;\mathsf{fma}\left(1 - y, z, \mathsf{fma}\left(1 - t, a, x\right)\right)\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if b < -6.8e8 or 2.2000000000000001e-7 < b Initial program 91.7%
Taylor expanded in z around 0
sub-negN/A
+-commutativeN/A
*-commutativeN/A
distribute-lft-neg-inN/A
mul-1-negN/A
lower-fma.f64N/A
mul-1-negN/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
distribute-neg-inN/A
metadata-evalN/A
sub-negN/A
lower--.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower--.f64N/A
+-commutativeN/A
lower-+.f6489.3
Applied rewrites89.3%
Taylor expanded in t around 0
Applied rewrites89.3%
if -6.8e8 < b < 2.2000000000000001e-7Initial program 99.2%
Taylor expanded in b around 0
+-commutativeN/A
associate--r+N/A
sub-negN/A
+-commutativeN/A
associate-+r-N/A
*-commutativeN/A
distribute-lft-neg-inN/A
mul-1-negN/A
lower-fma.f64N/A
sub-negN/A
metadata-evalN/A
distribute-lft-inN/A
metadata-evalN/A
+-commutativeN/A
neg-mul-1N/A
sub-negN/A
lower--.f64N/A
sub-negN/A
+-commutativeN/A
Applied rewrites91.5%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (* (- y 2.0) b)))
(if (<= t -1.75e+25)
(* b t)
(if (<= t -1.08e-305)
t_1
(if (<= t 7.2e-156)
(- x (- a))
(if (<= t 620000000000.0) t_1 (* (- a) t)))))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = (y - 2.0) * b;
double tmp;
if (t <= -1.75e+25) {
tmp = b * t;
} else if (t <= -1.08e-305) {
tmp = t_1;
} else if (t <= 7.2e-156) {
tmp = x - -a;
} else if (t <= 620000000000.0) {
tmp = t_1;
} else {
tmp = -a * t;
}
return tmp;
}
real(8) function code(x, y, z, t, a, b)
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) :: t_1
real(8) :: tmp
t_1 = (y - 2.0d0) * b
if (t <= (-1.75d+25)) then
tmp = b * t
else if (t <= (-1.08d-305)) then
tmp = t_1
else if (t <= 7.2d-156) then
tmp = x - -a
else if (t <= 620000000000.0d0) then
tmp = t_1
else
tmp = -a * t
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b) {
double t_1 = (y - 2.0) * b;
double tmp;
if (t <= -1.75e+25) {
tmp = b * t;
} else if (t <= -1.08e-305) {
tmp = t_1;
} else if (t <= 7.2e-156) {
tmp = x - -a;
} else if (t <= 620000000000.0) {
tmp = t_1;
} else {
tmp = -a * t;
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = (y - 2.0) * b tmp = 0 if t <= -1.75e+25: tmp = b * t elif t <= -1.08e-305: tmp = t_1 elif t <= 7.2e-156: tmp = x - -a elif t <= 620000000000.0: tmp = t_1 else: tmp = -a * t return tmp
function code(x, y, z, t, a, b) t_1 = Float64(Float64(y - 2.0) * b) tmp = 0.0 if (t <= -1.75e+25) tmp = Float64(b * t); elseif (t <= -1.08e-305) tmp = t_1; elseif (t <= 7.2e-156) tmp = Float64(x - Float64(-a)); elseif (t <= 620000000000.0) tmp = t_1; else tmp = Float64(Float64(-a) * t); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) t_1 = (y - 2.0) * b; tmp = 0.0; if (t <= -1.75e+25) tmp = b * t; elseif (t <= -1.08e-305) tmp = t_1; elseif (t <= 7.2e-156) tmp = x - -a; elseif (t <= 620000000000.0) tmp = t_1; else tmp = -a * t; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(N[(y - 2.0), $MachinePrecision] * b), $MachinePrecision]}, If[LessEqual[t, -1.75e+25], N[(b * t), $MachinePrecision], If[LessEqual[t, -1.08e-305], t$95$1, If[LessEqual[t, 7.2e-156], N[(x - (-a)), $MachinePrecision], If[LessEqual[t, 620000000000.0], t$95$1, N[((-a) * t), $MachinePrecision]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \left(y - 2\right) \cdot b\\
\mathbf{if}\;t \leq -1.75 \cdot 10^{+25}:\\
\;\;\;\;b \cdot t\\
\mathbf{elif}\;t \leq -1.08 \cdot 10^{-305}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;t \leq 7.2 \cdot 10^{-156}:\\
\;\;\;\;x - \left(-a\right)\\
\mathbf{elif}\;t \leq 620000000000:\\
\;\;\;\;t\_1\\
\mathbf{else}:\\
\;\;\;\;\left(-a\right) \cdot t\\
\end{array}
\end{array}
if t < -1.75e25Initial program 86.2%
Taylor expanded in b around inf
*-commutativeN/A
lower-*.f64N/A
lower--.f64N/A
+-commutativeN/A
lower-+.f6453.6
Applied rewrites53.6%
Taylor expanded in t around inf
Applied rewrites49.8%
if -1.75e25 < t < -1.08000000000000004e-305 or 7.19999999999999998e-156 < t < 6.2e11Initial program 99.0%
Taylor expanded in b around inf
*-commutativeN/A
lower-*.f64N/A
lower--.f64N/A
+-commutativeN/A
lower-+.f6444.9
Applied rewrites44.9%
Taylor expanded in t around 0
Applied rewrites43.6%
if -1.08000000000000004e-305 < t < 7.19999999999999998e-156Initial program 96.6%
lift-+.f64N/A
lift--.f64N/A
associate-+l-N/A
lift--.f64N/A
associate--l-N/A
lower--.f64N/A
lift-*.f64N/A
*-commutativeN/A
lower-fma.f64N/A
lower--.f6496.7
lift-*.f64N/A
*-commutativeN/A
lower-*.f6496.7
lift-*.f64N/A
*-commutativeN/A
lower-*.f6496.7
lift-+.f64N/A
+-commutativeN/A
lower-+.f6496.7
Applied rewrites96.7%
Taylor expanded in a around inf
*-commutativeN/A
lower-*.f64N/A
lower--.f6446.1
Applied rewrites46.1%
Taylor expanded in t around 0
Applied rewrites46.1%
if 6.2e11 < t Initial program 95.9%
Taylor expanded in t around inf
*-commutativeN/A
lower-*.f64N/A
lower--.f6465.8
Applied rewrites65.8%
Taylor expanded in b around 0
Applied rewrites39.7%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (fma (- (+ t y) 2.0) b x)))
(if (<= b -680000000.0)
t_1
(if (<= b 2.2e-7)
(- x (* z (- y 1.0)))
(if (<= b 6e+89) (+ (* (- b a) t) a) t_1)))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = fma(((t + y) - 2.0), b, x);
double tmp;
if (b <= -680000000.0) {
tmp = t_1;
} else if (b <= 2.2e-7) {
tmp = x - (z * (y - 1.0));
} else if (b <= 6e+89) {
tmp = ((b - a) * t) + a;
} else {
tmp = t_1;
}
return tmp;
}
function code(x, y, z, t, a, b) t_1 = fma(Float64(Float64(t + y) - 2.0), b, x) tmp = 0.0 if (b <= -680000000.0) tmp = t_1; elseif (b <= 2.2e-7) tmp = Float64(x - Float64(z * Float64(y - 1.0))); elseif (b <= 6e+89) tmp = Float64(Float64(Float64(b - a) * t) + a); else tmp = t_1; end return tmp end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(N[(N[(t + y), $MachinePrecision] - 2.0), $MachinePrecision] * b + x), $MachinePrecision]}, If[LessEqual[b, -680000000.0], t$95$1, If[LessEqual[b, 2.2e-7], N[(x - N[(z * N[(y - 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[b, 6e+89], N[(N[(N[(b - a), $MachinePrecision] * t), $MachinePrecision] + a), $MachinePrecision], t$95$1]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \mathsf{fma}\left(\left(t + y\right) - 2, b, x\right)\\
\mathbf{if}\;b \leq -680000000:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;b \leq 2.2 \cdot 10^{-7}:\\
\;\;\;\;x - z \cdot \left(y - 1\right)\\
\mathbf{elif}\;b \leq 6 \cdot 10^{+89}:\\
\;\;\;\;\left(b - a\right) \cdot t + a\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if b < -6.8e8 or 6.00000000000000025e89 < b Initial program 90.5%
Taylor expanded in z around 0
sub-negN/A
+-commutativeN/A
*-commutativeN/A
distribute-lft-neg-inN/A
mul-1-negN/A
lower-fma.f64N/A
mul-1-negN/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
distribute-neg-inN/A
metadata-evalN/A
sub-negN/A
lower--.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower--.f64N/A
+-commutativeN/A
lower-+.f6491.0
Applied rewrites91.0%
Taylor expanded in t around 0
Applied rewrites90.0%
Taylor expanded in a around 0
Applied rewrites85.5%
if -6.8e8 < b < 2.2000000000000001e-7Initial program 99.2%
lift-+.f64N/A
lift--.f64N/A
associate-+l-N/A
lift--.f64N/A
associate--l-N/A
lower--.f64N/A
lift-*.f64N/A
*-commutativeN/A
lower-fma.f64N/A
lower--.f6499.2
lift-*.f64N/A
*-commutativeN/A
lower-*.f6499.2
lift-*.f64N/A
*-commutativeN/A
lower-*.f6499.2
lift-+.f64N/A
+-commutativeN/A
lower-+.f6499.2
Applied rewrites99.2%
Taylor expanded in z around inf
*-commutativeN/A
lower-*.f64N/A
lower--.f6462.8
Applied rewrites62.8%
if 2.2000000000000001e-7 < b < 6.00000000000000025e89Initial program 96.4%
Taylor expanded in z around 0
sub-negN/A
+-commutativeN/A
*-commutativeN/A
distribute-lft-neg-inN/A
mul-1-negN/A
lower-fma.f64N/A
mul-1-negN/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
distribute-neg-inN/A
metadata-evalN/A
sub-negN/A
lower--.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower--.f64N/A
+-commutativeN/A
lower-+.f6483.0
Applied rewrites83.0%
Taylor expanded in t around 0
Applied rewrites86.6%
Taylor expanded in t around inf
Applied rewrites61.9%
Final simplification72.1%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (- x (* (- z b) y))))
(if (<= y -5.4e+62)
t_1
(if (<= y 1.12e+133) (fma (- 1.0 t) a (fma (- t 2.0) b x)) t_1))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = x - ((z - b) * y);
double tmp;
if (y <= -5.4e+62) {
tmp = t_1;
} else if (y <= 1.12e+133) {
tmp = fma((1.0 - t), a, fma((t - 2.0), b, x));
} else {
tmp = t_1;
}
return tmp;
}
function code(x, y, z, t, a, b) t_1 = Float64(x - Float64(Float64(z - b) * y)) tmp = 0.0 if (y <= -5.4e+62) tmp = t_1; elseif (y <= 1.12e+133) tmp = fma(Float64(1.0 - t), a, fma(Float64(t - 2.0), b, x)); else tmp = t_1; end return tmp end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(x - N[(N[(z - b), $MachinePrecision] * y), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, -5.4e+62], t$95$1, If[LessEqual[y, 1.12e+133], N[(N[(1.0 - t), $MachinePrecision] * a + N[(N[(t - 2.0), $MachinePrecision] * b + x), $MachinePrecision]), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x - \left(z - b\right) \cdot y\\
\mathbf{if}\;y \leq -5.4 \cdot 10^{+62}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y \leq 1.12 \cdot 10^{+133}:\\
\;\;\;\;\mathsf{fma}\left(1 - t, a, \mathsf{fma}\left(t - 2, b, x\right)\right)\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if y < -5.4e62 or 1.12e133 < y Initial program 90.9%
lift-+.f64N/A
lift--.f64N/A
associate-+l-N/A
lift--.f64N/A
associate--l-N/A
lower--.f64N/A
lift-*.f64N/A
*-commutativeN/A
lower-fma.f64N/A
lower--.f6492.0
lift-*.f64N/A
*-commutativeN/A
lower-*.f6492.0
lift-*.f64N/A
*-commutativeN/A
lower-*.f6492.0
lift-+.f64N/A
+-commutativeN/A
lower-+.f6492.0
Applied rewrites92.0%
Taylor expanded in y around inf
*-commutativeN/A
sub-negN/A
mul-1-negN/A
+-commutativeN/A
*-lft-identityN/A
metadata-evalN/A
cancel-sign-sub-invN/A
lower-*.f64N/A
cancel-sign-sub-invN/A
metadata-evalN/A
*-lft-identityN/A
+-commutativeN/A
mul-1-negN/A
sub-negN/A
lower--.f6485.3
Applied rewrites85.3%
if -5.4e62 < y < 1.12e133Initial program 97.6%
Taylor expanded in z around 0
sub-negN/A
+-commutativeN/A
*-commutativeN/A
distribute-lft-neg-inN/A
mul-1-negN/A
lower-fma.f64N/A
mul-1-negN/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
distribute-neg-inN/A
metadata-evalN/A
sub-negN/A
lower--.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower--.f64N/A
+-commutativeN/A
lower-+.f6476.7
Applied rewrites76.7%
Taylor expanded in y around 0
Applied rewrites73.8%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (* (- (+ t y) 2.0) b)))
(if (<= b -750000000.0)
t_1
(if (<= b 2.2e-7)
(- x (* z (- y 1.0)))
(if (<= b 7.4e+90) (+ (* (- b a) t) a) t_1)))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = ((t + y) - 2.0) * b;
double tmp;
if (b <= -750000000.0) {
tmp = t_1;
} else if (b <= 2.2e-7) {
tmp = x - (z * (y - 1.0));
} else if (b <= 7.4e+90) {
tmp = ((b - a) * t) + a;
} else {
tmp = t_1;
}
return tmp;
}
real(8) function code(x, y, z, t, a, b)
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) :: t_1
real(8) :: tmp
t_1 = ((t + y) - 2.0d0) * b
if (b <= (-750000000.0d0)) then
tmp = t_1
else if (b <= 2.2d-7) then
tmp = x - (z * (y - 1.0d0))
else if (b <= 7.4d+90) then
tmp = ((b - a) * t) + a
else
tmp = t_1
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b) {
double t_1 = ((t + y) - 2.0) * b;
double tmp;
if (b <= -750000000.0) {
tmp = t_1;
} else if (b <= 2.2e-7) {
tmp = x - (z * (y - 1.0));
} else if (b <= 7.4e+90) {
tmp = ((b - a) * t) + a;
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = ((t + y) - 2.0) * b tmp = 0 if b <= -750000000.0: tmp = t_1 elif b <= 2.2e-7: tmp = x - (z * (y - 1.0)) elif b <= 7.4e+90: tmp = ((b - a) * t) + a else: tmp = t_1 return tmp
function code(x, y, z, t, a, b) t_1 = Float64(Float64(Float64(t + y) - 2.0) * b) tmp = 0.0 if (b <= -750000000.0) tmp = t_1; elseif (b <= 2.2e-7) tmp = Float64(x - Float64(z * Float64(y - 1.0))); elseif (b <= 7.4e+90) tmp = Float64(Float64(Float64(b - a) * t) + a); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a, b) t_1 = ((t + y) - 2.0) * b; tmp = 0.0; if (b <= -750000000.0) tmp = t_1; elseif (b <= 2.2e-7) tmp = x - (z * (y - 1.0)); elseif (b <= 7.4e+90) tmp = ((b - a) * t) + a; else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(N[(N[(t + y), $MachinePrecision] - 2.0), $MachinePrecision] * b), $MachinePrecision]}, If[LessEqual[b, -750000000.0], t$95$1, If[LessEqual[b, 2.2e-7], N[(x - N[(z * N[(y - 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[b, 7.4e+90], N[(N[(N[(b - a), $MachinePrecision] * t), $MachinePrecision] + a), $MachinePrecision], t$95$1]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \left(\left(t + y\right) - 2\right) \cdot b\\
\mathbf{if}\;b \leq -750000000:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;b \leq 2.2 \cdot 10^{-7}:\\
\;\;\;\;x - z \cdot \left(y - 1\right)\\
\mathbf{elif}\;b \leq 7.4 \cdot 10^{+90}:\\
\;\;\;\;\left(b - a\right) \cdot t + a\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if b < -7.5e8 or 7.4e90 < b Initial program 90.5%
Taylor expanded in b around inf
*-commutativeN/A
lower-*.f64N/A
lower--.f64N/A
+-commutativeN/A
lower-+.f6478.6
Applied rewrites78.6%
if -7.5e8 < b < 2.2000000000000001e-7Initial program 99.2%
lift-+.f64N/A
lift--.f64N/A
associate-+l-N/A
lift--.f64N/A
associate--l-N/A
lower--.f64N/A
lift-*.f64N/A
*-commutativeN/A
lower-fma.f64N/A
lower--.f6499.2
lift-*.f64N/A
*-commutativeN/A
lower-*.f6499.2
lift-*.f64N/A
*-commutativeN/A
lower-*.f6499.2
lift-+.f64N/A
+-commutativeN/A
lower-+.f6499.2
Applied rewrites99.2%
Taylor expanded in z around inf
*-commutativeN/A
lower-*.f64N/A
lower--.f6462.8
Applied rewrites62.8%
if 2.2000000000000001e-7 < b < 7.4e90Initial program 96.4%
Taylor expanded in z around 0
sub-negN/A
+-commutativeN/A
*-commutativeN/A
distribute-lft-neg-inN/A
mul-1-negN/A
lower-fma.f64N/A
mul-1-negN/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
distribute-neg-inN/A
metadata-evalN/A
sub-negN/A
lower--.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower--.f64N/A
+-commutativeN/A
lower-+.f6483.0
Applied rewrites83.0%
Taylor expanded in t around 0
Applied rewrites86.6%
Taylor expanded in t around inf
Applied rewrites61.9%
Final simplification69.2%
(FPCore (x y z t a b)
:precision binary64
(if (<= b -2.4e+15)
(* (- y 2.0) b)
(if (<= b 1.65e-7)
(* (- 1.0 y) z)
(if (<= b 2.6e+60) (+ (* b y) a) (* (- t 2.0) b)))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (b <= -2.4e+15) {
tmp = (y - 2.0) * b;
} else if (b <= 1.65e-7) {
tmp = (1.0 - y) * z;
} else if (b <= 2.6e+60) {
tmp = (b * y) + a;
} else {
tmp = (t - 2.0) * b;
}
return tmp;
}
real(8) function code(x, y, z, t, a, b)
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) :: tmp
if (b <= (-2.4d+15)) then
tmp = (y - 2.0d0) * b
else if (b <= 1.65d-7) then
tmp = (1.0d0 - y) * z
else if (b <= 2.6d+60) then
tmp = (b * y) + a
else
tmp = (t - 2.0d0) * b
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (b <= -2.4e+15) {
tmp = (y - 2.0) * b;
} else if (b <= 1.65e-7) {
tmp = (1.0 - y) * z;
} else if (b <= 2.6e+60) {
tmp = (b * y) + a;
} else {
tmp = (t - 2.0) * b;
}
return tmp;
}
def code(x, y, z, t, a, b): tmp = 0 if b <= -2.4e+15: tmp = (y - 2.0) * b elif b <= 1.65e-7: tmp = (1.0 - y) * z elif b <= 2.6e+60: tmp = (b * y) + a else: tmp = (t - 2.0) * b return tmp
function code(x, y, z, t, a, b) tmp = 0.0 if (b <= -2.4e+15) tmp = Float64(Float64(y - 2.0) * b); elseif (b <= 1.65e-7) tmp = Float64(Float64(1.0 - y) * z); elseif (b <= 2.6e+60) tmp = Float64(Float64(b * y) + a); else tmp = Float64(Float64(t - 2.0) * b); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) tmp = 0.0; if (b <= -2.4e+15) tmp = (y - 2.0) * b; elseif (b <= 1.65e-7) tmp = (1.0 - y) * z; elseif (b <= 2.6e+60) tmp = (b * y) + a; else tmp = (t - 2.0) * b; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[b, -2.4e+15], N[(N[(y - 2.0), $MachinePrecision] * b), $MachinePrecision], If[LessEqual[b, 1.65e-7], N[(N[(1.0 - y), $MachinePrecision] * z), $MachinePrecision], If[LessEqual[b, 2.6e+60], N[(N[(b * y), $MachinePrecision] + a), $MachinePrecision], N[(N[(t - 2.0), $MachinePrecision] * b), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;b \leq -2.4 \cdot 10^{+15}:\\
\;\;\;\;\left(y - 2\right) \cdot b\\
\mathbf{elif}\;b \leq 1.65 \cdot 10^{-7}:\\
\;\;\;\;\left(1 - y\right) \cdot z\\
\mathbf{elif}\;b \leq 2.6 \cdot 10^{+60}:\\
\;\;\;\;b \cdot y + a\\
\mathbf{else}:\\
\;\;\;\;\left(t - 2\right) \cdot b\\
\end{array}
\end{array}
if b < -2.4e15Initial program 93.2%
Taylor expanded in b around inf
*-commutativeN/A
lower-*.f64N/A
lower--.f64N/A
+-commutativeN/A
lower-+.f6480.0
Applied rewrites80.0%
Taylor expanded in t around 0
Applied rewrites54.5%
if -2.4e15 < b < 1.6500000000000001e-7Initial program 99.2%
Taylor expanded in z around inf
*-commutativeN/A
sub-negN/A
metadata-evalN/A
distribute-neg-inN/A
+-commutativeN/A
metadata-evalN/A
sub-negN/A
mul-1-negN/A
lower-*.f64N/A
sub-negN/A
metadata-evalN/A
distribute-lft-inN/A
metadata-evalN/A
+-commutativeN/A
neg-mul-1N/A
sub-negN/A
lower--.f6445.8
Applied rewrites45.8%
if 1.6500000000000001e-7 < b < 2.60000000000000008e60Initial program 99.9%
Taylor expanded in z around 0
sub-negN/A
+-commutativeN/A
*-commutativeN/A
distribute-lft-neg-inN/A
mul-1-negN/A
lower-fma.f64N/A
mul-1-negN/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
distribute-neg-inN/A
metadata-evalN/A
sub-negN/A
lower--.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower--.f64N/A
+-commutativeN/A
lower-+.f6485.9
Applied rewrites85.9%
Taylor expanded in t around 0
Applied rewrites86.0%
Taylor expanded in y around inf
Applied rewrites52.0%
if 2.60000000000000008e60 < b Initial program 87.2%
Taylor expanded in b around inf
*-commutativeN/A
lower-*.f64N/A
lower--.f64N/A
+-commutativeN/A
lower-+.f6471.2
Applied rewrites71.2%
Taylor expanded in y around 0
Applied rewrites54.0%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (* (- 1.0 t) a)))
(if (<= a -8e+75)
t_1
(if (<= a 4.5e-87)
(* (- t 2.0) b)
(if (<= a 4.8e+186) (* (- y 2.0) b) t_1)))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = (1.0 - t) * a;
double tmp;
if (a <= -8e+75) {
tmp = t_1;
} else if (a <= 4.5e-87) {
tmp = (t - 2.0) * b;
} else if (a <= 4.8e+186) {
tmp = (y - 2.0) * b;
} else {
tmp = t_1;
}
return tmp;
}
real(8) function code(x, y, z, t, a, b)
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) :: t_1
real(8) :: tmp
t_1 = (1.0d0 - t) * a
if (a <= (-8d+75)) then
tmp = t_1
else if (a <= 4.5d-87) then
tmp = (t - 2.0d0) * b
else if (a <= 4.8d+186) then
tmp = (y - 2.0d0) * b
else
tmp = t_1
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b) {
double t_1 = (1.0 - t) * a;
double tmp;
if (a <= -8e+75) {
tmp = t_1;
} else if (a <= 4.5e-87) {
tmp = (t - 2.0) * b;
} else if (a <= 4.8e+186) {
tmp = (y - 2.0) * b;
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = (1.0 - t) * a tmp = 0 if a <= -8e+75: tmp = t_1 elif a <= 4.5e-87: tmp = (t - 2.0) * b elif a <= 4.8e+186: tmp = (y - 2.0) * b else: tmp = t_1 return tmp
function code(x, y, z, t, a, b) t_1 = Float64(Float64(1.0 - t) * a) tmp = 0.0 if (a <= -8e+75) tmp = t_1; elseif (a <= 4.5e-87) tmp = Float64(Float64(t - 2.0) * b); elseif (a <= 4.8e+186) tmp = Float64(Float64(y - 2.0) * b); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a, b) t_1 = (1.0 - t) * a; tmp = 0.0; if (a <= -8e+75) tmp = t_1; elseif (a <= 4.5e-87) tmp = (t - 2.0) * b; elseif (a <= 4.8e+186) tmp = (y - 2.0) * b; else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(N[(1.0 - t), $MachinePrecision] * a), $MachinePrecision]}, If[LessEqual[a, -8e+75], t$95$1, If[LessEqual[a, 4.5e-87], N[(N[(t - 2.0), $MachinePrecision] * b), $MachinePrecision], If[LessEqual[a, 4.8e+186], N[(N[(y - 2.0), $MachinePrecision] * b), $MachinePrecision], t$95$1]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \left(1 - t\right) \cdot a\\
\mathbf{if}\;a \leq -8 \cdot 10^{+75}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;a \leq 4.5 \cdot 10^{-87}:\\
\;\;\;\;\left(t - 2\right) \cdot b\\
\mathbf{elif}\;a \leq 4.8 \cdot 10^{+186}:\\
\;\;\;\;\left(y - 2\right) \cdot b\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if a < -7.99999999999999941e75 or 4.7999999999999999e186 < a Initial program 91.2%
Taylor expanded in a around inf
*-commutativeN/A
sub-negN/A
metadata-evalN/A
distribute-neg-inN/A
+-commutativeN/A
metadata-evalN/A
sub-negN/A
mul-1-negN/A
lower-*.f64N/A
mul-1-negN/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
distribute-neg-inN/A
metadata-evalN/A
sub-negN/A
lower--.f6462.8
Applied rewrites62.8%
if -7.99999999999999941e75 < a < 4.49999999999999958e-87Initial program 98.4%
Taylor expanded in b around inf
*-commutativeN/A
lower-*.f64N/A
lower--.f64N/A
+-commutativeN/A
lower-+.f6448.4
Applied rewrites48.4%
Taylor expanded in y around 0
Applied rewrites36.1%
if 4.49999999999999958e-87 < a < 4.7999999999999999e186Initial program 94.1%
Taylor expanded in b around inf
*-commutativeN/A
lower-*.f64N/A
lower--.f64N/A
+-commutativeN/A
lower-+.f6449.7
Applied rewrites49.7%
Taylor expanded in t around 0
Applied rewrites38.3%
(FPCore (x y z t a b)
:precision binary64
(if (<= t -1.75e+25)
(* b t)
(if (<= t 2.7e-238)
(+ (* b y) a)
(if (<= t 620000000000.0) (- x (- a)) (* (- a) t)))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (t <= -1.75e+25) {
tmp = b * t;
} else if (t <= 2.7e-238) {
tmp = (b * y) + a;
} else if (t <= 620000000000.0) {
tmp = x - -a;
} else {
tmp = -a * t;
}
return tmp;
}
real(8) function code(x, y, z, t, a, b)
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) :: tmp
if (t <= (-1.75d+25)) then
tmp = b * t
else if (t <= 2.7d-238) then
tmp = (b * y) + a
else if (t <= 620000000000.0d0) then
tmp = x - -a
else
tmp = -a * t
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (t <= -1.75e+25) {
tmp = b * t;
} else if (t <= 2.7e-238) {
tmp = (b * y) + a;
} else if (t <= 620000000000.0) {
tmp = x - -a;
} else {
tmp = -a * t;
}
return tmp;
}
def code(x, y, z, t, a, b): tmp = 0 if t <= -1.75e+25: tmp = b * t elif t <= 2.7e-238: tmp = (b * y) + a elif t <= 620000000000.0: tmp = x - -a else: tmp = -a * t return tmp
function code(x, y, z, t, a, b) tmp = 0.0 if (t <= -1.75e+25) tmp = Float64(b * t); elseif (t <= 2.7e-238) tmp = Float64(Float64(b * y) + a); elseif (t <= 620000000000.0) tmp = Float64(x - Float64(-a)); else tmp = Float64(Float64(-a) * t); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) tmp = 0.0; if (t <= -1.75e+25) tmp = b * t; elseif (t <= 2.7e-238) tmp = (b * y) + a; elseif (t <= 620000000000.0) tmp = x - -a; else tmp = -a * t; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[t, -1.75e+25], N[(b * t), $MachinePrecision], If[LessEqual[t, 2.7e-238], N[(N[(b * y), $MachinePrecision] + a), $MachinePrecision], If[LessEqual[t, 620000000000.0], N[(x - (-a)), $MachinePrecision], N[((-a) * t), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;t \leq -1.75 \cdot 10^{+25}:\\
\;\;\;\;b \cdot t\\
\mathbf{elif}\;t \leq 2.7 \cdot 10^{-238}:\\
\;\;\;\;b \cdot y + a\\
\mathbf{elif}\;t \leq 620000000000:\\
\;\;\;\;x - \left(-a\right)\\
\mathbf{else}:\\
\;\;\;\;\left(-a\right) \cdot t\\
\end{array}
\end{array}
if t < -1.75e25Initial program 86.2%
Taylor expanded in b around inf
*-commutativeN/A
lower-*.f64N/A
lower--.f64N/A
+-commutativeN/A
lower-+.f6453.6
Applied rewrites53.6%
Taylor expanded in t around inf
Applied rewrites49.8%
if -1.75e25 < t < 2.69999999999999991e-238Initial program 98.7%
Taylor expanded in z around 0
sub-negN/A
+-commutativeN/A
*-commutativeN/A
distribute-lft-neg-inN/A
mul-1-negN/A
lower-fma.f64N/A
mul-1-negN/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
distribute-neg-inN/A
metadata-evalN/A
sub-negN/A
lower--.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower--.f64N/A
+-commutativeN/A
lower-+.f6468.7
Applied rewrites68.7%
Taylor expanded in t around 0
Applied rewrites68.6%
Taylor expanded in y around inf
Applied rewrites44.0%
if 2.69999999999999991e-238 < t < 6.2e11Initial program 98.1%
lift-+.f64N/A
lift--.f64N/A
associate-+l-N/A
lift--.f64N/A
associate--l-N/A
lower--.f64N/A
lift-*.f64N/A
*-commutativeN/A
lower-fma.f64N/A
lower--.f6498.1
lift-*.f64N/A
*-commutativeN/A
lower-*.f6498.1
lift-*.f64N/A
*-commutativeN/A
lower-*.f6498.1
lift-+.f64N/A
+-commutativeN/A
lower-+.f6498.1
Applied rewrites98.1%
Taylor expanded in a around inf
*-commutativeN/A
lower-*.f64N/A
lower--.f6433.4
Applied rewrites33.4%
Taylor expanded in t around 0
Applied rewrites32.2%
if 6.2e11 < t Initial program 95.9%
Taylor expanded in t around inf
*-commutativeN/A
lower-*.f64N/A
lower--.f6465.8
Applied rewrites65.8%
Taylor expanded in b around 0
Applied rewrites39.7%
(FPCore (x y z t a b)
:precision binary64
(if (<= t -1.75e+25)
(* b t)
(if (<= t -3.15e-285)
(* b y)
(if (<= t 620000000000.0) (- x (- a)) (* (- a) t)))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (t <= -1.75e+25) {
tmp = b * t;
} else if (t <= -3.15e-285) {
tmp = b * y;
} else if (t <= 620000000000.0) {
tmp = x - -a;
} else {
tmp = -a * t;
}
return tmp;
}
real(8) function code(x, y, z, t, a, b)
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) :: tmp
if (t <= (-1.75d+25)) then
tmp = b * t
else if (t <= (-3.15d-285)) then
tmp = b * y
else if (t <= 620000000000.0d0) then
tmp = x - -a
else
tmp = -a * t
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (t <= -1.75e+25) {
tmp = b * t;
} else if (t <= -3.15e-285) {
tmp = b * y;
} else if (t <= 620000000000.0) {
tmp = x - -a;
} else {
tmp = -a * t;
}
return tmp;
}
def code(x, y, z, t, a, b): tmp = 0 if t <= -1.75e+25: tmp = b * t elif t <= -3.15e-285: tmp = b * y elif t <= 620000000000.0: tmp = x - -a else: tmp = -a * t return tmp
function code(x, y, z, t, a, b) tmp = 0.0 if (t <= -1.75e+25) tmp = Float64(b * t); elseif (t <= -3.15e-285) tmp = Float64(b * y); elseif (t <= 620000000000.0) tmp = Float64(x - Float64(-a)); else tmp = Float64(Float64(-a) * t); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) tmp = 0.0; if (t <= -1.75e+25) tmp = b * t; elseif (t <= -3.15e-285) tmp = b * y; elseif (t <= 620000000000.0) tmp = x - -a; else tmp = -a * t; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[t, -1.75e+25], N[(b * t), $MachinePrecision], If[LessEqual[t, -3.15e-285], N[(b * y), $MachinePrecision], If[LessEqual[t, 620000000000.0], N[(x - (-a)), $MachinePrecision], N[((-a) * t), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;t \leq -1.75 \cdot 10^{+25}:\\
\;\;\;\;b \cdot t\\
\mathbf{elif}\;t \leq -3.15 \cdot 10^{-285}:\\
\;\;\;\;b \cdot y\\
\mathbf{elif}\;t \leq 620000000000:\\
\;\;\;\;x - \left(-a\right)\\
\mathbf{else}:\\
\;\;\;\;\left(-a\right) \cdot t\\
\end{array}
\end{array}
if t < -1.75e25Initial program 86.2%
Taylor expanded in b around inf
*-commutativeN/A
lower-*.f64N/A
lower--.f64N/A
+-commutativeN/A
lower-+.f6453.6
Applied rewrites53.6%
Taylor expanded in t around inf
Applied rewrites49.8%
if -1.75e25 < t < -3.14999999999999994e-285Initial program 98.3%
Taylor expanded in z around 0
sub-negN/A
+-commutativeN/A
*-commutativeN/A
distribute-lft-neg-inN/A
mul-1-negN/A
lower-fma.f64N/A
mul-1-negN/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
distribute-neg-inN/A
metadata-evalN/A
sub-negN/A
lower--.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower--.f64N/A
+-commutativeN/A
lower-+.f6472.4
Applied rewrites72.4%
Taylor expanded in y around inf
Applied rewrites35.5%
if -3.14999999999999994e-285 < t < 6.2e11Initial program 98.6%
lift-+.f64N/A
lift--.f64N/A
associate-+l-N/A
lift--.f64N/A
associate--l-N/A
lower--.f64N/A
lift-*.f64N/A
*-commutativeN/A
lower-fma.f64N/A
lower--.f6498.6
lift-*.f64N/A
*-commutativeN/A
lower-*.f6498.6
lift-*.f64N/A
*-commutativeN/A
lower-*.f6498.6
lift-+.f64N/A
+-commutativeN/A
lower-+.f6498.6
Applied rewrites98.6%
Taylor expanded in a around inf
*-commutativeN/A
lower-*.f64N/A
lower--.f6436.3
Applied rewrites36.3%
Taylor expanded in t around 0
Applied rewrites35.5%
if 6.2e11 < t Initial program 95.9%
Taylor expanded in t around inf
*-commutativeN/A
lower-*.f64N/A
lower--.f6465.8
Applied rewrites65.8%
Taylor expanded in b around 0
Applied rewrites39.7%
(FPCore (x y z t a b)
:precision binary64
(let* ((t_1 (- x (* (- a b) t))))
(if (<= t -75000000.0)
t_1
(if (<= t 500000000000.0) (fma (- y 2.0) b (+ a x)) t_1))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = x - ((a - b) * t);
double tmp;
if (t <= -75000000.0) {
tmp = t_1;
} else if (t <= 500000000000.0) {
tmp = fma((y - 2.0), b, (a + x));
} else {
tmp = t_1;
}
return tmp;
}
function code(x, y, z, t, a, b) t_1 = Float64(x - Float64(Float64(a - b) * t)) tmp = 0.0 if (t <= -75000000.0) tmp = t_1; elseif (t <= 500000000000.0) tmp = fma(Float64(y - 2.0), b, Float64(a + x)); else tmp = t_1; end return tmp end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(x - N[(N[(a - b), $MachinePrecision] * t), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t, -75000000.0], t$95$1, If[LessEqual[t, 500000000000.0], N[(N[(y - 2.0), $MachinePrecision] * b + N[(a + x), $MachinePrecision]), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x - \left(a - b\right) \cdot t\\
\mathbf{if}\;t \leq -75000000:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;t \leq 500000000000:\\
\;\;\;\;\mathsf{fma}\left(y - 2, b, a + x\right)\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if t < -7.5e7 or 5e11 < t Initial program 92.3%
lift-+.f64N/A
lift--.f64N/A
associate-+l-N/A
lift--.f64N/A
associate--l-N/A
lower--.f64N/A
lift-*.f64N/A
*-commutativeN/A
lower-fma.f64N/A
lower--.f6493.1
lift-*.f64N/A
*-commutativeN/A
lower-*.f6493.1
lift-*.f64N/A
*-commutativeN/A
lower-*.f6493.1
lift-+.f64N/A
+-commutativeN/A
lower-+.f6493.1
Applied rewrites93.1%
Taylor expanded in t around inf
*-commutativeN/A
sub-negN/A
mul-1-negN/A
lower-*.f64N/A
mul-1-negN/A
sub-negN/A
lower--.f6474.5
Applied rewrites74.5%
if -7.5e7 < t < 5e11Initial program 98.4%
Taylor expanded in z around 0
sub-negN/A
+-commutativeN/A
*-commutativeN/A
distribute-lft-neg-inN/A
mul-1-negN/A
lower-fma.f64N/A
mul-1-negN/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
distribute-neg-inN/A
metadata-evalN/A
sub-negN/A
lower--.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower--.f64N/A
+-commutativeN/A
lower-+.f6467.1
Applied rewrites67.1%
Taylor expanded in t around 0
Applied rewrites65.9%
(FPCore (x y z t a b)
:precision binary64
(if (<= t -1.75e+25)
(* b t)
(if (<= t -3.15e-285)
(* b y)
(if (<= t 17000000000000.0) (- x (- a)) (* b t)))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (t <= -1.75e+25) {
tmp = b * t;
} else if (t <= -3.15e-285) {
tmp = b * y;
} else if (t <= 17000000000000.0) {
tmp = x - -a;
} else {
tmp = b * t;
}
return tmp;
}
real(8) function code(x, y, z, t, a, b)
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) :: tmp
if (t <= (-1.75d+25)) then
tmp = b * t
else if (t <= (-3.15d-285)) then
tmp = b * y
else if (t <= 17000000000000.0d0) then
tmp = x - -a
else
tmp = b * t
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (t <= -1.75e+25) {
tmp = b * t;
} else if (t <= -3.15e-285) {
tmp = b * y;
} else if (t <= 17000000000000.0) {
tmp = x - -a;
} else {
tmp = b * t;
}
return tmp;
}
def code(x, y, z, t, a, b): tmp = 0 if t <= -1.75e+25: tmp = b * t elif t <= -3.15e-285: tmp = b * y elif t <= 17000000000000.0: tmp = x - -a else: tmp = b * t return tmp
function code(x, y, z, t, a, b) tmp = 0.0 if (t <= -1.75e+25) tmp = Float64(b * t); elseif (t <= -3.15e-285) tmp = Float64(b * y); elseif (t <= 17000000000000.0) tmp = Float64(x - Float64(-a)); else tmp = Float64(b * t); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) tmp = 0.0; if (t <= -1.75e+25) tmp = b * t; elseif (t <= -3.15e-285) tmp = b * y; elseif (t <= 17000000000000.0) tmp = x - -a; else tmp = b * t; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[t, -1.75e+25], N[(b * t), $MachinePrecision], If[LessEqual[t, -3.15e-285], N[(b * y), $MachinePrecision], If[LessEqual[t, 17000000000000.0], N[(x - (-a)), $MachinePrecision], N[(b * t), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;t \leq -1.75 \cdot 10^{+25}:\\
\;\;\;\;b \cdot t\\
\mathbf{elif}\;t \leq -3.15 \cdot 10^{-285}:\\
\;\;\;\;b \cdot y\\
\mathbf{elif}\;t \leq 17000000000000:\\
\;\;\;\;x - \left(-a\right)\\
\mathbf{else}:\\
\;\;\;\;b \cdot t\\
\end{array}
\end{array}
if t < -1.75e25 or 1.7e13 < t Initial program 92.0%
Taylor expanded in b around inf
*-commutativeN/A
lower-*.f64N/A
lower--.f64N/A
+-commutativeN/A
lower-+.f6445.7
Applied rewrites45.7%
Taylor expanded in t around inf
Applied rewrites41.6%
if -1.75e25 < t < -3.14999999999999994e-285Initial program 98.3%
Taylor expanded in z around 0
sub-negN/A
+-commutativeN/A
*-commutativeN/A
distribute-lft-neg-inN/A
mul-1-negN/A
lower-fma.f64N/A
mul-1-negN/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
distribute-neg-inN/A
metadata-evalN/A
sub-negN/A
lower--.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower--.f64N/A
+-commutativeN/A
lower-+.f6472.4
Applied rewrites72.4%
Taylor expanded in y around inf
Applied rewrites35.5%
if -3.14999999999999994e-285 < t < 1.7e13Initial program 98.6%
lift-+.f64N/A
lift--.f64N/A
associate-+l-N/A
lift--.f64N/A
associate--l-N/A
lower--.f64N/A
lift-*.f64N/A
*-commutativeN/A
lower-fma.f64N/A
lower--.f6498.6
lift-*.f64N/A
*-commutativeN/A
lower-*.f6498.6
lift-*.f64N/A
*-commutativeN/A
lower-*.f6498.6
lift-+.f64N/A
+-commutativeN/A
lower-+.f6498.6
Applied rewrites98.6%
Taylor expanded in a around inf
*-commutativeN/A
lower-*.f64N/A
lower--.f6436.3
Applied rewrites36.3%
Taylor expanded in t around 0
Applied rewrites35.5%
(FPCore (x y z t a b) :precision binary64 (if (<= t -1.75e+25) (* b t) (if (<= t 1.3e-154) (* b y) (if (<= t 490000000000.0) (* -2.0 b) (* b t)))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (t <= -1.75e+25) {
tmp = b * t;
} else if (t <= 1.3e-154) {
tmp = b * y;
} else if (t <= 490000000000.0) {
tmp = -2.0 * b;
} else {
tmp = b * t;
}
return tmp;
}
real(8) function code(x, y, z, t, a, b)
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) :: tmp
if (t <= (-1.75d+25)) then
tmp = b * t
else if (t <= 1.3d-154) then
tmp = b * y
else if (t <= 490000000000.0d0) then
tmp = (-2.0d0) * b
else
tmp = b * t
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (t <= -1.75e+25) {
tmp = b * t;
} else if (t <= 1.3e-154) {
tmp = b * y;
} else if (t <= 490000000000.0) {
tmp = -2.0 * b;
} else {
tmp = b * t;
}
return tmp;
}
def code(x, y, z, t, a, b): tmp = 0 if t <= -1.75e+25: tmp = b * t elif t <= 1.3e-154: tmp = b * y elif t <= 490000000000.0: tmp = -2.0 * b else: tmp = b * t return tmp
function code(x, y, z, t, a, b) tmp = 0.0 if (t <= -1.75e+25) tmp = Float64(b * t); elseif (t <= 1.3e-154) tmp = Float64(b * y); elseif (t <= 490000000000.0) tmp = Float64(-2.0 * b); else tmp = Float64(b * t); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) tmp = 0.0; if (t <= -1.75e+25) tmp = b * t; elseif (t <= 1.3e-154) tmp = b * y; elseif (t <= 490000000000.0) tmp = -2.0 * b; else tmp = b * t; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[t, -1.75e+25], N[(b * t), $MachinePrecision], If[LessEqual[t, 1.3e-154], N[(b * y), $MachinePrecision], If[LessEqual[t, 490000000000.0], N[(-2.0 * b), $MachinePrecision], N[(b * t), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;t \leq -1.75 \cdot 10^{+25}:\\
\;\;\;\;b \cdot t\\
\mathbf{elif}\;t \leq 1.3 \cdot 10^{-154}:\\
\;\;\;\;b \cdot y\\
\mathbf{elif}\;t \leq 490000000000:\\
\;\;\;\;-2 \cdot b\\
\mathbf{else}:\\
\;\;\;\;b \cdot t\\
\end{array}
\end{array}
if t < -1.75e25 or 4.9e11 < t Initial program 92.0%
Taylor expanded in b around inf
*-commutativeN/A
lower-*.f64N/A
lower--.f64N/A
+-commutativeN/A
lower-+.f6445.7
Applied rewrites45.7%
Taylor expanded in t around inf
Applied rewrites41.6%
if -1.75e25 < t < 1.3e-154Initial program 97.8%
Taylor expanded in z around 0
sub-negN/A
+-commutativeN/A
*-commutativeN/A
distribute-lft-neg-inN/A
mul-1-negN/A
lower-fma.f64N/A
mul-1-negN/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
distribute-neg-inN/A
metadata-evalN/A
sub-negN/A
lower--.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower--.f64N/A
+-commutativeN/A
lower-+.f6466.9
Applied rewrites66.9%
Taylor expanded in y around inf
Applied rewrites27.2%
if 1.3e-154 < t < 4.9e11Initial program 100.0%
Taylor expanded in b around inf
*-commutativeN/A
lower-*.f64N/A
lower--.f64N/A
+-commutativeN/A
lower-+.f6438.9
Applied rewrites38.9%
Taylor expanded in t around 0
Applied rewrites38.9%
Taylor expanded in y around 0
Applied rewrites29.4%
(FPCore (x y z t a b) :precision binary64 (if (<= t -1.75e+25) (* b t) (if (<= t -5.9e-288) (* b y) (if (<= t 300000.0) (* 1.0 a) (* b t)))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (t <= -1.75e+25) {
tmp = b * t;
} else if (t <= -5.9e-288) {
tmp = b * y;
} else if (t <= 300000.0) {
tmp = 1.0 * a;
} else {
tmp = b * t;
}
return tmp;
}
real(8) function code(x, y, z, t, a, b)
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) :: tmp
if (t <= (-1.75d+25)) then
tmp = b * t
else if (t <= (-5.9d-288)) then
tmp = b * y
else if (t <= 300000.0d0) then
tmp = 1.0d0 * a
else
tmp = b * t
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (t <= -1.75e+25) {
tmp = b * t;
} else if (t <= -5.9e-288) {
tmp = b * y;
} else if (t <= 300000.0) {
tmp = 1.0 * a;
} else {
tmp = b * t;
}
return tmp;
}
def code(x, y, z, t, a, b): tmp = 0 if t <= -1.75e+25: tmp = b * t elif t <= -5.9e-288: tmp = b * y elif t <= 300000.0: tmp = 1.0 * a else: tmp = b * t return tmp
function code(x, y, z, t, a, b) tmp = 0.0 if (t <= -1.75e+25) tmp = Float64(b * t); elseif (t <= -5.9e-288) tmp = Float64(b * y); elseif (t <= 300000.0) tmp = Float64(1.0 * a); else tmp = Float64(b * t); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) tmp = 0.0; if (t <= -1.75e+25) tmp = b * t; elseif (t <= -5.9e-288) tmp = b * y; elseif (t <= 300000.0) tmp = 1.0 * a; else tmp = b * t; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[t, -1.75e+25], N[(b * t), $MachinePrecision], If[LessEqual[t, -5.9e-288], N[(b * y), $MachinePrecision], If[LessEqual[t, 300000.0], N[(1.0 * a), $MachinePrecision], N[(b * t), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;t \leq -1.75 \cdot 10^{+25}:\\
\;\;\;\;b \cdot t\\
\mathbf{elif}\;t \leq -5.9 \cdot 10^{-288}:\\
\;\;\;\;b \cdot y\\
\mathbf{elif}\;t \leq 300000:\\
\;\;\;\;1 \cdot a\\
\mathbf{else}:\\
\;\;\;\;b \cdot t\\
\end{array}
\end{array}
if t < -1.75e25 or 3e5 < t Initial program 92.1%
Taylor expanded in b around inf
*-commutativeN/A
lower-*.f64N/A
lower--.f64N/A
+-commutativeN/A
lower-+.f6446.1
Applied rewrites46.1%
Taylor expanded in t around inf
Applied rewrites41.2%
if -1.75e25 < t < -5.90000000000000014e-288Initial program 98.3%
Taylor expanded in z around 0
sub-negN/A
+-commutativeN/A
*-commutativeN/A
distribute-lft-neg-inN/A
mul-1-negN/A
lower-fma.f64N/A
mul-1-negN/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
distribute-neg-inN/A
metadata-evalN/A
sub-negN/A
lower--.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower--.f64N/A
+-commutativeN/A
lower-+.f6472.4
Applied rewrites72.4%
Taylor expanded in y around inf
Applied rewrites35.5%
if -5.90000000000000014e-288 < t < 3e5Initial program 98.5%
Taylor expanded in a around inf
*-commutativeN/A
sub-negN/A
metadata-evalN/A
distribute-neg-inN/A
+-commutativeN/A
metadata-evalN/A
sub-negN/A
mul-1-negN/A
lower-*.f64N/A
mul-1-negN/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
distribute-neg-inN/A
metadata-evalN/A
sub-negN/A
lower--.f6420.0
Applied rewrites20.0%
Taylor expanded in t around 0
Applied rewrites19.4%
(FPCore (x y z t a b) :precision binary64 (let* ((t_1 (* (- b z) y))) (if (<= y -3.7e+60) t_1 (if (<= y 5.5e+115) (- x (* (- a b) t)) t_1))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = (b - z) * y;
double tmp;
if (y <= -3.7e+60) {
tmp = t_1;
} else if (y <= 5.5e+115) {
tmp = x - ((a - b) * t);
} else {
tmp = t_1;
}
return tmp;
}
real(8) function code(x, y, z, t, a, b)
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) :: t_1
real(8) :: tmp
t_1 = (b - z) * y
if (y <= (-3.7d+60)) then
tmp = t_1
else if (y <= 5.5d+115) then
tmp = x - ((a - b) * t)
else
tmp = t_1
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b) {
double t_1 = (b - z) * y;
double tmp;
if (y <= -3.7e+60) {
tmp = t_1;
} else if (y <= 5.5e+115) {
tmp = x - ((a - b) * t);
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = (b - z) * y tmp = 0 if y <= -3.7e+60: tmp = t_1 elif y <= 5.5e+115: tmp = x - ((a - b) * t) else: tmp = t_1 return tmp
function code(x, y, z, t, a, b) t_1 = Float64(Float64(b - z) * y) tmp = 0.0 if (y <= -3.7e+60) tmp = t_1; elseif (y <= 5.5e+115) tmp = Float64(x - Float64(Float64(a - b) * t)); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a, b) t_1 = (b - z) * y; tmp = 0.0; if (y <= -3.7e+60) tmp = t_1; elseif (y <= 5.5e+115) tmp = x - ((a - b) * t); else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(N[(b - z), $MachinePrecision] * y), $MachinePrecision]}, If[LessEqual[y, -3.7e+60], t$95$1, If[LessEqual[y, 5.5e+115], N[(x - N[(N[(a - b), $MachinePrecision] * t), $MachinePrecision]), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \left(b - z\right) \cdot y\\
\mathbf{if}\;y \leq -3.7 \cdot 10^{+60}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y \leq 5.5 \cdot 10^{+115}:\\
\;\;\;\;x - \left(a - b\right) \cdot t\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if y < -3.69999999999999988e60 or 5.5e115 < y Initial program 91.5%
Taylor expanded in y around inf
*-commutativeN/A
lower-*.f64N/A
lower--.f6476.0
Applied rewrites76.0%
if -3.69999999999999988e60 < y < 5.5e115Initial program 97.5%
lift-+.f64N/A
lift--.f64N/A
associate-+l-N/A
lift--.f64N/A
associate--l-N/A
lower--.f64N/A
lift-*.f64N/A
*-commutativeN/A
lower-fma.f64N/A
lower--.f6497.5
lift-*.f64N/A
*-commutativeN/A
lower-*.f6497.5
lift-*.f64N/A
*-commutativeN/A
lower-*.f6497.5
lift-+.f64N/A
+-commutativeN/A
lower-+.f6497.5
Applied rewrites97.5%
Taylor expanded in t around inf
*-commutativeN/A
sub-negN/A
mul-1-negN/A
lower-*.f64N/A
mul-1-negN/A
sub-negN/A
lower--.f6455.7
Applied rewrites55.7%
(FPCore (x y z t a b) :precision binary64 (let* ((t_1 (* (- b z) y))) (if (<= y -5.4e+62) t_1 (if (<= y 3.3e+114) (+ (* (- b a) t) a) t_1))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = (b - z) * y;
double tmp;
if (y <= -5.4e+62) {
tmp = t_1;
} else if (y <= 3.3e+114) {
tmp = ((b - a) * t) + a;
} else {
tmp = t_1;
}
return tmp;
}
real(8) function code(x, y, z, t, a, b)
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) :: t_1
real(8) :: tmp
t_1 = (b - z) * y
if (y <= (-5.4d+62)) then
tmp = t_1
else if (y <= 3.3d+114) then
tmp = ((b - a) * t) + a
else
tmp = t_1
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b) {
double t_1 = (b - z) * y;
double tmp;
if (y <= -5.4e+62) {
tmp = t_1;
} else if (y <= 3.3e+114) {
tmp = ((b - a) * t) + a;
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = (b - z) * y tmp = 0 if y <= -5.4e+62: tmp = t_1 elif y <= 3.3e+114: tmp = ((b - a) * t) + a else: tmp = t_1 return tmp
function code(x, y, z, t, a, b) t_1 = Float64(Float64(b - z) * y) tmp = 0.0 if (y <= -5.4e+62) tmp = t_1; elseif (y <= 3.3e+114) tmp = Float64(Float64(Float64(b - a) * t) + a); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a, b) t_1 = (b - z) * y; tmp = 0.0; if (y <= -5.4e+62) tmp = t_1; elseif (y <= 3.3e+114) tmp = ((b - a) * t) + a; else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(N[(b - z), $MachinePrecision] * y), $MachinePrecision]}, If[LessEqual[y, -5.4e+62], t$95$1, If[LessEqual[y, 3.3e+114], N[(N[(N[(b - a), $MachinePrecision] * t), $MachinePrecision] + a), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \left(b - z\right) \cdot y\\
\mathbf{if}\;y \leq -5.4 \cdot 10^{+62}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y \leq 3.3 \cdot 10^{+114}:\\
\;\;\;\;\left(b - a\right) \cdot t + a\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if y < -5.4e62 or 3.3000000000000001e114 < y Initial program 91.5%
Taylor expanded in y around inf
*-commutativeN/A
lower-*.f64N/A
lower--.f6476.8
Applied rewrites76.8%
if -5.4e62 < y < 3.3000000000000001e114Initial program 97.5%
Taylor expanded in z around 0
sub-negN/A
+-commutativeN/A
*-commutativeN/A
distribute-lft-neg-inN/A
mul-1-negN/A
lower-fma.f64N/A
mul-1-negN/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
distribute-neg-inN/A
metadata-evalN/A
sub-negN/A
lower--.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower--.f64N/A
+-commutativeN/A
lower-+.f6477.1
Applied rewrites77.1%
Taylor expanded in t around 0
Applied rewrites78.3%
Taylor expanded in t around inf
Applied rewrites50.3%
(FPCore (x y z t a b) :precision binary64 (let* ((t_1 (* (- b z) y))) (if (<= y -5.4e+62) t_1 (if (<= y 5.1e+35) (fma (- 1.0 t) a x) t_1))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = (b - z) * y;
double tmp;
if (y <= -5.4e+62) {
tmp = t_1;
} else if (y <= 5.1e+35) {
tmp = fma((1.0 - t), a, x);
} else {
tmp = t_1;
}
return tmp;
}
function code(x, y, z, t, a, b) t_1 = Float64(Float64(b - z) * y) tmp = 0.0 if (y <= -5.4e+62) tmp = t_1; elseif (y <= 5.1e+35) tmp = fma(Float64(1.0 - t), a, x); else tmp = t_1; end return tmp end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(N[(b - z), $MachinePrecision] * y), $MachinePrecision]}, If[LessEqual[y, -5.4e+62], t$95$1, If[LessEqual[y, 5.1e+35], N[(N[(1.0 - t), $MachinePrecision] * a + x), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \left(b - z\right) \cdot y\\
\mathbf{if}\;y \leq -5.4 \cdot 10^{+62}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y \leq 5.1 \cdot 10^{+35}:\\
\;\;\;\;\mathsf{fma}\left(1 - t, a, x\right)\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if y < -5.4e62 or 5.10000000000000017e35 < y Initial program 92.9%
Taylor expanded in y around inf
*-commutativeN/A
lower-*.f64N/A
lower--.f6471.2
Applied rewrites71.2%
if -5.4e62 < y < 5.10000000000000017e35Initial program 97.2%
Taylor expanded in z around 0
sub-negN/A
+-commutativeN/A
*-commutativeN/A
distribute-lft-neg-inN/A
mul-1-negN/A
lower-fma.f64N/A
mul-1-negN/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
distribute-neg-inN/A
metadata-evalN/A
sub-negN/A
lower--.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower--.f64N/A
+-commutativeN/A
lower-+.f6477.4
Applied rewrites77.4%
Taylor expanded in b around 0
Applied rewrites46.0%
(FPCore (x y z t a b) :precision binary64 (let* ((t_1 (+ (* b y) a))) (if (<= y -0.0034) t_1 (if (<= y 1.35e+129) (* (- t 2.0) b) t_1))))
double code(double x, double y, double z, double t, double a, double b) {
double t_1 = (b * y) + a;
double tmp;
if (y <= -0.0034) {
tmp = t_1;
} else if (y <= 1.35e+129) {
tmp = (t - 2.0) * b;
} else {
tmp = t_1;
}
return tmp;
}
real(8) function code(x, y, z, t, a, b)
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) :: t_1
real(8) :: tmp
t_1 = (b * y) + a
if (y <= (-0.0034d0)) then
tmp = t_1
else if (y <= 1.35d+129) then
tmp = (t - 2.0d0) * b
else
tmp = t_1
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b) {
double t_1 = (b * y) + a;
double tmp;
if (y <= -0.0034) {
tmp = t_1;
} else if (y <= 1.35e+129) {
tmp = (t - 2.0) * b;
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b): t_1 = (b * y) + a tmp = 0 if y <= -0.0034: tmp = t_1 elif y <= 1.35e+129: tmp = (t - 2.0) * b else: tmp = t_1 return tmp
function code(x, y, z, t, a, b) t_1 = Float64(Float64(b * y) + a) tmp = 0.0 if (y <= -0.0034) tmp = t_1; elseif (y <= 1.35e+129) tmp = Float64(Float64(t - 2.0) * b); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a, b) t_1 = (b * y) + a; tmp = 0.0; if (y <= -0.0034) tmp = t_1; elseif (y <= 1.35e+129) tmp = (t - 2.0) * b; else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := Block[{t$95$1 = N[(N[(b * y), $MachinePrecision] + a), $MachinePrecision]}, If[LessEqual[y, -0.0034], t$95$1, If[LessEqual[y, 1.35e+129], N[(N[(t - 2.0), $MachinePrecision] * b), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := b \cdot y + a\\
\mathbf{if}\;y \leq -0.0034:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;y \leq 1.35 \cdot 10^{+129}:\\
\;\;\;\;\left(t - 2\right) \cdot b\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if y < -0.00339999999999999981 or 1.35e129 < y Initial program 92.2%
Taylor expanded in z around 0
sub-negN/A
+-commutativeN/A
*-commutativeN/A
distribute-lft-neg-inN/A
mul-1-negN/A
lower-fma.f64N/A
mul-1-negN/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
distribute-neg-inN/A
metadata-evalN/A
sub-negN/A
lower--.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower--.f64N/A
+-commutativeN/A
lower-+.f6466.1
Applied rewrites66.1%
Taylor expanded in t around 0
Applied rewrites64.2%
Taylor expanded in y around inf
Applied rewrites45.1%
if -0.00339999999999999981 < y < 1.35e129Initial program 97.4%
Taylor expanded in b around inf
*-commutativeN/A
lower-*.f64N/A
lower--.f64N/A
+-commutativeN/A
lower-+.f6440.0
Applied rewrites40.0%
Taylor expanded in y around 0
Applied rewrites37.4%
(FPCore (x y z t a b) :precision binary64 (if (<= t -1.75e+25) (* b t) (if (<= t 17500000000000.0) (* b y) (* b t))))
double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (t <= -1.75e+25) {
tmp = b * t;
} else if (t <= 17500000000000.0) {
tmp = b * y;
} else {
tmp = b * t;
}
return tmp;
}
real(8) function code(x, y, z, t, a, b)
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) :: tmp
if (t <= (-1.75d+25)) then
tmp = b * t
else if (t <= 17500000000000.0d0) then
tmp = b * y
else
tmp = b * t
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b) {
double tmp;
if (t <= -1.75e+25) {
tmp = b * t;
} else if (t <= 17500000000000.0) {
tmp = b * y;
} else {
tmp = b * t;
}
return tmp;
}
def code(x, y, z, t, a, b): tmp = 0 if t <= -1.75e+25: tmp = b * t elif t <= 17500000000000.0: tmp = b * y else: tmp = b * t return tmp
function code(x, y, z, t, a, b) tmp = 0.0 if (t <= -1.75e+25) tmp = Float64(b * t); elseif (t <= 17500000000000.0) tmp = Float64(b * y); else tmp = Float64(b * t); end return tmp end
function tmp_2 = code(x, y, z, t, a, b) tmp = 0.0; if (t <= -1.75e+25) tmp = b * t; elseif (t <= 17500000000000.0) tmp = b * y; else tmp = b * t; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_] := If[LessEqual[t, -1.75e+25], N[(b * t), $MachinePrecision], If[LessEqual[t, 17500000000000.0], N[(b * y), $MachinePrecision], N[(b * t), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;t \leq -1.75 \cdot 10^{+25}:\\
\;\;\;\;b \cdot t\\
\mathbf{elif}\;t \leq 17500000000000:\\
\;\;\;\;b \cdot y\\
\mathbf{else}:\\
\;\;\;\;b \cdot t\\
\end{array}
\end{array}
if t < -1.75e25 or 1.75e13 < t Initial program 92.0%
Taylor expanded in b around inf
*-commutativeN/A
lower-*.f64N/A
lower--.f64N/A
+-commutativeN/A
lower-+.f6445.7
Applied rewrites45.7%
Taylor expanded in t around inf
Applied rewrites41.6%
if -1.75e25 < t < 1.75e13Initial program 98.4%
Taylor expanded in z around 0
sub-negN/A
+-commutativeN/A
*-commutativeN/A
distribute-lft-neg-inN/A
mul-1-negN/A
lower-fma.f64N/A
mul-1-negN/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
distribute-neg-inN/A
metadata-evalN/A
sub-negN/A
lower--.f64N/A
+-commutativeN/A
*-commutativeN/A
lower-fma.f64N/A
lower--.f64N/A
+-commutativeN/A
lower-+.f6467.9
Applied rewrites67.9%
Taylor expanded in y around inf
Applied rewrites23.0%
(FPCore (x y z t a b) :precision binary64 (* b t))
double code(double x, double y, double z, double t, double a, double b) {
return b * t;
}
real(8) function code(x, y, z, t, a, b)
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
code = b * t
end function
public static double code(double x, double y, double z, double t, double a, double b) {
return b * t;
}
def code(x, y, z, t, a, b): return b * t
function code(x, y, z, t, a, b) return Float64(b * t) end
function tmp = code(x, y, z, t, a, b) tmp = b * t; end
code[x_, y_, z_, t_, a_, b_] := N[(b * t), $MachinePrecision]
\begin{array}{l}
\\
b \cdot t
\end{array}
Initial program 95.3%
Taylor expanded in b around inf
*-commutativeN/A
lower-*.f64N/A
lower--.f64N/A
+-commutativeN/A
lower-+.f6441.8
Applied rewrites41.8%
Taylor expanded in t around inf
Applied rewrites21.9%
herbie shell --seed 2024255
(FPCore (x y z t a b)
:name "Statistics.Distribution.Beta:$centropy from math-functions-0.1.5.2"
:precision binary64
(+ (- (- x (* (- y 1.0) z)) (* (- t 1.0) a)) (* (- (+ y t) 2.0) b)))