| 1870 | |
| 1871 | |
| 1872 | void cftf162(float *a) |
| 1873 | { |
| 1874 | float wn4r, wk1r, wk1i, wk2r, wk2i, wk3r, wk3i, |
| 1875 | x0r, x0i, x1r, x1i, x2r, x2i, |
| 1876 | y0r, y0i, y1r, y1i, y2r, y2i, y3r, y3i, |
| 1877 | y4r, y4i, y5r, y5i, y6r, y6i, y7r, y7i, |
| 1878 | y8r, y8i, y9r, y9i, y10r, y10i, y11r, y11i, |
| 1879 | y12r, y12i, y13r, y13i, y14r, y14i, y15r, y15i; |
| 1880 | |
| 1881 | wn4r = WR5000; |
| 1882 | wk1r = WR1250; |
| 1883 | wk1i = WI1250; |
| 1884 | wk2r = WR2500; |
| 1885 | wk2i = WI2500; |
| 1886 | wk3r = WR3750; |
| 1887 | wk3i = WI3750; |
| 1888 | x1r = a[0] - a[17]; |
| 1889 | x1i = a[1] + a[16]; |
| 1890 | x0r = a[8] - a[25]; |
| 1891 | x0i = a[9] + a[24]; |
| 1892 | x2r = wn4r * (x0r - x0i); |
| 1893 | x2i = wn4r * (x0i + x0r); |
| 1894 | y0r = x1r + x2r; |
| 1895 | y0i = x1i + x2i; |
| 1896 | y4r = x1r - x2r; |
| 1897 | y4i = x1i - x2i; |
| 1898 | x1r = a[0] + a[17]; |
| 1899 | x1i = a[1] - a[16]; |
| 1900 | x0r = a[8] + a[25]; |
| 1901 | x0i = a[9] - a[24]; |
| 1902 | x2r = wn4r * (x0r - x0i); |
| 1903 | x2i = wn4r * (x0i + x0r); |
| 1904 | y8r = x1r - x2i; |
| 1905 | y8i = x1i + x2r; |
| 1906 | y12r = x1r + x2i; |
| 1907 | y12i = x1i - x2r; |
| 1908 | x0r = a[2] - a[19]; |
| 1909 | x0i = a[3] + a[18]; |
| 1910 | x1r = wk1r * x0r - wk1i * x0i; |
| 1911 | x1i = wk1r * x0i + wk1i * x0r; |
| 1912 | x0r = a[10] - a[27]; |
| 1913 | x0i = a[11] + a[26]; |
| 1914 | x2r = wk3i * x0r - wk3r * x0i; |
| 1915 | x2i = wk3i * x0i + wk3r * x0r; |
| 1916 | y1r = x1r + x2r; |
| 1917 | y1i = x1i + x2i; |
| 1918 | y5r = x1r - x2r; |
| 1919 | y5i = x1i - x2i; |
| 1920 | x0r = a[2] + a[19]; |
| 1921 | x0i = a[3] - a[18]; |
| 1922 | x1r = wk3r * x0r - wk3i * x0i; |
| 1923 | x1i = wk3r * x0i + wk3i * x0r; |
| 1924 | x0r = a[10] + a[27]; |
| 1925 | x0i = a[11] - a[26]; |
| 1926 | x2r = wk1r * x0r + wk1i * x0i; |
| 1927 | x2i = wk1r * x0i - wk1i * x0r; |
| 1928 | y9r = x1r - x2r; |
| 1929 | y9i = x1i - x2i; |