indextables.cpp

普林斯顿开发的快速球面调和变换算法

  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};




indextables::indextables()
{

}

indextables::~indextables()
{

}

const int* indextables::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;
    }


}



void indextables::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);


}