table_xjj.c

同样是浮点型的fft算法一样为vc开发的

int tab_sqrt[101]=
{
 16384,16466,16547,16628,16708,16789,16868,16948,17027,17105,17184,17262,17339,17416,
 17493,17570,17646,17722,17798,17873,17948,18022,18097,18171,18244,18318,18391,18464,
 18536 ,18609,18681,18752,18824,18895,18966,19036 ,19107,19177,19247,19316,19386,19455,
 19524,19592,19661,19729,19797,19865,19932,19999,20066,20133,20200,20266,20332,20398,
 20464,20529,20594,20659,20724,20789,20853,20918,20982,21046,21109,21173,21236,21299,
 21362,21425,21487,21550,21612,21674,21736,21797,21859,21920,21981,22042,22103,22164,
 22224,22285,22345,22405,22465,22524,22584,22643,22702,22761,22820,22879,22938,22996,
 23054,23112,23170 
};

//分母1--2000

unsigned int daoshu[1001]={
0,65535,32768,21845,16384,13107,10922,9362,8192,7281,6553,
5957,5461,5041,4681,4369,4096,3855,3640,3449,3276,
3120,2978,2849,2730,2621,2520,2427,2340,2259,2184,
2114,2048,1985,1927,1872,1820,1771,1724,1680,1638,
1598,1560,1524,1489,1456,1424,1394,1365,1337,1310,
1285,1260,1236,1213,1191,1170,1149,1129,1110,1092,
1074,1057,1040,1024,1008,992,978,963,949,936,923,910,897,885,873,862,851,840,829,819,
809,799,789,780,771,762,753,744,736,728,720,712,704,697,689,682,675,668,661,655,
648,642,636,630,624,618,612,606,601,595,590,585,579,574,569,564,560,555,550,546,
541,537,532,528,524,520,516,512,508,504,500,496,492,489,485,481,478,474,471,468,
464,461,458,455,451,448,445,442,439,436,434,431,428,425,422,420,417,414,412,409,
407,404,402,399,397,394,392,390,387,385,383,381,378,376,374,372,370,368,366,364,
362,360,358,356,354,352,350,348,346,344,343,341,339,337,336,334,332,330,329,327,
326,324,322,321,319,318,316,315,313,312,310,309,307,306,304,303,302,300,299,297,
296,295,293,292,291,289,288,287,286,284,283,282,281,280,278,277,276,275,274,273,
271,270,269,268,267,266,265,264,263,262,261,260,259,258,257,256,255,254,253,252,
251,250,249,248,247,246,245,244,243,242,241,240,240,239,238,237,236,235,234,234,
233,232,231,230,229,229,228,227,226,225,225,224,223,222,222,221,220,219,219,218,
217,217,216,215,214,214,213,212,212,211,210,210,209,208,208,207,206,206,205,204,
204,203,202,202,201,201,200,199,199,198,197,197,196,196,195,195,194,193,193,192,
192,191,191,190,189,189,188,188,187,187,186,186,185,185,184,184,183,183,182,182,
181,181,180,180,179,179,178,178,177,177,176,176,175,175,174,174,173,173,172,172,
172,171,171,170,170,169,169,168,168,168,167,167,166,166,165,165,165,164,164,163,
163,163,162,162,161,161,161,160,160,159,159,159,158,158,157,157,157,156,156,156,
155,155,154,154,154,153,153,153,152,152,152,151,151,151,150,150,149,149,149,148,
148,148,147,147,147,146,146,146,145,145,145,144,144,144,144,143,143,143,142,142,
142,141,141,141,140,140,140,140,139,139,139,138,138,138,137,137,137,137,136,136,
136,135,135,135,135,134,134,134,134,133,133,133,132,132,132,132,131,131,131,131,
130,130,130,130,129,129,129,129,128,128,128,128,127,127,127,127,126,126,126,126,
125,125,125,125,124,124,124,124,123,123,123,123,122,122,122,122,122,121,121,121,
121,120,120,120,120,120,119,119,119,119,118,118,118,118,118,117,117,117,117,117,
116,116,116,116,115,115,115,115,115,114,114,114,114,114,113,113,113,113,113,112,
112,112,112,112,112,111,111,111,111,111,110,110,110,110,110,109,109,109,109,109,
109,108,108,108,108,108,107,107,107,107,107,107,106,106,106,106,106,106,105,105,
105,105,105,105,104,104,104,104,104,104,103,103,103,103,103,103,102,102,102,102,
102,102,101,101,101,101,101,101,100,100,100,100,100,100,100,99,99,99,99,99,
99,98,98,98,98,98,98,98,97,97,97,97,97,97,97,96,96,96,96,96,96,96,95,95,95,95,95,95,95,94,
94,94,94,94,94,94,94,93,93,93,93,93,93,93,92,92,92,92,92,92,92,92,91,91,91,91,91,91,91,91,
90,90,90,90,90,90,90,90,89,89,89,89,89,89,89,89,88,88,88,88,88,88,88,88,87,87,87,87,87,87,
87,87,87,86,86,86,86,86,86,86,86,86,85,85,85,85,85,85,85,85,85,84,84,84,84,84,84,84,84,84,
83,83,83,83,83,83,83,83,83,82,82,82,82,82,82,82,82,82,82,81,81,81,81,81,81,81,81,81,81,80,
80,80,80,80,80,80,80,80,80,79,79,79,79,79,79,79,79,79,79,78,78,78,78,78,78,78,78,78,78,78,
77,77,77,77,77,77,77,77,77,77,77,76,76,76,76,76,76,76,76,76,76,76,75,75,75,75,75,75,75,75,
75,75,75,74,74,74,74,74,74,74,74,74,74,74,74,73,73,73,73,73,73,73,73,73,73,73,73,72,72,72,
72,72,72,72,72,72,72,72,72,72,71,71,71,71,71,71,71,71,71,71,71,71,71,70,70,70,70,70,70,70,
70,70,70,70,70,70,69,69,69,69,69,69,69,69,69,69,69,69,69,68,68,68,68,68,68,68,68,68,68,68,
68,68,68,67,67,67,67,67,67,67,67,67,67,67,67,67,67,67,66,66,66,66,66,66,66,66,66,66,66,66,
66,66,65,65,65,65,65,65,65,65};

