newfct.cpp

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

  1743, 1328, 719, 816, 1231, 1840, 207, 176, 1871, 1200, 847, 688, 1359, 
  1712, 335, 432, 1615, 1456, 591, 944, 1103, 1968, 79, 112, 1935, 1136, 911, 
  624, 1423, 1648, 399, 368, 1679, 1392, 655, 880, 1167, 1904, 143, 240, 
  1807, 1264, 783, 752, 1295, 1776, 271, 496, 1551, 1520, 527, 1008, 1039, 
  2032, 15, 8, 2039, 1032, 1015, 520, 1527, 1544, 503, 264, 1783, 1288, 759, 
  776, 1271, 1800, 247, 136, 1911, 1160, 887, 648, 1399, 1672, 375, 392, 
  1655, 1416, 631, 904, 1143, 1928, 119, 72, 1975, 1096, 951, 584, 1463, 
  1608, 439, 328, 1719, 1352, 695, 840, 1207, 1864, 183, 200, 1847, 1224, 
  823, 712, 1335, 1736, 311, 456, 1591, 1480, 567, 968, 1079, 1992, 55, 40, 
  2007, 1064, 983, 552, 1495, 1576, 471, 296, 1751, 1320, 727, 808, 1239, 
  1832, 215, 168, 1879, 1192, 855, 680, 1367, 1704, 343, 424, 1623, 1448, 
  599, 936, 1111, 1960, 87, 104, 1943, 1128, 919, 616, 1431, 1640, 407, 360, 
  1687, 1384, 663, 872, 1175, 1896, 151, 232, 1815, 1256, 791, 744, 1303, 
  1768, 279, 488, 1559, 1512, 535, 1000, 1047, 2024, 23, 24, 2023, 1048, 999, 
  536, 1511, 1560, 487, 280, 1767, 1304, 743, 792, 1255, 1816, 231, 152, 
  1895, 1176, 871, 664, 1383, 1688, 359, 408, 1639, 1432, 615, 920, 1127, 
  1944, 103, 88, 1959, 1112, 935, 600, 1447, 1624, 423, 344, 1703, 1368, 679, 
  856, 1191, 1880, 167, 216, 1831, 1240, 807, 728, 1319, 1752, 295, 472, 
  1575, 1496, 551, 984, 1063, 2008, 39, 56, 1991, 1080, 967, 568, 1479, 1592, 
  455, 312, 1735, 1336, 711, 824, 1223, 1848, 199, 184, 1863, 1208, 839, 696, 
  1351, 1720, 327, 440, 1607, 1464, 583, 952, 1095, 1976, 71, 120, 1927, 
  1144, 903, 632, 1415, 1656, 391, 376, 1671, 1400, 647, 888, 1159, 1912, 
  135, 248, 1799, 1272, 775, 760, 1287, 1784, 263, 504, 1543, 1528, 519, 
  1016, 1031, 2040, 7, 4, 2043, 1028, 1019, 516, 1531, 1540, 507, 260, 1787, 
  1284, 763, 772, 1275, 1796, 251, 132, 1915, 1156, 891, 644, 1403, 1668, 
  379, 388, 1659, 1412, 635, 900, 1147, 1924, 123, 68, 1979, 1092, 955, 580, 
  1467, 1604, 443, 324, 1723, 1348, 699, 836, 1211, 1860, 187, 196, 1851, 
  1220, 827, 708, 1339, 1732, 315, 452, 1595, 1476, 571, 964, 1083, 1988, 59, 
  36, 2011, 1060, 987, 548, 1499, 1572, 475, 292, 1755, 1316, 731, 804, 1243, 
  1828, 219, 164, 1883, 1188, 859, 676, 1371, 1700, 347, 420, 1627, 1444, 
  603, 932, 1115, 1956, 91, 100, 1947, 1124, 923, 612, 1435, 1636, 411, 356, 
  1691, 1380, 667, 868, 1179, 1892, 155, 228, 1819, 1252, 795, 740, 1307, 
  1764, 283, 484, 1563, 1508, 539, 996, 1051, 2020, 27, 20, 2027, 1044, 1003, 
  532, 1515, 1556, 491, 276, 1771, 1300, 747, 788, 1259, 1812, 235, 148, 
  1899, 1172, 875, 660, 1387, 1684, 363, 404, 1643, 1428, 619, 916, 1131, 
  1940, 107, 84, 1963, 1108, 939, 596, 1451, 1620, 427, 340, 1707, 1364, 683, 
  852, 1195, 1876, 171, 212, 1835, 1236, 811, 724, 1323, 1748, 299, 468, 
  1579, 1492, 555, 980, 1067, 2004, 43, 52, 1995, 1076, 971, 564, 1483, 1588, 
  459, 308, 1739, 1332, 715, 820, 1227, 1844, 203, 180, 1867, 1204, 843, 692, 
  1355, 1716, 331, 436, 1611, 1460, 587, 948, 1099, 1972, 75, 116, 1931, 
  1140, 907, 628, 1419, 1652, 395, 372, 1675, 1396, 651, 884, 1163, 1908, 
  139, 244, 1803, 1268, 779, 756, 1291, 1780, 267, 500, 1547, 1524, 523, 
  1012, 1035, 2036, 11, 12, 2035, 1036, 1011, 524, 1523, 1548, 499, 268, 
  1779, 1292, 755, 780, 1267, 1804, 243, 140, 1907, 1164, 883, 652, 1395, 
  1676, 371, 396, 1651, 1420, 627, 908, 1139, 1932, 115, 76, 1971, 1100, 947, 
  588, 1459, 1612, 435, 332, 1715, 1356, 691, 844, 1203, 1868, 179, 204, 
  1843, 1228, 819, 716, 1331, 1740, 307, 460, 1587, 1484, 563, 972, 1075, 
  1996, 51, 44, 2003, 1068, 979, 556, 1491, 1580, 467, 300, 1747, 1324, 723, 
  812, 1235, 1836, 211, 172, 1875, 1196, 851, 684, 1363, 1708, 339, 428, 
  1619, 1452, 595, 940, 1107, 1964, 83, 108, 1939, 1132, 915, 620, 1427, 
  1644, 403, 364, 1683, 1388, 659, 876, 1171, 1900, 147, 236, 1811, 1260, 
  787, 748, 1299, 1772, 275, 492, 1555, 1516, 531, 1004, 1043, 2028, 19, 28, 
  2019, 1052, 995, 540, 1507, 1564, 483, 284, 1763, 1308, 739, 796, 1251, 
  1820, 227, 156, 1891, 1180, 867, 668, 1379, 1692, 355, 412, 1635, 1436, 
  611, 924, 1123, 1948, 99, 92, 1955, 1116, 931, 604, 1443, 1628, 419, 348, 
  1699, 1372, 675, 860, 1187, 1884, 163, 220, 1827, 1244, 803, 732, 1315, 
  1756, 291, 476, 1571, 1500, 547, 988, 1059, 2012, 35, 60, 1987, 1084, 963, 
  572, 1475, 1596, 451, 316, 1731, 1340, 707, 828, 1219, 1852, 195, 188, 
  1859, 1212, 835, 700, 1347, 1724, 323, 444, 1603, 1468, 579, 956, 1091, 
  1980, 67, 124, 1923, 1148, 899, 636, 1411, 1660, 387, 380, 1667, 1404, 643, 
  892, 1155, 1916, 131, 252, 1795, 1276, 771, 764, 1283, 1788, 259, 508, 
  1539, 1532, 515, 1020, 1027, 2044, 3, 2, 2045, 1026, 1021, 514, 1533, 1538, 
  509, 258, 1789, 1282, 765, 770, 1277, 1794, 253, 130, 1917, 1154, 893, 642, 
  1405, 1666, 381, 386, 1661, 1410, 637, 898, 1149, 1922, 125, 66, 1981, 
  1090, 957, 578, 1469, 1602, 445, 322, 1725, 1346, 701, 834, 1213, 1858, 
  189, 194, 1853, 1218, 829, 706, 1341, 1730, 317, 450, 1597, 1474, 573, 962, 
  1085, 1986, 61, 34, 2013, 1058, 989, 546, 1501, 1570, 477, 290, 1757, 1314, 
  733, 802, 1245, 1826, 221, 162, 1885, 1186, 861, 674, 1373, 1698, 349, 418, 
  1629, 1442, 605, 930, 1117, 1954, 93, 98, 1949, 1122, 925, 610, 1437, 1634, 
  413, 354, 1693, 1378, 669, 866, 1181, 1890, 157, 226, 1821, 1250, 797, 738, 
  1309, 1762, 285, 482, 1565, 1506, 541, 994, 1053, 2018, 29, 18, 2029, 1042, 
  1005, 530, 1517, 1554, 493, 274, 1773, 1298, 749, 786, 1261, 1810, 237, 
  146, 1901, 1170, 877, 658, 1389, 1682, 365, 402, 1645, 1426, 621, 914, 
  1133, 1938, 109, 82, 1965, 1106, 941, 594, 1453, 1618, 429, 338, 1709, 
  1362, 685, 850, 1197, 1874, 173, 210, 1837, 1234, 813, 722, 1325, 1746, 
  301, 466, 1581, 1490, 557, 978, 1069, 2002, 45, 50, 1997, 1074, 973, 562, 
  1485, 1586, 461, 306, 1741, 1330, 717, 818, 1229, 1842, 205, 178, 1869, 
  1202, 845, 690, 1357, 1714, 333, 434, 1613, 1458, 589, 946, 1101, 1970, 77, 
  114, 1933, 1138, 909, 626, 1421, 1650, 397, 370, 1677, 1394, 653, 882, 
  1165, 1906, 141, 242, 1805, 1266, 781, 754, 1293, 1778, 269, 498, 1549, 
  1522, 525, 1010, 1037, 2034, 13, 10, 2037, 1034, 1013, 522, 1525, 1546, 
  501, 266, 1781, 1290, 757, 778, 1269, 1802, 245, 138, 1909, 1162, 885, 650, 
  1397, 1674, 373, 394, 1653, 1418, 629, 906, 1141, 1930, 117, 74, 1973, 
  1098, 949, 586, 1461, 1610, 437, 330, 1717, 1354, 693, 842, 1205, 1866, 
  181, 202, 1845, 1226, 821, 714, 1333, 1738, 309, 458, 1589, 1482, 565, 970, 
  1077, 1994, 53, 42, 2005, 1066, 981, 554, 1493, 1578, 469, 298, 1749, 1322, 
  725, 810, 1237, 1834, 213, 170, 1877, 1194, 853, 682, 1365, 1706, 341, 426, 
  1621, 1450, 597, 938, 1109, 1962, 85, 106, 1941, 1130, 917, 618, 1429, 
  1642, 405, 362, 1685, 1386, 661, 874, 1173, 1898, 149, 234, 1813, 1258, 
  789, 746, 1301, 1770, 277, 490, 1557, 1514, 533, 1002, 1045, 2026, 21, 26, 
  2021, 1050, 997, 538, 1509, 1562, 485, 282, 1765, 1306, 741, 794, 1253, 
  1818, 229, 154, 1893, 1178, 869, 666, 1381, 1690, 357, 410, 1637, 1434, 
  613, 922, 1125, 1946, 101, 90, 1957, 1114, 933, 602, 1445, 1626, 421, 346, 
  1701, 1370, 677, 858, 1189, 1882, 165, 218, 1829, 1242, 805, 730, 1317, 
  1754, 293, 474, 1573, 1498, 549, 986, 1061, 2010, 37, 58, 1989, 1082, 965, 
  570, 1477, 1594, 453, 314, 1733, 1338, 709, 826, 1221, 1850, 197, 186, 
  1861, 1210, 837, 698, 1349, 1722, 325, 442, 1605, 1466, 581, 954, 1093, 
  1978, 69, 122, 1925, 1146, 901, 634, 1413, 1658, 389, 378, 1669, 1402, 645, 
  890, 1157, 1914, 133, 250, 1797, 1274, 773, 762, 1285, 1786, 261, 506, 
  1541, 1530, 517, 1018, 1029, 2042, 5, 6, 2041, 1030, 1017, 518, 1529, 1542, 
  505, 262, 1785, 1286, 761, 774, 1273, 1798, 249, 134, 1913, 1158, 889, 646, 
  1401, 1670, 377, 390, 1657, 1414, 633, 902, 1145, 1926, 121, 70, 1977, 
  1094, 953, 582, 1465, 1606, 441, 326, 1721, 1350, 697, 838, 1209, 1862, 
  185, 198, 1849, 1222, 825, 710, 1337, 1734, 313, 454, 1593, 1478, 569, 966, 
  1081, 1990, 57, 38, 2009, 1062, 985, 550, 1497, 1574, 473, 294, 1753, 1318, 
  729, 806, 1241, 1830, 217, 166, 1881, 1190, 857, 678, 1369, 1702, 345, 422, 
  1625, 1446, 601, 934, 1113, 1958, 89, 102, 1945, 1126, 921, 614, 1433, 
  1638, 409, 358, 1689, 1382, 665, 870, 1177, 1894, 153, 230, 1817, 1254, 
  793, 742, 1305, 1766, 281, 486, 1561, 1510, 537, 998, 1049, 2022, 25, 22, 
  2025, 1046, 1001, 534, 1513, 1558, 489, 278, 1769, 1302, 745, 790, 1257, 
  1814, 233, 150, 1897, 1174, 873, 662, 1385, 1686, 361, 406, 1641, 1430, 
  617, 918, 1129, 1942, 105, 86, 1961, 1110, 937, 598, 1449, 1622, 425, 342, 
  1705, 1366, 681, 854, 1193, 1878, 169, 214, 1833, 1238, 809, 726, 1321, 
  1750, 297, 470, 1577, 1494, 553, 982, 1065, 2006, 41, 54, 1993, 1078, 969, 
  566, 1481, 1590, 457, 310, 1737, 1334, 713, 822, 1225, 1846, 201, 182, 
  1865, 1206, 841, 694, 1353, 1718, 329, 438, 1609, 1462, 585, 950, 1097, 
  1974, 73, 118, 1929, 1142, 905, 630, 1417, 1654, 393, 374, 1673, 1398, 649, 
  886, 1161, 1910, 137, 246, 1801, 1270, 777, 758, 1289, 1782, 265, 502, 
  1545, 1526, 521, 1014, 1033, 2038, 9, 14, 2033, 1038, 1009, 526, 1521, 
  1550, 497, 270, 1777, 1294, 753, 782, 1265, 1806, 241, 142, 1905, 1166, 
  881, 654, 1393, 1678, 369, 398, 1649, 1422, 625, 910, 1137, 1934, 113, 78, 
  1969, 1102, 945, 590, 1457, 1614, 433, 334, 1713, 1358, 689, 846, 1201, 
  1870, 177, 206, 1841, 1230, 817, 718, 1329, 1742, 305, 462, 1585, 1486, 
  561, 974, 1073, 1998, 49, 46, 2001, 1070, 977, 558, 1489, 1582, 465, 302, 
  1745, 1326, 721, 814, 1233, 1838, 209, 174, 1873, 1198, 849, 686, 1361, 
  1710, 337, 430, 1617, 1454, 593, 942, 1105, 1966, 81, 110, 1937, 1134, 913, 
  622, 1425, 1646, 401, 366, 1681, 1390, 657, 878, 1169, 1902, 145, 238, 
  1809, 1262, 785, 750, 1297, 1774, 273, 494, 1553, 1518, 529, 1006, 1041, 
  2030, 17, 30, 2017, 1054, 993, 542, 1505, 1566, 481, 286, 1761, 1310, 737, 
  798, 1249, 1822, 225, 158, 1889, 1182, 865, 670, 1377, 1694, 353, 414, 
  1633, 1438, 609, 926, 1121, 1950, 97, 94, 1953, 1118, 929, 606, 1441, 1630, 
  417, 350, 1697, 1374, 673, 862, 1185, 1886, 161, 222, 1825, 1246, 801, 734, 
  1313, 1758, 289, 478, 1569, 1502, 545, 990, 1057, 2014, 33, 62, 1985, 1086, 
  961, 574, 1473, 1598, 449, 318, 1729, 1342, 705, 830, 1217, 1854, 193, 190, 
  1857, 1214, 833, 702, 1345, 1726, 321, 446, 1601, 1470, 577, 958, 1089, 
  1982, 65, 126, 1921, 1150, 897, 638, 1409, 1662, 385, 382, 1665, 1406, 641, 
  894, 1153, 1918, 129, 254, 1793, 1278, 769, 766, 1281, 1790, 257, 510, 
  1537, 1534, 513, 1022, 1025, 2046, 1};

