
(FPCore (x y z t a b c) :precision binary64 (+ (- (+ (* x y) (/ (* z t) 16.0)) (/ (* a b) 4.0)) c))
double code(double x, double y, double z, double t, double a, double b, double c) {
return (((x * y) + ((z * t) / 16.0)) - ((a * b) / 4.0)) + c;
}
real(8) function code(x, y, z, t, a, b, c)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8), intent (in) :: c
code = (((x * y) + ((z * t) / 16.0d0)) - ((a * b) / 4.0d0)) + c
end function
public static double code(double x, double y, double z, double t, double a, double b, double c) {
return (((x * y) + ((z * t) / 16.0)) - ((a * b) / 4.0)) + c;
}
def code(x, y, z, t, a, b, c): return (((x * y) + ((z * t) / 16.0)) - ((a * b) / 4.0)) + c
function code(x, y, z, t, a, b, c) return Float64(Float64(Float64(Float64(x * y) + Float64(Float64(z * t) / 16.0)) - Float64(Float64(a * b) / 4.0)) + c) end
function tmp = code(x, y, z, t, a, b, c) tmp = (((x * y) + ((z * t) / 16.0)) - ((a * b) / 4.0)) + c; end
code[x_, y_, z_, t_, a_, b_, c_] := N[(N[(N[(N[(x * y), $MachinePrecision] + N[(N[(z * t), $MachinePrecision] / 16.0), $MachinePrecision]), $MachinePrecision] - N[(N[(a * b), $MachinePrecision] / 4.0), $MachinePrecision]), $MachinePrecision] + c), $MachinePrecision]
\begin{array}{l}
\\
\left(\left(x \cdot y + \frac{z \cdot t}{16}\right) - \frac{a \cdot b}{4}\right) + c
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 14 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (x y z t a b c) :precision binary64 (+ (- (+ (* x y) (/ (* z t) 16.0)) (/ (* a b) 4.0)) c))
double code(double x, double y, double z, double t, double a, double b, double c) {
return (((x * y) + ((z * t) / 16.0)) - ((a * b) / 4.0)) + c;
}
real(8) function code(x, y, z, t, a, b, c)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8), intent (in) :: c
code = (((x * y) + ((z * t) / 16.0d0)) - ((a * b) / 4.0d0)) + c
end function
public static double code(double x, double y, double z, double t, double a, double b, double c) {
return (((x * y) + ((z * t) / 16.0)) - ((a * b) / 4.0)) + c;
}
def code(x, y, z, t, a, b, c): return (((x * y) + ((z * t) / 16.0)) - ((a * b) / 4.0)) + c
function code(x, y, z, t, a, b, c) return Float64(Float64(Float64(Float64(x * y) + Float64(Float64(z * t) / 16.0)) - Float64(Float64(a * b) / 4.0)) + c) end
function tmp = code(x, y, z, t, a, b, c) tmp = (((x * y) + ((z * t) / 16.0)) - ((a * b) / 4.0)) + c; end
code[x_, y_, z_, t_, a_, b_, c_] := N[(N[(N[(N[(x * y), $MachinePrecision] + N[(N[(z * t), $MachinePrecision] / 16.0), $MachinePrecision]), $MachinePrecision] - N[(N[(a * b), $MachinePrecision] / 4.0), $MachinePrecision]), $MachinePrecision] + c), $MachinePrecision]
\begin{array}{l}
\\
\left(\left(x \cdot y + \frac{z \cdot t}{16}\right) - \frac{a \cdot b}{4}\right) + c
\end{array}
(FPCore (x y z t a b c)
:precision binary64
(let* ((t_1 (- (+ (* x y) (/ (* z t) 16.0)) (/ (* a b) 4.0))))
(if (<= t_1 1e+285)
(+ t_1 c)
(* b (- (/ (+ (* x y) (- c (* z (* t -0.0625)))) b) (* a 0.25))))))
double code(double x, double y, double z, double t, double a, double b, double c) {
double t_1 = ((x * y) + ((z * t) / 16.0)) - ((a * b) / 4.0);
double tmp;
if (t_1 <= 1e+285) {
tmp = t_1 + c;
} else {
tmp = b * ((((x * y) + (c - (z * (t * -0.0625)))) / b) - (a * 0.25));
}
return tmp;
}
real(8) function code(x, y, z, t, a, b, c)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8), intent (in) :: c
real(8) :: t_1
real(8) :: tmp
t_1 = ((x * y) + ((z * t) / 16.0d0)) - ((a * b) / 4.0d0)
if (t_1 <= 1d+285) then
tmp = t_1 + c
else
tmp = b * ((((x * y) + (c - (z * (t * (-0.0625d0))))) / b) - (a * 0.25d0))
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b, double c) {
double t_1 = ((x * y) + ((z * t) / 16.0)) - ((a * b) / 4.0);
double tmp;
if (t_1 <= 1e+285) {
tmp = t_1 + c;
} else {
tmp = b * ((((x * y) + (c - (z * (t * -0.0625)))) / b) - (a * 0.25));
}
return tmp;
}
def code(x, y, z, t, a, b, c): t_1 = ((x * y) + ((z * t) / 16.0)) - ((a * b) / 4.0) tmp = 0 if t_1 <= 1e+285: tmp = t_1 + c else: tmp = b * ((((x * y) + (c - (z * (t * -0.0625)))) / b) - (a * 0.25)) return tmp
function code(x, y, z, t, a, b, c) t_1 = Float64(Float64(Float64(x * y) + Float64(Float64(z * t) / 16.0)) - Float64(Float64(a * b) / 4.0)) tmp = 0.0 if (t_1 <= 1e+285) tmp = Float64(t_1 + c); else tmp = Float64(b * Float64(Float64(Float64(Float64(x * y) + Float64(c - Float64(z * Float64(t * -0.0625)))) / b) - Float64(a * 0.25))); end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c) t_1 = ((x * y) + ((z * t) / 16.0)) - ((a * b) / 4.0); tmp = 0.0; if (t_1 <= 1e+285) tmp = t_1 + c; else tmp = b * ((((x * y) + (c - (z * (t * -0.0625)))) / b) - (a * 0.25)); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_] := Block[{t$95$1 = N[(N[(N[(x * y), $MachinePrecision] + N[(N[(z * t), $MachinePrecision] / 16.0), $MachinePrecision]), $MachinePrecision] - N[(N[(a * b), $MachinePrecision] / 4.0), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$1, 1e+285], N[(t$95$1 + c), $MachinePrecision], N[(b * N[(N[(N[(N[(x * y), $MachinePrecision] + N[(c - N[(z * N[(t * -0.0625), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / b), $MachinePrecision] - N[(a * 0.25), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \left(x \cdot y + \frac{z \cdot t}{16}\right) - \frac{a \cdot b}{4}\\
\mathbf{if}\;t\_1 \leq 10^{+285}:\\
\;\;\;\;t\_1 + c\\
\mathbf{else}:\\
\;\;\;\;b \cdot \left(\frac{x \cdot y + \left(c - z \cdot \left(t \cdot -0.0625\right)\right)}{b} - a \cdot 0.25\right)\\
\end{array}
\end{array}
if (-.f64 (+.f64 (*.f64 x y) (/.f64 (*.f64 z t) #s(literal 16 binary64))) (/.f64 (*.f64 a b) #s(literal 4 binary64))) < 9.9999999999999998e284Initial program 100.0%
if 9.9999999999999998e284 < (-.f64 (+.f64 (*.f64 x y) (/.f64 (*.f64 z t) #s(literal 16 binary64))) (/.f64 (*.f64 a b) #s(literal 4 binary64))) Initial program 83.9%
+-commutativeN/A
sub-negN/A
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f64N/A
distribute-neg-frac2N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
metadata-evalN/A
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f6483.9%
Simplified83.9%
Taylor expanded in b around -inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
Simplified92.7%
Final simplification98.4%
(FPCore (x y z t a b c)
:precision binary64
(let* ((t_1 (+ (* x y) c)))
(if (<= (* a b) -1e+152)
(* b (+ (* a -0.25) (/ c b)))
(if (<= (* a b) -5e-197)
t_1
(if (<= (* a b) 2e-301)
(+ (* x y) (* 0.0625 (* z t)))
(if (<= (* a b) 1e+46) t_1 (+ (* x y) (/ (* a b) -4.0))))))))
double code(double x, double y, double z, double t, double a, double b, double c) {
double t_1 = (x * y) + c;
double tmp;
if ((a * b) <= -1e+152) {
tmp = b * ((a * -0.25) + (c / b));
} else if ((a * b) <= -5e-197) {
tmp = t_1;
} else if ((a * b) <= 2e-301) {
tmp = (x * y) + (0.0625 * (z * t));
} else if ((a * b) <= 1e+46) {
tmp = t_1;
} else {
tmp = (x * y) + ((a * b) / -4.0);
}
return tmp;
}
real(8) function code(x, y, z, t, a, b, c)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8), intent (in) :: c
real(8) :: t_1
real(8) :: tmp
t_1 = (x * y) + c
if ((a * b) <= (-1d+152)) then
tmp = b * ((a * (-0.25d0)) + (c / b))
else if ((a * b) <= (-5d-197)) then
tmp = t_1
else if ((a * b) <= 2d-301) then
tmp = (x * y) + (0.0625d0 * (z * t))
else if ((a * b) <= 1d+46) then
tmp = t_1
else
tmp = (x * y) + ((a * b) / (-4.0d0))
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b, double c) {
double t_1 = (x * y) + c;
double tmp;
if ((a * b) <= -1e+152) {
tmp = b * ((a * -0.25) + (c / b));
} else if ((a * b) <= -5e-197) {
tmp = t_1;
} else if ((a * b) <= 2e-301) {
tmp = (x * y) + (0.0625 * (z * t));
} else if ((a * b) <= 1e+46) {
tmp = t_1;
} else {
tmp = (x * y) + ((a * b) / -4.0);
}
return tmp;
}
def code(x, y, z, t, a, b, c): t_1 = (x * y) + c tmp = 0 if (a * b) <= -1e+152: tmp = b * ((a * -0.25) + (c / b)) elif (a * b) <= -5e-197: tmp = t_1 elif (a * b) <= 2e-301: tmp = (x * y) + (0.0625 * (z * t)) elif (a * b) <= 1e+46: tmp = t_1 else: tmp = (x * y) + ((a * b) / -4.0) return tmp
function code(x, y, z, t, a, b, c) t_1 = Float64(Float64(x * y) + c) tmp = 0.0 if (Float64(a * b) <= -1e+152) tmp = Float64(b * Float64(Float64(a * -0.25) + Float64(c / b))); elseif (Float64(a * b) <= -5e-197) tmp = t_1; elseif (Float64(a * b) <= 2e-301) tmp = Float64(Float64(x * y) + Float64(0.0625 * Float64(z * t))); elseif (Float64(a * b) <= 1e+46) tmp = t_1; else tmp = Float64(Float64(x * y) + Float64(Float64(a * b) / -4.0)); end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c) t_1 = (x * y) + c; tmp = 0.0; if ((a * b) <= -1e+152) tmp = b * ((a * -0.25) + (c / b)); elseif ((a * b) <= -5e-197) tmp = t_1; elseif ((a * b) <= 2e-301) tmp = (x * y) + (0.0625 * (z * t)); elseif ((a * b) <= 1e+46) tmp = t_1; else tmp = (x * y) + ((a * b) / -4.0); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_] := Block[{t$95$1 = N[(N[(x * y), $MachinePrecision] + c), $MachinePrecision]}, If[LessEqual[N[(a * b), $MachinePrecision], -1e+152], N[(b * N[(N[(a * -0.25), $MachinePrecision] + N[(c / b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[N[(a * b), $MachinePrecision], -5e-197], t$95$1, If[LessEqual[N[(a * b), $MachinePrecision], 2e-301], N[(N[(x * y), $MachinePrecision] + N[(0.0625 * N[(z * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[N[(a * b), $MachinePrecision], 1e+46], t$95$1, N[(N[(x * y), $MachinePrecision] + N[(N[(a * b), $MachinePrecision] / -4.0), $MachinePrecision]), $MachinePrecision]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x \cdot y + c\\
\mathbf{if}\;a \cdot b \leq -1 \cdot 10^{+152}:\\
\;\;\;\;b \cdot \left(a \cdot -0.25 + \frac{c}{b}\right)\\
\mathbf{elif}\;a \cdot b \leq -5 \cdot 10^{-197}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;a \cdot b \leq 2 \cdot 10^{-301}:\\
\;\;\;\;x \cdot y + 0.0625 \cdot \left(z \cdot t\right)\\
\mathbf{elif}\;a \cdot b \leq 10^{+46}:\\
\;\;\;\;t\_1\\
\mathbf{else}:\\
\;\;\;\;x \cdot y + \frac{a \cdot b}{-4}\\
\end{array}
\end{array}
if (*.f64 a b) < -1e152Initial program 90.4%
+-commutativeN/A
sub-negN/A
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f64N/A
distribute-neg-frac2N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
metadata-evalN/A
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f6490.4%
Simplified90.4%
Taylor expanded in c around inf
Simplified79.0%
Taylor expanded in b around inf
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f6481.1%
Simplified81.1%
if -1e152 < (*.f64 a b) < -5.0000000000000002e-197 or 2.00000000000000013e-301 < (*.f64 a b) < 9.9999999999999999e45Initial program 100.0%
+-commutativeN/A
sub-negN/A
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f64N/A
distribute-neg-frac2N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
metadata-evalN/A
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64100.0%
Simplified100.0%
Taylor expanded in a around 0
+-commutativeN/A
associate-+l+N/A
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f6495.6%
Simplified95.6%
Taylor expanded in t around 0
+-lowering-+.f64N/A
*-lowering-*.f6473.2%
Simplified73.2%
if -5.0000000000000002e-197 < (*.f64 a b) < 2.00000000000000013e-301Initial program 100.0%
+-commutativeN/A
sub-negN/A
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f64N/A
distribute-neg-frac2N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
metadata-evalN/A
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64100.0%
Simplified100.0%
Taylor expanded in a around 0
+-commutativeN/A
associate-+l+N/A
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64100.0%
Simplified100.0%
Taylor expanded in c around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6479.8%
Simplified79.8%
if 9.9999999999999999e45 < (*.f64 a b) Initial program 90.1%
+-commutativeN/A
sub-negN/A
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f64N/A
distribute-neg-frac2N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
metadata-evalN/A
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f6490.1%
Simplified90.1%
Taylor expanded in x around inf
*-lowering-*.f6473.4%
Simplified73.4%
Final simplification75.8%
(FPCore (x y z t a b c)
:precision binary64
(let* ((t_1 (+ (* x y) c)) (t_2 (* b (+ (* a -0.25) (/ c b)))))
(if (<= (* a b) -1e+152)
t_2
(if (<= (* a b) -5e-197)
t_1
(if (<= (* a b) 2e-301)
(+ (* x y) (* 0.0625 (* z t)))
(if (<= (* a b) 1e+46) t_1 t_2))))))
double code(double x, double y, double z, double t, double a, double b, double c) {
double t_1 = (x * y) + c;
double t_2 = b * ((a * -0.25) + (c / b));
double tmp;
if ((a * b) <= -1e+152) {
tmp = t_2;
} else if ((a * b) <= -5e-197) {
tmp = t_1;
} else if ((a * b) <= 2e-301) {
tmp = (x * y) + (0.0625 * (z * t));
} else if ((a * b) <= 1e+46) {
tmp = t_1;
} else {
tmp = t_2;
}
return tmp;
}
real(8) function code(x, y, z, t, a, b, c)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8), intent (in) :: c
real(8) :: t_1
real(8) :: t_2
real(8) :: tmp
t_1 = (x * y) + c
t_2 = b * ((a * (-0.25d0)) + (c / b))
if ((a * b) <= (-1d+152)) then
tmp = t_2
else if ((a * b) <= (-5d-197)) then
tmp = t_1
else if ((a * b) <= 2d-301) then
tmp = (x * y) + (0.0625d0 * (z * t))
else if ((a * b) <= 1d+46) then
tmp = t_1
else
tmp = t_2
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b, double c) {
double t_1 = (x * y) + c;
double t_2 = b * ((a * -0.25) + (c / b));
double tmp;
if ((a * b) <= -1e+152) {
tmp = t_2;
} else if ((a * b) <= -5e-197) {
tmp = t_1;
} else if ((a * b) <= 2e-301) {
tmp = (x * y) + (0.0625 * (z * t));
} else if ((a * b) <= 1e+46) {
tmp = t_1;
} else {
tmp = t_2;
}
return tmp;
}
def code(x, y, z, t, a, b, c): t_1 = (x * y) + c t_2 = b * ((a * -0.25) + (c / b)) tmp = 0 if (a * b) <= -1e+152: tmp = t_2 elif (a * b) <= -5e-197: tmp = t_1 elif (a * b) <= 2e-301: tmp = (x * y) + (0.0625 * (z * t)) elif (a * b) <= 1e+46: tmp = t_1 else: tmp = t_2 return tmp
function code(x, y, z, t, a, b, c) t_1 = Float64(Float64(x * y) + c) t_2 = Float64(b * Float64(Float64(a * -0.25) + Float64(c / b))) tmp = 0.0 if (Float64(a * b) <= -1e+152) tmp = t_2; elseif (Float64(a * b) <= -5e-197) tmp = t_1; elseif (Float64(a * b) <= 2e-301) tmp = Float64(Float64(x * y) + Float64(0.0625 * Float64(z * t))); elseif (Float64(a * b) <= 1e+46) tmp = t_1; else tmp = t_2; end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c) t_1 = (x * y) + c; t_2 = b * ((a * -0.25) + (c / b)); tmp = 0.0; if ((a * b) <= -1e+152) tmp = t_2; elseif ((a * b) <= -5e-197) tmp = t_1; elseif ((a * b) <= 2e-301) tmp = (x * y) + (0.0625 * (z * t)); elseif ((a * b) <= 1e+46) tmp = t_1; else tmp = t_2; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_] := Block[{t$95$1 = N[(N[(x * y), $MachinePrecision] + c), $MachinePrecision]}, Block[{t$95$2 = N[(b * N[(N[(a * -0.25), $MachinePrecision] + N[(c / b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[N[(a * b), $MachinePrecision], -1e+152], t$95$2, If[LessEqual[N[(a * b), $MachinePrecision], -5e-197], t$95$1, If[LessEqual[N[(a * b), $MachinePrecision], 2e-301], N[(N[(x * y), $MachinePrecision] + N[(0.0625 * N[(z * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[N[(a * b), $MachinePrecision], 1e+46], t$95$1, t$95$2]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x \cdot y + c\\
t_2 := b \cdot \left(a \cdot -0.25 + \frac{c}{b}\right)\\
\mathbf{if}\;a \cdot b \leq -1 \cdot 10^{+152}:\\
\;\;\;\;t\_2\\
\mathbf{elif}\;a \cdot b \leq -5 \cdot 10^{-197}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;a \cdot b \leq 2 \cdot 10^{-301}:\\
\;\;\;\;x \cdot y + 0.0625 \cdot \left(z \cdot t\right)\\
\mathbf{elif}\;a \cdot b \leq 10^{+46}:\\
\;\;\;\;t\_1\\
\mathbf{else}:\\
\;\;\;\;t\_2\\
\end{array}
\end{array}
if (*.f64 a b) < -1e152 or 9.9999999999999999e45 < (*.f64 a b) Initial program 90.2%
+-commutativeN/A
sub-negN/A
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f64N/A
distribute-neg-frac2N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
metadata-evalN/A
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f6490.2%
Simplified90.2%
Taylor expanded in c around inf
Simplified75.6%
Taylor expanded in b around inf
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f6476.6%
Simplified76.6%
if -1e152 < (*.f64 a b) < -5.0000000000000002e-197 or 2.00000000000000013e-301 < (*.f64 a b) < 9.9999999999999999e45Initial program 100.0%
+-commutativeN/A
sub-negN/A
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f64N/A
distribute-neg-frac2N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
metadata-evalN/A
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64100.0%
Simplified100.0%
Taylor expanded in a around 0
+-commutativeN/A
associate-+l+N/A
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f6495.6%
Simplified95.6%
Taylor expanded in t around 0
+-lowering-+.f64N/A
*-lowering-*.f6473.2%
Simplified73.2%
if -5.0000000000000002e-197 < (*.f64 a b) < 2.00000000000000013e-301Initial program 100.0%
+-commutativeN/A
sub-negN/A
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f64N/A
distribute-neg-frac2N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
metadata-evalN/A
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64100.0%
Simplified100.0%
Taylor expanded in a around 0
+-commutativeN/A
associate-+l+N/A
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64100.0%
Simplified100.0%
Taylor expanded in c around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6479.8%
Simplified79.8%
Final simplification75.8%
(FPCore (x y z t a b c)
:precision binary64
(let* ((t_1 (/ (* z t) 16.0)))
(if (<= (- (+ (* x y) t_1) (/ (* a b) 4.0)) INFINITY)
(+ (/ a (/ -4.0 b)) (+ t_1 (+ (* x y) c)))
(* 0.0625 (* z t)))))
double code(double x, double y, double z, double t, double a, double b, double c) {
double t_1 = (z * t) / 16.0;
double tmp;
if ((((x * y) + t_1) - ((a * b) / 4.0)) <= ((double) INFINITY)) {
tmp = (a / (-4.0 / b)) + (t_1 + ((x * y) + c));
} else {
tmp = 0.0625 * (z * t);
}
return tmp;
}
public static double code(double x, double y, double z, double t, double a, double b, double c) {
double t_1 = (z * t) / 16.0;
double tmp;
if ((((x * y) + t_1) - ((a * b) / 4.0)) <= Double.POSITIVE_INFINITY) {
tmp = (a / (-4.0 / b)) + (t_1 + ((x * y) + c));
} else {
tmp = 0.0625 * (z * t);
}
return tmp;
}
def code(x, y, z, t, a, b, c): t_1 = (z * t) / 16.0 tmp = 0 if (((x * y) + t_1) - ((a * b) / 4.0)) <= math.inf: tmp = (a / (-4.0 / b)) + (t_1 + ((x * y) + c)) else: tmp = 0.0625 * (z * t) return tmp
function code(x, y, z, t, a, b, c) t_1 = Float64(Float64(z * t) / 16.0) tmp = 0.0 if (Float64(Float64(Float64(x * y) + t_1) - Float64(Float64(a * b) / 4.0)) <= Inf) tmp = Float64(Float64(a / Float64(-4.0 / b)) + Float64(t_1 + Float64(Float64(x * y) + c))); else tmp = Float64(0.0625 * Float64(z * t)); end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c) t_1 = (z * t) / 16.0; tmp = 0.0; if ((((x * y) + t_1) - ((a * b) / 4.0)) <= Inf) tmp = (a / (-4.0 / b)) + (t_1 + ((x * y) + c)); else tmp = 0.0625 * (z * t); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_] := Block[{t$95$1 = N[(N[(z * t), $MachinePrecision] / 16.0), $MachinePrecision]}, If[LessEqual[N[(N[(N[(x * y), $MachinePrecision] + t$95$1), $MachinePrecision] - N[(N[(a * b), $MachinePrecision] / 4.0), $MachinePrecision]), $MachinePrecision], Infinity], N[(N[(a / N[(-4.0 / b), $MachinePrecision]), $MachinePrecision] + N[(t$95$1 + N[(N[(x * y), $MachinePrecision] + c), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(0.0625 * N[(z * t), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{z \cdot t}{16}\\
\mathbf{if}\;\left(x \cdot y + t\_1\right) - \frac{a \cdot b}{4} \leq \infty:\\
\;\;\;\;\frac{a}{\frac{-4}{b}} + \left(t\_1 + \left(x \cdot y + c\right)\right)\\
\mathbf{else}:\\
\;\;\;\;0.0625 \cdot \left(z \cdot t\right)\\
\end{array}
\end{array}
if (-.f64 (+.f64 (*.f64 x y) (/.f64 (*.f64 z t) #s(literal 16 binary64))) (/.f64 (*.f64 a b) #s(literal 4 binary64))) < +inf.0Initial program 99.6%
+-commutativeN/A
sub-negN/A
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f64N/A
distribute-neg-frac2N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
metadata-evalN/A
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f6499.6%
Simplified99.6%
associate-/l*N/A
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
/-lowering-/.f6499.9%
Applied egg-rr99.9%
if +inf.0 < (-.f64 (+.f64 (*.f64 x y) (/.f64 (*.f64 z t) #s(literal 16 binary64))) (/.f64 (*.f64 a b) #s(literal 4 binary64))) Initial program 0.0%
+-commutativeN/A
sub-negN/A
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f64N/A
distribute-neg-frac2N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
metadata-evalN/A
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f640.0%
Simplified0.0%
Taylor expanded in z around inf
*-lowering-*.f64N/A
*-lowering-*.f6450.1%
Simplified50.1%
Final simplification98.4%
(FPCore (x y z t a b c) :precision binary64 (let* ((t_1 (- (+ (* x y) (/ (* z t) 16.0)) (/ (* a b) 4.0)))) (if (<= t_1 INFINITY) (+ t_1 c) (* 0.0625 (* z t)))))
double code(double x, double y, double z, double t, double a, double b, double c) {
double t_1 = ((x * y) + ((z * t) / 16.0)) - ((a * b) / 4.0);
double tmp;
if (t_1 <= ((double) INFINITY)) {
tmp = t_1 + c;
} else {
tmp = 0.0625 * (z * t);
}
return tmp;
}
public static double code(double x, double y, double z, double t, double a, double b, double c) {
double t_1 = ((x * y) + ((z * t) / 16.0)) - ((a * b) / 4.0);
double tmp;
if (t_1 <= Double.POSITIVE_INFINITY) {
tmp = t_1 + c;
} else {
tmp = 0.0625 * (z * t);
}
return tmp;
}
def code(x, y, z, t, a, b, c): t_1 = ((x * y) + ((z * t) / 16.0)) - ((a * b) / 4.0) tmp = 0 if t_1 <= math.inf: tmp = t_1 + c else: tmp = 0.0625 * (z * t) return tmp
function code(x, y, z, t, a, b, c) t_1 = Float64(Float64(Float64(x * y) + Float64(Float64(z * t) / 16.0)) - Float64(Float64(a * b) / 4.0)) tmp = 0.0 if (t_1 <= Inf) tmp = Float64(t_1 + c); else tmp = Float64(0.0625 * Float64(z * t)); end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c) t_1 = ((x * y) + ((z * t) / 16.0)) - ((a * b) / 4.0); tmp = 0.0; if (t_1 <= Inf) tmp = t_1 + c; else tmp = 0.0625 * (z * t); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_] := Block[{t$95$1 = N[(N[(N[(x * y), $MachinePrecision] + N[(N[(z * t), $MachinePrecision] / 16.0), $MachinePrecision]), $MachinePrecision] - N[(N[(a * b), $MachinePrecision] / 4.0), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$1, Infinity], N[(t$95$1 + c), $MachinePrecision], N[(0.0625 * N[(z * t), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \left(x \cdot y + \frac{z \cdot t}{16}\right) - \frac{a \cdot b}{4}\\
\mathbf{if}\;t\_1 \leq \infty:\\
\;\;\;\;t\_1 + c\\
\mathbf{else}:\\
\;\;\;\;0.0625 \cdot \left(z \cdot t\right)\\
\end{array}
\end{array}
if (-.f64 (+.f64 (*.f64 x y) (/.f64 (*.f64 z t) #s(literal 16 binary64))) (/.f64 (*.f64 a b) #s(literal 4 binary64))) < +inf.0Initial program 99.6%
if +inf.0 < (-.f64 (+.f64 (*.f64 x y) (/.f64 (*.f64 z t) #s(literal 16 binary64))) (/.f64 (*.f64 a b) #s(literal 4 binary64))) Initial program 0.0%
+-commutativeN/A
sub-negN/A
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f64N/A
distribute-neg-frac2N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
metadata-evalN/A
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f640.0%
Simplified0.0%
Taylor expanded in z around inf
*-lowering-*.f64N/A
*-lowering-*.f6450.1%
Simplified50.1%
Final simplification98.1%
(FPCore (x y z t a b c)
:precision binary64
(if (<= (* a b) -1e+152)
(* a (/ b -4.0))
(if (<= (* a b) 1e-140)
(* x y)
(if (<= (* a b) 1e+46) c (* (* a b) -0.25)))))
double code(double x, double y, double z, double t, double a, double b, double c) {
double tmp;
if ((a * b) <= -1e+152) {
tmp = a * (b / -4.0);
} else if ((a * b) <= 1e-140) {
tmp = x * y;
} else if ((a * b) <= 1e+46) {
tmp = c;
} else {
tmp = (a * b) * -0.25;
}
return tmp;
}
real(8) function code(x, y, z, t, a, b, c)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8), intent (in) :: c
real(8) :: tmp
if ((a * b) <= (-1d+152)) then
tmp = a * (b / (-4.0d0))
else if ((a * b) <= 1d-140) then
tmp = x * y
else if ((a * b) <= 1d+46) then
tmp = c
else
tmp = (a * b) * (-0.25d0)
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b, double c) {
double tmp;
if ((a * b) <= -1e+152) {
tmp = a * (b / -4.0);
} else if ((a * b) <= 1e-140) {
tmp = x * y;
} else if ((a * b) <= 1e+46) {
tmp = c;
} else {
tmp = (a * b) * -0.25;
}
return tmp;
}
def code(x, y, z, t, a, b, c): tmp = 0 if (a * b) <= -1e+152: tmp = a * (b / -4.0) elif (a * b) <= 1e-140: tmp = x * y elif (a * b) <= 1e+46: tmp = c else: tmp = (a * b) * -0.25 return tmp
function code(x, y, z, t, a, b, c) tmp = 0.0 if (Float64(a * b) <= -1e+152) tmp = Float64(a * Float64(b / -4.0)); elseif (Float64(a * b) <= 1e-140) tmp = Float64(x * y); elseif (Float64(a * b) <= 1e+46) tmp = c; else tmp = Float64(Float64(a * b) * -0.25); end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c) tmp = 0.0; if ((a * b) <= -1e+152) tmp = a * (b / -4.0); elseif ((a * b) <= 1e-140) tmp = x * y; elseif ((a * b) <= 1e+46) tmp = c; else tmp = (a * b) * -0.25; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_] := If[LessEqual[N[(a * b), $MachinePrecision], -1e+152], N[(a * N[(b / -4.0), $MachinePrecision]), $MachinePrecision], If[LessEqual[N[(a * b), $MachinePrecision], 1e-140], N[(x * y), $MachinePrecision], If[LessEqual[N[(a * b), $MachinePrecision], 1e+46], c, N[(N[(a * b), $MachinePrecision] * -0.25), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;a \cdot b \leq -1 \cdot 10^{+152}:\\
\;\;\;\;a \cdot \frac{b}{-4}\\
\mathbf{elif}\;a \cdot b \leq 10^{-140}:\\
\;\;\;\;x \cdot y\\
\mathbf{elif}\;a \cdot b \leq 10^{+46}:\\
\;\;\;\;c\\
\mathbf{else}:\\
\;\;\;\;\left(a \cdot b\right) \cdot -0.25\\
\end{array}
\end{array}
if (*.f64 a b) < -1e152Initial program 90.4%
+-commutativeN/A
sub-negN/A
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f64N/A
distribute-neg-frac2N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
metadata-evalN/A
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f6490.4%
Simplified90.4%
Taylor expanded in a around inf
*-lowering-*.f64N/A
*-lowering-*.f6469.8%
Simplified69.8%
*-commutativeN/A
associate-*r*N/A
metadata-evalN/A
associate-/r/N/A
clear-numN/A
*-lowering-*.f64N/A
/-lowering-/.f6471.9%
Applied egg-rr71.9%
if -1e152 < (*.f64 a b) < 9.9999999999999998e-141Initial program 100.0%
+-commutativeN/A
sub-negN/A
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f64N/A
distribute-neg-frac2N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
metadata-evalN/A
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64100.0%
Simplified100.0%
Taylor expanded in x around inf
*-lowering-*.f6442.1%
Simplified42.1%
if 9.9999999999999998e-141 < (*.f64 a b) < 9.9999999999999999e45Initial program 100.0%
+-commutativeN/A
sub-negN/A
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f64N/A
distribute-neg-frac2N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
metadata-evalN/A
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64100.0%
Simplified100.0%
Taylor expanded in c around inf
Simplified45.7%
if 9.9999999999999999e45 < (*.f64 a b) Initial program 90.1%
+-commutativeN/A
sub-negN/A
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f64N/A
distribute-neg-frac2N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
metadata-evalN/A
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f6490.1%
Simplified90.1%
Taylor expanded in a around inf
*-lowering-*.f64N/A
*-lowering-*.f6469.3%
Simplified69.3%
Final simplification52.8%
(FPCore (x y z t a b c)
:precision binary64
(let* ((t_1 (* (* a b) -0.25)))
(if (<= (* a b) -3.7e+137)
t_1
(if (<= (* a b) 1e-140) (* x y) (if (<= (* a b) 1e+52) c t_1)))))
double code(double x, double y, double z, double t, double a, double b, double c) {
double t_1 = (a * b) * -0.25;
double tmp;
if ((a * b) <= -3.7e+137) {
tmp = t_1;
} else if ((a * b) <= 1e-140) {
tmp = x * y;
} else if ((a * b) <= 1e+52) {
tmp = c;
} else {
tmp = t_1;
}
return tmp;
}
real(8) function code(x, y, z, t, a, b, c)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8), intent (in) :: c
real(8) :: t_1
real(8) :: tmp
t_1 = (a * b) * (-0.25d0)
if ((a * b) <= (-3.7d+137)) then
tmp = t_1
else if ((a * b) <= 1d-140) then
tmp = x * y
else if ((a * b) <= 1d+52) then
tmp = c
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 c) {
double t_1 = (a * b) * -0.25;
double tmp;
if ((a * b) <= -3.7e+137) {
tmp = t_1;
} else if ((a * b) <= 1e-140) {
tmp = x * y;
} else if ((a * b) <= 1e+52) {
tmp = c;
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b, c): t_1 = (a * b) * -0.25 tmp = 0 if (a * b) <= -3.7e+137: tmp = t_1 elif (a * b) <= 1e-140: tmp = x * y elif (a * b) <= 1e+52: tmp = c else: tmp = t_1 return tmp
function code(x, y, z, t, a, b, c) t_1 = Float64(Float64(a * b) * -0.25) tmp = 0.0 if (Float64(a * b) <= -3.7e+137) tmp = t_1; elseif (Float64(a * b) <= 1e-140) tmp = Float64(x * y); elseif (Float64(a * b) <= 1e+52) tmp = c; else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c) t_1 = (a * b) * -0.25; tmp = 0.0; if ((a * b) <= -3.7e+137) tmp = t_1; elseif ((a * b) <= 1e-140) tmp = x * y; elseif ((a * b) <= 1e+52) tmp = c; else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_] := Block[{t$95$1 = N[(N[(a * b), $MachinePrecision] * -0.25), $MachinePrecision]}, If[LessEqual[N[(a * b), $MachinePrecision], -3.7e+137], t$95$1, If[LessEqual[N[(a * b), $MachinePrecision], 1e-140], N[(x * y), $MachinePrecision], If[LessEqual[N[(a * b), $MachinePrecision], 1e+52], c, t$95$1]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \left(a \cdot b\right) \cdot -0.25\\
\mathbf{if}\;a \cdot b \leq -3.7 \cdot 10^{+137}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;a \cdot b \leq 10^{-140}:\\
\;\;\;\;x \cdot y\\
\mathbf{elif}\;a \cdot b \leq 10^{+52}:\\
\;\;\;\;c\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if (*.f64 a b) < -3.7000000000000002e137 or 9.9999999999999999e51 < (*.f64 a b) Initial program 90.2%
+-commutativeN/A
sub-negN/A
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f64N/A
distribute-neg-frac2N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
metadata-evalN/A
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f6490.2%
Simplified90.2%
Taylor expanded in a around inf
*-lowering-*.f64N/A
*-lowering-*.f6469.5%
Simplified69.5%
if -3.7000000000000002e137 < (*.f64 a b) < 9.9999999999999998e-141Initial program 100.0%
+-commutativeN/A
sub-negN/A
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f64N/A
distribute-neg-frac2N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
metadata-evalN/A
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64100.0%
Simplified100.0%
Taylor expanded in x around inf
*-lowering-*.f6442.1%
Simplified42.1%
if 9.9999999999999998e-141 < (*.f64 a b) < 9.9999999999999999e51Initial program 100.0%
+-commutativeN/A
sub-negN/A
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f64N/A
distribute-neg-frac2N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
metadata-evalN/A
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64100.0%
Simplified100.0%
Taylor expanded in c around inf
Simplified45.7%
Final simplification52.4%
(FPCore (x y z t a b c)
:precision binary64
(let* ((t_1 (* 0.0625 (* z t))))
(if (<= (* a b) -1e+203)
(* b (+ (* a -0.25) (/ c b)))
(if (<= (* a b) 1e+46) (+ t_1 (+ (* x y) c)) (+ t_1 (/ (* a b) -4.0))))))
double code(double x, double y, double z, double t, double a, double b, double c) {
double t_1 = 0.0625 * (z * t);
double tmp;
if ((a * b) <= -1e+203) {
tmp = b * ((a * -0.25) + (c / b));
} else if ((a * b) <= 1e+46) {
tmp = t_1 + ((x * y) + c);
} else {
tmp = t_1 + ((a * b) / -4.0);
}
return tmp;
}
real(8) function code(x, y, z, t, a, b, c)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8), intent (in) :: c
real(8) :: t_1
real(8) :: tmp
t_1 = 0.0625d0 * (z * t)
if ((a * b) <= (-1d+203)) then
tmp = b * ((a * (-0.25d0)) + (c / b))
else if ((a * b) <= 1d+46) then
tmp = t_1 + ((x * y) + c)
else
tmp = t_1 + ((a * b) / (-4.0d0))
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b, double c) {
double t_1 = 0.0625 * (z * t);
double tmp;
if ((a * b) <= -1e+203) {
tmp = b * ((a * -0.25) + (c / b));
} else if ((a * b) <= 1e+46) {
tmp = t_1 + ((x * y) + c);
} else {
tmp = t_1 + ((a * b) / -4.0);
}
return tmp;
}
def code(x, y, z, t, a, b, c): t_1 = 0.0625 * (z * t) tmp = 0 if (a * b) <= -1e+203: tmp = b * ((a * -0.25) + (c / b)) elif (a * b) <= 1e+46: tmp = t_1 + ((x * y) + c) else: tmp = t_1 + ((a * b) / -4.0) return tmp
function code(x, y, z, t, a, b, c) t_1 = Float64(0.0625 * Float64(z * t)) tmp = 0.0 if (Float64(a * b) <= -1e+203) tmp = Float64(b * Float64(Float64(a * -0.25) + Float64(c / b))); elseif (Float64(a * b) <= 1e+46) tmp = Float64(t_1 + Float64(Float64(x * y) + c)); else tmp = Float64(t_1 + Float64(Float64(a * b) / -4.0)); end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c) t_1 = 0.0625 * (z * t); tmp = 0.0; if ((a * b) <= -1e+203) tmp = b * ((a * -0.25) + (c / b)); elseif ((a * b) <= 1e+46) tmp = t_1 + ((x * y) + c); else tmp = t_1 + ((a * b) / -4.0); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_] := Block[{t$95$1 = N[(0.0625 * N[(z * t), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[N[(a * b), $MachinePrecision], -1e+203], N[(b * N[(N[(a * -0.25), $MachinePrecision] + N[(c / b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[N[(a * b), $MachinePrecision], 1e+46], N[(t$95$1 + N[(N[(x * y), $MachinePrecision] + c), $MachinePrecision]), $MachinePrecision], N[(t$95$1 + N[(N[(a * b), $MachinePrecision] / -4.0), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := 0.0625 \cdot \left(z \cdot t\right)\\
\mathbf{if}\;a \cdot b \leq -1 \cdot 10^{+203}:\\
\;\;\;\;b \cdot \left(a \cdot -0.25 + \frac{c}{b}\right)\\
\mathbf{elif}\;a \cdot b \leq 10^{+46}:\\
\;\;\;\;t\_1 + \left(x \cdot y + c\right)\\
\mathbf{else}:\\
\;\;\;\;t\_1 + \frac{a \cdot b}{-4}\\
\end{array}
\end{array}
if (*.f64 a b) < -9.9999999999999999e202Initial program 88.0%
+-commutativeN/A
sub-negN/A
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f64N/A
distribute-neg-frac2N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
metadata-evalN/A
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f6488.0%
Simplified88.0%
Taylor expanded in c around inf
Simplified82.6%
Taylor expanded in b around inf
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f6485.2%
Simplified85.2%
if -9.9999999999999999e202 < (*.f64 a b) < 9.9999999999999999e45Initial program 100.0%
+-commutativeN/A
sub-negN/A
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f64N/A
distribute-neg-frac2N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
metadata-evalN/A
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64100.0%
Simplified100.0%
Taylor expanded in a around 0
+-commutativeN/A
associate-+l+N/A
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f6495.7%
Simplified95.7%
if 9.9999999999999999e45 < (*.f64 a b) Initial program 90.1%
+-commutativeN/A
sub-negN/A
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f64N/A
distribute-neg-frac2N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
metadata-evalN/A
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f6490.1%
Simplified90.1%
Taylor expanded in z around inf
*-lowering-*.f64N/A
*-lowering-*.f6484.5%
Simplified84.5%
Final simplification92.2%
(FPCore (x y z t a b c)
:precision binary64
(if (<= (* a b) -1e+203)
(* b (+ (* a -0.25) (/ c b)))
(if (<= (* a b) 1e+56)
(+ (* 0.0625 (* z t)) (+ (* x y) c))
(+ (* x y) (/ (* a b) -4.0)))))
double code(double x, double y, double z, double t, double a, double b, double c) {
double tmp;
if ((a * b) <= -1e+203) {
tmp = b * ((a * -0.25) + (c / b));
} else if ((a * b) <= 1e+56) {
tmp = (0.0625 * (z * t)) + ((x * y) + c);
} else {
tmp = (x * y) + ((a * b) / -4.0);
}
return tmp;
}
real(8) function code(x, y, z, t, a, b, c)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8), intent (in) :: c
real(8) :: tmp
if ((a * b) <= (-1d+203)) then
tmp = b * ((a * (-0.25d0)) + (c / b))
else if ((a * b) <= 1d+56) then
tmp = (0.0625d0 * (z * t)) + ((x * y) + c)
else
tmp = (x * y) + ((a * b) / (-4.0d0))
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b, double c) {
double tmp;
if ((a * b) <= -1e+203) {
tmp = b * ((a * -0.25) + (c / b));
} else if ((a * b) <= 1e+56) {
tmp = (0.0625 * (z * t)) + ((x * y) + c);
} else {
tmp = (x * y) + ((a * b) / -4.0);
}
return tmp;
}
def code(x, y, z, t, a, b, c): tmp = 0 if (a * b) <= -1e+203: tmp = b * ((a * -0.25) + (c / b)) elif (a * b) <= 1e+56: tmp = (0.0625 * (z * t)) + ((x * y) + c) else: tmp = (x * y) + ((a * b) / -4.0) return tmp
function code(x, y, z, t, a, b, c) tmp = 0.0 if (Float64(a * b) <= -1e+203) tmp = Float64(b * Float64(Float64(a * -0.25) + Float64(c / b))); elseif (Float64(a * b) <= 1e+56) tmp = Float64(Float64(0.0625 * Float64(z * t)) + Float64(Float64(x * y) + c)); else tmp = Float64(Float64(x * y) + Float64(Float64(a * b) / -4.0)); end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c) tmp = 0.0; if ((a * b) <= -1e+203) tmp = b * ((a * -0.25) + (c / b)); elseif ((a * b) <= 1e+56) tmp = (0.0625 * (z * t)) + ((x * y) + c); else tmp = (x * y) + ((a * b) / -4.0); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_] := If[LessEqual[N[(a * b), $MachinePrecision], -1e+203], N[(b * N[(N[(a * -0.25), $MachinePrecision] + N[(c / b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[N[(a * b), $MachinePrecision], 1e+56], N[(N[(0.0625 * N[(z * t), $MachinePrecision]), $MachinePrecision] + N[(N[(x * y), $MachinePrecision] + c), $MachinePrecision]), $MachinePrecision], N[(N[(x * y), $MachinePrecision] + N[(N[(a * b), $MachinePrecision] / -4.0), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;a \cdot b \leq -1 \cdot 10^{+203}:\\
\;\;\;\;b \cdot \left(a \cdot -0.25 + \frac{c}{b}\right)\\
\mathbf{elif}\;a \cdot b \leq 10^{+56}:\\
\;\;\;\;0.0625 \cdot \left(z \cdot t\right) + \left(x \cdot y + c\right)\\
\mathbf{else}:\\
\;\;\;\;x \cdot y + \frac{a \cdot b}{-4}\\
\end{array}
\end{array}
if (*.f64 a b) < -9.9999999999999999e202Initial program 88.0%
+-commutativeN/A
sub-negN/A
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f64N/A
distribute-neg-frac2N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
metadata-evalN/A
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f6488.0%
Simplified88.0%
Taylor expanded in c around inf
Simplified82.6%
Taylor expanded in b around inf
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f6485.2%
Simplified85.2%
if -9.9999999999999999e202 < (*.f64 a b) < 1.00000000000000009e56Initial program 100.0%
+-commutativeN/A
sub-negN/A
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f64N/A
distribute-neg-frac2N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
metadata-evalN/A
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64100.0%
Simplified100.0%
Taylor expanded in a around 0
+-commutativeN/A
associate-+l+N/A
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f6495.7%
Simplified95.7%
if 1.00000000000000009e56 < (*.f64 a b) Initial program 89.9%
+-commutativeN/A
sub-negN/A
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f64N/A
distribute-neg-frac2N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
metadata-evalN/A
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f6489.9%
Simplified89.9%
Taylor expanded in x around inf
*-lowering-*.f6474.9%
Simplified74.9%
Final simplification90.3%
(FPCore (x y z t a b c) :precision binary64 (let* ((t_1 (* b (+ (* a -0.25) (/ c b))))) (if (<= (* a b) -1e+152) t_1 (if (<= (* a b) 1e+46) (+ (* x y) c) t_1))))
double code(double x, double y, double z, double t, double a, double b, double c) {
double t_1 = b * ((a * -0.25) + (c / b));
double tmp;
if ((a * b) <= -1e+152) {
tmp = t_1;
} else if ((a * b) <= 1e+46) {
tmp = (x * y) + c;
} else {
tmp = t_1;
}
return tmp;
}
real(8) function code(x, y, z, t, a, b, c)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8), intent (in) :: c
real(8) :: t_1
real(8) :: tmp
t_1 = b * ((a * (-0.25d0)) + (c / b))
if ((a * b) <= (-1d+152)) then
tmp = t_1
else if ((a * b) <= 1d+46) then
tmp = (x * y) + c
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 c) {
double t_1 = b * ((a * -0.25) + (c / b));
double tmp;
if ((a * b) <= -1e+152) {
tmp = t_1;
} else if ((a * b) <= 1e+46) {
tmp = (x * y) + c;
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b, c): t_1 = b * ((a * -0.25) + (c / b)) tmp = 0 if (a * b) <= -1e+152: tmp = t_1 elif (a * b) <= 1e+46: tmp = (x * y) + c else: tmp = t_1 return tmp
function code(x, y, z, t, a, b, c) t_1 = Float64(b * Float64(Float64(a * -0.25) + Float64(c / b))) tmp = 0.0 if (Float64(a * b) <= -1e+152) tmp = t_1; elseif (Float64(a * b) <= 1e+46) tmp = Float64(Float64(x * y) + c); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c) t_1 = b * ((a * -0.25) + (c / b)); tmp = 0.0; if ((a * b) <= -1e+152) tmp = t_1; elseif ((a * b) <= 1e+46) tmp = (x * y) + c; else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_] := Block[{t$95$1 = N[(b * N[(N[(a * -0.25), $MachinePrecision] + N[(c / b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[N[(a * b), $MachinePrecision], -1e+152], t$95$1, If[LessEqual[N[(a * b), $MachinePrecision], 1e+46], N[(N[(x * y), $MachinePrecision] + c), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := b \cdot \left(a \cdot -0.25 + \frac{c}{b}\right)\\
\mathbf{if}\;a \cdot b \leq -1 \cdot 10^{+152}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;a \cdot b \leq 10^{+46}:\\
\;\;\;\;x \cdot y + c\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if (*.f64 a b) < -1e152 or 9.9999999999999999e45 < (*.f64 a b) Initial program 90.2%
+-commutativeN/A
sub-negN/A
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f64N/A
distribute-neg-frac2N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
metadata-evalN/A
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f6490.2%
Simplified90.2%
Taylor expanded in c around inf
Simplified75.6%
Taylor expanded in b around inf
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f6476.6%
Simplified76.6%
if -1e152 < (*.f64 a b) < 9.9999999999999999e45Initial program 100.0%
+-commutativeN/A
sub-negN/A
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f64N/A
distribute-neg-frac2N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
metadata-evalN/A
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64100.0%
Simplified100.0%
Taylor expanded in a around 0
+-commutativeN/A
associate-+l+N/A
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f6497.0%
Simplified97.0%
Taylor expanded in t around 0
+-lowering-+.f64N/A
*-lowering-*.f6470.0%
Simplified70.0%
Final simplification72.3%
(FPCore (x y z t a b c) :precision binary64 (let* ((t_1 (+ c (/ (* a b) -4.0)))) (if (<= (* a b) -1e+152) t_1 (if (<= (* a b) 1e+44) (+ (* x y) c) t_1))))
double code(double x, double y, double z, double t, double a, double b, double c) {
double t_1 = c + ((a * b) / -4.0);
double tmp;
if ((a * b) <= -1e+152) {
tmp = t_1;
} else if ((a * b) <= 1e+44) {
tmp = (x * y) + c;
} else {
tmp = t_1;
}
return tmp;
}
real(8) function code(x, y, z, t, a, b, c)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8), intent (in) :: c
real(8) :: t_1
real(8) :: tmp
t_1 = c + ((a * b) / (-4.0d0))
if ((a * b) <= (-1d+152)) then
tmp = t_1
else if ((a * b) <= 1d+44) then
tmp = (x * y) + c
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 c) {
double t_1 = c + ((a * b) / -4.0);
double tmp;
if ((a * b) <= -1e+152) {
tmp = t_1;
} else if ((a * b) <= 1e+44) {
tmp = (x * y) + c;
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b, c): t_1 = c + ((a * b) / -4.0) tmp = 0 if (a * b) <= -1e+152: tmp = t_1 elif (a * b) <= 1e+44: tmp = (x * y) + c else: tmp = t_1 return tmp
function code(x, y, z, t, a, b, c) t_1 = Float64(c + Float64(Float64(a * b) / -4.0)) tmp = 0.0 if (Float64(a * b) <= -1e+152) tmp = t_1; elseif (Float64(a * b) <= 1e+44) tmp = Float64(Float64(x * y) + c); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c) t_1 = c + ((a * b) / -4.0); tmp = 0.0; if ((a * b) <= -1e+152) tmp = t_1; elseif ((a * b) <= 1e+44) tmp = (x * y) + c; else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_] := Block[{t$95$1 = N[(c + N[(N[(a * b), $MachinePrecision] / -4.0), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[N[(a * b), $MachinePrecision], -1e+152], t$95$1, If[LessEqual[N[(a * b), $MachinePrecision], 1e+44], N[(N[(x * y), $MachinePrecision] + c), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := c + \frac{a \cdot b}{-4}\\
\mathbf{if}\;a \cdot b \leq -1 \cdot 10^{+152}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;a \cdot b \leq 10^{+44}:\\
\;\;\;\;x \cdot y + c\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if (*.f64 a b) < -1e152 or 1.0000000000000001e44 < (*.f64 a b) Initial program 90.4%
+-commutativeN/A
sub-negN/A
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f64N/A
distribute-neg-frac2N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
metadata-evalN/A
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f6490.4%
Simplified90.4%
Taylor expanded in c around inf
Simplified75.9%
if -1e152 < (*.f64 a b) < 1.0000000000000001e44Initial program 100.0%
+-commutativeN/A
sub-negN/A
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f64N/A
distribute-neg-frac2N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
metadata-evalN/A
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64100.0%
Simplified100.0%
Taylor expanded in a around 0
+-commutativeN/A
associate-+l+N/A
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f6497.0%
Simplified97.0%
Taylor expanded in t around 0
+-lowering-+.f64N/A
*-lowering-*.f6469.8%
Simplified69.8%
Final simplification72.0%
(FPCore (x y z t a b c) :precision binary64 (if (<= (* a b) -1e+203) (* a (/ b -4.0)) (if (<= (* a b) 1e+46) (+ (* x y) c) (* (* a b) -0.25))))
double code(double x, double y, double z, double t, double a, double b, double c) {
double tmp;
if ((a * b) <= -1e+203) {
tmp = a * (b / -4.0);
} else if ((a * b) <= 1e+46) {
tmp = (x * y) + c;
} else {
tmp = (a * b) * -0.25;
}
return tmp;
}
real(8) function code(x, y, z, t, a, b, c)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8), intent (in) :: c
real(8) :: tmp
if ((a * b) <= (-1d+203)) then
tmp = a * (b / (-4.0d0))
else if ((a * b) <= 1d+46) then
tmp = (x * y) + c
else
tmp = (a * b) * (-0.25d0)
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b, double c) {
double tmp;
if ((a * b) <= -1e+203) {
tmp = a * (b / -4.0);
} else if ((a * b) <= 1e+46) {
tmp = (x * y) + c;
} else {
tmp = (a * b) * -0.25;
}
return tmp;
}
def code(x, y, z, t, a, b, c): tmp = 0 if (a * b) <= -1e+203: tmp = a * (b / -4.0) elif (a * b) <= 1e+46: tmp = (x * y) + c else: tmp = (a * b) * -0.25 return tmp
function code(x, y, z, t, a, b, c) tmp = 0.0 if (Float64(a * b) <= -1e+203) tmp = Float64(a * Float64(b / -4.0)); elseif (Float64(a * b) <= 1e+46) tmp = Float64(Float64(x * y) + c); else tmp = Float64(Float64(a * b) * -0.25); end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c) tmp = 0.0; if ((a * b) <= -1e+203) tmp = a * (b / -4.0); elseif ((a * b) <= 1e+46) tmp = (x * y) + c; else tmp = (a * b) * -0.25; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_] := If[LessEqual[N[(a * b), $MachinePrecision], -1e+203], N[(a * N[(b / -4.0), $MachinePrecision]), $MachinePrecision], If[LessEqual[N[(a * b), $MachinePrecision], 1e+46], N[(N[(x * y), $MachinePrecision] + c), $MachinePrecision], N[(N[(a * b), $MachinePrecision] * -0.25), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;a \cdot b \leq -1 \cdot 10^{+203}:\\
\;\;\;\;a \cdot \frac{b}{-4}\\
\mathbf{elif}\;a \cdot b \leq 10^{+46}:\\
\;\;\;\;x \cdot y + c\\
\mathbf{else}:\\
\;\;\;\;\left(a \cdot b\right) \cdot -0.25\\
\end{array}
\end{array}
if (*.f64 a b) < -9.9999999999999999e202Initial program 88.0%
+-commutativeN/A
sub-negN/A
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f64N/A
distribute-neg-frac2N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
metadata-evalN/A
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f6488.0%
Simplified88.0%
Taylor expanded in a around inf
*-lowering-*.f64N/A
*-lowering-*.f6479.7%
Simplified79.7%
*-commutativeN/A
associate-*r*N/A
metadata-evalN/A
associate-/r/N/A
clear-numN/A
*-lowering-*.f64N/A
/-lowering-/.f6482.3%
Applied egg-rr82.3%
if -9.9999999999999999e202 < (*.f64 a b) < 9.9999999999999999e45Initial program 100.0%
+-commutativeN/A
sub-negN/A
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f64N/A
distribute-neg-frac2N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
metadata-evalN/A
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64100.0%
Simplified100.0%
Taylor expanded in a around 0
+-commutativeN/A
associate-+l+N/A
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f6495.7%
Simplified95.7%
Taylor expanded in t around 0
+-lowering-+.f64N/A
*-lowering-*.f6468.7%
Simplified68.7%
if 9.9999999999999999e45 < (*.f64 a b) Initial program 90.1%
+-commutativeN/A
sub-negN/A
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f64N/A
distribute-neg-frac2N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
metadata-evalN/A
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f6490.1%
Simplified90.1%
Taylor expanded in a around inf
*-lowering-*.f64N/A
*-lowering-*.f6469.3%
Simplified69.3%
Final simplification70.6%
(FPCore (x y z t a b c) :precision binary64 (if (<= (* x y) -1.4e+120) (* x y) (if (<= (* x y) 1.4e+37) c (* x y))))
double code(double x, double y, double z, double t, double a, double b, double c) {
double tmp;
if ((x * y) <= -1.4e+120) {
tmp = x * y;
} else if ((x * y) <= 1.4e+37) {
tmp = c;
} else {
tmp = x * y;
}
return tmp;
}
real(8) function code(x, y, z, t, a, b, c)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8), intent (in) :: c
real(8) :: tmp
if ((x * y) <= (-1.4d+120)) then
tmp = x * y
else if ((x * y) <= 1.4d+37) then
tmp = c
else
tmp = x * y
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b, double c) {
double tmp;
if ((x * y) <= -1.4e+120) {
tmp = x * y;
} else if ((x * y) <= 1.4e+37) {
tmp = c;
} else {
tmp = x * y;
}
return tmp;
}
def code(x, y, z, t, a, b, c): tmp = 0 if (x * y) <= -1.4e+120: tmp = x * y elif (x * y) <= 1.4e+37: tmp = c else: tmp = x * y return tmp
function code(x, y, z, t, a, b, c) tmp = 0.0 if (Float64(x * y) <= -1.4e+120) tmp = Float64(x * y); elseif (Float64(x * y) <= 1.4e+37) tmp = c; else tmp = Float64(x * y); end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c) tmp = 0.0; if ((x * y) <= -1.4e+120) tmp = x * y; elseif ((x * y) <= 1.4e+37) tmp = c; else tmp = x * y; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_] := If[LessEqual[N[(x * y), $MachinePrecision], -1.4e+120], N[(x * y), $MachinePrecision], If[LessEqual[N[(x * y), $MachinePrecision], 1.4e+37], c, N[(x * y), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \cdot y \leq -1.4 \cdot 10^{+120}:\\
\;\;\;\;x \cdot y\\
\mathbf{elif}\;x \cdot y \leq 1.4 \cdot 10^{+37}:\\
\;\;\;\;c\\
\mathbf{else}:\\
\;\;\;\;x \cdot y\\
\end{array}
\end{array}
if (*.f64 x y) < -1.4e120 or 1.3999999999999999e37 < (*.f64 x y) Initial program 94.2%
+-commutativeN/A
sub-negN/A
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f64N/A
distribute-neg-frac2N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
metadata-evalN/A
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f6494.2%
Simplified94.2%
Taylor expanded in x around inf
*-lowering-*.f6465.3%
Simplified65.3%
if -1.4e120 < (*.f64 x y) < 1.3999999999999999e37Initial program 97.7%
+-commutativeN/A
sub-negN/A
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f64N/A
distribute-neg-frac2N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
metadata-evalN/A
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f6497.7%
Simplified97.7%
Taylor expanded in c around inf
Simplified33.3%
(FPCore (x y z t a b c) :precision binary64 c)
double code(double x, double y, double z, double t, double a, double b, double c) {
return c;
}
real(8) function code(x, y, z, t, a, b, c)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8), intent (in) :: a
real(8), intent (in) :: b
real(8), intent (in) :: c
code = c
end function
public static double code(double x, double y, double z, double t, double a, double b, double c) {
return c;
}
def code(x, y, z, t, a, b, c): return c
function code(x, y, z, t, a, b, c) return c end
function tmp = code(x, y, z, t, a, b, c) tmp = c; end
code[x_, y_, z_, t_, a_, b_, c_] := c
\begin{array}{l}
\\
c
\end{array}
Initial program 96.5%
+-commutativeN/A
sub-negN/A
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f64N/A
distribute-neg-frac2N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
metadata-evalN/A
associate-+r+N/A
+-commutativeN/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f6496.5%
Simplified96.5%
Taylor expanded in c around inf
Simplified24.9%
herbie shell --seed 2024161
(FPCore (x y z t a b c)
:name "Diagrams.Solve.Polynomial:quartForm from diagrams-solve-0.1, C"
:precision binary64
(+ (- (+ (* x y) (/ (* z t) 16.0)) (/ (* a b) 4.0)) c))