
(FPCore (x y) :precision binary64 (* (sin x) (/ (sinh y) y)))
double code(double x, double y) {
return sin(x) * (sinh(y) / y);
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
code = sin(x) * (sinh(y) / y)
end function
public static double code(double x, double y) {
return Math.sin(x) * (Math.sinh(y) / y);
}
def code(x, y): return math.sin(x) * (math.sinh(y) / y)
function code(x, y) return Float64(sin(x) * Float64(sinh(y) / y)) end
function tmp = code(x, y) tmp = sin(x) * (sinh(y) / y); end
code[x_, y_] := N[(N[Sin[x], $MachinePrecision] * N[(N[Sinh[y], $MachinePrecision] / y), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\sin x \cdot \frac{\sinh y}{y}
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 21 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (x y) :precision binary64 (* (sin x) (/ (sinh y) y)))
double code(double x, double y) {
return sin(x) * (sinh(y) / y);
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
code = sin(x) * (sinh(y) / y)
end function
public static double code(double x, double y) {
return Math.sin(x) * (Math.sinh(y) / y);
}
def code(x, y): return math.sin(x) * (math.sinh(y) / y)
function code(x, y) return Float64(sin(x) * Float64(sinh(y) / y)) end
function tmp = code(x, y) tmp = sin(x) * (sinh(y) / y); end
code[x_, y_] := N[(N[Sin[x], $MachinePrecision] * N[(N[Sinh[y], $MachinePrecision] / y), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\sin x \cdot \frac{\sinh y}{y}
\end{array}
(FPCore (x y) :precision binary64 (/ (sin x) (/ y (sinh y))))
double code(double x, double y) {
return sin(x) / (y / sinh(y));
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
code = sin(x) / (y / sinh(y))
end function
public static double code(double x, double y) {
return Math.sin(x) / (y / Math.sinh(y));
}
def code(x, y): return math.sin(x) / (y / math.sinh(y))
function code(x, y) return Float64(sin(x) / Float64(y / sinh(y))) end
function tmp = code(x, y) tmp = sin(x) / (y / sinh(y)); end
code[x_, y_] := N[(N[Sin[x], $MachinePrecision] / N[(y / N[Sinh[y], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{\sin x}{\frac{y}{\sinh y}}
\end{array}
Initial program 100.0%
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
sin-lowering-sin.f64N/A
/-lowering-/.f64N/A
sinh-lowering-sinh.f64100.0%
Applied egg-rr100.0%
(FPCore (x y) :precision binary64 (* (sin x) (/ (sinh y) y)))
double code(double x, double y) {
return sin(x) * (sinh(y) / y);
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
code = sin(x) * (sinh(y) / y)
end function
public static double code(double x, double y) {
return Math.sin(x) * (Math.sinh(y) / y);
}
def code(x, y): return math.sin(x) * (math.sinh(y) / y)
function code(x, y) return Float64(sin(x) * Float64(sinh(y) / y)) end
function tmp = code(x, y) tmp = sin(x) * (sinh(y) / y); end
code[x_, y_] := N[(N[Sin[x], $MachinePrecision] * N[(N[Sinh[y], $MachinePrecision] / y), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\sin x \cdot \frac{\sinh y}{y}
\end{array}
Initial program 100.0%
(FPCore (x y)
:precision binary64
(if (<= y 5.9e-19)
(sin x)
(if (<= y 3.8e+77)
(* (* x (/ 2.0 y)) (/ (sinh y) 2.0))
(* (sin x) (* y (* y (* (* y y) 0.008333333333333333)))))))
double code(double x, double y) {
double tmp;
if (y <= 5.9e-19) {
tmp = sin(x);
} else if (y <= 3.8e+77) {
tmp = (x * (2.0 / y)) * (sinh(y) / 2.0);
} else {
tmp = sin(x) * (y * (y * ((y * y) * 0.008333333333333333)));
}
return tmp;
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8) :: tmp
if (y <= 5.9d-19) then
tmp = sin(x)
else if (y <= 3.8d+77) then
tmp = (x * (2.0d0 / y)) * (sinh(y) / 2.0d0)
else
tmp = sin(x) * (y * (y * ((y * y) * 0.008333333333333333d0)))
end if
code = tmp
end function
public static double code(double x, double y) {
double tmp;
if (y <= 5.9e-19) {
tmp = Math.sin(x);
} else if (y <= 3.8e+77) {
tmp = (x * (2.0 / y)) * (Math.sinh(y) / 2.0);
} else {
tmp = Math.sin(x) * (y * (y * ((y * y) * 0.008333333333333333)));
}
return tmp;
}
def code(x, y): tmp = 0 if y <= 5.9e-19: tmp = math.sin(x) elif y <= 3.8e+77: tmp = (x * (2.0 / y)) * (math.sinh(y) / 2.0) else: tmp = math.sin(x) * (y * (y * ((y * y) * 0.008333333333333333))) return tmp
function code(x, y) tmp = 0.0 if (y <= 5.9e-19) tmp = sin(x); elseif (y <= 3.8e+77) tmp = Float64(Float64(x * Float64(2.0 / y)) * Float64(sinh(y) / 2.0)); else tmp = Float64(sin(x) * Float64(y * Float64(y * Float64(Float64(y * y) * 0.008333333333333333)))); end return tmp end
function tmp_2 = code(x, y) tmp = 0.0; if (y <= 5.9e-19) tmp = sin(x); elseif (y <= 3.8e+77) tmp = (x * (2.0 / y)) * (sinh(y) / 2.0); else tmp = sin(x) * (y * (y * ((y * y) * 0.008333333333333333))); end tmp_2 = tmp; end
code[x_, y_] := If[LessEqual[y, 5.9e-19], N[Sin[x], $MachinePrecision], If[LessEqual[y, 3.8e+77], N[(N[(x * N[(2.0 / y), $MachinePrecision]), $MachinePrecision] * N[(N[Sinh[y], $MachinePrecision] / 2.0), $MachinePrecision]), $MachinePrecision], N[(N[Sin[x], $MachinePrecision] * N[(y * N[(y * N[(N[(y * y), $MachinePrecision] * 0.008333333333333333), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq 5.9 \cdot 10^{-19}:\\
\;\;\;\;\sin x\\
\mathbf{elif}\;y \leq 3.8 \cdot 10^{+77}:\\
\;\;\;\;\left(x \cdot \frac{2}{y}\right) \cdot \frac{\sinh y}{2}\\
\mathbf{else}:\\
\;\;\;\;\sin x \cdot \left(y \cdot \left(y \cdot \left(\left(y \cdot y\right) \cdot 0.008333333333333333\right)\right)\right)\\
\end{array}
\end{array}
if y < 5.90000000000000038e-19Initial program 100.0%
Taylor expanded in y around 0
sin-lowering-sin.f6470.1%
Simplified70.1%
if 5.90000000000000038e-19 < y < 3.8000000000000001e77Initial program 99.8%
sinh-defN/A
associate-/l/N/A
sinh-undefN/A
times-fracN/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
sinh-lowering-sinh.f64100.0%
Applied egg-rr100.0%
Taylor expanded in x around 0
Simplified92.3%
if 3.8000000000000001e77 < y Initial program 100.0%
Taylor expanded in y around 0
distribute-rgt-inN/A
*-rgt-identityN/A
distribute-rgt-inN/A
*-commutativeN/A
associate-*r*N/A
distribute-rgt-outN/A
+-commutativeN/A
associate-*l*N/A
*-commutativeN/A
distribute-lft-inN/A
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
+-lowering-+.f64N/A
Simplified100.0%
Taylor expanded in y around inf
associate-*r*N/A
*-commutativeN/A
metadata-evalN/A
pow-sqrN/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
associate-*r*N/A
unpow2N/A
unpow3N/A
*-lowering-*.f64N/A
unpow3N/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64100.0%
Simplified100.0%
(FPCore (x y)
:precision binary64
(if (<= y 5.9e-19)
(sin x)
(if (<= y 4.5e+146)
(/ x (/ y (sinh y)))
(* (sin x) (+ 1.0 (* y (* y 0.16666666666666666)))))))
double code(double x, double y) {
double tmp;
if (y <= 5.9e-19) {
tmp = sin(x);
} else if (y <= 4.5e+146) {
tmp = x / (y / sinh(y));
} else {
tmp = sin(x) * (1.0 + (y * (y * 0.16666666666666666)));
}
return tmp;
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8) :: tmp
if (y <= 5.9d-19) then
tmp = sin(x)
else if (y <= 4.5d+146) then
tmp = x / (y / sinh(y))
else
tmp = sin(x) * (1.0d0 + (y * (y * 0.16666666666666666d0)))
end if
code = tmp
end function
public static double code(double x, double y) {
double tmp;
if (y <= 5.9e-19) {
tmp = Math.sin(x);
} else if (y <= 4.5e+146) {
tmp = x / (y / Math.sinh(y));
} else {
tmp = Math.sin(x) * (1.0 + (y * (y * 0.16666666666666666)));
}
return tmp;
}
def code(x, y): tmp = 0 if y <= 5.9e-19: tmp = math.sin(x) elif y <= 4.5e+146: tmp = x / (y / math.sinh(y)) else: tmp = math.sin(x) * (1.0 + (y * (y * 0.16666666666666666))) return tmp
function code(x, y) tmp = 0.0 if (y <= 5.9e-19) tmp = sin(x); elseif (y <= 4.5e+146) tmp = Float64(x / Float64(y / sinh(y))); else tmp = Float64(sin(x) * Float64(1.0 + Float64(y * Float64(y * 0.16666666666666666)))); end return tmp end
function tmp_2 = code(x, y) tmp = 0.0; if (y <= 5.9e-19) tmp = sin(x); elseif (y <= 4.5e+146) tmp = x / (y / sinh(y)); else tmp = sin(x) * (1.0 + (y * (y * 0.16666666666666666))); end tmp_2 = tmp; end
code[x_, y_] := If[LessEqual[y, 5.9e-19], N[Sin[x], $MachinePrecision], If[LessEqual[y, 4.5e+146], N[(x / N[(y / N[Sinh[y], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[Sin[x], $MachinePrecision] * N[(1.0 + N[(y * N[(y * 0.16666666666666666), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq 5.9 \cdot 10^{-19}:\\
\;\;\;\;\sin x\\
\mathbf{elif}\;y \leq 4.5 \cdot 10^{+146}:\\
\;\;\;\;\frac{x}{\frac{y}{\sinh y}}\\
\mathbf{else}:\\
\;\;\;\;\sin x \cdot \left(1 + y \cdot \left(y \cdot 0.16666666666666666\right)\right)\\
\end{array}
\end{array}
if y < 5.90000000000000038e-19Initial program 100.0%
Taylor expanded in y around 0
sin-lowering-sin.f6470.1%
Simplified70.1%
if 5.90000000000000038e-19 < y < 4.50000000000000026e146Initial program 99.9%
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
sin-lowering-sin.f64N/A
/-lowering-/.f64N/A
sinh-lowering-sinh.f6499.9%
Applied egg-rr99.9%
Taylor expanded in x around 0
Simplified87.0%
if 4.50000000000000026e146 < y Initial program 100.0%
Taylor expanded in y around 0
*-lft-identityN/A
associate-*r*N/A
distribute-rgt-inN/A
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f6497.8%
Simplified97.8%
(FPCore (x y)
:precision binary64
(if (<= y 5.9e-19)
(sin x)
(if (<= y 4.5e+146)
(/ x (/ y (sinh y)))
(* (sin x) (* (* y y) 0.16666666666666666)))))
double code(double x, double y) {
double tmp;
if (y <= 5.9e-19) {
tmp = sin(x);
} else if (y <= 4.5e+146) {
tmp = x / (y / sinh(y));
} else {
tmp = sin(x) * ((y * y) * 0.16666666666666666);
}
return tmp;
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8) :: tmp
if (y <= 5.9d-19) then
tmp = sin(x)
else if (y <= 4.5d+146) then
tmp = x / (y / sinh(y))
else
tmp = sin(x) * ((y * y) * 0.16666666666666666d0)
end if
code = tmp
end function
public static double code(double x, double y) {
double tmp;
if (y <= 5.9e-19) {
tmp = Math.sin(x);
} else if (y <= 4.5e+146) {
tmp = x / (y / Math.sinh(y));
} else {
tmp = Math.sin(x) * ((y * y) * 0.16666666666666666);
}
return tmp;
}
def code(x, y): tmp = 0 if y <= 5.9e-19: tmp = math.sin(x) elif y <= 4.5e+146: tmp = x / (y / math.sinh(y)) else: tmp = math.sin(x) * ((y * y) * 0.16666666666666666) return tmp
function code(x, y) tmp = 0.0 if (y <= 5.9e-19) tmp = sin(x); elseif (y <= 4.5e+146) tmp = Float64(x / Float64(y / sinh(y))); else tmp = Float64(sin(x) * Float64(Float64(y * y) * 0.16666666666666666)); end return tmp end
function tmp_2 = code(x, y) tmp = 0.0; if (y <= 5.9e-19) tmp = sin(x); elseif (y <= 4.5e+146) tmp = x / (y / sinh(y)); else tmp = sin(x) * ((y * y) * 0.16666666666666666); end tmp_2 = tmp; end
code[x_, y_] := If[LessEqual[y, 5.9e-19], N[Sin[x], $MachinePrecision], If[LessEqual[y, 4.5e+146], N[(x / N[(y / N[Sinh[y], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[Sin[x], $MachinePrecision] * N[(N[(y * y), $MachinePrecision] * 0.16666666666666666), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq 5.9 \cdot 10^{-19}:\\
\;\;\;\;\sin x\\
\mathbf{elif}\;y \leq 4.5 \cdot 10^{+146}:\\
\;\;\;\;\frac{x}{\frac{y}{\sinh y}}\\
\mathbf{else}:\\
\;\;\;\;\sin x \cdot \left(\left(y \cdot y\right) \cdot 0.16666666666666666\right)\\
\end{array}
\end{array}
if y < 5.90000000000000038e-19Initial program 100.0%
Taylor expanded in y around 0
sin-lowering-sin.f6470.1%
Simplified70.1%
if 5.90000000000000038e-19 < y < 4.50000000000000026e146Initial program 99.9%
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
sin-lowering-sin.f64N/A
/-lowering-/.f64N/A
sinh-lowering-sinh.f6499.9%
Applied egg-rr99.9%
Taylor expanded in x around 0
Simplified87.0%
if 4.50000000000000026e146 < y Initial program 100.0%
Taylor expanded in y around 0
*-lft-identityN/A
associate-*r*N/A
distribute-rgt-inN/A
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f6497.8%
Simplified97.8%
Taylor expanded in y around inf
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6497.8%
Simplified97.8%
Final simplification76.6%
(FPCore (x y)
:precision binary64
(if (<= y 5.9e-19)
(sin x)
(if (<= y 2e+254)
(/ x (/ y (sinh y)))
(*
(* y (* x (+ 0.008333333333333333 (* (* x x) -0.001388888888888889))))
(* y (* y y))))))
double code(double x, double y) {
double tmp;
if (y <= 5.9e-19) {
tmp = sin(x);
} else if (y <= 2e+254) {
tmp = x / (y / sinh(y));
} else {
tmp = (y * (x * (0.008333333333333333 + ((x * x) * -0.001388888888888889)))) * (y * (y * y));
}
return tmp;
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8) :: tmp
if (y <= 5.9d-19) then
tmp = sin(x)
else if (y <= 2d+254) then
tmp = x / (y / sinh(y))
else
tmp = (y * (x * (0.008333333333333333d0 + ((x * x) * (-0.001388888888888889d0))))) * (y * (y * y))
end if
code = tmp
end function
public static double code(double x, double y) {
double tmp;
if (y <= 5.9e-19) {
tmp = Math.sin(x);
} else if (y <= 2e+254) {
tmp = x / (y / Math.sinh(y));
} else {
tmp = (y * (x * (0.008333333333333333 + ((x * x) * -0.001388888888888889)))) * (y * (y * y));
}
return tmp;
}
def code(x, y): tmp = 0 if y <= 5.9e-19: tmp = math.sin(x) elif y <= 2e+254: tmp = x / (y / math.sinh(y)) else: tmp = (y * (x * (0.008333333333333333 + ((x * x) * -0.001388888888888889)))) * (y * (y * y)) return tmp
function code(x, y) tmp = 0.0 if (y <= 5.9e-19) tmp = sin(x); elseif (y <= 2e+254) tmp = Float64(x / Float64(y / sinh(y))); else tmp = Float64(Float64(y * Float64(x * Float64(0.008333333333333333 + Float64(Float64(x * x) * -0.001388888888888889)))) * Float64(y * Float64(y * y))); end return tmp end
function tmp_2 = code(x, y) tmp = 0.0; if (y <= 5.9e-19) tmp = sin(x); elseif (y <= 2e+254) tmp = x / (y / sinh(y)); else tmp = (y * (x * (0.008333333333333333 + ((x * x) * -0.001388888888888889)))) * (y * (y * y)); end tmp_2 = tmp; end
code[x_, y_] := If[LessEqual[y, 5.9e-19], N[Sin[x], $MachinePrecision], If[LessEqual[y, 2e+254], N[(x / N[(y / N[Sinh[y], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(y * N[(x * N[(0.008333333333333333 + N[(N[(x * x), $MachinePrecision] * -0.001388888888888889), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(y * N[(y * y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq 5.9 \cdot 10^{-19}:\\
\;\;\;\;\sin x\\
\mathbf{elif}\;y \leq 2 \cdot 10^{+254}:\\
\;\;\;\;\frac{x}{\frac{y}{\sinh y}}\\
\mathbf{else}:\\
\;\;\;\;\left(y \cdot \left(x \cdot \left(0.008333333333333333 + \left(x \cdot x\right) \cdot -0.001388888888888889\right)\right)\right) \cdot \left(y \cdot \left(y \cdot y\right)\right)\\
\end{array}
\end{array}
if y < 5.90000000000000038e-19Initial program 100.0%
Taylor expanded in y around 0
sin-lowering-sin.f6470.1%
Simplified70.1%
if 5.90000000000000038e-19 < y < 1.9999999999999999e254Initial program 100.0%
clear-numN/A
un-div-invN/A
/-lowering-/.f64N/A
sin-lowering-sin.f64N/A
/-lowering-/.f64N/A
sinh-lowering-sinh.f64100.0%
Applied egg-rr100.0%
Taylor expanded in x around 0
Simplified83.0%
if 1.9999999999999999e254 < y Initial program 100.0%
Taylor expanded in y around 0
distribute-rgt-inN/A
*-rgt-identityN/A
distribute-rgt-inN/A
*-commutativeN/A
associate-*r*N/A
distribute-rgt-outN/A
+-commutativeN/A
associate-*l*N/A
*-commutativeN/A
distribute-lft-inN/A
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
+-lowering-+.f64N/A
Simplified100.0%
Taylor expanded in y around inf
associate-*r*N/A
*-commutativeN/A
metadata-evalN/A
pow-sqrN/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
associate-*r*N/A
unpow2N/A
unpow3N/A
*-lowering-*.f64N/A
unpow3N/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64100.0%
Simplified100.0%
Taylor expanded in x around 0
distribute-lft-inN/A
associate-*r*N/A
associate-*r*N/A
associate-*r*N/A
*-commutativeN/A
distribute-rgt-outN/A
*-lowering-*.f64N/A
metadata-evalN/A
pow-sqrN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-commutativeN/A
distribute-lft-outN/A
+-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64100.0%
Simplified100.0%
*-commutativeN/A
associate-*l*N/A
cube-unmultN/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
cube-unmultN/A
*-lowering-*.f64N/A
*-lowering-*.f64100.0%
Applied egg-rr100.0%
Final simplification74.6%
(FPCore (x y)
:precision binary64
(if (<= y 5.9e-19)
(sin x)
(if (<= y 3.95e+254)
(* x (/ (sinh y) y))
(*
(* y (* x (+ 0.008333333333333333 (* (* x x) -0.001388888888888889))))
(* y (* y y))))))
double code(double x, double y) {
double tmp;
if (y <= 5.9e-19) {
tmp = sin(x);
} else if (y <= 3.95e+254) {
tmp = x * (sinh(y) / y);
} else {
tmp = (y * (x * (0.008333333333333333 + ((x * x) * -0.001388888888888889)))) * (y * (y * y));
}
return tmp;
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8) :: tmp
if (y <= 5.9d-19) then
tmp = sin(x)
else if (y <= 3.95d+254) then
tmp = x * (sinh(y) / y)
else
tmp = (y * (x * (0.008333333333333333d0 + ((x * x) * (-0.001388888888888889d0))))) * (y * (y * y))
end if
code = tmp
end function
public static double code(double x, double y) {
double tmp;
if (y <= 5.9e-19) {
tmp = Math.sin(x);
} else if (y <= 3.95e+254) {
tmp = x * (Math.sinh(y) / y);
} else {
tmp = (y * (x * (0.008333333333333333 + ((x * x) * -0.001388888888888889)))) * (y * (y * y));
}
return tmp;
}
def code(x, y): tmp = 0 if y <= 5.9e-19: tmp = math.sin(x) elif y <= 3.95e+254: tmp = x * (math.sinh(y) / y) else: tmp = (y * (x * (0.008333333333333333 + ((x * x) * -0.001388888888888889)))) * (y * (y * y)) return tmp
function code(x, y) tmp = 0.0 if (y <= 5.9e-19) tmp = sin(x); elseif (y <= 3.95e+254) tmp = Float64(x * Float64(sinh(y) / y)); else tmp = Float64(Float64(y * Float64(x * Float64(0.008333333333333333 + Float64(Float64(x * x) * -0.001388888888888889)))) * Float64(y * Float64(y * y))); end return tmp end
function tmp_2 = code(x, y) tmp = 0.0; if (y <= 5.9e-19) tmp = sin(x); elseif (y <= 3.95e+254) tmp = x * (sinh(y) / y); else tmp = (y * (x * (0.008333333333333333 + ((x * x) * -0.001388888888888889)))) * (y * (y * y)); end tmp_2 = tmp; end
code[x_, y_] := If[LessEqual[y, 5.9e-19], N[Sin[x], $MachinePrecision], If[LessEqual[y, 3.95e+254], N[(x * N[(N[Sinh[y], $MachinePrecision] / y), $MachinePrecision]), $MachinePrecision], N[(N[(y * N[(x * N[(0.008333333333333333 + N[(N[(x * x), $MachinePrecision] * -0.001388888888888889), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(y * N[(y * y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq 5.9 \cdot 10^{-19}:\\
\;\;\;\;\sin x\\
\mathbf{elif}\;y \leq 3.95 \cdot 10^{+254}:\\
\;\;\;\;x \cdot \frac{\sinh y}{y}\\
\mathbf{else}:\\
\;\;\;\;\left(y \cdot \left(x \cdot \left(0.008333333333333333 + \left(x \cdot x\right) \cdot -0.001388888888888889\right)\right)\right) \cdot \left(y \cdot \left(y \cdot y\right)\right)\\
\end{array}
\end{array}
if y < 5.90000000000000038e-19Initial program 100.0%
Taylor expanded in y around 0
sin-lowering-sin.f6470.1%
Simplified70.1%
if 5.90000000000000038e-19 < y < 3.9499999999999999e254Initial program 100.0%
Taylor expanded in x around 0
Simplified83.0%
if 3.9499999999999999e254 < y Initial program 100.0%
Taylor expanded in y around 0
distribute-rgt-inN/A
*-rgt-identityN/A
distribute-rgt-inN/A
*-commutativeN/A
associate-*r*N/A
distribute-rgt-outN/A
+-commutativeN/A
associate-*l*N/A
*-commutativeN/A
distribute-lft-inN/A
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
+-lowering-+.f64N/A
Simplified100.0%
Taylor expanded in y around inf
associate-*r*N/A
*-commutativeN/A
metadata-evalN/A
pow-sqrN/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
associate-*r*N/A
unpow2N/A
unpow3N/A
*-lowering-*.f64N/A
unpow3N/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64100.0%
Simplified100.0%
Taylor expanded in x around 0
distribute-lft-inN/A
associate-*r*N/A
associate-*r*N/A
associate-*r*N/A
*-commutativeN/A
distribute-rgt-outN/A
*-lowering-*.f64N/A
metadata-evalN/A
pow-sqrN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-commutativeN/A
distribute-lft-outN/A
+-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64100.0%
Simplified100.0%
*-commutativeN/A
associate-*l*N/A
cube-unmultN/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
cube-unmultN/A
*-lowering-*.f64N/A
*-lowering-*.f64100.0%
Applied egg-rr100.0%
Final simplification74.6%
(FPCore (x y)
:precision binary64
(if (<= y 0.00095)
(sin x)
(if (<= y 2e+252)
(/
(+
1.0
(*
(* y y)
(+
0.16666666666666666
(*
(* y y)
(+ 0.008333333333333333 (* y (* y 0.0001984126984126984)))))))
(/ 1.0 x))
(*
(* y (* x (+ 0.008333333333333333 (* (* x x) -0.001388888888888889))))
(* y (* y y))))))
double code(double x, double y) {
double tmp;
if (y <= 0.00095) {
tmp = sin(x);
} else if (y <= 2e+252) {
tmp = (1.0 + ((y * y) * (0.16666666666666666 + ((y * y) * (0.008333333333333333 + (y * (y * 0.0001984126984126984))))))) / (1.0 / x);
} else {
tmp = (y * (x * (0.008333333333333333 + ((x * x) * -0.001388888888888889)))) * (y * (y * y));
}
return tmp;
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8) :: tmp
if (y <= 0.00095d0) then
tmp = sin(x)
else if (y <= 2d+252) then
tmp = (1.0d0 + ((y * y) * (0.16666666666666666d0 + ((y * y) * (0.008333333333333333d0 + (y * (y * 0.0001984126984126984d0))))))) / (1.0d0 / x)
else
tmp = (y * (x * (0.008333333333333333d0 + ((x * x) * (-0.001388888888888889d0))))) * (y * (y * y))
end if
code = tmp
end function
public static double code(double x, double y) {
double tmp;
if (y <= 0.00095) {
tmp = Math.sin(x);
} else if (y <= 2e+252) {
tmp = (1.0 + ((y * y) * (0.16666666666666666 + ((y * y) * (0.008333333333333333 + (y * (y * 0.0001984126984126984))))))) / (1.0 / x);
} else {
tmp = (y * (x * (0.008333333333333333 + ((x * x) * -0.001388888888888889)))) * (y * (y * y));
}
return tmp;
}
def code(x, y): tmp = 0 if y <= 0.00095: tmp = math.sin(x) elif y <= 2e+252: tmp = (1.0 + ((y * y) * (0.16666666666666666 + ((y * y) * (0.008333333333333333 + (y * (y * 0.0001984126984126984))))))) / (1.0 / x) else: tmp = (y * (x * (0.008333333333333333 + ((x * x) * -0.001388888888888889)))) * (y * (y * y)) return tmp
function code(x, y) tmp = 0.0 if (y <= 0.00095) tmp = sin(x); elseif (y <= 2e+252) tmp = Float64(Float64(1.0 + Float64(Float64(y * y) * Float64(0.16666666666666666 + Float64(Float64(y * y) * Float64(0.008333333333333333 + Float64(y * Float64(y * 0.0001984126984126984))))))) / Float64(1.0 / x)); else tmp = Float64(Float64(y * Float64(x * Float64(0.008333333333333333 + Float64(Float64(x * x) * -0.001388888888888889)))) * Float64(y * Float64(y * y))); end return tmp end
function tmp_2 = code(x, y) tmp = 0.0; if (y <= 0.00095) tmp = sin(x); elseif (y <= 2e+252) tmp = (1.0 + ((y * y) * (0.16666666666666666 + ((y * y) * (0.008333333333333333 + (y * (y * 0.0001984126984126984))))))) / (1.0 / x); else tmp = (y * (x * (0.008333333333333333 + ((x * x) * -0.001388888888888889)))) * (y * (y * y)); end tmp_2 = tmp; end
code[x_, y_] := If[LessEqual[y, 0.00095], N[Sin[x], $MachinePrecision], If[LessEqual[y, 2e+252], N[(N[(1.0 + N[(N[(y * y), $MachinePrecision] * N[(0.16666666666666666 + N[(N[(y * y), $MachinePrecision] * N[(0.008333333333333333 + N[(y * N[(y * 0.0001984126984126984), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(1.0 / x), $MachinePrecision]), $MachinePrecision], N[(N[(y * N[(x * N[(0.008333333333333333 + N[(N[(x * x), $MachinePrecision] * -0.001388888888888889), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(y * N[(y * y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq 0.00095:\\
\;\;\;\;\sin x\\
\mathbf{elif}\;y \leq 2 \cdot 10^{+252}:\\
\;\;\;\;\frac{1 + \left(y \cdot y\right) \cdot \left(0.16666666666666666 + \left(y \cdot y\right) \cdot \left(0.008333333333333333 + y \cdot \left(y \cdot 0.0001984126984126984\right)\right)\right)}{\frac{1}{x}}\\
\mathbf{else}:\\
\;\;\;\;\left(y \cdot \left(x \cdot \left(0.008333333333333333 + \left(x \cdot x\right) \cdot -0.001388888888888889\right)\right)\right) \cdot \left(y \cdot \left(y \cdot y\right)\right)\\
\end{array}
\end{array}
if y < 9.49999999999999998e-4Initial program 100.0%
Taylor expanded in y around 0
sin-lowering-sin.f6470.3%
Simplified70.3%
if 9.49999999999999998e-4 < y < 2.0000000000000002e252Initial program 100.0%
Taylor expanded in y around 0
Simplified82.4%
*-commutativeN/A
remove-double-divN/A
un-div-invN/A
/-lowering-/.f64N/A
Applied egg-rr82.4%
Taylor expanded in x around 0
Simplified70.1%
if 2.0000000000000002e252 < y Initial program 100.0%
Taylor expanded in y around 0
distribute-rgt-inN/A
*-rgt-identityN/A
distribute-rgt-inN/A
*-commutativeN/A
associate-*r*N/A
distribute-rgt-outN/A
+-commutativeN/A
associate-*l*N/A
*-commutativeN/A
distribute-lft-inN/A
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
+-lowering-+.f64N/A
Simplified100.0%
Taylor expanded in y around inf
associate-*r*N/A
*-commutativeN/A
metadata-evalN/A
pow-sqrN/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
associate-*r*N/A
unpow2N/A
unpow3N/A
*-lowering-*.f64N/A
unpow3N/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64100.0%
Simplified100.0%
Taylor expanded in x around 0
distribute-lft-inN/A
associate-*r*N/A
associate-*r*N/A
associate-*r*N/A
*-commutativeN/A
distribute-rgt-outN/A
*-lowering-*.f64N/A
metadata-evalN/A
pow-sqrN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-commutativeN/A
distribute-lft-outN/A
+-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64100.0%
Simplified100.0%
*-commutativeN/A
associate-*l*N/A
cube-unmultN/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
cube-unmultN/A
*-lowering-*.f64N/A
*-lowering-*.f64100.0%
Applied egg-rr100.0%
Final simplification71.7%
(FPCore (x y)
:precision binary64
(let* ((t_0
(*
x
(+
1.0
(*
(* x x)
(+
-0.16666666666666666
(*
x
(*
x
(+
0.008333333333333333
(* (* x x) -0.0001984126984126984))))))))))
(if (<= x 5.1e+36)
(/
(+
1.0
(*
(* y y)
(+
0.16666666666666666
(*
(* y y)
(+ 0.008333333333333333 (* y (* y 0.0001984126984126984)))))))
(/ (+ 1.0 (* 0.16666666666666666 (* x x))) x))
(if (<= x 6.8e+106)
t_0
(if (<= x 1.45e+235)
(*
x
(+
1.0
(*
x
(* x (+ -0.16666666666666666 (* 0.008333333333333333 (* x x)))))))
t_0)))))
double code(double x, double y) {
double t_0 = x * (1.0 + ((x * x) * (-0.16666666666666666 + (x * (x * (0.008333333333333333 + ((x * x) * -0.0001984126984126984)))))));
double tmp;
if (x <= 5.1e+36) {
tmp = (1.0 + ((y * y) * (0.16666666666666666 + ((y * y) * (0.008333333333333333 + (y * (y * 0.0001984126984126984))))))) / ((1.0 + (0.16666666666666666 * (x * x))) / x);
} else if (x <= 6.8e+106) {
tmp = t_0;
} else if (x <= 1.45e+235) {
tmp = x * (1.0 + (x * (x * (-0.16666666666666666 + (0.008333333333333333 * (x * x))))));
} else {
tmp = t_0;
}
return tmp;
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8) :: t_0
real(8) :: tmp
t_0 = x * (1.0d0 + ((x * x) * ((-0.16666666666666666d0) + (x * (x * (0.008333333333333333d0 + ((x * x) * (-0.0001984126984126984d0))))))))
if (x <= 5.1d+36) then
tmp = (1.0d0 + ((y * y) * (0.16666666666666666d0 + ((y * y) * (0.008333333333333333d0 + (y * (y * 0.0001984126984126984d0))))))) / ((1.0d0 + (0.16666666666666666d0 * (x * x))) / x)
else if (x <= 6.8d+106) then
tmp = t_0
else if (x <= 1.45d+235) then
tmp = x * (1.0d0 + (x * (x * ((-0.16666666666666666d0) + (0.008333333333333333d0 * (x * x))))))
else
tmp = t_0
end if
code = tmp
end function
public static double code(double x, double y) {
double t_0 = x * (1.0 + ((x * x) * (-0.16666666666666666 + (x * (x * (0.008333333333333333 + ((x * x) * -0.0001984126984126984)))))));
double tmp;
if (x <= 5.1e+36) {
tmp = (1.0 + ((y * y) * (0.16666666666666666 + ((y * y) * (0.008333333333333333 + (y * (y * 0.0001984126984126984))))))) / ((1.0 + (0.16666666666666666 * (x * x))) / x);
} else if (x <= 6.8e+106) {
tmp = t_0;
} else if (x <= 1.45e+235) {
tmp = x * (1.0 + (x * (x * (-0.16666666666666666 + (0.008333333333333333 * (x * x))))));
} else {
tmp = t_0;
}
return tmp;
}
def code(x, y): t_0 = x * (1.0 + ((x * x) * (-0.16666666666666666 + (x * (x * (0.008333333333333333 + ((x * x) * -0.0001984126984126984))))))) tmp = 0 if x <= 5.1e+36: tmp = (1.0 + ((y * y) * (0.16666666666666666 + ((y * y) * (0.008333333333333333 + (y * (y * 0.0001984126984126984))))))) / ((1.0 + (0.16666666666666666 * (x * x))) / x) elif x <= 6.8e+106: tmp = t_0 elif x <= 1.45e+235: tmp = x * (1.0 + (x * (x * (-0.16666666666666666 + (0.008333333333333333 * (x * x)))))) else: tmp = t_0 return tmp
function code(x, y) t_0 = Float64(x * Float64(1.0 + Float64(Float64(x * x) * Float64(-0.16666666666666666 + Float64(x * Float64(x * Float64(0.008333333333333333 + Float64(Float64(x * x) * -0.0001984126984126984)))))))) tmp = 0.0 if (x <= 5.1e+36) tmp = Float64(Float64(1.0 + Float64(Float64(y * y) * Float64(0.16666666666666666 + Float64(Float64(y * y) * Float64(0.008333333333333333 + Float64(y * Float64(y * 0.0001984126984126984))))))) / Float64(Float64(1.0 + Float64(0.16666666666666666 * Float64(x * x))) / x)); elseif (x <= 6.8e+106) tmp = t_0; elseif (x <= 1.45e+235) tmp = Float64(x * Float64(1.0 + Float64(x * Float64(x * Float64(-0.16666666666666666 + Float64(0.008333333333333333 * Float64(x * x))))))); else tmp = t_0; end return tmp end
function tmp_2 = code(x, y) t_0 = x * (1.0 + ((x * x) * (-0.16666666666666666 + (x * (x * (0.008333333333333333 + ((x * x) * -0.0001984126984126984))))))); tmp = 0.0; if (x <= 5.1e+36) tmp = (1.0 + ((y * y) * (0.16666666666666666 + ((y * y) * (0.008333333333333333 + (y * (y * 0.0001984126984126984))))))) / ((1.0 + (0.16666666666666666 * (x * x))) / x); elseif (x <= 6.8e+106) tmp = t_0; elseif (x <= 1.45e+235) tmp = x * (1.0 + (x * (x * (-0.16666666666666666 + (0.008333333333333333 * (x * x)))))); else tmp = t_0; end tmp_2 = tmp; end
code[x_, y_] := Block[{t$95$0 = N[(x * N[(1.0 + N[(N[(x * x), $MachinePrecision] * N[(-0.16666666666666666 + N[(x * N[(x * N[(0.008333333333333333 + N[(N[(x * x), $MachinePrecision] * -0.0001984126984126984), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[x, 5.1e+36], N[(N[(1.0 + N[(N[(y * y), $MachinePrecision] * N[(0.16666666666666666 + N[(N[(y * y), $MachinePrecision] * N[(0.008333333333333333 + N[(y * N[(y * 0.0001984126984126984), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(N[(1.0 + N[(0.16666666666666666 * N[(x * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / x), $MachinePrecision]), $MachinePrecision], If[LessEqual[x, 6.8e+106], t$95$0, If[LessEqual[x, 1.45e+235], N[(x * N[(1.0 + N[(x * N[(x * N[(-0.16666666666666666 + N[(0.008333333333333333 * N[(x * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$0]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := x \cdot \left(1 + \left(x \cdot x\right) \cdot \left(-0.16666666666666666 + x \cdot \left(x \cdot \left(0.008333333333333333 + \left(x \cdot x\right) \cdot -0.0001984126984126984\right)\right)\right)\right)\\
\mathbf{if}\;x \leq 5.1 \cdot 10^{+36}:\\
\;\;\;\;\frac{1 + \left(y \cdot y\right) \cdot \left(0.16666666666666666 + \left(y \cdot y\right) \cdot \left(0.008333333333333333 + y \cdot \left(y \cdot 0.0001984126984126984\right)\right)\right)}{\frac{1 + 0.16666666666666666 \cdot \left(x \cdot x\right)}{x}}\\
\mathbf{elif}\;x \leq 6.8 \cdot 10^{+106}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;x \leq 1.45 \cdot 10^{+235}:\\
\;\;\;\;x \cdot \left(1 + x \cdot \left(x \cdot \left(-0.16666666666666666 + 0.008333333333333333 \cdot \left(x \cdot x\right)\right)\right)\right)\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
if x < 5.09999999999999973e36Initial program 100.0%
Taylor expanded in y around 0
Simplified92.7%
*-commutativeN/A
remove-double-divN/A
un-div-invN/A
/-lowering-/.f64N/A
Applied egg-rr92.6%
Taylor expanded in x around 0
/-lowering-/.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6467.5%
Simplified67.5%
if 5.09999999999999973e36 < x < 6.79999999999999989e106 or 1.45000000000000011e235 < x Initial program 100.0%
Taylor expanded in y around 0
sin-lowering-sin.f6444.6%
Simplified44.6%
Taylor expanded in x around 0
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f64N/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6447.2%
Simplified47.2%
if 6.79999999999999989e106 < x < 1.45000000000000011e235Initial program 100.0%
Taylor expanded in y around 0
sin-lowering-sin.f6438.4%
Simplified38.4%
Taylor expanded in x around 0
*-lowering-*.f64N/A
+-lowering-+.f64N/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6455.5%
Simplified55.5%
Final simplification64.2%
(FPCore (x y)
:precision binary64
(let* ((t_0
(*
x
(+
1.0
(*
(* x x)
(+
-0.16666666666666666
(*
x
(*
x
(+
0.008333333333333333
(* (* x x) -0.0001984126984126984))))))))))
(if (<= x 5.1e+36)
(*
x
(+
1.0
(*
y
(*
y
(+
0.16666666666666666
(*
y
(*
y
(+ 0.008333333333333333 (* (* y y) 0.0001984126984126984)))))))))
(if (<= x 6.8e+106)
t_0
(if (<= x 1.45e+235)
(*
x
(+
1.0
(*
x
(* x (+ -0.16666666666666666 (* 0.008333333333333333 (* x x)))))))
t_0)))))
double code(double x, double y) {
double t_0 = x * (1.0 + ((x * x) * (-0.16666666666666666 + (x * (x * (0.008333333333333333 + ((x * x) * -0.0001984126984126984)))))));
double tmp;
if (x <= 5.1e+36) {
tmp = x * (1.0 + (y * (y * (0.16666666666666666 + (y * (y * (0.008333333333333333 + ((y * y) * 0.0001984126984126984))))))));
} else if (x <= 6.8e+106) {
tmp = t_0;
} else if (x <= 1.45e+235) {
tmp = x * (1.0 + (x * (x * (-0.16666666666666666 + (0.008333333333333333 * (x * x))))));
} else {
tmp = t_0;
}
return tmp;
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8) :: t_0
real(8) :: tmp
t_0 = x * (1.0d0 + ((x * x) * ((-0.16666666666666666d0) + (x * (x * (0.008333333333333333d0 + ((x * x) * (-0.0001984126984126984d0))))))))
if (x <= 5.1d+36) then
tmp = x * (1.0d0 + (y * (y * (0.16666666666666666d0 + (y * (y * (0.008333333333333333d0 + ((y * y) * 0.0001984126984126984d0))))))))
else if (x <= 6.8d+106) then
tmp = t_0
else if (x <= 1.45d+235) then
tmp = x * (1.0d0 + (x * (x * ((-0.16666666666666666d0) + (0.008333333333333333d0 * (x * x))))))
else
tmp = t_0
end if
code = tmp
end function
public static double code(double x, double y) {
double t_0 = x * (1.0 + ((x * x) * (-0.16666666666666666 + (x * (x * (0.008333333333333333 + ((x * x) * -0.0001984126984126984)))))));
double tmp;
if (x <= 5.1e+36) {
tmp = x * (1.0 + (y * (y * (0.16666666666666666 + (y * (y * (0.008333333333333333 + ((y * y) * 0.0001984126984126984))))))));
} else if (x <= 6.8e+106) {
tmp = t_0;
} else if (x <= 1.45e+235) {
tmp = x * (1.0 + (x * (x * (-0.16666666666666666 + (0.008333333333333333 * (x * x))))));
} else {
tmp = t_0;
}
return tmp;
}
def code(x, y): t_0 = x * (1.0 + ((x * x) * (-0.16666666666666666 + (x * (x * (0.008333333333333333 + ((x * x) * -0.0001984126984126984))))))) tmp = 0 if x <= 5.1e+36: tmp = x * (1.0 + (y * (y * (0.16666666666666666 + (y * (y * (0.008333333333333333 + ((y * y) * 0.0001984126984126984)))))))) elif x <= 6.8e+106: tmp = t_0 elif x <= 1.45e+235: tmp = x * (1.0 + (x * (x * (-0.16666666666666666 + (0.008333333333333333 * (x * x)))))) else: tmp = t_0 return tmp
function code(x, y) t_0 = Float64(x * Float64(1.0 + Float64(Float64(x * x) * Float64(-0.16666666666666666 + Float64(x * Float64(x * Float64(0.008333333333333333 + Float64(Float64(x * x) * -0.0001984126984126984)))))))) tmp = 0.0 if (x <= 5.1e+36) tmp = Float64(x * Float64(1.0 + Float64(y * Float64(y * Float64(0.16666666666666666 + Float64(y * Float64(y * Float64(0.008333333333333333 + Float64(Float64(y * y) * 0.0001984126984126984))))))))); elseif (x <= 6.8e+106) tmp = t_0; elseif (x <= 1.45e+235) tmp = Float64(x * Float64(1.0 + Float64(x * Float64(x * Float64(-0.16666666666666666 + Float64(0.008333333333333333 * Float64(x * x))))))); else tmp = t_0; end return tmp end
function tmp_2 = code(x, y) t_0 = x * (1.0 + ((x * x) * (-0.16666666666666666 + (x * (x * (0.008333333333333333 + ((x * x) * -0.0001984126984126984))))))); tmp = 0.0; if (x <= 5.1e+36) tmp = x * (1.0 + (y * (y * (0.16666666666666666 + (y * (y * (0.008333333333333333 + ((y * y) * 0.0001984126984126984)))))))); elseif (x <= 6.8e+106) tmp = t_0; elseif (x <= 1.45e+235) tmp = x * (1.0 + (x * (x * (-0.16666666666666666 + (0.008333333333333333 * (x * x)))))); else tmp = t_0; end tmp_2 = tmp; end
code[x_, y_] := Block[{t$95$0 = N[(x * N[(1.0 + N[(N[(x * x), $MachinePrecision] * N[(-0.16666666666666666 + N[(x * N[(x * N[(0.008333333333333333 + N[(N[(x * x), $MachinePrecision] * -0.0001984126984126984), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[x, 5.1e+36], N[(x * N[(1.0 + N[(y * N[(y * N[(0.16666666666666666 + N[(y * N[(y * N[(0.008333333333333333 + N[(N[(y * y), $MachinePrecision] * 0.0001984126984126984), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[x, 6.8e+106], t$95$0, If[LessEqual[x, 1.45e+235], N[(x * N[(1.0 + N[(x * N[(x * N[(-0.16666666666666666 + N[(0.008333333333333333 * N[(x * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], t$95$0]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := x \cdot \left(1 + \left(x \cdot x\right) \cdot \left(-0.16666666666666666 + x \cdot \left(x \cdot \left(0.008333333333333333 + \left(x \cdot x\right) \cdot -0.0001984126984126984\right)\right)\right)\right)\\
\mathbf{if}\;x \leq 5.1 \cdot 10^{+36}:\\
\;\;\;\;x \cdot \left(1 + y \cdot \left(y \cdot \left(0.16666666666666666 + y \cdot \left(y \cdot \left(0.008333333333333333 + \left(y \cdot y\right) \cdot 0.0001984126984126984\right)\right)\right)\right)\right)\\
\mathbf{elif}\;x \leq 6.8 \cdot 10^{+106}:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;x \leq 1.45 \cdot 10^{+235}:\\
\;\;\;\;x \cdot \left(1 + x \cdot \left(x \cdot \left(-0.16666666666666666 + 0.008333333333333333 \cdot \left(x \cdot x\right)\right)\right)\right)\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
if x < 5.09999999999999973e36Initial program 100.0%
Taylor expanded in y around 0
Simplified92.7%
Taylor expanded in x around 0
*-lowering-*.f64N/A
+-lowering-+.f64N/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6471.0%
Simplified71.0%
if 5.09999999999999973e36 < x < 6.79999999999999989e106 or 1.45000000000000011e235 < x Initial program 100.0%
Taylor expanded in y around 0
sin-lowering-sin.f6444.6%
Simplified44.6%
Taylor expanded in x around 0
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f64N/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6447.2%
Simplified47.2%
if 6.79999999999999989e106 < x < 1.45000000000000011e235Initial program 100.0%
Taylor expanded in y around 0
sin-lowering-sin.f6438.4%
Simplified38.4%
Taylor expanded in x around 0
*-lowering-*.f64N/A
+-lowering-+.f64N/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6455.5%
Simplified55.5%
Final simplification67.0%
(FPCore (x y)
:precision binary64
(if (<= y 3.7e+252)
(*
x
(+
1.0
(*
y
(*
y
(+
0.16666666666666666
(*
y
(* y (+ 0.008333333333333333 (* (* y y) 0.0001984126984126984)))))))))
(*
(* y (* x (+ 0.008333333333333333 (* (* x x) -0.001388888888888889))))
(* y (* y y)))))
double code(double x, double y) {
double tmp;
if (y <= 3.7e+252) {
tmp = x * (1.0 + (y * (y * (0.16666666666666666 + (y * (y * (0.008333333333333333 + ((y * y) * 0.0001984126984126984))))))));
} else {
tmp = (y * (x * (0.008333333333333333 + ((x * x) * -0.001388888888888889)))) * (y * (y * y));
}
return tmp;
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8) :: tmp
if (y <= 3.7d+252) then
tmp = x * (1.0d0 + (y * (y * (0.16666666666666666d0 + (y * (y * (0.008333333333333333d0 + ((y * y) * 0.0001984126984126984d0))))))))
else
tmp = (y * (x * (0.008333333333333333d0 + ((x * x) * (-0.001388888888888889d0))))) * (y * (y * y))
end if
code = tmp
end function
public static double code(double x, double y) {
double tmp;
if (y <= 3.7e+252) {
tmp = x * (1.0 + (y * (y * (0.16666666666666666 + (y * (y * (0.008333333333333333 + ((y * y) * 0.0001984126984126984))))))));
} else {
tmp = (y * (x * (0.008333333333333333 + ((x * x) * -0.001388888888888889)))) * (y * (y * y));
}
return tmp;
}
def code(x, y): tmp = 0 if y <= 3.7e+252: tmp = x * (1.0 + (y * (y * (0.16666666666666666 + (y * (y * (0.008333333333333333 + ((y * y) * 0.0001984126984126984)))))))) else: tmp = (y * (x * (0.008333333333333333 + ((x * x) * -0.001388888888888889)))) * (y * (y * y)) return tmp
function code(x, y) tmp = 0.0 if (y <= 3.7e+252) tmp = Float64(x * Float64(1.0 + Float64(y * Float64(y * Float64(0.16666666666666666 + Float64(y * Float64(y * Float64(0.008333333333333333 + Float64(Float64(y * y) * 0.0001984126984126984))))))))); else tmp = Float64(Float64(y * Float64(x * Float64(0.008333333333333333 + Float64(Float64(x * x) * -0.001388888888888889)))) * Float64(y * Float64(y * y))); end return tmp end
function tmp_2 = code(x, y) tmp = 0.0; if (y <= 3.7e+252) tmp = x * (1.0 + (y * (y * (0.16666666666666666 + (y * (y * (0.008333333333333333 + ((y * y) * 0.0001984126984126984)))))))); else tmp = (y * (x * (0.008333333333333333 + ((x * x) * -0.001388888888888889)))) * (y * (y * y)); end tmp_2 = tmp; end
code[x_, y_] := If[LessEqual[y, 3.7e+252], N[(x * N[(1.0 + N[(y * N[(y * N[(0.16666666666666666 + N[(y * N[(y * N[(0.008333333333333333 + N[(N[(y * y), $MachinePrecision] * 0.0001984126984126984), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(y * N[(x * N[(0.008333333333333333 + N[(N[(x * x), $MachinePrecision] * -0.001388888888888889), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(y * N[(y * y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq 3.7 \cdot 10^{+252}:\\
\;\;\;\;x \cdot \left(1 + y \cdot \left(y \cdot \left(0.16666666666666666 + y \cdot \left(y \cdot \left(0.008333333333333333 + \left(y \cdot y\right) \cdot 0.0001984126984126984\right)\right)\right)\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\left(y \cdot \left(x \cdot \left(0.008333333333333333 + \left(x \cdot x\right) \cdot -0.001388888888888889\right)\right)\right) \cdot \left(y \cdot \left(y \cdot y\right)\right)\\
\end{array}
\end{array}
if y < 3.6999999999999998e252Initial program 100.0%
Taylor expanded in y around 0
Simplified91.9%
Taylor expanded in x around 0
*-lowering-*.f64N/A
+-lowering-+.f64N/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6462.9%
Simplified62.9%
if 3.6999999999999998e252 < y Initial program 100.0%
Taylor expanded in y around 0
distribute-rgt-inN/A
*-rgt-identityN/A
distribute-rgt-inN/A
*-commutativeN/A
associate-*r*N/A
distribute-rgt-outN/A
+-commutativeN/A
associate-*l*N/A
*-commutativeN/A
distribute-lft-inN/A
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
+-lowering-+.f64N/A
Simplified100.0%
Taylor expanded in y around inf
associate-*r*N/A
*-commutativeN/A
metadata-evalN/A
pow-sqrN/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
associate-*r*N/A
unpow2N/A
unpow3N/A
*-lowering-*.f64N/A
unpow3N/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64100.0%
Simplified100.0%
Taylor expanded in x around 0
distribute-lft-inN/A
associate-*r*N/A
associate-*r*N/A
associate-*r*N/A
*-commutativeN/A
distribute-rgt-outN/A
*-lowering-*.f64N/A
metadata-evalN/A
pow-sqrN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-commutativeN/A
distribute-lft-outN/A
+-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64100.0%
Simplified100.0%
*-commutativeN/A
associate-*l*N/A
cube-unmultN/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
cube-unmultN/A
*-lowering-*.f64N/A
*-lowering-*.f64100.0%
Applied egg-rr100.0%
Final simplification64.8%
(FPCore (x y)
:precision binary64
(if (<= y 2.3e+18)
(*
x
(+
1.0
(* x (* x (+ -0.16666666666666666 (* 0.008333333333333333 (* x x)))))))
(if (<= y 1e+237)
(* x (* y (* y (* (* y y) 0.008333333333333333))))
(*
(* y y)
(* x (+ 0.16666666666666666 (* (* x x) -0.027777777777777776)))))))
double code(double x, double y) {
double tmp;
if (y <= 2.3e+18) {
tmp = x * (1.0 + (x * (x * (-0.16666666666666666 + (0.008333333333333333 * (x * x))))));
} else if (y <= 1e+237) {
tmp = x * (y * (y * ((y * y) * 0.008333333333333333)));
} else {
tmp = (y * y) * (x * (0.16666666666666666 + ((x * x) * -0.027777777777777776)));
}
return tmp;
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8) :: tmp
if (y <= 2.3d+18) then
tmp = x * (1.0d0 + (x * (x * ((-0.16666666666666666d0) + (0.008333333333333333d0 * (x * x))))))
else if (y <= 1d+237) then
tmp = x * (y * (y * ((y * y) * 0.008333333333333333d0)))
else
tmp = (y * y) * (x * (0.16666666666666666d0 + ((x * x) * (-0.027777777777777776d0))))
end if
code = tmp
end function
public static double code(double x, double y) {
double tmp;
if (y <= 2.3e+18) {
tmp = x * (1.0 + (x * (x * (-0.16666666666666666 + (0.008333333333333333 * (x * x))))));
} else if (y <= 1e+237) {
tmp = x * (y * (y * ((y * y) * 0.008333333333333333)));
} else {
tmp = (y * y) * (x * (0.16666666666666666 + ((x * x) * -0.027777777777777776)));
}
return tmp;
}
def code(x, y): tmp = 0 if y <= 2.3e+18: tmp = x * (1.0 + (x * (x * (-0.16666666666666666 + (0.008333333333333333 * (x * x)))))) elif y <= 1e+237: tmp = x * (y * (y * ((y * y) * 0.008333333333333333))) else: tmp = (y * y) * (x * (0.16666666666666666 + ((x * x) * -0.027777777777777776))) return tmp
function code(x, y) tmp = 0.0 if (y <= 2.3e+18) tmp = Float64(x * Float64(1.0 + Float64(x * Float64(x * Float64(-0.16666666666666666 + Float64(0.008333333333333333 * Float64(x * x))))))); elseif (y <= 1e+237) tmp = Float64(x * Float64(y * Float64(y * Float64(Float64(y * y) * 0.008333333333333333)))); else tmp = Float64(Float64(y * y) * Float64(x * Float64(0.16666666666666666 + Float64(Float64(x * x) * -0.027777777777777776)))); end return tmp end
function tmp_2 = code(x, y) tmp = 0.0; if (y <= 2.3e+18) tmp = x * (1.0 + (x * (x * (-0.16666666666666666 + (0.008333333333333333 * (x * x)))))); elseif (y <= 1e+237) tmp = x * (y * (y * ((y * y) * 0.008333333333333333))); else tmp = (y * y) * (x * (0.16666666666666666 + ((x * x) * -0.027777777777777776))); end tmp_2 = tmp; end
code[x_, y_] := If[LessEqual[y, 2.3e+18], N[(x * N[(1.0 + N[(x * N[(x * N[(-0.16666666666666666 + N[(0.008333333333333333 * N[(x * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 1e+237], N[(x * N[(y * N[(y * N[(N[(y * y), $MachinePrecision] * 0.008333333333333333), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(y * y), $MachinePrecision] * N[(x * N[(0.16666666666666666 + N[(N[(x * x), $MachinePrecision] * -0.027777777777777776), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq 2.3 \cdot 10^{+18}:\\
\;\;\;\;x \cdot \left(1 + x \cdot \left(x \cdot \left(-0.16666666666666666 + 0.008333333333333333 \cdot \left(x \cdot x\right)\right)\right)\right)\\
\mathbf{elif}\;y \leq 10^{+237}:\\
\;\;\;\;x \cdot \left(y \cdot \left(y \cdot \left(\left(y \cdot y\right) \cdot 0.008333333333333333\right)\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\left(y \cdot y\right) \cdot \left(x \cdot \left(0.16666666666666666 + \left(x \cdot x\right) \cdot -0.027777777777777776\right)\right)\\
\end{array}
\end{array}
if y < 2.3e18Initial program 100.0%
Taylor expanded in y around 0
sin-lowering-sin.f6468.5%
Simplified68.5%
Taylor expanded in x around 0
*-lowering-*.f64N/A
+-lowering-+.f64N/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6446.2%
Simplified46.2%
if 2.3e18 < y < 9.9999999999999994e236Initial program 100.0%
Taylor expanded in y around 0
distribute-rgt-inN/A
*-rgt-identityN/A
distribute-rgt-inN/A
*-commutativeN/A
associate-*r*N/A
distribute-rgt-outN/A
+-commutativeN/A
associate-*l*N/A
*-commutativeN/A
distribute-lft-inN/A
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
+-lowering-+.f64N/A
Simplified87.2%
Taylor expanded in y around inf
associate-*r*N/A
*-commutativeN/A
metadata-evalN/A
pow-sqrN/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
associate-*r*N/A
unpow2N/A
unpow3N/A
*-lowering-*.f64N/A
unpow3N/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6487.2%
Simplified87.2%
Taylor expanded in x around 0
Simplified73.6%
if 9.9999999999999994e236 < y Initial program 100.0%
Taylor expanded in y around 0
*-lft-identityN/A
associate-*r*N/A
distribute-rgt-inN/A
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64100.0%
Simplified100.0%
Taylor expanded in y around inf
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64100.0%
Simplified100.0%
Taylor expanded in x around 0
distribute-lft-inN/A
associate-*r*N/A
associate-*r*N/A
associate-*r*N/A
*-commutativeN/A
distribute-rgt-outN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-commutativeN/A
distribute-lft-outN/A
*-lowering-*.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64100.0%
Simplified100.0%
Final simplification55.2%
(FPCore (x y)
:precision binary64
(if (<= y 7.5e+14)
(* x (+ 1.0 (* x (* x -0.16666666666666666))))
(if (<= y 6.2e+235)
(* x (* y (* y (* (* y y) 0.008333333333333333))))
(*
(* y y)
(* x (+ 0.16666666666666666 (* (* x x) -0.027777777777777776)))))))
double code(double x, double y) {
double tmp;
if (y <= 7.5e+14) {
tmp = x * (1.0 + (x * (x * -0.16666666666666666)));
} else if (y <= 6.2e+235) {
tmp = x * (y * (y * ((y * y) * 0.008333333333333333)));
} else {
tmp = (y * y) * (x * (0.16666666666666666 + ((x * x) * -0.027777777777777776)));
}
return tmp;
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8) :: tmp
if (y <= 7.5d+14) then
tmp = x * (1.0d0 + (x * (x * (-0.16666666666666666d0))))
else if (y <= 6.2d+235) then
tmp = x * (y * (y * ((y * y) * 0.008333333333333333d0)))
else
tmp = (y * y) * (x * (0.16666666666666666d0 + ((x * x) * (-0.027777777777777776d0))))
end if
code = tmp
end function
public static double code(double x, double y) {
double tmp;
if (y <= 7.5e+14) {
tmp = x * (1.0 + (x * (x * -0.16666666666666666)));
} else if (y <= 6.2e+235) {
tmp = x * (y * (y * ((y * y) * 0.008333333333333333)));
} else {
tmp = (y * y) * (x * (0.16666666666666666 + ((x * x) * -0.027777777777777776)));
}
return tmp;
}
def code(x, y): tmp = 0 if y <= 7.5e+14: tmp = x * (1.0 + (x * (x * -0.16666666666666666))) elif y <= 6.2e+235: tmp = x * (y * (y * ((y * y) * 0.008333333333333333))) else: tmp = (y * y) * (x * (0.16666666666666666 + ((x * x) * -0.027777777777777776))) return tmp
function code(x, y) tmp = 0.0 if (y <= 7.5e+14) tmp = Float64(x * Float64(1.0 + Float64(x * Float64(x * -0.16666666666666666)))); elseif (y <= 6.2e+235) tmp = Float64(x * Float64(y * Float64(y * Float64(Float64(y * y) * 0.008333333333333333)))); else tmp = Float64(Float64(y * y) * Float64(x * Float64(0.16666666666666666 + Float64(Float64(x * x) * -0.027777777777777776)))); end return tmp end
function tmp_2 = code(x, y) tmp = 0.0; if (y <= 7.5e+14) tmp = x * (1.0 + (x * (x * -0.16666666666666666))); elseif (y <= 6.2e+235) tmp = x * (y * (y * ((y * y) * 0.008333333333333333))); else tmp = (y * y) * (x * (0.16666666666666666 + ((x * x) * -0.027777777777777776))); end tmp_2 = tmp; end
code[x_, y_] := If[LessEqual[y, 7.5e+14], N[(x * N[(1.0 + N[(x * N[(x * -0.16666666666666666), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 6.2e+235], N[(x * N[(y * N[(y * N[(N[(y * y), $MachinePrecision] * 0.008333333333333333), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(y * y), $MachinePrecision] * N[(x * N[(0.16666666666666666 + N[(N[(x * x), $MachinePrecision] * -0.027777777777777776), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq 7.5 \cdot 10^{+14}:\\
\;\;\;\;x \cdot \left(1 + x \cdot \left(x \cdot -0.16666666666666666\right)\right)\\
\mathbf{elif}\;y \leq 6.2 \cdot 10^{+235}:\\
\;\;\;\;x \cdot \left(y \cdot \left(y \cdot \left(\left(y \cdot y\right) \cdot 0.008333333333333333\right)\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\left(y \cdot y\right) \cdot \left(x \cdot \left(0.16666666666666666 + \left(x \cdot x\right) \cdot -0.027777777777777776\right)\right)\\
\end{array}
\end{array}
if y < 7.5e14Initial program 100.0%
Taylor expanded in y around 0
sin-lowering-sin.f6468.5%
Simplified68.5%
Taylor expanded in x around 0
*-lowering-*.f64N/A
+-lowering-+.f64N/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6446.2%
Simplified46.2%
Taylor expanded in x around 0
*-commutativeN/A
*-lowering-*.f6444.1%
Simplified44.1%
if 7.5e14 < y < 6.20000000000000022e235Initial program 100.0%
Taylor expanded in y around 0
distribute-rgt-inN/A
*-rgt-identityN/A
distribute-rgt-inN/A
*-commutativeN/A
associate-*r*N/A
distribute-rgt-outN/A
+-commutativeN/A
associate-*l*N/A
*-commutativeN/A
distribute-lft-inN/A
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
+-lowering-+.f64N/A
Simplified87.2%
Taylor expanded in y around inf
associate-*r*N/A
*-commutativeN/A
metadata-evalN/A
pow-sqrN/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
associate-*r*N/A
unpow2N/A
unpow3N/A
*-lowering-*.f64N/A
unpow3N/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6487.2%
Simplified87.2%
Taylor expanded in x around 0
Simplified73.6%
if 6.20000000000000022e235 < y Initial program 100.0%
Taylor expanded in y around 0
*-lft-identityN/A
associate-*r*N/A
distribute-rgt-inN/A
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64100.0%
Simplified100.0%
Taylor expanded in y around inf
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64100.0%
Simplified100.0%
Taylor expanded in x around 0
distribute-lft-inN/A
associate-*r*N/A
associate-*r*N/A
associate-*r*N/A
*-commutativeN/A
distribute-rgt-outN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-commutativeN/A
distribute-lft-outN/A
*-lowering-*.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64100.0%
Simplified100.0%
(FPCore (x y)
:precision binary64
(if (<= y 1.6e+255)
(*
x
(+
1.0
(* (* y y) (+ (* (* y y) 0.008333333333333333) 0.16666666666666666))))
(*
(* y (* x (+ 0.008333333333333333 (* (* x x) -0.001388888888888889))))
(* y (* y y)))))
double code(double x, double y) {
double tmp;
if (y <= 1.6e+255) {
tmp = x * (1.0 + ((y * y) * (((y * y) * 0.008333333333333333) + 0.16666666666666666)));
} else {
tmp = (y * (x * (0.008333333333333333 + ((x * x) * -0.001388888888888889)))) * (y * (y * y));
}
return tmp;
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8) :: tmp
if (y <= 1.6d+255) then
tmp = x * (1.0d0 + ((y * y) * (((y * y) * 0.008333333333333333d0) + 0.16666666666666666d0)))
else
tmp = (y * (x * (0.008333333333333333d0 + ((x * x) * (-0.001388888888888889d0))))) * (y * (y * y))
end if
code = tmp
end function
public static double code(double x, double y) {
double tmp;
if (y <= 1.6e+255) {
tmp = x * (1.0 + ((y * y) * (((y * y) * 0.008333333333333333) + 0.16666666666666666)));
} else {
tmp = (y * (x * (0.008333333333333333 + ((x * x) * -0.001388888888888889)))) * (y * (y * y));
}
return tmp;
}
def code(x, y): tmp = 0 if y <= 1.6e+255: tmp = x * (1.0 + ((y * y) * (((y * y) * 0.008333333333333333) + 0.16666666666666666))) else: tmp = (y * (x * (0.008333333333333333 + ((x * x) * -0.001388888888888889)))) * (y * (y * y)) return tmp
function code(x, y) tmp = 0.0 if (y <= 1.6e+255) tmp = Float64(x * Float64(1.0 + Float64(Float64(y * y) * Float64(Float64(Float64(y * y) * 0.008333333333333333) + 0.16666666666666666)))); else tmp = Float64(Float64(y * Float64(x * Float64(0.008333333333333333 + Float64(Float64(x * x) * -0.001388888888888889)))) * Float64(y * Float64(y * y))); end return tmp end
function tmp_2 = code(x, y) tmp = 0.0; if (y <= 1.6e+255) tmp = x * (1.0 + ((y * y) * (((y * y) * 0.008333333333333333) + 0.16666666666666666))); else tmp = (y * (x * (0.008333333333333333 + ((x * x) * -0.001388888888888889)))) * (y * (y * y)); end tmp_2 = tmp; end
code[x_, y_] := If[LessEqual[y, 1.6e+255], N[(x * N[(1.0 + N[(N[(y * y), $MachinePrecision] * N[(N[(N[(y * y), $MachinePrecision] * 0.008333333333333333), $MachinePrecision] + 0.16666666666666666), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(y * N[(x * N[(0.008333333333333333 + N[(N[(x * x), $MachinePrecision] * -0.001388888888888889), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(y * N[(y * y), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq 1.6 \cdot 10^{+255}:\\
\;\;\;\;x \cdot \left(1 + \left(y \cdot y\right) \cdot \left(\left(y \cdot y\right) \cdot 0.008333333333333333 + 0.16666666666666666\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\left(y \cdot \left(x \cdot \left(0.008333333333333333 + \left(x \cdot x\right) \cdot -0.001388888888888889\right)\right)\right) \cdot \left(y \cdot \left(y \cdot y\right)\right)\\
\end{array}
\end{array}
if y < 1.5999999999999999e255Initial program 100.0%
Taylor expanded in y around 0
distribute-rgt-inN/A
*-rgt-identityN/A
distribute-rgt-inN/A
*-commutativeN/A
associate-*r*N/A
distribute-rgt-outN/A
+-commutativeN/A
associate-*l*N/A
*-commutativeN/A
distribute-lft-inN/A
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
+-lowering-+.f64N/A
Simplified90.6%
Taylor expanded in x around 0
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6462.1%
Simplified62.1%
if 1.5999999999999999e255 < y Initial program 100.0%
Taylor expanded in y around 0
distribute-rgt-inN/A
*-rgt-identityN/A
distribute-rgt-inN/A
*-commutativeN/A
associate-*r*N/A
distribute-rgt-outN/A
+-commutativeN/A
associate-*l*N/A
*-commutativeN/A
distribute-lft-inN/A
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
+-lowering-+.f64N/A
Simplified100.0%
Taylor expanded in y around inf
associate-*r*N/A
*-commutativeN/A
metadata-evalN/A
pow-sqrN/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
associate-*r*N/A
unpow2N/A
unpow3N/A
*-lowering-*.f64N/A
unpow3N/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64100.0%
Simplified100.0%
Taylor expanded in x around 0
distribute-lft-inN/A
associate-*r*N/A
associate-*r*N/A
associate-*r*N/A
*-commutativeN/A
distribute-rgt-outN/A
*-lowering-*.f64N/A
metadata-evalN/A
pow-sqrN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-commutativeN/A
distribute-lft-outN/A
+-commutativeN/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64100.0%
Simplified100.0%
*-commutativeN/A
associate-*l*N/A
cube-unmultN/A
associate-*r*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
cube-unmultN/A
*-lowering-*.f64N/A
*-lowering-*.f64100.0%
Applied egg-rr100.0%
Final simplification64.0%
(FPCore (x y)
:precision binary64
(if (<= x 6.5e+127)
(* x (+ 1.0 (* (* y y) 0.16666666666666666)))
(if (<= x 1.45e+235)
(* x (* 0.008333333333333333 (* (* x x) (* x x))))
(* x (+ 1.0 (* x (* x -0.16666666666666666)))))))
double code(double x, double y) {
double tmp;
if (x <= 6.5e+127) {
tmp = x * (1.0 + ((y * y) * 0.16666666666666666));
} else if (x <= 1.45e+235) {
tmp = x * (0.008333333333333333 * ((x * x) * (x * x)));
} else {
tmp = x * (1.0 + (x * (x * -0.16666666666666666)));
}
return tmp;
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8) :: tmp
if (x <= 6.5d+127) then
tmp = x * (1.0d0 + ((y * y) * 0.16666666666666666d0))
else if (x <= 1.45d+235) then
tmp = x * (0.008333333333333333d0 * ((x * x) * (x * x)))
else
tmp = x * (1.0d0 + (x * (x * (-0.16666666666666666d0))))
end if
code = tmp
end function
public static double code(double x, double y) {
double tmp;
if (x <= 6.5e+127) {
tmp = x * (1.0 + ((y * y) * 0.16666666666666666));
} else if (x <= 1.45e+235) {
tmp = x * (0.008333333333333333 * ((x * x) * (x * x)));
} else {
tmp = x * (1.0 + (x * (x * -0.16666666666666666)));
}
return tmp;
}
def code(x, y): tmp = 0 if x <= 6.5e+127: tmp = x * (1.0 + ((y * y) * 0.16666666666666666)) elif x <= 1.45e+235: tmp = x * (0.008333333333333333 * ((x * x) * (x * x))) else: tmp = x * (1.0 + (x * (x * -0.16666666666666666))) return tmp
function code(x, y) tmp = 0.0 if (x <= 6.5e+127) tmp = Float64(x * Float64(1.0 + Float64(Float64(y * y) * 0.16666666666666666))); elseif (x <= 1.45e+235) tmp = Float64(x * Float64(0.008333333333333333 * Float64(Float64(x * x) * Float64(x * x)))); else tmp = Float64(x * Float64(1.0 + Float64(x * Float64(x * -0.16666666666666666)))); end return tmp end
function tmp_2 = code(x, y) tmp = 0.0; if (x <= 6.5e+127) tmp = x * (1.0 + ((y * y) * 0.16666666666666666)); elseif (x <= 1.45e+235) tmp = x * (0.008333333333333333 * ((x * x) * (x * x))); else tmp = x * (1.0 + (x * (x * -0.16666666666666666))); end tmp_2 = tmp; end
code[x_, y_] := If[LessEqual[x, 6.5e+127], N[(x * N[(1.0 + N[(N[(y * y), $MachinePrecision] * 0.16666666666666666), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[x, 1.45e+235], N[(x * N[(0.008333333333333333 * N[(N[(x * x), $MachinePrecision] * N[(x * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x * N[(1.0 + N[(x * N[(x * -0.16666666666666666), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq 6.5 \cdot 10^{+127}:\\
\;\;\;\;x \cdot \left(1 + \left(y \cdot y\right) \cdot 0.16666666666666666\right)\\
\mathbf{elif}\;x \leq 1.45 \cdot 10^{+235}:\\
\;\;\;\;x \cdot \left(0.008333333333333333 \cdot \left(\left(x \cdot x\right) \cdot \left(x \cdot x\right)\right)\right)\\
\mathbf{else}:\\
\;\;\;\;x \cdot \left(1 + x \cdot \left(x \cdot -0.16666666666666666\right)\right)\\
\end{array}
\end{array}
if x < 6.5e127Initial program 100.0%
Taylor expanded in y around 0
*-lft-identityN/A
associate-*r*N/A
distribute-rgt-inN/A
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f6478.1%
Simplified78.1%
Taylor expanded in x around 0
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6457.0%
Simplified57.0%
if 6.5e127 < x < 1.45000000000000011e235Initial program 100.0%
Taylor expanded in y around 0
sin-lowering-sin.f6438.4%
Simplified38.4%
Taylor expanded in x around 0
*-lowering-*.f64N/A
+-lowering-+.f64N/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6455.5%
Simplified55.5%
Taylor expanded in x around inf
metadata-evalN/A
pow-sqrN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6455.5%
Simplified55.5%
if 1.45000000000000011e235 < x Initial program 100.0%
Taylor expanded in y around 0
sin-lowering-sin.f6448.3%
Simplified48.3%
Taylor expanded in x around 0
*-lowering-*.f64N/A
+-lowering-+.f64N/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f647.6%
Simplified7.6%
Taylor expanded in x around 0
*-commutativeN/A
*-lowering-*.f6447.5%
Simplified47.5%
Final simplification56.3%
(FPCore (x y)
:precision binary64
(if (<= y 4e+236)
(*
x
(+
1.0
(* (* y y) (+ (* (* y y) 0.008333333333333333) 0.16666666666666666))))
(*
(* y y)
(* x (+ 0.16666666666666666 (* (* x x) -0.027777777777777776))))))
double code(double x, double y) {
double tmp;
if (y <= 4e+236) {
tmp = x * (1.0 + ((y * y) * (((y * y) * 0.008333333333333333) + 0.16666666666666666)));
} else {
tmp = (y * y) * (x * (0.16666666666666666 + ((x * x) * -0.027777777777777776)));
}
return tmp;
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8) :: tmp
if (y <= 4d+236) then
tmp = x * (1.0d0 + ((y * y) * (((y * y) * 0.008333333333333333d0) + 0.16666666666666666d0)))
else
tmp = (y * y) * (x * (0.16666666666666666d0 + ((x * x) * (-0.027777777777777776d0))))
end if
code = tmp
end function
public static double code(double x, double y) {
double tmp;
if (y <= 4e+236) {
tmp = x * (1.0 + ((y * y) * (((y * y) * 0.008333333333333333) + 0.16666666666666666)));
} else {
tmp = (y * y) * (x * (0.16666666666666666 + ((x * x) * -0.027777777777777776)));
}
return tmp;
}
def code(x, y): tmp = 0 if y <= 4e+236: tmp = x * (1.0 + ((y * y) * (((y * y) * 0.008333333333333333) + 0.16666666666666666))) else: tmp = (y * y) * (x * (0.16666666666666666 + ((x * x) * -0.027777777777777776))) return tmp
function code(x, y) tmp = 0.0 if (y <= 4e+236) tmp = Float64(x * Float64(1.0 + Float64(Float64(y * y) * Float64(Float64(Float64(y * y) * 0.008333333333333333) + 0.16666666666666666)))); else tmp = Float64(Float64(y * y) * Float64(x * Float64(0.16666666666666666 + Float64(Float64(x * x) * -0.027777777777777776)))); end return tmp end
function tmp_2 = code(x, y) tmp = 0.0; if (y <= 4e+236) tmp = x * (1.0 + ((y * y) * (((y * y) * 0.008333333333333333) + 0.16666666666666666))); else tmp = (y * y) * (x * (0.16666666666666666 + ((x * x) * -0.027777777777777776))); end tmp_2 = tmp; end
code[x_, y_] := If[LessEqual[y, 4e+236], N[(x * N[(1.0 + N[(N[(y * y), $MachinePrecision] * N[(N[(N[(y * y), $MachinePrecision] * 0.008333333333333333), $MachinePrecision] + 0.16666666666666666), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(y * y), $MachinePrecision] * N[(x * N[(0.16666666666666666 + N[(N[(x * x), $MachinePrecision] * -0.027777777777777776), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq 4 \cdot 10^{+236}:\\
\;\;\;\;x \cdot \left(1 + \left(y \cdot y\right) \cdot \left(\left(y \cdot y\right) \cdot 0.008333333333333333 + 0.16666666666666666\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\left(y \cdot y\right) \cdot \left(x \cdot \left(0.16666666666666666 + \left(x \cdot x\right) \cdot -0.027777777777777776\right)\right)\\
\end{array}
\end{array}
if y < 4.00000000000000021e236Initial program 100.0%
Taylor expanded in y around 0
distribute-rgt-inN/A
*-rgt-identityN/A
distribute-rgt-inN/A
*-commutativeN/A
associate-*r*N/A
distribute-rgt-outN/A
+-commutativeN/A
associate-*l*N/A
*-commutativeN/A
distribute-lft-inN/A
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
+-lowering-+.f64N/A
Simplified90.3%
Taylor expanded in x around 0
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6461.4%
Simplified61.4%
if 4.00000000000000021e236 < y Initial program 100.0%
Taylor expanded in y around 0
*-lft-identityN/A
associate-*r*N/A
distribute-rgt-inN/A
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64100.0%
Simplified100.0%
Taylor expanded in y around inf
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64100.0%
Simplified100.0%
Taylor expanded in x around 0
distribute-lft-inN/A
associate-*r*N/A
associate-*r*N/A
associate-*r*N/A
*-commutativeN/A
distribute-rgt-outN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64N/A
*-commutativeN/A
distribute-lft-outN/A
*-lowering-*.f64N/A
+-commutativeN/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f64100.0%
Simplified100.0%
Final simplification64.4%
(FPCore (x y) :precision binary64 (if (<= y 6.5e+14) (* x (+ 1.0 (* x (* x -0.16666666666666666)))) (* x (* y (* y (* (* y y) 0.008333333333333333))))))
double code(double x, double y) {
double tmp;
if (y <= 6.5e+14) {
tmp = x * (1.0 + (x * (x * -0.16666666666666666)));
} else {
tmp = x * (y * (y * ((y * y) * 0.008333333333333333)));
}
return tmp;
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8) :: tmp
if (y <= 6.5d+14) then
tmp = x * (1.0d0 + (x * (x * (-0.16666666666666666d0))))
else
tmp = x * (y * (y * ((y * y) * 0.008333333333333333d0)))
end if
code = tmp
end function
public static double code(double x, double y) {
double tmp;
if (y <= 6.5e+14) {
tmp = x * (1.0 + (x * (x * -0.16666666666666666)));
} else {
tmp = x * (y * (y * ((y * y) * 0.008333333333333333)));
}
return tmp;
}
def code(x, y): tmp = 0 if y <= 6.5e+14: tmp = x * (1.0 + (x * (x * -0.16666666666666666))) else: tmp = x * (y * (y * ((y * y) * 0.008333333333333333))) return tmp
function code(x, y) tmp = 0.0 if (y <= 6.5e+14) tmp = Float64(x * Float64(1.0 + Float64(x * Float64(x * -0.16666666666666666)))); else tmp = Float64(x * Float64(y * Float64(y * Float64(Float64(y * y) * 0.008333333333333333)))); end return tmp end
function tmp_2 = code(x, y) tmp = 0.0; if (y <= 6.5e+14) tmp = x * (1.0 + (x * (x * -0.16666666666666666))); else tmp = x * (y * (y * ((y * y) * 0.008333333333333333))); end tmp_2 = tmp; end
code[x_, y_] := If[LessEqual[y, 6.5e+14], N[(x * N[(1.0 + N[(x * N[(x * -0.16666666666666666), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x * N[(y * N[(y * N[(N[(y * y), $MachinePrecision] * 0.008333333333333333), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq 6.5 \cdot 10^{+14}:\\
\;\;\;\;x \cdot \left(1 + x \cdot \left(x \cdot -0.16666666666666666\right)\right)\\
\mathbf{else}:\\
\;\;\;\;x \cdot \left(y \cdot \left(y \cdot \left(\left(y \cdot y\right) \cdot 0.008333333333333333\right)\right)\right)\\
\end{array}
\end{array}
if y < 6.5e14Initial program 100.0%
Taylor expanded in y around 0
sin-lowering-sin.f6468.5%
Simplified68.5%
Taylor expanded in x around 0
*-lowering-*.f64N/A
+-lowering-+.f64N/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6446.2%
Simplified46.2%
Taylor expanded in x around 0
*-commutativeN/A
*-lowering-*.f6444.1%
Simplified44.1%
if 6.5e14 < y Initial program 100.0%
Taylor expanded in y around 0
distribute-rgt-inN/A
*-rgt-identityN/A
distribute-rgt-inN/A
*-commutativeN/A
associate-*r*N/A
distribute-rgt-outN/A
+-commutativeN/A
associate-*l*N/A
*-commutativeN/A
distribute-lft-inN/A
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
+-lowering-+.f64N/A
Simplified91.1%
Taylor expanded in y around inf
associate-*r*N/A
*-commutativeN/A
metadata-evalN/A
pow-sqrN/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
associate-*r*N/A
unpow2N/A
unpow3N/A
*-lowering-*.f64N/A
unpow3N/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6491.1%
Simplified91.1%
Taylor expanded in x around 0
Simplified71.0%
(FPCore (x y) :precision binary64 (if (<= x 1.45e+235) (* x (+ 1.0 (* (* y y) 0.16666666666666666))) (* x (+ 1.0 (* x (* x -0.16666666666666666))))))
double code(double x, double y) {
double tmp;
if (x <= 1.45e+235) {
tmp = x * (1.0 + ((y * y) * 0.16666666666666666));
} else {
tmp = x * (1.0 + (x * (x * -0.16666666666666666)));
}
return tmp;
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8) :: tmp
if (x <= 1.45d+235) then
tmp = x * (1.0d0 + ((y * y) * 0.16666666666666666d0))
else
tmp = x * (1.0d0 + (x * (x * (-0.16666666666666666d0))))
end if
code = tmp
end function
public static double code(double x, double y) {
double tmp;
if (x <= 1.45e+235) {
tmp = x * (1.0 + ((y * y) * 0.16666666666666666));
} else {
tmp = x * (1.0 + (x * (x * -0.16666666666666666)));
}
return tmp;
}
def code(x, y): tmp = 0 if x <= 1.45e+235: tmp = x * (1.0 + ((y * y) * 0.16666666666666666)) else: tmp = x * (1.0 + (x * (x * -0.16666666666666666))) return tmp
function code(x, y) tmp = 0.0 if (x <= 1.45e+235) tmp = Float64(x * Float64(1.0 + Float64(Float64(y * y) * 0.16666666666666666))); else tmp = Float64(x * Float64(1.0 + Float64(x * Float64(x * -0.16666666666666666)))); end return tmp end
function tmp_2 = code(x, y) tmp = 0.0; if (x <= 1.45e+235) tmp = x * (1.0 + ((y * y) * 0.16666666666666666)); else tmp = x * (1.0 + (x * (x * -0.16666666666666666))); end tmp_2 = tmp; end
code[x_, y_] := If[LessEqual[x, 1.45e+235], N[(x * N[(1.0 + N[(N[(y * y), $MachinePrecision] * 0.16666666666666666), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(x * N[(1.0 + N[(x * N[(x * -0.16666666666666666), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq 1.45 \cdot 10^{+235}:\\
\;\;\;\;x \cdot \left(1 + \left(y \cdot y\right) \cdot 0.16666666666666666\right)\\
\mathbf{else}:\\
\;\;\;\;x \cdot \left(1 + x \cdot \left(x \cdot -0.16666666666666666\right)\right)\\
\end{array}
\end{array}
if x < 1.45000000000000011e235Initial program 100.0%
Taylor expanded in y around 0
*-lft-identityN/A
associate-*r*N/A
distribute-rgt-inN/A
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f6477.3%
Simplified77.3%
Taylor expanded in x around 0
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6456.1%
Simplified56.1%
if 1.45000000000000011e235 < x Initial program 100.0%
Taylor expanded in y around 0
sin-lowering-sin.f6448.3%
Simplified48.3%
Taylor expanded in x around 0
*-lowering-*.f64N/A
+-lowering-+.f64N/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
sub-negN/A
metadata-evalN/A
+-commutativeN/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f647.6%
Simplified7.6%
Taylor expanded in x around 0
*-commutativeN/A
*-lowering-*.f6447.5%
Simplified47.5%
Final simplification55.6%
(FPCore (x y) :precision binary64 (if (<= y 140.0) x (* x (* (* y y) 0.16666666666666666))))
double code(double x, double y) {
double tmp;
if (y <= 140.0) {
tmp = x;
} else {
tmp = x * ((y * y) * 0.16666666666666666);
}
return tmp;
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8) :: tmp
if (y <= 140.0d0) then
tmp = x
else
tmp = x * ((y * y) * 0.16666666666666666d0)
end if
code = tmp
end function
public static double code(double x, double y) {
double tmp;
if (y <= 140.0) {
tmp = x;
} else {
tmp = x * ((y * y) * 0.16666666666666666);
}
return tmp;
}
def code(x, y): tmp = 0 if y <= 140.0: tmp = x else: tmp = x * ((y * y) * 0.16666666666666666) return tmp
function code(x, y) tmp = 0.0 if (y <= 140.0) tmp = x; else tmp = Float64(x * Float64(Float64(y * y) * 0.16666666666666666)); end return tmp end
function tmp_2 = code(x, y) tmp = 0.0; if (y <= 140.0) tmp = x; else tmp = x * ((y * y) * 0.16666666666666666); end tmp_2 = tmp; end
code[x_, y_] := If[LessEqual[y, 140.0], x, N[(x * N[(N[(y * y), $MachinePrecision] * 0.16666666666666666), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq 140:\\
\;\;\;\;x\\
\mathbf{else}:\\
\;\;\;\;x \cdot \left(\left(y \cdot y\right) \cdot 0.16666666666666666\right)\\
\end{array}
\end{array}
if y < 140Initial program 100.0%
Taylor expanded in y around 0
sin-lowering-sin.f6469.9%
Simplified69.9%
Taylor expanded in x around 0
Simplified40.5%
if 140 < y Initial program 100.0%
Taylor expanded in y around 0
*-lft-identityN/A
associate-*r*N/A
distribute-rgt-inN/A
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f6459.7%
Simplified59.7%
Taylor expanded in y around inf
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6459.7%
Simplified59.7%
Taylor expanded in x around 0
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6448.9%
Simplified48.9%
Final simplification42.7%
(FPCore (x y) :precision binary64 (* x (+ 1.0 (* (* y y) 0.16666666666666666))))
double code(double x, double y) {
return x * (1.0 + ((y * y) * 0.16666666666666666));
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
code = x * (1.0d0 + ((y * y) * 0.16666666666666666d0))
end function
public static double code(double x, double y) {
return x * (1.0 + ((y * y) * 0.16666666666666666));
}
def code(x, y): return x * (1.0 + ((y * y) * 0.16666666666666666))
function code(x, y) return Float64(x * Float64(1.0 + Float64(Float64(y * y) * 0.16666666666666666))) end
function tmp = code(x, y) tmp = x * (1.0 + ((y * y) * 0.16666666666666666)); end
code[x_, y_] := N[(x * N[(1.0 + N[(N[(y * y), $MachinePrecision] * 0.16666666666666666), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
x \cdot \left(1 + \left(y \cdot y\right) \cdot 0.16666666666666666\right)
\end{array}
Initial program 100.0%
Taylor expanded in y around 0
*-lft-identityN/A
associate-*r*N/A
distribute-rgt-inN/A
*-lowering-*.f64N/A
sin-lowering-sin.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f6478.2%
Simplified78.2%
Taylor expanded in x around 0
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6453.2%
Simplified53.2%
Final simplification53.2%
(FPCore (x y) :precision binary64 x)
double code(double x, double y) {
return x;
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
code = x
end function
public static double code(double x, double y) {
return x;
}
def code(x, y): return x
function code(x, y) return x end
function tmp = code(x, y) tmp = x; end
code[x_, y_] := x
\begin{array}{l}
\\
x
\end{array}
Initial program 100.0%
Taylor expanded in y around 0
sin-lowering-sin.f6451.8%
Simplified51.8%
Taylor expanded in x around 0
Simplified30.2%
herbie shell --seed 2024192
(FPCore (x y)
:name "Linear.Quaternion:$ccos from linear-1.19.1.3"
:precision binary64
(* (sin x) (/ (sinh y) y)))