unsigned  int expi[640]={0,0,0,
40592,35280,31672,29020,26968,25320,23962,22817,21840,20988,20239,19575,18981,18446,17961,17518,17113,16740,
16396,16076,15779,15502,15243,15000,14771,14556,14353,14161,13979,13807,13643,13488,13340,13198,13063,12933,12810,12691,
12577,12468,12363,12263,12166,12072,11982,11895,11812,11731,11653,11577,11505,11434,11366,11299,11235,11173,11113,11054,
10997,10942,10889,10836,10786,10737,10689,10642,10597,10552,10510,10468,10427,10387,10348,10310,10273,10237,10202,10168,
10134,10101,10069,10038,10007,9977,9948,9919,9891,9864,9837,9811,9785,9760,9735,9711,9687,9664,9641,9619,
9597,9576,9555,9534,9514,9495,9475,9456,9437,9419,9401,9384,9366,9349,9333,9316,9300,9284,9269,9254,
9239,9224,9210,9196,9182,9168,9154,9141,9128,9115,9103,9090,9078,9066,9054,9043,9031,9020,9009,8998,
8988,8977,8967,8957,8947,8937,8927,8918,8908,8899,8890,8881,8872,8863,8855,8846,8838,8830,8822,8814,
8806,8798,8791,8783,8776,8768,8761,8754,8747,8740,8734,8727,8720,8714,8708,8701,8695,8689,8683,8677,
8671,8665,8660,8654,8648,8643,8638,8632,8627,8622,8617,8612,8607,8602,8597,8592,8587,8583,8578,8574,
8569,8565,8561,8556,8552,8548,8544,8540,8536,8532,8528,8524,8520,8516,8513,8509,8505,8502,8498,8495,
8491,8488,8485,8481,8478,8475,8472,8468,8465,8462,8459,8456,8453,8450,8447,8445,8442,8439,8436,8434,
8431,8428,8426,8423,8421,8418,8415,8413,8411,8408,8406,8403,8401,8399,8397,8394,8392,8390,8388,8386,
8384,8382,8379,8377,8375,8373,8372,8370,8368,8366,8364,8362,8360,8358,8357,8355,8353,8351,8350,8348,
8346,8345,8343,8341,8340,8338,8337,8335,8334,8332,8331,8329,8328,8326,8325,8324,8322,8321,8319,8318,
8317,8315,8314,8313,8312,8310,8309,8308,8307,8306,8304,8303,8302,8301,8300,8299,8297,8296,8295,8294,
8293,8292,8291,8290,8289,8288,8287,8286,8285,8284,8283,8282,8281,8280,8280,8279,8278,8277,8276,8275,
8274,8273,8273,8272,8271,8270,8269,8269,8268,8267,8266,8265,8265,8264,8263,8263,8262,8261,8260,8260,
8259,8258,8258,8257,8256,8256,8255,8254,8254,8253,8253,8252,8251,8251,8250,8250,8249,8248,8248,8247,
8247,8246,8246,8245,8245,8244,8244,8243,8243,8242,8242,8241,8241,8240,8240,8239,8239,8238,8238,8237,
8237,8236,8236,8235,8235,8235,8234,8234,8233,8233,8233,8232,8232,8231,8231,8231,8230,8230,8229,8229,
8229,8228,8228,8228,8227,8227,8227,8226,8226,8226,8225,8225,8225,8224,8224,8224,8223,8223,8223,8222,
8222,8222,8222,8221,8221,8221,8220,8220,8220,8220,8219,8219,8219,8218,8218,8218,8218,8217,8217,8217,
8217,8216,8216,8216,8216,8216,8215,8215,8215,8215,8214,8214,8214,8214,8214,8213,8213,8213,8213,8212,
8212,8212,8212,8212,8212,8211,8211,8211,8211,8211,8210,8210,8210,8210,8210,8210,8209,8209,8209,8209,
8209,8209,8208,8208,8208,8208,8208,8208,8207,8207,8207,8207,8207,8207,8207,8206,8206,8206,8206,8206,
8206,8206,8205,8205,8205,8205,8205,8205,8205,8205,8204,8204,8204,8204,8204,8204,8204,8204,8203,8203,
8203,8203,8203,8203,8203,8203,8203,8203,8202,8202,8202,8202,8202,8202,8202,8202,8202,8202,8201,8201,
8201,8201,8201,8201,8201,8201,8201,8201,8201,8200,8200,8200,8200,8200,8200,8200,8200,8200,8200,8200,
8200,8200,8199,8199,8199,8199,8199,8199,8199,8199,8199,8199,8199,8199,8199,8199,8198,8198,8198,8198,
8198,8198,8198,8198,8198,8198,8198,8198,8198,8198,8198,8198,8198,8197,8197,8197,8197,8197,8197,8197,
8197,8197,8197,8197,8197,8197,8197,8197,8197,8197,8197,8197,8197,8196,8196,8196,8196,8196,8196};