static const int OUR4096[4096] = 
{0, 4095, 2048, 2047, 1024, 3071, 3072, 1023, 512, 3583, 2560, 1535, 1536, 
  2559, 3584, 511, 256, 3839, 2304, 1791, 1280, 2815, 3328, 767, 768, 3327, 
  2816, 1279, 1792, 2303, 3840, 255, 128, 3967, 2176, 1919, 1152, 2943, 3200, 
  895, 640, 3455, 2688, 1407, 1664, 2431, 3712, 383, 384, 3711, 2432, 1663, 
  1408, 2687, 3456, 639, 896, 3199, 2944, 1151, 1920, 2175, 3968, 127, 64, 
  4031, 2112, 1983, 1088, 3007, 3136, 959, 576, 3519, 2624, 1471, 1600, 2495, 
  3648, 447, 320, 3775, 2368, 1727, 1344, 2751, 3392, 703, 832, 3263, 2880, 
  1215, 1856, 2239, 3904, 191, 192, 3903, 2240, 1855, 1216, 2879, 3264, 831, 
  704, 3391, 2752, 1343, 1728, 2367, 3776, 319, 448, 3647, 2496, 1599, 1472, 
  2623, 3520, 575, 960, 3135, 3008, 1087, 1984, 2111, 4032, 63, 32, 4063, 
  2080, 2015, 1056, 3039, 3104, 991, 544, 3551, 2592, 1503, 1568, 2527, 3616, 
  479, 288, 3807, 2336, 1759, 1312, 2783, 3360, 735, 800, 3295, 2848, 1247, 
  1824, 2271, 3872, 223, 160, 3935, 2208, 1887, 1184, 2911, 3232, 863, 672, 
  3423, 2720, 1375, 1696, 2399, 3744, 351, 416, 3679, 2464, 1631, 1440, 2655, 
  3488, 607, 928, 3167, 2976, 1119, 1952, 2143, 4000, 95, 96, 3999, 2144, 
  1951, 1120, 2975, 3168, 927, 608, 3487, 2656, 1439, 1632, 2463, 3680, 415, 
  352, 3743, 2400, 1695, 1376, 2719, 3424, 671, 864, 3231, 2912, 1183, 1888, 
  2207, 3936, 159, 224, 3871, 2272, 1823, 1248, 2847, 3296, 799, 736, 3359, 
  2784, 1311, 1760, 2335, 3808, 287, 480, 3615, 2528, 1567, 1504, 2591, 3552, 
  543, 992, 3103, 3040, 1055, 2016, 2079, 4064, 31, 16, 4079, 2064, 2031, 
  1040, 3055, 3088, 1007, 528, 3567, 2576, 1519, 1552, 2543, 3600, 495, 272, 
  3823, 2320, 1775, 1296, 2799, 3344, 751, 784, 3311, 2832, 1263, 1808, 2287, 
  3856, 239, 144, 3951, 2192, 1903, 1168, 2927, 3216, 879, 656, 3439, 2704, 
  1391, 1680, 2415, 3728, 367, 400, 3695, 2448, 1647, 1424, 2671, 3472, 623, 
  912, 3183, 2960, 1135, 1936, 2159, 3984, 111, 80, 4015, 2128, 1967, 1104, 
  2991, 3152, 943, 592, 3503, 2640, 1455, 1616, 2479, 3664, 431, 336, 3759, 
  2384, 1711, 1360, 2735, 3408, 687, 848, 3247, 2896, 1199, 1872, 2223, 3920, 
  175, 208, 3887, 2256, 1839, 1232, 2863, 3280, 815, 720, 3375, 2768, 1327, 
  1744, 2351, 3792, 303, 464, 3631, 2512, 1583, 1488, 2607, 3536, 559, 976, 
  3119, 3024, 1071, 2000, 2095, 4048, 47, 48, 4047, 2096, 1999, 1072, 3023, 
  3120, 975, 560, 3535, 2608, 1487, 1584, 2511, 3632, 463, 304, 3791, 2352, 
  1743, 1328, 2767, 3376, 719, 816, 3279, 2864, 1231, 1840, 2255, 3888, 207, 
  176, 3919, 2224, 1871, 1200, 2895, 3248, 847, 688, 3407, 2736, 1359, 1712, 
  2383, 3760, 335, 432, 3663, 2480, 1615, 1456, 2639, 3504, 591, 944, 3151, 
  2992, 1103, 1968, 2127, 4016, 79, 112, 3983, 2160, 1935, 1136, 2959, 3184, 
  911, 624, 3471, 2672, 1423, 1648, 2447, 3696, 399, 368, 3727, 2416, 1679, 
  1392, 2703, 3440, 655, 880, 3215, 2928, 1167, 1904, 2191, 3952, 143, 240, 
  3855, 2288, 1807, 1264, 2831, 3312, 783, 752, 3343, 2800, 1295, 1776, 2319, 
  3824, 271, 496, 3599, 2544, 1551, 1520, 2575, 3568, 527, 1008, 3087, 3056, 
  1039, 2032, 2063, 4080, 15, 8, 4087, 2056, 2039, 1032, 3063, 3080, 1015, 
  520, 3575, 2568, 1527, 1544, 2551, 3592, 503, 264, 3831, 2312, 1783, 1288, 
  2807, 3336, 759, 776, 3319, 2824, 1271, 1800, 2295, 3848, 247, 136, 3959, 
  2184, 1911, 1160, 2935, 3208, 887, 648, 3447, 2696, 1399, 1672, 2423, 3720, 
  375, 392, 3703, 2440, 1655, 1416, 2679, 3464, 631, 904, 3191, 2952, 1143, 
  1928, 2167, 3976, 119, 72, 4023, 2120, 1975, 1096, 2999, 3144, 951, 584, 
  3511, 2632, 1463, 1608, 2487, 3656, 439, 328, 3767, 2376, 1719, 1352, 2743, 
  3400, 695, 840, 3255, 2888, 1207, 1864, 2231, 3912, 183, 200, 3895, 2248, 
  1847, 1224, 2871, 3272, 823, 712, 3383, 2760, 1335, 1736, 2359, 3784, 311,