trigtbl.h
1 
2 /*
3 ** Please only distribute this with it's associated FHT routine.
4 ** This algorithm is apparently patented(!) and the code copyrighted.
5 ** See the comment with the fht routine for more info.
6 ** -Thanks,
7 ** Ron Mayer
8 */
9 
10 #ifdef __cplusplus
11 extern "C" {
12 #endif
13 
14 #ifdef REAL
15 #else
16 #define REAL double
17 #endif
18 
19 #ifdef GOOD_TRIG
20 #else
21 #define FAST_TRIG
22 #endif
23 
24 #if defined(GOOD_TRIG)
25 #define FHT_SWAP(a,b,t) {(t)=(a);(a)=(b);(b)=(t);}
26 #define TRIG_VARS \
27  int t_lam=0;
28 #define TRIG_INIT(k,c,s) \
29  { \
30  int i; \
31  for (i=2 ; i<=k ; i++) \
32  {coswrk[i]=costab[i];sinwrk[i]=sintab[i];} \
33  t_lam = 0; \
34  c = 1; \
35  s = 0; \
36  }
37 #define TRIG_NEXT(k,c,s) \
38  { \
39  int i,j; \
40  (t_lam)++; \
41  for (i=0 ; !((1<<i)&t_lam) ; i++); \
42  i = k-i; \
43  s = sinwrk[i]; \
44  c = coswrk[i]; \
45  if (i>1) \
46  { \
47  for (j=k-i+2 ; (1<<j)&t_lam ; j++); \
48  j = k - j; \
49  sinwrk[i] = halsec[i] * (sinwrk[i-1] + sinwrk[j]); \
50  coswrk[i] = halsec[i] * (coswrk[i-1] + coswrk[j]); \
51  } \
52  }
53 #define TRIG_RESET(k,c,s)
54 #endif
55 
56 #if defined(FAST_TRIG)
57 #define TRIG_VARS \
58  REAL t_c,t_s;
59 #define TRIG_INIT(k,c,s) \
60  { \
61  t_c = costab[k]; \
62  t_s = sintab[k]; \
63  c = 1; \
64  s = 0; \
65  }
66 #define TRIG_NEXT(k,c,s) \
67  { \
68  REAL t = c; \
69  c = t*t_c - s*t_s; \
70  s = t*t_s + s*t_c; \
71  }
72 #define TRIG_RESET(k,c,s)
73 #endif
74 
75 static REAL halsec[20]=
76  {
77  0,
78  0,
79  .54119610014619698439972320536638942006107206337801,
80  .50979557910415916894193980398784391368261849190893,
81  .50241928618815570551167011928012092247859337193963,
82  .50060299823519630134550410676638239611758632599591,
83  .50015063602065098821477101271097658495974913010340,
84  .50003765191554772296778139077905492847503165398345,
85  .50000941253588775676512870469186533538523133757983,
86  .50000235310628608051401267171204408939326297376426,
87  .50000058827484117879868526730916804925780637276181,
88  .50000014706860214875463798283871198206179118093251,
89  .50000003676714377807315864400643020315103490883972,
90  .50000000919178552207366560348853455333939112569380,
91  .50000000229794635411562887767906868558991922348920,
92  .50000000057448658687873302235147272458812263401372
93  };
94 static REAL costab[20]=
95  {
96  .00000000000000000000000000000000000000000000000000,
97  .70710678118654752440084436210484903928483593768847,
98  .92387953251128675612818318939678828682241662586364,
99  .98078528040323044912618223613423903697393373089333,
100  .99518472667219688624483695310947992157547486872985,
101  .99879545620517239271477160475910069444320361470461,
102  .99969881869620422011576564966617219685006108125772,
103  .99992470183914454092164649119638322435060646880221,
104  .99998117528260114265699043772856771617391725094433,
105  .99999529380957617151158012570011989955298763362218,
106  .99999882345170190992902571017152601904826792288976,
107  .99999970586288221916022821773876567711626389934930,
108  .99999992646571785114473148070738785694820115568892,
109  .99999998161642929380834691540290971450507605124278,
110  .99999999540410731289097193313960614895889430318945,
111  .99999999885102682756267330779455410840053741619428
112  };
113 static REAL sintab[20]=
114  {
115  1.0000000000000000000000000000000000000000000000000,
116  .70710678118654752440084436210484903928483593768846,
117  .38268343236508977172845998403039886676134456248561,
118  .19509032201612826784828486847702224092769161775195,
119  .09801714032956060199419556388864184586113667316749,
120  .04906767432741801425495497694268265831474536302574,
121  .02454122852291228803173452945928292506546611923944,
122  .01227153828571992607940826195100321214037231959176,
123  .00613588464915447535964023459037258091705788631738,
124  .00306795676296597627014536549091984251894461021344,
125  .00153398018628476561230369715026407907995486457522,
126  .00076699031874270452693856835794857664314091945205,
127  .00038349518757139558907246168118138126339502603495,
128  .00019174759731070330743990956198900093346887403385,
129  .00009587379909597734587051721097647635118706561284,
130  .00004793689960306688454900399049465887274686668768
131  };
132 static REAL coswrk[20]=
133  {
134  .00000000000000000000000000000000000000000000000000,
135  .70710678118654752440084436210484903928483593768847,
136  .92387953251128675612818318939678828682241662586364,
137  .98078528040323044912618223613423903697393373089333,
138  .99518472667219688624483695310947992157547486872985,
139  .99879545620517239271477160475910069444320361470461,
140  .99969881869620422011576564966617219685006108125772,
141  .99992470183914454092164649119638322435060646880221,
142  .99998117528260114265699043772856771617391725094433,
143  .99999529380957617151158012570011989955298763362218,
144  .99999882345170190992902571017152601904826792288976,
145  .99999970586288221916022821773876567711626389934930,
146  .99999992646571785114473148070738785694820115568892,
147  .99999998161642929380834691540290971450507605124278,
148  .99999999540410731289097193313960614895889430318945,
149  .99999999885102682756267330779455410840053741619428
150  };
151 static REAL sinwrk[20]=
152  {
153  1.0000000000000000000000000000000000000000000000000,
154  .70710678118654752440084436210484903928483593768846,
155  .38268343236508977172845998403039886676134456248561,
156  .19509032201612826784828486847702224092769161775195,
157  .09801714032956060199419556388864184586113667316749,
158  .04906767432741801425495497694268265831474536302574,
159  .02454122852291228803173452945928292506546611923944,
160  .01227153828571992607940826195100321214037231959176,
161  .00613588464915447535964023459037258091705788631738,
162  .00306795676296597627014536549091984251894461021344,
163  .00153398018628476561230369715026407907995486457522,
164  .00076699031874270452693856835794857664314091945205,
165  .00038349518757139558907246168118138126339502603495,
166  .00019174759731070330743990956198900093346887403385,
167  .00009587379909597734587051721097647635118706561284,
168  .00004793689960306688454900399049465887274686668768
169  };
170 
171 #ifdef __cplusplus
172 }
173 #endif