
(FPCore (x y z t a b c)
:precision binary64
(/
x
(+
x
(*
y
(exp
(*
2.0
(-
(/ (* z (sqrt (+ t a))) t)
(* (- b c) (- (+ a (/ 5.0 6.0)) (/ 2.0 (* t 3.0)))))))))))
double code(double x, double y, double z, double t, double a, double b, double c) {
return x / (x + (y * exp((2.0 * (((z * sqrt((t + a))) / t) - ((b - c) * ((a + (5.0 / 6.0)) - (2.0 / (t * 3.0)))))))));
}
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 / (x + (y * exp((2.0d0 * (((z * sqrt((t + a))) / t) - ((b - c) * ((a + (5.0d0 / 6.0d0)) - (2.0d0 / (t * 3.0d0)))))))))
end function
public static double code(double x, double y, double z, double t, double a, double b, double c) {
return x / (x + (y * Math.exp((2.0 * (((z * Math.sqrt((t + a))) / t) - ((b - c) * ((a + (5.0 / 6.0)) - (2.0 / (t * 3.0)))))))));
}
def code(x, y, z, t, a, b, c): return x / (x + (y * math.exp((2.0 * (((z * math.sqrt((t + a))) / t) - ((b - c) * ((a + (5.0 / 6.0)) - (2.0 / (t * 3.0)))))))))
function code(x, y, z, t, a, b, c) return Float64(x / Float64(x + Float64(y * exp(Float64(2.0 * Float64(Float64(Float64(z * sqrt(Float64(t + a))) / t) - Float64(Float64(b - c) * Float64(Float64(a + Float64(5.0 / 6.0)) - Float64(2.0 / Float64(t * 3.0)))))))))) end
function tmp = code(x, y, z, t, a, b, c) tmp = x / (x + (y * exp((2.0 * (((z * sqrt((t + a))) / t) - ((b - c) * ((a + (5.0 / 6.0)) - (2.0 / (t * 3.0))))))))); end
code[x_, y_, z_, t_, a_, b_, c_] := N[(x / N[(x + N[(y * N[Exp[N[(2.0 * N[(N[(N[(z * N[Sqrt[N[(t + a), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] / t), $MachinePrecision] - N[(N[(b - c), $MachinePrecision] * N[(N[(a + N[(5.0 / 6.0), $MachinePrecision]), $MachinePrecision] - N[(2.0 / N[(t * 3.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{x}{x + y \cdot e^{2 \cdot \left(\frac{z \cdot \sqrt{t + a}}{t} - \left(b - c\right) \cdot \left(\left(a + \frac{5}{6}\right) - \frac{2}{t \cdot 3}\right)\right)}}
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 19 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (x y z t a b c)
:precision binary64
(/
x
(+
x
(*
y
(exp
(*
2.0
(-
(/ (* z (sqrt (+ t a))) t)
(* (- b c) (- (+ a (/ 5.0 6.0)) (/ 2.0 (* t 3.0)))))))))))
double code(double x, double y, double z, double t, double a, double b, double c) {
return x / (x + (y * exp((2.0 * (((z * sqrt((t + a))) / t) - ((b - c) * ((a + (5.0 / 6.0)) - (2.0 / (t * 3.0)))))))));
}
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 / (x + (y * exp((2.0d0 * (((z * sqrt((t + a))) / t) - ((b - c) * ((a + (5.0d0 / 6.0d0)) - (2.0d0 / (t * 3.0d0)))))))))
end function
public static double code(double x, double y, double z, double t, double a, double b, double c) {
return x / (x + (y * Math.exp((2.0 * (((z * Math.sqrt((t + a))) / t) - ((b - c) * ((a + (5.0 / 6.0)) - (2.0 / (t * 3.0)))))))));
}
def code(x, y, z, t, a, b, c): return x / (x + (y * math.exp((2.0 * (((z * math.sqrt((t + a))) / t) - ((b - c) * ((a + (5.0 / 6.0)) - (2.0 / (t * 3.0)))))))))
function code(x, y, z, t, a, b, c) return Float64(x / Float64(x + Float64(y * exp(Float64(2.0 * Float64(Float64(Float64(z * sqrt(Float64(t + a))) / t) - Float64(Float64(b - c) * Float64(Float64(a + Float64(5.0 / 6.0)) - Float64(2.0 / Float64(t * 3.0)))))))))) end
function tmp = code(x, y, z, t, a, b, c) tmp = x / (x + (y * exp((2.0 * (((z * sqrt((t + a))) / t) - ((b - c) * ((a + (5.0 / 6.0)) - (2.0 / (t * 3.0))))))))); end
code[x_, y_, z_, t_, a_, b_, c_] := N[(x / N[(x + N[(y * N[Exp[N[(2.0 * N[(N[(N[(z * N[Sqrt[N[(t + a), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] / t), $MachinePrecision] - N[(N[(b - c), $MachinePrecision] * N[(N[(a + N[(5.0 / 6.0), $MachinePrecision]), $MachinePrecision] - N[(2.0 / N[(t * 3.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{x}{x + y \cdot e^{2 \cdot \left(\frac{z \cdot \sqrt{t + a}}{t} - \left(b - c\right) \cdot \left(\left(a + \frac{5}{6}\right) - \frac{2}{t \cdot 3}\right)\right)}}
\end{array}
(FPCore (x y z t a b c)
:precision binary64
(let* ((t_1
(+
(/ (* z (sqrt (+ t a))) t)
(* (- b c) (- (/ 2.0 (* t 3.0)) (+ a (/ 5.0 6.0)))))))
(if (<= t_1 INFINITY)
(/ x (+ x (* y (exp (* 2.0 t_1)))))
(/
x
(+
x
(*
y
(exp
(*
2.0
(* c (+ (- a (/ 0.6666666666666666 t)) 0.8333333333333334))))))))))
double code(double x, double y, double z, double t, double a, double b, double c) {
double t_1 = ((z * sqrt((t + a))) / t) + ((b - c) * ((2.0 / (t * 3.0)) - (a + (5.0 / 6.0))));
double tmp;
if (t_1 <= ((double) INFINITY)) {
tmp = x / (x + (y * exp((2.0 * t_1))));
} else {
tmp = x / (x + (y * exp((2.0 * (c * ((a - (0.6666666666666666 / t)) + 0.8333333333333334))))));
}
return tmp;
}
public static double code(double x, double y, double z, double t, double a, double b, double c) {
double t_1 = ((z * Math.sqrt((t + a))) / t) + ((b - c) * ((2.0 / (t * 3.0)) - (a + (5.0 / 6.0))));
double tmp;
if (t_1 <= Double.POSITIVE_INFINITY) {
tmp = x / (x + (y * Math.exp((2.0 * t_1))));
} else {
tmp = x / (x + (y * Math.exp((2.0 * (c * ((a - (0.6666666666666666 / t)) + 0.8333333333333334))))));
}
return tmp;
}
def code(x, y, z, t, a, b, c): t_1 = ((z * math.sqrt((t + a))) / t) + ((b - c) * ((2.0 / (t * 3.0)) - (a + (5.0 / 6.0)))) tmp = 0 if t_1 <= math.inf: tmp = x / (x + (y * math.exp((2.0 * t_1)))) else: tmp = x / (x + (y * math.exp((2.0 * (c * ((a - (0.6666666666666666 / t)) + 0.8333333333333334)))))) return tmp
function code(x, y, z, t, a, b, c) t_1 = Float64(Float64(Float64(z * sqrt(Float64(t + a))) / t) + Float64(Float64(b - c) * Float64(Float64(2.0 / Float64(t * 3.0)) - Float64(a + Float64(5.0 / 6.0))))) tmp = 0.0 if (t_1 <= Inf) tmp = Float64(x / Float64(x + Float64(y * exp(Float64(2.0 * t_1))))); else tmp = Float64(x / Float64(x + Float64(y * exp(Float64(2.0 * Float64(c * Float64(Float64(a - Float64(0.6666666666666666 / t)) + 0.8333333333333334))))))); end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c) t_1 = ((z * sqrt((t + a))) / t) + ((b - c) * ((2.0 / (t * 3.0)) - (a + (5.0 / 6.0)))); tmp = 0.0; if (t_1 <= Inf) tmp = x / (x + (y * exp((2.0 * t_1)))); else tmp = x / (x + (y * exp((2.0 * (c * ((a - (0.6666666666666666 / t)) + 0.8333333333333334)))))); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_] := Block[{t$95$1 = N[(N[(N[(z * N[Sqrt[N[(t + a), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] / t), $MachinePrecision] + N[(N[(b - c), $MachinePrecision] * N[(N[(2.0 / N[(t * 3.0), $MachinePrecision]), $MachinePrecision] - N[(a + N[(5.0 / 6.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$1, Infinity], N[(x / N[(x + N[(y * N[Exp[N[(2.0 * t$95$1), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x / N[(x + N[(y * N[Exp[N[(2.0 * N[(c * N[(N[(a - N[(0.6666666666666666 / t), $MachinePrecision]), $MachinePrecision] + 0.8333333333333334), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{z \cdot \sqrt{t + a}}{t} + \left(b - c\right) \cdot \left(\frac{2}{t \cdot 3} - \left(a + \frac{5}{6}\right)\right)\\
\mathbf{if}\;t\_1 \leq \infty:\\
\;\;\;\;\frac{x}{x + y \cdot e^{2 \cdot t\_1}}\\
\mathbf{else}:\\
\;\;\;\;\frac{x}{x + y \cdot e^{2 \cdot \left(c \cdot \left(\left(a - \frac{0.6666666666666666}{t}\right) + 0.8333333333333334\right)\right)}}\\
\end{array}
\end{array}
if (-.f64 (/.f64 (*.f64 z (sqrt.f64 (+.f64 t a))) t) (*.f64 (-.f64 b c) (-.f64 (+.f64 a (/.f64 #s(literal 5 binary64) #s(literal 6 binary64))) (/.f64 #s(literal 2 binary64) (*.f64 t #s(literal 3 binary64)))))) < +inf.0Initial program 99.2%
if +inf.0 < (-.f64 (/.f64 (*.f64 z (sqrt.f64 (+.f64 t a))) t) (*.f64 (-.f64 b c) (-.f64 (+.f64 a (/.f64 #s(literal 5 binary64) #s(literal 6 binary64))) (/.f64 #s(literal 2 binary64) (*.f64 t #s(literal 3 binary64)))))) Initial program 0.0%
Taylor expanded in c around inf
*-lowering-*.f64N/A
associate--l+N/A
+-lowering-+.f64N/A
--lowering--.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f6474.3%
Simplified74.3%
Final simplification97.7%
(FPCore (x y z t a b c)
:precision binary64
(let* ((t_1
(/
x
(+
x
(*
y
(exp
(*
2.0
(*
c
(+ (- a (/ 0.6666666666666666 t)) 0.8333333333333334)))))))))
(if (<= c -1.4e+25)
t_1
(if (<= c -1.25e-298)
(/
x
(+
x
(*
y
(exp
(*
2.0
(* b (- (/ 0.6666666666666666 t) (+ a 0.8333333333333334))))))))
(if (<= c 620000000.0)
(/ x (+ x (* y (exp (* 2.0 (* (sqrt (+ t a)) (/ z t)))))))
t_1)))))
double code(double x, double y, double z, double t, double a, double b, double c) {
double t_1 = x / (x + (y * exp((2.0 * (c * ((a - (0.6666666666666666 / t)) + 0.8333333333333334))))));
double tmp;
if (c <= -1.4e+25) {
tmp = t_1;
} else if (c <= -1.25e-298) {
tmp = x / (x + (y * exp((2.0 * (b * ((0.6666666666666666 / t) - (a + 0.8333333333333334)))))));
} else if (c <= 620000000.0) {
tmp = x / (x + (y * exp((2.0 * (sqrt((t + a)) * (z / t))))));
} 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 = x / (x + (y * exp((2.0d0 * (c * ((a - (0.6666666666666666d0 / t)) + 0.8333333333333334d0))))))
if (c <= (-1.4d+25)) then
tmp = t_1
else if (c <= (-1.25d-298)) then
tmp = x / (x + (y * exp((2.0d0 * (b * ((0.6666666666666666d0 / t) - (a + 0.8333333333333334d0)))))))
else if (c <= 620000000.0d0) then
tmp = x / (x + (y * exp((2.0d0 * (sqrt((t + a)) * (z / t))))))
else
tmp = t_1
end if
code = tmp
end function
public static double code(double x, double y, double z, double t, double a, double b, double c) {
double t_1 = x / (x + (y * Math.exp((2.0 * (c * ((a - (0.6666666666666666 / t)) + 0.8333333333333334))))));
double tmp;
if (c <= -1.4e+25) {
tmp = t_1;
} else if (c <= -1.25e-298) {
tmp = x / (x + (y * Math.exp((2.0 * (b * ((0.6666666666666666 / t) - (a + 0.8333333333333334)))))));
} else if (c <= 620000000.0) {
tmp = x / (x + (y * Math.exp((2.0 * (Math.sqrt((t + a)) * (z / t))))));
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b, c): t_1 = x / (x + (y * math.exp((2.0 * (c * ((a - (0.6666666666666666 / t)) + 0.8333333333333334)))))) tmp = 0 if c <= -1.4e+25: tmp = t_1 elif c <= -1.25e-298: tmp = x / (x + (y * math.exp((2.0 * (b * ((0.6666666666666666 / t) - (a + 0.8333333333333334))))))) elif c <= 620000000.0: tmp = x / (x + (y * math.exp((2.0 * (math.sqrt((t + a)) * (z / t)))))) else: tmp = t_1 return tmp
function code(x, y, z, t, a, b, c) t_1 = Float64(x / Float64(x + Float64(y * exp(Float64(2.0 * Float64(c * Float64(Float64(a - Float64(0.6666666666666666 / t)) + 0.8333333333333334))))))) tmp = 0.0 if (c <= -1.4e+25) tmp = t_1; elseif (c <= -1.25e-298) tmp = Float64(x / Float64(x + Float64(y * exp(Float64(2.0 * Float64(b * Float64(Float64(0.6666666666666666 / t) - Float64(a + 0.8333333333333334)))))))); elseif (c <= 620000000.0) tmp = Float64(x / Float64(x + Float64(y * exp(Float64(2.0 * Float64(sqrt(Float64(t + a)) * Float64(z / t))))))); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c) t_1 = x / (x + (y * exp((2.0 * (c * ((a - (0.6666666666666666 / t)) + 0.8333333333333334)))))); tmp = 0.0; if (c <= -1.4e+25) tmp = t_1; elseif (c <= -1.25e-298) tmp = x / (x + (y * exp((2.0 * (b * ((0.6666666666666666 / t) - (a + 0.8333333333333334))))))); elseif (c <= 620000000.0) tmp = x / (x + (y * exp((2.0 * (sqrt((t + a)) * (z / t)))))); else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_] := Block[{t$95$1 = N[(x / N[(x + N[(y * N[Exp[N[(2.0 * N[(c * N[(N[(a - N[(0.6666666666666666 / t), $MachinePrecision]), $MachinePrecision] + 0.8333333333333334), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[c, -1.4e+25], t$95$1, If[LessEqual[c, -1.25e-298], N[(x / N[(x + N[(y * N[Exp[N[(2.0 * N[(b * N[(N[(0.6666666666666666 / t), $MachinePrecision] - N[(a + 0.8333333333333334), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[c, 620000000.0], N[(x / N[(x + N[(y * N[Exp[N[(2.0 * N[(N[Sqrt[N[(t + a), $MachinePrecision]], $MachinePrecision] * N[(z / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$1]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{x}{x + y \cdot e^{2 \cdot \left(c \cdot \left(\left(a - \frac{0.6666666666666666}{t}\right) + 0.8333333333333334\right)\right)}}\\
\mathbf{if}\;c \leq -1.4 \cdot 10^{+25}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;c \leq -1.25 \cdot 10^{-298}:\\
\;\;\;\;\frac{x}{x + y \cdot e^{2 \cdot \left(b \cdot \left(\frac{0.6666666666666666}{t} - \left(a + 0.8333333333333334\right)\right)\right)}}\\
\mathbf{elif}\;c \leq 620000000:\\
\;\;\;\;\frac{x}{x + y \cdot e^{2 \cdot \left(\sqrt{t + a} \cdot \frac{z}{t}\right)}}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if c < -1.4000000000000001e25 or 6.2e8 < c Initial program 93.8%
Taylor expanded in c around inf
*-lowering-*.f64N/A
associate--l+N/A
+-lowering-+.f64N/A
--lowering--.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f6494.8%
Simplified94.8%
if -1.4000000000000001e25 < c < -1.25000000000000005e-298Initial program 92.9%
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
count-2N/A
exp-lowering-exp.f64N/A
Simplified91.4%
Taylor expanded in b around inf
*-lowering-*.f64N/A
--lowering--.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f64N/A
+-lowering-+.f6479.8%
Simplified79.8%
if -1.25000000000000005e-298 < c < 6.2e8Initial program 93.0%
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
count-2N/A
exp-lowering-exp.f64N/A
Simplified96.4%
Taylor expanded in z around inf
*-lowering-*.f64N/A
/-lowering-/.f64N/A
sqrt-lowering-sqrt.f64N/A
+-commutativeN/A
+-lowering-+.f6481.1%
Simplified81.1%
Final simplification87.7%
(FPCore (x y z t a b c)
:precision binary64
(/
x
(+
x
(*
y
(exp
(*
2.0
(+
(* (sqrt (+ t a)) (/ z t))
(*
(- b c)
(+ (/ 0.6666666666666666 t) (- -0.8333333333333334 a))))))))))
double code(double x, double y, double z, double t, double a, double b, double c) {
return x / (x + (y * exp((2.0 * ((sqrt((t + a)) * (z / t)) + ((b - c) * ((0.6666666666666666 / t) + (-0.8333333333333334 - a))))))));
}
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 / (x + (y * exp((2.0d0 * ((sqrt((t + a)) * (z / t)) + ((b - c) * ((0.6666666666666666d0 / t) + ((-0.8333333333333334d0) - a))))))))
end function
public static double code(double x, double y, double z, double t, double a, double b, double c) {
return x / (x + (y * Math.exp((2.0 * ((Math.sqrt((t + a)) * (z / t)) + ((b - c) * ((0.6666666666666666 / t) + (-0.8333333333333334 - a))))))));
}
def code(x, y, z, t, a, b, c): return x / (x + (y * math.exp((2.0 * ((math.sqrt((t + a)) * (z / t)) + ((b - c) * ((0.6666666666666666 / t) + (-0.8333333333333334 - a))))))))
function code(x, y, z, t, a, b, c) return Float64(x / Float64(x + Float64(y * exp(Float64(2.0 * Float64(Float64(sqrt(Float64(t + a)) * Float64(z / t)) + Float64(Float64(b - c) * Float64(Float64(0.6666666666666666 / t) + Float64(-0.8333333333333334 - a))))))))) end
function tmp = code(x, y, z, t, a, b, c) tmp = x / (x + (y * exp((2.0 * ((sqrt((t + a)) * (z / t)) + ((b - c) * ((0.6666666666666666 / t) + (-0.8333333333333334 - a)))))))); end
code[x_, y_, z_, t_, a_, b_, c_] := N[(x / N[(x + N[(y * N[Exp[N[(2.0 * N[(N[(N[Sqrt[N[(t + a), $MachinePrecision]], $MachinePrecision] * N[(z / t), $MachinePrecision]), $MachinePrecision] + N[(N[(b - c), $MachinePrecision] * N[(N[(0.6666666666666666 / t), $MachinePrecision] + N[(-0.8333333333333334 - a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{x}{x + y \cdot e^{2 \cdot \left(\sqrt{t + a} \cdot \frac{z}{t} + \left(b - c\right) \cdot \left(\frac{0.6666666666666666}{t} + \left(-0.8333333333333334 - a\right)\right)\right)}}
\end{array}
Initial program 93.4%
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
count-2N/A
exp-lowering-exp.f64N/A
Simplified94.5%
(FPCore (x y z t a b c)
:precision binary64
(let* ((t_1
(/
x
(+
x
(*
y
(exp
(*
2.0
(*
c
(+ (- a (/ 0.6666666666666666 t)) 0.8333333333333334)))))))))
(if (<= c -1e+26)
t_1
(if (<= c 440000000.0)
(/
x
(+
x
(*
y
(exp
(*
2.0
(* b (- (/ 0.6666666666666666 t) (+ a 0.8333333333333334))))))))
t_1))))
double code(double x, double y, double z, double t, double a, double b, double c) {
double t_1 = x / (x + (y * exp((2.0 * (c * ((a - (0.6666666666666666 / t)) + 0.8333333333333334))))));
double tmp;
if (c <= -1e+26) {
tmp = t_1;
} else if (c <= 440000000.0) {
tmp = x / (x + (y * exp((2.0 * (b * ((0.6666666666666666 / t) - (a + 0.8333333333333334)))))));
} 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 = x / (x + (y * exp((2.0d0 * (c * ((a - (0.6666666666666666d0 / t)) + 0.8333333333333334d0))))))
if (c <= (-1d+26)) then
tmp = t_1
else if (c <= 440000000.0d0) then
tmp = x / (x + (y * exp((2.0d0 * (b * ((0.6666666666666666d0 / t) - (a + 0.8333333333333334d0)))))))
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 = x / (x + (y * Math.exp((2.0 * (c * ((a - (0.6666666666666666 / t)) + 0.8333333333333334))))));
double tmp;
if (c <= -1e+26) {
tmp = t_1;
} else if (c <= 440000000.0) {
tmp = x / (x + (y * Math.exp((2.0 * (b * ((0.6666666666666666 / t) - (a + 0.8333333333333334)))))));
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b, c): t_1 = x / (x + (y * math.exp((2.0 * (c * ((a - (0.6666666666666666 / t)) + 0.8333333333333334)))))) tmp = 0 if c <= -1e+26: tmp = t_1 elif c <= 440000000.0: tmp = x / (x + (y * math.exp((2.0 * (b * ((0.6666666666666666 / t) - (a + 0.8333333333333334))))))) else: tmp = t_1 return tmp
function code(x, y, z, t, a, b, c) t_1 = Float64(x / Float64(x + Float64(y * exp(Float64(2.0 * Float64(c * Float64(Float64(a - Float64(0.6666666666666666 / t)) + 0.8333333333333334))))))) tmp = 0.0 if (c <= -1e+26) tmp = t_1; elseif (c <= 440000000.0) tmp = Float64(x / Float64(x + Float64(y * exp(Float64(2.0 * Float64(b * Float64(Float64(0.6666666666666666 / t) - Float64(a + 0.8333333333333334)))))))); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c) t_1 = x / (x + (y * exp((2.0 * (c * ((a - (0.6666666666666666 / t)) + 0.8333333333333334)))))); tmp = 0.0; if (c <= -1e+26) tmp = t_1; elseif (c <= 440000000.0) tmp = x / (x + (y * exp((2.0 * (b * ((0.6666666666666666 / t) - (a + 0.8333333333333334))))))); else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_] := Block[{t$95$1 = N[(x / N[(x + N[(y * N[Exp[N[(2.0 * N[(c * N[(N[(a - N[(0.6666666666666666 / t), $MachinePrecision]), $MachinePrecision] + 0.8333333333333334), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[c, -1e+26], t$95$1, If[LessEqual[c, 440000000.0], N[(x / N[(x + N[(y * N[Exp[N[(2.0 * N[(b * N[(N[(0.6666666666666666 / t), $MachinePrecision] - N[(a + 0.8333333333333334), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{x}{x + y \cdot e^{2 \cdot \left(c \cdot \left(\left(a - \frac{0.6666666666666666}{t}\right) + 0.8333333333333334\right)\right)}}\\
\mathbf{if}\;c \leq -1 \cdot 10^{+26}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;c \leq 440000000:\\
\;\;\;\;\frac{x}{x + y \cdot e^{2 \cdot \left(b \cdot \left(\frac{0.6666666666666666}{t} - \left(a + 0.8333333333333334\right)\right)\right)}}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if c < -1.00000000000000005e26 or 4.4e8 < c Initial program 93.8%
Taylor expanded in c around inf
*-lowering-*.f64N/A
associate--l+N/A
+-lowering-+.f64N/A
--lowering--.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f6494.8%
Simplified94.8%
if -1.00000000000000005e26 < c < 4.4e8Initial program 92.9%
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
count-2N/A
exp-lowering-exp.f64N/A
Simplified93.7%
Taylor expanded in b around inf
*-lowering-*.f64N/A
--lowering--.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f64N/A
+-lowering-+.f6477.0%
Simplified77.0%
Final simplification86.0%
(FPCore (x y z t a b c)
:precision binary64
(let* ((t_1
(/
x
(+
x
(*
y
(exp
(*
2.0
(*
b
(- (/ 0.6666666666666666 t) (+ a 0.8333333333333334))))))))))
(if (<= b -3.3e+18)
t_1
(if (<= b 1.4e+26)
(/ x (+ x (* y (exp (* 2.0 (* c (+ a 0.8333333333333334)))))))
t_1))))
double code(double x, double y, double z, double t, double a, double b, double c) {
double t_1 = x / (x + (y * exp((2.0 * (b * ((0.6666666666666666 / t) - (a + 0.8333333333333334)))))));
double tmp;
if (b <= -3.3e+18) {
tmp = t_1;
} else if (b <= 1.4e+26) {
tmp = x / (x + (y * exp((2.0 * (c * (a + 0.8333333333333334))))));
} 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 = x / (x + (y * exp((2.0d0 * (b * ((0.6666666666666666d0 / t) - (a + 0.8333333333333334d0)))))))
if (b <= (-3.3d+18)) then
tmp = t_1
else if (b <= 1.4d+26) then
tmp = x / (x + (y * exp((2.0d0 * (c * (a + 0.8333333333333334d0))))))
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 = x / (x + (y * Math.exp((2.0 * (b * ((0.6666666666666666 / t) - (a + 0.8333333333333334)))))));
double tmp;
if (b <= -3.3e+18) {
tmp = t_1;
} else if (b <= 1.4e+26) {
tmp = x / (x + (y * Math.exp((2.0 * (c * (a + 0.8333333333333334))))));
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b, c): t_1 = x / (x + (y * math.exp((2.0 * (b * ((0.6666666666666666 / t) - (a + 0.8333333333333334))))))) tmp = 0 if b <= -3.3e+18: tmp = t_1 elif b <= 1.4e+26: tmp = x / (x + (y * math.exp((2.0 * (c * (a + 0.8333333333333334)))))) else: tmp = t_1 return tmp
function code(x, y, z, t, a, b, c) t_1 = Float64(x / Float64(x + Float64(y * exp(Float64(2.0 * Float64(b * Float64(Float64(0.6666666666666666 / t) - Float64(a + 0.8333333333333334)))))))) tmp = 0.0 if (b <= -3.3e+18) tmp = t_1; elseif (b <= 1.4e+26) tmp = Float64(x / Float64(x + Float64(y * exp(Float64(2.0 * Float64(c * Float64(a + 0.8333333333333334))))))); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c) t_1 = x / (x + (y * exp((2.0 * (b * ((0.6666666666666666 / t) - (a + 0.8333333333333334))))))); tmp = 0.0; if (b <= -3.3e+18) tmp = t_1; elseif (b <= 1.4e+26) tmp = x / (x + (y * exp((2.0 * (c * (a + 0.8333333333333334)))))); else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_] := Block[{t$95$1 = N[(x / N[(x + N[(y * N[Exp[N[(2.0 * N[(b * N[(N[(0.6666666666666666 / t), $MachinePrecision] - N[(a + 0.8333333333333334), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[b, -3.3e+18], t$95$1, If[LessEqual[b, 1.4e+26], N[(x / N[(x + N[(y * N[Exp[N[(2.0 * N[(c * N[(a + 0.8333333333333334), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{x}{x + y \cdot e^{2 \cdot \left(b \cdot \left(\frac{0.6666666666666666}{t} - \left(a + 0.8333333333333334\right)\right)\right)}}\\
\mathbf{if}\;b \leq -3.3 \cdot 10^{+18}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;b \leq 1.4 \cdot 10^{+26}:\\
\;\;\;\;\frac{x}{x + y \cdot e^{2 \cdot \left(c \cdot \left(a + 0.8333333333333334\right)\right)}}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if b < -3.3e18 or 1.4e26 < b Initial program 90.7%
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
count-2N/A
exp-lowering-exp.f64N/A
Simplified89.7%
Taylor expanded in b around inf
*-lowering-*.f64N/A
--lowering--.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f64N/A
+-lowering-+.f6489.1%
Simplified89.1%
if -3.3e18 < b < 1.4e26Initial program 95.3%
Taylor expanded in c around inf
*-lowering-*.f64N/A
associate--l+N/A
+-lowering-+.f64N/A
--lowering--.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f6481.0%
Simplified81.0%
Taylor expanded in t around inf
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
exp-lowering-exp.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f6468.7%
Simplified68.7%
Final simplification77.2%
(FPCore (x y z t a b c)
:precision binary64
(let* ((t_1 (/ x (+ x (* y (exp (* -2.0 (* b (+ a 0.8333333333333334)))))))))
(if (<= b -1e+42)
t_1
(if (<= b 1.25e+26)
(/ x (+ x (* y (exp (* 2.0 (* c (+ a 0.8333333333333334)))))))
t_1))))
double code(double x, double y, double z, double t, double a, double b, double c) {
double t_1 = x / (x + (y * exp((-2.0 * (b * (a + 0.8333333333333334))))));
double tmp;
if (b <= -1e+42) {
tmp = t_1;
} else if (b <= 1.25e+26) {
tmp = x / (x + (y * exp((2.0 * (c * (a + 0.8333333333333334))))));
} 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 = x / (x + (y * exp(((-2.0d0) * (b * (a + 0.8333333333333334d0))))))
if (b <= (-1d+42)) then
tmp = t_1
else if (b <= 1.25d+26) then
tmp = x / (x + (y * exp((2.0d0 * (c * (a + 0.8333333333333334d0))))))
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 = x / (x + (y * Math.exp((-2.0 * (b * (a + 0.8333333333333334))))));
double tmp;
if (b <= -1e+42) {
tmp = t_1;
} else if (b <= 1.25e+26) {
tmp = x / (x + (y * Math.exp((2.0 * (c * (a + 0.8333333333333334))))));
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b, c): t_1 = x / (x + (y * math.exp((-2.0 * (b * (a + 0.8333333333333334)))))) tmp = 0 if b <= -1e+42: tmp = t_1 elif b <= 1.25e+26: tmp = x / (x + (y * math.exp((2.0 * (c * (a + 0.8333333333333334)))))) else: tmp = t_1 return tmp
function code(x, y, z, t, a, b, c) t_1 = Float64(x / Float64(x + Float64(y * exp(Float64(-2.0 * Float64(b * Float64(a + 0.8333333333333334))))))) tmp = 0.0 if (b <= -1e+42) tmp = t_1; elseif (b <= 1.25e+26) tmp = Float64(x / Float64(x + Float64(y * exp(Float64(2.0 * Float64(c * Float64(a + 0.8333333333333334))))))); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c) t_1 = x / (x + (y * exp((-2.0 * (b * (a + 0.8333333333333334)))))); tmp = 0.0; if (b <= -1e+42) tmp = t_1; elseif (b <= 1.25e+26) tmp = x / (x + (y * exp((2.0 * (c * (a + 0.8333333333333334)))))); else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_] := Block[{t$95$1 = N[(x / N[(x + N[(y * N[Exp[N[(-2.0 * N[(b * N[(a + 0.8333333333333334), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[b, -1e+42], t$95$1, If[LessEqual[b, 1.25e+26], N[(x / N[(x + N[(y * N[Exp[N[(2.0 * N[(c * N[(a + 0.8333333333333334), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{x}{x + y \cdot e^{-2 \cdot \left(b \cdot \left(a + 0.8333333333333334\right)\right)}}\\
\mathbf{if}\;b \leq -1 \cdot 10^{+42}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;b \leq 1.25 \cdot 10^{+26}:\\
\;\;\;\;\frac{x}{x + y \cdot e^{2 \cdot \left(c \cdot \left(a + 0.8333333333333334\right)\right)}}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if b < -1.00000000000000004e42 or 1.25e26 < b Initial program 90.3%
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
count-2N/A
exp-lowering-exp.f64N/A
Simplified89.3%
Taylor expanded in b around inf
*-lowering-*.f64N/A
--lowering--.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f64N/A
+-lowering-+.f6489.7%
Simplified89.7%
Taylor expanded in t around inf
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
exp-lowering-exp.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f6475.5%
Simplified75.5%
if -1.00000000000000004e42 < b < 1.25e26Initial program 95.4%
Taylor expanded in c around inf
*-lowering-*.f64N/A
associate--l+N/A
+-lowering-+.f64N/A
--lowering--.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f6480.2%
Simplified80.2%
Taylor expanded in t around inf
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
exp-lowering-exp.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f6468.2%
Simplified68.2%
Final simplification71.2%
(FPCore (x y z t a b c)
:precision binary64
(let* ((t_1 (/ x (+ x (* y (exp (* -2.0 (* b (+ a 0.8333333333333334)))))))))
(if (<= b -1.45e+40)
t_1
(if (<= b 1.45e-25) (/ x (+ x (* y (exp (* 2.0 (* a c)))))) t_1))))
double code(double x, double y, double z, double t, double a, double b, double c) {
double t_1 = x / (x + (y * exp((-2.0 * (b * (a + 0.8333333333333334))))));
double tmp;
if (b <= -1.45e+40) {
tmp = t_1;
} else if (b <= 1.45e-25) {
tmp = x / (x + (y * exp((2.0 * (a * 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 = x / (x + (y * exp(((-2.0d0) * (b * (a + 0.8333333333333334d0))))))
if (b <= (-1.45d+40)) then
tmp = t_1
else if (b <= 1.45d-25) then
tmp = x / (x + (y * exp((2.0d0 * (a * 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 = x / (x + (y * Math.exp((-2.0 * (b * (a + 0.8333333333333334))))));
double tmp;
if (b <= -1.45e+40) {
tmp = t_1;
} else if (b <= 1.45e-25) {
tmp = x / (x + (y * Math.exp((2.0 * (a * c)))));
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t, a, b, c): t_1 = x / (x + (y * math.exp((-2.0 * (b * (a + 0.8333333333333334)))))) tmp = 0 if b <= -1.45e+40: tmp = t_1 elif b <= 1.45e-25: tmp = x / (x + (y * math.exp((2.0 * (a * c))))) else: tmp = t_1 return tmp
function code(x, y, z, t, a, b, c) t_1 = Float64(x / Float64(x + Float64(y * exp(Float64(-2.0 * Float64(b * Float64(a + 0.8333333333333334))))))) tmp = 0.0 if (b <= -1.45e+40) tmp = t_1; elseif (b <= 1.45e-25) tmp = Float64(x / Float64(x + Float64(y * exp(Float64(2.0 * Float64(a * c)))))); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c) t_1 = x / (x + (y * exp((-2.0 * (b * (a + 0.8333333333333334)))))); tmp = 0.0; if (b <= -1.45e+40) tmp = t_1; elseif (b <= 1.45e-25) tmp = x / (x + (y * exp((2.0 * (a * c))))); else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_] := Block[{t$95$1 = N[(x / N[(x + N[(y * N[Exp[N[(-2.0 * N[(b * N[(a + 0.8333333333333334), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[b, -1.45e+40], t$95$1, If[LessEqual[b, 1.45e-25], N[(x / N[(x + N[(y * N[Exp[N[(2.0 * N[(a * c), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$1]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{x}{x + y \cdot e^{-2 \cdot \left(b \cdot \left(a + 0.8333333333333334\right)\right)}}\\
\mathbf{if}\;b \leq -1.45 \cdot 10^{+40}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;b \leq 1.45 \cdot 10^{-25}:\\
\;\;\;\;\frac{x}{x + y \cdot e^{2 \cdot \left(a \cdot c\right)}}\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if b < -1.45000000000000009e40 or 1.45e-25 < b Initial program 91.4%
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
count-2N/A
exp-lowering-exp.f64N/A
Simplified90.5%
Taylor expanded in b around inf
*-lowering-*.f64N/A
--lowering--.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f64N/A
+-lowering-+.f6486.7%
Simplified86.7%
Taylor expanded in t around inf
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
exp-lowering-exp.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f6475.0%
Simplified75.0%
if -1.45000000000000009e40 < b < 1.45e-25Initial program 95.0%
Taylor expanded in c around inf
*-lowering-*.f64N/A
associate--l+N/A
+-lowering-+.f64N/A
--lowering--.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f6478.4%
Simplified78.4%
Taylor expanded in a around inf
Simplified63.2%
Final simplification68.6%
(FPCore (x y z t a b c)
:precision binary64
(let* ((t_1 (- (/ 0.6666666666666666 t) (+ a 0.8333333333333334))))
(if (<= b -1.2e+113)
(/ (* x 0.75) (* (* y (* b (* b b))) (* t_1 (* t_1 t_1))))
(if (<= b 2.5e+18) (/ x (+ x (* y (exp (* 2.0 (* a c)))))) 1.0))))
double code(double x, double y, double z, double t, double a, double b, double c) {
double t_1 = (0.6666666666666666 / t) - (a + 0.8333333333333334);
double tmp;
if (b <= -1.2e+113) {
tmp = (x * 0.75) / ((y * (b * (b * b))) * (t_1 * (t_1 * t_1)));
} else if (b <= 2.5e+18) {
tmp = x / (x + (y * exp((2.0 * (a * c)))));
} else {
tmp = 1.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.6666666666666666d0 / t) - (a + 0.8333333333333334d0)
if (b <= (-1.2d+113)) then
tmp = (x * 0.75d0) / ((y * (b * (b * b))) * (t_1 * (t_1 * t_1)))
else if (b <= 2.5d+18) then
tmp = x / (x + (y * exp((2.0d0 * (a * c)))))
else
tmp = 1.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.6666666666666666 / t) - (a + 0.8333333333333334);
double tmp;
if (b <= -1.2e+113) {
tmp = (x * 0.75) / ((y * (b * (b * b))) * (t_1 * (t_1 * t_1)));
} else if (b <= 2.5e+18) {
tmp = x / (x + (y * Math.exp((2.0 * (a * c)))));
} else {
tmp = 1.0;
}
return tmp;
}
def code(x, y, z, t, a, b, c): t_1 = (0.6666666666666666 / t) - (a + 0.8333333333333334) tmp = 0 if b <= -1.2e+113: tmp = (x * 0.75) / ((y * (b * (b * b))) * (t_1 * (t_1 * t_1))) elif b <= 2.5e+18: tmp = x / (x + (y * math.exp((2.0 * (a * c))))) else: tmp = 1.0 return tmp
function code(x, y, z, t, a, b, c) t_1 = Float64(Float64(0.6666666666666666 / t) - Float64(a + 0.8333333333333334)) tmp = 0.0 if (b <= -1.2e+113) tmp = Float64(Float64(x * 0.75) / Float64(Float64(y * Float64(b * Float64(b * b))) * Float64(t_1 * Float64(t_1 * t_1)))); elseif (b <= 2.5e+18) tmp = Float64(x / Float64(x + Float64(y * exp(Float64(2.0 * Float64(a * c)))))); else tmp = 1.0; end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c) t_1 = (0.6666666666666666 / t) - (a + 0.8333333333333334); tmp = 0.0; if (b <= -1.2e+113) tmp = (x * 0.75) / ((y * (b * (b * b))) * (t_1 * (t_1 * t_1))); elseif (b <= 2.5e+18) tmp = x / (x + (y * exp((2.0 * (a * c))))); else tmp = 1.0; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_] := Block[{t$95$1 = N[(N[(0.6666666666666666 / t), $MachinePrecision] - N[(a + 0.8333333333333334), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[b, -1.2e+113], N[(N[(x * 0.75), $MachinePrecision] / N[(N[(y * N[(b * N[(b * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(t$95$1 * N[(t$95$1 * t$95$1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[b, 2.5e+18], N[(x / N[(x + N[(y * N[Exp[N[(2.0 * N[(a * c), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 1.0]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{0.6666666666666666}{t} - \left(a + 0.8333333333333334\right)\\
\mathbf{if}\;b \leq -1.2 \cdot 10^{+113}:\\
\;\;\;\;\frac{x \cdot 0.75}{\left(y \cdot \left(b \cdot \left(b \cdot b\right)\right)\right) \cdot \left(t\_1 \cdot \left(t\_1 \cdot t\_1\right)\right)}\\
\mathbf{elif}\;b \leq 2.5 \cdot 10^{+18}:\\
\;\;\;\;\frac{x}{x + y \cdot e^{2 \cdot \left(a \cdot c\right)}}\\
\mathbf{else}:\\
\;\;\;\;1\\
\end{array}
\end{array}
if b < -1.19999999999999992e113Initial program 90.0%
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
count-2N/A
exp-lowering-exp.f64N/A
Simplified90.0%
Taylor expanded in b around inf
*-lowering-*.f64N/A
--lowering--.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f64N/A
+-lowering-+.f6495.2%
Simplified95.2%
Taylor expanded in b around 0
+-lowering-+.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
Simplified63.0%
Taylor expanded in b around inf
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
Simplified68.5%
if -1.19999999999999992e113 < b < 2.5e18Initial program 94.7%
Taylor expanded in c around inf
*-lowering-*.f64N/A
associate--l+N/A
+-lowering-+.f64N/A
--lowering--.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f6479.3%
Simplified79.3%
Taylor expanded in a around inf
Simplified63.9%
if 2.5e18 < b Initial program 91.5%
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
count-2N/A
exp-lowering-exp.f64N/A
Simplified87.2%
Taylor expanded in x around inf
Simplified67.0%
Final simplification65.2%
(FPCore (x y z t a b c)
:precision binary64
(let* ((t_1 (+ (- a (/ 0.6666666666666666 t)) 0.8333333333333334))
(t_2 (- (/ 0.6666666666666666 t) (+ a 0.8333333333333334)))
(t_3 (* t_2 t_2))
(t_4 (* t_1 t_1)))
(if (<= c -1e-209)
1.0
(if (<= c 31.0)
(/
x
(+
x
(*
y
(+
(*
b
(+
(* 2.0 t_2)
(* b (+ (* (* t_2 t_3) (* b 1.3333333333333333)) (* 2.0 t_3)))))
1.0))))
(if (<= c 1.05e+65)
1.0
(/
x
(+
x
(*
y
(+
(*
c
(+
(* 2.0 t_1)
(* c (+ (* (* c 1.3333333333333333) (* t_1 t_4)) (* 2.0 t_4)))))
1.0)))))))))
double code(double x, double y, double z, double t, double a, double b, double c) {
double t_1 = (a - (0.6666666666666666 / t)) + 0.8333333333333334;
double t_2 = (0.6666666666666666 / t) - (a + 0.8333333333333334);
double t_3 = t_2 * t_2;
double t_4 = t_1 * t_1;
double tmp;
if (c <= -1e-209) {
tmp = 1.0;
} else if (c <= 31.0) {
tmp = x / (x + (y * ((b * ((2.0 * t_2) + (b * (((t_2 * t_3) * (b * 1.3333333333333333)) + (2.0 * t_3))))) + 1.0)));
} else if (c <= 1.05e+65) {
tmp = 1.0;
} else {
tmp = x / (x + (y * ((c * ((2.0 * t_1) + (c * (((c * 1.3333333333333333) * (t_1 * t_4)) + (2.0 * t_4))))) + 1.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) :: t_2
real(8) :: t_3
real(8) :: t_4
real(8) :: tmp
t_1 = (a - (0.6666666666666666d0 / t)) + 0.8333333333333334d0
t_2 = (0.6666666666666666d0 / t) - (a + 0.8333333333333334d0)
t_3 = t_2 * t_2
t_4 = t_1 * t_1
if (c <= (-1d-209)) then
tmp = 1.0d0
else if (c <= 31.0d0) then
tmp = x / (x + (y * ((b * ((2.0d0 * t_2) + (b * (((t_2 * t_3) * (b * 1.3333333333333333d0)) + (2.0d0 * t_3))))) + 1.0d0)))
else if (c <= 1.05d+65) then
tmp = 1.0d0
else
tmp = x / (x + (y * ((c * ((2.0d0 * t_1) + (c * (((c * 1.3333333333333333d0) * (t_1 * t_4)) + (2.0d0 * t_4))))) + 1.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 = (a - (0.6666666666666666 / t)) + 0.8333333333333334;
double t_2 = (0.6666666666666666 / t) - (a + 0.8333333333333334);
double t_3 = t_2 * t_2;
double t_4 = t_1 * t_1;
double tmp;
if (c <= -1e-209) {
tmp = 1.0;
} else if (c <= 31.0) {
tmp = x / (x + (y * ((b * ((2.0 * t_2) + (b * (((t_2 * t_3) * (b * 1.3333333333333333)) + (2.0 * t_3))))) + 1.0)));
} else if (c <= 1.05e+65) {
tmp = 1.0;
} else {
tmp = x / (x + (y * ((c * ((2.0 * t_1) + (c * (((c * 1.3333333333333333) * (t_1 * t_4)) + (2.0 * t_4))))) + 1.0)));
}
return tmp;
}
def code(x, y, z, t, a, b, c): t_1 = (a - (0.6666666666666666 / t)) + 0.8333333333333334 t_2 = (0.6666666666666666 / t) - (a + 0.8333333333333334) t_3 = t_2 * t_2 t_4 = t_1 * t_1 tmp = 0 if c <= -1e-209: tmp = 1.0 elif c <= 31.0: tmp = x / (x + (y * ((b * ((2.0 * t_2) + (b * (((t_2 * t_3) * (b * 1.3333333333333333)) + (2.0 * t_3))))) + 1.0))) elif c <= 1.05e+65: tmp = 1.0 else: tmp = x / (x + (y * ((c * ((2.0 * t_1) + (c * (((c * 1.3333333333333333) * (t_1 * t_4)) + (2.0 * t_4))))) + 1.0))) return tmp
function code(x, y, z, t, a, b, c) t_1 = Float64(Float64(a - Float64(0.6666666666666666 / t)) + 0.8333333333333334) t_2 = Float64(Float64(0.6666666666666666 / t) - Float64(a + 0.8333333333333334)) t_3 = Float64(t_2 * t_2) t_4 = Float64(t_1 * t_1) tmp = 0.0 if (c <= -1e-209) tmp = 1.0; elseif (c <= 31.0) tmp = Float64(x / Float64(x + Float64(y * Float64(Float64(b * Float64(Float64(2.0 * t_2) + Float64(b * Float64(Float64(Float64(t_2 * t_3) * Float64(b * 1.3333333333333333)) + Float64(2.0 * t_3))))) + 1.0)))); elseif (c <= 1.05e+65) tmp = 1.0; else tmp = Float64(x / Float64(x + Float64(y * Float64(Float64(c * Float64(Float64(2.0 * t_1) + Float64(c * Float64(Float64(Float64(c * 1.3333333333333333) * Float64(t_1 * t_4)) + Float64(2.0 * t_4))))) + 1.0)))); end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c) t_1 = (a - (0.6666666666666666 / t)) + 0.8333333333333334; t_2 = (0.6666666666666666 / t) - (a + 0.8333333333333334); t_3 = t_2 * t_2; t_4 = t_1 * t_1; tmp = 0.0; if (c <= -1e-209) tmp = 1.0; elseif (c <= 31.0) tmp = x / (x + (y * ((b * ((2.0 * t_2) + (b * (((t_2 * t_3) * (b * 1.3333333333333333)) + (2.0 * t_3))))) + 1.0))); elseif (c <= 1.05e+65) tmp = 1.0; else tmp = x / (x + (y * ((c * ((2.0 * t_1) + (c * (((c * 1.3333333333333333) * (t_1 * t_4)) + (2.0 * t_4))))) + 1.0))); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_] := Block[{t$95$1 = N[(N[(a - N[(0.6666666666666666 / t), $MachinePrecision]), $MachinePrecision] + 0.8333333333333334), $MachinePrecision]}, Block[{t$95$2 = N[(N[(0.6666666666666666 / t), $MachinePrecision] - N[(a + 0.8333333333333334), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$3 = N[(t$95$2 * t$95$2), $MachinePrecision]}, Block[{t$95$4 = N[(t$95$1 * t$95$1), $MachinePrecision]}, If[LessEqual[c, -1e-209], 1.0, If[LessEqual[c, 31.0], N[(x / N[(x + N[(y * N[(N[(b * N[(N[(2.0 * t$95$2), $MachinePrecision] + N[(b * N[(N[(N[(t$95$2 * t$95$3), $MachinePrecision] * N[(b * 1.3333333333333333), $MachinePrecision]), $MachinePrecision] + N[(2.0 * t$95$3), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[c, 1.05e+65], 1.0, N[(x / N[(x + N[(y * N[(N[(c * N[(N[(2.0 * t$95$1), $MachinePrecision] + N[(c * N[(N[(N[(c * 1.3333333333333333), $MachinePrecision] * N[(t$95$1 * t$95$4), $MachinePrecision]), $MachinePrecision] + N[(2.0 * t$95$4), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \left(a - \frac{0.6666666666666666}{t}\right) + 0.8333333333333334\\
t_2 := \frac{0.6666666666666666}{t} - \left(a + 0.8333333333333334\right)\\
t_3 := t\_2 \cdot t\_2\\
t_4 := t\_1 \cdot t\_1\\
\mathbf{if}\;c \leq -1 \cdot 10^{-209}:\\
\;\;\;\;1\\
\mathbf{elif}\;c \leq 31:\\
\;\;\;\;\frac{x}{x + y \cdot \left(b \cdot \left(2 \cdot t\_2 + b \cdot \left(\left(t\_2 \cdot t\_3\right) \cdot \left(b \cdot 1.3333333333333333\right) + 2 \cdot t\_3\right)\right) + 1\right)}\\
\mathbf{elif}\;c \leq 1.05 \cdot 10^{+65}:\\
\;\;\;\;1\\
\mathbf{else}:\\
\;\;\;\;\frac{x}{x + y \cdot \left(c \cdot \left(2 \cdot t\_1 + c \cdot \left(\left(c \cdot 1.3333333333333333\right) \cdot \left(t\_1 \cdot t\_4\right) + 2 \cdot t\_4\right)\right) + 1\right)}\\
\end{array}
\end{array}
if c < -1e-209 or 31 < c < 1.04999999999999996e65Initial program 90.9%
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
count-2N/A
exp-lowering-exp.f64N/A
Simplified91.7%
Taylor expanded in x around inf
Simplified63.2%
if -1e-209 < c < 31Initial program 94.8%
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
count-2N/A
exp-lowering-exp.f64N/A
Simplified96.1%
Taylor expanded in b around inf
*-lowering-*.f64N/A
--lowering--.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f64N/A
+-lowering-+.f6478.4%
Simplified78.4%
Taylor expanded in b around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
Simplified61.4%
if 1.04999999999999996e65 < c Initial program 96.6%
Taylor expanded in c around inf
*-lowering-*.f64N/A
associate--l+N/A
+-lowering-+.f64N/A
--lowering--.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f6495.1%
Simplified95.1%
Taylor expanded in c around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
Simplified60.4%
Final simplification62.0%
(FPCore (x y z t a b c)
:precision binary64
(let* ((t_1 (- (/ 0.6666666666666666 t) (+ a 0.8333333333333334)))
(t_2 (* t_1 t_1)))
(if (<= c -3.3e-209)
1.0
(if (<= c 2.4e-272)
(/
x
(+
x
(*
y
(+
(*
b
(+
(* 2.0 t_1)
(* b (+ (* (* t_1 t_2) (* b 1.3333333333333333)) (* 2.0 t_2)))))
1.0))))
(if (<= c 6.5e-143)
(/ x (+ x (* y (* c (* a 2.0)))))
(if (<= c 1.12e+113)
(/ x (+ x (* y (+ (* b (* 2.0 (+ t_1 (* b t_2)))) 1.0))))
(/
x
(+
x
(*
y
(-
1.0
(/
(+
(* c 1.3333333333333333)
(* (/ (* c c) t) -0.8888888888888888))
t)))))))))))
double code(double x, double y, double z, double t, double a, double b, double c) {
double t_1 = (0.6666666666666666 / t) - (a + 0.8333333333333334);
double t_2 = t_1 * t_1;
double tmp;
if (c <= -3.3e-209) {
tmp = 1.0;
} else if (c <= 2.4e-272) {
tmp = x / (x + (y * ((b * ((2.0 * t_1) + (b * (((t_1 * t_2) * (b * 1.3333333333333333)) + (2.0 * t_2))))) + 1.0)));
} else if (c <= 6.5e-143) {
tmp = x / (x + (y * (c * (a * 2.0))));
} else if (c <= 1.12e+113) {
tmp = x / (x + (y * ((b * (2.0 * (t_1 + (b * t_2)))) + 1.0)));
} else {
tmp = x / (x + (y * (1.0 - (((c * 1.3333333333333333) + (((c * c) / t) * -0.8888888888888888)) / t))));
}
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 = (0.6666666666666666d0 / t) - (a + 0.8333333333333334d0)
t_2 = t_1 * t_1
if (c <= (-3.3d-209)) then
tmp = 1.0d0
else if (c <= 2.4d-272) then
tmp = x / (x + (y * ((b * ((2.0d0 * t_1) + (b * (((t_1 * t_2) * (b * 1.3333333333333333d0)) + (2.0d0 * t_2))))) + 1.0d0)))
else if (c <= 6.5d-143) then
tmp = x / (x + (y * (c * (a * 2.0d0))))
else if (c <= 1.12d+113) then
tmp = x / (x + (y * ((b * (2.0d0 * (t_1 + (b * t_2)))) + 1.0d0)))
else
tmp = x / (x + (y * (1.0d0 - (((c * 1.3333333333333333d0) + (((c * c) / t) * (-0.8888888888888888d0))) / t))))
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.6666666666666666 / t) - (a + 0.8333333333333334);
double t_2 = t_1 * t_1;
double tmp;
if (c <= -3.3e-209) {
tmp = 1.0;
} else if (c <= 2.4e-272) {
tmp = x / (x + (y * ((b * ((2.0 * t_1) + (b * (((t_1 * t_2) * (b * 1.3333333333333333)) + (2.0 * t_2))))) + 1.0)));
} else if (c <= 6.5e-143) {
tmp = x / (x + (y * (c * (a * 2.0))));
} else if (c <= 1.12e+113) {
tmp = x / (x + (y * ((b * (2.0 * (t_1 + (b * t_2)))) + 1.0)));
} else {
tmp = x / (x + (y * (1.0 - (((c * 1.3333333333333333) + (((c * c) / t) * -0.8888888888888888)) / t))));
}
return tmp;
}
def code(x, y, z, t, a, b, c): t_1 = (0.6666666666666666 / t) - (a + 0.8333333333333334) t_2 = t_1 * t_1 tmp = 0 if c <= -3.3e-209: tmp = 1.0 elif c <= 2.4e-272: tmp = x / (x + (y * ((b * ((2.0 * t_1) + (b * (((t_1 * t_2) * (b * 1.3333333333333333)) + (2.0 * t_2))))) + 1.0))) elif c <= 6.5e-143: tmp = x / (x + (y * (c * (a * 2.0)))) elif c <= 1.12e+113: tmp = x / (x + (y * ((b * (2.0 * (t_1 + (b * t_2)))) + 1.0))) else: tmp = x / (x + (y * (1.0 - (((c * 1.3333333333333333) + (((c * c) / t) * -0.8888888888888888)) / t)))) return tmp
function code(x, y, z, t, a, b, c) t_1 = Float64(Float64(0.6666666666666666 / t) - Float64(a + 0.8333333333333334)) t_2 = Float64(t_1 * t_1) tmp = 0.0 if (c <= -3.3e-209) tmp = 1.0; elseif (c <= 2.4e-272) tmp = Float64(x / Float64(x + Float64(y * Float64(Float64(b * Float64(Float64(2.0 * t_1) + Float64(b * Float64(Float64(Float64(t_1 * t_2) * Float64(b * 1.3333333333333333)) + Float64(2.0 * t_2))))) + 1.0)))); elseif (c <= 6.5e-143) tmp = Float64(x / Float64(x + Float64(y * Float64(c * Float64(a * 2.0))))); elseif (c <= 1.12e+113) tmp = Float64(x / Float64(x + Float64(y * Float64(Float64(b * Float64(2.0 * Float64(t_1 + Float64(b * t_2)))) + 1.0)))); else tmp = Float64(x / Float64(x + Float64(y * Float64(1.0 - Float64(Float64(Float64(c * 1.3333333333333333) + Float64(Float64(Float64(c * c) / t) * -0.8888888888888888)) / t))))); end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c) t_1 = (0.6666666666666666 / t) - (a + 0.8333333333333334); t_2 = t_1 * t_1; tmp = 0.0; if (c <= -3.3e-209) tmp = 1.0; elseif (c <= 2.4e-272) tmp = x / (x + (y * ((b * ((2.0 * t_1) + (b * (((t_1 * t_2) * (b * 1.3333333333333333)) + (2.0 * t_2))))) + 1.0))); elseif (c <= 6.5e-143) tmp = x / (x + (y * (c * (a * 2.0)))); elseif (c <= 1.12e+113) tmp = x / (x + (y * ((b * (2.0 * (t_1 + (b * t_2)))) + 1.0))); else tmp = x / (x + (y * (1.0 - (((c * 1.3333333333333333) + (((c * c) / t) * -0.8888888888888888)) / t)))); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_] := Block[{t$95$1 = N[(N[(0.6666666666666666 / t), $MachinePrecision] - N[(a + 0.8333333333333334), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(t$95$1 * t$95$1), $MachinePrecision]}, If[LessEqual[c, -3.3e-209], 1.0, If[LessEqual[c, 2.4e-272], N[(x / N[(x + N[(y * N[(N[(b * N[(N[(2.0 * t$95$1), $MachinePrecision] + N[(b * N[(N[(N[(t$95$1 * t$95$2), $MachinePrecision] * N[(b * 1.3333333333333333), $MachinePrecision]), $MachinePrecision] + N[(2.0 * t$95$2), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[c, 6.5e-143], N[(x / N[(x + N[(y * N[(c * N[(a * 2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[c, 1.12e+113], N[(x / N[(x + N[(y * N[(N[(b * N[(2.0 * N[(t$95$1 + N[(b * t$95$2), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x / N[(x + N[(y * N[(1.0 - N[(N[(N[(c * 1.3333333333333333), $MachinePrecision] + N[(N[(N[(c * c), $MachinePrecision] / t), $MachinePrecision] * -0.8888888888888888), $MachinePrecision]), $MachinePrecision] / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{0.6666666666666666}{t} - \left(a + 0.8333333333333334\right)\\
t_2 := t\_1 \cdot t\_1\\
\mathbf{if}\;c \leq -3.3 \cdot 10^{-209}:\\
\;\;\;\;1\\
\mathbf{elif}\;c \leq 2.4 \cdot 10^{-272}:\\
\;\;\;\;\frac{x}{x + y \cdot \left(b \cdot \left(2 \cdot t\_1 + b \cdot \left(\left(t\_1 \cdot t\_2\right) \cdot \left(b \cdot 1.3333333333333333\right) + 2 \cdot t\_2\right)\right) + 1\right)}\\
\mathbf{elif}\;c \leq 6.5 \cdot 10^{-143}:\\
\;\;\;\;\frac{x}{x + y \cdot \left(c \cdot \left(a \cdot 2\right)\right)}\\
\mathbf{elif}\;c \leq 1.12 \cdot 10^{+113}:\\
\;\;\;\;\frac{x}{x + y \cdot \left(b \cdot \left(2 \cdot \left(t\_1 + b \cdot t\_2\right)\right) + 1\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{x}{x + y \cdot \left(1 - \frac{c \cdot 1.3333333333333333 + \frac{c \cdot c}{t} \cdot -0.8888888888888888}{t}\right)}\\
\end{array}
\end{array}
if c < -3.29999999999999974e-209Initial program 90.9%
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
count-2N/A
exp-lowering-exp.f64N/A
Simplified91.8%
Taylor expanded in x around inf
Simplified62.1%
if -3.29999999999999974e-209 < c < 2.3999999999999999e-272Initial program 97.1%
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
count-2N/A
exp-lowering-exp.f64N/A
Simplified94.1%
Taylor expanded in b around inf
*-lowering-*.f64N/A
--lowering--.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f64N/A
+-lowering-+.f6476.8%
Simplified76.8%
Taylor expanded in b around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
Simplified71.1%
if 2.3999999999999999e-272 < c < 6.4999999999999999e-143Initial program 86.0%
Taylor expanded in c around inf
*-lowering-*.f64N/A
associate--l+N/A
+-lowering-+.f64N/A
--lowering--.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f6450.9%
Simplified50.9%
Taylor expanded in c around 0
+-lowering-+.f64N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
associate--l+N/A
+-lowering-+.f64N/A
--lowering--.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f6446.4%
Simplified46.4%
Taylor expanded in a around inf
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f6458.9%
Simplified58.9%
if 6.4999999999999999e-143 < c < 1.1200000000000001e113Initial program 95.5%
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
count-2N/A
exp-lowering-exp.f64N/A
Simplified97.7%
Taylor expanded in b around inf
*-lowering-*.f64N/A
--lowering--.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f64N/A
+-lowering-+.f6464.8%
Simplified64.8%
Taylor expanded in b around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
distribute-lft-outN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
Simplified60.6%
if 1.1200000000000001e113 < c Initial program 97.9%
Taylor expanded in c around inf
*-lowering-*.f64N/A
associate--l+N/A
+-lowering-+.f64N/A
--lowering--.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f6493.8%
Simplified93.8%
Taylor expanded in t around 0
associate-*r/N/A
/-lowering-/.f64N/A
*-commutativeN/A
*-lowering-*.f6462.9%
Simplified62.9%
Taylor expanded in t around -inf
+-lowering-+.f64N/A
mul-1-negN/A
neg-lowering-neg.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-lowering-*.f6465.0%
Simplified65.0%
Final simplification63.3%
(FPCore (x y z t a b c)
:precision binary64
(let* ((t_1 (- (/ 0.6666666666666666 t) (+ a 0.8333333333333334)))
(t_2 (+ (- a (/ 0.6666666666666666 t)) 0.8333333333333334)))
(if (<= b -8.5e-16)
(/ (* x 0.75) (* (* y (* b (* b b))) (* t_1 (* t_1 t_1))))
(if (<= b 1.8e-225)
1.0
(if (<= b 4.2e-152)
(/ x (+ x (* y (+ (* c (* 2.0 (+ t_2 (* c (* t_2 t_2))))) 1.0))))
1.0)))))
double code(double x, double y, double z, double t, double a, double b, double c) {
double t_1 = (0.6666666666666666 / t) - (a + 0.8333333333333334);
double t_2 = (a - (0.6666666666666666 / t)) + 0.8333333333333334;
double tmp;
if (b <= -8.5e-16) {
tmp = (x * 0.75) / ((y * (b * (b * b))) * (t_1 * (t_1 * t_1)));
} else if (b <= 1.8e-225) {
tmp = 1.0;
} else if (b <= 4.2e-152) {
tmp = x / (x + (y * ((c * (2.0 * (t_2 + (c * (t_2 * t_2))))) + 1.0)));
} else {
tmp = 1.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) :: t_2
real(8) :: tmp
t_1 = (0.6666666666666666d0 / t) - (a + 0.8333333333333334d0)
t_2 = (a - (0.6666666666666666d0 / t)) + 0.8333333333333334d0
if (b <= (-8.5d-16)) then
tmp = (x * 0.75d0) / ((y * (b * (b * b))) * (t_1 * (t_1 * t_1)))
else if (b <= 1.8d-225) then
tmp = 1.0d0
else if (b <= 4.2d-152) then
tmp = x / (x + (y * ((c * (2.0d0 * (t_2 + (c * (t_2 * t_2))))) + 1.0d0)))
else
tmp = 1.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.6666666666666666 / t) - (a + 0.8333333333333334);
double t_2 = (a - (0.6666666666666666 / t)) + 0.8333333333333334;
double tmp;
if (b <= -8.5e-16) {
tmp = (x * 0.75) / ((y * (b * (b * b))) * (t_1 * (t_1 * t_1)));
} else if (b <= 1.8e-225) {
tmp = 1.0;
} else if (b <= 4.2e-152) {
tmp = x / (x + (y * ((c * (2.0 * (t_2 + (c * (t_2 * t_2))))) + 1.0)));
} else {
tmp = 1.0;
}
return tmp;
}
def code(x, y, z, t, a, b, c): t_1 = (0.6666666666666666 / t) - (a + 0.8333333333333334) t_2 = (a - (0.6666666666666666 / t)) + 0.8333333333333334 tmp = 0 if b <= -8.5e-16: tmp = (x * 0.75) / ((y * (b * (b * b))) * (t_1 * (t_1 * t_1))) elif b <= 1.8e-225: tmp = 1.0 elif b <= 4.2e-152: tmp = x / (x + (y * ((c * (2.0 * (t_2 + (c * (t_2 * t_2))))) + 1.0))) else: tmp = 1.0 return tmp
function code(x, y, z, t, a, b, c) t_1 = Float64(Float64(0.6666666666666666 / t) - Float64(a + 0.8333333333333334)) t_2 = Float64(Float64(a - Float64(0.6666666666666666 / t)) + 0.8333333333333334) tmp = 0.0 if (b <= -8.5e-16) tmp = Float64(Float64(x * 0.75) / Float64(Float64(y * Float64(b * Float64(b * b))) * Float64(t_1 * Float64(t_1 * t_1)))); elseif (b <= 1.8e-225) tmp = 1.0; elseif (b <= 4.2e-152) tmp = Float64(x / Float64(x + Float64(y * Float64(Float64(c * Float64(2.0 * Float64(t_2 + Float64(c * Float64(t_2 * t_2))))) + 1.0)))); else tmp = 1.0; end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c) t_1 = (0.6666666666666666 / t) - (a + 0.8333333333333334); t_2 = (a - (0.6666666666666666 / t)) + 0.8333333333333334; tmp = 0.0; if (b <= -8.5e-16) tmp = (x * 0.75) / ((y * (b * (b * b))) * (t_1 * (t_1 * t_1))); elseif (b <= 1.8e-225) tmp = 1.0; elseif (b <= 4.2e-152) tmp = x / (x + (y * ((c * (2.0 * (t_2 + (c * (t_2 * t_2))))) + 1.0))); else tmp = 1.0; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_] := Block[{t$95$1 = N[(N[(0.6666666666666666 / t), $MachinePrecision] - N[(a + 0.8333333333333334), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(N[(a - N[(0.6666666666666666 / t), $MachinePrecision]), $MachinePrecision] + 0.8333333333333334), $MachinePrecision]}, If[LessEqual[b, -8.5e-16], N[(N[(x * 0.75), $MachinePrecision] / N[(N[(y * N[(b * N[(b * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(t$95$1 * N[(t$95$1 * t$95$1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[b, 1.8e-225], 1.0, If[LessEqual[b, 4.2e-152], N[(x / N[(x + N[(y * N[(N[(c * N[(2.0 * N[(t$95$2 + N[(c * N[(t$95$2 * t$95$2), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 1.0]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{0.6666666666666666}{t} - \left(a + 0.8333333333333334\right)\\
t_2 := \left(a - \frac{0.6666666666666666}{t}\right) + 0.8333333333333334\\
\mathbf{if}\;b \leq -8.5 \cdot 10^{-16}:\\
\;\;\;\;\frac{x \cdot 0.75}{\left(y \cdot \left(b \cdot \left(b \cdot b\right)\right)\right) \cdot \left(t\_1 \cdot \left(t\_1 \cdot t\_1\right)\right)}\\
\mathbf{elif}\;b \leq 1.8 \cdot 10^{-225}:\\
\;\;\;\;1\\
\mathbf{elif}\;b \leq 4.2 \cdot 10^{-152}:\\
\;\;\;\;\frac{x}{x + y \cdot \left(c \cdot \left(2 \cdot \left(t\_2 + c \cdot \left(t\_2 \cdot t\_2\right)\right)\right) + 1\right)}\\
\mathbf{else}:\\
\;\;\;\;1\\
\end{array}
\end{array}
if b < -8.5000000000000001e-16Initial program 91.2%
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
count-2N/A
exp-lowering-exp.f64N/A
Simplified92.6%
Taylor expanded in b around inf
*-lowering-*.f64N/A
--lowering--.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f64N/A
+-lowering-+.f6485.8%
Simplified85.8%
Taylor expanded in b around 0
+-lowering-+.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
Simplified55.1%
Taylor expanded in b around inf
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
Simplified63.0%
if -8.5000000000000001e-16 < b < 1.80000000000000005e-225 or 4.19999999999999998e-152 < b Initial program 93.8%
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
count-2N/A
exp-lowering-exp.f64N/A
Simplified94.9%
Taylor expanded in x around inf
Simplified58.9%
if 1.80000000000000005e-225 < b < 4.19999999999999998e-152Initial program 100.0%
Taylor expanded in c around inf
*-lowering-*.f64N/A
associate--l+N/A
+-lowering-+.f64N/A
--lowering--.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f6491.2%
Simplified91.2%
Taylor expanded in c around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
distribute-lft-outN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
Simplified91.3%
Final simplification61.4%
(FPCore (x y z t a b c)
:precision binary64
(let* ((t_1 (- (/ 0.6666666666666666 t) (+ a 0.8333333333333334))))
(if (<= b -1.75e-14)
(/ (* x 0.75) (* (* y (* b (* b b))) (* t_1 (* t_1 t_1))))
1.0)))
double code(double x, double y, double z, double t, double a, double b, double c) {
double t_1 = (0.6666666666666666 / t) - (a + 0.8333333333333334);
double tmp;
if (b <= -1.75e-14) {
tmp = (x * 0.75) / ((y * (b * (b * b))) * (t_1 * (t_1 * t_1)));
} else {
tmp = 1.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.6666666666666666d0 / t) - (a + 0.8333333333333334d0)
if (b <= (-1.75d-14)) then
tmp = (x * 0.75d0) / ((y * (b * (b * b))) * (t_1 * (t_1 * t_1)))
else
tmp = 1.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.6666666666666666 / t) - (a + 0.8333333333333334);
double tmp;
if (b <= -1.75e-14) {
tmp = (x * 0.75) / ((y * (b * (b * b))) * (t_1 * (t_1 * t_1)));
} else {
tmp = 1.0;
}
return tmp;
}
def code(x, y, z, t, a, b, c): t_1 = (0.6666666666666666 / t) - (a + 0.8333333333333334) tmp = 0 if b <= -1.75e-14: tmp = (x * 0.75) / ((y * (b * (b * b))) * (t_1 * (t_1 * t_1))) else: tmp = 1.0 return tmp
function code(x, y, z, t, a, b, c) t_1 = Float64(Float64(0.6666666666666666 / t) - Float64(a + 0.8333333333333334)) tmp = 0.0 if (b <= -1.75e-14) tmp = Float64(Float64(x * 0.75) / Float64(Float64(y * Float64(b * Float64(b * b))) * Float64(t_1 * Float64(t_1 * t_1)))); else tmp = 1.0; end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c) t_1 = (0.6666666666666666 / t) - (a + 0.8333333333333334); tmp = 0.0; if (b <= -1.75e-14) tmp = (x * 0.75) / ((y * (b * (b * b))) * (t_1 * (t_1 * t_1))); else tmp = 1.0; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_] := Block[{t$95$1 = N[(N[(0.6666666666666666 / t), $MachinePrecision] - N[(a + 0.8333333333333334), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[b, -1.75e-14], N[(N[(x * 0.75), $MachinePrecision] / N[(N[(y * N[(b * N[(b * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(t$95$1 * N[(t$95$1 * t$95$1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 1.0]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \frac{0.6666666666666666}{t} - \left(a + 0.8333333333333334\right)\\
\mathbf{if}\;b \leq -1.75 \cdot 10^{-14}:\\
\;\;\;\;\frac{x \cdot 0.75}{\left(y \cdot \left(b \cdot \left(b \cdot b\right)\right)\right) \cdot \left(t\_1 \cdot \left(t\_1 \cdot t\_1\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;1\\
\end{array}
\end{array}
if b < -1.7500000000000001e-14Initial program 91.2%
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
count-2N/A
exp-lowering-exp.f64N/A
Simplified92.6%
Taylor expanded in b around inf
*-lowering-*.f64N/A
--lowering--.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f64N/A
+-lowering-+.f6485.8%
Simplified85.8%
Taylor expanded in b around 0
+-lowering-+.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
Simplified55.1%
Taylor expanded in b around inf
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
cube-multN/A
unpow2N/A
*-lowering-*.f64N/A
Simplified63.0%
if -1.7500000000000001e-14 < b Initial program 94.2%
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
count-2N/A
exp-lowering-exp.f64N/A
Simplified95.2%
Taylor expanded in x around inf
Simplified56.2%
Final simplification58.0%
(FPCore (x y z t a b c)
:precision binary64
(if (<= c -1.75e-209)
1.0
(if (<= c 1.1e+113)
(/
x
(+
x
(*
y
(+
(*
2.0
(*
b
(-
(* b (* (+ a 0.8333333333333334) (+ a 0.8333333333333334)))
(+ a 0.8333333333333334))))
1.0))))
(/
x
(+
x
(*
y
(-
1.0
(/
(+ (* c 1.3333333333333333) (* (/ (* c c) t) -0.8888888888888888))
t))))))))
double code(double x, double y, double z, double t, double a, double b, double c) {
double tmp;
if (c <= -1.75e-209) {
tmp = 1.0;
} else if (c <= 1.1e+113) {
tmp = x / (x + (y * ((2.0 * (b * ((b * ((a + 0.8333333333333334) * (a + 0.8333333333333334))) - (a + 0.8333333333333334)))) + 1.0)));
} else {
tmp = x / (x + (y * (1.0 - (((c * 1.3333333333333333) + (((c * c) / t) * -0.8888888888888888)) / t))));
}
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 (c <= (-1.75d-209)) then
tmp = 1.0d0
else if (c <= 1.1d+113) then
tmp = x / (x + (y * ((2.0d0 * (b * ((b * ((a + 0.8333333333333334d0) * (a + 0.8333333333333334d0))) - (a + 0.8333333333333334d0)))) + 1.0d0)))
else
tmp = x / (x + (y * (1.0d0 - (((c * 1.3333333333333333d0) + (((c * c) / t) * (-0.8888888888888888d0))) / t))))
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 (c <= -1.75e-209) {
tmp = 1.0;
} else if (c <= 1.1e+113) {
tmp = x / (x + (y * ((2.0 * (b * ((b * ((a + 0.8333333333333334) * (a + 0.8333333333333334))) - (a + 0.8333333333333334)))) + 1.0)));
} else {
tmp = x / (x + (y * (1.0 - (((c * 1.3333333333333333) + (((c * c) / t) * -0.8888888888888888)) / t))));
}
return tmp;
}
def code(x, y, z, t, a, b, c): tmp = 0 if c <= -1.75e-209: tmp = 1.0 elif c <= 1.1e+113: tmp = x / (x + (y * ((2.0 * (b * ((b * ((a + 0.8333333333333334) * (a + 0.8333333333333334))) - (a + 0.8333333333333334)))) + 1.0))) else: tmp = x / (x + (y * (1.0 - (((c * 1.3333333333333333) + (((c * c) / t) * -0.8888888888888888)) / t)))) return tmp
function code(x, y, z, t, a, b, c) tmp = 0.0 if (c <= -1.75e-209) tmp = 1.0; elseif (c <= 1.1e+113) tmp = Float64(x / Float64(x + Float64(y * Float64(Float64(2.0 * Float64(b * Float64(Float64(b * Float64(Float64(a + 0.8333333333333334) * Float64(a + 0.8333333333333334))) - Float64(a + 0.8333333333333334)))) + 1.0)))); else tmp = Float64(x / Float64(x + Float64(y * Float64(1.0 - Float64(Float64(Float64(c * 1.3333333333333333) + Float64(Float64(Float64(c * c) / t) * -0.8888888888888888)) / t))))); end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c) tmp = 0.0; if (c <= -1.75e-209) tmp = 1.0; elseif (c <= 1.1e+113) tmp = x / (x + (y * ((2.0 * (b * ((b * ((a + 0.8333333333333334) * (a + 0.8333333333333334))) - (a + 0.8333333333333334)))) + 1.0))); else tmp = x / (x + (y * (1.0 - (((c * 1.3333333333333333) + (((c * c) / t) * -0.8888888888888888)) / t)))); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_] := If[LessEqual[c, -1.75e-209], 1.0, If[LessEqual[c, 1.1e+113], N[(x / N[(x + N[(y * N[(N[(2.0 * N[(b * N[(N[(b * N[(N[(a + 0.8333333333333334), $MachinePrecision] * N[(a + 0.8333333333333334), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - N[(a + 0.8333333333333334), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x / N[(x + N[(y * N[(1.0 - N[(N[(N[(c * 1.3333333333333333), $MachinePrecision] + N[(N[(N[(c * c), $MachinePrecision] / t), $MachinePrecision] * -0.8888888888888888), $MachinePrecision]), $MachinePrecision] / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;c \leq -1.75 \cdot 10^{-209}:\\
\;\;\;\;1\\
\mathbf{elif}\;c \leq 1.1 \cdot 10^{+113}:\\
\;\;\;\;\frac{x}{x + y \cdot \left(2 \cdot \left(b \cdot \left(b \cdot \left(\left(a + 0.8333333333333334\right) \cdot \left(a + 0.8333333333333334\right)\right) - \left(a + 0.8333333333333334\right)\right)\right) + 1\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{x}{x + y \cdot \left(1 - \frac{c \cdot 1.3333333333333333 + \frac{c \cdot c}{t} \cdot -0.8888888888888888}{t}\right)}\\
\end{array}
\end{array}
if c < -1.75000000000000001e-209Initial program 90.9%
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
count-2N/A
exp-lowering-exp.f64N/A
Simplified91.8%
Taylor expanded in x around inf
Simplified62.1%
if -1.75000000000000001e-209 < c < 1.10000000000000005e113Initial program 94.0%
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
count-2N/A
exp-lowering-exp.f64N/A
Simplified96.0%
Taylor expanded in b around inf
*-lowering-*.f64N/A
--lowering--.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f64N/A
+-lowering-+.f6469.7%
Simplified69.7%
Taylor expanded in b around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
distribute-lft-outN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
Simplified52.9%
Taylor expanded in t around inf
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
--lowering--.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
+-lowering-+.f6452.0%
Simplified52.0%
if 1.10000000000000005e113 < c Initial program 97.9%
Taylor expanded in c around inf
*-lowering-*.f64N/A
associate--l+N/A
+-lowering-+.f64N/A
--lowering--.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f6493.8%
Simplified93.8%
Taylor expanded in t around 0
associate-*r/N/A
/-lowering-/.f64N/A
*-commutativeN/A
*-lowering-*.f6462.9%
Simplified62.9%
Taylor expanded in t around -inf
+-lowering-+.f64N/A
mul-1-negN/A
neg-lowering-neg.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-lowering-*.f6465.0%
Simplified65.0%
Final simplification58.7%
(FPCore (x y z t a b c)
:precision binary64
(if (<= a 5e-305)
1.0
(if (<= a 2.2e-96)
(/
x
(+
x
(*
y
(+
(*
c
(- (* 0.8888888888888888 (/ c (* t t))) (/ 1.3333333333333333 t)))
1.0))))
(if (<= a 5e+240)
1.0
(/ x (- x (* y (- -1.0 (* c (* 2.0 (+ a (* c (* a a)))))))))))))
double code(double x, double y, double z, double t, double a, double b, double c) {
double tmp;
if (a <= 5e-305) {
tmp = 1.0;
} else if (a <= 2.2e-96) {
tmp = x / (x + (y * ((c * ((0.8888888888888888 * (c / (t * t))) - (1.3333333333333333 / t))) + 1.0)));
} else if (a <= 5e+240) {
tmp = 1.0;
} else {
tmp = x / (x - (y * (-1.0 - (c * (2.0 * (a + (c * (a * a))))))));
}
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 <= 5d-305) then
tmp = 1.0d0
else if (a <= 2.2d-96) then
tmp = x / (x + (y * ((c * ((0.8888888888888888d0 * (c / (t * t))) - (1.3333333333333333d0 / t))) + 1.0d0)))
else if (a <= 5d+240) then
tmp = 1.0d0
else
tmp = x / (x - (y * ((-1.0d0) - (c * (2.0d0 * (a + (c * (a * a))))))))
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 <= 5e-305) {
tmp = 1.0;
} else if (a <= 2.2e-96) {
tmp = x / (x + (y * ((c * ((0.8888888888888888 * (c / (t * t))) - (1.3333333333333333 / t))) + 1.0)));
} else if (a <= 5e+240) {
tmp = 1.0;
} else {
tmp = x / (x - (y * (-1.0 - (c * (2.0 * (a + (c * (a * a))))))));
}
return tmp;
}
def code(x, y, z, t, a, b, c): tmp = 0 if a <= 5e-305: tmp = 1.0 elif a <= 2.2e-96: tmp = x / (x + (y * ((c * ((0.8888888888888888 * (c / (t * t))) - (1.3333333333333333 / t))) + 1.0))) elif a <= 5e+240: tmp = 1.0 else: tmp = x / (x - (y * (-1.0 - (c * (2.0 * (a + (c * (a * a)))))))) return tmp
function code(x, y, z, t, a, b, c) tmp = 0.0 if (a <= 5e-305) tmp = 1.0; elseif (a <= 2.2e-96) tmp = Float64(x / Float64(x + Float64(y * Float64(Float64(c * Float64(Float64(0.8888888888888888 * Float64(c / Float64(t * t))) - Float64(1.3333333333333333 / t))) + 1.0)))); elseif (a <= 5e+240) tmp = 1.0; else tmp = Float64(x / Float64(x - Float64(y * Float64(-1.0 - Float64(c * Float64(2.0 * Float64(a + Float64(c * Float64(a * a))))))))); end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c) tmp = 0.0; if (a <= 5e-305) tmp = 1.0; elseif (a <= 2.2e-96) tmp = x / (x + (y * ((c * ((0.8888888888888888 * (c / (t * t))) - (1.3333333333333333 / t))) + 1.0))); elseif (a <= 5e+240) tmp = 1.0; else tmp = x / (x - (y * (-1.0 - (c * (2.0 * (a + (c * (a * a)))))))); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_] := If[LessEqual[a, 5e-305], 1.0, If[LessEqual[a, 2.2e-96], N[(x / N[(x + N[(y * N[(N[(c * N[(N[(0.8888888888888888 * N[(c / N[(t * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - N[(1.3333333333333333 / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[a, 5e+240], 1.0, N[(x / N[(x - N[(y * N[(-1.0 - N[(c * N[(2.0 * N[(a + N[(c * N[(a * a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;a \leq 5 \cdot 10^{-305}:\\
\;\;\;\;1\\
\mathbf{elif}\;a \leq 2.2 \cdot 10^{-96}:\\
\;\;\;\;\frac{x}{x + y \cdot \left(c \cdot \left(0.8888888888888888 \cdot \frac{c}{t \cdot t} - \frac{1.3333333333333333}{t}\right) + 1\right)}\\
\mathbf{elif}\;a \leq 5 \cdot 10^{+240}:\\
\;\;\;\;1\\
\mathbf{else}:\\
\;\;\;\;\frac{x}{x - y \cdot \left(-1 - c \cdot \left(2 \cdot \left(a + c \cdot \left(a \cdot a\right)\right)\right)\right)}\\
\end{array}
\end{array}
if a < 4.99999999999999985e-305 or 2.19999999999999979e-96 < a < 5.0000000000000003e240Initial program 96.0%
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
count-2N/A
exp-lowering-exp.f64N/A
Simplified97.1%
Taylor expanded in x around inf
Simplified59.1%
if 4.99999999999999985e-305 < a < 2.19999999999999979e-96Initial program 92.3%
Taylor expanded in c around inf
*-lowering-*.f64N/A
associate--l+N/A
+-lowering-+.f64N/A
--lowering--.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f6470.4%
Simplified70.4%
Taylor expanded in t around 0
associate-*r/N/A
/-lowering-/.f64N/A
*-commutativeN/A
*-lowering-*.f6453.7%
Simplified53.7%
Taylor expanded in c around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
--lowering--.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
unpow2N/A
*-lowering-*.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f6459.3%
Simplified59.3%
if 5.0000000000000003e240 < a Initial program 80.6%
Taylor expanded in c around inf
*-lowering-*.f64N/A
associate--l+N/A
+-lowering-+.f64N/A
--lowering--.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f6475.0%
Simplified75.0%
Taylor expanded in a around inf
Simplified75.0%
Taylor expanded in c around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
distribute-lft-outN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6462.5%
Simplified62.5%
Final simplification59.6%
(FPCore (x y z t a b c)
:precision binary64
(if (<= a 1.4e-306)
1.0
(if (<= a 7.6e-68)
(/ x (* c (- (/ (+ x y) c) (* 1.3333333333333333 (/ y t)))))
(if (<= a 3e+241)
1.0
(/ x (- x (* y (- -1.0 (* c (* 2.0 (+ a (* c (* a a)))))))))))))
double code(double x, double y, double z, double t, double a, double b, double c) {
double tmp;
if (a <= 1.4e-306) {
tmp = 1.0;
} else if (a <= 7.6e-68) {
tmp = x / (c * (((x + y) / c) - (1.3333333333333333 * (y / t))));
} else if (a <= 3e+241) {
tmp = 1.0;
} else {
tmp = x / (x - (y * (-1.0 - (c * (2.0 * (a + (c * (a * a))))))));
}
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 <= 1.4d-306) then
tmp = 1.0d0
else if (a <= 7.6d-68) then
tmp = x / (c * (((x + y) / c) - (1.3333333333333333d0 * (y / t))))
else if (a <= 3d+241) then
tmp = 1.0d0
else
tmp = x / (x - (y * ((-1.0d0) - (c * (2.0d0 * (a + (c * (a * a))))))))
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 <= 1.4e-306) {
tmp = 1.0;
} else if (a <= 7.6e-68) {
tmp = x / (c * (((x + y) / c) - (1.3333333333333333 * (y / t))));
} else if (a <= 3e+241) {
tmp = 1.0;
} else {
tmp = x / (x - (y * (-1.0 - (c * (2.0 * (a + (c * (a * a))))))));
}
return tmp;
}
def code(x, y, z, t, a, b, c): tmp = 0 if a <= 1.4e-306: tmp = 1.0 elif a <= 7.6e-68: tmp = x / (c * (((x + y) / c) - (1.3333333333333333 * (y / t)))) elif a <= 3e+241: tmp = 1.0 else: tmp = x / (x - (y * (-1.0 - (c * (2.0 * (a + (c * (a * a)))))))) return tmp
function code(x, y, z, t, a, b, c) tmp = 0.0 if (a <= 1.4e-306) tmp = 1.0; elseif (a <= 7.6e-68) tmp = Float64(x / Float64(c * Float64(Float64(Float64(x + y) / c) - Float64(1.3333333333333333 * Float64(y / t))))); elseif (a <= 3e+241) tmp = 1.0; else tmp = Float64(x / Float64(x - Float64(y * Float64(-1.0 - Float64(c * Float64(2.0 * Float64(a + Float64(c * Float64(a * a))))))))); end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c) tmp = 0.0; if (a <= 1.4e-306) tmp = 1.0; elseif (a <= 7.6e-68) tmp = x / (c * (((x + y) / c) - (1.3333333333333333 * (y / t)))); elseif (a <= 3e+241) tmp = 1.0; else tmp = x / (x - (y * (-1.0 - (c * (2.0 * (a + (c * (a * a)))))))); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_] := If[LessEqual[a, 1.4e-306], 1.0, If[LessEqual[a, 7.6e-68], N[(x / N[(c * N[(N[(N[(x + y), $MachinePrecision] / c), $MachinePrecision] - N[(1.3333333333333333 * N[(y / t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[a, 3e+241], 1.0, N[(x / N[(x - N[(y * N[(-1.0 - N[(c * N[(2.0 * N[(a + N[(c * N[(a * a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;a \leq 1.4 \cdot 10^{-306}:\\
\;\;\;\;1\\
\mathbf{elif}\;a \leq 7.6 \cdot 10^{-68}:\\
\;\;\;\;\frac{x}{c \cdot \left(\frac{x + y}{c} - 1.3333333333333333 \cdot \frac{y}{t}\right)}\\
\mathbf{elif}\;a \leq 3 \cdot 10^{+241}:\\
\;\;\;\;1\\
\mathbf{else}:\\
\;\;\;\;\frac{x}{x - y \cdot \left(-1 - c \cdot \left(2 \cdot \left(a + c \cdot \left(a \cdot a\right)\right)\right)\right)}\\
\end{array}
\end{array}
if a < 1.4000000000000001e-306 or 7.60000000000000075e-68 < a < 3.00000000000000015e241Initial program 95.6%
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
count-2N/A
exp-lowering-exp.f64N/A
Simplified96.9%
Taylor expanded in x around inf
Simplified61.0%
if 1.4000000000000001e-306 < a < 7.60000000000000075e-68Initial program 93.9%
Taylor expanded in c around inf
*-lowering-*.f64N/A
associate--l+N/A
+-lowering-+.f64N/A
--lowering--.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f6470.2%
Simplified70.2%
Taylor expanded in t around 0
associate-*r/N/A
/-lowering-/.f64N/A
*-commutativeN/A
*-lowering-*.f6449.7%
Simplified49.7%
Taylor expanded in c around 0
associate-+r+N/A
+-lowering-+.f64N/A
+-lowering-+.f64N/A
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f6441.1%
Simplified41.1%
Taylor expanded in c around -inf
mul-1-negN/A
neg-sub0N/A
--lowering--.f64N/A
*-lowering-*.f64N/A
+-commutativeN/A
mul-1-negN/A
unsub-negN/A
--lowering--.f64N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
+-lowering-+.f6447.4%
Simplified47.4%
if 3.00000000000000015e241 < a Initial program 80.6%
Taylor expanded in c around inf
*-lowering-*.f64N/A
associate--l+N/A
+-lowering-+.f64N/A
--lowering--.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f6475.0%
Simplified75.0%
Taylor expanded in a around inf
Simplified75.0%
Taylor expanded in c around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
distribute-lft-outN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6462.5%
Simplified62.5%
Final simplification57.7%
(FPCore (x y z t a b c)
:precision binary64
(if (<= a 3.2e-306)
1.0
(if (<= a 1.7e-97)
(/ x (+ x (* y (+ (/ (* c -1.3333333333333333) t) 1.0))))
1.0)))
double code(double x, double y, double z, double t, double a, double b, double c) {
double tmp;
if (a <= 3.2e-306) {
tmp = 1.0;
} else if (a <= 1.7e-97) {
tmp = x / (x + (y * (((c * -1.3333333333333333) / t) + 1.0)));
} else {
tmp = 1.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 <= 3.2d-306) then
tmp = 1.0d0
else if (a <= 1.7d-97) then
tmp = x / (x + (y * (((c * (-1.3333333333333333d0)) / t) + 1.0d0)))
else
tmp = 1.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 <= 3.2e-306) {
tmp = 1.0;
} else if (a <= 1.7e-97) {
tmp = x / (x + (y * (((c * -1.3333333333333333) / t) + 1.0)));
} else {
tmp = 1.0;
}
return tmp;
}
def code(x, y, z, t, a, b, c): tmp = 0 if a <= 3.2e-306: tmp = 1.0 elif a <= 1.7e-97: tmp = x / (x + (y * (((c * -1.3333333333333333) / t) + 1.0))) else: tmp = 1.0 return tmp
function code(x, y, z, t, a, b, c) tmp = 0.0 if (a <= 3.2e-306) tmp = 1.0; elseif (a <= 1.7e-97) tmp = Float64(x / Float64(x + Float64(y * Float64(Float64(Float64(c * -1.3333333333333333) / t) + 1.0)))); else tmp = 1.0; end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c) tmp = 0.0; if (a <= 3.2e-306) tmp = 1.0; elseif (a <= 1.7e-97) tmp = x / (x + (y * (((c * -1.3333333333333333) / t) + 1.0))); else tmp = 1.0; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_] := If[LessEqual[a, 3.2e-306], 1.0, If[LessEqual[a, 1.7e-97], N[(x / N[(x + N[(y * N[(N[(N[(c * -1.3333333333333333), $MachinePrecision] / t), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 1.0]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;a \leq 3.2 \cdot 10^{-306}:\\
\;\;\;\;1\\
\mathbf{elif}\;a \leq 1.7 \cdot 10^{-97}:\\
\;\;\;\;\frac{x}{x + y \cdot \left(\frac{c \cdot -1.3333333333333333}{t} + 1\right)}\\
\mathbf{else}:\\
\;\;\;\;1\\
\end{array}
\end{array}
if a < 3.19999999999999971e-306 or 1.6999999999999999e-97 < a Initial program 93.7%
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
count-2N/A
exp-lowering-exp.f64N/A
Simplified96.1%
Taylor expanded in x around inf
Simplified56.3%
if 3.19999999999999971e-306 < a < 1.6999999999999999e-97Initial program 92.3%
Taylor expanded in c around inf
*-lowering-*.f64N/A
associate--l+N/A
+-lowering-+.f64N/A
--lowering--.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f6470.4%
Simplified70.4%
Taylor expanded in t around 0
associate-*r/N/A
/-lowering-/.f64N/A
*-commutativeN/A
*-lowering-*.f6453.7%
Simplified53.7%
Taylor expanded in c around 0
+-lowering-+.f64N/A
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f6446.4%
Simplified46.4%
Final simplification54.3%
(FPCore (x y z t a b c) :precision binary64 (if (<= a 1.6e-205) 1.0 (if (<= a 1.45e-15) (* (/ x (- (* x x) (* y y))) (- x y)) 1.0)))
double code(double x, double y, double z, double t, double a, double b, double c) {
double tmp;
if (a <= 1.6e-205) {
tmp = 1.0;
} else if (a <= 1.45e-15) {
tmp = (x / ((x * x) - (y * y))) * (x - y);
} else {
tmp = 1.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 <= 1.6d-205) then
tmp = 1.0d0
else if (a <= 1.45d-15) then
tmp = (x / ((x * x) - (y * y))) * (x - y)
else
tmp = 1.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 <= 1.6e-205) {
tmp = 1.0;
} else if (a <= 1.45e-15) {
tmp = (x / ((x * x) - (y * y))) * (x - y);
} else {
tmp = 1.0;
}
return tmp;
}
def code(x, y, z, t, a, b, c): tmp = 0 if a <= 1.6e-205: tmp = 1.0 elif a <= 1.45e-15: tmp = (x / ((x * x) - (y * y))) * (x - y) else: tmp = 1.0 return tmp
function code(x, y, z, t, a, b, c) tmp = 0.0 if (a <= 1.6e-205) tmp = 1.0; elseif (a <= 1.45e-15) tmp = Float64(Float64(x / Float64(Float64(x * x) - Float64(y * y))) * Float64(x - y)); else tmp = 1.0; end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c) tmp = 0.0; if (a <= 1.6e-205) tmp = 1.0; elseif (a <= 1.45e-15) tmp = (x / ((x * x) - (y * y))) * (x - y); else tmp = 1.0; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_] := If[LessEqual[a, 1.6e-205], 1.0, If[LessEqual[a, 1.45e-15], N[(N[(x / N[(N[(x * x), $MachinePrecision] - N[(y * y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(x - y), $MachinePrecision]), $MachinePrecision], 1.0]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;a \leq 1.6 \cdot 10^{-205}:\\
\;\;\;\;1\\
\mathbf{elif}\;a \leq 1.45 \cdot 10^{-15}:\\
\;\;\;\;\frac{x}{x \cdot x - y \cdot y} \cdot \left(x - y\right)\\
\mathbf{else}:\\
\;\;\;\;1\\
\end{array}
\end{array}
if a < 1.60000000000000005e-205 or 1.45000000000000009e-15 < a Initial program 94.4%
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
count-2N/A
exp-lowering-exp.f64N/A
Simplified96.9%
Taylor expanded in x around inf
Simplified56.1%
if 1.60000000000000005e-205 < a < 1.45000000000000009e-15Initial program 90.3%
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
count-2N/A
exp-lowering-exp.f64N/A
Simplified87.1%
Taylor expanded in b around inf
*-lowering-*.f64N/A
--lowering--.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f64N/A
+-lowering-+.f6460.3%
Simplified60.3%
Taylor expanded in b around 0
/-lowering-/.f64N/A
+-lowering-+.f6431.6%
Simplified31.6%
flip-+N/A
associate-/r/N/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
--lowering--.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
--lowering--.f6447.8%
Applied egg-rr47.8%
(FPCore (x y z t a b c) :precision binary64 (if (<= b -3e+91) (/ x (+ x (* y (* (* b b) (* 2.0 (* a a)))))) 1.0))
double code(double x, double y, double z, double t, double a, double b, double c) {
double tmp;
if (b <= -3e+91) {
tmp = x / (x + (y * ((b * b) * (2.0 * (a * a)))));
} else {
tmp = 1.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 (b <= (-3d+91)) then
tmp = x / (x + (y * ((b * b) * (2.0d0 * (a * a)))))
else
tmp = 1.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 (b <= -3e+91) {
tmp = x / (x + (y * ((b * b) * (2.0 * (a * a)))));
} else {
tmp = 1.0;
}
return tmp;
}
def code(x, y, z, t, a, b, c): tmp = 0 if b <= -3e+91: tmp = x / (x + (y * ((b * b) * (2.0 * (a * a))))) else: tmp = 1.0 return tmp
function code(x, y, z, t, a, b, c) tmp = 0.0 if (b <= -3e+91) tmp = Float64(x / Float64(x + Float64(y * Float64(Float64(b * b) * Float64(2.0 * Float64(a * a)))))); else tmp = 1.0; end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c) tmp = 0.0; if (b <= -3e+91) tmp = x / (x + (y * ((b * b) * (2.0 * (a * a))))); else tmp = 1.0; end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_] := If[LessEqual[b, -3e+91], N[(x / N[(x + N[(y * N[(N[(b * b), $MachinePrecision] * N[(2.0 * N[(a * a), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], 1.0]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;b \leq -3 \cdot 10^{+91}:\\
\;\;\;\;\frac{x}{x + y \cdot \left(\left(b \cdot b\right) \cdot \left(2 \cdot \left(a \cdot a\right)\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;1\\
\end{array}
\end{array}
if b < -3.00000000000000006e91Initial program 91.5%
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
count-2N/A
exp-lowering-exp.f64N/A
Simplified91.5%
Taylor expanded in b around inf
*-lowering-*.f64N/A
--lowering--.f64N/A
associate-*r/N/A
metadata-evalN/A
/-lowering-/.f64N/A
+-lowering-+.f6495.9%
Simplified95.9%
Taylor expanded in b around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
distribute-lft-outN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
Simplified61.2%
Taylor expanded in a around inf
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6450.4%
Simplified50.4%
if -3.00000000000000006e91 < b Initial program 93.8%
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
count-2N/A
exp-lowering-exp.f64N/A
Simplified95.2%
Taylor expanded in x around inf
Simplified54.6%
Final simplification53.8%
(FPCore (x y z t a b c) :precision binary64 1.0)
double code(double x, double y, double z, double t, double a, double b, double c) {
return 1.0;
}
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 = 1.0d0
end function
public static double code(double x, double y, double z, double t, double a, double b, double c) {
return 1.0;
}
def code(x, y, z, t, a, b, c): return 1.0
function code(x, y, z, t, a, b, c) return 1.0 end
function tmp = code(x, y, z, t, a, b, c) tmp = 1.0; end
code[x_, y_, z_, t_, a_, b_, c_] := 1.0
\begin{array}{l}
\\
1
\end{array}
Initial program 93.4%
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
count-2N/A
exp-lowering-exp.f64N/A
Simplified94.5%
Taylor expanded in x around inf
Simplified50.8%
(FPCore (x y z t a b c)
:precision binary64
(let* ((t_1 (* z (sqrt (+ t a)))) (t_2 (- a (/ 5.0 6.0))))
(if (< t -2.118326644891581e-50)
(/
x
(+
x
(* y (exp (* 2.0 (- (+ (* a c) (* 0.8333333333333334 c)) (* a b)))))))
(if (< t 5.196588770651547e-123)
(/
x
(+
x
(*
y
(exp
(*
2.0
(/
(-
(* t_1 (* (* 3.0 t) t_2))
(*
(- (* (+ (/ 5.0 6.0) a) (* 3.0 t)) 2.0)
(* t_2 (* (- b c) t))))
(* (* (* t t) 3.0) t_2)))))))
(/
x
(+
x
(*
y
(exp
(*
2.0
(-
(/ t_1 t)
(* (- b c) (- (+ a (/ 5.0 6.0)) (/ 2.0 (* t 3.0))))))))))))))
double code(double x, double y, double z, double t, double a, double b, double c) {
double t_1 = z * sqrt((t + a));
double t_2 = a - (5.0 / 6.0);
double tmp;
if (t < -2.118326644891581e-50) {
tmp = x / (x + (y * exp((2.0 * (((a * c) + (0.8333333333333334 * c)) - (a * b))))));
} else if (t < 5.196588770651547e-123) {
tmp = x / (x + (y * exp((2.0 * (((t_1 * ((3.0 * t) * t_2)) - (((((5.0 / 6.0) + a) * (3.0 * t)) - 2.0) * (t_2 * ((b - c) * t)))) / (((t * t) * 3.0) * t_2))))));
} else {
tmp = x / (x + (y * exp((2.0 * ((t_1 / t) - ((b - c) * ((a + (5.0 / 6.0)) - (2.0 / (t * 3.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) :: t_2
real(8) :: tmp
t_1 = z * sqrt((t + a))
t_2 = a - (5.0d0 / 6.0d0)
if (t < (-2.118326644891581d-50)) then
tmp = x / (x + (y * exp((2.0d0 * (((a * c) + (0.8333333333333334d0 * c)) - (a * b))))))
else if (t < 5.196588770651547d-123) then
tmp = x / (x + (y * exp((2.0d0 * (((t_1 * ((3.0d0 * t) * t_2)) - (((((5.0d0 / 6.0d0) + a) * (3.0d0 * t)) - 2.0d0) * (t_2 * ((b - c) * t)))) / (((t * t) * 3.0d0) * t_2))))))
else
tmp = x / (x + (y * exp((2.0d0 * ((t_1 / t) - ((b - c) * ((a + (5.0d0 / 6.0d0)) - (2.0d0 / (t * 3.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 = z * Math.sqrt((t + a));
double t_2 = a - (5.0 / 6.0);
double tmp;
if (t < -2.118326644891581e-50) {
tmp = x / (x + (y * Math.exp((2.0 * (((a * c) + (0.8333333333333334 * c)) - (a * b))))));
} else if (t < 5.196588770651547e-123) {
tmp = x / (x + (y * Math.exp((2.0 * (((t_1 * ((3.0 * t) * t_2)) - (((((5.0 / 6.0) + a) * (3.0 * t)) - 2.0) * (t_2 * ((b - c) * t)))) / (((t * t) * 3.0) * t_2))))));
} else {
tmp = x / (x + (y * Math.exp((2.0 * ((t_1 / t) - ((b - c) * ((a + (5.0 / 6.0)) - (2.0 / (t * 3.0)))))))));
}
return tmp;
}
def code(x, y, z, t, a, b, c): t_1 = z * math.sqrt((t + a)) t_2 = a - (5.0 / 6.0) tmp = 0 if t < -2.118326644891581e-50: tmp = x / (x + (y * math.exp((2.0 * (((a * c) + (0.8333333333333334 * c)) - (a * b)))))) elif t < 5.196588770651547e-123: tmp = x / (x + (y * math.exp((2.0 * (((t_1 * ((3.0 * t) * t_2)) - (((((5.0 / 6.0) + a) * (3.0 * t)) - 2.0) * (t_2 * ((b - c) * t)))) / (((t * t) * 3.0) * t_2)))))) else: tmp = x / (x + (y * math.exp((2.0 * ((t_1 / t) - ((b - c) * ((a + (5.0 / 6.0)) - (2.0 / (t * 3.0))))))))) return tmp
function code(x, y, z, t, a, b, c) t_1 = Float64(z * sqrt(Float64(t + a))) t_2 = Float64(a - Float64(5.0 / 6.0)) tmp = 0.0 if (t < -2.118326644891581e-50) tmp = Float64(x / Float64(x + Float64(y * exp(Float64(2.0 * Float64(Float64(Float64(a * c) + Float64(0.8333333333333334 * c)) - Float64(a * b))))))); elseif (t < 5.196588770651547e-123) tmp = Float64(x / Float64(x + Float64(y * exp(Float64(2.0 * Float64(Float64(Float64(t_1 * Float64(Float64(3.0 * t) * t_2)) - Float64(Float64(Float64(Float64(Float64(5.0 / 6.0) + a) * Float64(3.0 * t)) - 2.0) * Float64(t_2 * Float64(Float64(b - c) * t)))) / Float64(Float64(Float64(t * t) * 3.0) * t_2))))))); else tmp = Float64(x / Float64(x + Float64(y * exp(Float64(2.0 * Float64(Float64(t_1 / t) - Float64(Float64(b - c) * Float64(Float64(a + Float64(5.0 / 6.0)) - Float64(2.0 / Float64(t * 3.0)))))))))); end return tmp end
function tmp_2 = code(x, y, z, t, a, b, c) t_1 = z * sqrt((t + a)); t_2 = a - (5.0 / 6.0); tmp = 0.0; if (t < -2.118326644891581e-50) tmp = x / (x + (y * exp((2.0 * (((a * c) + (0.8333333333333334 * c)) - (a * b)))))); elseif (t < 5.196588770651547e-123) tmp = x / (x + (y * exp((2.0 * (((t_1 * ((3.0 * t) * t_2)) - (((((5.0 / 6.0) + a) * (3.0 * t)) - 2.0) * (t_2 * ((b - c) * t)))) / (((t * t) * 3.0) * t_2)))))); else tmp = x / (x + (y * exp((2.0 * ((t_1 / t) - ((b - c) * ((a + (5.0 / 6.0)) - (2.0 / (t * 3.0))))))))); end tmp_2 = tmp; end
code[x_, y_, z_, t_, a_, b_, c_] := Block[{t$95$1 = N[(z * N[Sqrt[N[(t + a), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(a - N[(5.0 / 6.0), $MachinePrecision]), $MachinePrecision]}, If[Less[t, -2.118326644891581e-50], N[(x / N[(x + N[(y * N[Exp[N[(2.0 * N[(N[(N[(a * c), $MachinePrecision] + N[(0.8333333333333334 * c), $MachinePrecision]), $MachinePrecision] - N[(a * b), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[Less[t, 5.196588770651547e-123], N[(x / N[(x + N[(y * N[Exp[N[(2.0 * N[(N[(N[(t$95$1 * N[(N[(3.0 * t), $MachinePrecision] * t$95$2), $MachinePrecision]), $MachinePrecision] - N[(N[(N[(N[(N[(5.0 / 6.0), $MachinePrecision] + a), $MachinePrecision] * N[(3.0 * t), $MachinePrecision]), $MachinePrecision] - 2.0), $MachinePrecision] * N[(t$95$2 * N[(N[(b - c), $MachinePrecision] * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(N[(N[(t * t), $MachinePrecision] * 3.0), $MachinePrecision] * t$95$2), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x / N[(x + N[(y * N[Exp[N[(2.0 * N[(N[(t$95$1 / t), $MachinePrecision] - N[(N[(b - c), $MachinePrecision] * N[(N[(a + N[(5.0 / 6.0), $MachinePrecision]), $MachinePrecision] - N[(2.0 / N[(t * 3.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := z \cdot \sqrt{t + a}\\
t_2 := a - \frac{5}{6}\\
\mathbf{if}\;t < -2.118326644891581 \cdot 10^{-50}:\\
\;\;\;\;\frac{x}{x + y \cdot e^{2 \cdot \left(\left(a \cdot c + 0.8333333333333334 \cdot c\right) - a \cdot b\right)}}\\
\mathbf{elif}\;t < 5.196588770651547 \cdot 10^{-123}:\\
\;\;\;\;\frac{x}{x + y \cdot e^{2 \cdot \frac{t\_1 \cdot \left(\left(3 \cdot t\right) \cdot t\_2\right) - \left(\left(\frac{5}{6} + a\right) \cdot \left(3 \cdot t\right) - 2\right) \cdot \left(t\_2 \cdot \left(\left(b - c\right) \cdot t\right)\right)}{\left(\left(t \cdot t\right) \cdot 3\right) \cdot t\_2}}}\\
\mathbf{else}:\\
\;\;\;\;\frac{x}{x + y \cdot e^{2 \cdot \left(\frac{t\_1}{t} - \left(b - c\right) \cdot \left(\left(a + \frac{5}{6}\right) - \frac{2}{t \cdot 3}\right)\right)}}\\
\end{array}
\end{array}
herbie shell --seed 2024138
(FPCore (x y z t a b c)
:name "Numeric.SpecFunctions:invIncompleteBetaWorker from math-functions-0.1.5.2, I"
:precision binary64
:alt
(! :herbie-platform default (if (< t -2118326644891581/100000000000000000000000000000000000000000000000000000000000000000) (/ x (+ x (* y (exp (* 2 (- (+ (* a c) (* 4166666666666667/5000000000000000 c)) (* a b))))))) (if (< t 5196588770651547/1000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000) (/ x (+ x (* y (exp (* 2 (/ (- (* (* z (sqrt (+ t a))) (* (* 3 t) (- a (/ 5 6)))) (* (- (* (+ (/ 5 6) a) (* 3 t)) 2) (* (- a (/ 5 6)) (* (- b c) t)))) (* (* (* t t) 3) (- a (/ 5 6))))))))) (/ x (+ x (* y (exp (* 2 (- (/ (* z (sqrt (+ t a))) t) (* (- b c) (- (+ a (/ 5 6)) (/ 2 (* t 3)))))))))))))
(/ x (+ x (* y (exp (* 2.0 (- (/ (* z (sqrt (+ t a))) t) (* (- b c) (- (+ a (/ 5.0 6.0)) (/ 2.0 (* t 3.0)))))))))))