
(FPCore (x y z)
:precision binary64
(/
(*
(- x 2.0)
(+
(*
(+ (* (+ (* (+ (* x 4.16438922228) 78.6994924154) x) 137.519416416) x) y)
x)
z))
(+
(* (+ (* (+ (* (+ x 43.3400022514) x) 263.505074721) x) 313.399215894) x)
47.066876606)))
double code(double x, double y, double z) {
return ((x - 2.0) * ((((((((x * 4.16438922228) + 78.6994924154) * x) + 137.519416416) * x) + y) * x) + z)) / (((((((x + 43.3400022514) * x) + 263.505074721) * x) + 313.399215894) * x) + 47.066876606);
}
real(8) function code(x, y, z)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
code = ((x - 2.0d0) * ((((((((x * 4.16438922228d0) + 78.6994924154d0) * x) + 137.519416416d0) * x) + y) * x) + z)) / (((((((x + 43.3400022514d0) * x) + 263.505074721d0) * x) + 313.399215894d0) * x) + 47.066876606d0)
end function
public static double code(double x, double y, double z) {
return ((x - 2.0) * ((((((((x * 4.16438922228) + 78.6994924154) * x) + 137.519416416) * x) + y) * x) + z)) / (((((((x + 43.3400022514) * x) + 263.505074721) * x) + 313.399215894) * x) + 47.066876606);
}
def code(x, y, z): return ((x - 2.0) * ((((((((x * 4.16438922228) + 78.6994924154) * x) + 137.519416416) * x) + y) * x) + z)) / (((((((x + 43.3400022514) * x) + 263.505074721) * x) + 313.399215894) * x) + 47.066876606)
function code(x, y, z) return Float64(Float64(Float64(x - 2.0) * Float64(Float64(Float64(Float64(Float64(Float64(Float64(Float64(x * 4.16438922228) + 78.6994924154) * x) + 137.519416416) * x) + y) * x) + z)) / Float64(Float64(Float64(Float64(Float64(Float64(Float64(x + 43.3400022514) * x) + 263.505074721) * x) + 313.399215894) * x) + 47.066876606)) end
function tmp = code(x, y, z) tmp = ((x - 2.0) * ((((((((x * 4.16438922228) + 78.6994924154) * x) + 137.519416416) * x) + y) * x) + z)) / (((((((x + 43.3400022514) * x) + 263.505074721) * x) + 313.399215894) * x) + 47.066876606); end
code[x_, y_, z_] := N[(N[(N[(x - 2.0), $MachinePrecision] * N[(N[(N[(N[(N[(N[(N[(N[(x * 4.16438922228), $MachinePrecision] + 78.6994924154), $MachinePrecision] * x), $MachinePrecision] + 137.519416416), $MachinePrecision] * x), $MachinePrecision] + y), $MachinePrecision] * x), $MachinePrecision] + z), $MachinePrecision]), $MachinePrecision] / N[(N[(N[(N[(N[(N[(N[(x + 43.3400022514), $MachinePrecision] * x), $MachinePrecision] + 263.505074721), $MachinePrecision] * x), $MachinePrecision] + 313.399215894), $MachinePrecision] * x), $MachinePrecision] + 47.066876606), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{\left(x - 2\right) \cdot \left(\left(\left(\left(x \cdot 4.16438922228 + 78.6994924154\right) \cdot x + 137.519416416\right) \cdot x + y\right) \cdot x + z\right)}{\left(\left(\left(x + 43.3400022514\right) \cdot x + 263.505074721\right) \cdot x + 313.399215894\right) \cdot x + 47.066876606}
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 23 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (x y z)
:precision binary64
(/
(*
(- x 2.0)
(+
(*
(+ (* (+ (* (+ (* x 4.16438922228) 78.6994924154) x) 137.519416416) x) y)
x)
z))
(+
(* (+ (* (+ (* (+ x 43.3400022514) x) 263.505074721) x) 313.399215894) x)
47.066876606)))
double code(double x, double y, double z) {
return ((x - 2.0) * ((((((((x * 4.16438922228) + 78.6994924154) * x) + 137.519416416) * x) + y) * x) + z)) / (((((((x + 43.3400022514) * x) + 263.505074721) * x) + 313.399215894) * x) + 47.066876606);
}
real(8) function code(x, y, z)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
code = ((x - 2.0d0) * ((((((((x * 4.16438922228d0) + 78.6994924154d0) * x) + 137.519416416d0) * x) + y) * x) + z)) / (((((((x + 43.3400022514d0) * x) + 263.505074721d0) * x) + 313.399215894d0) * x) + 47.066876606d0)
end function
public static double code(double x, double y, double z) {
return ((x - 2.0) * ((((((((x * 4.16438922228) + 78.6994924154) * x) + 137.519416416) * x) + y) * x) + z)) / (((((((x + 43.3400022514) * x) + 263.505074721) * x) + 313.399215894) * x) + 47.066876606);
}
def code(x, y, z): return ((x - 2.0) * ((((((((x * 4.16438922228) + 78.6994924154) * x) + 137.519416416) * x) + y) * x) + z)) / (((((((x + 43.3400022514) * x) + 263.505074721) * x) + 313.399215894) * x) + 47.066876606)
function code(x, y, z) return Float64(Float64(Float64(x - 2.0) * Float64(Float64(Float64(Float64(Float64(Float64(Float64(Float64(x * 4.16438922228) + 78.6994924154) * x) + 137.519416416) * x) + y) * x) + z)) / Float64(Float64(Float64(Float64(Float64(Float64(Float64(x + 43.3400022514) * x) + 263.505074721) * x) + 313.399215894) * x) + 47.066876606)) end
function tmp = code(x, y, z) tmp = ((x - 2.0) * ((((((((x * 4.16438922228) + 78.6994924154) * x) + 137.519416416) * x) + y) * x) + z)) / (((((((x + 43.3400022514) * x) + 263.505074721) * x) + 313.399215894) * x) + 47.066876606); end
code[x_, y_, z_] := N[(N[(N[(x - 2.0), $MachinePrecision] * N[(N[(N[(N[(N[(N[(N[(N[(x * 4.16438922228), $MachinePrecision] + 78.6994924154), $MachinePrecision] * x), $MachinePrecision] + 137.519416416), $MachinePrecision] * x), $MachinePrecision] + y), $MachinePrecision] * x), $MachinePrecision] + z), $MachinePrecision]), $MachinePrecision] / N[(N[(N[(N[(N[(N[(N[(x + 43.3400022514), $MachinePrecision] * x), $MachinePrecision] + 263.505074721), $MachinePrecision] * x), $MachinePrecision] + 313.399215894), $MachinePrecision] * x), $MachinePrecision] + 47.066876606), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{\left(x - 2\right) \cdot \left(\left(\left(\left(x \cdot 4.16438922228 + 78.6994924154\right) \cdot x + 137.519416416\right) \cdot x + y\right) \cdot x + z\right)}{\left(\left(\left(x + 43.3400022514\right) \cdot x + 263.505074721\right) \cdot x + 313.399215894\right) \cdot x + 47.066876606}
\end{array}
(FPCore (x y z)
:precision binary64
(let* ((t_0
(+
(*
x
(+ (* x (+ (* x (+ x 43.3400022514)) 263.505074721)) 313.399215894))
47.066876606))
(t_1
(*
x
(+
(* x (+ (* x (+ (* x 4.16438922228) 78.6994924154)) 137.519416416))
y))))
(if (<= (/ (* (- x 2.0) (+ t_1 z)) t_0) INFINITY)
(* (+ x -2.0) (+ (/ t_1 t_0) (/ z t_0)))
(+
(+
(+ (* x 4.16438922228) (/ 3655.1204654076414 x))
(/ (- y 130977.50649958357) (* x x)))
-110.1139242984811))))
double code(double x, double y, double z) {
double t_0 = (x * ((x * ((x * (x + 43.3400022514)) + 263.505074721)) + 313.399215894)) + 47.066876606;
double t_1 = x * ((x * ((x * ((x * 4.16438922228) + 78.6994924154)) + 137.519416416)) + y);
double tmp;
if ((((x - 2.0) * (t_1 + z)) / t_0) <= ((double) INFINITY)) {
tmp = (x + -2.0) * ((t_1 / t_0) + (z / t_0));
} else {
tmp = (((x * 4.16438922228) + (3655.1204654076414 / x)) + ((y - 130977.50649958357) / (x * x))) + -110.1139242984811;
}
return tmp;
}
public static double code(double x, double y, double z) {
double t_0 = (x * ((x * ((x * (x + 43.3400022514)) + 263.505074721)) + 313.399215894)) + 47.066876606;
double t_1 = x * ((x * ((x * ((x * 4.16438922228) + 78.6994924154)) + 137.519416416)) + y);
double tmp;
if ((((x - 2.0) * (t_1 + z)) / t_0) <= Double.POSITIVE_INFINITY) {
tmp = (x + -2.0) * ((t_1 / t_0) + (z / t_0));
} else {
tmp = (((x * 4.16438922228) + (3655.1204654076414 / x)) + ((y - 130977.50649958357) / (x * x))) + -110.1139242984811;
}
return tmp;
}
def code(x, y, z): t_0 = (x * ((x * ((x * (x + 43.3400022514)) + 263.505074721)) + 313.399215894)) + 47.066876606 t_1 = x * ((x * ((x * ((x * 4.16438922228) + 78.6994924154)) + 137.519416416)) + y) tmp = 0 if (((x - 2.0) * (t_1 + z)) / t_0) <= math.inf: tmp = (x + -2.0) * ((t_1 / t_0) + (z / t_0)) else: tmp = (((x * 4.16438922228) + (3655.1204654076414 / x)) + ((y - 130977.50649958357) / (x * x))) + -110.1139242984811 return tmp
function code(x, y, z) t_0 = Float64(Float64(x * Float64(Float64(x * Float64(Float64(x * Float64(x + 43.3400022514)) + 263.505074721)) + 313.399215894)) + 47.066876606) t_1 = Float64(x * Float64(Float64(x * Float64(Float64(x * Float64(Float64(x * 4.16438922228) + 78.6994924154)) + 137.519416416)) + y)) tmp = 0.0 if (Float64(Float64(Float64(x - 2.0) * Float64(t_1 + z)) / t_0) <= Inf) tmp = Float64(Float64(x + -2.0) * Float64(Float64(t_1 / t_0) + Float64(z / t_0))); else tmp = Float64(Float64(Float64(Float64(x * 4.16438922228) + Float64(3655.1204654076414 / x)) + Float64(Float64(y - 130977.50649958357) / Float64(x * x))) + -110.1139242984811); end return tmp end
function tmp_2 = code(x, y, z) t_0 = (x * ((x * ((x * (x + 43.3400022514)) + 263.505074721)) + 313.399215894)) + 47.066876606; t_1 = x * ((x * ((x * ((x * 4.16438922228) + 78.6994924154)) + 137.519416416)) + y); tmp = 0.0; if ((((x - 2.0) * (t_1 + z)) / t_0) <= Inf) tmp = (x + -2.0) * ((t_1 / t_0) + (z / t_0)); else tmp = (((x * 4.16438922228) + (3655.1204654076414 / x)) + ((y - 130977.50649958357) / (x * x))) + -110.1139242984811; end tmp_2 = tmp; end
code[x_, y_, z_] := Block[{t$95$0 = N[(N[(x * N[(N[(x * N[(N[(x * N[(x + 43.3400022514), $MachinePrecision]), $MachinePrecision] + 263.505074721), $MachinePrecision]), $MachinePrecision] + 313.399215894), $MachinePrecision]), $MachinePrecision] + 47.066876606), $MachinePrecision]}, Block[{t$95$1 = N[(x * N[(N[(x * N[(N[(x * N[(N[(x * 4.16438922228), $MachinePrecision] + 78.6994924154), $MachinePrecision]), $MachinePrecision] + 137.519416416), $MachinePrecision]), $MachinePrecision] + y), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[N[(N[(N[(x - 2.0), $MachinePrecision] * N[(t$95$1 + z), $MachinePrecision]), $MachinePrecision] / t$95$0), $MachinePrecision], Infinity], N[(N[(x + -2.0), $MachinePrecision] * N[(N[(t$95$1 / t$95$0), $MachinePrecision] + N[(z / t$95$0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(N[(N[(x * 4.16438922228), $MachinePrecision] + N[(3655.1204654076414 / x), $MachinePrecision]), $MachinePrecision] + N[(N[(y - 130977.50649958357), $MachinePrecision] / N[(x * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + -110.1139242984811), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := x \cdot \left(x \cdot \left(x \cdot \left(x + 43.3400022514\right) + 263.505074721\right) + 313.399215894\right) + 47.066876606\\
t_1 := x \cdot \left(x \cdot \left(x \cdot \left(x \cdot 4.16438922228 + 78.6994924154\right) + 137.519416416\right) + y\right)\\
\mathbf{if}\;\frac{\left(x - 2\right) \cdot \left(t_1 + z\right)}{t_0} \leq \infty:\\
\;\;\;\;\left(x + -2\right) \cdot \left(\frac{t_1}{t_0} + \frac{z}{t_0}\right)\\
\mathbf{else}:\\
\;\;\;\;\left(\left(x \cdot 4.16438922228 + \frac{3655.1204654076414}{x}\right) + \frac{y - 130977.50649958357}{x \cdot x}\right) + -110.1139242984811\\
\end{array}
\end{array}
if (/.f64 (*.f64 (-.f64 x 2) (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 x 104109730557/25000000000) 393497462077/5000000000) x) 4297481763/31250000) x) y) x) z)) (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 x 216700011257/5000000000) x) 263505074721/1000000000) x) 156699607947/500000000) x) 23533438303/500000000)) < +inf.0Initial program 94.4%
associate-*r/98.8%
sub-neg98.8%
metadata-eval98.8%
*-commutative98.8%
fma-def98.8%
*-commutative98.8%
fma-def98.8%
*-commutative98.8%
fma-def98.8%
fma-def98.8%
*-commutative98.8%
Simplified98.8%
Taylor expanded in z around 0 98.8%
if +inf.0 < (/.f64 (*.f64 (-.f64 x 2) (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 x 104109730557/25000000000) 393497462077/5000000000) x) 4297481763/31250000) x) y) x) z)) (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 x 216700011257/5000000000) x) 263505074721/1000000000) x) 156699607947/500000000) x) 23533438303/500000000)) Initial program 0.0%
associate-*r/0.0%
sub-neg0.0%
metadata-eval0.0%
*-commutative0.0%
fma-def0.0%
*-commutative0.0%
fma-def0.0%
*-commutative0.0%
fma-def0.0%
fma-def0.0%
*-commutative0.0%
Simplified0.0%
Taylor expanded in x around -inf 99.1%
sub-neg99.1%
+-commutative99.1%
mul-1-neg99.1%
unsub-neg99.1%
*-commutative99.1%
fma-def99.1%
associate-*r/99.1%
metadata-eval99.1%
mul-1-neg99.1%
unsub-neg99.1%
unpow299.1%
metadata-eval99.1%
Simplified99.1%
fma-udef99.1%
*-commutative99.1%
Applied egg-rr99.1%
Final simplification98.9%
(FPCore (x y z)
:precision binary64
(let* ((t_0
(/
(*
(- x 2.0)
(+
(*
x
(+
(*
x
(+ (* x (+ (* x 4.16438922228) 78.6994924154)) 137.519416416))
y))
z))
(+
(*
x
(+
(* x (+ (* x (+ x 43.3400022514)) 263.505074721))
313.399215894))
47.066876606))))
(if (<= t_0 2e+293)
t_0
(+
(+
(+ (* x 4.16438922228) (/ 3655.1204654076414 x))
(/ (- y 130977.50649958357) (* x x)))
-110.1139242984811))))
double code(double x, double y, double z) {
double t_0 = ((x - 2.0) * ((x * ((x * ((x * ((x * 4.16438922228) + 78.6994924154)) + 137.519416416)) + y)) + z)) / ((x * ((x * ((x * (x + 43.3400022514)) + 263.505074721)) + 313.399215894)) + 47.066876606);
double tmp;
if (t_0 <= 2e+293) {
tmp = t_0;
} else {
tmp = (((x * 4.16438922228) + (3655.1204654076414 / x)) + ((y - 130977.50649958357) / (x * x))) + -110.1139242984811;
}
return tmp;
}
real(8) function code(x, y, z)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8) :: t_0
real(8) :: tmp
t_0 = ((x - 2.0d0) * ((x * ((x * ((x * ((x * 4.16438922228d0) + 78.6994924154d0)) + 137.519416416d0)) + y)) + z)) / ((x * ((x * ((x * (x + 43.3400022514d0)) + 263.505074721d0)) + 313.399215894d0)) + 47.066876606d0)
if (t_0 <= 2d+293) then
tmp = t_0
else
tmp = (((x * 4.16438922228d0) + (3655.1204654076414d0 / x)) + ((y - 130977.50649958357d0) / (x * x))) + (-110.1139242984811d0)
end if
code = tmp
end function
public static double code(double x, double y, double z) {
double t_0 = ((x - 2.0) * ((x * ((x * ((x * ((x * 4.16438922228) + 78.6994924154)) + 137.519416416)) + y)) + z)) / ((x * ((x * ((x * (x + 43.3400022514)) + 263.505074721)) + 313.399215894)) + 47.066876606);
double tmp;
if (t_0 <= 2e+293) {
tmp = t_0;
} else {
tmp = (((x * 4.16438922228) + (3655.1204654076414 / x)) + ((y - 130977.50649958357) / (x * x))) + -110.1139242984811;
}
return tmp;
}
def code(x, y, z): t_0 = ((x - 2.0) * ((x * ((x * ((x * ((x * 4.16438922228) + 78.6994924154)) + 137.519416416)) + y)) + z)) / ((x * ((x * ((x * (x + 43.3400022514)) + 263.505074721)) + 313.399215894)) + 47.066876606) tmp = 0 if t_0 <= 2e+293: tmp = t_0 else: tmp = (((x * 4.16438922228) + (3655.1204654076414 / x)) + ((y - 130977.50649958357) / (x * x))) + -110.1139242984811 return tmp
function code(x, y, z) t_0 = Float64(Float64(Float64(x - 2.0) * Float64(Float64(x * Float64(Float64(x * Float64(Float64(x * Float64(Float64(x * 4.16438922228) + 78.6994924154)) + 137.519416416)) + y)) + z)) / Float64(Float64(x * Float64(Float64(x * Float64(Float64(x * Float64(x + 43.3400022514)) + 263.505074721)) + 313.399215894)) + 47.066876606)) tmp = 0.0 if (t_0 <= 2e+293) tmp = t_0; else tmp = Float64(Float64(Float64(Float64(x * 4.16438922228) + Float64(3655.1204654076414 / x)) + Float64(Float64(y - 130977.50649958357) / Float64(x * x))) + -110.1139242984811); end return tmp end
function tmp_2 = code(x, y, z) t_0 = ((x - 2.0) * ((x * ((x * ((x * ((x * 4.16438922228) + 78.6994924154)) + 137.519416416)) + y)) + z)) / ((x * ((x * ((x * (x + 43.3400022514)) + 263.505074721)) + 313.399215894)) + 47.066876606); tmp = 0.0; if (t_0 <= 2e+293) tmp = t_0; else tmp = (((x * 4.16438922228) + (3655.1204654076414 / x)) + ((y - 130977.50649958357) / (x * x))) + -110.1139242984811; end tmp_2 = tmp; end
code[x_, y_, z_] := Block[{t$95$0 = N[(N[(N[(x - 2.0), $MachinePrecision] * N[(N[(x * N[(N[(x * N[(N[(x * N[(N[(x * 4.16438922228), $MachinePrecision] + 78.6994924154), $MachinePrecision]), $MachinePrecision] + 137.519416416), $MachinePrecision]), $MachinePrecision] + y), $MachinePrecision]), $MachinePrecision] + z), $MachinePrecision]), $MachinePrecision] / N[(N[(x * N[(N[(x * N[(N[(x * N[(x + 43.3400022514), $MachinePrecision]), $MachinePrecision] + 263.505074721), $MachinePrecision]), $MachinePrecision] + 313.399215894), $MachinePrecision]), $MachinePrecision] + 47.066876606), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[t$95$0, 2e+293], t$95$0, N[(N[(N[(N[(x * 4.16438922228), $MachinePrecision] + N[(3655.1204654076414 / x), $MachinePrecision]), $MachinePrecision] + N[(N[(y - 130977.50649958357), $MachinePrecision] / N[(x * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + -110.1139242984811), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \frac{\left(x - 2\right) \cdot \left(x \cdot \left(x \cdot \left(x \cdot \left(x \cdot 4.16438922228 + 78.6994924154\right) + 137.519416416\right) + y\right) + z\right)}{x \cdot \left(x \cdot \left(x \cdot \left(x + 43.3400022514\right) + 263.505074721\right) + 313.399215894\right) + 47.066876606}\\
\mathbf{if}\;t_0 \leq 2 \cdot 10^{+293}:\\
\;\;\;\;t_0\\
\mathbf{else}:\\
\;\;\;\;\left(\left(x \cdot 4.16438922228 + \frac{3655.1204654076414}{x}\right) + \frac{y - 130977.50649958357}{x \cdot x}\right) + -110.1139242984811\\
\end{array}
\end{array}
if (/.f64 (*.f64 (-.f64 x 2) (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 x 104109730557/25000000000) 393497462077/5000000000) x) 4297481763/31250000) x) y) x) z)) (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 x 216700011257/5000000000) x) 263505074721/1000000000) x) 156699607947/500000000) x) 23533438303/500000000)) < 1.9999999999999998e293Initial program 97.5%
if 1.9999999999999998e293 < (/.f64 (*.f64 (-.f64 x 2) (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 x 104109730557/25000000000) 393497462077/5000000000) x) 4297481763/31250000) x) y) x) z)) (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 (*.f64 (+.f64 x 216700011257/5000000000) x) 263505074721/1000000000) x) 156699607947/500000000) x) 23533438303/500000000)) Initial program 0.2%
associate-*r/4.5%
sub-neg4.5%
metadata-eval4.5%
*-commutative4.5%
fma-def4.5%
*-commutative4.5%
fma-def4.5%
*-commutative4.5%
fma-def4.5%
fma-def4.5%
*-commutative4.5%
Simplified4.5%
Taylor expanded in x around -inf 98.3%
sub-neg98.3%
+-commutative98.3%
mul-1-neg98.3%
unsub-neg98.3%
*-commutative98.3%
fma-def98.3%
associate-*r/98.3%
metadata-eval98.3%
mul-1-neg98.3%
unsub-neg98.3%
unpow298.3%
metadata-eval98.3%
Simplified98.3%
fma-udef98.3%
*-commutative98.3%
Applied egg-rr98.3%
Final simplification97.9%
(FPCore (x y z)
:precision binary64
(if (<= x -82000000000000.0)
(*
(+ x -2.0)
(+
4.16438922228
(/ z (+ 47.066876606 (* x (+ 313.399215894 (* x (* x x))))))))
(if (<= x 800.0)
(/
(* (- x 2.0) (+ z (* x (+ y (* x 137.519416416)))))
(+
(*
x
(+ (* x (+ (* x (+ x 43.3400022514)) 263.505074721)) 313.399215894))
47.066876606))
(+
(+
(+ (* x 4.16438922228) (/ 3655.1204654076414 x))
(/ (- y 130977.50649958357) (* x x)))
-110.1139242984811))))
double code(double x, double y, double z) {
double tmp;
if (x <= -82000000000000.0) {
tmp = (x + -2.0) * (4.16438922228 + (z / (47.066876606 + (x * (313.399215894 + (x * (x * x)))))));
} else if (x <= 800.0) {
tmp = ((x - 2.0) * (z + (x * (y + (x * 137.519416416))))) / ((x * ((x * ((x * (x + 43.3400022514)) + 263.505074721)) + 313.399215894)) + 47.066876606);
} else {
tmp = (((x * 4.16438922228) + (3655.1204654076414 / x)) + ((y - 130977.50649958357) / (x * x))) + -110.1139242984811;
}
return tmp;
}
real(8) function code(x, y, z)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8) :: tmp
if (x <= (-82000000000000.0d0)) then
tmp = (x + (-2.0d0)) * (4.16438922228d0 + (z / (47.066876606d0 + (x * (313.399215894d0 + (x * (x * x)))))))
else if (x <= 800.0d0) then
tmp = ((x - 2.0d0) * (z + (x * (y + (x * 137.519416416d0))))) / ((x * ((x * ((x * (x + 43.3400022514d0)) + 263.505074721d0)) + 313.399215894d0)) + 47.066876606d0)
else
tmp = (((x * 4.16438922228d0) + (3655.1204654076414d0 / x)) + ((y - 130977.50649958357d0) / (x * x))) + (-110.1139242984811d0)
end if
code = tmp
end function
public static double code(double x, double y, double z) {
double tmp;
if (x <= -82000000000000.0) {
tmp = (x + -2.0) * (4.16438922228 + (z / (47.066876606 + (x * (313.399215894 + (x * (x * x)))))));
} else if (x <= 800.0) {
tmp = ((x - 2.0) * (z + (x * (y + (x * 137.519416416))))) / ((x * ((x * ((x * (x + 43.3400022514)) + 263.505074721)) + 313.399215894)) + 47.066876606);
} else {
tmp = (((x * 4.16438922228) + (3655.1204654076414 / x)) + ((y - 130977.50649958357) / (x * x))) + -110.1139242984811;
}
return tmp;
}
def code(x, y, z): tmp = 0 if x <= -82000000000000.0: tmp = (x + -2.0) * (4.16438922228 + (z / (47.066876606 + (x * (313.399215894 + (x * (x * x))))))) elif x <= 800.0: tmp = ((x - 2.0) * (z + (x * (y + (x * 137.519416416))))) / ((x * ((x * ((x * (x + 43.3400022514)) + 263.505074721)) + 313.399215894)) + 47.066876606) else: tmp = (((x * 4.16438922228) + (3655.1204654076414 / x)) + ((y - 130977.50649958357) / (x * x))) + -110.1139242984811 return tmp
function code(x, y, z) tmp = 0.0 if (x <= -82000000000000.0) tmp = Float64(Float64(x + -2.0) * Float64(4.16438922228 + Float64(z / Float64(47.066876606 + Float64(x * Float64(313.399215894 + Float64(x * Float64(x * x)))))))); elseif (x <= 800.0) tmp = Float64(Float64(Float64(x - 2.0) * Float64(z + Float64(x * Float64(y + Float64(x * 137.519416416))))) / Float64(Float64(x * Float64(Float64(x * Float64(Float64(x * Float64(x + 43.3400022514)) + 263.505074721)) + 313.399215894)) + 47.066876606)); else tmp = Float64(Float64(Float64(Float64(x * 4.16438922228) + Float64(3655.1204654076414 / x)) + Float64(Float64(y - 130977.50649958357) / Float64(x * x))) + -110.1139242984811); end return tmp end
function tmp_2 = code(x, y, z) tmp = 0.0; if (x <= -82000000000000.0) tmp = (x + -2.0) * (4.16438922228 + (z / (47.066876606 + (x * (313.399215894 + (x * (x * x))))))); elseif (x <= 800.0) tmp = ((x - 2.0) * (z + (x * (y + (x * 137.519416416))))) / ((x * ((x * ((x * (x + 43.3400022514)) + 263.505074721)) + 313.399215894)) + 47.066876606); else tmp = (((x * 4.16438922228) + (3655.1204654076414 / x)) + ((y - 130977.50649958357) / (x * x))) + -110.1139242984811; end tmp_2 = tmp; end
code[x_, y_, z_] := If[LessEqual[x, -82000000000000.0], N[(N[(x + -2.0), $MachinePrecision] * N[(4.16438922228 + N[(z / N[(47.066876606 + N[(x * N[(313.399215894 + N[(x * N[(x * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[x, 800.0], N[(N[(N[(x - 2.0), $MachinePrecision] * N[(z + N[(x * N[(y + N[(x * 137.519416416), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(N[(x * N[(N[(x * N[(N[(x * N[(x + 43.3400022514), $MachinePrecision]), $MachinePrecision] + 263.505074721), $MachinePrecision]), $MachinePrecision] + 313.399215894), $MachinePrecision]), $MachinePrecision] + 47.066876606), $MachinePrecision]), $MachinePrecision], N[(N[(N[(N[(x * 4.16438922228), $MachinePrecision] + N[(3655.1204654076414 / x), $MachinePrecision]), $MachinePrecision] + N[(N[(y - 130977.50649958357), $MachinePrecision] / N[(x * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + -110.1139242984811), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq -82000000000000:\\
\;\;\;\;\left(x + -2\right) \cdot \left(4.16438922228 + \frac{z}{47.066876606 + x \cdot \left(313.399215894 + x \cdot \left(x \cdot x\right)\right)}\right)\\
\mathbf{elif}\;x \leq 800:\\
\;\;\;\;\frac{\left(x - 2\right) \cdot \left(z + x \cdot \left(y + x \cdot 137.519416416\right)\right)}{x \cdot \left(x \cdot \left(x \cdot \left(x + 43.3400022514\right) + 263.505074721\right) + 313.399215894\right) + 47.066876606}\\
\mathbf{else}:\\
\;\;\;\;\left(\left(x \cdot 4.16438922228 + \frac{3655.1204654076414}{x}\right) + \frac{y - 130977.50649958357}{x \cdot x}\right) + -110.1139242984811\\
\end{array}
\end{array}
if x < -8.2e13Initial program 12.8%
associate-*r/17.3%
sub-neg17.3%
metadata-eval17.3%
*-commutative17.3%
fma-def17.3%
*-commutative17.3%
fma-def17.3%
*-commutative17.3%
fma-def17.3%
fma-def17.3%
*-commutative17.3%
Simplified17.3%
Taylor expanded in z around 0 17.3%
Taylor expanded in x around inf 97.3%
Taylor expanded in x around inf 42.3%
cube-mult42.3%
unpow242.3%
distribute-rgt-out97.3%
unpow297.3%
associate-*r*97.3%
+-commutative97.3%
Simplified97.3%
Taylor expanded in x around inf 97.3%
unpow297.3%
Simplified97.3%
if -8.2e13 < x < 800Initial program 99.6%
Taylor expanded in x around 0 98.3%
*-commutative98.3%
Simplified98.3%
if 800 < x Initial program 17.0%
associate-*r/22.3%
sub-neg22.3%
metadata-eval22.3%
*-commutative22.3%
fma-def22.3%
*-commutative22.3%
fma-def22.3%
*-commutative22.3%
fma-def22.3%
fma-def22.3%
*-commutative22.3%
Simplified22.3%
Taylor expanded in x around -inf 95.8%
sub-neg95.8%
+-commutative95.8%
mul-1-neg95.8%
unsub-neg95.8%
*-commutative95.8%
fma-def95.8%
associate-*r/95.8%
metadata-eval95.8%
mul-1-neg95.8%
unsub-neg95.8%
unpow295.8%
metadata-eval95.8%
Simplified95.8%
fma-udef95.8%
*-commutative95.8%
Applied egg-rr95.8%
Final simplification97.4%
(FPCore (x y z)
:precision binary64
(if (<= x -0.075)
(*
(+ x -2.0)
(+
4.16438922228
(/
z
(+
(*
x
(+ (* x (+ (* x (+ x 43.3400022514)) 263.505074721)) 313.399215894))
47.066876606))))
(if (<= x 95.0)
(/
(* (- x 2.0) (+ z (* x (+ y (* x 137.519416416)))))
(+ 47.066876606 (* x (+ 313.399215894 (* x 263.505074721)))))
(+
(+
(+ (* x 4.16438922228) (/ 3655.1204654076414 x))
(/ (- y 130977.50649958357) (* x x)))
-110.1139242984811))))
double code(double x, double y, double z) {
double tmp;
if (x <= -0.075) {
tmp = (x + -2.0) * (4.16438922228 + (z / ((x * ((x * ((x * (x + 43.3400022514)) + 263.505074721)) + 313.399215894)) + 47.066876606)));
} else if (x <= 95.0) {
tmp = ((x - 2.0) * (z + (x * (y + (x * 137.519416416))))) / (47.066876606 + (x * (313.399215894 + (x * 263.505074721))));
} else {
tmp = (((x * 4.16438922228) + (3655.1204654076414 / x)) + ((y - 130977.50649958357) / (x * x))) + -110.1139242984811;
}
return tmp;
}
real(8) function code(x, y, z)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8) :: tmp
if (x <= (-0.075d0)) then
tmp = (x + (-2.0d0)) * (4.16438922228d0 + (z / ((x * ((x * ((x * (x + 43.3400022514d0)) + 263.505074721d0)) + 313.399215894d0)) + 47.066876606d0)))
else if (x <= 95.0d0) then
tmp = ((x - 2.0d0) * (z + (x * (y + (x * 137.519416416d0))))) / (47.066876606d0 + (x * (313.399215894d0 + (x * 263.505074721d0))))
else
tmp = (((x * 4.16438922228d0) + (3655.1204654076414d0 / x)) + ((y - 130977.50649958357d0) / (x * x))) + (-110.1139242984811d0)
end if
code = tmp
end function
public static double code(double x, double y, double z) {
double tmp;
if (x <= -0.075) {
tmp = (x + -2.0) * (4.16438922228 + (z / ((x * ((x * ((x * (x + 43.3400022514)) + 263.505074721)) + 313.399215894)) + 47.066876606)));
} else if (x <= 95.0) {
tmp = ((x - 2.0) * (z + (x * (y + (x * 137.519416416))))) / (47.066876606 + (x * (313.399215894 + (x * 263.505074721))));
} else {
tmp = (((x * 4.16438922228) + (3655.1204654076414 / x)) + ((y - 130977.50649958357) / (x * x))) + -110.1139242984811;
}
return tmp;
}
def code(x, y, z): tmp = 0 if x <= -0.075: tmp = (x + -2.0) * (4.16438922228 + (z / ((x * ((x * ((x * (x + 43.3400022514)) + 263.505074721)) + 313.399215894)) + 47.066876606))) elif x <= 95.0: tmp = ((x - 2.0) * (z + (x * (y + (x * 137.519416416))))) / (47.066876606 + (x * (313.399215894 + (x * 263.505074721)))) else: tmp = (((x * 4.16438922228) + (3655.1204654076414 / x)) + ((y - 130977.50649958357) / (x * x))) + -110.1139242984811 return tmp
function code(x, y, z) tmp = 0.0 if (x <= -0.075) tmp = Float64(Float64(x + -2.0) * Float64(4.16438922228 + Float64(z / Float64(Float64(x * Float64(Float64(x * Float64(Float64(x * Float64(x + 43.3400022514)) + 263.505074721)) + 313.399215894)) + 47.066876606)))); elseif (x <= 95.0) tmp = Float64(Float64(Float64(x - 2.0) * Float64(z + Float64(x * Float64(y + Float64(x * 137.519416416))))) / Float64(47.066876606 + Float64(x * Float64(313.399215894 + Float64(x * 263.505074721))))); else tmp = Float64(Float64(Float64(Float64(x * 4.16438922228) + Float64(3655.1204654076414 / x)) + Float64(Float64(y - 130977.50649958357) / Float64(x * x))) + -110.1139242984811); end return tmp end
function tmp_2 = code(x, y, z) tmp = 0.0; if (x <= -0.075) tmp = (x + -2.0) * (4.16438922228 + (z / ((x * ((x * ((x * (x + 43.3400022514)) + 263.505074721)) + 313.399215894)) + 47.066876606))); elseif (x <= 95.0) tmp = ((x - 2.0) * (z + (x * (y + (x * 137.519416416))))) / (47.066876606 + (x * (313.399215894 + (x * 263.505074721)))); else tmp = (((x * 4.16438922228) + (3655.1204654076414 / x)) + ((y - 130977.50649958357) / (x * x))) + -110.1139242984811; end tmp_2 = tmp; end
code[x_, y_, z_] := If[LessEqual[x, -0.075], N[(N[(x + -2.0), $MachinePrecision] * N[(4.16438922228 + N[(z / N[(N[(x * N[(N[(x * N[(N[(x * N[(x + 43.3400022514), $MachinePrecision]), $MachinePrecision] + 263.505074721), $MachinePrecision]), $MachinePrecision] + 313.399215894), $MachinePrecision]), $MachinePrecision] + 47.066876606), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[x, 95.0], N[(N[(N[(x - 2.0), $MachinePrecision] * N[(z + N[(x * N[(y + N[(x * 137.519416416), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(47.066876606 + N[(x * N[(313.399215894 + N[(x * 263.505074721), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(N[(N[(x * 4.16438922228), $MachinePrecision] + N[(3655.1204654076414 / x), $MachinePrecision]), $MachinePrecision] + N[(N[(y - 130977.50649958357), $MachinePrecision] / N[(x * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + -110.1139242984811), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq -0.075:\\
\;\;\;\;\left(x + -2\right) \cdot \left(4.16438922228 + \frac{z}{x \cdot \left(x \cdot \left(x \cdot \left(x + 43.3400022514\right) + 263.505074721\right) + 313.399215894\right) + 47.066876606}\right)\\
\mathbf{elif}\;x \leq 95:\\
\;\;\;\;\frac{\left(x - 2\right) \cdot \left(z + x \cdot \left(y + x \cdot 137.519416416\right)\right)}{47.066876606 + x \cdot \left(313.399215894 + x \cdot 263.505074721\right)}\\
\mathbf{else}:\\
\;\;\;\;\left(\left(x \cdot 4.16438922228 + \frac{3655.1204654076414}{x}\right) + \frac{y - 130977.50649958357}{x \cdot x}\right) + -110.1139242984811\\
\end{array}
\end{array}
if x < -0.0749999999999999972Initial program 16.7%
associate-*r/21.0%
sub-neg21.0%
metadata-eval21.0%
*-commutative21.0%
fma-def21.0%
*-commutative21.0%
fma-def21.0%
*-commutative21.0%
fma-def21.0%
fma-def21.0%
*-commutative21.0%
Simplified20.9%
Taylor expanded in z around 0 21.0%
Taylor expanded in x around inf 95.2%
if -0.0749999999999999972 < x < 95Initial program 99.6%
Taylor expanded in x around 0 99.6%
*-commutative99.6%
Simplified99.6%
Taylor expanded in x around 0 97.4%
*-commutative97.4%
Simplified97.4%
if 95 < x Initial program 18.1%
associate-*r/23.4%
sub-neg23.4%
metadata-eval23.4%
*-commutative23.4%
fma-def23.4%
*-commutative23.4%
fma-def23.4%
*-commutative23.4%
fma-def23.4%
fma-def23.4%
*-commutative23.4%
Simplified23.4%
Taylor expanded in x around -inf 94.8%
sub-neg94.8%
+-commutative94.8%
mul-1-neg94.8%
unsub-neg94.8%
*-commutative94.8%
fma-def94.8%
associate-*r/94.8%
metadata-eval94.8%
mul-1-neg94.8%
unsub-neg94.8%
unpow294.8%
metadata-eval94.8%
Simplified94.8%
fma-udef94.8%
*-commutative94.8%
Applied egg-rr94.8%
Final simplification96.1%
(FPCore (x y z)
:precision binary64
(if (<= x -0.0012)
(*
(+ x -2.0)
(+
4.16438922228
(/
z
(+
(*
x
(+ (* x (+ (* x (+ x 43.3400022514)) 263.505074721)) 313.399215894))
47.066876606))))
(if (<= x 50.0)
(/
(* (- x 2.0) (+ z (* x (+ y (* x 137.519416416)))))
(+ 47.066876606 (* x 313.399215894)))
(+
(+
(+ (* x 4.16438922228) (/ 3655.1204654076414 x))
(/ (- y 130977.50649958357) (* x x)))
-110.1139242984811))))
double code(double x, double y, double z) {
double tmp;
if (x <= -0.0012) {
tmp = (x + -2.0) * (4.16438922228 + (z / ((x * ((x * ((x * (x + 43.3400022514)) + 263.505074721)) + 313.399215894)) + 47.066876606)));
} else if (x <= 50.0) {
tmp = ((x - 2.0) * (z + (x * (y + (x * 137.519416416))))) / (47.066876606 + (x * 313.399215894));
} else {
tmp = (((x * 4.16438922228) + (3655.1204654076414 / x)) + ((y - 130977.50649958357) / (x * x))) + -110.1139242984811;
}
return tmp;
}
real(8) function code(x, y, z)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8) :: tmp
if (x <= (-0.0012d0)) then
tmp = (x + (-2.0d0)) * (4.16438922228d0 + (z / ((x * ((x * ((x * (x + 43.3400022514d0)) + 263.505074721d0)) + 313.399215894d0)) + 47.066876606d0)))
else if (x <= 50.0d0) then
tmp = ((x - 2.0d0) * (z + (x * (y + (x * 137.519416416d0))))) / (47.066876606d0 + (x * 313.399215894d0))
else
tmp = (((x * 4.16438922228d0) + (3655.1204654076414d0 / x)) + ((y - 130977.50649958357d0) / (x * x))) + (-110.1139242984811d0)
end if
code = tmp
end function
public static double code(double x, double y, double z) {
double tmp;
if (x <= -0.0012) {
tmp = (x + -2.0) * (4.16438922228 + (z / ((x * ((x * ((x * (x + 43.3400022514)) + 263.505074721)) + 313.399215894)) + 47.066876606)));
} else if (x <= 50.0) {
tmp = ((x - 2.0) * (z + (x * (y + (x * 137.519416416))))) / (47.066876606 + (x * 313.399215894));
} else {
tmp = (((x * 4.16438922228) + (3655.1204654076414 / x)) + ((y - 130977.50649958357) / (x * x))) + -110.1139242984811;
}
return tmp;
}
def code(x, y, z): tmp = 0 if x <= -0.0012: tmp = (x + -2.0) * (4.16438922228 + (z / ((x * ((x * ((x * (x + 43.3400022514)) + 263.505074721)) + 313.399215894)) + 47.066876606))) elif x <= 50.0: tmp = ((x - 2.0) * (z + (x * (y + (x * 137.519416416))))) / (47.066876606 + (x * 313.399215894)) else: tmp = (((x * 4.16438922228) + (3655.1204654076414 / x)) + ((y - 130977.50649958357) / (x * x))) + -110.1139242984811 return tmp
function code(x, y, z) tmp = 0.0 if (x <= -0.0012) tmp = Float64(Float64(x + -2.0) * Float64(4.16438922228 + Float64(z / Float64(Float64(x * Float64(Float64(x * Float64(Float64(x * Float64(x + 43.3400022514)) + 263.505074721)) + 313.399215894)) + 47.066876606)))); elseif (x <= 50.0) tmp = Float64(Float64(Float64(x - 2.0) * Float64(z + Float64(x * Float64(y + Float64(x * 137.519416416))))) / Float64(47.066876606 + Float64(x * 313.399215894))); else tmp = Float64(Float64(Float64(Float64(x * 4.16438922228) + Float64(3655.1204654076414 / x)) + Float64(Float64(y - 130977.50649958357) / Float64(x * x))) + -110.1139242984811); end return tmp end
function tmp_2 = code(x, y, z) tmp = 0.0; if (x <= -0.0012) tmp = (x + -2.0) * (4.16438922228 + (z / ((x * ((x * ((x * (x + 43.3400022514)) + 263.505074721)) + 313.399215894)) + 47.066876606))); elseif (x <= 50.0) tmp = ((x - 2.0) * (z + (x * (y + (x * 137.519416416))))) / (47.066876606 + (x * 313.399215894)); else tmp = (((x * 4.16438922228) + (3655.1204654076414 / x)) + ((y - 130977.50649958357) / (x * x))) + -110.1139242984811; end tmp_2 = tmp; end
code[x_, y_, z_] := If[LessEqual[x, -0.0012], N[(N[(x + -2.0), $MachinePrecision] * N[(4.16438922228 + N[(z / N[(N[(x * N[(N[(x * N[(N[(x * N[(x + 43.3400022514), $MachinePrecision]), $MachinePrecision] + 263.505074721), $MachinePrecision]), $MachinePrecision] + 313.399215894), $MachinePrecision]), $MachinePrecision] + 47.066876606), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[x, 50.0], N[(N[(N[(x - 2.0), $MachinePrecision] * N[(z + N[(x * N[(y + N[(x * 137.519416416), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(47.066876606 + N[(x * 313.399215894), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(N[(N[(x * 4.16438922228), $MachinePrecision] + N[(3655.1204654076414 / x), $MachinePrecision]), $MachinePrecision] + N[(N[(y - 130977.50649958357), $MachinePrecision] / N[(x * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + -110.1139242984811), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq -0.0012:\\
\;\;\;\;\left(x + -2\right) \cdot \left(4.16438922228 + \frac{z}{x \cdot \left(x \cdot \left(x \cdot \left(x + 43.3400022514\right) + 263.505074721\right) + 313.399215894\right) + 47.066876606}\right)\\
\mathbf{elif}\;x \leq 50:\\
\;\;\;\;\frac{\left(x - 2\right) \cdot \left(z + x \cdot \left(y + x \cdot 137.519416416\right)\right)}{47.066876606 + x \cdot 313.399215894}\\
\mathbf{else}:\\
\;\;\;\;\left(\left(x \cdot 4.16438922228 + \frac{3655.1204654076414}{x}\right) + \frac{y - 130977.50649958357}{x \cdot x}\right) + -110.1139242984811\\
\end{array}
\end{array}
if x < -0.00119999999999999989Initial program 16.7%
associate-*r/21.0%
sub-neg21.0%
metadata-eval21.0%
*-commutative21.0%
fma-def21.0%
*-commutative21.0%
fma-def21.0%
*-commutative21.0%
fma-def21.0%
fma-def21.0%
*-commutative21.0%
Simplified20.9%
Taylor expanded in z around 0 21.0%
Taylor expanded in x around inf 95.2%
if -0.00119999999999999989 < x < 50Initial program 99.6%
Taylor expanded in x around 0 99.6%
*-commutative99.6%
Simplified99.6%
Taylor expanded in x around 0 96.7%
*-commutative96.7%
Simplified96.7%
if 50 < x Initial program 18.1%
associate-*r/23.4%
sub-neg23.4%
metadata-eval23.4%
*-commutative23.4%
fma-def23.4%
*-commutative23.4%
fma-def23.4%
*-commutative23.4%
fma-def23.4%
fma-def23.4%
*-commutative23.4%
Simplified23.4%
Taylor expanded in x around -inf 94.8%
sub-neg94.8%
+-commutative94.8%
mul-1-neg94.8%
unsub-neg94.8%
*-commutative94.8%
fma-def94.8%
associate-*r/94.8%
metadata-eval94.8%
mul-1-neg94.8%
unsub-neg94.8%
unpow294.8%
metadata-eval94.8%
Simplified94.8%
fma-udef94.8%
*-commutative94.8%
Applied egg-rr94.8%
Final simplification95.8%
(FPCore (x y z)
:precision binary64
(if (<= x -35.0)
(*
(+ x -2.0)
(+
4.16438922228
(/ z (+ 47.066876606 (* x (+ 313.399215894 (* x (* x x))))))))
(if (<= x 45.0)
(/
(* (- x 2.0) (+ z (* x (+ y (* x 137.519416416)))))
(+ 47.066876606 (* x 313.399215894)))
(+
(+
(+ (* x 4.16438922228) (/ 3655.1204654076414 x))
(/ (- y 130977.50649958357) (* x x)))
-110.1139242984811))))
double code(double x, double y, double z) {
double tmp;
if (x <= -35.0) {
tmp = (x + -2.0) * (4.16438922228 + (z / (47.066876606 + (x * (313.399215894 + (x * (x * x)))))));
} else if (x <= 45.0) {
tmp = ((x - 2.0) * (z + (x * (y + (x * 137.519416416))))) / (47.066876606 + (x * 313.399215894));
} else {
tmp = (((x * 4.16438922228) + (3655.1204654076414 / x)) + ((y - 130977.50649958357) / (x * x))) + -110.1139242984811;
}
return tmp;
}
real(8) function code(x, y, z)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8) :: tmp
if (x <= (-35.0d0)) then
tmp = (x + (-2.0d0)) * (4.16438922228d0 + (z / (47.066876606d0 + (x * (313.399215894d0 + (x * (x * x)))))))
else if (x <= 45.0d0) then
tmp = ((x - 2.0d0) * (z + (x * (y + (x * 137.519416416d0))))) / (47.066876606d0 + (x * 313.399215894d0))
else
tmp = (((x * 4.16438922228d0) + (3655.1204654076414d0 / x)) + ((y - 130977.50649958357d0) / (x * x))) + (-110.1139242984811d0)
end if
code = tmp
end function
public static double code(double x, double y, double z) {
double tmp;
if (x <= -35.0) {
tmp = (x + -2.0) * (4.16438922228 + (z / (47.066876606 + (x * (313.399215894 + (x * (x * x)))))));
} else if (x <= 45.0) {
tmp = ((x - 2.0) * (z + (x * (y + (x * 137.519416416))))) / (47.066876606 + (x * 313.399215894));
} else {
tmp = (((x * 4.16438922228) + (3655.1204654076414 / x)) + ((y - 130977.50649958357) / (x * x))) + -110.1139242984811;
}
return tmp;
}
def code(x, y, z): tmp = 0 if x <= -35.0: tmp = (x + -2.0) * (4.16438922228 + (z / (47.066876606 + (x * (313.399215894 + (x * (x * x))))))) elif x <= 45.0: tmp = ((x - 2.0) * (z + (x * (y + (x * 137.519416416))))) / (47.066876606 + (x * 313.399215894)) else: tmp = (((x * 4.16438922228) + (3655.1204654076414 / x)) + ((y - 130977.50649958357) / (x * x))) + -110.1139242984811 return tmp
function code(x, y, z) tmp = 0.0 if (x <= -35.0) tmp = Float64(Float64(x + -2.0) * Float64(4.16438922228 + Float64(z / Float64(47.066876606 + Float64(x * Float64(313.399215894 + Float64(x * Float64(x * x)))))))); elseif (x <= 45.0) tmp = Float64(Float64(Float64(x - 2.0) * Float64(z + Float64(x * Float64(y + Float64(x * 137.519416416))))) / Float64(47.066876606 + Float64(x * 313.399215894))); else tmp = Float64(Float64(Float64(Float64(x * 4.16438922228) + Float64(3655.1204654076414 / x)) + Float64(Float64(y - 130977.50649958357) / Float64(x * x))) + -110.1139242984811); end return tmp end
function tmp_2 = code(x, y, z) tmp = 0.0; if (x <= -35.0) tmp = (x + -2.0) * (4.16438922228 + (z / (47.066876606 + (x * (313.399215894 + (x * (x * x))))))); elseif (x <= 45.0) tmp = ((x - 2.0) * (z + (x * (y + (x * 137.519416416))))) / (47.066876606 + (x * 313.399215894)); else tmp = (((x * 4.16438922228) + (3655.1204654076414 / x)) + ((y - 130977.50649958357) / (x * x))) + -110.1139242984811; end tmp_2 = tmp; end
code[x_, y_, z_] := If[LessEqual[x, -35.0], N[(N[(x + -2.0), $MachinePrecision] * N[(4.16438922228 + N[(z / N[(47.066876606 + N[(x * N[(313.399215894 + N[(x * N[(x * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[x, 45.0], N[(N[(N[(x - 2.0), $MachinePrecision] * N[(z + N[(x * N[(y + N[(x * 137.519416416), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(47.066876606 + N[(x * 313.399215894), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(N[(N[(x * 4.16438922228), $MachinePrecision] + N[(3655.1204654076414 / x), $MachinePrecision]), $MachinePrecision] + N[(N[(y - 130977.50649958357), $MachinePrecision] / N[(x * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + -110.1139242984811), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq -35:\\
\;\;\;\;\left(x + -2\right) \cdot \left(4.16438922228 + \frac{z}{47.066876606 + x \cdot \left(313.399215894 + x \cdot \left(x \cdot x\right)\right)}\right)\\
\mathbf{elif}\;x \leq 45:\\
\;\;\;\;\frac{\left(x - 2\right) \cdot \left(z + x \cdot \left(y + x \cdot 137.519416416\right)\right)}{47.066876606 + x \cdot 313.399215894}\\
\mathbf{else}:\\
\;\;\;\;\left(\left(x \cdot 4.16438922228 + \frac{3655.1204654076414}{x}\right) + \frac{y - 130977.50649958357}{x \cdot x}\right) + -110.1139242984811\\
\end{array}
\end{array}
if x < -35Initial program 16.7%
associate-*r/21.0%
sub-neg21.0%
metadata-eval21.0%
*-commutative21.0%
fma-def21.0%
*-commutative21.0%
fma-def21.0%
*-commutative21.0%
fma-def21.0%
fma-def21.0%
*-commutative21.0%
Simplified20.9%
Taylor expanded in z around 0 21.0%
Taylor expanded in x around inf 95.2%
Taylor expanded in x around inf 41.8%
cube-mult41.8%
unpow241.8%
distribute-rgt-out94.3%
unpow294.3%
associate-*r*94.3%
+-commutative94.3%
Simplified94.3%
Taylor expanded in x around inf 94.1%
unpow294.1%
Simplified94.1%
if -35 < x < 45Initial program 99.6%
Taylor expanded in x around 0 99.6%
*-commutative99.6%
Simplified99.6%
Taylor expanded in x around 0 96.7%
*-commutative96.7%
Simplified96.7%
if 45 < x Initial program 18.1%
associate-*r/23.4%
sub-neg23.4%
metadata-eval23.4%
*-commutative23.4%
fma-def23.4%
*-commutative23.4%
fma-def23.4%
*-commutative23.4%
fma-def23.4%
fma-def23.4%
*-commutative23.4%
Simplified23.4%
Taylor expanded in x around -inf 94.8%
sub-neg94.8%
+-commutative94.8%
mul-1-neg94.8%
unsub-neg94.8%
*-commutative94.8%
fma-def94.8%
associate-*r/94.8%
metadata-eval94.8%
mul-1-neg94.8%
unsub-neg94.8%
unpow294.8%
metadata-eval94.8%
Simplified94.8%
fma-udef94.8%
*-commutative94.8%
Applied egg-rr94.8%
Final simplification95.5%
(FPCore (x y z)
:precision binary64
(if (<= x -35.0)
(*
(+ x -2.0)
(+
4.16438922228
(/
z
(+
47.066876606
(* x (+ 313.399215894 (* x (* x (+ x 43.3400022514)))))))))
(if (<= x 56.0)
(/
(* (- x 2.0) (+ z (* x (+ y (* x 137.519416416)))))
(+ 47.066876606 (* x 313.399215894)))
(+
(+
(+ (* x 4.16438922228) (/ 3655.1204654076414 x))
(/ (- y 130977.50649958357) (* x x)))
-110.1139242984811))))
double code(double x, double y, double z) {
double tmp;
if (x <= -35.0) {
tmp = (x + -2.0) * (4.16438922228 + (z / (47.066876606 + (x * (313.399215894 + (x * (x * (x + 43.3400022514))))))));
} else if (x <= 56.0) {
tmp = ((x - 2.0) * (z + (x * (y + (x * 137.519416416))))) / (47.066876606 + (x * 313.399215894));
} else {
tmp = (((x * 4.16438922228) + (3655.1204654076414 / x)) + ((y - 130977.50649958357) / (x * x))) + -110.1139242984811;
}
return tmp;
}
real(8) function code(x, y, z)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8) :: tmp
if (x <= (-35.0d0)) then
tmp = (x + (-2.0d0)) * (4.16438922228d0 + (z / (47.066876606d0 + (x * (313.399215894d0 + (x * (x * (x + 43.3400022514d0))))))))
else if (x <= 56.0d0) then
tmp = ((x - 2.0d0) * (z + (x * (y + (x * 137.519416416d0))))) / (47.066876606d0 + (x * 313.399215894d0))
else
tmp = (((x * 4.16438922228d0) + (3655.1204654076414d0 / x)) + ((y - 130977.50649958357d0) / (x * x))) + (-110.1139242984811d0)
end if
code = tmp
end function
public static double code(double x, double y, double z) {
double tmp;
if (x <= -35.0) {
tmp = (x + -2.0) * (4.16438922228 + (z / (47.066876606 + (x * (313.399215894 + (x * (x * (x + 43.3400022514))))))));
} else if (x <= 56.0) {
tmp = ((x - 2.0) * (z + (x * (y + (x * 137.519416416))))) / (47.066876606 + (x * 313.399215894));
} else {
tmp = (((x * 4.16438922228) + (3655.1204654076414 / x)) + ((y - 130977.50649958357) / (x * x))) + -110.1139242984811;
}
return tmp;
}
def code(x, y, z): tmp = 0 if x <= -35.0: tmp = (x + -2.0) * (4.16438922228 + (z / (47.066876606 + (x * (313.399215894 + (x * (x * (x + 43.3400022514)))))))) elif x <= 56.0: tmp = ((x - 2.0) * (z + (x * (y + (x * 137.519416416))))) / (47.066876606 + (x * 313.399215894)) else: tmp = (((x * 4.16438922228) + (3655.1204654076414 / x)) + ((y - 130977.50649958357) / (x * x))) + -110.1139242984811 return tmp
function code(x, y, z) tmp = 0.0 if (x <= -35.0) tmp = Float64(Float64(x + -2.0) * Float64(4.16438922228 + Float64(z / Float64(47.066876606 + Float64(x * Float64(313.399215894 + Float64(x * Float64(x * Float64(x + 43.3400022514))))))))); elseif (x <= 56.0) tmp = Float64(Float64(Float64(x - 2.0) * Float64(z + Float64(x * Float64(y + Float64(x * 137.519416416))))) / Float64(47.066876606 + Float64(x * 313.399215894))); else tmp = Float64(Float64(Float64(Float64(x * 4.16438922228) + Float64(3655.1204654076414 / x)) + Float64(Float64(y - 130977.50649958357) / Float64(x * x))) + -110.1139242984811); end return tmp end
function tmp_2 = code(x, y, z) tmp = 0.0; if (x <= -35.0) tmp = (x + -2.0) * (4.16438922228 + (z / (47.066876606 + (x * (313.399215894 + (x * (x * (x + 43.3400022514)))))))); elseif (x <= 56.0) tmp = ((x - 2.0) * (z + (x * (y + (x * 137.519416416))))) / (47.066876606 + (x * 313.399215894)); else tmp = (((x * 4.16438922228) + (3655.1204654076414 / x)) + ((y - 130977.50649958357) / (x * x))) + -110.1139242984811; end tmp_2 = tmp; end
code[x_, y_, z_] := If[LessEqual[x, -35.0], N[(N[(x + -2.0), $MachinePrecision] * N[(4.16438922228 + N[(z / N[(47.066876606 + N[(x * N[(313.399215894 + N[(x * N[(x * N[(x + 43.3400022514), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[x, 56.0], N[(N[(N[(x - 2.0), $MachinePrecision] * N[(z + N[(x * N[(y + N[(x * 137.519416416), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(47.066876606 + N[(x * 313.399215894), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(N[(N[(x * 4.16438922228), $MachinePrecision] + N[(3655.1204654076414 / x), $MachinePrecision]), $MachinePrecision] + N[(N[(y - 130977.50649958357), $MachinePrecision] / N[(x * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + -110.1139242984811), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq -35:\\
\;\;\;\;\left(x + -2\right) \cdot \left(4.16438922228 + \frac{z}{47.066876606 + x \cdot \left(313.399215894 + x \cdot \left(x \cdot \left(x + 43.3400022514\right)\right)\right)}\right)\\
\mathbf{elif}\;x \leq 56:\\
\;\;\;\;\frac{\left(x - 2\right) \cdot \left(z + x \cdot \left(y + x \cdot 137.519416416\right)\right)}{47.066876606 + x \cdot 313.399215894}\\
\mathbf{else}:\\
\;\;\;\;\left(\left(x \cdot 4.16438922228 + \frac{3655.1204654076414}{x}\right) + \frac{y - 130977.50649958357}{x \cdot x}\right) + -110.1139242984811\\
\end{array}
\end{array}
if x < -35Initial program 16.7%
associate-*r/21.0%
sub-neg21.0%
metadata-eval21.0%
*-commutative21.0%
fma-def21.0%
*-commutative21.0%
fma-def21.0%
*-commutative21.0%
fma-def21.0%
fma-def21.0%
*-commutative21.0%
Simplified20.9%
Taylor expanded in z around 0 21.0%
Taylor expanded in x around inf 95.2%
Taylor expanded in x around inf 41.8%
cube-mult41.8%
unpow241.8%
distribute-rgt-out94.3%
unpow294.3%
associate-*r*94.3%
+-commutative94.3%
Simplified94.3%
if -35 < x < 56Initial program 99.6%
Taylor expanded in x around 0 99.6%
*-commutative99.6%
Simplified99.6%
Taylor expanded in x around 0 96.7%
*-commutative96.7%
Simplified96.7%
if 56 < x Initial program 18.1%
associate-*r/23.4%
sub-neg23.4%
metadata-eval23.4%
*-commutative23.4%
fma-def23.4%
*-commutative23.4%
fma-def23.4%
*-commutative23.4%
fma-def23.4%
fma-def23.4%
*-commutative23.4%
Simplified23.4%
Taylor expanded in x around -inf 94.8%
sub-neg94.8%
+-commutative94.8%
mul-1-neg94.8%
unsub-neg94.8%
*-commutative94.8%
fma-def94.8%
associate-*r/94.8%
metadata-eval94.8%
mul-1-neg94.8%
unsub-neg94.8%
unpow294.8%
metadata-eval94.8%
Simplified94.8%
fma-udef94.8%
*-commutative94.8%
Applied egg-rr94.8%
Final simplification95.6%
(FPCore (x y z)
:precision binary64
(let* ((t_0 (/ (+ x -2.0) (+ (/ 5.86923874282773 x) 0.24013125253755718)))
(t_1 (* (+ x -2.0) (/ z (+ 47.066876606 (* x 313.399215894))))))
(if (<= x -0.175)
t_0
(if (<= x -2.6e-135)
(* (* x y) -0.0424927283095952)
(if (<= x -5.2e-191)
t_1
(if (<= x -3.1e-212)
(* x (* y -0.0424927283095952))
(if (<= x 3.7) t_1 t_0)))))))
double code(double x, double y, double z) {
double t_0 = (x + -2.0) / ((5.86923874282773 / x) + 0.24013125253755718);
double t_1 = (x + -2.0) * (z / (47.066876606 + (x * 313.399215894)));
double tmp;
if (x <= -0.175) {
tmp = t_0;
} else if (x <= -2.6e-135) {
tmp = (x * y) * -0.0424927283095952;
} else if (x <= -5.2e-191) {
tmp = t_1;
} else if (x <= -3.1e-212) {
tmp = x * (y * -0.0424927283095952);
} else if (x <= 3.7) {
tmp = t_1;
} else {
tmp = t_0;
}
return tmp;
}
real(8) function code(x, y, z)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8) :: t_0
real(8) :: t_1
real(8) :: tmp
t_0 = (x + (-2.0d0)) / ((5.86923874282773d0 / x) + 0.24013125253755718d0)
t_1 = (x + (-2.0d0)) * (z / (47.066876606d0 + (x * 313.399215894d0)))
if (x <= (-0.175d0)) then
tmp = t_0
else if (x <= (-2.6d-135)) then
tmp = (x * y) * (-0.0424927283095952d0)
else if (x <= (-5.2d-191)) then
tmp = t_1
else if (x <= (-3.1d-212)) then
tmp = x * (y * (-0.0424927283095952d0))
else if (x <= 3.7d0) then
tmp = t_1
else
tmp = t_0
end if
code = tmp
end function
public static double code(double x, double y, double z) {
double t_0 = (x + -2.0) / ((5.86923874282773 / x) + 0.24013125253755718);
double t_1 = (x + -2.0) * (z / (47.066876606 + (x * 313.399215894)));
double tmp;
if (x <= -0.175) {
tmp = t_0;
} else if (x <= -2.6e-135) {
tmp = (x * y) * -0.0424927283095952;
} else if (x <= -5.2e-191) {
tmp = t_1;
} else if (x <= -3.1e-212) {
tmp = x * (y * -0.0424927283095952);
} else if (x <= 3.7) {
tmp = t_1;
} else {
tmp = t_0;
}
return tmp;
}
def code(x, y, z): t_0 = (x + -2.0) / ((5.86923874282773 / x) + 0.24013125253755718) t_1 = (x + -2.0) * (z / (47.066876606 + (x * 313.399215894))) tmp = 0 if x <= -0.175: tmp = t_0 elif x <= -2.6e-135: tmp = (x * y) * -0.0424927283095952 elif x <= -5.2e-191: tmp = t_1 elif x <= -3.1e-212: tmp = x * (y * -0.0424927283095952) elif x <= 3.7: tmp = t_1 else: tmp = t_0 return tmp
function code(x, y, z) t_0 = Float64(Float64(x + -2.0) / Float64(Float64(5.86923874282773 / x) + 0.24013125253755718)) t_1 = Float64(Float64(x + -2.0) * Float64(z / Float64(47.066876606 + Float64(x * 313.399215894)))) tmp = 0.0 if (x <= -0.175) tmp = t_0; elseif (x <= -2.6e-135) tmp = Float64(Float64(x * y) * -0.0424927283095952); elseif (x <= -5.2e-191) tmp = t_1; elseif (x <= -3.1e-212) tmp = Float64(x * Float64(y * -0.0424927283095952)); elseif (x <= 3.7) tmp = t_1; else tmp = t_0; end return tmp end
function tmp_2 = code(x, y, z) t_0 = (x + -2.0) / ((5.86923874282773 / x) + 0.24013125253755718); t_1 = (x + -2.0) * (z / (47.066876606 + (x * 313.399215894))); tmp = 0.0; if (x <= -0.175) tmp = t_0; elseif (x <= -2.6e-135) tmp = (x * y) * -0.0424927283095952; elseif (x <= -5.2e-191) tmp = t_1; elseif (x <= -3.1e-212) tmp = x * (y * -0.0424927283095952); elseif (x <= 3.7) tmp = t_1; else tmp = t_0; end tmp_2 = tmp; end
code[x_, y_, z_] := Block[{t$95$0 = N[(N[(x + -2.0), $MachinePrecision] / N[(N[(5.86923874282773 / x), $MachinePrecision] + 0.24013125253755718), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(N[(x + -2.0), $MachinePrecision] * N[(z / N[(47.066876606 + N[(x * 313.399215894), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[x, -0.175], t$95$0, If[LessEqual[x, -2.6e-135], N[(N[(x * y), $MachinePrecision] * -0.0424927283095952), $MachinePrecision], If[LessEqual[x, -5.2e-191], t$95$1, If[LessEqual[x, -3.1e-212], N[(x * N[(y * -0.0424927283095952), $MachinePrecision]), $MachinePrecision], If[LessEqual[x, 3.7], t$95$1, t$95$0]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \frac{x + -2}{\frac{5.86923874282773}{x} + 0.24013125253755718}\\
t_1 := \left(x + -2\right) \cdot \frac{z}{47.066876606 + x \cdot 313.399215894}\\
\mathbf{if}\;x \leq -0.175:\\
\;\;\;\;t_0\\
\mathbf{elif}\;x \leq -2.6 \cdot 10^{-135}:\\
\;\;\;\;\left(x \cdot y\right) \cdot -0.0424927283095952\\
\mathbf{elif}\;x \leq -5.2 \cdot 10^{-191}:\\
\;\;\;\;t_1\\
\mathbf{elif}\;x \leq -3.1 \cdot 10^{-212}:\\
\;\;\;\;x \cdot \left(y \cdot -0.0424927283095952\right)\\
\mathbf{elif}\;x \leq 3.7:\\
\;\;\;\;t_1\\
\mathbf{else}:\\
\;\;\;\;t_0\\
\end{array}
\end{array}
if x < -0.17499999999999999 or 3.7000000000000002 < x Initial program 18.0%
associate-/l*22.8%
sub-neg22.8%
metadata-eval22.8%
fma-def22.8%
fma-def22.8%
fma-def22.8%
fma-def22.8%
fma-def22.8%
fma-def22.8%
fma-def22.8%
Simplified22.8%
Taylor expanded in x around inf 88.6%
associate-*r/88.6%
metadata-eval88.6%
Simplified88.6%
if -0.17499999999999999 < x < -2.60000000000000004e-135Initial program 99.5%
*-commutative99.5%
associate-*r/99.3%
*-commutative99.3%
fma-def99.3%
*-commutative99.3%
fma-def99.3%
*-commutative99.3%
fma-def99.3%
fma-def99.3%
Simplified99.3%
Taylor expanded in z around 0 65.4%
Taylor expanded in x around 0 53.3%
if -2.60000000000000004e-135 < x < -5.19999999999999972e-191 or -3.10000000000000006e-212 < x < 3.7000000000000002Initial program 99.6%
associate-*r/99.7%
sub-neg99.7%
metadata-eval99.7%
*-commutative99.7%
fma-def99.7%
*-commutative99.7%
fma-def99.7%
*-commutative99.7%
fma-def99.7%
fma-def99.7%
*-commutative99.7%
Simplified99.7%
Taylor expanded in z around inf 68.3%
Taylor expanded in x around 0 66.8%
*-commutative66.8%
Simplified66.8%
if -5.19999999999999972e-191 < x < -3.10000000000000006e-212Initial program 99.5%
*-commutative99.5%
associate-*r/99.7%
*-commutative99.7%
fma-def99.7%
*-commutative99.7%
fma-def99.7%
*-commutative99.7%
fma-def99.7%
fma-def99.7%
Simplified99.7%
Taylor expanded in z around 0 70.4%
Taylor expanded in x around 0 70.7%
associate-*r*70.9%
*-commutative70.9%
Simplified70.9%
Final simplification77.7%
(FPCore (x y z)
:precision binary64
(if (or (<= x -7500.0) (not (<= x 620.0)))
(/
(+ x -2.0)
(+
(/ 5.86923874282773 x)
(- 0.24013125253755718 (/ 55.572073733743466 (* x x)))))
(*
(+ x -2.0)
(-
(* z 0.0212463641547976)
(* x (- (* z 0.14147091005106402) (* y 0.0212463641547976)))))))
double code(double x, double y, double z) {
double tmp;
if ((x <= -7500.0) || !(x <= 620.0)) {
tmp = (x + -2.0) / ((5.86923874282773 / x) + (0.24013125253755718 - (55.572073733743466 / (x * x))));
} else {
tmp = (x + -2.0) * ((z * 0.0212463641547976) - (x * ((z * 0.14147091005106402) - (y * 0.0212463641547976))));
}
return tmp;
}
real(8) function code(x, y, z)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8) :: tmp
if ((x <= (-7500.0d0)) .or. (.not. (x <= 620.0d0))) then
tmp = (x + (-2.0d0)) / ((5.86923874282773d0 / x) + (0.24013125253755718d0 - (55.572073733743466d0 / (x * x))))
else
tmp = (x + (-2.0d0)) * ((z * 0.0212463641547976d0) - (x * ((z * 0.14147091005106402d0) - (y * 0.0212463641547976d0))))
end if
code = tmp
end function
public static double code(double x, double y, double z) {
double tmp;
if ((x <= -7500.0) || !(x <= 620.0)) {
tmp = (x + -2.0) / ((5.86923874282773 / x) + (0.24013125253755718 - (55.572073733743466 / (x * x))));
} else {
tmp = (x + -2.0) * ((z * 0.0212463641547976) - (x * ((z * 0.14147091005106402) - (y * 0.0212463641547976))));
}
return tmp;
}
def code(x, y, z): tmp = 0 if (x <= -7500.0) or not (x <= 620.0): tmp = (x + -2.0) / ((5.86923874282773 / x) + (0.24013125253755718 - (55.572073733743466 / (x * x)))) else: tmp = (x + -2.0) * ((z * 0.0212463641547976) - (x * ((z * 0.14147091005106402) - (y * 0.0212463641547976)))) return tmp
function code(x, y, z) tmp = 0.0 if ((x <= -7500.0) || !(x <= 620.0)) tmp = Float64(Float64(x + -2.0) / Float64(Float64(5.86923874282773 / x) + Float64(0.24013125253755718 - Float64(55.572073733743466 / Float64(x * x))))); else tmp = Float64(Float64(x + -2.0) * Float64(Float64(z * 0.0212463641547976) - Float64(x * Float64(Float64(z * 0.14147091005106402) - Float64(y * 0.0212463641547976))))); end return tmp end
function tmp_2 = code(x, y, z) tmp = 0.0; if ((x <= -7500.0) || ~((x <= 620.0))) tmp = (x + -2.0) / ((5.86923874282773 / x) + (0.24013125253755718 - (55.572073733743466 / (x * x)))); else tmp = (x + -2.0) * ((z * 0.0212463641547976) - (x * ((z * 0.14147091005106402) - (y * 0.0212463641547976)))); end tmp_2 = tmp; end
code[x_, y_, z_] := If[Or[LessEqual[x, -7500.0], N[Not[LessEqual[x, 620.0]], $MachinePrecision]], N[(N[(x + -2.0), $MachinePrecision] / N[(N[(5.86923874282773 / x), $MachinePrecision] + N[(0.24013125253755718 - N[(55.572073733743466 / N[(x * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(x + -2.0), $MachinePrecision] * N[(N[(z * 0.0212463641547976), $MachinePrecision] - N[(x * N[(N[(z * 0.14147091005106402), $MachinePrecision] - N[(y * 0.0212463641547976), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq -7500 \lor \neg \left(x \leq 620\right):\\
\;\;\;\;\frac{x + -2}{\frac{5.86923874282773}{x} + \left(0.24013125253755718 - \frac{55.572073733743466}{x \cdot x}\right)}\\
\mathbf{else}:\\
\;\;\;\;\left(x + -2\right) \cdot \left(z \cdot 0.0212463641547976 - x \cdot \left(z \cdot 0.14147091005106402 - y \cdot 0.0212463641547976\right)\right)\\
\end{array}
\end{array}
if x < -7500 or 620 < x Initial program 16.2%
associate-/l*21.1%
sub-neg21.1%
metadata-eval21.1%
fma-def21.1%
fma-def21.1%
fma-def21.1%
fma-def21.1%
fma-def21.1%
fma-def21.1%
fma-def21.1%
Simplified21.1%
Taylor expanded in x around inf 90.7%
+-commutative90.7%
associate--l+90.7%
associate-*r/90.7%
metadata-eval90.7%
associate-*r/90.7%
metadata-eval90.7%
unpow290.7%
Simplified90.7%
if -7500 < x < 620Initial program 99.6%
associate-*r/99.6%
sub-neg99.6%
metadata-eval99.6%
*-commutative99.6%
fma-def99.6%
*-commutative99.6%
fma-def99.6%
*-commutative99.6%
fma-def99.6%
fma-def99.6%
*-commutative99.6%
Simplified99.6%
Taylor expanded in x around 0 91.5%
Final simplification91.1%
(FPCore (x y z)
:precision binary64
(if (or (<= x -2600.0) (not (<= x 1.15)))
(+
(+
(+ (* x 4.16438922228) (/ 3655.1204654076414 x))
(/ (- y 130977.50649958357) (* x x)))
-110.1139242984811)
(*
(+ x -2.0)
(-
(* z 0.0212463641547976)
(* x (- (* z 0.14147091005106402) (* y 0.0212463641547976)))))))
double code(double x, double y, double z) {
double tmp;
if ((x <= -2600.0) || !(x <= 1.15)) {
tmp = (((x * 4.16438922228) + (3655.1204654076414 / x)) + ((y - 130977.50649958357) / (x * x))) + -110.1139242984811;
} else {
tmp = (x + -2.0) * ((z * 0.0212463641547976) - (x * ((z * 0.14147091005106402) - (y * 0.0212463641547976))));
}
return tmp;
}
real(8) function code(x, y, z)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8) :: tmp
if ((x <= (-2600.0d0)) .or. (.not. (x <= 1.15d0))) then
tmp = (((x * 4.16438922228d0) + (3655.1204654076414d0 / x)) + ((y - 130977.50649958357d0) / (x * x))) + (-110.1139242984811d0)
else
tmp = (x + (-2.0d0)) * ((z * 0.0212463641547976d0) - (x * ((z * 0.14147091005106402d0) - (y * 0.0212463641547976d0))))
end if
code = tmp
end function
public static double code(double x, double y, double z) {
double tmp;
if ((x <= -2600.0) || !(x <= 1.15)) {
tmp = (((x * 4.16438922228) + (3655.1204654076414 / x)) + ((y - 130977.50649958357) / (x * x))) + -110.1139242984811;
} else {
tmp = (x + -2.0) * ((z * 0.0212463641547976) - (x * ((z * 0.14147091005106402) - (y * 0.0212463641547976))));
}
return tmp;
}
def code(x, y, z): tmp = 0 if (x <= -2600.0) or not (x <= 1.15): tmp = (((x * 4.16438922228) + (3655.1204654076414 / x)) + ((y - 130977.50649958357) / (x * x))) + -110.1139242984811 else: tmp = (x + -2.0) * ((z * 0.0212463641547976) - (x * ((z * 0.14147091005106402) - (y * 0.0212463641547976)))) return tmp
function code(x, y, z) tmp = 0.0 if ((x <= -2600.0) || !(x <= 1.15)) tmp = Float64(Float64(Float64(Float64(x * 4.16438922228) + Float64(3655.1204654076414 / x)) + Float64(Float64(y - 130977.50649958357) / Float64(x * x))) + -110.1139242984811); else tmp = Float64(Float64(x + -2.0) * Float64(Float64(z * 0.0212463641547976) - Float64(x * Float64(Float64(z * 0.14147091005106402) - Float64(y * 0.0212463641547976))))); end return tmp end
function tmp_2 = code(x, y, z) tmp = 0.0; if ((x <= -2600.0) || ~((x <= 1.15))) tmp = (((x * 4.16438922228) + (3655.1204654076414 / x)) + ((y - 130977.50649958357) / (x * x))) + -110.1139242984811; else tmp = (x + -2.0) * ((z * 0.0212463641547976) - (x * ((z * 0.14147091005106402) - (y * 0.0212463641547976)))); end tmp_2 = tmp; end
code[x_, y_, z_] := If[Or[LessEqual[x, -2600.0], N[Not[LessEqual[x, 1.15]], $MachinePrecision]], N[(N[(N[(N[(x * 4.16438922228), $MachinePrecision] + N[(3655.1204654076414 / x), $MachinePrecision]), $MachinePrecision] + N[(N[(y - 130977.50649958357), $MachinePrecision] / N[(x * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + -110.1139242984811), $MachinePrecision], N[(N[(x + -2.0), $MachinePrecision] * N[(N[(z * 0.0212463641547976), $MachinePrecision] - N[(x * N[(N[(z * 0.14147091005106402), $MachinePrecision] - N[(y * 0.0212463641547976), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq -2600 \lor \neg \left(x \leq 1.15\right):\\
\;\;\;\;\left(\left(x \cdot 4.16438922228 + \frac{3655.1204654076414}{x}\right) + \frac{y - 130977.50649958357}{x \cdot x}\right) + -110.1139242984811\\
\mathbf{else}:\\
\;\;\;\;\left(x + -2\right) \cdot \left(z \cdot 0.0212463641547976 - x \cdot \left(z \cdot 0.14147091005106402 - y \cdot 0.0212463641547976\right)\right)\\
\end{array}
\end{array}
if x < -2600 or 1.1499999999999999 < x Initial program 18.0%
associate-*r/22.8%
sub-neg22.8%
metadata-eval22.8%
*-commutative22.8%
fma-def22.8%
*-commutative22.8%
fma-def22.8%
*-commutative22.8%
fma-def22.8%
fma-def22.8%
*-commutative22.8%
Simplified22.8%
Taylor expanded in x around -inf 93.2%
sub-neg93.2%
+-commutative93.2%
mul-1-neg93.2%
unsub-neg93.2%
*-commutative93.2%
fma-def93.2%
associate-*r/93.2%
metadata-eval93.2%
mul-1-neg93.2%
unsub-neg93.2%
unpow293.2%
metadata-eval93.2%
Simplified93.2%
fma-udef93.2%
*-commutative93.2%
Applied egg-rr93.2%
if -2600 < x < 1.1499999999999999Initial program 99.6%
associate-*r/99.6%
sub-neg99.6%
metadata-eval99.6%
*-commutative99.6%
fma-def99.6%
*-commutative99.6%
fma-def99.6%
*-commutative99.6%
fma-def99.6%
fma-def99.6%
*-commutative99.6%
Simplified99.6%
Taylor expanded in x around 0 93.6%
Final simplification93.4%
(FPCore (x y z)
:precision binary64
(if (<= x -0.175)
(*
(+ x -2.0)
(+
4.16438922228
(/ z (+ 47.066876606 (* x (+ 313.399215894 (* x (* x x))))))))
(if (<= x 1.15)
(*
(+ x -2.0)
(-
(* z 0.0212463641547976)
(* x (- (* z 0.14147091005106402) (* y 0.0212463641547976)))))
(+
(+
(+ (* x 4.16438922228) (/ 3655.1204654076414 x))
(/ (- y 130977.50649958357) (* x x)))
-110.1139242984811))))
double code(double x, double y, double z) {
double tmp;
if (x <= -0.175) {
tmp = (x + -2.0) * (4.16438922228 + (z / (47.066876606 + (x * (313.399215894 + (x * (x * x)))))));
} else if (x <= 1.15) {
tmp = (x + -2.0) * ((z * 0.0212463641547976) - (x * ((z * 0.14147091005106402) - (y * 0.0212463641547976))));
} else {
tmp = (((x * 4.16438922228) + (3655.1204654076414 / x)) + ((y - 130977.50649958357) / (x * x))) + -110.1139242984811;
}
return tmp;
}
real(8) function code(x, y, z)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8) :: tmp
if (x <= (-0.175d0)) then
tmp = (x + (-2.0d0)) * (4.16438922228d0 + (z / (47.066876606d0 + (x * (313.399215894d0 + (x * (x * x)))))))
else if (x <= 1.15d0) then
tmp = (x + (-2.0d0)) * ((z * 0.0212463641547976d0) - (x * ((z * 0.14147091005106402d0) - (y * 0.0212463641547976d0))))
else
tmp = (((x * 4.16438922228d0) + (3655.1204654076414d0 / x)) + ((y - 130977.50649958357d0) / (x * x))) + (-110.1139242984811d0)
end if
code = tmp
end function
public static double code(double x, double y, double z) {
double tmp;
if (x <= -0.175) {
tmp = (x + -2.0) * (4.16438922228 + (z / (47.066876606 + (x * (313.399215894 + (x * (x * x)))))));
} else if (x <= 1.15) {
tmp = (x + -2.0) * ((z * 0.0212463641547976) - (x * ((z * 0.14147091005106402) - (y * 0.0212463641547976))));
} else {
tmp = (((x * 4.16438922228) + (3655.1204654076414 / x)) + ((y - 130977.50649958357) / (x * x))) + -110.1139242984811;
}
return tmp;
}
def code(x, y, z): tmp = 0 if x <= -0.175: tmp = (x + -2.0) * (4.16438922228 + (z / (47.066876606 + (x * (313.399215894 + (x * (x * x))))))) elif x <= 1.15: tmp = (x + -2.0) * ((z * 0.0212463641547976) - (x * ((z * 0.14147091005106402) - (y * 0.0212463641547976)))) else: tmp = (((x * 4.16438922228) + (3655.1204654076414 / x)) + ((y - 130977.50649958357) / (x * x))) + -110.1139242984811 return tmp
function code(x, y, z) tmp = 0.0 if (x <= -0.175) tmp = Float64(Float64(x + -2.0) * Float64(4.16438922228 + Float64(z / Float64(47.066876606 + Float64(x * Float64(313.399215894 + Float64(x * Float64(x * x)))))))); elseif (x <= 1.15) tmp = Float64(Float64(x + -2.0) * Float64(Float64(z * 0.0212463641547976) - Float64(x * Float64(Float64(z * 0.14147091005106402) - Float64(y * 0.0212463641547976))))); else tmp = Float64(Float64(Float64(Float64(x * 4.16438922228) + Float64(3655.1204654076414 / x)) + Float64(Float64(y - 130977.50649958357) / Float64(x * x))) + -110.1139242984811); end return tmp end
function tmp_2 = code(x, y, z) tmp = 0.0; if (x <= -0.175) tmp = (x + -2.0) * (4.16438922228 + (z / (47.066876606 + (x * (313.399215894 + (x * (x * x))))))); elseif (x <= 1.15) tmp = (x + -2.0) * ((z * 0.0212463641547976) - (x * ((z * 0.14147091005106402) - (y * 0.0212463641547976)))); else tmp = (((x * 4.16438922228) + (3655.1204654076414 / x)) + ((y - 130977.50649958357) / (x * x))) + -110.1139242984811; end tmp_2 = tmp; end
code[x_, y_, z_] := If[LessEqual[x, -0.175], N[(N[(x + -2.0), $MachinePrecision] * N[(4.16438922228 + N[(z / N[(47.066876606 + N[(x * N[(313.399215894 + N[(x * N[(x * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[x, 1.15], N[(N[(x + -2.0), $MachinePrecision] * N[(N[(z * 0.0212463641547976), $MachinePrecision] - N[(x * N[(N[(z * 0.14147091005106402), $MachinePrecision] - N[(y * 0.0212463641547976), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(N[(N[(x * 4.16438922228), $MachinePrecision] + N[(3655.1204654076414 / x), $MachinePrecision]), $MachinePrecision] + N[(N[(y - 130977.50649958357), $MachinePrecision] / N[(x * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + -110.1139242984811), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq -0.175:\\
\;\;\;\;\left(x + -2\right) \cdot \left(4.16438922228 + \frac{z}{47.066876606 + x \cdot \left(313.399215894 + x \cdot \left(x \cdot x\right)\right)}\right)\\
\mathbf{elif}\;x \leq 1.15:\\
\;\;\;\;\left(x + -2\right) \cdot \left(z \cdot 0.0212463641547976 - x \cdot \left(z \cdot 0.14147091005106402 - y \cdot 0.0212463641547976\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\left(\left(x \cdot 4.16438922228 + \frac{3655.1204654076414}{x}\right) + \frac{y - 130977.50649958357}{x \cdot x}\right) + -110.1139242984811\\
\end{array}
\end{array}
if x < -0.17499999999999999Initial program 16.7%
associate-*r/21.0%
sub-neg21.0%
metadata-eval21.0%
*-commutative21.0%
fma-def21.0%
*-commutative21.0%
fma-def21.0%
*-commutative21.0%
fma-def21.0%
fma-def21.0%
*-commutative21.0%
Simplified20.9%
Taylor expanded in z around 0 21.0%
Taylor expanded in x around inf 95.2%
Taylor expanded in x around inf 41.8%
cube-mult41.8%
unpow241.8%
distribute-rgt-out94.3%
unpow294.3%
associate-*r*94.3%
+-commutative94.3%
Simplified94.3%
Taylor expanded in x around inf 94.1%
unpow294.1%
Simplified94.1%
if -0.17499999999999999 < x < 1.1499999999999999Initial program 99.6%
associate-*r/99.6%
sub-neg99.6%
metadata-eval99.6%
*-commutative99.6%
fma-def99.6%
*-commutative99.6%
fma-def99.6%
*-commutative99.6%
fma-def99.6%
fma-def99.6%
*-commutative99.6%
Simplified99.6%
Taylor expanded in x around 0 94.4%
if 1.1499999999999999 < x Initial program 20.3%
associate-*r/25.4%
sub-neg25.4%
metadata-eval25.4%
*-commutative25.4%
fma-def25.4%
*-commutative25.4%
fma-def25.4%
*-commutative25.4%
fma-def25.4%
fma-def25.4%
*-commutative25.4%
Simplified25.4%
Taylor expanded in x around -inf 92.5%
sub-neg92.5%
+-commutative92.5%
mul-1-neg92.5%
unsub-neg92.5%
*-commutative92.5%
fma-def92.5%
associate-*r/92.5%
metadata-eval92.5%
mul-1-neg92.5%
unsub-neg92.5%
unpow292.5%
metadata-eval92.5%
Simplified92.5%
fma-udef92.5%
*-commutative92.5%
Applied egg-rr92.5%
Final simplification93.7%
(FPCore (x y z)
:precision binary64
(let* ((t_0 (/ (+ x -2.0) (+ (/ 5.86923874282773 x) 0.24013125253755718))))
(if (<= x -0.175)
t_0
(if (<= x -2.6e-135)
(* (* x y) -0.0424927283095952)
(if (<= x -5.2e-191)
(* z -0.0424927283095952)
(if (<= x -3.1e-212)
(* x (* y -0.0424927283095952))
(if (<= x 5.2) (/ (+ x -2.0) (/ 47.066876606 z)) t_0)))))))
double code(double x, double y, double z) {
double t_0 = (x + -2.0) / ((5.86923874282773 / x) + 0.24013125253755718);
double tmp;
if (x <= -0.175) {
tmp = t_0;
} else if (x <= -2.6e-135) {
tmp = (x * y) * -0.0424927283095952;
} else if (x <= -5.2e-191) {
tmp = z * -0.0424927283095952;
} else if (x <= -3.1e-212) {
tmp = x * (y * -0.0424927283095952);
} else if (x <= 5.2) {
tmp = (x + -2.0) / (47.066876606 / z);
} else {
tmp = t_0;
}
return tmp;
}
real(8) function code(x, y, z)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8) :: t_0
real(8) :: tmp
t_0 = (x + (-2.0d0)) / ((5.86923874282773d0 / x) + 0.24013125253755718d0)
if (x <= (-0.175d0)) then
tmp = t_0
else if (x <= (-2.6d-135)) then
tmp = (x * y) * (-0.0424927283095952d0)
else if (x <= (-5.2d-191)) then
tmp = z * (-0.0424927283095952d0)
else if (x <= (-3.1d-212)) then
tmp = x * (y * (-0.0424927283095952d0))
else if (x <= 5.2d0) then
tmp = (x + (-2.0d0)) / (47.066876606d0 / z)
else
tmp = t_0
end if
code = tmp
end function
public static double code(double x, double y, double z) {
double t_0 = (x + -2.0) / ((5.86923874282773 / x) + 0.24013125253755718);
double tmp;
if (x <= -0.175) {
tmp = t_0;
} else if (x <= -2.6e-135) {
tmp = (x * y) * -0.0424927283095952;
} else if (x <= -5.2e-191) {
tmp = z * -0.0424927283095952;
} else if (x <= -3.1e-212) {
tmp = x * (y * -0.0424927283095952);
} else if (x <= 5.2) {
tmp = (x + -2.0) / (47.066876606 / z);
} else {
tmp = t_0;
}
return tmp;
}
def code(x, y, z): t_0 = (x + -2.0) / ((5.86923874282773 / x) + 0.24013125253755718) tmp = 0 if x <= -0.175: tmp = t_0 elif x <= -2.6e-135: tmp = (x * y) * -0.0424927283095952 elif x <= -5.2e-191: tmp = z * -0.0424927283095952 elif x <= -3.1e-212: tmp = x * (y * -0.0424927283095952) elif x <= 5.2: tmp = (x + -2.0) / (47.066876606 / z) else: tmp = t_0 return tmp
function code(x, y, z) t_0 = Float64(Float64(x + -2.0) / Float64(Float64(5.86923874282773 / x) + 0.24013125253755718)) tmp = 0.0 if (x <= -0.175) tmp = t_0; elseif (x <= -2.6e-135) tmp = Float64(Float64(x * y) * -0.0424927283095952); elseif (x <= -5.2e-191) tmp = Float64(z * -0.0424927283095952); elseif (x <= -3.1e-212) tmp = Float64(x * Float64(y * -0.0424927283095952)); elseif (x <= 5.2) tmp = Float64(Float64(x + -2.0) / Float64(47.066876606 / z)); else tmp = t_0; end return tmp end
function tmp_2 = code(x, y, z) t_0 = (x + -2.0) / ((5.86923874282773 / x) + 0.24013125253755718); tmp = 0.0; if (x <= -0.175) tmp = t_0; elseif (x <= -2.6e-135) tmp = (x * y) * -0.0424927283095952; elseif (x <= -5.2e-191) tmp = z * -0.0424927283095952; elseif (x <= -3.1e-212) tmp = x * (y * -0.0424927283095952); elseif (x <= 5.2) tmp = (x + -2.0) / (47.066876606 / z); else tmp = t_0; end tmp_2 = tmp; end
code[x_, y_, z_] := Block[{t$95$0 = N[(N[(x + -2.0), $MachinePrecision] / N[(N[(5.86923874282773 / x), $MachinePrecision] + 0.24013125253755718), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[x, -0.175], t$95$0, If[LessEqual[x, -2.6e-135], N[(N[(x * y), $MachinePrecision] * -0.0424927283095952), $MachinePrecision], If[LessEqual[x, -5.2e-191], N[(z * -0.0424927283095952), $MachinePrecision], If[LessEqual[x, -3.1e-212], N[(x * N[(y * -0.0424927283095952), $MachinePrecision]), $MachinePrecision], If[LessEqual[x, 5.2], N[(N[(x + -2.0), $MachinePrecision] / N[(47.066876606 / z), $MachinePrecision]), $MachinePrecision], t$95$0]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \frac{x + -2}{\frac{5.86923874282773}{x} + 0.24013125253755718}\\
\mathbf{if}\;x \leq -0.175:\\
\;\;\;\;t_0\\
\mathbf{elif}\;x \leq -2.6 \cdot 10^{-135}:\\
\;\;\;\;\left(x \cdot y\right) \cdot -0.0424927283095952\\
\mathbf{elif}\;x \leq -5.2 \cdot 10^{-191}:\\
\;\;\;\;z \cdot -0.0424927283095952\\
\mathbf{elif}\;x \leq -3.1 \cdot 10^{-212}:\\
\;\;\;\;x \cdot \left(y \cdot -0.0424927283095952\right)\\
\mathbf{elif}\;x \leq 5.2:\\
\;\;\;\;\frac{x + -2}{\frac{47.066876606}{z}}\\
\mathbf{else}:\\
\;\;\;\;t_0\\
\end{array}
\end{array}
if x < -0.17499999999999999 or 5.20000000000000018 < x Initial program 18.0%
associate-/l*22.8%
sub-neg22.8%
metadata-eval22.8%
fma-def22.8%
fma-def22.8%
fma-def22.8%
fma-def22.8%
fma-def22.8%
fma-def22.8%
fma-def22.8%
Simplified22.8%
Taylor expanded in x around inf 88.6%
associate-*r/88.6%
metadata-eval88.6%
Simplified88.6%
if -0.17499999999999999 < x < -2.60000000000000004e-135Initial program 99.5%
*-commutative99.5%
associate-*r/99.3%
*-commutative99.3%
fma-def99.3%
*-commutative99.3%
fma-def99.3%
*-commutative99.3%
fma-def99.3%
fma-def99.3%
Simplified99.3%
Taylor expanded in z around 0 65.4%
Taylor expanded in x around 0 53.3%
if -2.60000000000000004e-135 < x < -5.19999999999999972e-191Initial program 99.8%
associate-*r/99.8%
sub-neg99.8%
metadata-eval99.8%
*-commutative99.8%
fma-def99.8%
*-commutative99.8%
fma-def99.8%
*-commutative99.8%
fma-def99.8%
fma-def99.8%
*-commutative99.8%
Simplified99.8%
Taylor expanded in x around 0 87.9%
*-commutative87.9%
Simplified87.9%
if -5.19999999999999972e-191 < x < -3.10000000000000006e-212Initial program 99.5%
*-commutative99.5%
associate-*r/99.7%
*-commutative99.7%
fma-def99.7%
*-commutative99.7%
fma-def99.7%
*-commutative99.7%
fma-def99.7%
fma-def99.7%
Simplified99.7%
Taylor expanded in z around 0 70.4%
Taylor expanded in x around 0 70.7%
associate-*r*70.9%
*-commutative70.9%
Simplified70.9%
if -3.10000000000000006e-212 < x < 5.20000000000000018Initial program 99.6%
associate-/l*99.5%
sub-neg99.5%
metadata-eval99.5%
fma-def99.5%
fma-def99.5%
fma-def99.5%
fma-def99.5%
fma-def99.5%
fma-def99.5%
fma-def99.5%
Simplified99.5%
Taylor expanded in x around 0 64.2%
Final simplification77.5%
(FPCore (x y z)
:precision binary64
(if (or (<= x -3000.0) (not (<= x 620.0)))
(/
(+ x -2.0)
(+
(/ 5.86923874282773 x)
(- 0.24013125253755718 (/ 55.572073733743466 (* x x)))))
(* (+ x -2.0) (+ (* z 0.0212463641547976) (* 0.0212463641547976 (* x y))))))
double code(double x, double y, double z) {
double tmp;
if ((x <= -3000.0) || !(x <= 620.0)) {
tmp = (x + -2.0) / ((5.86923874282773 / x) + (0.24013125253755718 - (55.572073733743466 / (x * x))));
} else {
tmp = (x + -2.0) * ((z * 0.0212463641547976) + (0.0212463641547976 * (x * y)));
}
return tmp;
}
real(8) function code(x, y, z)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8) :: tmp
if ((x <= (-3000.0d0)) .or. (.not. (x <= 620.0d0))) then
tmp = (x + (-2.0d0)) / ((5.86923874282773d0 / x) + (0.24013125253755718d0 - (55.572073733743466d0 / (x * x))))
else
tmp = (x + (-2.0d0)) * ((z * 0.0212463641547976d0) + (0.0212463641547976d0 * (x * y)))
end if
code = tmp
end function
public static double code(double x, double y, double z) {
double tmp;
if ((x <= -3000.0) || !(x <= 620.0)) {
tmp = (x + -2.0) / ((5.86923874282773 / x) + (0.24013125253755718 - (55.572073733743466 / (x * x))));
} else {
tmp = (x + -2.0) * ((z * 0.0212463641547976) + (0.0212463641547976 * (x * y)));
}
return tmp;
}
def code(x, y, z): tmp = 0 if (x <= -3000.0) or not (x <= 620.0): tmp = (x + -2.0) / ((5.86923874282773 / x) + (0.24013125253755718 - (55.572073733743466 / (x * x)))) else: tmp = (x + -2.0) * ((z * 0.0212463641547976) + (0.0212463641547976 * (x * y))) return tmp
function code(x, y, z) tmp = 0.0 if ((x <= -3000.0) || !(x <= 620.0)) tmp = Float64(Float64(x + -2.0) / Float64(Float64(5.86923874282773 / x) + Float64(0.24013125253755718 - Float64(55.572073733743466 / Float64(x * x))))); else tmp = Float64(Float64(x + -2.0) * Float64(Float64(z * 0.0212463641547976) + Float64(0.0212463641547976 * Float64(x * y)))); end return tmp end
function tmp_2 = code(x, y, z) tmp = 0.0; if ((x <= -3000.0) || ~((x <= 620.0))) tmp = (x + -2.0) / ((5.86923874282773 / x) + (0.24013125253755718 - (55.572073733743466 / (x * x)))); else tmp = (x + -2.0) * ((z * 0.0212463641547976) + (0.0212463641547976 * (x * y))); end tmp_2 = tmp; end
code[x_, y_, z_] := If[Or[LessEqual[x, -3000.0], N[Not[LessEqual[x, 620.0]], $MachinePrecision]], N[(N[(x + -2.0), $MachinePrecision] / N[(N[(5.86923874282773 / x), $MachinePrecision] + N[(0.24013125253755718 - N[(55.572073733743466 / N[(x * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(x + -2.0), $MachinePrecision] * N[(N[(z * 0.0212463641547976), $MachinePrecision] + N[(0.0212463641547976 * N[(x * y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq -3000 \lor \neg \left(x \leq 620\right):\\
\;\;\;\;\frac{x + -2}{\frac{5.86923874282773}{x} + \left(0.24013125253755718 - \frac{55.572073733743466}{x \cdot x}\right)}\\
\mathbf{else}:\\
\;\;\;\;\left(x + -2\right) \cdot \left(z \cdot 0.0212463641547976 + 0.0212463641547976 \cdot \left(x \cdot y\right)\right)\\
\end{array}
\end{array}
if x < -3e3 or 620 < x Initial program 16.2%
associate-/l*21.1%
sub-neg21.1%
metadata-eval21.1%
fma-def21.1%
fma-def21.1%
fma-def21.1%
fma-def21.1%
fma-def21.1%
fma-def21.1%
fma-def21.1%
Simplified21.1%
Taylor expanded in x around inf 90.7%
+-commutative90.7%
associate--l+90.7%
associate-*r/90.7%
metadata-eval90.7%
associate-*r/90.7%
metadata-eval90.7%
unpow290.7%
Simplified90.7%
if -3e3 < x < 620Initial program 99.6%
associate-*r/99.6%
sub-neg99.6%
metadata-eval99.6%
*-commutative99.6%
fma-def99.6%
*-commutative99.6%
fma-def99.6%
*-commutative99.6%
fma-def99.6%
fma-def99.6%
*-commutative99.6%
Simplified99.6%
Taylor expanded in x around 0 91.5%
Taylor expanded in y around inf 91.4%
Final simplification91.0%
(FPCore (x y z)
:precision binary64
(if (<= x -0.175)
(- (* x 4.16438922228) 110.1139242984811)
(if (<= x -2.6e-135)
(* (* x y) -0.0424927283095952)
(if (<= x -6e-191)
(* z -0.0424927283095952)
(if (<= x -3.1e-212)
(* x (* y -0.0424927283095952))
(if (<= x 4.8)
(/ (+ x -2.0) (/ 47.066876606 z))
(/ (+ x -2.0) 0.24013125253755718)))))))
double code(double x, double y, double z) {
double tmp;
if (x <= -0.175) {
tmp = (x * 4.16438922228) - 110.1139242984811;
} else if (x <= -2.6e-135) {
tmp = (x * y) * -0.0424927283095952;
} else if (x <= -6e-191) {
tmp = z * -0.0424927283095952;
} else if (x <= -3.1e-212) {
tmp = x * (y * -0.0424927283095952);
} else if (x <= 4.8) {
tmp = (x + -2.0) / (47.066876606 / z);
} else {
tmp = (x + -2.0) / 0.24013125253755718;
}
return tmp;
}
real(8) function code(x, y, z)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8) :: tmp
if (x <= (-0.175d0)) then
tmp = (x * 4.16438922228d0) - 110.1139242984811d0
else if (x <= (-2.6d-135)) then
tmp = (x * y) * (-0.0424927283095952d0)
else if (x <= (-6d-191)) then
tmp = z * (-0.0424927283095952d0)
else if (x <= (-3.1d-212)) then
tmp = x * (y * (-0.0424927283095952d0))
else if (x <= 4.8d0) then
tmp = (x + (-2.0d0)) / (47.066876606d0 / z)
else
tmp = (x + (-2.0d0)) / 0.24013125253755718d0
end if
code = tmp
end function
public static double code(double x, double y, double z) {
double tmp;
if (x <= -0.175) {
tmp = (x * 4.16438922228) - 110.1139242984811;
} else if (x <= -2.6e-135) {
tmp = (x * y) * -0.0424927283095952;
} else if (x <= -6e-191) {
tmp = z * -0.0424927283095952;
} else if (x <= -3.1e-212) {
tmp = x * (y * -0.0424927283095952);
} else if (x <= 4.8) {
tmp = (x + -2.0) / (47.066876606 / z);
} else {
tmp = (x + -2.0) / 0.24013125253755718;
}
return tmp;
}
def code(x, y, z): tmp = 0 if x <= -0.175: tmp = (x * 4.16438922228) - 110.1139242984811 elif x <= -2.6e-135: tmp = (x * y) * -0.0424927283095952 elif x <= -6e-191: tmp = z * -0.0424927283095952 elif x <= -3.1e-212: tmp = x * (y * -0.0424927283095952) elif x <= 4.8: tmp = (x + -2.0) / (47.066876606 / z) else: tmp = (x + -2.0) / 0.24013125253755718 return tmp
function code(x, y, z) tmp = 0.0 if (x <= -0.175) tmp = Float64(Float64(x * 4.16438922228) - 110.1139242984811); elseif (x <= -2.6e-135) tmp = Float64(Float64(x * y) * -0.0424927283095952); elseif (x <= -6e-191) tmp = Float64(z * -0.0424927283095952); elseif (x <= -3.1e-212) tmp = Float64(x * Float64(y * -0.0424927283095952)); elseif (x <= 4.8) tmp = Float64(Float64(x + -2.0) / Float64(47.066876606 / z)); else tmp = Float64(Float64(x + -2.0) / 0.24013125253755718); end return tmp end
function tmp_2 = code(x, y, z) tmp = 0.0; if (x <= -0.175) tmp = (x * 4.16438922228) - 110.1139242984811; elseif (x <= -2.6e-135) tmp = (x * y) * -0.0424927283095952; elseif (x <= -6e-191) tmp = z * -0.0424927283095952; elseif (x <= -3.1e-212) tmp = x * (y * -0.0424927283095952); elseif (x <= 4.8) tmp = (x + -2.0) / (47.066876606 / z); else tmp = (x + -2.0) / 0.24013125253755718; end tmp_2 = tmp; end
code[x_, y_, z_] := If[LessEqual[x, -0.175], N[(N[(x * 4.16438922228), $MachinePrecision] - 110.1139242984811), $MachinePrecision], If[LessEqual[x, -2.6e-135], N[(N[(x * y), $MachinePrecision] * -0.0424927283095952), $MachinePrecision], If[LessEqual[x, -6e-191], N[(z * -0.0424927283095952), $MachinePrecision], If[LessEqual[x, -3.1e-212], N[(x * N[(y * -0.0424927283095952), $MachinePrecision]), $MachinePrecision], If[LessEqual[x, 4.8], N[(N[(x + -2.0), $MachinePrecision] / N[(47.066876606 / z), $MachinePrecision]), $MachinePrecision], N[(N[(x + -2.0), $MachinePrecision] / 0.24013125253755718), $MachinePrecision]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq -0.175:\\
\;\;\;\;x \cdot 4.16438922228 - 110.1139242984811\\
\mathbf{elif}\;x \leq -2.6 \cdot 10^{-135}:\\
\;\;\;\;\left(x \cdot y\right) \cdot -0.0424927283095952\\
\mathbf{elif}\;x \leq -6 \cdot 10^{-191}:\\
\;\;\;\;z \cdot -0.0424927283095952\\
\mathbf{elif}\;x \leq -3.1 \cdot 10^{-212}:\\
\;\;\;\;x \cdot \left(y \cdot -0.0424927283095952\right)\\
\mathbf{elif}\;x \leq 4.8:\\
\;\;\;\;\frac{x + -2}{\frac{47.066876606}{z}}\\
\mathbf{else}:\\
\;\;\;\;\frac{x + -2}{0.24013125253755718}\\
\end{array}
\end{array}
if x < -0.17499999999999999Initial program 16.7%
associate-*r/21.0%
sub-neg21.0%
metadata-eval21.0%
*-commutative21.0%
fma-def21.0%
*-commutative21.0%
fma-def21.0%
*-commutative21.0%
fma-def21.0%
fma-def21.0%
*-commutative21.0%
Simplified20.9%
Taylor expanded in x around inf 89.7%
if -0.17499999999999999 < x < -2.60000000000000004e-135Initial program 99.5%
*-commutative99.5%
associate-*r/99.3%
*-commutative99.3%
fma-def99.3%
*-commutative99.3%
fma-def99.3%
*-commutative99.3%
fma-def99.3%
fma-def99.3%
Simplified99.3%
Taylor expanded in z around 0 65.4%
Taylor expanded in x around 0 53.3%
if -2.60000000000000004e-135 < x < -6.0000000000000001e-191Initial program 99.8%
associate-*r/99.8%
sub-neg99.8%
metadata-eval99.8%
*-commutative99.8%
fma-def99.8%
*-commutative99.8%
fma-def99.8%
*-commutative99.8%
fma-def99.8%
fma-def99.8%
*-commutative99.8%
Simplified99.8%
Taylor expanded in x around 0 87.9%
*-commutative87.9%
Simplified87.9%
if -6.0000000000000001e-191 < x < -3.10000000000000006e-212Initial program 99.5%
*-commutative99.5%
associate-*r/99.7%
*-commutative99.7%
fma-def99.7%
*-commutative99.7%
fma-def99.7%
*-commutative99.7%
fma-def99.7%
fma-def99.7%
Simplified99.7%
Taylor expanded in z around 0 70.4%
Taylor expanded in x around 0 70.7%
associate-*r*70.9%
*-commutative70.9%
Simplified70.9%
if -3.10000000000000006e-212 < x < 4.79999999999999982Initial program 99.6%
associate-/l*99.5%
sub-neg99.5%
metadata-eval99.5%
fma-def99.5%
fma-def99.5%
fma-def99.5%
fma-def99.5%
fma-def99.5%
fma-def99.5%
fma-def99.5%
Simplified99.5%
Taylor expanded in x around 0 64.2%
if 4.79999999999999982 < x Initial program 19.2%
associate-/l*24.5%
sub-neg24.5%
metadata-eval24.5%
fma-def24.5%
fma-def24.5%
fma-def24.4%
fma-def24.4%
fma-def24.4%
fma-def24.4%
fma-def24.4%
Simplified24.4%
Taylor expanded in x around inf 86.9%
Final simplification77.3%
(FPCore (x y z)
:precision binary64
(if (<= x -0.072)
(* (+ x -2.0) (- 4.16438922228 (/ 101.7851458539211 x)))
(if (<= x -2.6e-135)
(* (* x y) -0.0424927283095952)
(if (<= x -5.2e-191)
(* z -0.0424927283095952)
(if (<= x -3.1e-212)
(* x (* y -0.0424927283095952))
(if (<= x 3.5)
(/ (+ x -2.0) (/ 47.066876606 z))
(/ (+ x -2.0) 0.24013125253755718)))))))
double code(double x, double y, double z) {
double tmp;
if (x <= -0.072) {
tmp = (x + -2.0) * (4.16438922228 - (101.7851458539211 / x));
} else if (x <= -2.6e-135) {
tmp = (x * y) * -0.0424927283095952;
} else if (x <= -5.2e-191) {
tmp = z * -0.0424927283095952;
} else if (x <= -3.1e-212) {
tmp = x * (y * -0.0424927283095952);
} else if (x <= 3.5) {
tmp = (x + -2.0) / (47.066876606 / z);
} else {
tmp = (x + -2.0) / 0.24013125253755718;
}
return tmp;
}
real(8) function code(x, y, z)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8) :: tmp
if (x <= (-0.072d0)) then
tmp = (x + (-2.0d0)) * (4.16438922228d0 - (101.7851458539211d0 / x))
else if (x <= (-2.6d-135)) then
tmp = (x * y) * (-0.0424927283095952d0)
else if (x <= (-5.2d-191)) then
tmp = z * (-0.0424927283095952d0)
else if (x <= (-3.1d-212)) then
tmp = x * (y * (-0.0424927283095952d0))
else if (x <= 3.5d0) then
tmp = (x + (-2.0d0)) / (47.066876606d0 / z)
else
tmp = (x + (-2.0d0)) / 0.24013125253755718d0
end if
code = tmp
end function
public static double code(double x, double y, double z) {
double tmp;
if (x <= -0.072) {
tmp = (x + -2.0) * (4.16438922228 - (101.7851458539211 / x));
} else if (x <= -2.6e-135) {
tmp = (x * y) * -0.0424927283095952;
} else if (x <= -5.2e-191) {
tmp = z * -0.0424927283095952;
} else if (x <= -3.1e-212) {
tmp = x * (y * -0.0424927283095952);
} else if (x <= 3.5) {
tmp = (x + -2.0) / (47.066876606 / z);
} else {
tmp = (x + -2.0) / 0.24013125253755718;
}
return tmp;
}
def code(x, y, z): tmp = 0 if x <= -0.072: tmp = (x + -2.0) * (4.16438922228 - (101.7851458539211 / x)) elif x <= -2.6e-135: tmp = (x * y) * -0.0424927283095952 elif x <= -5.2e-191: tmp = z * -0.0424927283095952 elif x <= -3.1e-212: tmp = x * (y * -0.0424927283095952) elif x <= 3.5: tmp = (x + -2.0) / (47.066876606 / z) else: tmp = (x + -2.0) / 0.24013125253755718 return tmp
function code(x, y, z) tmp = 0.0 if (x <= -0.072) tmp = Float64(Float64(x + -2.0) * Float64(4.16438922228 - Float64(101.7851458539211 / x))); elseif (x <= -2.6e-135) tmp = Float64(Float64(x * y) * -0.0424927283095952); elseif (x <= -5.2e-191) tmp = Float64(z * -0.0424927283095952); elseif (x <= -3.1e-212) tmp = Float64(x * Float64(y * -0.0424927283095952)); elseif (x <= 3.5) tmp = Float64(Float64(x + -2.0) / Float64(47.066876606 / z)); else tmp = Float64(Float64(x + -2.0) / 0.24013125253755718); end return tmp end
function tmp_2 = code(x, y, z) tmp = 0.0; if (x <= -0.072) tmp = (x + -2.0) * (4.16438922228 - (101.7851458539211 / x)); elseif (x <= -2.6e-135) tmp = (x * y) * -0.0424927283095952; elseif (x <= -5.2e-191) tmp = z * -0.0424927283095952; elseif (x <= -3.1e-212) tmp = x * (y * -0.0424927283095952); elseif (x <= 3.5) tmp = (x + -2.0) / (47.066876606 / z); else tmp = (x + -2.0) / 0.24013125253755718; end tmp_2 = tmp; end
code[x_, y_, z_] := If[LessEqual[x, -0.072], N[(N[(x + -2.0), $MachinePrecision] * N[(4.16438922228 - N[(101.7851458539211 / x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[x, -2.6e-135], N[(N[(x * y), $MachinePrecision] * -0.0424927283095952), $MachinePrecision], If[LessEqual[x, -5.2e-191], N[(z * -0.0424927283095952), $MachinePrecision], If[LessEqual[x, -3.1e-212], N[(x * N[(y * -0.0424927283095952), $MachinePrecision]), $MachinePrecision], If[LessEqual[x, 3.5], N[(N[(x + -2.0), $MachinePrecision] / N[(47.066876606 / z), $MachinePrecision]), $MachinePrecision], N[(N[(x + -2.0), $MachinePrecision] / 0.24013125253755718), $MachinePrecision]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq -0.072:\\
\;\;\;\;\left(x + -2\right) \cdot \left(4.16438922228 - \frac{101.7851458539211}{x}\right)\\
\mathbf{elif}\;x \leq -2.6 \cdot 10^{-135}:\\
\;\;\;\;\left(x \cdot y\right) \cdot -0.0424927283095952\\
\mathbf{elif}\;x \leq -5.2 \cdot 10^{-191}:\\
\;\;\;\;z \cdot -0.0424927283095952\\
\mathbf{elif}\;x \leq -3.1 \cdot 10^{-212}:\\
\;\;\;\;x \cdot \left(y \cdot -0.0424927283095952\right)\\
\mathbf{elif}\;x \leq 3.5:\\
\;\;\;\;\frac{x + -2}{\frac{47.066876606}{z}}\\
\mathbf{else}:\\
\;\;\;\;\frac{x + -2}{0.24013125253755718}\\
\end{array}
\end{array}
if x < -0.0719999999999999946Initial program 16.7%
associate-*r/21.0%
sub-neg21.0%
metadata-eval21.0%
*-commutative21.0%
fma-def21.0%
*-commutative21.0%
fma-def21.0%
*-commutative21.0%
fma-def21.0%
fma-def21.0%
*-commutative21.0%
Simplified20.9%
Taylor expanded in x around inf 89.7%
associate-*r/89.7%
metadata-eval89.7%
Simplified89.7%
if -0.0719999999999999946 < x < -2.60000000000000004e-135Initial program 99.5%
*-commutative99.5%
associate-*r/99.3%
*-commutative99.3%
fma-def99.3%
*-commutative99.3%
fma-def99.3%
*-commutative99.3%
fma-def99.3%
fma-def99.3%
Simplified99.3%
Taylor expanded in z around 0 65.4%
Taylor expanded in x around 0 53.3%
if -2.60000000000000004e-135 < x < -5.19999999999999972e-191Initial program 99.8%
associate-*r/99.8%
sub-neg99.8%
metadata-eval99.8%
*-commutative99.8%
fma-def99.8%
*-commutative99.8%
fma-def99.8%
*-commutative99.8%
fma-def99.8%
fma-def99.8%
*-commutative99.8%
Simplified99.8%
Taylor expanded in x around 0 87.9%
*-commutative87.9%
Simplified87.9%
if -5.19999999999999972e-191 < x < -3.10000000000000006e-212Initial program 99.5%
*-commutative99.5%
associate-*r/99.7%
*-commutative99.7%
fma-def99.7%
*-commutative99.7%
fma-def99.7%
*-commutative99.7%
fma-def99.7%
fma-def99.7%
Simplified99.7%
Taylor expanded in z around 0 70.4%
Taylor expanded in x around 0 70.7%
associate-*r*70.9%
*-commutative70.9%
Simplified70.9%
if -3.10000000000000006e-212 < x < 3.5Initial program 99.6%
associate-/l*99.5%
sub-neg99.5%
metadata-eval99.5%
fma-def99.5%
fma-def99.5%
fma-def99.5%
fma-def99.5%
fma-def99.5%
fma-def99.5%
fma-def99.5%
Simplified99.5%
Taylor expanded in x around 0 64.2%
if 3.5 < x Initial program 19.2%
associate-/l*24.5%
sub-neg24.5%
metadata-eval24.5%
fma-def24.5%
fma-def24.5%
fma-def24.4%
fma-def24.4%
fma-def24.4%
fma-def24.4%
fma-def24.4%
Simplified24.4%
Taylor expanded in x around inf 86.9%
Final simplification77.3%
(FPCore (x y z) :precision binary64 (if (or (<= x -22000.0) (not (<= x 800.0))) (/ (+ x -2.0) (+ (/ 5.86923874282773 x) 0.24013125253755718)) (* (+ x -2.0) (+ (* z 0.0212463641547976) (* 0.0212463641547976 (* x y))))))
double code(double x, double y, double z) {
double tmp;
if ((x <= -22000.0) || !(x <= 800.0)) {
tmp = (x + -2.0) / ((5.86923874282773 / x) + 0.24013125253755718);
} else {
tmp = (x + -2.0) * ((z * 0.0212463641547976) + (0.0212463641547976 * (x * y)));
}
return tmp;
}
real(8) function code(x, y, z)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8) :: tmp
if ((x <= (-22000.0d0)) .or. (.not. (x <= 800.0d0))) then
tmp = (x + (-2.0d0)) / ((5.86923874282773d0 / x) + 0.24013125253755718d0)
else
tmp = (x + (-2.0d0)) * ((z * 0.0212463641547976d0) + (0.0212463641547976d0 * (x * y)))
end if
code = tmp
end function
public static double code(double x, double y, double z) {
double tmp;
if ((x <= -22000.0) || !(x <= 800.0)) {
tmp = (x + -2.0) / ((5.86923874282773 / x) + 0.24013125253755718);
} else {
tmp = (x + -2.0) * ((z * 0.0212463641547976) + (0.0212463641547976 * (x * y)));
}
return tmp;
}
def code(x, y, z): tmp = 0 if (x <= -22000.0) or not (x <= 800.0): tmp = (x + -2.0) / ((5.86923874282773 / x) + 0.24013125253755718) else: tmp = (x + -2.0) * ((z * 0.0212463641547976) + (0.0212463641547976 * (x * y))) return tmp
function code(x, y, z) tmp = 0.0 if ((x <= -22000.0) || !(x <= 800.0)) tmp = Float64(Float64(x + -2.0) / Float64(Float64(5.86923874282773 / x) + 0.24013125253755718)); else tmp = Float64(Float64(x + -2.0) * Float64(Float64(z * 0.0212463641547976) + Float64(0.0212463641547976 * Float64(x * y)))); end return tmp end
function tmp_2 = code(x, y, z) tmp = 0.0; if ((x <= -22000.0) || ~((x <= 800.0))) tmp = (x + -2.0) / ((5.86923874282773 / x) + 0.24013125253755718); else tmp = (x + -2.0) * ((z * 0.0212463641547976) + (0.0212463641547976 * (x * y))); end tmp_2 = tmp; end
code[x_, y_, z_] := If[Or[LessEqual[x, -22000.0], N[Not[LessEqual[x, 800.0]], $MachinePrecision]], N[(N[(x + -2.0), $MachinePrecision] / N[(N[(5.86923874282773 / x), $MachinePrecision] + 0.24013125253755718), $MachinePrecision]), $MachinePrecision], N[(N[(x + -2.0), $MachinePrecision] * N[(N[(z * 0.0212463641547976), $MachinePrecision] + N[(0.0212463641547976 * N[(x * y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq -22000 \lor \neg \left(x \leq 800\right):\\
\;\;\;\;\frac{x + -2}{\frac{5.86923874282773}{x} + 0.24013125253755718}\\
\mathbf{else}:\\
\;\;\;\;\left(x + -2\right) \cdot \left(z \cdot 0.0212463641547976 + 0.0212463641547976 \cdot \left(x \cdot y\right)\right)\\
\end{array}
\end{array}
if x < -22000 or 800 < x Initial program 16.2%
associate-/l*21.1%
sub-neg21.1%
metadata-eval21.1%
fma-def21.1%
fma-def21.1%
fma-def21.1%
fma-def21.1%
fma-def21.1%
fma-def21.1%
fma-def21.1%
Simplified21.1%
Taylor expanded in x around inf 90.5%
associate-*r/90.5%
metadata-eval90.5%
Simplified90.5%
if -22000 < x < 800Initial program 99.6%
associate-*r/99.6%
sub-neg99.6%
metadata-eval99.6%
*-commutative99.6%
fma-def99.6%
*-commutative99.6%
fma-def99.6%
*-commutative99.6%
fma-def99.6%
fma-def99.6%
*-commutative99.6%
Simplified99.6%
Taylor expanded in x around 0 91.5%
Taylor expanded in y around inf 91.4%
Final simplification90.9%
(FPCore (x y z)
:precision binary64
(let* ((t_0 (- (* x 4.16438922228) 110.1139242984811)))
(if (<= x -0.14)
t_0
(if (<= x -5.6e-136)
(* (* x y) -0.0424927283095952)
(if (<= x -5.2e-191)
(* z -0.0424927283095952)
(if (<= x -3.1e-212)
(* x (* y -0.0424927283095952))
(if (<= x 620.0) (* z -0.0424927283095952) t_0)))))))
double code(double x, double y, double z) {
double t_0 = (x * 4.16438922228) - 110.1139242984811;
double tmp;
if (x <= -0.14) {
tmp = t_0;
} else if (x <= -5.6e-136) {
tmp = (x * y) * -0.0424927283095952;
} else if (x <= -5.2e-191) {
tmp = z * -0.0424927283095952;
} else if (x <= -3.1e-212) {
tmp = x * (y * -0.0424927283095952);
} else if (x <= 620.0) {
tmp = z * -0.0424927283095952;
} else {
tmp = t_0;
}
return tmp;
}
real(8) function code(x, y, z)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8) :: t_0
real(8) :: tmp
t_0 = (x * 4.16438922228d0) - 110.1139242984811d0
if (x <= (-0.14d0)) then
tmp = t_0
else if (x <= (-5.6d-136)) then
tmp = (x * y) * (-0.0424927283095952d0)
else if (x <= (-5.2d-191)) then
tmp = z * (-0.0424927283095952d0)
else if (x <= (-3.1d-212)) then
tmp = x * (y * (-0.0424927283095952d0))
else if (x <= 620.0d0) then
tmp = z * (-0.0424927283095952d0)
else
tmp = t_0
end if
code = tmp
end function
public static double code(double x, double y, double z) {
double t_0 = (x * 4.16438922228) - 110.1139242984811;
double tmp;
if (x <= -0.14) {
tmp = t_0;
} else if (x <= -5.6e-136) {
tmp = (x * y) * -0.0424927283095952;
} else if (x <= -5.2e-191) {
tmp = z * -0.0424927283095952;
} else if (x <= -3.1e-212) {
tmp = x * (y * -0.0424927283095952);
} else if (x <= 620.0) {
tmp = z * -0.0424927283095952;
} else {
tmp = t_0;
}
return tmp;
}
def code(x, y, z): t_0 = (x * 4.16438922228) - 110.1139242984811 tmp = 0 if x <= -0.14: tmp = t_0 elif x <= -5.6e-136: tmp = (x * y) * -0.0424927283095952 elif x <= -5.2e-191: tmp = z * -0.0424927283095952 elif x <= -3.1e-212: tmp = x * (y * -0.0424927283095952) elif x <= 620.0: tmp = z * -0.0424927283095952 else: tmp = t_0 return tmp
function code(x, y, z) t_0 = Float64(Float64(x * 4.16438922228) - 110.1139242984811) tmp = 0.0 if (x <= -0.14) tmp = t_0; elseif (x <= -5.6e-136) tmp = Float64(Float64(x * y) * -0.0424927283095952); elseif (x <= -5.2e-191) tmp = Float64(z * -0.0424927283095952); elseif (x <= -3.1e-212) tmp = Float64(x * Float64(y * -0.0424927283095952)); elseif (x <= 620.0) tmp = Float64(z * -0.0424927283095952); else tmp = t_0; end return tmp end
function tmp_2 = code(x, y, z) t_0 = (x * 4.16438922228) - 110.1139242984811; tmp = 0.0; if (x <= -0.14) tmp = t_0; elseif (x <= -5.6e-136) tmp = (x * y) * -0.0424927283095952; elseif (x <= -5.2e-191) tmp = z * -0.0424927283095952; elseif (x <= -3.1e-212) tmp = x * (y * -0.0424927283095952); elseif (x <= 620.0) tmp = z * -0.0424927283095952; else tmp = t_0; end tmp_2 = tmp; end
code[x_, y_, z_] := Block[{t$95$0 = N[(N[(x * 4.16438922228), $MachinePrecision] - 110.1139242984811), $MachinePrecision]}, If[LessEqual[x, -0.14], t$95$0, If[LessEqual[x, -5.6e-136], N[(N[(x * y), $MachinePrecision] * -0.0424927283095952), $MachinePrecision], If[LessEqual[x, -5.2e-191], N[(z * -0.0424927283095952), $MachinePrecision], If[LessEqual[x, -3.1e-212], N[(x * N[(y * -0.0424927283095952), $MachinePrecision]), $MachinePrecision], If[LessEqual[x, 620.0], N[(z * -0.0424927283095952), $MachinePrecision], t$95$0]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := x \cdot 4.16438922228 - 110.1139242984811\\
\mathbf{if}\;x \leq -0.14:\\
\;\;\;\;t_0\\
\mathbf{elif}\;x \leq -5.6 \cdot 10^{-136}:\\
\;\;\;\;\left(x \cdot y\right) \cdot -0.0424927283095952\\
\mathbf{elif}\;x \leq -5.2 \cdot 10^{-191}:\\
\;\;\;\;z \cdot -0.0424927283095952\\
\mathbf{elif}\;x \leq -3.1 \cdot 10^{-212}:\\
\;\;\;\;x \cdot \left(y \cdot -0.0424927283095952\right)\\
\mathbf{elif}\;x \leq 620:\\
\;\;\;\;z \cdot -0.0424927283095952\\
\mathbf{else}:\\
\;\;\;\;t_0\\
\end{array}
\end{array}
if x < -0.14000000000000001 or 620 < x Initial program 16.8%
associate-*r/21.7%
sub-neg21.7%
metadata-eval21.7%
*-commutative21.7%
fma-def21.7%
*-commutative21.7%
fma-def21.7%
*-commutative21.7%
fma-def21.7%
fma-def21.7%
*-commutative21.7%
Simplified21.6%
Taylor expanded in x around inf 89.3%
if -0.14000000000000001 < x < -5.6000000000000002e-136Initial program 99.5%
*-commutative99.5%
associate-*r/99.3%
*-commutative99.3%
fma-def99.3%
*-commutative99.3%
fma-def99.3%
*-commutative99.3%
fma-def99.3%
fma-def99.3%
Simplified99.3%
Taylor expanded in z around 0 65.4%
Taylor expanded in x around 0 53.3%
if -5.6000000000000002e-136 < x < -5.19999999999999972e-191 or -3.10000000000000006e-212 < x < 620Initial program 99.6%
associate-*r/99.6%
sub-neg99.6%
metadata-eval99.6%
*-commutative99.6%
fma-def99.6%
*-commutative99.6%
fma-def99.6%
*-commutative99.6%
fma-def99.6%
fma-def99.6%
*-commutative99.6%
Simplified99.7%
Taylor expanded in x around 0 64.8%
*-commutative64.8%
Simplified64.8%
if -5.19999999999999972e-191 < x < -3.10000000000000006e-212Initial program 99.5%
*-commutative99.5%
associate-*r/99.7%
*-commutative99.7%
fma-def99.7%
*-commutative99.7%
fma-def99.7%
*-commutative99.7%
fma-def99.7%
fma-def99.7%
Simplified99.7%
Taylor expanded in z around 0 70.4%
Taylor expanded in x around 0 70.7%
associate-*r*70.9%
*-commutative70.9%
Simplified70.9%
Final simplification77.2%
(FPCore (x y z)
:precision binary64
(if (<= x -5.5)
(- (* x 4.16438922228) 110.1139242984811)
(if (<= x -2.6e-135)
(* (* x y) -0.0424927283095952)
(if (<= x -5.3e-191)
(* z -0.0424927283095952)
(if (<= x -7.2e-213)
(* x (* y -0.0424927283095952))
(if (<= x 620.0)
(* z -0.0424927283095952)
(/ (+ x -2.0) 0.24013125253755718)))))))
double code(double x, double y, double z) {
double tmp;
if (x <= -5.5) {
tmp = (x * 4.16438922228) - 110.1139242984811;
} else if (x <= -2.6e-135) {
tmp = (x * y) * -0.0424927283095952;
} else if (x <= -5.3e-191) {
tmp = z * -0.0424927283095952;
} else if (x <= -7.2e-213) {
tmp = x * (y * -0.0424927283095952);
} else if (x <= 620.0) {
tmp = z * -0.0424927283095952;
} else {
tmp = (x + -2.0) / 0.24013125253755718;
}
return tmp;
}
real(8) function code(x, y, z)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8) :: tmp
if (x <= (-5.5d0)) then
tmp = (x * 4.16438922228d0) - 110.1139242984811d0
else if (x <= (-2.6d-135)) then
tmp = (x * y) * (-0.0424927283095952d0)
else if (x <= (-5.3d-191)) then
tmp = z * (-0.0424927283095952d0)
else if (x <= (-7.2d-213)) then
tmp = x * (y * (-0.0424927283095952d0))
else if (x <= 620.0d0) then
tmp = z * (-0.0424927283095952d0)
else
tmp = (x + (-2.0d0)) / 0.24013125253755718d0
end if
code = tmp
end function
public static double code(double x, double y, double z) {
double tmp;
if (x <= -5.5) {
tmp = (x * 4.16438922228) - 110.1139242984811;
} else if (x <= -2.6e-135) {
tmp = (x * y) * -0.0424927283095952;
} else if (x <= -5.3e-191) {
tmp = z * -0.0424927283095952;
} else if (x <= -7.2e-213) {
tmp = x * (y * -0.0424927283095952);
} else if (x <= 620.0) {
tmp = z * -0.0424927283095952;
} else {
tmp = (x + -2.0) / 0.24013125253755718;
}
return tmp;
}
def code(x, y, z): tmp = 0 if x <= -5.5: tmp = (x * 4.16438922228) - 110.1139242984811 elif x <= -2.6e-135: tmp = (x * y) * -0.0424927283095952 elif x <= -5.3e-191: tmp = z * -0.0424927283095952 elif x <= -7.2e-213: tmp = x * (y * -0.0424927283095952) elif x <= 620.0: tmp = z * -0.0424927283095952 else: tmp = (x + -2.0) / 0.24013125253755718 return tmp
function code(x, y, z) tmp = 0.0 if (x <= -5.5) tmp = Float64(Float64(x * 4.16438922228) - 110.1139242984811); elseif (x <= -2.6e-135) tmp = Float64(Float64(x * y) * -0.0424927283095952); elseif (x <= -5.3e-191) tmp = Float64(z * -0.0424927283095952); elseif (x <= -7.2e-213) tmp = Float64(x * Float64(y * -0.0424927283095952)); elseif (x <= 620.0) tmp = Float64(z * -0.0424927283095952); else tmp = Float64(Float64(x + -2.0) / 0.24013125253755718); end return tmp end
function tmp_2 = code(x, y, z) tmp = 0.0; if (x <= -5.5) tmp = (x * 4.16438922228) - 110.1139242984811; elseif (x <= -2.6e-135) tmp = (x * y) * -0.0424927283095952; elseif (x <= -5.3e-191) tmp = z * -0.0424927283095952; elseif (x <= -7.2e-213) tmp = x * (y * -0.0424927283095952); elseif (x <= 620.0) tmp = z * -0.0424927283095952; else tmp = (x + -2.0) / 0.24013125253755718; end tmp_2 = tmp; end
code[x_, y_, z_] := If[LessEqual[x, -5.5], N[(N[(x * 4.16438922228), $MachinePrecision] - 110.1139242984811), $MachinePrecision], If[LessEqual[x, -2.6e-135], N[(N[(x * y), $MachinePrecision] * -0.0424927283095952), $MachinePrecision], If[LessEqual[x, -5.3e-191], N[(z * -0.0424927283095952), $MachinePrecision], If[LessEqual[x, -7.2e-213], N[(x * N[(y * -0.0424927283095952), $MachinePrecision]), $MachinePrecision], If[LessEqual[x, 620.0], N[(z * -0.0424927283095952), $MachinePrecision], N[(N[(x + -2.0), $MachinePrecision] / 0.24013125253755718), $MachinePrecision]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq -5.5:\\
\;\;\;\;x \cdot 4.16438922228 - 110.1139242984811\\
\mathbf{elif}\;x \leq -2.6 \cdot 10^{-135}:\\
\;\;\;\;\left(x \cdot y\right) \cdot -0.0424927283095952\\
\mathbf{elif}\;x \leq -5.3 \cdot 10^{-191}:\\
\;\;\;\;z \cdot -0.0424927283095952\\
\mathbf{elif}\;x \leq -7.2 \cdot 10^{-213}:\\
\;\;\;\;x \cdot \left(y \cdot -0.0424927283095952\right)\\
\mathbf{elif}\;x \leq 620:\\
\;\;\;\;z \cdot -0.0424927283095952\\
\mathbf{else}:\\
\;\;\;\;\frac{x + -2}{0.24013125253755718}\\
\end{array}
\end{array}
if x < -5.5Initial program 16.7%
associate-*r/21.0%
sub-neg21.0%
metadata-eval21.0%
*-commutative21.0%
fma-def21.0%
*-commutative21.0%
fma-def21.0%
*-commutative21.0%
fma-def21.0%
fma-def21.0%
*-commutative21.0%
Simplified20.9%
Taylor expanded in x around inf 89.7%
if -5.5 < x < -2.60000000000000004e-135Initial program 99.5%
*-commutative99.5%
associate-*r/99.3%
*-commutative99.3%
fma-def99.3%
*-commutative99.3%
fma-def99.3%
*-commutative99.3%
fma-def99.3%
fma-def99.3%
Simplified99.3%
Taylor expanded in z around 0 65.4%
Taylor expanded in x around 0 53.3%
if -2.60000000000000004e-135 < x < -5.29999999999999985e-191 or -7.2000000000000002e-213 < x < 620Initial program 99.6%
associate-*r/99.6%
sub-neg99.6%
metadata-eval99.6%
*-commutative99.6%
fma-def99.6%
*-commutative99.6%
fma-def99.6%
*-commutative99.6%
fma-def99.6%
fma-def99.6%
*-commutative99.6%
Simplified99.7%
Taylor expanded in x around 0 64.8%
*-commutative64.8%
Simplified64.8%
if -5.29999999999999985e-191 < x < -7.2000000000000002e-213Initial program 99.5%
*-commutative99.5%
associate-*r/99.7%
*-commutative99.7%
fma-def99.7%
*-commutative99.7%
fma-def99.7%
*-commutative99.7%
fma-def99.7%
fma-def99.7%
Simplified99.7%
Taylor expanded in z around 0 70.4%
Taylor expanded in x around 0 70.7%
associate-*r*70.9%
*-commutative70.9%
Simplified70.9%
if 620 < x Initial program 17.0%
associate-/l*22.4%
sub-neg22.4%
metadata-eval22.4%
fma-def22.4%
fma-def22.4%
fma-def22.3%
fma-def22.3%
fma-def22.3%
fma-def22.3%
fma-def22.3%
Simplified22.3%
Taylor expanded in x around inf 89.3%
Final simplification77.3%
(FPCore (x y z)
:precision binary64
(let* ((t_0 (* (* x y) -0.0424927283095952)))
(if (<= x -5.5)
(* x 4.16438922228)
(if (<= x -2.45e-135)
t_0
(if (<= x -5.2e-191)
(* z -0.0424927283095952)
(if (<= x -3.1e-212)
t_0
(if (<= x 620.0)
(* z -0.0424927283095952)
(* x 4.16438922228))))))))
double code(double x, double y, double z) {
double t_0 = (x * y) * -0.0424927283095952;
double tmp;
if (x <= -5.5) {
tmp = x * 4.16438922228;
} else if (x <= -2.45e-135) {
tmp = t_0;
} else if (x <= -5.2e-191) {
tmp = z * -0.0424927283095952;
} else if (x <= -3.1e-212) {
tmp = t_0;
} else if (x <= 620.0) {
tmp = z * -0.0424927283095952;
} else {
tmp = x * 4.16438922228;
}
return tmp;
}
real(8) function code(x, y, z)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8) :: t_0
real(8) :: tmp
t_0 = (x * y) * (-0.0424927283095952d0)
if (x <= (-5.5d0)) then
tmp = x * 4.16438922228d0
else if (x <= (-2.45d-135)) then
tmp = t_0
else if (x <= (-5.2d-191)) then
tmp = z * (-0.0424927283095952d0)
else if (x <= (-3.1d-212)) then
tmp = t_0
else if (x <= 620.0d0) then
tmp = z * (-0.0424927283095952d0)
else
tmp = x * 4.16438922228d0
end if
code = tmp
end function
public static double code(double x, double y, double z) {
double t_0 = (x * y) * -0.0424927283095952;
double tmp;
if (x <= -5.5) {
tmp = x * 4.16438922228;
} else if (x <= -2.45e-135) {
tmp = t_0;
} else if (x <= -5.2e-191) {
tmp = z * -0.0424927283095952;
} else if (x <= -3.1e-212) {
tmp = t_0;
} else if (x <= 620.0) {
tmp = z * -0.0424927283095952;
} else {
tmp = x * 4.16438922228;
}
return tmp;
}
def code(x, y, z): t_0 = (x * y) * -0.0424927283095952 tmp = 0 if x <= -5.5: tmp = x * 4.16438922228 elif x <= -2.45e-135: tmp = t_0 elif x <= -5.2e-191: tmp = z * -0.0424927283095952 elif x <= -3.1e-212: tmp = t_0 elif x <= 620.0: tmp = z * -0.0424927283095952 else: tmp = x * 4.16438922228 return tmp
function code(x, y, z) t_0 = Float64(Float64(x * y) * -0.0424927283095952) tmp = 0.0 if (x <= -5.5) tmp = Float64(x * 4.16438922228); elseif (x <= -2.45e-135) tmp = t_0; elseif (x <= -5.2e-191) tmp = Float64(z * -0.0424927283095952); elseif (x <= -3.1e-212) tmp = t_0; elseif (x <= 620.0) tmp = Float64(z * -0.0424927283095952); else tmp = Float64(x * 4.16438922228); end return tmp end
function tmp_2 = code(x, y, z) t_0 = (x * y) * -0.0424927283095952; tmp = 0.0; if (x <= -5.5) tmp = x * 4.16438922228; elseif (x <= -2.45e-135) tmp = t_0; elseif (x <= -5.2e-191) tmp = z * -0.0424927283095952; elseif (x <= -3.1e-212) tmp = t_0; elseif (x <= 620.0) tmp = z * -0.0424927283095952; else tmp = x * 4.16438922228; end tmp_2 = tmp; end
code[x_, y_, z_] := Block[{t$95$0 = N[(N[(x * y), $MachinePrecision] * -0.0424927283095952), $MachinePrecision]}, If[LessEqual[x, -5.5], N[(x * 4.16438922228), $MachinePrecision], If[LessEqual[x, -2.45e-135], t$95$0, If[LessEqual[x, -5.2e-191], N[(z * -0.0424927283095952), $MachinePrecision], If[LessEqual[x, -3.1e-212], t$95$0, If[LessEqual[x, 620.0], N[(z * -0.0424927283095952), $MachinePrecision], N[(x * 4.16438922228), $MachinePrecision]]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \left(x \cdot y\right) \cdot -0.0424927283095952\\
\mathbf{if}\;x \leq -5.5:\\
\;\;\;\;x \cdot 4.16438922228\\
\mathbf{elif}\;x \leq -2.45 \cdot 10^{-135}:\\
\;\;\;\;t_0\\
\mathbf{elif}\;x \leq -5.2 \cdot 10^{-191}:\\
\;\;\;\;z \cdot -0.0424927283095952\\
\mathbf{elif}\;x \leq -3.1 \cdot 10^{-212}:\\
\;\;\;\;t_0\\
\mathbf{elif}\;x \leq 620:\\
\;\;\;\;z \cdot -0.0424927283095952\\
\mathbf{else}:\\
\;\;\;\;x \cdot 4.16438922228\\
\end{array}
\end{array}
if x < -5.5 or 620 < x Initial program 16.8%
associate-*r/21.7%
sub-neg21.7%
metadata-eval21.7%
*-commutative21.7%
fma-def21.7%
*-commutative21.7%
fma-def21.7%
*-commutative21.7%
fma-def21.7%
fma-def21.7%
*-commutative21.7%
Simplified21.6%
Taylor expanded in x around inf 88.9%
*-commutative88.9%
Simplified88.9%
if -5.5 < x < -2.4500000000000001e-135 or -5.19999999999999972e-191 < x < -3.10000000000000006e-212Initial program 99.5%
*-commutative99.5%
associate-*r/99.4%
*-commutative99.4%
fma-def99.4%
*-commutative99.4%
fma-def99.4%
*-commutative99.4%
fma-def99.4%
fma-def99.4%
Simplified99.4%
Taylor expanded in z around 0 66.9%
Taylor expanded in x around 0 58.5%
if -2.4500000000000001e-135 < x < -5.19999999999999972e-191 or -3.10000000000000006e-212 < x < 620Initial program 99.6%
associate-*r/99.6%
sub-neg99.6%
metadata-eval99.6%
*-commutative99.6%
fma-def99.6%
*-commutative99.6%
fma-def99.6%
*-commutative99.6%
fma-def99.6%
fma-def99.6%
*-commutative99.6%
Simplified99.7%
Taylor expanded in x around 0 64.8%
*-commutative64.8%
Simplified64.8%
Final simplification77.0%
(FPCore (x y z)
:precision binary64
(if (<= x -0.2)
(* x 4.16438922228)
(if (<= x -2.6e-135)
(* (* x y) -0.0424927283095952)
(if (<= x -5.2e-191)
(* z -0.0424927283095952)
(if (<= x -3.1e-212)
(* x (* y -0.0424927283095952))
(if (<= x 720.0) (* z -0.0424927283095952) (* x 4.16438922228)))))))
double code(double x, double y, double z) {
double tmp;
if (x <= -0.2) {
tmp = x * 4.16438922228;
} else if (x <= -2.6e-135) {
tmp = (x * y) * -0.0424927283095952;
} else if (x <= -5.2e-191) {
tmp = z * -0.0424927283095952;
} else if (x <= -3.1e-212) {
tmp = x * (y * -0.0424927283095952);
} else if (x <= 720.0) {
tmp = z * -0.0424927283095952;
} else {
tmp = x * 4.16438922228;
}
return tmp;
}
real(8) function code(x, y, z)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8) :: tmp
if (x <= (-0.2d0)) then
tmp = x * 4.16438922228d0
else if (x <= (-2.6d-135)) then
tmp = (x * y) * (-0.0424927283095952d0)
else if (x <= (-5.2d-191)) then
tmp = z * (-0.0424927283095952d0)
else if (x <= (-3.1d-212)) then
tmp = x * (y * (-0.0424927283095952d0))
else if (x <= 720.0d0) then
tmp = z * (-0.0424927283095952d0)
else
tmp = x * 4.16438922228d0
end if
code = tmp
end function
public static double code(double x, double y, double z) {
double tmp;
if (x <= -0.2) {
tmp = x * 4.16438922228;
} else if (x <= -2.6e-135) {
tmp = (x * y) * -0.0424927283095952;
} else if (x <= -5.2e-191) {
tmp = z * -0.0424927283095952;
} else if (x <= -3.1e-212) {
tmp = x * (y * -0.0424927283095952);
} else if (x <= 720.0) {
tmp = z * -0.0424927283095952;
} else {
tmp = x * 4.16438922228;
}
return tmp;
}
def code(x, y, z): tmp = 0 if x <= -0.2: tmp = x * 4.16438922228 elif x <= -2.6e-135: tmp = (x * y) * -0.0424927283095952 elif x <= -5.2e-191: tmp = z * -0.0424927283095952 elif x <= -3.1e-212: tmp = x * (y * -0.0424927283095952) elif x <= 720.0: tmp = z * -0.0424927283095952 else: tmp = x * 4.16438922228 return tmp
function code(x, y, z) tmp = 0.0 if (x <= -0.2) tmp = Float64(x * 4.16438922228); elseif (x <= -2.6e-135) tmp = Float64(Float64(x * y) * -0.0424927283095952); elseif (x <= -5.2e-191) tmp = Float64(z * -0.0424927283095952); elseif (x <= -3.1e-212) tmp = Float64(x * Float64(y * -0.0424927283095952)); elseif (x <= 720.0) tmp = Float64(z * -0.0424927283095952); else tmp = Float64(x * 4.16438922228); end return tmp end
function tmp_2 = code(x, y, z) tmp = 0.0; if (x <= -0.2) tmp = x * 4.16438922228; elseif (x <= -2.6e-135) tmp = (x * y) * -0.0424927283095952; elseif (x <= -5.2e-191) tmp = z * -0.0424927283095952; elseif (x <= -3.1e-212) tmp = x * (y * -0.0424927283095952); elseif (x <= 720.0) tmp = z * -0.0424927283095952; else tmp = x * 4.16438922228; end tmp_2 = tmp; end
code[x_, y_, z_] := If[LessEqual[x, -0.2], N[(x * 4.16438922228), $MachinePrecision], If[LessEqual[x, -2.6e-135], N[(N[(x * y), $MachinePrecision] * -0.0424927283095952), $MachinePrecision], If[LessEqual[x, -5.2e-191], N[(z * -0.0424927283095952), $MachinePrecision], If[LessEqual[x, -3.1e-212], N[(x * N[(y * -0.0424927283095952), $MachinePrecision]), $MachinePrecision], If[LessEqual[x, 720.0], N[(z * -0.0424927283095952), $MachinePrecision], N[(x * 4.16438922228), $MachinePrecision]]]]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq -0.2:\\
\;\;\;\;x \cdot 4.16438922228\\
\mathbf{elif}\;x \leq -2.6 \cdot 10^{-135}:\\
\;\;\;\;\left(x \cdot y\right) \cdot -0.0424927283095952\\
\mathbf{elif}\;x \leq -5.2 \cdot 10^{-191}:\\
\;\;\;\;z \cdot -0.0424927283095952\\
\mathbf{elif}\;x \leq -3.1 \cdot 10^{-212}:\\
\;\;\;\;x \cdot \left(y \cdot -0.0424927283095952\right)\\
\mathbf{elif}\;x \leq 720:\\
\;\;\;\;z \cdot -0.0424927283095952\\
\mathbf{else}:\\
\;\;\;\;x \cdot 4.16438922228\\
\end{array}
\end{array}
if x < -0.20000000000000001 or 720 < x Initial program 16.8%
associate-*r/21.7%
sub-neg21.7%
metadata-eval21.7%
*-commutative21.7%
fma-def21.7%
*-commutative21.7%
fma-def21.7%
*-commutative21.7%
fma-def21.7%
fma-def21.7%
*-commutative21.7%
Simplified21.6%
Taylor expanded in x around inf 88.9%
*-commutative88.9%
Simplified88.9%
if -0.20000000000000001 < x < -2.60000000000000004e-135Initial program 99.5%
*-commutative99.5%
associate-*r/99.3%
*-commutative99.3%
fma-def99.3%
*-commutative99.3%
fma-def99.3%
*-commutative99.3%
fma-def99.3%
fma-def99.3%
Simplified99.3%
Taylor expanded in z around 0 65.4%
Taylor expanded in x around 0 53.3%
if -2.60000000000000004e-135 < x < -5.19999999999999972e-191 or -3.10000000000000006e-212 < x < 720Initial program 99.6%
associate-*r/99.6%
sub-neg99.6%
metadata-eval99.6%
*-commutative99.6%
fma-def99.6%
*-commutative99.6%
fma-def99.6%
*-commutative99.6%
fma-def99.6%
fma-def99.6%
*-commutative99.6%
Simplified99.7%
Taylor expanded in x around 0 64.8%
*-commutative64.8%
Simplified64.8%
if -5.19999999999999972e-191 < x < -3.10000000000000006e-212Initial program 99.5%
*-commutative99.5%
associate-*r/99.7%
*-commutative99.7%
fma-def99.7%
*-commutative99.7%
fma-def99.7%
*-commutative99.7%
fma-def99.7%
fma-def99.7%
Simplified99.7%
Taylor expanded in z around 0 70.4%
Taylor expanded in x around 0 70.7%
associate-*r*70.9%
*-commutative70.9%
Simplified70.9%
Final simplification77.0%
(FPCore (x y z) :precision binary64 (if (<= x -5500.0) (* x 4.16438922228) (if (<= x 620.0) (* z -0.0424927283095952) (* x 4.16438922228))))
double code(double x, double y, double z) {
double tmp;
if (x <= -5500.0) {
tmp = x * 4.16438922228;
} else if (x <= 620.0) {
tmp = z * -0.0424927283095952;
} else {
tmp = x * 4.16438922228;
}
return tmp;
}
real(8) function code(x, y, z)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8) :: tmp
if (x <= (-5500.0d0)) then
tmp = x * 4.16438922228d0
else if (x <= 620.0d0) then
tmp = z * (-0.0424927283095952d0)
else
tmp = x * 4.16438922228d0
end if
code = tmp
end function
public static double code(double x, double y, double z) {
double tmp;
if (x <= -5500.0) {
tmp = x * 4.16438922228;
} else if (x <= 620.0) {
tmp = z * -0.0424927283095952;
} else {
tmp = x * 4.16438922228;
}
return tmp;
}
def code(x, y, z): tmp = 0 if x <= -5500.0: tmp = x * 4.16438922228 elif x <= 620.0: tmp = z * -0.0424927283095952 else: tmp = x * 4.16438922228 return tmp
function code(x, y, z) tmp = 0.0 if (x <= -5500.0) tmp = Float64(x * 4.16438922228); elseif (x <= 620.0) tmp = Float64(z * -0.0424927283095952); else tmp = Float64(x * 4.16438922228); end return tmp end
function tmp_2 = code(x, y, z) tmp = 0.0; if (x <= -5500.0) tmp = x * 4.16438922228; elseif (x <= 620.0) tmp = z * -0.0424927283095952; else tmp = x * 4.16438922228; end tmp_2 = tmp; end
code[x_, y_, z_] := If[LessEqual[x, -5500.0], N[(x * 4.16438922228), $MachinePrecision], If[LessEqual[x, 620.0], N[(z * -0.0424927283095952), $MachinePrecision], N[(x * 4.16438922228), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq -5500:\\
\;\;\;\;x \cdot 4.16438922228\\
\mathbf{elif}\;x \leq 620:\\
\;\;\;\;z \cdot -0.0424927283095952\\
\mathbf{else}:\\
\;\;\;\;x \cdot 4.16438922228\\
\end{array}
\end{array}
if x < -5500 or 620 < x Initial program 16.2%
associate-*r/21.1%
sub-neg21.1%
metadata-eval21.1%
*-commutative21.1%
fma-def21.1%
*-commutative21.1%
fma-def21.1%
*-commutative21.1%
fma-def21.1%
fma-def21.1%
*-commutative21.1%
Simplified21.1%
Taylor expanded in x around inf 89.5%
*-commutative89.5%
Simplified89.5%
if -5500 < x < 620Initial program 99.6%
associate-*r/99.6%
sub-neg99.6%
metadata-eval99.6%
*-commutative99.6%
fma-def99.6%
*-commutative99.6%
fma-def99.6%
*-commutative99.6%
fma-def99.6%
fma-def99.6%
*-commutative99.6%
Simplified99.6%
Taylor expanded in x around 0 55.5%
*-commutative55.5%
Simplified55.5%
Final simplification73.6%
(FPCore (x y z) :precision binary64 (* x -0.3407596943375357))
double code(double x, double y, double z) {
return x * -0.3407596943375357;
}
real(8) function code(x, y, z)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
code = x * (-0.3407596943375357d0)
end function
public static double code(double x, double y, double z) {
return x * -0.3407596943375357;
}
def code(x, y, z): return x * -0.3407596943375357
function code(x, y, z) return Float64(x * -0.3407596943375357) end
function tmp = code(x, y, z) tmp = x * -0.3407596943375357; end
code[x_, y_, z_] := N[(x * -0.3407596943375357), $MachinePrecision]
\begin{array}{l}
\\
x \cdot -0.3407596943375357
\end{array}
Initial program 55.3%
associate-/l*57.8%
sub-neg57.8%
metadata-eval57.8%
fma-def57.8%
fma-def57.8%
fma-def57.8%
fma-def57.8%
fma-def57.8%
fma-def57.8%
fma-def57.8%
Simplified57.8%
Taylor expanded in x around inf 49.6%
associate-*r/49.6%
metadata-eval49.6%
Simplified49.6%
Taylor expanded in x around 0 1.9%
*-commutative1.9%
Simplified1.9%
Final simplification1.9%
(FPCore (x y z) :precision binary64 (* x 4.16438922228))
double code(double x, double y, double z) {
return x * 4.16438922228;
}
real(8) function code(x, y, z)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
code = x * 4.16438922228d0
end function
public static double code(double x, double y, double z) {
return x * 4.16438922228;
}
def code(x, y, z): return x * 4.16438922228
function code(x, y, z) return Float64(x * 4.16438922228) end
function tmp = code(x, y, z) tmp = x * 4.16438922228; end
code[x_, y_, z_] := N[(x * 4.16438922228), $MachinePrecision]
\begin{array}{l}
\\
x \cdot 4.16438922228
\end{array}
Initial program 55.3%
associate-*r/57.9%
sub-neg57.9%
metadata-eval57.9%
*-commutative57.9%
fma-def57.9%
*-commutative57.9%
fma-def57.9%
*-commutative57.9%
fma-def57.9%
fma-def57.9%
*-commutative57.9%
Simplified57.9%
Taylor expanded in x around inf 49.2%
*-commutative49.2%
Simplified49.2%
Final simplification49.2%
(FPCore (x y z)
:precision binary64
(let* ((t_0 (- (+ (/ y (* x x)) (* 4.16438922228 x)) 110.1139242984811)))
(if (< x -3.326128725870005e+62)
t_0
(if (< x 9.429991714554673e+55)
(*
(/ (- x 2.0) 1.0)
(/
(+
(*
(+
(* (+ (* (+ (* x 4.16438922228) 78.6994924154) x) 137.519416416) x)
y)
x)
z)
(+
(*
(+
(+ (* 263.505074721 x) (+ (* 43.3400022514 (* x x)) (* x (* x x))))
313.399215894)
x)
47.066876606)))
t_0))))
double code(double x, double y, double z) {
double t_0 = ((y / (x * x)) + (4.16438922228 * x)) - 110.1139242984811;
double tmp;
if (x < -3.326128725870005e+62) {
tmp = t_0;
} else if (x < 9.429991714554673e+55) {
tmp = ((x - 2.0) / 1.0) * (((((((((x * 4.16438922228) + 78.6994924154) * x) + 137.519416416) * x) + y) * x) + z) / (((((263.505074721 * x) + ((43.3400022514 * (x * x)) + (x * (x * x)))) + 313.399215894) * x) + 47.066876606));
} else {
tmp = t_0;
}
return tmp;
}
real(8) function code(x, y, z)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8), intent (in) :: z
real(8) :: t_0
real(8) :: tmp
t_0 = ((y / (x * x)) + (4.16438922228d0 * x)) - 110.1139242984811d0
if (x < (-3.326128725870005d+62)) then
tmp = t_0
else if (x < 9.429991714554673d+55) then
tmp = ((x - 2.0d0) / 1.0d0) * (((((((((x * 4.16438922228d0) + 78.6994924154d0) * x) + 137.519416416d0) * x) + y) * x) + z) / (((((263.505074721d0 * x) + ((43.3400022514d0 * (x * x)) + (x * (x * x)))) + 313.399215894d0) * x) + 47.066876606d0))
else
tmp = t_0
end if
code = tmp
end function
public static double code(double x, double y, double z) {
double t_0 = ((y / (x * x)) + (4.16438922228 * x)) - 110.1139242984811;
double tmp;
if (x < -3.326128725870005e+62) {
tmp = t_0;
} else if (x < 9.429991714554673e+55) {
tmp = ((x - 2.0) / 1.0) * (((((((((x * 4.16438922228) + 78.6994924154) * x) + 137.519416416) * x) + y) * x) + z) / (((((263.505074721 * x) + ((43.3400022514 * (x * x)) + (x * (x * x)))) + 313.399215894) * x) + 47.066876606));
} else {
tmp = t_0;
}
return tmp;
}
def code(x, y, z): t_0 = ((y / (x * x)) + (4.16438922228 * x)) - 110.1139242984811 tmp = 0 if x < -3.326128725870005e+62: tmp = t_0 elif x < 9.429991714554673e+55: tmp = ((x - 2.0) / 1.0) * (((((((((x * 4.16438922228) + 78.6994924154) * x) + 137.519416416) * x) + y) * x) + z) / (((((263.505074721 * x) + ((43.3400022514 * (x * x)) + (x * (x * x)))) + 313.399215894) * x) + 47.066876606)) else: tmp = t_0 return tmp
function code(x, y, z) t_0 = Float64(Float64(Float64(y / Float64(x * x)) + Float64(4.16438922228 * x)) - 110.1139242984811) tmp = 0.0 if (x < -3.326128725870005e+62) tmp = t_0; elseif (x < 9.429991714554673e+55) tmp = Float64(Float64(Float64(x - 2.0) / 1.0) * Float64(Float64(Float64(Float64(Float64(Float64(Float64(Float64(Float64(x * 4.16438922228) + 78.6994924154) * x) + 137.519416416) * x) + y) * x) + z) / Float64(Float64(Float64(Float64(Float64(263.505074721 * x) + Float64(Float64(43.3400022514 * Float64(x * x)) + Float64(x * Float64(x * x)))) + 313.399215894) * x) + 47.066876606))); else tmp = t_0; end return tmp end
function tmp_2 = code(x, y, z) t_0 = ((y / (x * x)) + (4.16438922228 * x)) - 110.1139242984811; tmp = 0.0; if (x < -3.326128725870005e+62) tmp = t_0; elseif (x < 9.429991714554673e+55) tmp = ((x - 2.0) / 1.0) * (((((((((x * 4.16438922228) + 78.6994924154) * x) + 137.519416416) * x) + y) * x) + z) / (((((263.505074721 * x) + ((43.3400022514 * (x * x)) + (x * (x * x)))) + 313.399215894) * x) + 47.066876606)); else tmp = t_0; end tmp_2 = tmp; end
code[x_, y_, z_] := Block[{t$95$0 = N[(N[(N[(y / N[(x * x), $MachinePrecision]), $MachinePrecision] + N[(4.16438922228 * x), $MachinePrecision]), $MachinePrecision] - 110.1139242984811), $MachinePrecision]}, If[Less[x, -3.326128725870005e+62], t$95$0, If[Less[x, 9.429991714554673e+55], N[(N[(N[(x - 2.0), $MachinePrecision] / 1.0), $MachinePrecision] * N[(N[(N[(N[(N[(N[(N[(N[(N[(x * 4.16438922228), $MachinePrecision] + 78.6994924154), $MachinePrecision] * x), $MachinePrecision] + 137.519416416), $MachinePrecision] * x), $MachinePrecision] + y), $MachinePrecision] * x), $MachinePrecision] + z), $MachinePrecision] / N[(N[(N[(N[(N[(263.505074721 * x), $MachinePrecision] + N[(N[(43.3400022514 * N[(x * x), $MachinePrecision]), $MachinePrecision] + N[(x * N[(x * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + 313.399215894), $MachinePrecision] * x), $MachinePrecision] + 47.066876606), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$0]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \left(\frac{y}{x \cdot x} + 4.16438922228 \cdot x\right) - 110.1139242984811\\
\mathbf{if}\;x < -3.326128725870005 \cdot 10^{+62}:\\
\;\;\;\;t_0\\
\mathbf{elif}\;x < 9.429991714554673 \cdot 10^{+55}:\\
\;\;\;\;\frac{x - 2}{1} \cdot \frac{\left(\left(\left(x \cdot 4.16438922228 + 78.6994924154\right) \cdot x + 137.519416416\right) \cdot x + y\right) \cdot x + z}{\left(\left(263.505074721 \cdot x + \left(43.3400022514 \cdot \left(x \cdot x\right) + x \cdot \left(x \cdot x\right)\right)\right) + 313.399215894\right) \cdot x + 47.066876606}\\
\mathbf{else}:\\
\;\;\;\;t_0\\
\end{array}
\end{array}
herbie shell --seed 2023199
(FPCore (x y z)
:name "Numeric.SpecFunctions:logGamma from math-functions-0.1.5.2, C"
:precision binary64
:herbie-target
(if (< x -3.326128725870005e+62) (- (+ (/ y (* x x)) (* 4.16438922228 x)) 110.1139242984811) (if (< x 9.429991714554673e+55) (* (/ (- x 2.0) 1.0) (/ (+ (* (+ (* (+ (* (+ (* x 4.16438922228) 78.6994924154) x) 137.519416416) x) y) x) z) (+ (* (+ (+ (* 263.505074721 x) (+ (* 43.3400022514 (* x x)) (* x (* x x)))) 313.399215894) x) 47.066876606))) (- (+ (/ y (* x x)) (* 4.16438922228 x)) 110.1139242984811)))
(/ (* (- x 2.0) (+ (* (+ (* (+ (* (+ (* x 4.16438922228) 78.6994924154) x) 137.519416416) x) y) x) z)) (+ (* (+ (* (+ (* (+ x 43.3400022514) x) 263.505074721) x) 313.399215894) x) 47.066876606)))