| 1711 | } |
| 1712 | |
| 1713 | void cftf161(float *a) |
| 1714 | { |
| 1715 | float wn4r, wk1r, wk1i, |
| 1716 | x0r, x0i, x1r, x1i, x2r, x2i, x3r, x3i, |
| 1717 | y0r, y0i, y1r, y1i, y2r, y2i, y3r, y3i, |
| 1718 | y4r, y4i, y5r, y5i, y6r, y6i, y7r, y7i, |
| 1719 | y8r, y8i, y9r, y9i, y10r, y10i, y11r, y11i, |
| 1720 | y12r, y12i, y13r, y13i, y14r, y14i, y15r, y15i; |
| 1721 | |
| 1722 | wn4r = WR5000; |
| 1723 | wk1r = WR2500; |
| 1724 | wk1i = WI2500; |
| 1725 | x0r = a[0] + a[16]; |
| 1726 | x0i = a[1] + a[17]; |
| 1727 | x1r = a[0] - a[16]; |
| 1728 | x1i = a[1] - a[17]; |
| 1729 | x2r = a[8] + a[24]; |
| 1730 | x2i = a[9] + a[25]; |
| 1731 | x3r = a[8] - a[24]; |
| 1732 | x3i = a[9] - a[25]; |
| 1733 | y0r = x0r + x2r; |
| 1734 | y0i = x0i + x2i; |
| 1735 | y4r = x0r - x2r; |
| 1736 | y4i = x0i - x2i; |
| 1737 | y8r = x1r - x3i; |
| 1738 | y8i = x1i + x3r; |
| 1739 | y12r = x1r + x3i; |
| 1740 | y12i = x1i - x3r; |
| 1741 | x0r = a[2] + a[18]; |
| 1742 | x0i = a[3] + a[19]; |
| 1743 | x1r = a[2] - a[18]; |
| 1744 | x1i = a[3] - a[19]; |
| 1745 | x2r = a[10] + a[26]; |
| 1746 | x2i = a[11] + a[27]; |
| 1747 | x3r = a[10] - a[26]; |
| 1748 | x3i = a[11] - a[27]; |
| 1749 | y1r = x0r + x2r; |
| 1750 | y1i = x0i + x2i; |
| 1751 | y5r = x0r - x2r; |
| 1752 | y5i = x0i - x2i; |
| 1753 | x0r = x1r - x3i; |
| 1754 | x0i = x1i + x3r; |
| 1755 | y9r = wk1r * x0r - wk1i * x0i; |
| 1756 | y9i = wk1r * x0i + wk1i * x0r; |
| 1757 | x0r = x1r + x3i; |
| 1758 | x0i = x1i - x3r; |
| 1759 | y13r = wk1i * x0r - wk1r * x0i; |
| 1760 | y13i = wk1i * x0i + wk1r * x0r; |
| 1761 | x0r = a[4] + a[20]; |
| 1762 | x0i = a[5] + a[21]; |
| 1763 | x1r = a[4] - a[20]; |
| 1764 | x1i = a[5] - a[21]; |
| 1765 | x2r = a[12] + a[28]; |
| 1766 | x2i = a[13] + a[29]; |
| 1767 | x3r = a[12] - a[28]; |
| 1768 | x3i = a[13] - a[29]; |
| 1769 | y2r = x0r + x2r; |
| 1770 | y2i = x0i + x2i; |