
(FPCore (x y) :precision binary64 (* (cos x) (/ (sinh y) y)))
double code(double x, double y) {
return cos(x) * (sinh(y) / y);
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
code = cos(x) * (sinh(y) / y)
end function
public static double code(double x, double y) {
return Math.cos(x) * (Math.sinh(y) / y);
}
def code(x, y): return math.cos(x) * (math.sinh(y) / y)
function code(x, y) return Float64(cos(x) * Float64(sinh(y) / y)) end
function tmp = code(x, y) tmp = cos(x) * (sinh(y) / y); end
code[x_, y_] := N[(N[Cos[x], $MachinePrecision] * N[(N[Sinh[y], $MachinePrecision] / y), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\cos x \cdot \frac{\sinh y}{y}
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 20 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (x y) :precision binary64 (* (cos x) (/ (sinh y) y)))
double code(double x, double y) {
return cos(x) * (sinh(y) / y);
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
code = cos(x) * (sinh(y) / y)
end function
public static double code(double x, double y) {
return Math.cos(x) * (Math.sinh(y) / y);
}
def code(x, y): return math.cos(x) * (math.sinh(y) / y)
function code(x, y) return Float64(cos(x) * Float64(sinh(y) / y)) end
function tmp = code(x, y) tmp = cos(x) * (sinh(y) / y); end
code[x_, y_] := N[(N[Cos[x], $MachinePrecision] * N[(N[Sinh[y], $MachinePrecision] / y), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\cos x \cdot \frac{\sinh y}{y}
\end{array}
(FPCore (x y) :precision binary64 (* (cos x) (/ (sinh y) y)))
double code(double x, double y) {
return cos(x) * (sinh(y) / y);
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
code = cos(x) * (sinh(y) / y)
end function
public static double code(double x, double y) {
return Math.cos(x) * (Math.sinh(y) / y);
}
def code(x, y): return math.cos(x) * (math.sinh(y) / y)
function code(x, y) return Float64(cos(x) * Float64(sinh(y) / y)) end
function tmp = code(x, y) tmp = cos(x) * (sinh(y) / y); end
code[x_, y_] := N[(N[Cos[x], $MachinePrecision] * N[(N[Sinh[y], $MachinePrecision] / y), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\cos x \cdot \frac{\sinh y}{y}
\end{array}
Initial program 100.0%
(FPCore (x y)
:precision binary64
(let* ((t_0 (+ 0.16666666666666666 (* (* y y) 0.008333333333333333))))
(if (<= y 0.215)
(* (cos x) (+ 1.0 (* (* y y) t_0)))
(if (<= y 1.15e+62)
(/ (sinh y) y)
(/ (* (cos x) (* y (+ 1.0 (* y (* y t_0))))) y)))))
double code(double x, double y) {
double t_0 = 0.16666666666666666 + ((y * y) * 0.008333333333333333);
double tmp;
if (y <= 0.215) {
tmp = cos(x) * (1.0 + ((y * y) * t_0));
} else if (y <= 1.15e+62) {
tmp = sinh(y) / y;
} else {
tmp = (cos(x) * (y * (1.0 + (y * (y * t_0))))) / y;
}
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 = 0.16666666666666666d0 + ((y * y) * 0.008333333333333333d0)
if (y <= 0.215d0) then
tmp = cos(x) * (1.0d0 + ((y * y) * t_0))
else if (y <= 1.15d+62) then
tmp = sinh(y) / y
else
tmp = (cos(x) * (y * (1.0d0 + (y * (y * t_0))))) / y
end if
code = tmp
end function
public static double code(double x, double y) {
double t_0 = 0.16666666666666666 + ((y * y) * 0.008333333333333333);
double tmp;
if (y <= 0.215) {
tmp = Math.cos(x) * (1.0 + ((y * y) * t_0));
} else if (y <= 1.15e+62) {
tmp = Math.sinh(y) / y;
} else {
tmp = (Math.cos(x) * (y * (1.0 + (y * (y * t_0))))) / y;
}
return tmp;
}
def code(x, y): t_0 = 0.16666666666666666 + ((y * y) * 0.008333333333333333) tmp = 0 if y <= 0.215: tmp = math.cos(x) * (1.0 + ((y * y) * t_0)) elif y <= 1.15e+62: tmp = math.sinh(y) / y else: tmp = (math.cos(x) * (y * (1.0 + (y * (y * t_0))))) / y return tmp
function code(x, y) t_0 = Float64(0.16666666666666666 + Float64(Float64(y * y) * 0.008333333333333333)) tmp = 0.0 if (y <= 0.215) tmp = Float64(cos(x) * Float64(1.0 + Float64(Float64(y * y) * t_0))); elseif (y <= 1.15e+62) tmp = Float64(sinh(y) / y); else tmp = Float64(Float64(cos(x) * Float64(y * Float64(1.0 + Float64(y * Float64(y * t_0))))) / y); end return tmp end
function tmp_2 = code(x, y) t_0 = 0.16666666666666666 + ((y * y) * 0.008333333333333333); tmp = 0.0; if (y <= 0.215) tmp = cos(x) * (1.0 + ((y * y) * t_0)); elseif (y <= 1.15e+62) tmp = sinh(y) / y; else tmp = (cos(x) * (y * (1.0 + (y * (y * t_0))))) / y; end tmp_2 = tmp; end
code[x_, y_] := Block[{t$95$0 = N[(0.16666666666666666 + N[(N[(y * y), $MachinePrecision] * 0.008333333333333333), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, 0.215], N[(N[Cos[x], $MachinePrecision] * N[(1.0 + N[(N[(y * y), $MachinePrecision] * t$95$0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[y, 1.15e+62], N[(N[Sinh[y], $MachinePrecision] / y), $MachinePrecision], N[(N[(N[Cos[x], $MachinePrecision] * N[(y * N[(1.0 + N[(y * N[(y * t$95$0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := 0.16666666666666666 + \left(y \cdot y\right) \cdot 0.008333333333333333\\
\mathbf{if}\;y \leq 0.215:\\
\;\;\;\;\cos x \cdot \left(1 + \left(y \cdot y\right) \cdot t\_0\right)\\
\mathbf{elif}\;y \leq 1.15 \cdot 10^{+62}:\\
\;\;\;\;\frac{\sinh y}{y}\\
\mathbf{else}:\\
\;\;\;\;\frac{\cos x \cdot \left(y \cdot \left(1 + y \cdot \left(y \cdot t\_0\right)\right)\right)}{y}\\
\end{array}
\end{array}
if y < 0.214999999999999997Initial 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
cos-lowering-cos.f64N/A
+-lowering-+.f64N/A
Simplified91.0%
if 0.214999999999999997 < y < 1.14999999999999992e62Initial program 100.0%
Taylor expanded in x around 0
Simplified87.5%
*-lft-identityN/A
/-lowering-/.f64N/A
sinh-lowering-sinh.f6487.5%
Applied egg-rr87.5%
if 1.14999999999999992e62 < y Initial program 100.0%
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
cos-lowering-cos.f64N/A
sinh-lowering-sinh.f64100.0%
Applied egg-rr100.0%
Taylor expanded in y around 0
*-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-*.f64100.0%
Simplified100.0%
(FPCore (x y)
:precision binary64
(let* ((t_0
(*
(cos x)
(+
1.0
(*
(* y y)
(+ 0.16666666666666666 (* (* y y) 0.008333333333333333)))))))
(if (<= y 0.4) t_0 (if (<= y 3.8e+77) (/ (sinh y) y) t_0))))
double code(double x, double y) {
double t_0 = cos(x) * (1.0 + ((y * y) * (0.16666666666666666 + ((y * y) * 0.008333333333333333))));
double tmp;
if (y <= 0.4) {
tmp = t_0;
} else if (y <= 3.8e+77) {
tmp = sinh(y) / y;
} 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 = cos(x) * (1.0d0 + ((y * y) * (0.16666666666666666d0 + ((y * y) * 0.008333333333333333d0))))
if (y <= 0.4d0) then
tmp = t_0
else if (y <= 3.8d+77) then
tmp = sinh(y) / y
else
tmp = t_0
end if
code = tmp
end function
public static double code(double x, double y) {
double t_0 = Math.cos(x) * (1.0 + ((y * y) * (0.16666666666666666 + ((y * y) * 0.008333333333333333))));
double tmp;
if (y <= 0.4) {
tmp = t_0;
} else if (y <= 3.8e+77) {
tmp = Math.sinh(y) / y;
} else {
tmp = t_0;
}
return tmp;
}
def code(x, y): t_0 = math.cos(x) * (1.0 + ((y * y) * (0.16666666666666666 + ((y * y) * 0.008333333333333333)))) tmp = 0 if y <= 0.4: tmp = t_0 elif y <= 3.8e+77: tmp = math.sinh(y) / y else: tmp = t_0 return tmp
function code(x, y) t_0 = Float64(cos(x) * Float64(1.0 + Float64(Float64(y * y) * Float64(0.16666666666666666 + Float64(Float64(y * y) * 0.008333333333333333))))) tmp = 0.0 if (y <= 0.4) tmp = t_0; elseif (y <= 3.8e+77) tmp = Float64(sinh(y) / y); else tmp = t_0; end return tmp end
function tmp_2 = code(x, y) t_0 = cos(x) * (1.0 + ((y * y) * (0.16666666666666666 + ((y * y) * 0.008333333333333333)))); tmp = 0.0; if (y <= 0.4) tmp = t_0; elseif (y <= 3.8e+77) tmp = sinh(y) / y; else tmp = t_0; end tmp_2 = tmp; end
code[x_, y_] := Block[{t$95$0 = N[(N[Cos[x], $MachinePrecision] * N[(1.0 + N[(N[(y * y), $MachinePrecision] * N[(0.16666666666666666 + N[(N[(y * y), $MachinePrecision] * 0.008333333333333333), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, 0.4], t$95$0, If[LessEqual[y, 3.8e+77], N[(N[Sinh[y], $MachinePrecision] / y), $MachinePrecision], t$95$0]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \cos x \cdot \left(1 + \left(y \cdot y\right) \cdot \left(0.16666666666666666 + \left(y \cdot y\right) \cdot 0.008333333333333333\right)\right)\\
\mathbf{if}\;y \leq 0.4:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;y \leq 3.8 \cdot 10^{+77}:\\
\;\;\;\;\frac{\sinh y}{y}\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
if y < 0.40000000000000002 or 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
cos-lowering-cos.f64N/A
+-lowering-+.f64N/A
Simplified93.0%
if 0.40000000000000002 < y < 3.8000000000000001e77Initial program 100.0%
Taylor expanded in x around 0
Simplified84.6%
*-lft-identityN/A
/-lowering-/.f64N/A
sinh-lowering-sinh.f6484.6%
Applied egg-rr84.6%
(FPCore (x y)
:precision binary64
(let* ((t_0 (+ 1.0 (* (* y y) 0.16666666666666666))))
(if (<= y 0.0135)
(* (cos x) t_0)
(if (<= y 1.3e+103) (/ (sinh y) y) (/ (* t_0 (* (cos x) y)) y)))))
double code(double x, double y) {
double t_0 = 1.0 + ((y * y) * 0.16666666666666666);
double tmp;
if (y <= 0.0135) {
tmp = cos(x) * t_0;
} else if (y <= 1.3e+103) {
tmp = sinh(y) / y;
} else {
tmp = (t_0 * (cos(x) * y)) / y;
}
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 = 1.0d0 + ((y * y) * 0.16666666666666666d0)
if (y <= 0.0135d0) then
tmp = cos(x) * t_0
else if (y <= 1.3d+103) then
tmp = sinh(y) / y
else
tmp = (t_0 * (cos(x) * y)) / y
end if
code = tmp
end function
public static double code(double x, double y) {
double t_0 = 1.0 + ((y * y) * 0.16666666666666666);
double tmp;
if (y <= 0.0135) {
tmp = Math.cos(x) * t_0;
} else if (y <= 1.3e+103) {
tmp = Math.sinh(y) / y;
} else {
tmp = (t_0 * (Math.cos(x) * y)) / y;
}
return tmp;
}
def code(x, y): t_0 = 1.0 + ((y * y) * 0.16666666666666666) tmp = 0 if y <= 0.0135: tmp = math.cos(x) * t_0 elif y <= 1.3e+103: tmp = math.sinh(y) / y else: tmp = (t_0 * (math.cos(x) * y)) / y return tmp
function code(x, y) t_0 = Float64(1.0 + Float64(Float64(y * y) * 0.16666666666666666)) tmp = 0.0 if (y <= 0.0135) tmp = Float64(cos(x) * t_0); elseif (y <= 1.3e+103) tmp = Float64(sinh(y) / y); else tmp = Float64(Float64(t_0 * Float64(cos(x) * y)) / y); end return tmp end
function tmp_2 = code(x, y) t_0 = 1.0 + ((y * y) * 0.16666666666666666); tmp = 0.0; if (y <= 0.0135) tmp = cos(x) * t_0; elseif (y <= 1.3e+103) tmp = sinh(y) / y; else tmp = (t_0 * (cos(x) * y)) / y; end tmp_2 = tmp; end
code[x_, y_] := Block[{t$95$0 = N[(1.0 + N[(N[(y * y), $MachinePrecision] * 0.16666666666666666), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, 0.0135], N[(N[Cos[x], $MachinePrecision] * t$95$0), $MachinePrecision], If[LessEqual[y, 1.3e+103], N[(N[Sinh[y], $MachinePrecision] / y), $MachinePrecision], N[(N[(t$95$0 * N[(N[Cos[x], $MachinePrecision] * y), $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := 1 + \left(y \cdot y\right) \cdot 0.16666666666666666\\
\mathbf{if}\;y \leq 0.0135:\\
\;\;\;\;\cos x \cdot t\_0\\
\mathbf{elif}\;y \leq 1.3 \cdot 10^{+103}:\\
\;\;\;\;\frac{\sinh y}{y}\\
\mathbf{else}:\\
\;\;\;\;\frac{t\_0 \cdot \left(\cos x \cdot y\right)}{y}\\
\end{array}
\end{array}
if y < 0.0134999999999999998Initial program 100.0%
Taylor expanded in y around 0
associate-*r*N/A
*-lft-identityN/A
distribute-rgt-inN/A
*-lowering-*.f64N/A
cos-lowering-cos.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6481.8%
Simplified81.8%
if 0.0134999999999999998 < y < 1.3000000000000001e103Initial program 100.0%
Taylor expanded in x around 0
Simplified83.3%
*-lft-identityN/A
/-lowering-/.f64N/A
sinh-lowering-sinh.f6483.3%
Applied egg-rr83.3%
if 1.3000000000000001e103 < y Initial program 100.0%
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
cos-lowering-cos.f64N/A
sinh-lowering-sinh.f64100.0%
Applied egg-rr100.0%
Taylor expanded in y around 0
+-commutativeN/A
*-commutativeN/A
associate-*r*N/A
*-commutativeN/A
distribute-rgt-inN/A
associate-*r*N/A
*-commutativeN/A
associate-*l*N/A
*-commutativeN/A
*-commutativeN/A
distribute-lft1-inN/A
+-commutativeN/A
*-lowering-*.f64N/A
Simplified100.0%
Final simplification85.3%
(FPCore (x y) :precision binary64 (let* ((t_0 (* (cos x) (+ 1.0 (* (* y y) 0.16666666666666666))))) (if (<= y 0.03) t_0 (if (<= y 1.1e+154) (/ (sinh y) y) t_0))))
double code(double x, double y) {
double t_0 = cos(x) * (1.0 + ((y * y) * 0.16666666666666666));
double tmp;
if (y <= 0.03) {
tmp = t_0;
} else if (y <= 1.1e+154) {
tmp = sinh(y) / y;
} 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 = cos(x) * (1.0d0 + ((y * y) * 0.16666666666666666d0))
if (y <= 0.03d0) then
tmp = t_0
else if (y <= 1.1d+154) then
tmp = sinh(y) / y
else
tmp = t_0
end if
code = tmp
end function
public static double code(double x, double y) {
double t_0 = Math.cos(x) * (1.0 + ((y * y) * 0.16666666666666666));
double tmp;
if (y <= 0.03) {
tmp = t_0;
} else if (y <= 1.1e+154) {
tmp = Math.sinh(y) / y;
} else {
tmp = t_0;
}
return tmp;
}
def code(x, y): t_0 = math.cos(x) * (1.0 + ((y * y) * 0.16666666666666666)) tmp = 0 if y <= 0.03: tmp = t_0 elif y <= 1.1e+154: tmp = math.sinh(y) / y else: tmp = t_0 return tmp
function code(x, y) t_0 = Float64(cos(x) * Float64(1.0 + Float64(Float64(y * y) * 0.16666666666666666))) tmp = 0.0 if (y <= 0.03) tmp = t_0; elseif (y <= 1.1e+154) tmp = Float64(sinh(y) / y); else tmp = t_0; end return tmp end
function tmp_2 = code(x, y) t_0 = cos(x) * (1.0 + ((y * y) * 0.16666666666666666)); tmp = 0.0; if (y <= 0.03) tmp = t_0; elseif (y <= 1.1e+154) tmp = sinh(y) / y; else tmp = t_0; end tmp_2 = tmp; end
code[x_, y_] := Block[{t$95$0 = N[(N[Cos[x], $MachinePrecision] * N[(1.0 + N[(N[(y * y), $MachinePrecision] * 0.16666666666666666), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, 0.03], t$95$0, If[LessEqual[y, 1.1e+154], N[(N[Sinh[y], $MachinePrecision] / y), $MachinePrecision], t$95$0]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \cos x \cdot \left(1 + \left(y \cdot y\right) \cdot 0.16666666666666666\right)\\
\mathbf{if}\;y \leq 0.03:\\
\;\;\;\;t\_0\\
\mathbf{elif}\;y \leq 1.1 \cdot 10^{+154}:\\
\;\;\;\;\frac{\sinh y}{y}\\
\mathbf{else}:\\
\;\;\;\;t\_0\\
\end{array}
\end{array}
if y < 0.029999999999999999 or 1.1000000000000001e154 < y Initial program 100.0%
Taylor expanded in y around 0
associate-*r*N/A
*-lft-identityN/A
distribute-rgt-inN/A
*-lowering-*.f64N/A
cos-lowering-cos.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6484.2%
Simplified84.2%
if 0.029999999999999999 < y < 1.1000000000000001e154Initial program 100.0%
Taylor expanded in x around 0
Simplified78.1%
*-lft-identityN/A
/-lowering-/.f64N/A
sinh-lowering-sinh.f6478.1%
Applied egg-rr78.1%
Final simplification83.4%
(FPCore (x y) :precision binary64 (if (<= y 0.0076) (cos x) (/ (sinh y) y)))
double code(double x, double y) {
double tmp;
if (y <= 0.0076) {
tmp = cos(x);
} else {
tmp = sinh(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.0076d0) then
tmp = cos(x)
else
tmp = sinh(y) / y
end if
code = tmp
end function
public static double code(double x, double y) {
double tmp;
if (y <= 0.0076) {
tmp = Math.cos(x);
} else {
tmp = Math.sinh(y) / y;
}
return tmp;
}
def code(x, y): tmp = 0 if y <= 0.0076: tmp = math.cos(x) else: tmp = math.sinh(y) / y return tmp
function code(x, y) tmp = 0.0 if (y <= 0.0076) tmp = cos(x); else tmp = Float64(sinh(y) / y); end return tmp end
function tmp_2 = code(x, y) tmp = 0.0; if (y <= 0.0076) tmp = cos(x); else tmp = sinh(y) / y; end tmp_2 = tmp; end
code[x_, y_] := If[LessEqual[y, 0.0076], N[Cos[x], $MachinePrecision], N[(N[Sinh[y], $MachinePrecision] / y), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq 0.0076:\\
\;\;\;\;\cos x\\
\mathbf{else}:\\
\;\;\;\;\frac{\sinh y}{y}\\
\end{array}
\end{array}
if y < 0.00759999999999999998Initial program 100.0%
Taylor expanded in y around 0
cos-lowering-cos.f6464.4%
Simplified64.4%
if 0.00759999999999999998 < y Initial program 100.0%
Taylor expanded in x around 0
Simplified74.2%
*-lft-identityN/A
/-lowering-/.f64N/A
sinh-lowering-sinh.f6474.2%
Applied egg-rr74.2%
(FPCore (x y)
:precision binary64
(let* ((t_0
(*
(* y y)
(+
0.16666666666666666
(*
(* y y)
(+ 0.008333333333333333 (* (* y y) 0.0001984126984126984)))))))
(if (<= y 5.8)
(cos x)
(if (<= y 1.35e+106)
(/ 1.0 (/ (+ 1.0 (* y (* y -0.16666666666666666))) (- 1.0 (* t_0 t_0))))
(*
(+
1.0
(* y (* y (+ 0.16666666666666666 (* (* y y) 0.008333333333333333)))))
(+ 1.0 (* x (* x -0.5))))))))
double code(double x, double y) {
double t_0 = (y * y) * (0.16666666666666666 + ((y * y) * (0.008333333333333333 + ((y * y) * 0.0001984126984126984))));
double tmp;
if (y <= 5.8) {
tmp = cos(x);
} else if (y <= 1.35e+106) {
tmp = 1.0 / ((1.0 + (y * (y * -0.16666666666666666))) / (1.0 - (t_0 * t_0)));
} else {
tmp = (1.0 + (y * (y * (0.16666666666666666 + ((y * y) * 0.008333333333333333))))) * (1.0 + (x * (x * -0.5)));
}
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 = (y * y) * (0.16666666666666666d0 + ((y * y) * (0.008333333333333333d0 + ((y * y) * 0.0001984126984126984d0))))
if (y <= 5.8d0) then
tmp = cos(x)
else if (y <= 1.35d+106) then
tmp = 1.0d0 / ((1.0d0 + (y * (y * (-0.16666666666666666d0)))) / (1.0d0 - (t_0 * t_0)))
else
tmp = (1.0d0 + (y * (y * (0.16666666666666666d0 + ((y * y) * 0.008333333333333333d0))))) * (1.0d0 + (x * (x * (-0.5d0))))
end if
code = tmp
end function
public static double code(double x, double y) {
double t_0 = (y * y) * (0.16666666666666666 + ((y * y) * (0.008333333333333333 + ((y * y) * 0.0001984126984126984))));
double tmp;
if (y <= 5.8) {
tmp = Math.cos(x);
} else if (y <= 1.35e+106) {
tmp = 1.0 / ((1.0 + (y * (y * -0.16666666666666666))) / (1.0 - (t_0 * t_0)));
} else {
tmp = (1.0 + (y * (y * (0.16666666666666666 + ((y * y) * 0.008333333333333333))))) * (1.0 + (x * (x * -0.5)));
}
return tmp;
}
def code(x, y): t_0 = (y * y) * (0.16666666666666666 + ((y * y) * (0.008333333333333333 + ((y * y) * 0.0001984126984126984)))) tmp = 0 if y <= 5.8: tmp = math.cos(x) elif y <= 1.35e+106: tmp = 1.0 / ((1.0 + (y * (y * -0.16666666666666666))) / (1.0 - (t_0 * t_0))) else: tmp = (1.0 + (y * (y * (0.16666666666666666 + ((y * y) * 0.008333333333333333))))) * (1.0 + (x * (x * -0.5))) return tmp
function code(x, y) t_0 = Float64(Float64(y * y) * Float64(0.16666666666666666 + Float64(Float64(y * y) * Float64(0.008333333333333333 + Float64(Float64(y * y) * 0.0001984126984126984))))) tmp = 0.0 if (y <= 5.8) tmp = cos(x); elseif (y <= 1.35e+106) tmp = Float64(1.0 / Float64(Float64(1.0 + Float64(y * Float64(y * -0.16666666666666666))) / Float64(1.0 - Float64(t_0 * t_0)))); else tmp = Float64(Float64(1.0 + Float64(y * Float64(y * Float64(0.16666666666666666 + Float64(Float64(y * y) * 0.008333333333333333))))) * Float64(1.0 + Float64(x * Float64(x * -0.5)))); end return tmp end
function tmp_2 = code(x, y) t_0 = (y * y) * (0.16666666666666666 + ((y * y) * (0.008333333333333333 + ((y * y) * 0.0001984126984126984)))); tmp = 0.0; if (y <= 5.8) tmp = cos(x); elseif (y <= 1.35e+106) tmp = 1.0 / ((1.0 + (y * (y * -0.16666666666666666))) / (1.0 - (t_0 * t_0))); else tmp = (1.0 + (y * (y * (0.16666666666666666 + ((y * y) * 0.008333333333333333))))) * (1.0 + (x * (x * -0.5))); end tmp_2 = tmp; end
code[x_, y_] := Block[{t$95$0 = N[(N[(y * y), $MachinePrecision] * N[(0.16666666666666666 + N[(N[(y * y), $MachinePrecision] * N[(0.008333333333333333 + N[(N[(y * y), $MachinePrecision] * 0.0001984126984126984), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, 5.8], N[Cos[x], $MachinePrecision], If[LessEqual[y, 1.35e+106], N[(1.0 / N[(N[(1.0 + N[(y * N[(y * -0.16666666666666666), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(1.0 - N[(t$95$0 * t$95$0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(1.0 + N[(y * N[(y * N[(0.16666666666666666 + N[(N[(y * y), $MachinePrecision] * 0.008333333333333333), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(1.0 + N[(x * N[(x * -0.5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \left(y \cdot y\right) \cdot \left(0.16666666666666666 + \left(y \cdot y\right) \cdot \left(0.008333333333333333 + \left(y \cdot y\right) \cdot 0.0001984126984126984\right)\right)\\
\mathbf{if}\;y \leq 5.8:\\
\;\;\;\;\cos x\\
\mathbf{elif}\;y \leq 1.35 \cdot 10^{+106}:\\
\;\;\;\;\frac{1}{\frac{1 + y \cdot \left(y \cdot -0.16666666666666666\right)}{1 - t\_0 \cdot t\_0}}\\
\mathbf{else}:\\
\;\;\;\;\left(1 + y \cdot \left(y \cdot \left(0.16666666666666666 + \left(y \cdot y\right) \cdot 0.008333333333333333\right)\right)\right) \cdot \left(1 + x \cdot \left(x \cdot -0.5\right)\right)\\
\end{array}
\end{array}
if y < 5.79999999999999982Initial program 100.0%
Taylor expanded in y around 0
cos-lowering-cos.f6464.4%
Simplified64.4%
if 5.79999999999999982 < y < 1.35000000000000003e106Initial program 100.0%
Taylor expanded in x around 0
Simplified85.0%
Taylor expanded in y around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/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-*.f6457.1%
Simplified57.1%
flip-+N/A
clear-numN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
Applied egg-rr6.9%
Taylor expanded in y around 0
+-lowering-+.f64N/A
*-commutativeN/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f6461.8%
Simplified61.8%
if 1.35000000000000003e106 < y Initial program 100.0%
Taylor expanded in x around 0
+-lowering-+.f64N/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f6482.6%
Simplified82.6%
Taylor expanded in y around 0
+-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-*.f6482.6%
Simplified82.6%
Final simplification67.5%
(FPCore (x y)
:precision binary64
(let* ((t_0
(*
(* y y)
(+
0.16666666666666666
(*
(* y y)
(+ 0.008333333333333333 (* (* y y) 0.0001984126984126984)))))))
(if (<= y 7e+106)
(/ 1.0 (/ (+ 1.0 (* y (* y -0.16666666666666666))) (- 1.0 (* t_0 t_0))))
(*
(+
1.0
(* y (* y (+ 0.16666666666666666 (* (* y y) 0.008333333333333333)))))
(+ 1.0 (* x (* x -0.5)))))))
double code(double x, double y) {
double t_0 = (y * y) * (0.16666666666666666 + ((y * y) * (0.008333333333333333 + ((y * y) * 0.0001984126984126984))));
double tmp;
if (y <= 7e+106) {
tmp = 1.0 / ((1.0 + (y * (y * -0.16666666666666666))) / (1.0 - (t_0 * t_0)));
} else {
tmp = (1.0 + (y * (y * (0.16666666666666666 + ((y * y) * 0.008333333333333333))))) * (1.0 + (x * (x * -0.5)));
}
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 = (y * y) * (0.16666666666666666d0 + ((y * y) * (0.008333333333333333d0 + ((y * y) * 0.0001984126984126984d0))))
if (y <= 7d+106) then
tmp = 1.0d0 / ((1.0d0 + (y * (y * (-0.16666666666666666d0)))) / (1.0d0 - (t_0 * t_0)))
else
tmp = (1.0d0 + (y * (y * (0.16666666666666666d0 + ((y * y) * 0.008333333333333333d0))))) * (1.0d0 + (x * (x * (-0.5d0))))
end if
code = tmp
end function
public static double code(double x, double y) {
double t_0 = (y * y) * (0.16666666666666666 + ((y * y) * (0.008333333333333333 + ((y * y) * 0.0001984126984126984))));
double tmp;
if (y <= 7e+106) {
tmp = 1.0 / ((1.0 + (y * (y * -0.16666666666666666))) / (1.0 - (t_0 * t_0)));
} else {
tmp = (1.0 + (y * (y * (0.16666666666666666 + ((y * y) * 0.008333333333333333))))) * (1.0 + (x * (x * -0.5)));
}
return tmp;
}
def code(x, y): t_0 = (y * y) * (0.16666666666666666 + ((y * y) * (0.008333333333333333 + ((y * y) * 0.0001984126984126984)))) tmp = 0 if y <= 7e+106: tmp = 1.0 / ((1.0 + (y * (y * -0.16666666666666666))) / (1.0 - (t_0 * t_0))) else: tmp = (1.0 + (y * (y * (0.16666666666666666 + ((y * y) * 0.008333333333333333))))) * (1.0 + (x * (x * -0.5))) return tmp
function code(x, y) t_0 = Float64(Float64(y * y) * Float64(0.16666666666666666 + Float64(Float64(y * y) * Float64(0.008333333333333333 + Float64(Float64(y * y) * 0.0001984126984126984))))) tmp = 0.0 if (y <= 7e+106) tmp = Float64(1.0 / Float64(Float64(1.0 + Float64(y * Float64(y * -0.16666666666666666))) / Float64(1.0 - Float64(t_0 * t_0)))); else tmp = Float64(Float64(1.0 + Float64(y * Float64(y * Float64(0.16666666666666666 + Float64(Float64(y * y) * 0.008333333333333333))))) * Float64(1.0 + Float64(x * Float64(x * -0.5)))); end return tmp end
function tmp_2 = code(x, y) t_0 = (y * y) * (0.16666666666666666 + ((y * y) * (0.008333333333333333 + ((y * y) * 0.0001984126984126984)))); tmp = 0.0; if (y <= 7e+106) tmp = 1.0 / ((1.0 + (y * (y * -0.16666666666666666))) / (1.0 - (t_0 * t_0))); else tmp = (1.0 + (y * (y * (0.16666666666666666 + ((y * y) * 0.008333333333333333))))) * (1.0 + (x * (x * -0.5))); end tmp_2 = tmp; end
code[x_, y_] := Block[{t$95$0 = N[(N[(y * y), $MachinePrecision] * N[(0.16666666666666666 + N[(N[(y * y), $MachinePrecision] * N[(0.008333333333333333 + N[(N[(y * y), $MachinePrecision] * 0.0001984126984126984), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[y, 7e+106], N[(1.0 / N[(N[(1.0 + N[(y * N[(y * -0.16666666666666666), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(1.0 - N[(t$95$0 * t$95$0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(1.0 + N[(y * N[(y * N[(0.16666666666666666 + N[(N[(y * y), $MachinePrecision] * 0.008333333333333333), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(1.0 + N[(x * N[(x * -0.5), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \left(y \cdot y\right) \cdot \left(0.16666666666666666 + \left(y \cdot y\right) \cdot \left(0.008333333333333333 + \left(y \cdot y\right) \cdot 0.0001984126984126984\right)\right)\\
\mathbf{if}\;y \leq 7 \cdot 10^{+106}:\\
\;\;\;\;\frac{1}{\frac{1 + y \cdot \left(y \cdot -0.16666666666666666\right)}{1 - t\_0 \cdot t\_0}}\\
\mathbf{else}:\\
\;\;\;\;\left(1 + y \cdot \left(y \cdot \left(0.16666666666666666 + \left(y \cdot y\right) \cdot 0.008333333333333333\right)\right)\right) \cdot \left(1 + x \cdot \left(x \cdot -0.5\right)\right)\\
\end{array}
\end{array}
if y < 6.99999999999999962e106Initial program 100.0%
Taylor expanded in x around 0
Simplified71.2%
Taylor expanded in y around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/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-*.f6463.6%
Simplified63.6%
flip-+N/A
clear-numN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
Applied egg-rr40.5%
Taylor expanded in y around 0
+-lowering-+.f64N/A
*-commutativeN/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f6453.8%
Simplified53.8%
if 6.99999999999999962e106 < y Initial program 100.0%
Taylor expanded in x around 0
+-lowering-+.f64N/A
unpow2N/A
associate-*r*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f6482.6%
Simplified82.6%
Taylor expanded in y around 0
+-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-*.f6482.6%
Simplified82.6%
Final simplification59.0%
(FPCore (x y)
:precision binary64
(let* ((t_0 (+ 0.008333333333333333 (* (* y y) 0.0001984126984126984)))
(t_1 (* (* y y) (+ 0.16666666666666666 (* (* y y) t_0)))))
(if (<= x 8.5e+258)
(/ (* y (+ 1.0 (* (* y y) (+ 0.16666666666666666 (* y (* y t_0)))))) y)
(/ 1.0 (/ 1.0 (- 1.0 (* t_1 t_1)))))))
double code(double x, double y) {
double t_0 = 0.008333333333333333 + ((y * y) * 0.0001984126984126984);
double t_1 = (y * y) * (0.16666666666666666 + ((y * y) * t_0));
double tmp;
if (x <= 8.5e+258) {
tmp = (y * (1.0 + ((y * y) * (0.16666666666666666 + (y * (y * t_0)))))) / y;
} else {
tmp = 1.0 / (1.0 / (1.0 - (t_1 * t_1)));
}
return tmp;
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8) :: t_0
real(8) :: t_1
real(8) :: tmp
t_0 = 0.008333333333333333d0 + ((y * y) * 0.0001984126984126984d0)
t_1 = (y * y) * (0.16666666666666666d0 + ((y * y) * t_0))
if (x <= 8.5d+258) then
tmp = (y * (1.0d0 + ((y * y) * (0.16666666666666666d0 + (y * (y * t_0)))))) / y
else
tmp = 1.0d0 / (1.0d0 / (1.0d0 - (t_1 * t_1)))
end if
code = tmp
end function
public static double code(double x, double y) {
double t_0 = 0.008333333333333333 + ((y * y) * 0.0001984126984126984);
double t_1 = (y * y) * (0.16666666666666666 + ((y * y) * t_0));
double tmp;
if (x <= 8.5e+258) {
tmp = (y * (1.0 + ((y * y) * (0.16666666666666666 + (y * (y * t_0)))))) / y;
} else {
tmp = 1.0 / (1.0 / (1.0 - (t_1 * t_1)));
}
return tmp;
}
def code(x, y): t_0 = 0.008333333333333333 + ((y * y) * 0.0001984126984126984) t_1 = (y * y) * (0.16666666666666666 + ((y * y) * t_0)) tmp = 0 if x <= 8.5e+258: tmp = (y * (1.0 + ((y * y) * (0.16666666666666666 + (y * (y * t_0)))))) / y else: tmp = 1.0 / (1.0 / (1.0 - (t_1 * t_1))) return tmp
function code(x, y) t_0 = Float64(0.008333333333333333 + Float64(Float64(y * y) * 0.0001984126984126984)) t_1 = Float64(Float64(y * y) * Float64(0.16666666666666666 + Float64(Float64(y * y) * t_0))) tmp = 0.0 if (x <= 8.5e+258) tmp = Float64(Float64(y * Float64(1.0 + Float64(Float64(y * y) * Float64(0.16666666666666666 + Float64(y * Float64(y * t_0)))))) / y); else tmp = Float64(1.0 / Float64(1.0 / Float64(1.0 - Float64(t_1 * t_1)))); end return tmp end
function tmp_2 = code(x, y) t_0 = 0.008333333333333333 + ((y * y) * 0.0001984126984126984); t_1 = (y * y) * (0.16666666666666666 + ((y * y) * t_0)); tmp = 0.0; if (x <= 8.5e+258) tmp = (y * (1.0 + ((y * y) * (0.16666666666666666 + (y * (y * t_0)))))) / y; else tmp = 1.0 / (1.0 / (1.0 - (t_1 * t_1))); end tmp_2 = tmp; end
code[x_, y_] := Block[{t$95$0 = N[(0.008333333333333333 + N[(N[(y * y), $MachinePrecision] * 0.0001984126984126984), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(N[(y * y), $MachinePrecision] * N[(0.16666666666666666 + N[(N[(y * y), $MachinePrecision] * t$95$0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[x, 8.5e+258], N[(N[(y * N[(1.0 + N[(N[(y * y), $MachinePrecision] * N[(0.16666666666666666 + N[(y * N[(y * t$95$0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / y), $MachinePrecision], N[(1.0 / N[(1.0 / N[(1.0 - N[(t$95$1 * t$95$1), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := 0.008333333333333333 + \left(y \cdot y\right) \cdot 0.0001984126984126984\\
t_1 := \left(y \cdot y\right) \cdot \left(0.16666666666666666 + \left(y \cdot y\right) \cdot t\_0\right)\\
\mathbf{if}\;x \leq 8.5 \cdot 10^{+258}:\\
\;\;\;\;\frac{y \cdot \left(1 + \left(y \cdot y\right) \cdot \left(0.16666666666666666 + y \cdot \left(y \cdot t\_0\right)\right)\right)}{y}\\
\mathbf{else}:\\
\;\;\;\;\frac{1}{\frac{1}{1 - t\_1 \cdot t\_1}}\\
\end{array}
\end{array}
if x < 8.49999999999999974e258Initial program 100.0%
Taylor expanded in x around 0
Simplified73.4%
Taylor expanded in y around 0
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/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-*.f6468.8%
Simplified68.8%
if 8.49999999999999974e258 < x Initial program 100.0%
Taylor expanded in x around 0
Simplified23.3%
Taylor expanded in y around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/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-*.f6415.9%
Simplified15.9%
flip-+N/A
clear-numN/A
/-lowering-/.f64N/A
/-lowering-/.f64N/A
Applied egg-rr8.2%
Taylor expanded in y around 0
Simplified46.4%
Final simplification67.7%
(FPCore (x y)
:precision binary64
(if (<= x 5.8e+261)
(/
(*
y
(+
1.0
(*
(* y y)
(+
0.16666666666666666
(*
y
(* y (+ 0.008333333333333333 (* (* y y) 0.0001984126984126984))))))))
y)
(+ 1.0 (* -0.5 (* x x)))))
double code(double x, double y) {
double tmp;
if (x <= 5.8e+261) {
tmp = (y * (1.0 + ((y * y) * (0.16666666666666666 + (y * (y * (0.008333333333333333 + ((y * y) * 0.0001984126984126984)))))))) / y;
} else {
tmp = 1.0 + (-0.5 * (x * x));
}
return tmp;
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8) :: tmp
if (x <= 5.8d+261) then
tmp = (y * (1.0d0 + ((y * y) * (0.16666666666666666d0 + (y * (y * (0.008333333333333333d0 + ((y * y) * 0.0001984126984126984d0)))))))) / y
else
tmp = 1.0d0 + ((-0.5d0) * (x * x))
end if
code = tmp
end function
public static double code(double x, double y) {
double tmp;
if (x <= 5.8e+261) {
tmp = (y * (1.0 + ((y * y) * (0.16666666666666666 + (y * (y * (0.008333333333333333 + ((y * y) * 0.0001984126984126984)))))))) / y;
} else {
tmp = 1.0 + (-0.5 * (x * x));
}
return tmp;
}
def code(x, y): tmp = 0 if x <= 5.8e+261: tmp = (y * (1.0 + ((y * y) * (0.16666666666666666 + (y * (y * (0.008333333333333333 + ((y * y) * 0.0001984126984126984)))))))) / y else: tmp = 1.0 + (-0.5 * (x * x)) return tmp
function code(x, y) tmp = 0.0 if (x <= 5.8e+261) tmp = Float64(Float64(y * Float64(1.0 + Float64(Float64(y * y) * Float64(0.16666666666666666 + Float64(y * Float64(y * Float64(0.008333333333333333 + Float64(Float64(y * y) * 0.0001984126984126984)))))))) / y); else tmp = Float64(1.0 + Float64(-0.5 * Float64(x * x))); end return tmp end
function tmp_2 = code(x, y) tmp = 0.0; if (x <= 5.8e+261) tmp = (y * (1.0 + ((y * y) * (0.16666666666666666 + (y * (y * (0.008333333333333333 + ((y * y) * 0.0001984126984126984)))))))) / y; else tmp = 1.0 + (-0.5 * (x * x)); end tmp_2 = tmp; end
code[x_, y_] := If[LessEqual[x, 5.8e+261], N[(N[(y * N[(1.0 + N[(N[(y * y), $MachinePrecision] * 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] / y), $MachinePrecision], N[(1.0 + N[(-0.5 * N[(x * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq 5.8 \cdot 10^{+261}:\\
\;\;\;\;\frac{y \cdot \left(1 + \left(y \cdot y\right) \cdot \left(0.16666666666666666 + y \cdot \left(y \cdot \left(0.008333333333333333 + \left(y \cdot y\right) \cdot 0.0001984126984126984\right)\right)\right)\right)}{y}\\
\mathbf{else}:\\
\;\;\;\;1 + -0.5 \cdot \left(x \cdot x\right)\\
\end{array}
\end{array}
if x < 5.8e261Initial program 100.0%
Taylor expanded in x around 0
Simplified73.2%
Taylor expanded in y around 0
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/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-*.f6468.6%
Simplified68.6%
if 5.8e261 < x Initial program 100.0%
Taylor expanded in y around 0
cos-lowering-cos.f6443.5%
Simplified43.5%
Taylor expanded in x around 0
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6442.4%
Simplified42.4%
Final simplification67.4%
(FPCore (x y)
:precision binary64
(if (<= x 5.8e+261)
(+
1.0
(*
(* y y)
(+
0.16666666666666666
(* y (* y (+ 0.008333333333333333 (* (* y y) 0.0001984126984126984)))))))
(+ 1.0 (* -0.5 (* x x)))))
double code(double x, double y) {
double tmp;
if (x <= 5.8e+261) {
tmp = 1.0 + ((y * y) * (0.16666666666666666 + (y * (y * (0.008333333333333333 + ((y * y) * 0.0001984126984126984))))));
} else {
tmp = 1.0 + (-0.5 * (x * x));
}
return tmp;
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8) :: tmp
if (x <= 5.8d+261) then
tmp = 1.0d0 + ((y * y) * (0.16666666666666666d0 + (y * (y * (0.008333333333333333d0 + ((y * y) * 0.0001984126984126984d0))))))
else
tmp = 1.0d0 + ((-0.5d0) * (x * x))
end if
code = tmp
end function
public static double code(double x, double y) {
double tmp;
if (x <= 5.8e+261) {
tmp = 1.0 + ((y * y) * (0.16666666666666666 + (y * (y * (0.008333333333333333 + ((y * y) * 0.0001984126984126984))))));
} else {
tmp = 1.0 + (-0.5 * (x * x));
}
return tmp;
}
def code(x, y): tmp = 0 if x <= 5.8e+261: tmp = 1.0 + ((y * y) * (0.16666666666666666 + (y * (y * (0.008333333333333333 + ((y * y) * 0.0001984126984126984)))))) else: tmp = 1.0 + (-0.5 * (x * x)) return tmp
function code(x, y) tmp = 0.0 if (x <= 5.8e+261) tmp = Float64(1.0 + Float64(Float64(y * y) * Float64(0.16666666666666666 + Float64(y * Float64(y * Float64(0.008333333333333333 + Float64(Float64(y * y) * 0.0001984126984126984))))))); else tmp = Float64(1.0 + Float64(-0.5 * Float64(x * x))); end return tmp end
function tmp_2 = code(x, y) tmp = 0.0; if (x <= 5.8e+261) tmp = 1.0 + ((y * y) * (0.16666666666666666 + (y * (y * (0.008333333333333333 + ((y * y) * 0.0001984126984126984)))))); else tmp = 1.0 + (-0.5 * (x * x)); end tmp_2 = tmp; end
code[x_, y_] := If[LessEqual[x, 5.8e+261], N[(1.0 + N[(N[(y * y), $MachinePrecision] * N[(0.16666666666666666 + N[(y * N[(y * N[(0.008333333333333333 + N[(N[(y * y), $MachinePrecision] * 0.0001984126984126984), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(1.0 + N[(-0.5 * N[(x * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq 5.8 \cdot 10^{+261}:\\
\;\;\;\;1 + \left(y \cdot y\right) \cdot \left(0.16666666666666666 + y \cdot \left(y \cdot \left(0.008333333333333333 + \left(y \cdot y\right) \cdot 0.0001984126984126984\right)\right)\right)\\
\mathbf{else}:\\
\;\;\;\;1 + -0.5 \cdot \left(x \cdot x\right)\\
\end{array}
\end{array}
if x < 5.8e261Initial program 100.0%
Taylor expanded in x around 0
Simplified73.2%
Taylor expanded in y around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/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-*.f6467.1%
Simplified67.1%
if 5.8e261 < x Initial program 100.0%
Taylor expanded in y around 0
cos-lowering-cos.f6443.5%
Simplified43.5%
Taylor expanded in x around 0
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6442.4%
Simplified42.4%
Final simplification66.0%
(FPCore (x y)
:precision binary64
(if (<= x 5.8e+261)
(+
1.0
(*
(* y y)
(+ 0.16666666666666666 (* y (* y (* (* y y) 0.0001984126984126984))))))
(+ 1.0 (* -0.5 (* x x)))))
double code(double x, double y) {
double tmp;
if (x <= 5.8e+261) {
tmp = 1.0 + ((y * y) * (0.16666666666666666 + (y * (y * ((y * y) * 0.0001984126984126984)))));
} else {
tmp = 1.0 + (-0.5 * (x * x));
}
return tmp;
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8) :: tmp
if (x <= 5.8d+261) then
tmp = 1.0d0 + ((y * y) * (0.16666666666666666d0 + (y * (y * ((y * y) * 0.0001984126984126984d0)))))
else
tmp = 1.0d0 + ((-0.5d0) * (x * x))
end if
code = tmp
end function
public static double code(double x, double y) {
double tmp;
if (x <= 5.8e+261) {
tmp = 1.0 + ((y * y) * (0.16666666666666666 + (y * (y * ((y * y) * 0.0001984126984126984)))));
} else {
tmp = 1.0 + (-0.5 * (x * x));
}
return tmp;
}
def code(x, y): tmp = 0 if x <= 5.8e+261: tmp = 1.0 + ((y * y) * (0.16666666666666666 + (y * (y * ((y * y) * 0.0001984126984126984))))) else: tmp = 1.0 + (-0.5 * (x * x)) return tmp
function code(x, y) tmp = 0.0 if (x <= 5.8e+261) tmp = Float64(1.0 + Float64(Float64(y * y) * Float64(0.16666666666666666 + Float64(y * Float64(y * Float64(Float64(y * y) * 0.0001984126984126984)))))); else tmp = Float64(1.0 + Float64(-0.5 * Float64(x * x))); end return tmp end
function tmp_2 = code(x, y) tmp = 0.0; if (x <= 5.8e+261) tmp = 1.0 + ((y * y) * (0.16666666666666666 + (y * (y * ((y * y) * 0.0001984126984126984))))); else tmp = 1.0 + (-0.5 * (x * x)); end tmp_2 = tmp; end
code[x_, y_] := If[LessEqual[x, 5.8e+261], N[(1.0 + N[(N[(y * y), $MachinePrecision] * N[(0.16666666666666666 + N[(y * N[(y * N[(N[(y * y), $MachinePrecision] * 0.0001984126984126984), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(1.0 + N[(-0.5 * N[(x * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq 5.8 \cdot 10^{+261}:\\
\;\;\;\;1 + \left(y \cdot y\right) \cdot \left(0.16666666666666666 + y \cdot \left(y \cdot \left(\left(y \cdot y\right) \cdot 0.0001984126984126984\right)\right)\right)\\
\mathbf{else}:\\
\;\;\;\;1 + -0.5 \cdot \left(x \cdot x\right)\\
\end{array}
\end{array}
if x < 5.8e261Initial program 100.0%
Taylor expanded in x around 0
Simplified73.2%
Taylor expanded in y around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/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-*.f6467.1%
Simplified67.1%
associate-*r*N/A
+-commutativeN/A
distribute-rgt-inN/A
*-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f6467.1%
Applied egg-rr67.1%
Taylor expanded in y around inf
metadata-evalN/A
pow-sqrN/A
associate-*l*N/A
*-commutativeN/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6467.1%
Simplified67.1%
if 5.8e261 < x Initial program 100.0%
Taylor expanded in y around 0
cos-lowering-cos.f6443.5%
Simplified43.5%
Taylor expanded in x around 0
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6442.4%
Simplified42.4%
Final simplification66.0%
(FPCore (x y) :precision binary64 (if (<= x 5.8e+261) (+ 1.0 (* (* y y) (* y (* y (* (* y y) 0.0001984126984126984))))) (+ 1.0 (* -0.5 (* x x)))))
double code(double x, double y) {
double tmp;
if (x <= 5.8e+261) {
tmp = 1.0 + ((y * y) * (y * (y * ((y * y) * 0.0001984126984126984))));
} else {
tmp = 1.0 + (-0.5 * (x * x));
}
return tmp;
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8) :: tmp
if (x <= 5.8d+261) then
tmp = 1.0d0 + ((y * y) * (y * (y * ((y * y) * 0.0001984126984126984d0))))
else
tmp = 1.0d0 + ((-0.5d0) * (x * x))
end if
code = tmp
end function
public static double code(double x, double y) {
double tmp;
if (x <= 5.8e+261) {
tmp = 1.0 + ((y * y) * (y * (y * ((y * y) * 0.0001984126984126984))));
} else {
tmp = 1.0 + (-0.5 * (x * x));
}
return tmp;
}
def code(x, y): tmp = 0 if x <= 5.8e+261: tmp = 1.0 + ((y * y) * (y * (y * ((y * y) * 0.0001984126984126984)))) else: tmp = 1.0 + (-0.5 * (x * x)) return tmp
function code(x, y) tmp = 0.0 if (x <= 5.8e+261) tmp = Float64(1.0 + Float64(Float64(y * y) * Float64(y * Float64(y * Float64(Float64(y * y) * 0.0001984126984126984))))); else tmp = Float64(1.0 + Float64(-0.5 * Float64(x * x))); end return tmp end
function tmp_2 = code(x, y) tmp = 0.0; if (x <= 5.8e+261) tmp = 1.0 + ((y * y) * (y * (y * ((y * y) * 0.0001984126984126984)))); else tmp = 1.0 + (-0.5 * (x * x)); end tmp_2 = tmp; end
code[x_, y_] := If[LessEqual[x, 5.8e+261], N[(1.0 + N[(N[(y * y), $MachinePrecision] * N[(y * N[(y * N[(N[(y * y), $MachinePrecision] * 0.0001984126984126984), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(1.0 + N[(-0.5 * N[(x * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq 5.8 \cdot 10^{+261}:\\
\;\;\;\;1 + \left(y \cdot y\right) \cdot \left(y \cdot \left(y \cdot \left(\left(y \cdot y\right) \cdot 0.0001984126984126984\right)\right)\right)\\
\mathbf{else}:\\
\;\;\;\;1 + -0.5 \cdot \left(x \cdot x\right)\\
\end{array}
\end{array}
if x < 5.8e261Initial program 100.0%
Taylor expanded in x around 0
Simplified73.2%
Taylor expanded in y around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/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-*.f6467.1%
Simplified67.1%
associate-*r*N/A
+-commutativeN/A
distribute-rgt-inN/A
*-commutativeN/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f6467.1%
Applied egg-rr67.1%
Taylor expanded in y around inf
metadata-evalN/A
pow-sqrN/A
associate-*l*N/A
*-commutativeN/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6466.9%
Simplified66.9%
if 5.8e261 < x Initial program 100.0%
Taylor expanded in y around 0
cos-lowering-cos.f6443.5%
Simplified43.5%
Taylor expanded in x around 0
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6442.4%
Simplified42.4%
Final simplification65.8%
(FPCore (x y) :precision binary64 (if (<= x 5.8e+261) (+ 1.0 (* y (* y (+ 0.16666666666666666 (* (* y y) 0.008333333333333333))))) (+ 1.0 (* -0.5 (* x x)))))
double code(double x, double y) {
double tmp;
if (x <= 5.8e+261) {
tmp = 1.0 + (y * (y * (0.16666666666666666 + ((y * y) * 0.008333333333333333))));
} else {
tmp = 1.0 + (-0.5 * (x * x));
}
return tmp;
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8) :: tmp
if (x <= 5.8d+261) then
tmp = 1.0d0 + (y * (y * (0.16666666666666666d0 + ((y * y) * 0.008333333333333333d0))))
else
tmp = 1.0d0 + ((-0.5d0) * (x * x))
end if
code = tmp
end function
public static double code(double x, double y) {
double tmp;
if (x <= 5.8e+261) {
tmp = 1.0 + (y * (y * (0.16666666666666666 + ((y * y) * 0.008333333333333333))));
} else {
tmp = 1.0 + (-0.5 * (x * x));
}
return tmp;
}
def code(x, y): tmp = 0 if x <= 5.8e+261: tmp = 1.0 + (y * (y * (0.16666666666666666 + ((y * y) * 0.008333333333333333)))) else: tmp = 1.0 + (-0.5 * (x * x)) return tmp
function code(x, y) tmp = 0.0 if (x <= 5.8e+261) tmp = Float64(1.0 + Float64(y * Float64(y * Float64(0.16666666666666666 + Float64(Float64(y * y) * 0.008333333333333333))))); else tmp = Float64(1.0 + Float64(-0.5 * Float64(x * x))); end return tmp end
function tmp_2 = code(x, y) tmp = 0.0; if (x <= 5.8e+261) tmp = 1.0 + (y * (y * (0.16666666666666666 + ((y * y) * 0.008333333333333333)))); else tmp = 1.0 + (-0.5 * (x * x)); end tmp_2 = tmp; end
code[x_, y_] := If[LessEqual[x, 5.8e+261], N[(1.0 + N[(y * N[(y * N[(0.16666666666666666 + N[(N[(y * y), $MachinePrecision] * 0.008333333333333333), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(1.0 + N[(-0.5 * N[(x * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq 5.8 \cdot 10^{+261}:\\
\;\;\;\;1 + y \cdot \left(y \cdot \left(0.16666666666666666 + \left(y \cdot y\right) \cdot 0.008333333333333333\right)\right)\\
\mathbf{else}:\\
\;\;\;\;1 + -0.5 \cdot \left(x \cdot x\right)\\
\end{array}
\end{array}
if x < 5.8e261Initial program 100.0%
Taylor expanded in x around 0
Simplified73.2%
Taylor expanded in y around 0
+-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-*.f6463.2%
Simplified63.2%
if 5.8e261 < x Initial program 100.0%
Taylor expanded in y around 0
cos-lowering-cos.f6443.5%
Simplified43.5%
Taylor expanded in x around 0
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6442.4%
Simplified42.4%
Final simplification62.3%
(FPCore (x y) :precision binary64 (if (<= x 5.8e+261) (+ 1.0 (* y (* y (* (* y y) 0.008333333333333333)))) (+ 1.0 (* -0.5 (* x x)))))
double code(double x, double y) {
double tmp;
if (x <= 5.8e+261) {
tmp = 1.0 + (y * (y * ((y * y) * 0.008333333333333333)));
} else {
tmp = 1.0 + (-0.5 * (x * x));
}
return tmp;
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8) :: tmp
if (x <= 5.8d+261) then
tmp = 1.0d0 + (y * (y * ((y * y) * 0.008333333333333333d0)))
else
tmp = 1.0d0 + ((-0.5d0) * (x * x))
end if
code = tmp
end function
public static double code(double x, double y) {
double tmp;
if (x <= 5.8e+261) {
tmp = 1.0 + (y * (y * ((y * y) * 0.008333333333333333)));
} else {
tmp = 1.0 + (-0.5 * (x * x));
}
return tmp;
}
def code(x, y): tmp = 0 if x <= 5.8e+261: tmp = 1.0 + (y * (y * ((y * y) * 0.008333333333333333))) else: tmp = 1.0 + (-0.5 * (x * x)) return tmp
function code(x, y) tmp = 0.0 if (x <= 5.8e+261) tmp = Float64(1.0 + Float64(y * Float64(y * Float64(Float64(y * y) * 0.008333333333333333)))); else tmp = Float64(1.0 + Float64(-0.5 * Float64(x * x))); end return tmp end
function tmp_2 = code(x, y) tmp = 0.0; if (x <= 5.8e+261) tmp = 1.0 + (y * (y * ((y * y) * 0.008333333333333333))); else tmp = 1.0 + (-0.5 * (x * x)); end tmp_2 = tmp; end
code[x_, y_] := If[LessEqual[x, 5.8e+261], N[(1.0 + N[(y * N[(y * N[(N[(y * y), $MachinePrecision] * 0.008333333333333333), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(1.0 + N[(-0.5 * N[(x * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq 5.8 \cdot 10^{+261}:\\
\;\;\;\;1 + y \cdot \left(y \cdot \left(\left(y \cdot y\right) \cdot 0.008333333333333333\right)\right)\\
\mathbf{else}:\\
\;\;\;\;1 + -0.5 \cdot \left(x \cdot x\right)\\
\end{array}
\end{array}
if x < 5.8e261Initial program 100.0%
Taylor expanded in x around 0
Simplified73.2%
Taylor expanded in y around 0
+-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-*.f6463.2%
Simplified63.2%
Taylor expanded in y around inf
metadata-evalN/A
pow-sqrN/A
associate-*l*N/A
*-commutativeN/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6463.1%
Simplified63.1%
if 5.8e261 < x Initial program 100.0%
Taylor expanded in y around 0
cos-lowering-cos.f6443.5%
Simplified43.5%
Taylor expanded in x around 0
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6442.4%
Simplified42.4%
Final simplification62.1%
(FPCore (x y) :precision binary64 (if (<= y 0.0075) (+ 1.0 (* (* y y) 0.16666666666666666)) (* y (* y (* (* y y) 0.008333333333333333)))))
double code(double x, double y) {
double tmp;
if (y <= 0.0075) {
tmp = 1.0 + ((y * y) * 0.16666666666666666);
} else {
tmp = 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 <= 0.0075d0) then
tmp = 1.0d0 + ((y * y) * 0.16666666666666666d0)
else
tmp = y * (y * ((y * y) * 0.008333333333333333d0))
end if
code = tmp
end function
public static double code(double x, double y) {
double tmp;
if (y <= 0.0075) {
tmp = 1.0 + ((y * y) * 0.16666666666666666);
} else {
tmp = y * (y * ((y * y) * 0.008333333333333333));
}
return tmp;
}
def code(x, y): tmp = 0 if y <= 0.0075: tmp = 1.0 + ((y * y) * 0.16666666666666666) else: tmp = y * (y * ((y * y) * 0.008333333333333333)) return tmp
function code(x, y) tmp = 0.0 if (y <= 0.0075) tmp = Float64(1.0 + Float64(Float64(y * y) * 0.16666666666666666)); else tmp = 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 <= 0.0075) tmp = 1.0 + ((y * y) * 0.16666666666666666); else tmp = y * (y * ((y * y) * 0.008333333333333333)); end tmp_2 = tmp; end
code[x_, y_] := If[LessEqual[y, 0.0075], N[(1.0 + N[(N[(y * y), $MachinePrecision] * 0.16666666666666666), $MachinePrecision]), $MachinePrecision], N[(y * N[(y * N[(N[(y * y), $MachinePrecision] * 0.008333333333333333), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq 0.0075:\\
\;\;\;\;1 + \left(y \cdot y\right) \cdot 0.16666666666666666\\
\mathbf{else}:\\
\;\;\;\;y \cdot \left(y \cdot \left(\left(y \cdot y\right) \cdot 0.008333333333333333\right)\right)\\
\end{array}
\end{array}
if y < 0.0074999999999999997Initial program 100.0%
Taylor expanded in x around 0
Simplified70.1%
Taylor expanded in y around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6456.0%
Simplified56.0%
if 0.0074999999999999997 < y Initial program 100.0%
Taylor expanded in x around 0
Simplified73.2%
Taylor expanded in y around 0
+-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-*.f6457.9%
Simplified57.9%
Taylor expanded in y around inf
metadata-evalN/A
pow-sqrN/A
associate-*l*N/A
*-commutativeN/A
unpow2N/A
associate-*l*N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6457.9%
Simplified57.9%
Final simplification56.5%
(FPCore (x y) :precision binary64 (if (<= x 9.8e+163) (+ 1.0 (* (* y y) 0.16666666666666666)) (+ 1.0 (* -0.5 (* x x)))))
double code(double x, double y) {
double tmp;
if (x <= 9.8e+163) {
tmp = 1.0 + ((y * y) * 0.16666666666666666);
} else {
tmp = 1.0 + (-0.5 * (x * x));
}
return tmp;
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
real(8) :: tmp
if (x <= 9.8d+163) then
tmp = 1.0d0 + ((y * y) * 0.16666666666666666d0)
else
tmp = 1.0d0 + ((-0.5d0) * (x * x))
end if
code = tmp
end function
public static double code(double x, double y) {
double tmp;
if (x <= 9.8e+163) {
tmp = 1.0 + ((y * y) * 0.16666666666666666);
} else {
tmp = 1.0 + (-0.5 * (x * x));
}
return tmp;
}
def code(x, y): tmp = 0 if x <= 9.8e+163: tmp = 1.0 + ((y * y) * 0.16666666666666666) else: tmp = 1.0 + (-0.5 * (x * x)) return tmp
function code(x, y) tmp = 0.0 if (x <= 9.8e+163) tmp = Float64(1.0 + Float64(Float64(y * y) * 0.16666666666666666)); else tmp = Float64(1.0 + Float64(-0.5 * Float64(x * x))); end return tmp end
function tmp_2 = code(x, y) tmp = 0.0; if (x <= 9.8e+163) tmp = 1.0 + ((y * y) * 0.16666666666666666); else tmp = 1.0 + (-0.5 * (x * x)); end tmp_2 = tmp; end
code[x_, y_] := If[LessEqual[x, 9.8e+163], N[(1.0 + N[(N[(y * y), $MachinePrecision] * 0.16666666666666666), $MachinePrecision]), $MachinePrecision], N[(1.0 + N[(-0.5 * N[(x * x), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;x \leq 9.8 \cdot 10^{+163}:\\
\;\;\;\;1 + \left(y \cdot y\right) \cdot 0.16666666666666666\\
\mathbf{else}:\\
\;\;\;\;1 + -0.5 \cdot \left(x \cdot x\right)\\
\end{array}
\end{array}
if x < 9.8e163Initial program 100.0%
Taylor expanded in x around 0
Simplified76.8%
Taylor expanded in y around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6456.4%
Simplified56.4%
if 9.8e163 < x Initial program 99.9%
Taylor expanded in y around 0
cos-lowering-cos.f6443.4%
Simplified43.4%
Taylor expanded in x around 0
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6437.4%
Simplified37.4%
Final simplification54.1%
(FPCore (x y) :precision binary64 (if (<= y 0.0075) 1.0 (* (* y y) 0.16666666666666666)))
double code(double x, double y) {
double tmp;
if (y <= 0.0075) {
tmp = 1.0;
} else {
tmp = (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 <= 0.0075d0) then
tmp = 1.0d0
else
tmp = (y * y) * 0.16666666666666666d0
end if
code = tmp
end function
public static double code(double x, double y) {
double tmp;
if (y <= 0.0075) {
tmp = 1.0;
} else {
tmp = (y * y) * 0.16666666666666666;
}
return tmp;
}
def code(x, y): tmp = 0 if y <= 0.0075: tmp = 1.0 else: tmp = (y * y) * 0.16666666666666666 return tmp
function code(x, y) tmp = 0.0 if (y <= 0.0075) tmp = 1.0; else tmp = Float64(Float64(y * y) * 0.16666666666666666); end return tmp end
function tmp_2 = code(x, y) tmp = 0.0; if (y <= 0.0075) tmp = 1.0; else tmp = (y * y) * 0.16666666666666666; end tmp_2 = tmp; end
code[x_, y_] := If[LessEqual[y, 0.0075], 1.0, N[(N[(y * y), $MachinePrecision] * 0.16666666666666666), $MachinePrecision]]
\begin{array}{l}
\\
\begin{array}{l}
\mathbf{if}\;y \leq 0.0075:\\
\;\;\;\;1\\
\mathbf{else}:\\
\;\;\;\;\left(y \cdot y\right) \cdot 0.16666666666666666\\
\end{array}
\end{array}
if y < 0.0074999999999999997Initial program 100.0%
Taylor expanded in y around 0
cos-lowering-cos.f6464.5%
Simplified64.5%
Taylor expanded in x around 0
Simplified41.2%
if 0.0074999999999999997 < y Initial program 100.0%
Taylor expanded in x around 0
Simplified73.2%
Taylor expanded in y around 0
+-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-*.f6457.9%
Simplified57.9%
Taylor expanded in y around 0
+-lowering-+.f64N/A
*-commutativeN/A
unpow2N/A
associate-*l*N/A
*-commutativeN/A
*-lowering-*.f64N/A
*-commutativeN/A
*-lowering-*.f6438.2%
Simplified38.2%
Taylor expanded in y around inf
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6438.2%
Simplified38.2%
Final simplification40.4%
(FPCore (x y) :precision binary64 (+ 1.0 (* (* y y) 0.16666666666666666)))
double code(double x, double y) {
return 1.0 + ((y * y) * 0.16666666666666666);
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
code = 1.0d0 + ((y * y) * 0.16666666666666666d0)
end function
public static double code(double x, double y) {
return 1.0 + ((y * y) * 0.16666666666666666);
}
def code(x, y): return 1.0 + ((y * y) * 0.16666666666666666)
function code(x, y) return Float64(1.0 + Float64(Float64(y * y) * 0.16666666666666666)) end
function tmp = code(x, y) tmp = 1.0 + ((y * y) * 0.16666666666666666); end
code[x_, y_] := N[(1.0 + N[(N[(y * y), $MachinePrecision] * 0.16666666666666666), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
1 + \left(y \cdot y\right) \cdot 0.16666666666666666
\end{array}
Initial program 100.0%
Taylor expanded in x around 0
Simplified70.9%
Taylor expanded in y around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
unpow2N/A
*-lowering-*.f6451.3%
Simplified51.3%
Final simplification51.3%
(FPCore (x y) :precision binary64 1.0)
double code(double x, double y) {
return 1.0;
}
real(8) function code(x, y)
real(8), intent (in) :: x
real(8), intent (in) :: y
code = 1.0d0
end function
public static double code(double x, double y) {
return 1.0;
}
def code(x, y): return 1.0
function code(x, y) return 1.0 end
function tmp = code(x, y) tmp = 1.0; end
code[x_, y_] := 1.0
\begin{array}{l}
\\
1
\end{array}
Initial program 100.0%
Taylor expanded in y around 0
cos-lowering-cos.f6448.6%
Simplified48.6%
Taylor expanded in x around 0
Simplified31.1%
herbie shell --seed 2024155
(FPCore (x y)
:name "Linear.Quaternion:$csin from linear-1.19.1.3"
:precision binary64
(* (cos x) (/ (sinh y) y)))