3369, 809, 2857, 1833, 3881, 169, 2217, 1193, 3241, 681, 2729, 1705, 3753,   425, 2473, 1449, 3497, 937, 2985, 1961, 4009, 105, 2153, 1129, 3177, 617,   2665, 1641, 3689, 361, 2409, 1385, 3433, 873, 2921, 1897, 3945, 233, 2281,   1257, 3305, 745, 2793, 1769, 3817, 489, 2537, 1513, 3561, 1001, 3049, 2025,   4073, 25, 2073, 1049, 3097, 537, 2585, 1561, 3609, 281, 2329, 1305, 3353,   793, 2841, 1817, 3865, 153, 2201, 1177, 3225, 665, 2713, 1689, 3737, 409,   2457, 1433, 3481, 921, 2969, 1945, 3993, 89, 2137, 1113, 3161, 601, 2649,   1625, 3673, 345, 2393, 1369, 3417, 857, 2905, 1881, 3929, 217, 2265, 1241,   3289, 729, 2777, 1753, 3801, 473, 2521, 1497, 3545, 985, 3033, 2009, 4057,   57, 2105, 1081, 3129, 569, 2617, 1593, 3641, 313, 2361, 1337, 3385, 825,   2873, 1849, 3897, 185, 2233, 1209, 3257, 697, 2745, 1721, 3769, 441, 2489,   1465, 3513, 953, 3001, 1977, 4025, 121, 2169, 1145, 3193, 633, 2681, 1657,   3705, 377, 2425, 1401, 3449, 889, 2937, 1913, 3961, 249, 2297, 1273, 3321,   761, 2809, 1785, 3833, 505, 2553, 1529, 3577, 1017, 3065, 2041, 4089, 5,   2053, 1029, 3077, 517, 2565, 1541, 3589, 261, 2309, 1285, 3333, 773, 2821,   1797, 3845, 133, 2181, 1157, 3205, 645, 2693, 1669, 3717, 389, 2437, 1413,   3461, 901, 2949, 1925, 3973, 69, 2117, 1093, 3141, 581, 2629, 1605, 3653,   325, 2373, 1349, 3397, 837, 2885, 1861, 3909, 197, 2245, 1221, 3269, 709,   2757, 1733, 3781, 453, 2501, 1477, 3525, 965, 3013, 1989, 4037, 37, 2085,   1061, 3109, 549, 2597, 1573, 3621, 293, 2341, 1317, 3365, 805, 2853, 1829,   3877, 165, 2213, 1189, 3237, 677, 2725, 1701, 3749, 421, 2469, 1445, 3493,   933, 2981, 1957, 4005, 101, 2149, 1125, 3173, 613, 2661, 1637, 3685, 357,   2405, 1381, 3429, 869, 2917, 1893, 3941, 229, 2277, 1253, 3301, 741, 2789,   1765, 3813, 485, 2533, 1509, 3557, 997, 3045, 2021, 4069, 21, 2069, 1045,   3093, 533, 2581, 1557, 3605, 277, 2325, 1301, 3349, 789, 2837, 1813, 3861,   149, 2197, 1173, 3221, 661, 2709, 1685, 3733, 405, 2453, 1429, 3477, 917,   2965, 1941, 3989, 85, 2133, 1109, 3157, 597, 2645, 1621, 3669, 341, 2389,   1365, 3413, 853, 2901, 1877, 3925, 213, 2261, 1237, 3285, 725, 2773, 1749,   3797, 469, 2517, 1493, 3541, 981, 3029, 2005, 4053, 53, 2101, 1077, 3125,   565, 2613, 1589, 3637, 309, 2357, 1333, 3381, 821, 2869, 1845, 3893, 181,   2229, 1205, 3253, 693, 2741, 1717, 3765, 437, 2485, 1461, 3509, 949, 2997,   1973, 4021, 117, 2165, 1141, 3189, 629, 2677, 1653, 3701, 373, 2421, 1397,   3445, 885, 2933, 1909, 3957, 245, 2293, 1269, 3317, 757, 2805, 1781, 3829,   501, 2549, 1525, 3573, 1013, 3061, 2037, 4085, 13, 2061, 1037, 3085, 525,   2573, 1549, 3597, 269, 2317, 1293, 3341, 781, 2829, 1805, 3853, 141, 2189,   1165, 3213, 653, 2701, 1677, 3725, 397, 2445, 1421, 3469, 909, 2957, 1933,   3981, 77, 2125, 1101, 3149, 589, 2637, 1613, 3661, 333, 2381, 1357, 3405,   845, 2893, 1869, 3917, 205, 2253, 1229, 3277, 717, 2765, 1741, 3789, 461,   2509, 1485, 3533, 973, 3021, 1997, 4045, 45, 2093, 1069, 3117, 557, 2605,   1581, 3629, 301, 2349, 1325, 3373, 813, 2861, 1837, 3885, 173, 2221, 1197,   3245, 685, 2733, 1709, 3757, 429, 2477, 1453, 3501, 941, 2989, 1965, 4013,   109, 2157, 1133, 3181, 621, 2669, 1645, 3693, 365, 2413, 1389, 3437, 877,   2925, 1901, 3949, 237, 2285, 1261, 3309, 749, 2797, 1773, 3821, 493, 2541,   1517, 3565, 1005, 3053, 2029, 4077, 29, 2077, 1053, 3101, 541, 2589, 1565,   3613, 285, 2333, 1309, 3357, 797, 2845, 1821, 3869, 157, 2205, 1181, 3229,   669, 2717, 1693, 3741, 413, 2461, 1437, 3485, 925, 2973, 1949, 3997, 93,   2141, 1117, 3165, 605, 2653, 1629, 3677, 349, 2397, 1373, 3421, 861, 2909,   1885, 3933, 221, 2269, 1245, 3293, 733, 2781, 1757, 3805, 477, 2525, 1501,   3549, 989, 3037, 2013, 4061, 61, 2109, 1085, 3133, 573, 2621, 1597, 3645,   317, 2365, 1341, 3389, 829, 2877, 1853, 3901, 189, 2237, 1213, 3261, 701,   2749, 1725, 3773, 445, 2493, 1469, 3517, 957, 3005, 1981, 4029, 125, 2173,   1149, 3197, 637, 2685, 1661, 3709, 381, 2429, 1405, 3453, 893, 2941, 1917,   3965, 253, 2301, 1277, 3325, 765, 2813, 1789, 3837, 509, 2557, 1533, 3581,   1021, 3069, 2045, 4093, 3, 2051, 1027, 3075, 515, 2563, 1539, 3587, 259,   2307, 1283, 3331, 771, 2819, 1795, 3843, 131, 2179, 1155, 3203, 643, 2691,   1667, 3715, 387, 2435, 1411, 3459, 899, 2947, 1923, 3971, 67, 2115, 1091,   3139, 579, 2627, 1603, 3651, 323, 2371, 1347, 3395, 835, 2883, 1859, 3907,   195, 2243, 1219, 3267, 707, 2755, 1731, 3779, 451, 2499, 1475, 3523, 963,   3011, 1987, 4035, 35, 2083, 1059, 3107, 547, 2595, 1571, 3619, 291, 2339,   1315, 3363, 803, 2851, 1827, 3875, 163, 2211, 1187, 3235, 675, 2723, 1699,   3747, 419, 2467, 1443, 3491, 931, 2979, 1955, 4003, 99, 2147, 1123, 3171,   611, 2659, 1635, 3683, 355, 2403, 1379, 3427, 867, 2915, 1891, 3939, 227,   2275, 1251, 3299, 739, 2787, 1763, 3811, 483, 2531, 1507, 3555, 995, 3043,   2019, 4067, 19, 2067, 1043, 3091, 531, 2579, 1555, 3603, 275, 2323, 1299,   3347, 787, 2835, 1811, 3859, 147, 2195, 1171, 3219, 659, 2707, 1683, 3731,   403, 2451, 1427, 3475, 915, 2963, 1939, 3987, 83, 2131, 1107, 3155, 595,   2643, 1619, 3667, 339, 2387, 1363, 3411, 851, 2899, 1875, 3923, 211, 2259,   1235, 3283, 723, 2771, 1747, 3795, 467, 2515, 1491, 3539, 979, 3027, 2003,   4051, 51, 2099, 1075, 3123, 563, 2611, 1587, 3635, 307, 2355, 1331, 3379,   819, 2867, 1843, 3891, 179, 2227, 1203, 3251, 691, 2739, 1715, 3763, 435,   2483, 1459, 3507, 947, 2995, 1971, 4019, 115, 2163, 1139, 3187, 627, 2675,   1651, 3699, 371, 2419, 1395, 3443, 883, 2931, 1907, 3955, 243, 2291, 1267,   3315, 755, 2803, 1779, 3827, 499, 2547, 1523, 3571, 1011, 3059, 2035, 4083,   11, 2059, 1035, 3083, 523, 2571, 1547, 3595, 267, 2315, 1291, 3339, 779,   2827, 1803, 3851, 139, 2187, 1163, 3211, 651, 2699, 1675, 3723, 395, 2443,   1419, 3467, 907, 2955, 1931, 3979, 75, 2123, 1099, 3147, 587, 2635, 1611,   3659, 331, 2379, 1355, 3403, 843, 2891, 1867, 3915, 203, 2251, 1227, 3275,   715, 2763, 1739, 3787, 459, 2507, 1483, 3531, 971, 3019, 1995, 4043, 43,   2091, 1067, 3115, 555, 2603, 1579, 3627, 299, 2347, 1323, 3371, 811, 2859,   1835, 3883, 171, 2219, 1195, 3243, 683, 2731, 1707, 3755, 427, 2475, 1451,   3499, 939, 2987, 1963, 4011, 107, 2155, 1131, 3179, 619, 2667, 1643, 3691,   363, 2411, 1387, 3435, 875, 2923, 1899, 3947, 235, 2283, 1259, 3307, 747,   2795, 1771, 3819, 491, 2539, 1515, 3563, 1003, 3051, 2027, 4075, 27, 2075,   1051, 3099, 539, 2587, 1563, 3611, 283, 2331, 1307, 3355, 795, 2843, 1819,   3867, 155, 2203, 1179, 3227, 667, 2715, 1691, 3739, 411, 2459, 1435, 3483,   923, 2971, 1947, 3995, 91, 2139, 1115, 3163, 603, 2651, 1627, 3675, 347,   2395, 1371, 3419, 859, 2907, 1883, 3931, 219, 2267, 1243, 3291, 731, 2779,   1755, 3803, 475, 2523, 1499, 3547, 987, 3035, 2011, 4059, 59, 2107, 1083,   3131, 571, 2619, 1595, 3643, 315, 2363, 1339, 3387, 827, 2875, 1851, 3899,   187, 2235, 1211, 3259, 699, 2747, 1723, 3771, 443, 2491, 1467, 3515, 955,   3003, 1979, 4027, 123, 2171, 1147, 3195, 635, 2683, 1659, 3707, 379, 2427,   1403, 3451, 891, 2939, 1915, 3963, 251, 2299, 1275, 3323, 763, 2811, 1787,   3835, 507, 2555, 1531, 3579, 1019, 3067, 2043, 4091, 7, 2055, 1031, 3079,   519, 2567, 1543, 3591, 263, 2311, 1287, 3335, 775, 2823, 1799, 3847, 135,   2183, 1159, 3207, 647, 2695, 1671, 3719, 391, 2439, 1415, 3463, 903, 2951,   1927, 3975, 71, 2119, 1095, 3143, 583, 2631, 1607, 3655, 327, 2375, 1351,   3399, 839, 2887, 1863, 3911, 199, 2247, 1223, 3271, 711, 2759, 1735, 3783,   455, 2503, 1479, 3527, 967, 3015, 1991, 4039, 39, 2087, 1063, 3111, 551,   2599, 1575, 3623, 295, 2343, 1319, 3367, 807, 2855, 1831, 3879, 167, 2215,   1191, 3239, 679, 2727, 1703, 3751, 423, 2471, 1447, 3495, 935, 2983, 1959,   4007, 103, 2151, 1127, 3175, 615, 2663, 1639, 3687, 359, 2407, 1383, 3431,   871, 2919, 1895, 3943, 231, 2279, 1255, 3303, 743, 2791, 1767, 3815, 487,   2535, 1511, 3559, 999, 3047, 2023, 4071, 23, 2071, 1047, 3095, 535, 2583,   1559, 3607, 279, 2327, 1303, 3351, 791, 2839, 1815, 3863, 151, 2199, 1175,   3223, 663, 2711, 1687, 3735, 407, 2455, 1431, 3479, 919, 2967, 1943, 3991,   87, 2135, 1111, 3159, 599, 2647, 1623, 3671, 343, 2391, 1367, 3415, 855,   2903, 1879, 3927, 215, 2263, 1239, 3287, 727, 2775, 1751, 3799, 471, 2519,   1495, 3543, 983, 3031, 2007, 4055, 55, 2103, 1079, 3127, 567, 2615, 1591,   3639, 311, 2359, 1335, 3383, 823, 2871, 1847, 3895, 183, 2231, 1207, 3255,   695, 2743, 1719, 3767, 439, 2487, 1463, 3511, 951, 2999, 1975, 4023, 119,   2167, 1143, 3191, 631, 2679, 1655, 3703, 375, 2423, 1399, 3447, 887, 2935,   1911, 3959, 247, 2295, 1271, 3319, 759, 2807, 1783, 3831, 503, 2551, 1527,   3575, 1015, 3063, 2039, 4087, 15, 2063, 1039, 3087, 527, 2575, 1551, 3599,   271, 2319, 1295, 3343, 783, 2831, 1807, 3855, 143, 2191, 1167, 3215, 655,   2703, 1679, 3727, 399, 2447, 1423, 3471, 911, 2959, 1935, 3983, 79, 2127,   1103, 3151, 591, 2639, 1615, 3663, 335, 2383, 1359, 3407, 847, 2895, 1871,   3919, 207, 2255, 1231, 3279, 719, 2767, 1743, 3791, 463, 2511, 1487, 3535,   975, 3023, 1999, 4047, 47, 2095, 1071, 3119, 559, 2607, 1583, 3631, 303,   2351, 1327, 3375, 815, 2863, 1839, 3887, 175, 2223, 1199, 3247, 687, 2735,   1711, 3759, 431, 2479, 1455, 3503, 943, 2991, 1967, 4015, 111, 2159, 1135,   3183, 623, 2671, 1647, 3695, 367, 2415, 1391, 3439, 879, 2927, 1903, 3951,   239, 2287, 1263, 3311, 751, 2799, 1775, 3823, 495, 2543, 1519, 3567, 1007,   3055, 2031, 4079, 31, 2079, 1055, 3103, 543, 2591, 1567, 3615, 287, 2335,   1311, 3359, 799, 2847, 1823, 3871, 159, 2207, 1183, 3231, 671, 2719, 1695,   3743, 415, 2463, 1439, 3487, 927, 2975, 1951, 3999, 95, 2143, 1119, 3167,   607, 2655, 1631, 3679, 351, 2399, 1375, 3423, 863, 2911, 1887, 3935, 223,   2271, 1247, 3295, 735, 2783, 1759, 3807, 479, 2527, 1503, 3551, 991, 3039,   2015, 4063, 63, 2111, 1087, 3135, 575, 2623, 1599, 3647, 319, 2367, 1343,   3391, 831, 2879, 1855, 3903, 191, 2239, 1215, 3263, 703, 2751, 1727, 3775,   447, 2495, 1471, 3519, 959, 3007, 1983, 4031, 127, 2175, 1151, 3199, 639,   2687, 1663, 3711, 383, 2431, 1407, 3455, 895, 2943, 1919, 3967, 255, 2303,   1279, 3327, 767, 2815, 1791, 3839, 511, 2559, 1535, 3583, 1023, 3071, 2047,   4095};/************************************************************************/static const int *ptable(int n){  switch (n)    {    case 2: return p2;    case 4: return p4;    case 8: return p8;    case 16: return p16;    case 32: return p32;    case 64: return p64;    case 128: return p128;    case 256: return p256;    case 512: return p512;    case 1024: return p1024;    case 2048: return p2048;    case 4096: return p4096;    default: return 0;    }}/************************************************************************//* performs a bit-reversal on the data, a double array of size n *//* Needs workspace of size n */void bitreverse(double *data,		int n,		double *workspace){    int i;  const int *perm;  perm = ptable(n);    for (i=0; i<n; i++)    workspace[i] = data[perm[i]];  memcpy(data, workspace, sizeof(double) * n);}