
(FPCore (x y z t) :precision binary64 (* (* (- (* x 0.5) y) (sqrt (* z 2.0))) (exp (/ (* t t) 2.0))))
double code(double x, double y, double z, double t) {
return (((x * 0.5) - y) * sqrt((z * 2.0))) * exp(((t * t) / 2.0));
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
code = (((x * 0.5d0) - y) * sqrt((z * 2.0d0))) * exp(((t * t) / 2.0d0))
end function
public static double code(double x, double y, double z, double t) {
return (((x * 0.5) - y) * Math.sqrt((z * 2.0))) * Math.exp(((t * t) / 2.0));
}
def code(x, y, z, t): return (((x * 0.5) - y) * math.sqrt((z * 2.0))) * math.exp(((t * t) / 2.0))
function code(x, y, z, t) return Float64(Float64(Float64(Float64(x * 0.5) - y) * sqrt(Float64(z * 2.0))) * exp(Float64(Float64(t * t) / 2.0))) end
function tmp = code(x, y, z, t) tmp = (((x * 0.5) - y) * sqrt((z * 2.0))) * exp(((t * t) / 2.0)); end
code[x_, y_, z_, t_] := N[(N[(N[(N[(x * 0.5), $MachinePrecision] - y), $MachinePrecision] * N[Sqrt[N[(z * 2.0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] * N[Exp[N[(N[(t * t), $MachinePrecision] / 2.0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\left(\left(x \cdot 0.5 - y\right) \cdot \sqrt{z \cdot 2}\right) \cdot e^{\frac{t \cdot t}{2}}
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 24 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (x y z t) :precision binary64 (* (* (- (* x 0.5) y) (sqrt (* z 2.0))) (exp (/ (* t t) 2.0))))
double code(double x, double y, double z, double t) {
return (((x * 0.5) - y) * sqrt((z * 2.0))) * exp(((t * t) / 2.0));
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
code = (((x * 0.5d0) - y) * sqrt((z * 2.0d0))) * exp(((t * t) / 2.0d0))
end function
public static double code(double x, double y, double z, double t) {
return (((x * 0.5) - y) * Math.sqrt((z * 2.0))) * Math.exp(((t * t) / 2.0));
}
def code(x, y, z, t): return (((x * 0.5) - y) * math.sqrt((z * 2.0))) * math.exp(((t * t) / 2.0))
function code(x, y, z, t) return Float64(Float64(Float64(Float64(x * 0.5) - y) * sqrt(Float64(z * 2.0))) * exp(Float64(Float64(t * t) / 2.0))) end
function tmp = code(x, y, z, t) tmp = (((x * 0.5) - y) * sqrt((z * 2.0))) * exp(((t * t) / 2.0)); end
code[x_, y_, z_, t_] := N[(N[(N[(N[(x * 0.5), $MachinePrecision] - y), $MachinePrecision] * N[Sqrt[N[(z * 2.0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] * N[Exp[N[(N[(t * t), $MachinePrecision] / 2.0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\left(\left(x \cdot 0.5 - y\right) \cdot \sqrt{z \cdot 2}\right) \cdot e^{\frac{t \cdot t}{2}}
\end{array}
(FPCore (x y z t) :precision binary64 (* (- (* x 0.5) y) (sqrt (* z (* 2.0 (exp (* t t)))))))
double code(double x, double y, double z, double t) {
return ((x * 0.5) - y) * sqrt((z * (2.0 * exp((t * t)))));
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
code = ((x * 0.5d0) - y) * sqrt((z * (2.0d0 * exp((t * t)))))
end function
public static double code(double x, double y, double z, double t) {
return ((x * 0.5) - y) * Math.sqrt((z * (2.0 * Math.exp((t * t)))));
}
def code(x, y, z, t): return ((x * 0.5) - y) * math.sqrt((z * (2.0 * math.exp((t * t)))))
function code(x, y, z, t) return Float64(Float64(Float64(x * 0.5) - y) * sqrt(Float64(z * Float64(2.0 * exp(Float64(t * t)))))) end
function tmp = code(x, y, z, t) tmp = ((x * 0.5) - y) * sqrt((z * (2.0 * exp((t * t))))); end
code[x_, y_, z_, t_] := N[(N[(N[(x * 0.5), $MachinePrecision] - y), $MachinePrecision] * N[Sqrt[N[(z * N[(2.0 * N[Exp[N[(t * t), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\left(x \cdot 0.5 - y\right) \cdot \sqrt{z \cdot \left(2 \cdot e^{t \cdot t}\right)}
\end{array}
Initial program 99.8%
associate-*l*N/A
*-lowering-*.f64N/A
--lowering--.f64N/A
*-lowering-*.f64N/A
exp-sqrtN/A
sqrt-unprodN/A
sqrt-lowering-sqrt.f64N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
exp-lowering-exp.f64N/A
*-lowering-*.f6499.8%
Applied egg-rr99.8%
(FPCore (x y z t)
:precision binary64
(let* ((t_1 (- (* x 0.5) y)) (t_2 (sqrt (* z 2.0))))
(if (<= (* t t) 2e-20)
(* t_1 t_2)
(if (<= (* t t) 1e+77)
(* (sqrt (* z (* 2.0 (exp (* t t))))) (* x 0.5))
(*
t_2
(*
t_1
(+
1.0
(*
t
(*
t
(+
0.5
(* (* t t) (+ 0.125 (* (* t t) 0.020833333333333332)))))))))))))
double code(double x, double y, double z, double t) {
double t_1 = (x * 0.5) - y;
double t_2 = sqrt((z * 2.0));
double tmp;
if ((t * t) <= 2e-20) {
tmp = t_1 * t_2;
} else if ((t * t) <= 1e+77) {
tmp = sqrt((z * (2.0 * exp((t * t))))) * (x * 0.5);
} else {
tmp = t_2 * (t_1 * (1.0 + (t * (t * (0.5 + ((t * t) * (0.125 + ((t * t) * 0.020833333333333332))))))));
}
return tmp;
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8) :: t_1
real(8) :: t_2
real(8) :: tmp
t_1 = (x * 0.5d0) - y
t_2 = sqrt((z * 2.0d0))
if ((t * t) <= 2d-20) then
tmp = t_1 * t_2
else if ((t * t) <= 1d+77) then
tmp = sqrt((z * (2.0d0 * exp((t * t))))) * (x * 0.5d0)
else
tmp = t_2 * (t_1 * (1.0d0 + (t * (t * (0.5d0 + ((t * t) * (0.125d0 + ((t * t) * 0.020833333333333332d0))))))))
end if
code = tmp
end function
public static double code(double x, double y, double z, double t) {
double t_1 = (x * 0.5) - y;
double t_2 = Math.sqrt((z * 2.0));
double tmp;
if ((t * t) <= 2e-20) {
tmp = t_1 * t_2;
} else if ((t * t) <= 1e+77) {
tmp = Math.sqrt((z * (2.0 * Math.exp((t * t))))) * (x * 0.5);
} else {
tmp = t_2 * (t_1 * (1.0 + (t * (t * (0.5 + ((t * t) * (0.125 + ((t * t) * 0.020833333333333332))))))));
}
return tmp;
}
def code(x, y, z, t): t_1 = (x * 0.5) - y t_2 = math.sqrt((z * 2.0)) tmp = 0 if (t * t) <= 2e-20: tmp = t_1 * t_2 elif (t * t) <= 1e+77: tmp = math.sqrt((z * (2.0 * math.exp((t * t))))) * (x * 0.5) else: tmp = t_2 * (t_1 * (1.0 + (t * (t * (0.5 + ((t * t) * (0.125 + ((t * t) * 0.020833333333333332)))))))) return tmp
function code(x, y, z, t) t_1 = Float64(Float64(x * 0.5) - y) t_2 = sqrt(Float64(z * 2.0)) tmp = 0.0 if (Float64(t * t) <= 2e-20) tmp = Float64(t_1 * t_2); elseif (Float64(t * t) <= 1e+77) tmp = Float64(sqrt(Float64(z * Float64(2.0 * exp(Float64(t * t))))) * Float64(x * 0.5)); else tmp = Float64(t_2 * Float64(t_1 * Float64(1.0 + Float64(t * Float64(t * Float64(0.5 + Float64(Float64(t * t) * Float64(0.125 + Float64(Float64(t * t) * 0.020833333333333332))))))))); end return tmp end
function tmp_2 = code(x, y, z, t) t_1 = (x * 0.5) - y; t_2 = sqrt((z * 2.0)); tmp = 0.0; if ((t * t) <= 2e-20) tmp = t_1 * t_2; elseif ((t * t) <= 1e+77) tmp = sqrt((z * (2.0 * exp((t * t))))) * (x * 0.5); else tmp = t_2 * (t_1 * (1.0 + (t * (t * (0.5 + ((t * t) * (0.125 + ((t * t) * 0.020833333333333332)))))))); end tmp_2 = tmp; end
code[x_, y_, z_, t_] := Block[{t$95$1 = N[(N[(x * 0.5), $MachinePrecision] - y), $MachinePrecision]}, Block[{t$95$2 = N[Sqrt[N[(z * 2.0), $MachinePrecision]], $MachinePrecision]}, If[LessEqual[N[(t * t), $MachinePrecision], 2e-20], N[(t$95$1 * t$95$2), $MachinePrecision], If[LessEqual[N[(t * t), $MachinePrecision], 1e+77], N[(N[Sqrt[N[(z * N[(2.0 * N[Exp[N[(t * t), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * N[(x * 0.5), $MachinePrecision]), $MachinePrecision], N[(t$95$2 * N[(t$95$1 * N[(1.0 + N[(t * N[(t * N[(0.5 + N[(N[(t * t), $MachinePrecision] * N[(0.125 + N[(N[(t * t), $MachinePrecision] * 0.020833333333333332), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x \cdot 0.5 - y\\
t_2 := \sqrt{z \cdot 2}\\
\mathbf{if}\;t \cdot t \leq 2 \cdot 10^{-20}:\\
\;\;\;\;t\_1 \cdot t\_2\\
\mathbf{elif}\;t \cdot t \leq 10^{+77}:\\
\;\;\;\;\sqrt{z \cdot \left(2 \cdot e^{t \cdot t}\right)} \cdot \left(x \cdot 0.5\right)\\
\mathbf{else}:\\
\;\;\;\;t\_2 \cdot \left(t\_1 \cdot \left(1 + t \cdot \left(t \cdot \left(0.5 + \left(t \cdot t\right) \cdot \left(0.125 + \left(t \cdot t\right) \cdot 0.020833333333333332\right)\right)\right)\right)\right)\\
\end{array}
\end{array}
if (*.f64 t t) < 1.99999999999999989e-20Initial program 99.6%
Taylor expanded in t around 0
Simplified99.6%
*-rgt-identityN/A
*-lowering-*.f64N/A
--lowering--.f64N/A
*-lowering-*.f64N/A
sqrt-lowering-sqrt.f64N/A
*-lowering-*.f6499.6%
Applied egg-rr99.6%
if 1.99999999999999989e-20 < (*.f64 t t) < 9.99999999999999983e76Initial program 99.8%
associate-*l*N/A
*-lowering-*.f64N/A
--lowering--.f64N/A
*-lowering-*.f64N/A
exp-sqrtN/A
sqrt-unprodN/A
sqrt-lowering-sqrt.f64N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
exp-lowering-exp.f64N/A
*-lowering-*.f64100.0%
Applied egg-rr100.0%
Taylor expanded in x around inf
*-lowering-*.f6494.9%
Simplified94.9%
if 9.99999999999999983e76 < (*.f64 t t) Initial program 100.0%
Taylor expanded in t around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f6499.1%
Simplified99.1%
*-commutativeN/A
associate-*r*N/A
*-lowering-*.f64N/A
Applied egg-rr99.1%
Final simplification99.1%
(FPCore (x y z t)
:precision binary64
(let* ((t_1
(* t (+ 0.5 (* (* t t) (+ 0.125 (* (* t t) 0.020833333333333332))))))
(t_2 (- (* x 0.5) y))
(t_3 (* t t_1))
(t_4 (+ y (* x 0.5))))
(if (<= (* t t) 2e+47)
(/
(* (* t_2 (sqrt (* z 2.0))) (+ 1.0 (* t_3 (* (* t t) (* t_1 t_1)))))
(+ 1.0 (* t_3 (+ t_3 -1.0))))
(* t_2 (* t_4 (/ (/ (+ 1.0 t_3) (pow (* z 2.0) -0.5)) t_4))))))
double code(double x, double y, double z, double t) {
double t_1 = t * (0.5 + ((t * t) * (0.125 + ((t * t) * 0.020833333333333332))));
double t_2 = (x * 0.5) - y;
double t_3 = t * t_1;
double t_4 = y + (x * 0.5);
double tmp;
if ((t * t) <= 2e+47) {
tmp = ((t_2 * sqrt((z * 2.0))) * (1.0 + (t_3 * ((t * t) * (t_1 * t_1))))) / (1.0 + (t_3 * (t_3 + -1.0)));
} else {
tmp = t_2 * (t_4 * (((1.0 + t_3) / pow((z * 2.0), -0.5)) / t_4));
}
return tmp;
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8) :: t_1
real(8) :: t_2
real(8) :: t_3
real(8) :: t_4
real(8) :: tmp
t_1 = t * (0.5d0 + ((t * t) * (0.125d0 + ((t * t) * 0.020833333333333332d0))))
t_2 = (x * 0.5d0) - y
t_3 = t * t_1
t_4 = y + (x * 0.5d0)
if ((t * t) <= 2d+47) then
tmp = ((t_2 * sqrt((z * 2.0d0))) * (1.0d0 + (t_3 * ((t * t) * (t_1 * t_1))))) / (1.0d0 + (t_3 * (t_3 + (-1.0d0))))
else
tmp = t_2 * (t_4 * (((1.0d0 + t_3) / ((z * 2.0d0) ** (-0.5d0))) / t_4))
end if
code = tmp
end function
public static double code(double x, double y, double z, double t) {
double t_1 = t * (0.5 + ((t * t) * (0.125 + ((t * t) * 0.020833333333333332))));
double t_2 = (x * 0.5) - y;
double t_3 = t * t_1;
double t_4 = y + (x * 0.5);
double tmp;
if ((t * t) <= 2e+47) {
tmp = ((t_2 * Math.sqrt((z * 2.0))) * (1.0 + (t_3 * ((t * t) * (t_1 * t_1))))) / (1.0 + (t_3 * (t_3 + -1.0)));
} else {
tmp = t_2 * (t_4 * (((1.0 + t_3) / Math.pow((z * 2.0), -0.5)) / t_4));
}
return tmp;
}
def code(x, y, z, t): t_1 = t * (0.5 + ((t * t) * (0.125 + ((t * t) * 0.020833333333333332)))) t_2 = (x * 0.5) - y t_3 = t * t_1 t_4 = y + (x * 0.5) tmp = 0 if (t * t) <= 2e+47: tmp = ((t_2 * math.sqrt((z * 2.0))) * (1.0 + (t_3 * ((t * t) * (t_1 * t_1))))) / (1.0 + (t_3 * (t_3 + -1.0))) else: tmp = t_2 * (t_4 * (((1.0 + t_3) / math.pow((z * 2.0), -0.5)) / t_4)) return tmp
function code(x, y, z, t) t_1 = Float64(t * Float64(0.5 + Float64(Float64(t * t) * Float64(0.125 + Float64(Float64(t * t) * 0.020833333333333332))))) t_2 = Float64(Float64(x * 0.5) - y) t_3 = Float64(t * t_1) t_4 = Float64(y + Float64(x * 0.5)) tmp = 0.0 if (Float64(t * t) <= 2e+47) tmp = Float64(Float64(Float64(t_2 * sqrt(Float64(z * 2.0))) * Float64(1.0 + Float64(t_3 * Float64(Float64(t * t) * Float64(t_1 * t_1))))) / Float64(1.0 + Float64(t_3 * Float64(t_3 + -1.0)))); else tmp = Float64(t_2 * Float64(t_4 * Float64(Float64(Float64(1.0 + t_3) / (Float64(z * 2.0) ^ -0.5)) / t_4))); end return tmp end
function tmp_2 = code(x, y, z, t) t_1 = t * (0.5 + ((t * t) * (0.125 + ((t * t) * 0.020833333333333332)))); t_2 = (x * 0.5) - y; t_3 = t * t_1; t_4 = y + (x * 0.5); tmp = 0.0; if ((t * t) <= 2e+47) tmp = ((t_2 * sqrt((z * 2.0))) * (1.0 + (t_3 * ((t * t) * (t_1 * t_1))))) / (1.0 + (t_3 * (t_3 + -1.0))); else tmp = t_2 * (t_4 * (((1.0 + t_3) / ((z * 2.0) ^ -0.5)) / t_4)); end tmp_2 = tmp; end
code[x_, y_, z_, t_] := Block[{t$95$1 = N[(t * N[(0.5 + N[(N[(t * t), $MachinePrecision] * N[(0.125 + N[(N[(t * t), $MachinePrecision] * 0.020833333333333332), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(N[(x * 0.5), $MachinePrecision] - y), $MachinePrecision]}, Block[{t$95$3 = N[(t * t$95$1), $MachinePrecision]}, Block[{t$95$4 = N[(y + N[(x * 0.5), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[N[(t * t), $MachinePrecision], 2e+47], N[(N[(N[(t$95$2 * N[Sqrt[N[(z * 2.0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] * N[(1.0 + N[(t$95$3 * N[(N[(t * t), $MachinePrecision] * N[(t$95$1 * t$95$1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(1.0 + N[(t$95$3 * N[(t$95$3 + -1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(t$95$2 * N[(t$95$4 * N[(N[(N[(1.0 + t$95$3), $MachinePrecision] / N[Power[N[(z * 2.0), $MachinePrecision], -0.5], $MachinePrecision]), $MachinePrecision] / t$95$4), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := t \cdot \left(0.5 + \left(t \cdot t\right) \cdot \left(0.125 + \left(t \cdot t\right) \cdot 0.020833333333333332\right)\right)\\
t_2 := x \cdot 0.5 - y\\
t_3 := t \cdot t\_1\\
t_4 := y + x \cdot 0.5\\
\mathbf{if}\;t \cdot t \leq 2 \cdot 10^{+47}:\\
\;\;\;\;\frac{\left(t\_2 \cdot \sqrt{z \cdot 2}\right) \cdot \left(1 + t\_3 \cdot \left(\left(t \cdot t\right) \cdot \left(t\_1 \cdot t\_1\right)\right)\right)}{1 + t\_3 \cdot \left(t\_3 + -1\right)}\\
\mathbf{else}:\\
\;\;\;\;t\_2 \cdot \left(t\_4 \cdot \frac{\frac{1 + t\_3}{{\left(z \cdot 2\right)}^{-0.5}}}{t\_4}\right)\\
\end{array}
\end{array}
if (*.f64 t t) < 2.0000000000000001e47Initial program 99.6%
Taylor expanded in t around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f6492.7%
Simplified92.7%
Applied egg-rr96.1%
if 2.0000000000000001e47 < (*.f64 t t) Initial program 100.0%
Taylor expanded in t around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f6496.0%
Simplified96.0%
associate-*l*N/A
*-lowering-*.f64N/A
--lowering--.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
sqrt-lowering-sqrt.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
Applied egg-rr96.8%
/-rgt-identityN/A
clear-numN/A
/-lowering-/.f64N/A
pow1/2N/A
*-commutativeN/A
pow-flipN/A
metadata-evalN/A
metadata-evalN/A
metadata-evalN/A
pow-lowering-pow.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
metadata-evalN/A
metadata-eval96.8%
Applied egg-rr96.8%
Applied egg-rr97.6%
Final simplification96.8%
(FPCore (x y z t)
:precision binary64
(let* ((t_1
(* t (+ 0.5 (* (* t t) (+ 0.125 (* (* t t) 0.020833333333333332))))))
(t_2 (- (* x 0.5) y))
(t_3 (sqrt (* z 2.0)))
(t_4 (* t t_1)))
(if (<= (* t t) 1e+93)
(/ (* (* t_2 t_3) (- 1.0 (* (* t t) (* t_1 t_1)))) (- 1.0 t_4))
(* t_3 (* t_2 (+ 1.0 t_4))))))
double code(double x, double y, double z, double t) {
double t_1 = t * (0.5 + ((t * t) * (0.125 + ((t * t) * 0.020833333333333332))));
double t_2 = (x * 0.5) - y;
double t_3 = sqrt((z * 2.0));
double t_4 = t * t_1;
double tmp;
if ((t * t) <= 1e+93) {
tmp = ((t_2 * t_3) * (1.0 - ((t * t) * (t_1 * t_1)))) / (1.0 - t_4);
} else {
tmp = t_3 * (t_2 * (1.0 + t_4));
}
return tmp;
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8) :: t_1
real(8) :: t_2
real(8) :: t_3
real(8) :: t_4
real(8) :: tmp
t_1 = t * (0.5d0 + ((t * t) * (0.125d0 + ((t * t) * 0.020833333333333332d0))))
t_2 = (x * 0.5d0) - y
t_3 = sqrt((z * 2.0d0))
t_4 = t * t_1
if ((t * t) <= 1d+93) then
tmp = ((t_2 * t_3) * (1.0d0 - ((t * t) * (t_1 * t_1)))) / (1.0d0 - t_4)
else
tmp = t_3 * (t_2 * (1.0d0 + t_4))
end if
code = tmp
end function
public static double code(double x, double y, double z, double t) {
double t_1 = t * (0.5 + ((t * t) * (0.125 + ((t * t) * 0.020833333333333332))));
double t_2 = (x * 0.5) - y;
double t_3 = Math.sqrt((z * 2.0));
double t_4 = t * t_1;
double tmp;
if ((t * t) <= 1e+93) {
tmp = ((t_2 * t_3) * (1.0 - ((t * t) * (t_1 * t_1)))) / (1.0 - t_4);
} else {
tmp = t_3 * (t_2 * (1.0 + t_4));
}
return tmp;
}
def code(x, y, z, t): t_1 = t * (0.5 + ((t * t) * (0.125 + ((t * t) * 0.020833333333333332)))) t_2 = (x * 0.5) - y t_3 = math.sqrt((z * 2.0)) t_4 = t * t_1 tmp = 0 if (t * t) <= 1e+93: tmp = ((t_2 * t_3) * (1.0 - ((t * t) * (t_1 * t_1)))) / (1.0 - t_4) else: tmp = t_3 * (t_2 * (1.0 + t_4)) return tmp
function code(x, y, z, t) t_1 = Float64(t * Float64(0.5 + Float64(Float64(t * t) * Float64(0.125 + Float64(Float64(t * t) * 0.020833333333333332))))) t_2 = Float64(Float64(x * 0.5) - y) t_3 = sqrt(Float64(z * 2.0)) t_4 = Float64(t * t_1) tmp = 0.0 if (Float64(t * t) <= 1e+93) tmp = Float64(Float64(Float64(t_2 * t_3) * Float64(1.0 - Float64(Float64(t * t) * Float64(t_1 * t_1)))) / Float64(1.0 - t_4)); else tmp = Float64(t_3 * Float64(t_2 * Float64(1.0 + t_4))); end return tmp end
function tmp_2 = code(x, y, z, t) t_1 = t * (0.5 + ((t * t) * (0.125 + ((t * t) * 0.020833333333333332)))); t_2 = (x * 0.5) - y; t_3 = sqrt((z * 2.0)); t_4 = t * t_1; tmp = 0.0; if ((t * t) <= 1e+93) tmp = ((t_2 * t_3) * (1.0 - ((t * t) * (t_1 * t_1)))) / (1.0 - t_4); else tmp = t_3 * (t_2 * (1.0 + t_4)); end tmp_2 = tmp; end
code[x_, y_, z_, t_] := Block[{t$95$1 = N[(t * N[(0.5 + N[(N[(t * t), $MachinePrecision] * N[(0.125 + N[(N[(t * t), $MachinePrecision] * 0.020833333333333332), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(N[(x * 0.5), $MachinePrecision] - y), $MachinePrecision]}, Block[{t$95$3 = N[Sqrt[N[(z * 2.0), $MachinePrecision]], $MachinePrecision]}, Block[{t$95$4 = N[(t * t$95$1), $MachinePrecision]}, If[LessEqual[N[(t * t), $MachinePrecision], 1e+93], N[(N[(N[(t$95$2 * t$95$3), $MachinePrecision] * N[(1.0 - N[(N[(t * t), $MachinePrecision] * N[(t$95$1 * t$95$1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(1.0 - t$95$4), $MachinePrecision]), $MachinePrecision], N[(t$95$3 * N[(t$95$2 * N[(1.0 + t$95$4), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := t \cdot \left(0.5 + \left(t \cdot t\right) \cdot \left(0.125 + \left(t \cdot t\right) \cdot 0.020833333333333332\right)\right)\\
t_2 := x \cdot 0.5 - y\\
t_3 := \sqrt{z \cdot 2}\\
t_4 := t \cdot t\_1\\
\mathbf{if}\;t \cdot t \leq 10^{+93}:\\
\;\;\;\;\frac{\left(t\_2 \cdot t\_3\right) \cdot \left(1 - \left(t \cdot t\right) \cdot \left(t\_1 \cdot t\_1\right)\right)}{1 - t\_4}\\
\mathbf{else}:\\
\;\;\;\;t\_3 \cdot \left(t\_2 \cdot \left(1 + t\_4\right)\right)\\
\end{array}
\end{array}
if (*.f64 t t) < 1.00000000000000004e93Initial program 99.6%
Taylor expanded in t around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f6489.8%
Simplified89.8%
*-commutativeN/A
flip-+N/A
associate-*l/N/A
/-lowering-/.f64N/A
Applied egg-rr93.6%
if 1.00000000000000004e93 < (*.f64 t t) Initial program 100.0%
Taylor expanded in t around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64100.0%
Simplified100.0%
*-commutativeN/A
associate-*r*N/A
*-lowering-*.f64N/A
Applied egg-rr100.0%
Final simplification96.4%
(FPCore (x y z t)
:precision binary64
(if (<= (* t t) 1e-14)
(* (- (* x 0.5) y) (sqrt (* z (* 2.0 (+ (* t t) 1.0)))))
(*
y
(*
(sqrt (* z (+ 2.0 (* t (* t (+ 2.0 (* t t)))))))
(+ (/ (* x 0.5) y) -1.0)))))
double code(double x, double y, double z, double t) {
double tmp;
if ((t * t) <= 1e-14) {
tmp = ((x * 0.5) - y) * sqrt((z * (2.0 * ((t * t) + 1.0))));
} else {
tmp = y * (sqrt((z * (2.0 + (t * (t * (2.0 + (t * t))))))) * (((x * 0.5) / y) + -1.0));
}
return tmp;
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8) :: tmp
if ((t * t) <= 1d-14) then
tmp = ((x * 0.5d0) - y) * sqrt((z * (2.0d0 * ((t * t) + 1.0d0))))
else
tmp = y * (sqrt((z * (2.0d0 + (t * (t * (2.0d0 + (t * t))))))) * (((x * 0.5d0) / y) + (-1.0d0)))
end if
code = tmp
end function
public static double code(double x, double y, double z, double t) {
double tmp;
if ((t * t) <= 1e-14) {
tmp = ((x * 0.5) - y) * Math.sqrt((z * (2.0 * ((t * t) + 1.0))));
} else {
tmp = y * (Math.sqrt((z * (2.0 + (t * (t * (2.0 + (t * t))))))) * (((x * 0.5) / y) + -1.0));
}
return tmp;
}
def code(x, y, z, t): tmp = 0 if (t * t) <= 1e-14: tmp = ((x * 0.5) - y) * math.sqrt((z * (2.0 * ((t * t) + 1.0)))) else: tmp = y * (math.sqrt((z * (2.0 + (t * (t * (2.0 + (t * t))))))) * (((x * 0.5) / y) + -1.0)) return tmp
function code(x, y, z, t) tmp = 0.0 if (Float64(t * t) <= 1e-14) tmp = Float64(Float64(Float64(x * 0.5) - y) * sqrt(Float64(z * Float64(2.0 * Float64(Float64(t * t) + 1.0))))); else tmp = Float64(y * Float64(sqrt(Float64(z * Float64(2.0 + Float64(t * Float64(t * Float64(2.0 + Float64(t * t))))))) * Float64(Float64(Float64(x * 0.5) / y) + -1.0))); end return tmp end
function tmp_2 = code(x, y, z, t) tmp = 0.0; if ((t * t) <= 1e-14) tmp = ((x * 0.5) - y) * sqrt((z * (2.0 * ((t * t) + 1.0)))); else tmp = y * (sqrt((z * (2.0 + (t * (t * (2.0 + (t * t))))))) * (((x * 0.5) / y) + -1.0)); end tmp_2 = tmp; end
code[x_, y_, z_, t_] := If[LessEqual[N[(t * t), $MachinePrecision], 1e-14], N[(N[(N[(x * 0.5), $MachinePrecision] - y), $MachinePrecision] * N[Sqrt[N[(z * N[(2.0 * N[(N[(t * t), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], N[(y * N[(N[Sqrt[N[(z * N[(2.0 + N[(t * N[(t * N[(2.0 + N[(t * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * N[(N[(N[(x * 0.5), $MachinePrecision] / y), $MachinePrecision] + -1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;t \cdot t \leq 10^{-14}:\\
\;\;\;\;\left(x \cdot 0.5 - y\right) \cdot \sqrt{z \cdot \left(2 \cdot \left(t \cdot t + 1\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;y \cdot \left(\sqrt{z \cdot \left(2 + t \cdot \left(t \cdot \left(2 + t \cdot t\right)\right)\right)} \cdot \left(\frac{x \cdot 0.5}{y} + -1\right)\right)\\
\end{array}
\end{array}
if (*.f64 t t) < 9.99999999999999999e-15Initial program 99.6%
associate-*l*N/A
*-lowering-*.f64N/A
--lowering--.f64N/A
*-lowering-*.f64N/A
exp-sqrtN/A
sqrt-unprodN/A
sqrt-lowering-sqrt.f64N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
exp-lowering-exp.f64N/A
*-lowering-*.f6499.6%
Applied egg-rr99.6%
Taylor expanded in t around 0
+-commutativeN/A
+-lowering-+.f64N/A
unpow2N/A
*-lowering-*.f6499.6%
Simplified99.6%
if 9.99999999999999999e-15 < (*.f64 t t) Initial program 100.0%
associate-*l*N/A
*-lowering-*.f64N/A
--lowering--.f64N/A
*-lowering-*.f64N/A
exp-sqrtN/A
sqrt-unprodN/A
sqrt-lowering-sqrt.f64N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
exp-lowering-exp.f64N/A
*-lowering-*.f64100.0%
Applied egg-rr100.0%
Taylor expanded in t around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
distribute-rgt-outN/A
associate-*r*N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
unpow2N/A
*-lowering-*.f6481.3%
Simplified81.3%
Taylor expanded in y around inf
*-lowering-*.f64N/A
+-commutativeN/A
associate-*r*N/A
distribute-rgt-outN/A
metadata-evalN/A
sub-negN/A
*-lowering-*.f64N/A
Simplified85.7%
Final simplification92.5%
(FPCore (x y z t)
:precision binary64
(*
(- (* x 0.5) y)
(*
(sqrt (* z 2.0))
(+
1.0
(*
t
(* t (+ 0.5 (* (* t t) (+ 0.125 (* (* t t) 0.020833333333333332))))))))))
double code(double x, double y, double z, double t) {
return ((x * 0.5) - y) * (sqrt((z * 2.0)) * (1.0 + (t * (t * (0.5 + ((t * t) * (0.125 + ((t * t) * 0.020833333333333332))))))));
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
code = ((x * 0.5d0) - y) * (sqrt((z * 2.0d0)) * (1.0d0 + (t * (t * (0.5d0 + ((t * t) * (0.125d0 + ((t * t) * 0.020833333333333332d0))))))))
end function
public static double code(double x, double y, double z, double t) {
return ((x * 0.5) - y) * (Math.sqrt((z * 2.0)) * (1.0 + (t * (t * (0.5 + ((t * t) * (0.125 + ((t * t) * 0.020833333333333332))))))));
}
def code(x, y, z, t): return ((x * 0.5) - y) * (math.sqrt((z * 2.0)) * (1.0 + (t * (t * (0.5 + ((t * t) * (0.125 + ((t * t) * 0.020833333333333332))))))))
function code(x, y, z, t) return Float64(Float64(Float64(x * 0.5) - y) * Float64(sqrt(Float64(z * 2.0)) * Float64(1.0 + Float64(t * Float64(t * Float64(0.5 + Float64(Float64(t * t) * Float64(0.125 + Float64(Float64(t * t) * 0.020833333333333332))))))))) end
function tmp = code(x, y, z, t) tmp = ((x * 0.5) - y) * (sqrt((z * 2.0)) * (1.0 + (t * (t * (0.5 + ((t * t) * (0.125 + ((t * t) * 0.020833333333333332)))))))); end
code[x_, y_, z_, t_] := N[(N[(N[(x * 0.5), $MachinePrecision] - y), $MachinePrecision] * N[(N[Sqrt[N[(z * 2.0), $MachinePrecision]], $MachinePrecision] * N[(1.0 + N[(t * N[(t * N[(0.5 + N[(N[(t * t), $MachinePrecision] * N[(0.125 + N[(N[(t * t), $MachinePrecision] * 0.020833333333333332), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\left(x \cdot 0.5 - y\right) \cdot \left(\sqrt{z \cdot 2} \cdot \left(1 + t \cdot \left(t \cdot \left(0.5 + \left(t \cdot t\right) \cdot \left(0.125 + \left(t \cdot t\right) \cdot 0.020833333333333332\right)\right)\right)\right)\right)
\end{array}
Initial program 99.8%
Taylor expanded in t around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f6494.2%
Simplified94.2%
associate-*l*N/A
*-lowering-*.f64N/A
--lowering--.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
sqrt-lowering-sqrt.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
Applied egg-rr94.5%
(FPCore (x y z t)
:precision binary64
(*
(* (- (* x 0.5) y) (sqrt (* z 2.0)))
(+
1.0
(*
(* t t)
(+ 0.5 (* t (* t (+ 0.125 (* t (* t 0.020833333333333332))))))))))
double code(double x, double y, double z, double t) {
return (((x * 0.5) - y) * sqrt((z * 2.0))) * (1.0 + ((t * t) * (0.5 + (t * (t * (0.125 + (t * (t * 0.020833333333333332))))))));
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
code = (((x * 0.5d0) - y) * sqrt((z * 2.0d0))) * (1.0d0 + ((t * t) * (0.5d0 + (t * (t * (0.125d0 + (t * (t * 0.020833333333333332d0))))))))
end function
public static double code(double x, double y, double z, double t) {
return (((x * 0.5) - y) * Math.sqrt((z * 2.0))) * (1.0 + ((t * t) * (0.5 + (t * (t * (0.125 + (t * (t * 0.020833333333333332))))))));
}
def code(x, y, z, t): return (((x * 0.5) - y) * math.sqrt((z * 2.0))) * (1.0 + ((t * t) * (0.5 + (t * (t * (0.125 + (t * (t * 0.020833333333333332))))))))
function code(x, y, z, t) return Float64(Float64(Float64(Float64(x * 0.5) - y) * sqrt(Float64(z * 2.0))) * Float64(1.0 + Float64(Float64(t * t) * Float64(0.5 + Float64(t * Float64(t * Float64(0.125 + Float64(t * Float64(t * 0.020833333333333332))))))))) end
function tmp = code(x, y, z, t) tmp = (((x * 0.5) - y) * sqrt((z * 2.0))) * (1.0 + ((t * t) * (0.5 + (t * (t * (0.125 + (t * (t * 0.020833333333333332)))))))); end
code[x_, y_, z_, t_] := N[(N[(N[(N[(x * 0.5), $MachinePrecision] - y), $MachinePrecision] * N[Sqrt[N[(z * 2.0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] * N[(1.0 + N[(N[(t * t), $MachinePrecision] * N[(0.5 + N[(t * N[(t * N[(0.125 + N[(t * N[(t * 0.020833333333333332), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\left(\left(x \cdot 0.5 - y\right) \cdot \sqrt{z \cdot 2}\right) \cdot \left(1 + \left(t \cdot t\right) \cdot \left(0.5 + t \cdot \left(t \cdot \left(0.125 + t \cdot \left(t \cdot 0.020833333333333332\right)\right)\right)\right)\right)
\end{array}
Initial program 99.8%
Taylor expanded in t around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f6494.2%
Simplified94.2%
(FPCore (x y z t)
:precision binary64
(*
(sqrt (* z 2.0))
(*
(- (* x 0.5) y)
(+
1.0
(*
t
(* t (+ 0.5 (* (* t t) (+ 0.125 (* (* t t) 0.020833333333333332))))))))))
double code(double x, double y, double z, double t) {
return sqrt((z * 2.0)) * (((x * 0.5) - y) * (1.0 + (t * (t * (0.5 + ((t * t) * (0.125 + ((t * t) * 0.020833333333333332))))))));
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
code = sqrt((z * 2.0d0)) * (((x * 0.5d0) - y) * (1.0d0 + (t * (t * (0.5d0 + ((t * t) * (0.125d0 + ((t * t) * 0.020833333333333332d0))))))))
end function
public static double code(double x, double y, double z, double t) {
return Math.sqrt((z * 2.0)) * (((x * 0.5) - y) * (1.0 + (t * (t * (0.5 + ((t * t) * (0.125 + ((t * t) * 0.020833333333333332))))))));
}
def code(x, y, z, t): return math.sqrt((z * 2.0)) * (((x * 0.5) - y) * (1.0 + (t * (t * (0.5 + ((t * t) * (0.125 + ((t * t) * 0.020833333333333332))))))))
function code(x, y, z, t) return Float64(sqrt(Float64(z * 2.0)) * Float64(Float64(Float64(x * 0.5) - y) * Float64(1.0 + Float64(t * Float64(t * Float64(0.5 + Float64(Float64(t * t) * Float64(0.125 + Float64(Float64(t * t) * 0.020833333333333332))))))))) end
function tmp = code(x, y, z, t) tmp = sqrt((z * 2.0)) * (((x * 0.5) - y) * (1.0 + (t * (t * (0.5 + ((t * t) * (0.125 + ((t * t) * 0.020833333333333332)))))))); end
code[x_, y_, z_, t_] := N[(N[Sqrt[N[(z * 2.0), $MachinePrecision]], $MachinePrecision] * N[(N[(N[(x * 0.5), $MachinePrecision] - y), $MachinePrecision] * N[(1.0 + N[(t * N[(t * N[(0.5 + N[(N[(t * t), $MachinePrecision] * N[(0.125 + N[(N[(t * t), $MachinePrecision] * 0.020833333333333332), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\sqrt{z \cdot 2} \cdot \left(\left(x \cdot 0.5 - y\right) \cdot \left(1 + t \cdot \left(t \cdot \left(0.5 + \left(t \cdot t\right) \cdot \left(0.125 + \left(t \cdot t\right) \cdot 0.020833333333333332\right)\right)\right)\right)\right)
\end{array}
Initial program 99.8%
Taylor expanded in t around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f6494.2%
Simplified94.2%
*-commutativeN/A
associate-*r*N/A
*-lowering-*.f64N/A
Applied egg-rr94.2%
Final simplification94.2%
(FPCore (x y z t)
:precision binary64
(let* ((t_1 (+ 2.0 (* t t))))
(if (<= (* t t) 5000000.0)
(* (- (* x 0.5) y) (sqrt (+ (* z 2.0) (* t_1 (* z (* t t))))))
(* x (* (sqrt (* z (+ 2.0 (* t (* t t_1))))) (- 0.5 (/ y x)))))))
double code(double x, double y, double z, double t) {
double t_1 = 2.0 + (t * t);
double tmp;
if ((t * t) <= 5000000.0) {
tmp = ((x * 0.5) - y) * sqrt(((z * 2.0) + (t_1 * (z * (t * t)))));
} else {
tmp = x * (sqrt((z * (2.0 + (t * (t * t_1))))) * (0.5 - (y / x)));
}
return tmp;
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8) :: t_1
real(8) :: tmp
t_1 = 2.0d0 + (t * t)
if ((t * t) <= 5000000.0d0) then
tmp = ((x * 0.5d0) - y) * sqrt(((z * 2.0d0) + (t_1 * (z * (t * t)))))
else
tmp = x * (sqrt((z * (2.0d0 + (t * (t * t_1))))) * (0.5d0 - (y / x)))
end if
code = tmp
end function
public static double code(double x, double y, double z, double t) {
double t_1 = 2.0 + (t * t);
double tmp;
if ((t * t) <= 5000000.0) {
tmp = ((x * 0.5) - y) * Math.sqrt(((z * 2.0) + (t_1 * (z * (t * t)))));
} else {
tmp = x * (Math.sqrt((z * (2.0 + (t * (t * t_1))))) * (0.5 - (y / x)));
}
return tmp;
}
def code(x, y, z, t): t_1 = 2.0 + (t * t) tmp = 0 if (t * t) <= 5000000.0: tmp = ((x * 0.5) - y) * math.sqrt(((z * 2.0) + (t_1 * (z * (t * t))))) else: tmp = x * (math.sqrt((z * (2.0 + (t * (t * t_1))))) * (0.5 - (y / x))) return tmp
function code(x, y, z, t) t_1 = Float64(2.0 + Float64(t * t)) tmp = 0.0 if (Float64(t * t) <= 5000000.0) tmp = Float64(Float64(Float64(x * 0.5) - y) * sqrt(Float64(Float64(z * 2.0) + Float64(t_1 * Float64(z * Float64(t * t)))))); else tmp = Float64(x * Float64(sqrt(Float64(z * Float64(2.0 + Float64(t * Float64(t * t_1))))) * Float64(0.5 - Float64(y / x)))); end return tmp end
function tmp_2 = code(x, y, z, t) t_1 = 2.0 + (t * t); tmp = 0.0; if ((t * t) <= 5000000.0) tmp = ((x * 0.5) - y) * sqrt(((z * 2.0) + (t_1 * (z * (t * t))))); else tmp = x * (sqrt((z * (2.0 + (t * (t * t_1))))) * (0.5 - (y / x))); end tmp_2 = tmp; end
code[x_, y_, z_, t_] := Block[{t$95$1 = N[(2.0 + N[(t * t), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[N[(t * t), $MachinePrecision], 5000000.0], N[(N[(N[(x * 0.5), $MachinePrecision] - y), $MachinePrecision] * N[Sqrt[N[(N[(z * 2.0), $MachinePrecision] + N[(t$95$1 * N[(z * N[(t * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], N[(x * N[(N[Sqrt[N[(z * N[(2.0 + N[(t * N[(t * t$95$1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * N[(0.5 - N[(y / x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := 2 + t \cdot t\\
\mathbf{if}\;t \cdot t \leq 5000000:\\
\;\;\;\;\left(x \cdot 0.5 - y\right) \cdot \sqrt{z \cdot 2 + t\_1 \cdot \left(z \cdot \left(t \cdot t\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;x \cdot \left(\sqrt{z \cdot \left(2 + t \cdot \left(t \cdot t\_1\right)\right)} \cdot \left(0.5 - \frac{y}{x}\right)\right)\\
\end{array}
\end{array}
if (*.f64 t t) < 5e6Initial program 99.6%
associate-*l*N/A
*-lowering-*.f64N/A
--lowering--.f64N/A
*-lowering-*.f64N/A
exp-sqrtN/A
sqrt-unprodN/A
sqrt-lowering-sqrt.f64N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
exp-lowering-exp.f64N/A
*-lowering-*.f6499.6%
Applied egg-rr99.6%
Taylor expanded in t around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
distribute-rgt-outN/A
associate-*r*N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
unpow2N/A
*-lowering-*.f6496.4%
Simplified96.4%
if 5e6 < (*.f64 t t) Initial program 100.0%
associate-*l*N/A
*-lowering-*.f64N/A
--lowering--.f64N/A
*-lowering-*.f64N/A
exp-sqrtN/A
sqrt-unprodN/A
sqrt-lowering-sqrt.f64N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
exp-lowering-exp.f64N/A
*-lowering-*.f64100.0%
Applied egg-rr100.0%
Taylor expanded in t around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
distribute-rgt-outN/A
associate-*r*N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
unpow2N/A
*-lowering-*.f6483.9%
Simplified83.9%
Taylor expanded in x around inf
*-lowering-*.f64N/A
associate-*r*N/A
distribute-rgt-outN/A
+-commutativeN/A
*-lowering-*.f64N/A
Simplified88.4%
Final simplification92.5%
(FPCore (x y z t) :precision binary64 (if (<= (* t t) 0.4) (* (* (- (* x 0.5) y) (sqrt (* z 2.0))) (+ 1.0 (* 0.5 (* t t)))) (* x (* (sqrt (* z (+ 2.0 (* t (* t (+ 2.0 (* t t))))))) (- 0.5 (/ y x))))))
double code(double x, double y, double z, double t) {
double tmp;
if ((t * t) <= 0.4) {
tmp = (((x * 0.5) - y) * sqrt((z * 2.0))) * (1.0 + (0.5 * (t * t)));
} else {
tmp = x * (sqrt((z * (2.0 + (t * (t * (2.0 + (t * t))))))) * (0.5 - (y / x)));
}
return tmp;
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8) :: tmp
if ((t * t) <= 0.4d0) then
tmp = (((x * 0.5d0) - y) * sqrt((z * 2.0d0))) * (1.0d0 + (0.5d0 * (t * t)))
else
tmp = x * (sqrt((z * (2.0d0 + (t * (t * (2.0d0 + (t * t))))))) * (0.5d0 - (y / x)))
end if
code = tmp
end function
public static double code(double x, double y, double z, double t) {
double tmp;
if ((t * t) <= 0.4) {
tmp = (((x * 0.5) - y) * Math.sqrt((z * 2.0))) * (1.0 + (0.5 * (t * t)));
} else {
tmp = x * (Math.sqrt((z * (2.0 + (t * (t * (2.0 + (t * t))))))) * (0.5 - (y / x)));
}
return tmp;
}
def code(x, y, z, t): tmp = 0 if (t * t) <= 0.4: tmp = (((x * 0.5) - y) * math.sqrt((z * 2.0))) * (1.0 + (0.5 * (t * t))) else: tmp = x * (math.sqrt((z * (2.0 + (t * (t * (2.0 + (t * t))))))) * (0.5 - (y / x))) return tmp
function code(x, y, z, t) tmp = 0.0 if (Float64(t * t) <= 0.4) tmp = Float64(Float64(Float64(Float64(x * 0.5) - y) * sqrt(Float64(z * 2.0))) * Float64(1.0 + Float64(0.5 * Float64(t * t)))); else tmp = Float64(x * Float64(sqrt(Float64(z * Float64(2.0 + Float64(t * Float64(t * Float64(2.0 + Float64(t * t))))))) * Float64(0.5 - Float64(y / x)))); end return tmp end
function tmp_2 = code(x, y, z, t) tmp = 0.0; if ((t * t) <= 0.4) tmp = (((x * 0.5) - y) * sqrt((z * 2.0))) * (1.0 + (0.5 * (t * t))); else tmp = x * (sqrt((z * (2.0 + (t * (t * (2.0 + (t * t))))))) * (0.5 - (y / x))); end tmp_2 = tmp; end
code[x_, y_, z_, t_] := If[LessEqual[N[(t * t), $MachinePrecision], 0.4], N[(N[(N[(N[(x * 0.5), $MachinePrecision] - y), $MachinePrecision] * N[Sqrt[N[(z * 2.0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] * N[(1.0 + N[(0.5 * N[(t * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x * N[(N[Sqrt[N[(z * N[(2.0 + N[(t * N[(t * N[(2.0 + N[(t * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * N[(0.5 - N[(y / x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;t \cdot t \leq 0.4:\\
\;\;\;\;\left(\left(x \cdot 0.5 - y\right) \cdot \sqrt{z \cdot 2}\right) \cdot \left(1 + 0.5 \cdot \left(t \cdot t\right)\right)\\
\mathbf{else}:\\
\;\;\;\;x \cdot \left(\sqrt{z \cdot \left(2 + t \cdot \left(t \cdot \left(2 + t \cdot t\right)\right)\right)} \cdot \left(0.5 - \frac{y}{x}\right)\right)\\
\end{array}
\end{array}
if (*.f64 t t) < 0.40000000000000002Initial program 99.6%
Taylor expanded in t around 0
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6498.5%
Simplified98.5%
if 0.40000000000000002 < (*.f64 t t) Initial program 100.0%
associate-*l*N/A
*-lowering-*.f64N/A
--lowering--.f64N/A
*-lowering-*.f64N/A
exp-sqrtN/A
sqrt-unprodN/A
sqrt-lowering-sqrt.f64N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
exp-lowering-exp.f64N/A
*-lowering-*.f64100.0%
Applied egg-rr100.0%
Taylor expanded in t around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
distribute-rgt-outN/A
associate-*r*N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
unpow2N/A
*-lowering-*.f6482.0%
Simplified82.0%
Taylor expanded in x around inf
*-lowering-*.f64N/A
associate-*r*N/A
distribute-rgt-outN/A
+-commutativeN/A
*-lowering-*.f64N/A
Simplified86.4%
Final simplification92.5%
(FPCore (x y z t) :precision binary64 (if (<= (* t t) 0.4) (* (* (- (* x 0.5) y) (sqrt (* z 2.0))) (+ 1.0 (* 0.5 (* t t)))) (* (* y (+ (/ (* x 0.5) y) -1.0)) (sqrt (* t (* t (* z (* t t))))))))
double code(double x, double y, double z, double t) {
double tmp;
if ((t * t) <= 0.4) {
tmp = (((x * 0.5) - y) * sqrt((z * 2.0))) * (1.0 + (0.5 * (t * t)));
} else {
tmp = (y * (((x * 0.5) / y) + -1.0)) * sqrt((t * (t * (z * (t * t)))));
}
return tmp;
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8) :: tmp
if ((t * t) <= 0.4d0) then
tmp = (((x * 0.5d0) - y) * sqrt((z * 2.0d0))) * (1.0d0 + (0.5d0 * (t * t)))
else
tmp = (y * (((x * 0.5d0) / y) + (-1.0d0))) * sqrt((t * (t * (z * (t * t)))))
end if
code = tmp
end function
public static double code(double x, double y, double z, double t) {
double tmp;
if ((t * t) <= 0.4) {
tmp = (((x * 0.5) - y) * Math.sqrt((z * 2.0))) * (1.0 + (0.5 * (t * t)));
} else {
tmp = (y * (((x * 0.5) / y) + -1.0)) * Math.sqrt((t * (t * (z * (t * t)))));
}
return tmp;
}
def code(x, y, z, t): tmp = 0 if (t * t) <= 0.4: tmp = (((x * 0.5) - y) * math.sqrt((z * 2.0))) * (1.0 + (0.5 * (t * t))) else: tmp = (y * (((x * 0.5) / y) + -1.0)) * math.sqrt((t * (t * (z * (t * t))))) return tmp
function code(x, y, z, t) tmp = 0.0 if (Float64(t * t) <= 0.4) tmp = Float64(Float64(Float64(Float64(x * 0.5) - y) * sqrt(Float64(z * 2.0))) * Float64(1.0 + Float64(0.5 * Float64(t * t)))); else tmp = Float64(Float64(y * Float64(Float64(Float64(x * 0.5) / y) + -1.0)) * sqrt(Float64(t * Float64(t * Float64(z * Float64(t * t)))))); end return tmp end
function tmp_2 = code(x, y, z, t) tmp = 0.0; if ((t * t) <= 0.4) tmp = (((x * 0.5) - y) * sqrt((z * 2.0))) * (1.0 + (0.5 * (t * t))); else tmp = (y * (((x * 0.5) / y) + -1.0)) * sqrt((t * (t * (z * (t * t))))); end tmp_2 = tmp; end
code[x_, y_, z_, t_] := If[LessEqual[N[(t * t), $MachinePrecision], 0.4], N[(N[(N[(N[(x * 0.5), $MachinePrecision] - y), $MachinePrecision] * N[Sqrt[N[(z * 2.0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] * N[(1.0 + N[(0.5 * N[(t * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(y * N[(N[(N[(x * 0.5), $MachinePrecision] / y), $MachinePrecision] + -1.0), $MachinePrecision]), $MachinePrecision] * N[Sqrt[N[(t * N[(t * N[(z * N[(t * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;t \cdot t \leq 0.4:\\
\;\;\;\;\left(\left(x \cdot 0.5 - y\right) \cdot \sqrt{z \cdot 2}\right) \cdot \left(1 + 0.5 \cdot \left(t \cdot t\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\left(y \cdot \left(\frac{x \cdot 0.5}{y} + -1\right)\right) \cdot \sqrt{t \cdot \left(t \cdot \left(z \cdot \left(t \cdot t\right)\right)\right)}\\
\end{array}
\end{array}
if (*.f64 t t) < 0.40000000000000002Initial program 99.6%
Taylor expanded in t around 0
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6498.5%
Simplified98.5%
if 0.40000000000000002 < (*.f64 t t) Initial program 100.0%
associate-*l*N/A
*-lowering-*.f64N/A
--lowering--.f64N/A
*-lowering-*.f64N/A
exp-sqrtN/A
sqrt-unprodN/A
sqrt-lowering-sqrt.f64N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
exp-lowering-exp.f64N/A
*-lowering-*.f64100.0%
Applied egg-rr100.0%
Taylor expanded in t around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
distribute-rgt-outN/A
associate-*r*N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
unpow2N/A
*-lowering-*.f6482.0%
Simplified82.0%
Taylor expanded in t around inf
metadata-evalN/A
pow-sqrN/A
associate-*r*N/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6482.0%
Simplified82.0%
Taylor expanded in y around inf
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f64N/A
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f6482.8%
Simplified82.8%
Final simplification90.7%
(FPCore (x y z t)
:precision binary64
(let* ((t_1 (* (sqrt z) (- 0.0 (* y (* t t))))))
(if (<= t 220000000000.0)
(/ (* x 0.5) (pow (* z 2.0) -0.5))
(if (<= t 9.4e+148)
t_1
(if (<= t 2e+262) (* (sqrt z) (* x (* 0.5 (* t t)))) t_1)))))
double code(double x, double y, double z, double t) {
double t_1 = sqrt(z) * (0.0 - (y * (t * t)));
double tmp;
if (t <= 220000000000.0) {
tmp = (x * 0.5) / pow((z * 2.0), -0.5);
} else if (t <= 9.4e+148) {
tmp = t_1;
} else if (t <= 2e+262) {
tmp = sqrt(z) * (x * (0.5 * (t * t)));
} else {
tmp = t_1;
}
return tmp;
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8) :: t_1
real(8) :: tmp
t_1 = sqrt(z) * (0.0d0 - (y * (t * t)))
if (t <= 220000000000.0d0) then
tmp = (x * 0.5d0) / ((z * 2.0d0) ** (-0.5d0))
else if (t <= 9.4d+148) then
tmp = t_1
else if (t <= 2d+262) then
tmp = sqrt(z) * (x * (0.5d0 * (t * t)))
else
tmp = t_1
end if
code = tmp
end function
public static double code(double x, double y, double z, double t) {
double t_1 = Math.sqrt(z) * (0.0 - (y * (t * t)));
double tmp;
if (t <= 220000000000.0) {
tmp = (x * 0.5) / Math.pow((z * 2.0), -0.5);
} else if (t <= 9.4e+148) {
tmp = t_1;
} else if (t <= 2e+262) {
tmp = Math.sqrt(z) * (x * (0.5 * (t * t)));
} else {
tmp = t_1;
}
return tmp;
}
def code(x, y, z, t): t_1 = math.sqrt(z) * (0.0 - (y * (t * t))) tmp = 0 if t <= 220000000000.0: tmp = (x * 0.5) / math.pow((z * 2.0), -0.5) elif t <= 9.4e+148: tmp = t_1 elif t <= 2e+262: tmp = math.sqrt(z) * (x * (0.5 * (t * t))) else: tmp = t_1 return tmp
function code(x, y, z, t) t_1 = Float64(sqrt(z) * Float64(0.0 - Float64(y * Float64(t * t)))) tmp = 0.0 if (t <= 220000000000.0) tmp = Float64(Float64(x * 0.5) / (Float64(z * 2.0) ^ -0.5)); elseif (t <= 9.4e+148) tmp = t_1; elseif (t <= 2e+262) tmp = Float64(sqrt(z) * Float64(x * Float64(0.5 * Float64(t * t)))); else tmp = t_1; end return tmp end
function tmp_2 = code(x, y, z, t) t_1 = sqrt(z) * (0.0 - (y * (t * t))); tmp = 0.0; if (t <= 220000000000.0) tmp = (x * 0.5) / ((z * 2.0) ^ -0.5); elseif (t <= 9.4e+148) tmp = t_1; elseif (t <= 2e+262) tmp = sqrt(z) * (x * (0.5 * (t * t))); else tmp = t_1; end tmp_2 = tmp; end
code[x_, y_, z_, t_] := Block[{t$95$1 = N[(N[Sqrt[z], $MachinePrecision] * N[(0.0 - N[(y * N[(t * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t, 220000000000.0], N[(N[(x * 0.5), $MachinePrecision] / N[Power[N[(z * 2.0), $MachinePrecision], -0.5], $MachinePrecision]), $MachinePrecision], If[LessEqual[t, 9.4e+148], t$95$1, If[LessEqual[t, 2e+262], N[(N[Sqrt[z], $MachinePrecision] * N[(x * N[(0.5 * N[(t * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$1]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := \sqrt{z} \cdot \left(0 - y \cdot \left(t \cdot t\right)\right)\\
\mathbf{if}\;t \leq 220000000000:\\
\;\;\;\;\frac{x \cdot 0.5}{{\left(z \cdot 2\right)}^{-0.5}}\\
\mathbf{elif}\;t \leq 9.4 \cdot 10^{+148}:\\
\;\;\;\;t\_1\\
\mathbf{elif}\;t \leq 2 \cdot 10^{+262}:\\
\;\;\;\;\sqrt{z} \cdot \left(x \cdot \left(0.5 \cdot \left(t \cdot t\right)\right)\right)\\
\mathbf{else}:\\
\;\;\;\;t\_1\\
\end{array}
\end{array}
if t < 2.2e11Initial program 99.7%
Taylor expanded in t around 0
Simplified70.6%
Taylor expanded in x around inf
*-lowering-*.f6437.0%
Simplified37.0%
*-rgt-identityN/A
*-commutativeN/A
*-commutativeN/A
*-commutativeN/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
sqrt-lowering-sqrt.f64N/A
*-commutativeN/A
*-lowering-*.f6437.0%
Applied egg-rr37.0%
*-commutativeN/A
*-commutativeN/A
associate-*r*N/A
pow1/2N/A
metadata-evalN/A
pow-flipN/A
un-div-invN/A
/-lowering-/.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
pow-lowering-pow.f64N/A
*-lowering-*.f6437.0%
Applied egg-rr37.0%
if 2.2e11 < t < 9.3999999999999994e148 or 2e262 < t Initial program 100.0%
associate-*l*N/A
*-lowering-*.f64N/A
--lowering--.f64N/A
*-lowering-*.f64N/A
exp-sqrtN/A
sqrt-unprodN/A
sqrt-lowering-sqrt.f64N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
exp-lowering-exp.f64N/A
*-lowering-*.f64100.0%
Applied egg-rr100.0%
Taylor expanded in t around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
distribute-rgt-outN/A
associate-*r*N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
unpow2N/A
*-lowering-*.f6464.0%
Simplified64.0%
Taylor expanded in t around inf
metadata-evalN/A
pow-sqrN/A
associate-*r*N/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6464.0%
Simplified64.0%
Taylor expanded in x around 0
mul-1-negN/A
*-commutativeN/A
distribute-rgt-neg-inN/A
mul-1-negN/A
*-lowering-*.f64N/A
sqrt-lowering-sqrt.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6445.4%
Simplified45.4%
if 9.3999999999999994e148 < t < 2e262Initial program 100.0%
associate-*l*N/A
*-lowering-*.f64N/A
--lowering--.f64N/A
*-lowering-*.f64N/A
exp-sqrtN/A
sqrt-unprodN/A
sqrt-lowering-sqrt.f64N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
exp-lowering-exp.f64N/A
*-lowering-*.f64100.0%
Applied egg-rr100.0%
Taylor expanded in t around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
distribute-rgt-outN/A
associate-*r*N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
unpow2N/A
*-lowering-*.f64100.0%
Simplified100.0%
Taylor expanded in t around inf
metadata-evalN/A
pow-sqrN/A
associate-*r*N/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64100.0%
Simplified100.0%
Taylor expanded in x around inf
*-commutativeN/A
*-commutativeN/A
associate-*l*N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
sqrt-lowering-sqrt.f64N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6478.6%
Simplified78.6%
Final simplification42.6%
(FPCore (x y z t)
:precision binary64
(let* ((t_1 (- (* x 0.5) y)))
(if (<= (* t t) 0.4)
(* (* t_1 (sqrt (* z 2.0))) (+ 1.0 (* 0.5 (* t t))))
(* t_1 (sqrt (* t (* z (* t (* t t)))))))))
double code(double x, double y, double z, double t) {
double t_1 = (x * 0.5) - y;
double tmp;
if ((t * t) <= 0.4) {
tmp = (t_1 * sqrt((z * 2.0))) * (1.0 + (0.5 * (t * t)));
} else {
tmp = t_1 * sqrt((t * (z * (t * (t * t)))));
}
return tmp;
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8) :: t_1
real(8) :: tmp
t_1 = (x * 0.5d0) - y
if ((t * t) <= 0.4d0) then
tmp = (t_1 * sqrt((z * 2.0d0))) * (1.0d0 + (0.5d0 * (t * t)))
else
tmp = t_1 * sqrt((t * (z * (t * (t * t)))))
end if
code = tmp
end function
public static double code(double x, double y, double z, double t) {
double t_1 = (x * 0.5) - y;
double tmp;
if ((t * t) <= 0.4) {
tmp = (t_1 * Math.sqrt((z * 2.0))) * (1.0 + (0.5 * (t * t)));
} else {
tmp = t_1 * Math.sqrt((t * (z * (t * (t * t)))));
}
return tmp;
}
def code(x, y, z, t): t_1 = (x * 0.5) - y tmp = 0 if (t * t) <= 0.4: tmp = (t_1 * math.sqrt((z * 2.0))) * (1.0 + (0.5 * (t * t))) else: tmp = t_1 * math.sqrt((t * (z * (t * (t * t))))) return tmp
function code(x, y, z, t) t_1 = Float64(Float64(x * 0.5) - y) tmp = 0.0 if (Float64(t * t) <= 0.4) tmp = Float64(Float64(t_1 * sqrt(Float64(z * 2.0))) * Float64(1.0 + Float64(0.5 * Float64(t * t)))); else tmp = Float64(t_1 * sqrt(Float64(t * Float64(z * Float64(t * Float64(t * t)))))); end return tmp end
function tmp_2 = code(x, y, z, t) t_1 = (x * 0.5) - y; tmp = 0.0; if ((t * t) <= 0.4) tmp = (t_1 * sqrt((z * 2.0))) * (1.0 + (0.5 * (t * t))); else tmp = t_1 * sqrt((t * (z * (t * (t * t))))); end tmp_2 = tmp; end
code[x_, y_, z_, t_] := Block[{t$95$1 = N[(N[(x * 0.5), $MachinePrecision] - y), $MachinePrecision]}, If[LessEqual[N[(t * t), $MachinePrecision], 0.4], N[(N[(t$95$1 * N[Sqrt[N[(z * 2.0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] * N[(1.0 + N[(0.5 * N[(t * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(t$95$1 * N[Sqrt[N[(t * N[(z * N[(t * N[(t * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x \cdot 0.5 - y\\
\mathbf{if}\;t \cdot t \leq 0.4:\\
\;\;\;\;\left(t\_1 \cdot \sqrt{z \cdot 2}\right) \cdot \left(1 + 0.5 \cdot \left(t \cdot t\right)\right)\\
\mathbf{else}:\\
\;\;\;\;t\_1 \cdot \sqrt{t \cdot \left(z \cdot \left(t \cdot \left(t \cdot t\right)\right)\right)}\\
\end{array}
\end{array}
if (*.f64 t t) < 0.40000000000000002Initial program 99.6%
Taylor expanded in t around 0
+-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6498.5%
Simplified98.5%
if 0.40000000000000002 < (*.f64 t t) Initial program 100.0%
associate-*l*N/A
*-lowering-*.f64N/A
--lowering--.f64N/A
*-lowering-*.f64N/A
exp-sqrtN/A
sqrt-unprodN/A
sqrt-lowering-sqrt.f64N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
exp-lowering-exp.f64N/A
*-lowering-*.f64100.0%
Applied egg-rr100.0%
Taylor expanded in t around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
distribute-rgt-outN/A
associate-*r*N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
unpow2N/A
*-lowering-*.f6482.0%
Simplified82.0%
Taylor expanded in t around inf
metadata-evalN/A
pow-sqrN/A
associate-*r*N/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6482.0%
Simplified82.0%
associate-*r*N/A
cube-unmultN/A
*-lowering-*.f64N/A
cube-unmultN/A
*-lowering-*.f64N/A
*-lowering-*.f6482.8%
Applied egg-rr82.8%
Final simplification90.7%
(FPCore (x y z t) :precision binary64 (* (* (- (* x 0.5) y) (sqrt (* z 2.0))) (+ 1.0 (* (* t t) (+ 0.5 (* (* t t) 0.125))))))
double code(double x, double y, double z, double t) {
return (((x * 0.5) - y) * sqrt((z * 2.0))) * (1.0 + ((t * t) * (0.5 + ((t * t) * 0.125))));
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
code = (((x * 0.5d0) - y) * sqrt((z * 2.0d0))) * (1.0d0 + ((t * t) * (0.5d0 + ((t * t) * 0.125d0))))
end function
public static double code(double x, double y, double z, double t) {
return (((x * 0.5) - y) * Math.sqrt((z * 2.0))) * (1.0 + ((t * t) * (0.5 + ((t * t) * 0.125))));
}
def code(x, y, z, t): return (((x * 0.5) - y) * math.sqrt((z * 2.0))) * (1.0 + ((t * t) * (0.5 + ((t * t) * 0.125))))
function code(x, y, z, t) return Float64(Float64(Float64(Float64(x * 0.5) - y) * sqrt(Float64(z * 2.0))) * Float64(1.0 + Float64(Float64(t * t) * Float64(0.5 + Float64(Float64(t * t) * 0.125))))) end
function tmp = code(x, y, z, t) tmp = (((x * 0.5) - y) * sqrt((z * 2.0))) * (1.0 + ((t * t) * (0.5 + ((t * t) * 0.125)))); end
code[x_, y_, z_, t_] := N[(N[(N[(N[(x * 0.5), $MachinePrecision] - y), $MachinePrecision] * N[Sqrt[N[(z * 2.0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] * N[(1.0 + N[(N[(t * t), $MachinePrecision] * N[(0.5 + N[(N[(t * t), $MachinePrecision] * 0.125), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\left(\left(x \cdot 0.5 - y\right) \cdot \sqrt{z \cdot 2}\right) \cdot \left(1 + \left(t \cdot t\right) \cdot \left(0.5 + \left(t \cdot t\right) \cdot 0.125\right)\right)
\end{array}
Initial program 99.8%
Taylor expanded in t around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6491.5%
Simplified91.5%
(FPCore (x y z t)
:precision binary64
(let* ((t_1 (- (* x 0.5) y)))
(if (<= (* t t) 0.4)
(* t_1 (sqrt (* z (* 2.0 (+ (* t t) 1.0)))))
(* t_1 (sqrt (* t (* z (* t (* t t)))))))))
double code(double x, double y, double z, double t) {
double t_1 = (x * 0.5) - y;
double tmp;
if ((t * t) <= 0.4) {
tmp = t_1 * sqrt((z * (2.0 * ((t * t) + 1.0))));
} else {
tmp = t_1 * sqrt((t * (z * (t * (t * t)))));
}
return tmp;
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8) :: t_1
real(8) :: tmp
t_1 = (x * 0.5d0) - y
if ((t * t) <= 0.4d0) then
tmp = t_1 * sqrt((z * (2.0d0 * ((t * t) + 1.0d0))))
else
tmp = t_1 * sqrt((t * (z * (t * (t * t)))))
end if
code = tmp
end function
public static double code(double x, double y, double z, double t) {
double t_1 = (x * 0.5) - y;
double tmp;
if ((t * t) <= 0.4) {
tmp = t_1 * Math.sqrt((z * (2.0 * ((t * t) + 1.0))));
} else {
tmp = t_1 * Math.sqrt((t * (z * (t * (t * t)))));
}
return tmp;
}
def code(x, y, z, t): t_1 = (x * 0.5) - y tmp = 0 if (t * t) <= 0.4: tmp = t_1 * math.sqrt((z * (2.0 * ((t * t) + 1.0)))) else: tmp = t_1 * math.sqrt((t * (z * (t * (t * t))))) return tmp
function code(x, y, z, t) t_1 = Float64(Float64(x * 0.5) - y) tmp = 0.0 if (Float64(t * t) <= 0.4) tmp = Float64(t_1 * sqrt(Float64(z * Float64(2.0 * Float64(Float64(t * t) + 1.0))))); else tmp = Float64(t_1 * sqrt(Float64(t * Float64(z * Float64(t * Float64(t * t)))))); end return tmp end
function tmp_2 = code(x, y, z, t) t_1 = (x * 0.5) - y; tmp = 0.0; if ((t * t) <= 0.4) tmp = t_1 * sqrt((z * (2.0 * ((t * t) + 1.0)))); else tmp = t_1 * sqrt((t * (z * (t * (t * t))))); end tmp_2 = tmp; end
code[x_, y_, z_, t_] := Block[{t$95$1 = N[(N[(x * 0.5), $MachinePrecision] - y), $MachinePrecision]}, If[LessEqual[N[(t * t), $MachinePrecision], 0.4], N[(t$95$1 * N[Sqrt[N[(z * N[(2.0 * N[(N[(t * t), $MachinePrecision] + 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], N[(t$95$1 * N[Sqrt[N[(t * N[(z * N[(t * N[(t * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x \cdot 0.5 - y\\
\mathbf{if}\;t \cdot t \leq 0.4:\\
\;\;\;\;t\_1 \cdot \sqrt{z \cdot \left(2 \cdot \left(t \cdot t + 1\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;t\_1 \cdot \sqrt{t \cdot \left(z \cdot \left(t \cdot \left(t \cdot t\right)\right)\right)}\\
\end{array}
\end{array}
if (*.f64 t t) < 0.40000000000000002Initial program 99.6%
associate-*l*N/A
*-lowering-*.f64N/A
--lowering--.f64N/A
*-lowering-*.f64N/A
exp-sqrtN/A
sqrt-unprodN/A
sqrt-lowering-sqrt.f64N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
exp-lowering-exp.f64N/A
*-lowering-*.f6499.6%
Applied egg-rr99.6%
Taylor expanded in t around 0
+-commutativeN/A
+-lowering-+.f64N/A
unpow2N/A
*-lowering-*.f6498.5%
Simplified98.5%
if 0.40000000000000002 < (*.f64 t t) Initial program 100.0%
associate-*l*N/A
*-lowering-*.f64N/A
--lowering--.f64N/A
*-lowering-*.f64N/A
exp-sqrtN/A
sqrt-unprodN/A
sqrt-lowering-sqrt.f64N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
exp-lowering-exp.f64N/A
*-lowering-*.f64100.0%
Applied egg-rr100.0%
Taylor expanded in t around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
distribute-rgt-outN/A
associate-*r*N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
unpow2N/A
*-lowering-*.f6482.0%
Simplified82.0%
Taylor expanded in t around inf
metadata-evalN/A
pow-sqrN/A
associate-*r*N/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6482.0%
Simplified82.0%
associate-*r*N/A
cube-unmultN/A
*-lowering-*.f64N/A
cube-unmultN/A
*-lowering-*.f64N/A
*-lowering-*.f6482.8%
Applied egg-rr82.8%
Final simplification90.6%
(FPCore (x y z t)
:precision binary64
(let* ((t_1 (- (* x 0.5) y)))
(if (<= (* t t) 0.4)
(* t_1 (sqrt (* z 2.0)))
(* t_1 (sqrt (* t (* z (* t (* t t)))))))))
double code(double x, double y, double z, double t) {
double t_1 = (x * 0.5) - y;
double tmp;
if ((t * t) <= 0.4) {
tmp = t_1 * sqrt((z * 2.0));
} else {
tmp = t_1 * sqrt((t * (z * (t * (t * t)))));
}
return tmp;
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8) :: t_1
real(8) :: tmp
t_1 = (x * 0.5d0) - y
if ((t * t) <= 0.4d0) then
tmp = t_1 * sqrt((z * 2.0d0))
else
tmp = t_1 * sqrt((t * (z * (t * (t * t)))))
end if
code = tmp
end function
public static double code(double x, double y, double z, double t) {
double t_1 = (x * 0.5) - y;
double tmp;
if ((t * t) <= 0.4) {
tmp = t_1 * Math.sqrt((z * 2.0));
} else {
tmp = t_1 * Math.sqrt((t * (z * (t * (t * t)))));
}
return tmp;
}
def code(x, y, z, t): t_1 = (x * 0.5) - y tmp = 0 if (t * t) <= 0.4: tmp = t_1 * math.sqrt((z * 2.0)) else: tmp = t_1 * math.sqrt((t * (z * (t * (t * t))))) return tmp
function code(x, y, z, t) t_1 = Float64(Float64(x * 0.5) - y) tmp = 0.0 if (Float64(t * t) <= 0.4) tmp = Float64(t_1 * sqrt(Float64(z * 2.0))); else tmp = Float64(t_1 * sqrt(Float64(t * Float64(z * Float64(t * Float64(t * t)))))); end return tmp end
function tmp_2 = code(x, y, z, t) t_1 = (x * 0.5) - y; tmp = 0.0; if ((t * t) <= 0.4) tmp = t_1 * sqrt((z * 2.0)); else tmp = t_1 * sqrt((t * (z * (t * (t * t))))); end tmp_2 = tmp; end
code[x_, y_, z_, t_] := Block[{t$95$1 = N[(N[(x * 0.5), $MachinePrecision] - y), $MachinePrecision]}, If[LessEqual[N[(t * t), $MachinePrecision], 0.4], N[(t$95$1 * N[Sqrt[N[(z * 2.0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], N[(t$95$1 * N[Sqrt[N[(t * N[(z * N[(t * N[(t * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x \cdot 0.5 - y\\
\mathbf{if}\;t \cdot t \leq 0.4:\\
\;\;\;\;t\_1 \cdot \sqrt{z \cdot 2}\\
\mathbf{else}:\\
\;\;\;\;t\_1 \cdot \sqrt{t \cdot \left(z \cdot \left(t \cdot \left(t \cdot t\right)\right)\right)}\\
\end{array}
\end{array}
if (*.f64 t t) < 0.40000000000000002Initial program 99.6%
Taylor expanded in t around 0
Simplified98.4%
*-rgt-identityN/A
*-lowering-*.f64N/A
--lowering--.f64N/A
*-lowering-*.f64N/A
sqrt-lowering-sqrt.f64N/A
*-lowering-*.f6498.4%
Applied egg-rr98.4%
if 0.40000000000000002 < (*.f64 t t) Initial program 100.0%
associate-*l*N/A
*-lowering-*.f64N/A
--lowering--.f64N/A
*-lowering-*.f64N/A
exp-sqrtN/A
sqrt-unprodN/A
sqrt-lowering-sqrt.f64N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
exp-lowering-exp.f64N/A
*-lowering-*.f64100.0%
Applied egg-rr100.0%
Taylor expanded in t around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
distribute-rgt-outN/A
associate-*r*N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
unpow2N/A
*-lowering-*.f6482.0%
Simplified82.0%
Taylor expanded in t around inf
metadata-evalN/A
pow-sqrN/A
associate-*r*N/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6482.0%
Simplified82.0%
associate-*r*N/A
cube-unmultN/A
*-lowering-*.f64N/A
cube-unmultN/A
*-lowering-*.f64N/A
*-lowering-*.f6482.8%
Applied egg-rr82.8%
Final simplification90.6%
(FPCore (x y z t)
:precision binary64
(let* ((t_1 (- (* x 0.5) y)))
(if (<= (* t t) 0.4)
(* t_1 (sqrt (* z 2.0)))
(* t_1 (sqrt (* t (* t (* z (* t t)))))))))
double code(double x, double y, double z, double t) {
double t_1 = (x * 0.5) - y;
double tmp;
if ((t * t) <= 0.4) {
tmp = t_1 * sqrt((z * 2.0));
} else {
tmp = t_1 * sqrt((t * (t * (z * (t * t)))));
}
return tmp;
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8) :: t_1
real(8) :: tmp
t_1 = (x * 0.5d0) - y
if ((t * t) <= 0.4d0) then
tmp = t_1 * sqrt((z * 2.0d0))
else
tmp = t_1 * sqrt((t * (t * (z * (t * t)))))
end if
code = tmp
end function
public static double code(double x, double y, double z, double t) {
double t_1 = (x * 0.5) - y;
double tmp;
if ((t * t) <= 0.4) {
tmp = t_1 * Math.sqrt((z * 2.0));
} else {
tmp = t_1 * Math.sqrt((t * (t * (z * (t * t)))));
}
return tmp;
}
def code(x, y, z, t): t_1 = (x * 0.5) - y tmp = 0 if (t * t) <= 0.4: tmp = t_1 * math.sqrt((z * 2.0)) else: tmp = t_1 * math.sqrt((t * (t * (z * (t * t))))) return tmp
function code(x, y, z, t) t_1 = Float64(Float64(x * 0.5) - y) tmp = 0.0 if (Float64(t * t) <= 0.4) tmp = Float64(t_1 * sqrt(Float64(z * 2.0))); else tmp = Float64(t_1 * sqrt(Float64(t * Float64(t * Float64(z * Float64(t * t)))))); end return tmp end
function tmp_2 = code(x, y, z, t) t_1 = (x * 0.5) - y; tmp = 0.0; if ((t * t) <= 0.4) tmp = t_1 * sqrt((z * 2.0)); else tmp = t_1 * sqrt((t * (t * (z * (t * t))))); end tmp_2 = tmp; end
code[x_, y_, z_, t_] := Block[{t$95$1 = N[(N[(x * 0.5), $MachinePrecision] - y), $MachinePrecision]}, If[LessEqual[N[(t * t), $MachinePrecision], 0.4], N[(t$95$1 * N[Sqrt[N[(z * 2.0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], N[(t$95$1 * N[Sqrt[N[(t * N[(t * N[(z * N[(t * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x \cdot 0.5 - y\\
\mathbf{if}\;t \cdot t \leq 0.4:\\
\;\;\;\;t\_1 \cdot \sqrt{z \cdot 2}\\
\mathbf{else}:\\
\;\;\;\;t\_1 \cdot \sqrt{t \cdot \left(t \cdot \left(z \cdot \left(t \cdot t\right)\right)\right)}\\
\end{array}
\end{array}
if (*.f64 t t) < 0.40000000000000002Initial program 99.6%
Taylor expanded in t around 0
Simplified98.4%
*-rgt-identityN/A
*-lowering-*.f64N/A
--lowering--.f64N/A
*-lowering-*.f64N/A
sqrt-lowering-sqrt.f64N/A
*-lowering-*.f6498.4%
Applied egg-rr98.4%
if 0.40000000000000002 < (*.f64 t t) Initial program 100.0%
associate-*l*N/A
*-lowering-*.f64N/A
--lowering--.f64N/A
*-lowering-*.f64N/A
exp-sqrtN/A
sqrt-unprodN/A
sqrt-lowering-sqrt.f64N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
exp-lowering-exp.f64N/A
*-lowering-*.f64100.0%
Applied egg-rr100.0%
Taylor expanded in t around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
distribute-rgt-outN/A
associate-*r*N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
unpow2N/A
*-lowering-*.f6482.0%
Simplified82.0%
Taylor expanded in t around inf
metadata-evalN/A
pow-sqrN/A
associate-*r*N/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6482.0%
Simplified82.0%
Final simplification90.2%
(FPCore (x y z t)
:precision binary64
(if (<= t 9e+95)
(* (- (* x 0.5) y) (sqrt (* z 2.0)))
(if (<= t 2e+262)
(* (sqrt z) (* x (* 0.5 (* t t))))
(* (sqrt z) (- 0.0 (* y (* t t)))))))
double code(double x, double y, double z, double t) {
double tmp;
if (t <= 9e+95) {
tmp = ((x * 0.5) - y) * sqrt((z * 2.0));
} else if (t <= 2e+262) {
tmp = sqrt(z) * (x * (0.5 * (t * t)));
} else {
tmp = sqrt(z) * (0.0 - (y * (t * t)));
}
return tmp;
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8) :: tmp
if (t <= 9d+95) then
tmp = ((x * 0.5d0) - y) * sqrt((z * 2.0d0))
else if (t <= 2d+262) then
tmp = sqrt(z) * (x * (0.5d0 * (t * t)))
else
tmp = sqrt(z) * (0.0d0 - (y * (t * t)))
end if
code = tmp
end function
public static double code(double x, double y, double z, double t) {
double tmp;
if (t <= 9e+95) {
tmp = ((x * 0.5) - y) * Math.sqrt((z * 2.0));
} else if (t <= 2e+262) {
tmp = Math.sqrt(z) * (x * (0.5 * (t * t)));
} else {
tmp = Math.sqrt(z) * (0.0 - (y * (t * t)));
}
return tmp;
}
def code(x, y, z, t): tmp = 0 if t <= 9e+95: tmp = ((x * 0.5) - y) * math.sqrt((z * 2.0)) elif t <= 2e+262: tmp = math.sqrt(z) * (x * (0.5 * (t * t))) else: tmp = math.sqrt(z) * (0.0 - (y * (t * t))) return tmp
function code(x, y, z, t) tmp = 0.0 if (t <= 9e+95) tmp = Float64(Float64(Float64(x * 0.5) - y) * sqrt(Float64(z * 2.0))); elseif (t <= 2e+262) tmp = Float64(sqrt(z) * Float64(x * Float64(0.5 * Float64(t * t)))); else tmp = Float64(sqrt(z) * Float64(0.0 - Float64(y * Float64(t * t)))); end return tmp end
function tmp_2 = code(x, y, z, t) tmp = 0.0; if (t <= 9e+95) tmp = ((x * 0.5) - y) * sqrt((z * 2.0)); elseif (t <= 2e+262) tmp = sqrt(z) * (x * (0.5 * (t * t))); else tmp = sqrt(z) * (0.0 - (y * (t * t))); end tmp_2 = tmp; end
code[x_, y_, z_, t_] := If[LessEqual[t, 9e+95], N[(N[(N[(x * 0.5), $MachinePrecision] - y), $MachinePrecision] * N[Sqrt[N[(z * 2.0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], If[LessEqual[t, 2e+262], N[(N[Sqrt[z], $MachinePrecision] * N[(x * N[(0.5 * N[(t * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[Sqrt[z], $MachinePrecision] * N[(0.0 - N[(y * N[(t * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;t \leq 9 \cdot 10^{+95}:\\
\;\;\;\;\left(x \cdot 0.5 - y\right) \cdot \sqrt{z \cdot 2}\\
\mathbf{elif}\;t \leq 2 \cdot 10^{+262}:\\
\;\;\;\;\sqrt{z} \cdot \left(x \cdot \left(0.5 \cdot \left(t \cdot t\right)\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\sqrt{z} \cdot \left(0 - y \cdot \left(t \cdot t\right)\right)\\
\end{array}
\end{array}
if t < 9.00000000000000033e95Initial program 99.8%
Taylor expanded in t around 0
Simplified66.2%
*-rgt-identityN/A
*-lowering-*.f64N/A
--lowering--.f64N/A
*-lowering-*.f64N/A
sqrt-lowering-sqrt.f64N/A
*-lowering-*.f6466.2%
Applied egg-rr66.2%
if 9.00000000000000033e95 < t < 2e262Initial program 100.0%
associate-*l*N/A
*-lowering-*.f64N/A
--lowering--.f64N/A
*-lowering-*.f64N/A
exp-sqrtN/A
sqrt-unprodN/A
sqrt-lowering-sqrt.f64N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
exp-lowering-exp.f64N/A
*-lowering-*.f64100.0%
Applied egg-rr100.0%
Taylor expanded in t around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
distribute-rgt-outN/A
associate-*r*N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
unpow2N/A
*-lowering-*.f6497.3%
Simplified97.3%
Taylor expanded in t around inf
metadata-evalN/A
pow-sqrN/A
associate-*r*N/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6497.3%
Simplified97.3%
Taylor expanded in x around inf
*-commutativeN/A
*-commutativeN/A
associate-*l*N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
sqrt-lowering-sqrt.f64N/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6471.8%
Simplified71.8%
if 2e262 < t Initial program 100.0%
associate-*l*N/A
*-lowering-*.f64N/A
--lowering--.f64N/A
*-lowering-*.f64N/A
exp-sqrtN/A
sqrt-unprodN/A
sqrt-lowering-sqrt.f64N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
exp-lowering-exp.f64N/A
*-lowering-*.f64100.0%
Applied egg-rr100.0%
Taylor expanded in t around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
distribute-rgt-outN/A
associate-*r*N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
unpow2N/A
*-lowering-*.f64100.0%
Simplified100.0%
Taylor expanded in t around inf
metadata-evalN/A
pow-sqrN/A
associate-*r*N/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64100.0%
Simplified100.0%
Taylor expanded in x around 0
mul-1-negN/A
*-commutativeN/A
distribute-rgt-neg-inN/A
mul-1-negN/A
*-lowering-*.f64N/A
sqrt-lowering-sqrt.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64100.0%
Simplified100.0%
Final simplification68.0%
(FPCore (x y z t)
:precision binary64
(let* ((t_1 (- (* x 0.5) y)))
(if (<= (* t t) 0.4)
(* t_1 (sqrt (* z 2.0)))
(* (sqrt z) (* t_1 (* t t))))))
double code(double x, double y, double z, double t) {
double t_1 = (x * 0.5) - y;
double tmp;
if ((t * t) <= 0.4) {
tmp = t_1 * sqrt((z * 2.0));
} else {
tmp = sqrt(z) * (t_1 * (t * t));
}
return tmp;
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8) :: t_1
real(8) :: tmp
t_1 = (x * 0.5d0) - y
if ((t * t) <= 0.4d0) then
tmp = t_1 * sqrt((z * 2.0d0))
else
tmp = sqrt(z) * (t_1 * (t * t))
end if
code = tmp
end function
public static double code(double x, double y, double z, double t) {
double t_1 = (x * 0.5) - y;
double tmp;
if ((t * t) <= 0.4) {
tmp = t_1 * Math.sqrt((z * 2.0));
} else {
tmp = Math.sqrt(z) * (t_1 * (t * t));
}
return tmp;
}
def code(x, y, z, t): t_1 = (x * 0.5) - y tmp = 0 if (t * t) <= 0.4: tmp = t_1 * math.sqrt((z * 2.0)) else: tmp = math.sqrt(z) * (t_1 * (t * t)) return tmp
function code(x, y, z, t) t_1 = Float64(Float64(x * 0.5) - y) tmp = 0.0 if (Float64(t * t) <= 0.4) tmp = Float64(t_1 * sqrt(Float64(z * 2.0))); else tmp = Float64(sqrt(z) * Float64(t_1 * Float64(t * t))); end return tmp end
function tmp_2 = code(x, y, z, t) t_1 = (x * 0.5) - y; tmp = 0.0; if ((t * t) <= 0.4) tmp = t_1 * sqrt((z * 2.0)); else tmp = sqrt(z) * (t_1 * (t * t)); end tmp_2 = tmp; end
code[x_, y_, z_, t_] := Block[{t$95$1 = N[(N[(x * 0.5), $MachinePrecision] - y), $MachinePrecision]}, If[LessEqual[N[(t * t), $MachinePrecision], 0.4], N[(t$95$1 * N[Sqrt[N[(z * 2.0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], N[(N[Sqrt[z], $MachinePrecision] * N[(t$95$1 * N[(t * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x \cdot 0.5 - y\\
\mathbf{if}\;t \cdot t \leq 0.4:\\
\;\;\;\;t\_1 \cdot \sqrt{z \cdot 2}\\
\mathbf{else}:\\
\;\;\;\;\sqrt{z} \cdot \left(t\_1 \cdot \left(t \cdot t\right)\right)\\
\end{array}
\end{array}
if (*.f64 t t) < 0.40000000000000002Initial program 99.6%
Taylor expanded in t around 0
Simplified98.4%
*-rgt-identityN/A
*-lowering-*.f64N/A
--lowering--.f64N/A
*-lowering-*.f64N/A
sqrt-lowering-sqrt.f64N/A
*-lowering-*.f6498.4%
Applied egg-rr98.4%
if 0.40000000000000002 < (*.f64 t t) Initial program 100.0%
associate-*l*N/A
*-lowering-*.f64N/A
--lowering--.f64N/A
*-lowering-*.f64N/A
exp-sqrtN/A
sqrt-unprodN/A
sqrt-lowering-sqrt.f64N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
exp-lowering-exp.f64N/A
*-lowering-*.f64100.0%
Applied egg-rr100.0%
Taylor expanded in t around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
distribute-rgt-outN/A
associate-*r*N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
unpow2N/A
*-lowering-*.f6482.0%
Simplified82.0%
Taylor expanded in t around inf
metadata-evalN/A
pow-sqrN/A
associate-*r*N/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6482.0%
Simplified82.0%
Taylor expanded in z around -inf
mul-1-negN/A
*-commutativeN/A
distribute-rgt-neg-inN/A
mul-1-negN/A
*-lowering-*.f64N/A
sqrt-lowering-sqrt.f64N/A
associate-*r*N/A
associate-*r*N/A
*-lowering-*.f64N/A
Simplified77.0%
Final simplification87.7%
(FPCore (x y z t)
:precision binary64
(let* ((t_1 (- (* x 0.5) y)))
(if (<= (* t t) 0.4)
(* t_1 (sqrt (* z 2.0)))
(* t_1 (* (* t t) (sqrt z))))))
double code(double x, double y, double z, double t) {
double t_1 = (x * 0.5) - y;
double tmp;
if ((t * t) <= 0.4) {
tmp = t_1 * sqrt((z * 2.0));
} else {
tmp = t_1 * ((t * t) * sqrt(z));
}
return tmp;
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8) :: t_1
real(8) :: tmp
t_1 = (x * 0.5d0) - y
if ((t * t) <= 0.4d0) then
tmp = t_1 * sqrt((z * 2.0d0))
else
tmp = t_1 * ((t * t) * sqrt(z))
end if
code = tmp
end function
public static double code(double x, double y, double z, double t) {
double t_1 = (x * 0.5) - y;
double tmp;
if ((t * t) <= 0.4) {
tmp = t_1 * Math.sqrt((z * 2.0));
} else {
tmp = t_1 * ((t * t) * Math.sqrt(z));
}
return tmp;
}
def code(x, y, z, t): t_1 = (x * 0.5) - y tmp = 0 if (t * t) <= 0.4: tmp = t_1 * math.sqrt((z * 2.0)) else: tmp = t_1 * ((t * t) * math.sqrt(z)) return tmp
function code(x, y, z, t) t_1 = Float64(Float64(x * 0.5) - y) tmp = 0.0 if (Float64(t * t) <= 0.4) tmp = Float64(t_1 * sqrt(Float64(z * 2.0))); else tmp = Float64(t_1 * Float64(Float64(t * t) * sqrt(z))); end return tmp end
function tmp_2 = code(x, y, z, t) t_1 = (x * 0.5) - y; tmp = 0.0; if ((t * t) <= 0.4) tmp = t_1 * sqrt((z * 2.0)); else tmp = t_1 * ((t * t) * sqrt(z)); end tmp_2 = tmp; end
code[x_, y_, z_, t_] := Block[{t$95$1 = N[(N[(x * 0.5), $MachinePrecision] - y), $MachinePrecision]}, If[LessEqual[N[(t * t), $MachinePrecision], 0.4], N[(t$95$1 * N[Sqrt[N[(z * 2.0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], N[(t$95$1 * N[(N[(t * t), $MachinePrecision] * N[Sqrt[z], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x \cdot 0.5 - y\\
\mathbf{if}\;t \cdot t \leq 0.4:\\
\;\;\;\;t\_1 \cdot \sqrt{z \cdot 2}\\
\mathbf{else}:\\
\;\;\;\;t\_1 \cdot \left(\left(t \cdot t\right) \cdot \sqrt{z}\right)\\
\end{array}
\end{array}
if (*.f64 t t) < 0.40000000000000002Initial program 99.6%
Taylor expanded in t around 0
Simplified98.4%
*-rgt-identityN/A
*-lowering-*.f64N/A
--lowering--.f64N/A
*-lowering-*.f64N/A
sqrt-lowering-sqrt.f64N/A
*-lowering-*.f6498.4%
Applied egg-rr98.4%
if 0.40000000000000002 < (*.f64 t t) Initial program 100.0%
associate-*l*N/A
*-lowering-*.f64N/A
--lowering--.f64N/A
*-lowering-*.f64N/A
exp-sqrtN/A
sqrt-unprodN/A
sqrt-lowering-sqrt.f64N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
exp-lowering-exp.f64N/A
*-lowering-*.f64100.0%
Applied egg-rr100.0%
Taylor expanded in t around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
distribute-rgt-outN/A
associate-*r*N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
unpow2N/A
*-lowering-*.f6482.0%
Simplified82.0%
Taylor expanded in t around inf
metadata-evalN/A
pow-sqrN/A
associate-*r*N/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6482.0%
Simplified82.0%
Taylor expanded in t around 0
*-commutativeN/A
*-lowering-*.f64N/A
sqrt-lowering-sqrt.f64N/A
unpow2N/A
*-lowering-*.f6476.3%
Simplified76.3%
Final simplification87.3%
(FPCore (x y z t)
:precision binary64
(let* ((t_1 (- (* x 0.5) y)))
(if (<= (* t t) 0.4)
(* t_1 (sqrt (* z 2.0)))
(* (* t t) (* t_1 (sqrt z))))))
double code(double x, double y, double z, double t) {
double t_1 = (x * 0.5) - y;
double tmp;
if ((t * t) <= 0.4) {
tmp = t_1 * sqrt((z * 2.0));
} else {
tmp = (t * t) * (t_1 * sqrt(z));
}
return tmp;
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8) :: t_1
real(8) :: tmp
t_1 = (x * 0.5d0) - y
if ((t * t) <= 0.4d0) then
tmp = t_1 * sqrt((z * 2.0d0))
else
tmp = (t * t) * (t_1 * sqrt(z))
end if
code = tmp
end function
public static double code(double x, double y, double z, double t) {
double t_1 = (x * 0.5) - y;
double tmp;
if ((t * t) <= 0.4) {
tmp = t_1 * Math.sqrt((z * 2.0));
} else {
tmp = (t * t) * (t_1 * Math.sqrt(z));
}
return tmp;
}
def code(x, y, z, t): t_1 = (x * 0.5) - y tmp = 0 if (t * t) <= 0.4: tmp = t_1 * math.sqrt((z * 2.0)) else: tmp = (t * t) * (t_1 * math.sqrt(z)) return tmp
function code(x, y, z, t) t_1 = Float64(Float64(x * 0.5) - y) tmp = 0.0 if (Float64(t * t) <= 0.4) tmp = Float64(t_1 * sqrt(Float64(z * 2.0))); else tmp = Float64(Float64(t * t) * Float64(t_1 * sqrt(z))); end return tmp end
function tmp_2 = code(x, y, z, t) t_1 = (x * 0.5) - y; tmp = 0.0; if ((t * t) <= 0.4) tmp = t_1 * sqrt((z * 2.0)); else tmp = (t * t) * (t_1 * sqrt(z)); end tmp_2 = tmp; end
code[x_, y_, z_, t_] := Block[{t$95$1 = N[(N[(x * 0.5), $MachinePrecision] - y), $MachinePrecision]}, If[LessEqual[N[(t * t), $MachinePrecision], 0.4], N[(t$95$1 * N[Sqrt[N[(z * 2.0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision], N[(N[(t * t), $MachinePrecision] * N[(t$95$1 * N[Sqrt[z], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_1 := x \cdot 0.5 - y\\
\mathbf{if}\;t \cdot t \leq 0.4:\\
\;\;\;\;t\_1 \cdot \sqrt{z \cdot 2}\\
\mathbf{else}:\\
\;\;\;\;\left(t \cdot t\right) \cdot \left(t\_1 \cdot \sqrt{z}\right)\\
\end{array}
\end{array}
if (*.f64 t t) < 0.40000000000000002Initial program 99.6%
Taylor expanded in t around 0
Simplified98.4%
*-rgt-identityN/A
*-lowering-*.f64N/A
--lowering--.f64N/A
*-lowering-*.f64N/A
sqrt-lowering-sqrt.f64N/A
*-lowering-*.f6498.4%
Applied egg-rr98.4%
if 0.40000000000000002 < (*.f64 t t) Initial program 100.0%
associate-*l*N/A
*-lowering-*.f64N/A
--lowering--.f64N/A
*-lowering-*.f64N/A
exp-sqrtN/A
sqrt-unprodN/A
sqrt-lowering-sqrt.f64N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
exp-lowering-exp.f64N/A
*-lowering-*.f64100.0%
Applied egg-rr100.0%
Taylor expanded in t around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
distribute-rgt-outN/A
associate-*r*N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
unpow2N/A
*-lowering-*.f6482.0%
Simplified82.0%
Taylor expanded in t around inf
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
sqrt-lowering-sqrt.f64N/A
--lowering--.f64N/A
*-lowering-*.f6474.2%
Simplified74.2%
Final simplification86.3%
(FPCore (x y z t) :precision binary64 (if (<= t 150000000000.0) (/ (* x 0.5) (pow (* z 2.0) -0.5)) (* (sqrt z) (- 0.0 (* y (* t t))))))
double code(double x, double y, double z, double t) {
double tmp;
if (t <= 150000000000.0) {
tmp = (x * 0.5) / pow((z * 2.0), -0.5);
} else {
tmp = sqrt(z) * (0.0 - (y * (t * t)));
}
return tmp;
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
real(8) :: tmp
if (t <= 150000000000.0d0) then
tmp = (x * 0.5d0) / ((z * 2.0d0) ** (-0.5d0))
else
tmp = sqrt(z) * (0.0d0 - (y * (t * t)))
end if
code = tmp
end function
public static double code(double x, double y, double z, double t) {
double tmp;
if (t <= 150000000000.0) {
tmp = (x * 0.5) / Math.pow((z * 2.0), -0.5);
} else {
tmp = Math.sqrt(z) * (0.0 - (y * (t * t)));
}
return tmp;
}
def code(x, y, z, t): tmp = 0 if t <= 150000000000.0: tmp = (x * 0.5) / math.pow((z * 2.0), -0.5) else: tmp = math.sqrt(z) * (0.0 - (y * (t * t))) return tmp
function code(x, y, z, t) tmp = 0.0 if (t <= 150000000000.0) tmp = Float64(Float64(x * 0.5) / (Float64(z * 2.0) ^ -0.5)); else tmp = Float64(sqrt(z) * Float64(0.0 - Float64(y * Float64(t * t)))); end return tmp end
function tmp_2 = code(x, y, z, t) tmp = 0.0; if (t <= 150000000000.0) tmp = (x * 0.5) / ((z * 2.0) ^ -0.5); else tmp = sqrt(z) * (0.0 - (y * (t * t))); end tmp_2 = tmp; end
code[x_, y_, z_, t_] := If[LessEqual[t, 150000000000.0], N[(N[(x * 0.5), $MachinePrecision] / N[Power[N[(z * 2.0), $MachinePrecision], -0.5], $MachinePrecision]), $MachinePrecision], N[(N[Sqrt[z], $MachinePrecision] * N[(0.0 - N[(y * N[(t * t), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;t \leq 150000000000:\\
\;\;\;\;\frac{x \cdot 0.5}{{\left(z \cdot 2\right)}^{-0.5}}\\
\mathbf{else}:\\
\;\;\;\;\sqrt{z} \cdot \left(0 - y \cdot \left(t \cdot t\right)\right)\\
\end{array}
\end{array}
if t < 1.5e11Initial program 99.7%
Taylor expanded in t around 0
Simplified70.6%
Taylor expanded in x around inf
*-lowering-*.f6437.0%
Simplified37.0%
*-rgt-identityN/A
*-commutativeN/A
*-commutativeN/A
*-commutativeN/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
sqrt-lowering-sqrt.f64N/A
*-commutativeN/A
*-lowering-*.f6437.0%
Applied egg-rr37.0%
*-commutativeN/A
*-commutativeN/A
associate-*r*N/A
pow1/2N/A
metadata-evalN/A
pow-flipN/A
un-div-invN/A
/-lowering-/.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
pow-lowering-pow.f64N/A
*-lowering-*.f6437.0%
Applied egg-rr37.0%
if 1.5e11 < t Initial program 100.0%
associate-*l*N/A
*-lowering-*.f64N/A
--lowering--.f64N/A
*-lowering-*.f64N/A
exp-sqrtN/A
sqrt-unprodN/A
sqrt-lowering-sqrt.f64N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
exp-lowering-exp.f64N/A
*-lowering-*.f64100.0%
Applied egg-rr100.0%
Taylor expanded in t around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
distribute-rgt-outN/A
associate-*r*N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
unpow2N/A
*-lowering-*.f6480.8%
Simplified80.8%
Taylor expanded in t around inf
metadata-evalN/A
pow-sqrN/A
associate-*r*N/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6480.8%
Simplified80.8%
Taylor expanded in x around 0
mul-1-negN/A
*-commutativeN/A
distribute-rgt-neg-inN/A
mul-1-negN/A
*-lowering-*.f64N/A
sqrt-lowering-sqrt.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6461.0%
Simplified61.0%
Final simplification42.6%
(FPCore (x y z t) :precision binary64 (/ (* x 0.5) (pow (* z 2.0) -0.5)))
double code(double x, double y, double z, double t) {
return (x * 0.5) / pow((z * 2.0), -0.5);
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
code = (x * 0.5d0) / ((z * 2.0d0) ** (-0.5d0))
end function
public static double code(double x, double y, double z, double t) {
return (x * 0.5) / Math.pow((z * 2.0), -0.5);
}
def code(x, y, z, t): return (x * 0.5) / math.pow((z * 2.0), -0.5)
function code(x, y, z, t) return Float64(Float64(x * 0.5) / (Float64(z * 2.0) ^ -0.5)) end
function tmp = code(x, y, z, t) tmp = (x * 0.5) / ((z * 2.0) ^ -0.5); end
code[x_, y_, z_, t_] := N[(N[(x * 0.5), $MachinePrecision] / N[Power[N[(z * 2.0), $MachinePrecision], -0.5], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{x \cdot 0.5}{{\left(z \cdot 2\right)}^{-0.5}}
\end{array}
Initial program 99.8%
Taylor expanded in t around 0
Simplified57.3%
Taylor expanded in x around inf
*-lowering-*.f6430.8%
Simplified30.8%
*-rgt-identityN/A
*-commutativeN/A
*-commutativeN/A
*-commutativeN/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
sqrt-lowering-sqrt.f64N/A
*-commutativeN/A
*-lowering-*.f6430.8%
Applied egg-rr30.8%
*-commutativeN/A
*-commutativeN/A
associate-*r*N/A
pow1/2N/A
metadata-evalN/A
pow-flipN/A
un-div-invN/A
/-lowering-/.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
pow-lowering-pow.f64N/A
*-lowering-*.f6430.9%
Applied egg-rr30.9%
Final simplification30.9%
(FPCore (x y z t) :precision binary64 (* x (* 0.5 (sqrt (* z 2.0)))))
double code(double x, double y, double z, double t) {
return x * (0.5 * sqrt((z * 2.0)));
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
code = x * (0.5d0 * sqrt((z * 2.0d0)))
end function
public static double code(double x, double y, double z, double t) {
return x * (0.5 * Math.sqrt((z * 2.0)));
}
def code(x, y, z, t): return x * (0.5 * math.sqrt((z * 2.0)))
function code(x, y, z, t) return Float64(x * Float64(0.5 * sqrt(Float64(z * 2.0)))) end
function tmp = code(x, y, z, t) tmp = x * (0.5 * sqrt((z * 2.0))); end
code[x_, y_, z_, t_] := N[(x * N[(0.5 * N[Sqrt[N[(z * 2.0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
x \cdot \left(0.5 \cdot \sqrt{z \cdot 2}\right)
\end{array}
Initial program 99.8%
Taylor expanded in t around 0
Simplified57.3%
Taylor expanded in x around inf
*-lowering-*.f6430.8%
Simplified30.8%
*-rgt-identityN/A
*-commutativeN/A
*-commutativeN/A
*-commutativeN/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
sqrt-lowering-sqrt.f64N/A
*-commutativeN/A
*-lowering-*.f6430.8%
Applied egg-rr30.8%
Final simplification30.8%
(FPCore (x y z t) :precision binary64 (* (* (- (* x 0.5) y) (sqrt (* z 2.0))) (pow (exp 1.0) (/ (* t t) 2.0))))
double code(double x, double y, double z, double t) {
return (((x * 0.5) - y) * sqrt((z * 2.0))) * pow(exp(1.0), ((t * t) / 2.0));
}
real(8) function code(x, y, z, t)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8), intent (in) :: t
code = (((x * 0.5d0) - y) * sqrt((z * 2.0d0))) * (exp(1.0d0) ** ((t * t) / 2.0d0))
end function
public static double code(double x, double y, double z, double t) {
return (((x * 0.5) - y) * Math.sqrt((z * 2.0))) * Math.pow(Math.exp(1.0), ((t * t) / 2.0));
}
def code(x, y, z, t): return (((x * 0.5) - y) * math.sqrt((z * 2.0))) * math.pow(math.exp(1.0), ((t * t) / 2.0))
function code(x, y, z, t) return Float64(Float64(Float64(Float64(x * 0.5) - y) * sqrt(Float64(z * 2.0))) * (exp(1.0) ^ Float64(Float64(t * t) / 2.0))) end
function tmp = code(x, y, z, t) tmp = (((x * 0.5) - y) * sqrt((z * 2.0))) * (exp(1.0) ^ ((t * t) / 2.0)); end
code[x_, y_, z_, t_] := N[(N[(N[(N[(x * 0.5), $MachinePrecision] - y), $MachinePrecision] * N[Sqrt[N[(z * 2.0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] * N[Power[N[Exp[1.0], $MachinePrecision], N[(N[(t * t), $MachinePrecision] / 2.0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\left(\left(x \cdot 0.5 - y\right) \cdot \sqrt{z \cdot 2}\right) \cdot {\left(e^{1}\right)}^{\left(\frac{t \cdot t}{2}\right)}
\end{array}
herbie shell --seed 2024161
(FPCore (x y z t)
:name "Data.Number.Erf:$cinvnormcdf from erf-2.0.0.0, A"
:precision binary64
:alt
(! :herbie-platform default (* (* (- (* x 1/2) y) (sqrt (* z 2))) (pow (exp 1) (/ (* t t) 2))))
(* (* (- (* x 0.5) y) (sqrt (* z 2.0))) (exp (/ (* t t) 2.0))))