Logo Coherent WaveBurst  
Reference Guide 		Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
Meyer.cc
Go to the documentation of this file.
1 // Wavelet Analysis Tool
2 //--------------------------------------------------------------------
3 // Implementation of
4 // Meyer wavelets using Fast Wavelet Transform
5 // References:
6 //--------------------------------------------------------------------
7 //$Id: Meyer.hh,v 0.2 2001/08/06 19:37:00 klimenko Exp $
8 
9 #define MEYER_CC
10 
11 #include "Meyer.hh"
12 #include <iostream>
13 #include "stdio.h"
14 
15 //namespace datacondAPI {
16 //namespace wat {
17 
18 extern const double mey[16]=
19 { 7.4458559231880628e-01, 4.4459300275757724e-01, -3.5087555656258346e-02,
20  -1.3284520043622938e-01, 3.0655091960824263e-02, 6.3739024322801596e-02,
21  -2.4348745906078023e-02, -3.2130793990211758e-02, 1.7423434103729693e-02,
22  1.5270015130934803e-02, -1.1061496392513451e-02, -6.3877183184971563e-03,
23  6.0458140973233040e-03, 2.2025341009110021e-03, -2.7046721246437250e-03,
24  -6.0115023435160925e-04};
25 
26 extern const double mey004[32]=
27 { 7.43750353e-01, 4.44094660e-01, -3.50482247e-02, -1.32696489e-01,
28  3.06212119e-02, 6.36672446e-02, -2.43217736e-02, -3.20940567e-02,
29  1.74038765e-02, 1.52509309e-02, -1.10446159e-02, -6.38673997e-03,
30  6.04547967e-03, 2.19476459e-03, -2.70215388e-03, -5.78163043e-04,
31  8.59511073e-04, 1.61178251e-04, -9.35363005e-05, -1.39924765e-04,
32  -7.55811503e-05, 1.49736854e-04, 2.44738777e-05, -1.02790685e-04,
33  3.73859490e-05, 3.23295640e-05, -4.64169099e-05, 1.50170027e-05,
34  2.01182002e-05, -2.44497410e-05, 6.16393792e-06, 1.04921553e-05};
35 
36 extern const double mey008[48]=
37 { 7.33609639e-01, 4.46956642e-01, -2.58537823e-02, -1.40705138e-01,
38  2.39953217e-02, 7.52666880e-02, -2.11749933e-02, -4.52098615e-02,
39  1.77478786e-02, 2.78607152e-02, -1.41068440e-02, -1.69993874e-02,
40  1.06121463e-02, 1.00922602e-02, -7.53699278e-03, -5.78259209e-03,
41  5.03950041e-03, 3.19710122e-03, -3.16307359e-03, -1.72272472e-03,
42  1.85945070e-03, 9.25826653e-04, -1.02391321e-03, -5.12999856e-04,
43  5.31379384e-04, 3.00631034e-04, -2.64764748e-04, -1.85201368e-04,
44  1.31470288e-04, 1.15431643e-04, -6.82864331e-05, -6.94283079e-05,
45  3.79892424e-05, 3.88121845e-05, -2.19363081e-05, -1.97440301e-05,
46  1.23471763e-05, 9.18041839e-06, -6.45040920e-06, -4.13481950e-06,
47  3.11352396e-06, 2.05800311e-06, -1.49196971e-06, -1.22615994e-06,
48  8.12801410e-07, 7.67372574e-07, -5.15136932e-07, -4.06233507e-07};
49 
50 extern const double mey016[64]=
51 { 7.26060766e-01, 4.48510741e-01, -1.87063583e-02, -1.45182691e-01,
52  1.79824715e-02, 8.21488422e-02, -1.68369560e-02, -5.37407521e-02,
53  1.53534264e-02, 3.71819869e-02, -1.36345964e-02, -2.62900015e-02,
54  1.17907408e-02, 1.86830184e-02, -9.92827064e-03, -1.32245702e-02,
55  8.14003237e-03, 9.27736602e-03, -6.49849436e-03, -6.43454314e-03,
56  5.05236309e-03, 4.41002788e-03, -3.82653735e-03, -2.99022990e-03,
57  2.82480333e-03, 2.01133358e-03, -2.03436735e-03, -1.34748651e-03,
58  1.43122726e-03, 9.03470490e-04, -9.85466744e-04, -6.09099247e-04,
59  6.65769880e-04, 4.14298193e-04, -4.42743822e-04, -2.84604394e-04,
60  2.90941691e-04, 1.97138048e-04, -1.89729736e-04, -1.37177063e-04,
61  1.23292062e-04, 9.54358756e-05, -8.01004803e-05, -6.60692637e-05,
62  5.21278478e-05, 4.53290749e-05, -3.40029351e-05, -3.07304146e-05,
63  2.22340207e-05, 2.05603112e-05, -1.45771375e-05, -1.35893549e-05,
64  9.58175075e-06, 8.90067120e-06, -6.30288149e-06, -5.79549494e-06,
65  4.13268817e-06, 3.75088503e-06, -2.69342158e-06, -2.40111942e-06,
66  1.75266393e-06, 1.51368900e-06, -1.15405288e-06, -9.47428710e-07};
67 
68 extern const double mey032[96]=
69 { 7.20585503e-01, 4.49322176e-01, -1.33876233e-02, -1.47564003e-01,
70  1.31180419e-02, 8.59468746e-02, -1.26808775e-02, -5.87187863e-02,
71  1.20935038e-02, 4.30446119e-02, -1.13786310e-02, -3.27126251e-02,
72  1.05628353e-02, 2.53419265e-02, -9.67489952e-03, -1.98236902e-02,
73  8.74412541e-03, 1.55687636e-02, -7.79876613e-03, -1.22309446e-02,
74  6.86468095e-03, 9.58904815e-03, -5.96427568e-03, -7.49122252e-03,
75  5.11576757e-03, 5.82666639e-03, -4.33280026e-03, -4.51048065e-03,
76  3.62441235e-03, 3.47522836e-03, -2.99532338e-03, -2.66608068e-03,
77  2.44645878e-03, 2.03792774e-03, -1.97561400e-03, -1.55355878e-03,
78  1.57816957e-03, 1.18240002e-03, -1.24780015e-03, -8.99529934e-04,
79  9.77145013e-04, 6.84842091e-04, -7.58411779e-04, -5.22310303e-04,
80  5.83878332e-04, 3.99340552e-04, -4.46261029e-04, -3.06191991e-04,
81  3.38941838e-04, 2.35446834e-04, -2.56079543e-04, -1.81522688e-04,
82  1.92647355e-04, 1.40242548e-04, -1.44429583e-04, -1.08485956e-04,
83  1.07985946e-04, 8.39294256e-05, -8.05779149e-05, -6.48578157e-05,
84  6.00598599e-05, 5.00147104e-05, -4.47575881e-05, -3.84709859e-05,
85  3.33650554e-05, 2.95145229e-05, -2.48751954e-05, -2.25766867e-05,
86  1.85348434e-05, 1.72006556e-05, -1.38002625e-05, -1.30344747e-05,
87  1.02805058e-05, 9.82123371e-06, -7.67887314e-06, -7.37168717e-06,
88  5.75432885e-06, 5.52911299e-06, -4.31302644e-06, -4.14866747e-06,
89  3.21670349e-06, 3.10218711e-06, -2.38316834e-06, -2.29613316e-06,
90  1.76612899e-06, 1.67878805e-06, -1.32540443e-06, -1.22471389e-06,
91  1.01052828e-06, 9.08116537e-07, -7.69114071e-07, -6.88574605e-07,
92  5.67889970e-07, 5.20825336e-07, -4.02604215e-07, -3.76740644e-07};
93 
94 extern const double mey064[128]=
95 { 7.16664739e-01, 4.49737013e-01, -9.52509311e-03, -1.48794259e-01,
96  9.42718353e-03, 8.79506343e-02, -9.26626565e-03, -6.14289433e-02,
97  9.04568038e-03, 4.63727568e-02, -8.76996656e-03, -3.65540285e-02,
98  8.44469228e-03, 2.95812182e-02, -8.07624332e-03, -2.43407973e-02,
99  7.67160099e-03, 2.02446765e-02, -7.23813107e-03, -1.69529251e-02,
100  6.78338198e-03, 1.42549718e-02, -6.31487811e-03, -1.20128486e-02,
101  5.83990528e-03, 1.01317354e-02, -5.36530749e-03, -8.54364745e-03,
102  4.89732421e-03, 7.19792811e-03, -4.44148387e-03, -6.05549731e-03,
103  4.00254383e-03, 5.08527746e-03, -3.58445343e-03, -4.26192058e-03,
104  3.19032657e-03, 3.56431865e-03, -2.82243293e-03, -2.97459250e-03,
105  2.48222706e-03, 2.47739105e-03, -2.17042247e-03, -2.05941104e-03,
106  1.88709354e-03, 1.70907625e-03, -1.63177647e-03, -1.41631659e-03,
107  1.40355229e-03, 1.17239254e-03, -1.20111920e-03, -9.69732697e-04,
108  1.02287367e-03, 8.01783161e-04, -8.67009056e-04, -6.62884201e-04,
109  7.31617418e-04, 5.48181040e-04, -6.14769429e-04, -4.53554313e-04,
110  5.14559732e-04, 3.75545294e-04, -4.29129354e-04, -3.11263665e-04,
111  3.56689096e-04, 2.58289569e-04, -2.95556490e-04, -2.14593685e-04,
112  2.44195619e-04, 1.78487296e-04, -2.01237488e-04, -1.48590646e-04,
113  1.65470505e-04, 1.23795972e-04, -1.35814182e-04, -1.03213263e-04,
114  1.11300869e-04, 8.61101801e-05, -9.10784432e-05, -7.18692981e-05,
115  7.44231322e-05, 5.99739223e-05, -6.07397066e-05, -5.00101480e-05,
116  4.95379001e-05, 4.16610140e-05, -4.03973174e-05, -3.46804093e-05,
117  3.29447468e-05, 2.88579911e-05, -2.68558874e-05, -2.39982035e-05,
118  2.18698321e-05, 1.99247507e-05, -1.77927861e-05, -1.64983720e-05,
119  1.44792262e-05, 1.36240347e-05, -1.18022353e-05, -1.12355047e-05,
120  9.63651184e-06, 9.26887294e-06, -7.86574972e-06, -7.64844191e-06,
121  6.40251358e-06, 6.29662555e-06, -5.19694529e-06, -5.15597698e-06,
122  4.22242785e-06, 4.19969775e-06, -3.44991009e-06, -3.41878746e-06,
123  2.83439116e-06, 2.79757727e-06, -2.32528820e-06, -2.30118524e-06,
124  1.88884962e-06, 1.88665321e-06, -1.51899775e-06, -1.52596025e-06,
125  1.22477622e-06, 1.21732860e-06, -1.00615445e-06, -9.73026670e-07,
126  8.41733183e-07, 7.95446067e-07, -7.00115852e-07, -6.65015290e-07};
127 
128 extern const double mey128[192]=
129 { 7.13874857e-01, 4.49946782e-01, -6.75634361e-03, -1.49419865e-01,
130  6.72126119e-03, 8.89810886e-02, -6.66318178e-03, -6.28463092e-02,
131  6.58271363e-03, 4.81526973e-02, -6.48072078e-03, -3.86665398e-02,
132  6.35829589e-03, 3.19915158e-02, -6.21671757e-03, -2.70103180e-02,
133  6.05741630e-03, 2.31321538e-02, -5.88196124e-03, -2.00154955e-02,
134  5.69205682e-03, 1.74492695e-02, -5.48952611e-03, -1.52960492e-02,
135  5.27626983e-03, 1.34625138e-02, -5.05421361e-03, -1.18830206e-02,
136  4.82526746e-03, 1.05099617e-02, -4.59131003e-03, -9.30786174e-03,
137  4.35418637e-03, 8.24964702e-03, -4.11569602e-03, -7.31421627e-03,
138  3.87756010e-03, 6.48480242e-03, -3.64137952e-03, -5.74782559e-03,
139  3.40860795e-03, 5.09207359e-03, -3.18055166e-03, -4.50812305e-03,
140  2.95838430e-03, 3.98794154e-03, -2.74315314e-03, -3.52461175e-03,
141  2.53576474e-03, 3.11212215e-03, -2.33696174e-03, -2.74519036e-03,
142  2.14731404e-03, 2.41911564e-03, -1.96723589e-03, -2.12967297e-03,
143  1.79701675e-03, 1.87305269e-03, -1.63684203e-03, -1.64582811e-03,
144  1.48679163e-03, 1.44492316e-03, -1.34682759e-03, -1.26756542e-03,
145  1.21679416e-03, 1.11123371e-03, -1.09644184e-03, -9.73622283e-04,
146  9.85463097e-04, 8.52632074e-04, -8.83516198e-04, -7.46376303e-04,
147  7.90224945e-04, 6.53176176e-04, -7.05166027e-04, -5.71534445e-04,
148  6.27867363e-04, 5.00098286e-04, -5.57829147e-04, -4.37634821e-04,
149  4.94555708e-04, 3.83030948e-04, -4.37574603e-04, -3.35305600e-04,
150  3.86431130e-04, 2.93610816e-04, -3.40670064e-04, -2.57209815e-04,
151  2.99828230e-04, 2.25443881e-04, -2.63449743e-04, -1.97711592e-04,
152  2.31112201e-04, 1.73472226e-04, -2.02440347e-04, -1.52261530e-04,
153  1.77095480e-04, 1.33696306e-04, -1.54752496e-04, -1.17456008e-04,
154  1.35088212e-04, 1.03253143e-04, -1.17792828e-04, -9.08160149e-05,
155  1.02592849e-04, 7.98955829e-05, -8.92619378e-05, -7.02846186e-05,
156  7.76079954e-05, 6.18255623e-05, -6.74483122e-05, -5.43952605e-05,
157  5.85964038e-05, 4.78783358e-05, -5.08723664e-05, -4.21527247e-05,
158  4.41250033e-05, 3.70991295e-05, -3.82421809e-05, -3.26225654e-05,
159  3.31376531e-05, 2.86623967e-05, -2.87261562e-05, -2.51790476e-05,
160  2.49103579e-05, 2.21291417e-05, -2.15914552e-05, -1.94526151e-05,
161  1.86916430e-05, 1.70835602e-05, -1.61648846e-05, -1.49729930e-05,
162  1.39841985e-05, 1.30999537e-05, -1.21172445e-05, -1.14591373e-05,
163  1.05137757e-05, 1.00368267e-05, -9.11673804e-06, -8.79868318e-06,
164  7.88522928e-06, 7.70116876e-06, -6.80573957e-06, -6.71480408e-06,
165  5.87988418e-06, 5.83568375e-06, -5.10040617e-06, -5.07345765e-06,
166  4.43901459e-06, 4.42755830e-06, -3.85783858e-06, -3.87523967e-06,
167  3.33270839e-06, 3.38320811e-06, -2.86468813e-06, -2.93105633e-06,
168  2.46807056e-06, 2.52292752e-06, -2.14661427e-06, -2.17562281e-06,
169  1.88159042e-06, 1.89493599e-06, -1.64342181e-06, -1.66378420e-06,
170  1.41512696e-06, 1.45391789e-06, -1.20399650e-06, -1.24942327e-06,
171  1.02971383e-06, 1.05844635e-06, -9.00702688e-07, -9.01352903e-07,
172  8.02270160e-07, 7.87109450e-07, -7.08329725e-07, -7.01454992e-07,
173  6.04918371e-07, 6.18647449e-07, -5.01936480e-07, -5.24998941e-07,
174  4.21323580e-07, 4.30629196e-07, -3.73453754e-07, -3.57652076e-07,
175  3.45320152e-07, 3.16580472e-07, -3.12293890e-07, -2.94519438e-07,
176  2.61672192e-07, 2.66932032e-07, -2.04445070e-07, -2.21192051e-07};
177 
178 extern const double mey256[256]=
179 { 7.11895915e-01, 4.50052264e-01, -4.78497187e-03, -1.49735366e-01,
180  4.77249670e-03, 8.95037939e-02, -4.75175254e-03, -6.35715771e-02,
181  4.72283975e-03, 4.90741318e-02, -4.68592571e-03, -3.97760905e-02,
182  4.64123155e-03, 3.32796024e-02, -4.58900682e-03, -2.84659765e-02,
183  4.52951548e-03, 2.47432011e-02, -4.46304525e-03, -2.17687091e-02,
184  4.38992835e-03, 1.93305783e-02, -4.31055029e-03, -1.72907313e-02,
185  4.22533483e-03, 1.55553974e-02, -4.13471695e-03, -1.40586841e-02,
186  4.03912748e-03, 1.27529294e-02, -3.93900127e-03, -1.16027913e-02,
187  3.83479696e-03, 1.05815055e-02, -3.72700512e-03, -9.66844238e-03,
188  3.61613270e-03, 8.84745277e-03, -3.50267588e-03, -8.10570245e-03,
189  3.38710476e-03, 7.43283362e-03, -3.26987177e-03, -6.82036590e-03,
190  3.15143203e-03, 6.26127812e-03, -3.03225220e-03, -5.74971233e-03,
191  2.91279589e-03, 5.28074500e-03, -2.79349761e-03, -4.85019185e-03,
192  2.67474873e-03, 4.45444318e-03, -2.55690733e-03, -4.09034238e-03,
193  2.44032007e-03, 3.75511180e-03, -2.32533262e-03, -3.44630841e-03,
194  2.21227676e-03, 3.16178175e-03, -2.10144605e-03, -2.89961896e-03,
195  1.99308348e-03, 2.65808659e-03, -1.88739305e-03, -2.43559083e-03,
196  1.78456332e-03, 2.23066657e-03, -1.68477945e-03, -2.04198243e-03,
197  1.58821189e-03, 1.86833741e-03, -1.49499352e-03, -1.70863655e-03,
198  1.40520867e-03, 1.56185713e-03, -1.31890598e-03, -1.42702822e-03,
199  1.23612311e-03, 1.30323532e-03, -1.15689982e-03, -1.18963786e-03,
200  1.08126754e-03, 1.08547586e-03, -1.00922727e-03, -9.90053688e-04,
201  9.40739287e-04, 9.02712742e-04, -8.75736458e-04, -8.22816264e-04,
202  8.14149362e-04, 7.49758197e-04, -7.55919615e-04, -6.82984093e-04,
203  7.00989599e-04, 6.22000499e-04, -6.49280388e-04, -5.66360968e-04,
204  6.00681390e-04, 5.15640447e-04, -5.55063515e-04, -4.69421588e-04,
205  5.12304059e-04, 4.27304738e-04, -4.72299723e-04, -3.88929801e-04,
206  4.34955986e-04, 3.53986405e-04, -4.00164607e-04, -3.22200558e-04,
207  3.67792818e-04, 2.93309586e-04, -3.37696007e-04, -2.67048920e-04,
208  3.09742089e-04, 2.43162506e-04, -2.83824007e-04, -2.21425047e-04,
209  2.59848581e-04, 2.01652528e-04, -2.37713474e-04, -1.83689204e-04,
210  2.17295879e-04, 1.67382873e-04, -1.98464686e-04, -1.52571976e-04,
211  1.81104372e-04, 1.39096329e-04, -1.65127035e-04, -1.26819686e-04,
212  1.50460798e-04, 1.15640574e-04, -1.37026367e-04, -1.05479612e-04,
213  1.24725331e-04, 9.62550897e-05, -1.13451977e-04, -8.78703983e-05,
214  1.03116857e-04, 8.02250687e-05, -9.36585297e-05, -7.32376525e-05,
215  8.50317076e-05, 6.68568635e-05, -7.71835855e-05, -6.10491963e-05,
216  7.00419362e-05, 5.57748051e-05, -6.35267576e-05, -5.09752111e-05,
217  5.75736906e-05, 4.65846222e-05, -5.21456430e-05, -4.25530751e-05,
218  4.72208342e-05, 3.88578298e-05, -4.27690495e-05, -3.54912276e-05,
219  3.87396764e-05, 3.24367969e-05, -3.50733108e-05, -2.96571742e-05,
220  3.17251496e-05, 2.71056237e-05, -2.86766020e-05, -2.47493665e-05,
221  2.59233336e-05, 2.25811494e-05, -2.34515310e-05, -2.06072631e-05,
222  2.12259555e-05, 1.88237971e-05, -1.92015749e-05, -1.72046971e-05,
223  1.73469870e-05, 1.57134092e-05, -1.56560663e-05, -1.43263233e-05,
224  1.41360474e-05, 1.30444452e-05, -1.27838306e-05, -1.18815095e-05,
225  1.15740798e-05, 1.08403188e-05, -1.04708976e-05, -9.90087789e-06,
226  9.45129366e-06, 9.03210675e-06, -8.51687574e-06, -8.21534571e-06,
227  7.68197656e-06, 7.45608240e-06, -6.95000216e-06, -6.77211429e-06,
228  6.30157083e-06, 6.16993176e-06, -5.70622823e-06, -5.63289129e-06,
229  5.14595364e-06, 5.13296307e-06, -4.62688686e-06, -4.65426709e-06,
230  4.16748961e-06, 4.20482738e-06, -3.77492600e-06, -3.80475880e-06,
231  3.43323180e-06, 3.46267165e-06, -3.11505790e-06, -3.16386391e-06,
232  2.80520299e-06, 2.88208567e-06, -2.51236489e-06, -2.60309036e-06,
233  2.25732477e-06, 2.33640219e-06, -2.04934863e-06, -2.10351500e-06,
234  1.87437601e-06, 1.91430744e-06, -1.70678134e-06, -1.75524474e-06,
235  1.53292300e-06, 1.60115165e-06, -1.36290960e-06, -1.43877127e-06,
236  1.21879802e-06, 1.27853951e-06, -1.11100740e-06, -1.14279008e-06,
237  1.02651944e-06, 1.04217548e-06, -9.40650171e-07, -9.63873785e-07,
238  8.40608522e-07, 8.83366419e-07, -7.37271313e-07, -7.88001215e-07,
239  6.53389032e-07, 6.88770812e-07, -6.00007052e-07, -6.08521347e-07,
240  5.64672395e-07, 5.58376709e-07, -5.23211431e-07, -5.25948504e-07,
241  4.63293760e-07, 4.87116590e-07, -3.96211976e-07, -4.29594843e-07,
242  3.45091793e-07, 3.64711991e-07, -3.21317204e-07, -3.15622573e-07};
243 
244 extern const double mey384[384]=
245 { 7.11017584e-01, 4.50087409e-01, -3.90864520e-03, -1.49840731e-01,
246  3.90206191e-03, 8.96790649e-02, -3.89084227e-03, -6.38159860e-02,
247  3.87490944e-03, 4.93863601e-02, -3.85443835e-03, -4.01543853e-02,
248  3.82973748e-03, 3.37220018e-02, -3.80101268e-03, -2.89704313e-02,
249  3.76824883e-03, 2.53074870e-02, -3.73132705e-03, -2.23902353e-02,
250  3.69025970e-03, 2.00063033e-02, -3.64530741e-03, -1.80172933e-02,
251  3.59686042e-03, 1.63293518e-02, -3.54520195e-03, -1.48766360e-02,
252  3.49038937e-03, 1.36114551e-02, -3.43237106e-03, -1.24982537e-02,
253  3.37122106e-03, 1.15099762e-02, -3.30725584e-03, -1.06258357e-02,
254  3.24091535e-03, 9.82976528e-03, -3.17252619e-03, -9.10914602e-03,
255  3.10218270e-03, 8.45375526e-03, -3.02986363e-03, -7.85506037e-03,
256  2.95566681e-03, 7.30590531e-03, -2.87992593e-03, -6.80042517e-03,
257  2.80309157e-03, 6.33392168e-03, -2.72549458e-03, -5.90256016e-03,
258  2.64722721e-03, 5.50299075e-03, -2.56826006e-03, -5.13211877e-03,
259  2.48867691e-03, 4.78713456e-03, -2.40879175e-03, -4.46567920e-03,
260  2.32903023e-03, 4.16590632e-03, -2.24969322e-03, -3.88631909e-03,
261  2.17083836e-03, 3.62549780e-03, -2.09239740e-03, -3.38195207e-03,
262  2.01441139e-03, 3.15421411e-03, -1.93714812e-03, -2.94105418e-03,
263  1.86098398e-03, 2.74158159e-03, -1.78616796e-03, -2.55511292e-03,
264  1.71270364e-03, 2.38092463e-03, -1.64046686e-03, -2.21812548e-03,
265  1.56944142e-03, 2.06576607e-03, -1.49983691e-03, -1.92306771e-03,
266  1.43197091e-03, 1.78953434e-03, -1.36603342e-03, -1.66482974e-03,
267  1.30196917e-03, 1.54853763e-03, -1.23959572e-03, -1.44004013e-03,
268  1.17883935e-03, 1.33863258e-03, -1.11985330e-03, -1.24375656e-03,
269  1.06290017e-03, 1.15511541e-03, -1.00811665e-03, -1.07255448e-03,
270  9.55396053e-04, 9.95823582e-04, -9.04506526e-04, -9.24457717e-04,
271  8.55327235e-04, 8.57893516e-04, -8.07966638e-04, -7.95703775e-04,
272  7.62645097e-04, 7.37714244e-04, -7.19459630e-04, -6.83884826e-04,
273  6.78266519e-04, 6.34073000e-04, -6.38799615e-04, -5.87915177e-04,
274  6.00906492e-04, 5.44943812e-04, -5.64666746e-04, -5.04822423e-04,
275  5.30274583e-04, 4.67462815e-04, -4.97803536e-04, -4.32906634e-04,
276  4.67089022e-04, 4.01089094e-04, -4.37846574e-04, -3.71720589e-04,
277  4.09907918e-04, 3.44404015e-04, -3.83339184e-04, -3.18869979e-04,
278  3.58323383e-04, 2.95094182e-04, -3.34925025e-04, -2.73179107e-04,
279  3.12972551e-04, 2.53117874e-04, -2.92176448e-04, -2.34675958e-04,
280  2.72365201e-04, 2.17508613e-04, -2.53603369e-04, -2.01396171e-04,
281  2.36073945e-04, 1.86361512e-04, -2.19842838e-04, -1.72551841e-04,
282  2.04741179e-04, 1.60002621e-04, -1.90483310e-04, -1.48519365e-04,
283  1.76902612e-04, 1.37795143e-04, -1.64069468e-04, -1.27645963e-04,
284  1.52173503e-04, 1.18128308e-04, -1.41287958e-04, -1.09420999e-04,
285  1.31251892e-04, 1.01589203e-04, -1.21788072e-04, -9.44663026e-05,
286  1.12738710e-04, 8.77714836e-05, -1.04183332e-04, -8.13451855e-05,
287  9.63209473e-05, 7.52667260e-05, -8.92343436e-05, -6.97362337e-05,
288  8.27722214e-05, 6.48387355e-05, -7.66670221e-05, -6.04261010e-05,
289  7.07706056e-05, 5.62346999e-05, -6.51720769e-05, -5.21209285e-05,
290  6.00799054e-05, 4.81789030e-05, -5.55861831e-05, -4.46224615e-05,
291  5.15487245e-05, 4.15493180e-05, -4.77088656e-05, -3.88230703e-05,
292  4.39271165e-05, 3.61909193e-05, -4.03009685e-05, -3.35192549e-05,
293  3.70469977e-05, 3.09114062e-05, -3.42651102e-05, -2.85896985e-05,
294  3.18206326e-05, 2.66596591e-05, -2.94621038e-05, -2.49920735e-05,
295  2.70569244e-05, 2.33407512e-05, -2.47091627e-05, -2.15781546e-05,
296  2.26416598e-05, 1.98131916e-05, -2.09602777e-05, -1.82733110e-05,
297  1.95359940e-05, 1.70687447e-05, -1.81226984e-05, -1.60745977e-05,
298  1.65928502e-05, 1.50486440e-05, -1.50552921e-05, -1.38669855e-05,
299  1.37373646e-05, 1.26418718e-05, -1.27491643e-05, -1.16038210e-05,
300  1.19656476e-05, 1.08658854e-05, -1.11444373e-05, -1.03057618e-05,
301  1.01614907e-05, 9.68360712e-06, -9.12892437e-06, -8.87771523e-06,
302  8.27714049e-06, 8.00235344e-06, -7.71909503e-06, -7.28989912e-06,
303  7.33241173e-06, 6.85511870e-06, -6.87719924e-06, -6.57729069e-06,
304  6.23172852e-06, 6.21803201e-06, -5.51026882e-06, -5.65698534e-06,
305  4.94522316e-06, 5.00966549e-06, -4.65140639e-06, -4.50960504e-06,
306  4.50816730e-06, 4.27264119e-06, -4.27721386e-06, -4.17904541e-06,
307  3.83829877e-06, 3.99135588e-06, -3.30706353e-06, -3.59007062e-06,
308  2.91718177e-06, 3.09149959e-06, -2.78464422e-06, -2.72991556e-06,
309  2.78988692e-06, 2.62184526e-06, -2.69562226e-06, -2.64820266e-06,
310  2.38253026e-06, 2.57212448e-06, -1.96710869e-06, -2.27466668e-06,
311  1.68382189e-06, 1.87265988e-06, -1.64939057e-06, -1.60086337e-06,
312  1.74492415e-06, 1.57625873e-06, -1.73375546e-06, -1.68018495e-06,
313  1.49713605e-06, 1.67617565e-06, -1.15209033e-06, -1.44565738e-06,
314  9.33568272e-07, 1.10580727e-06, -9.58738388e-07, -8.91709022e-07,
315  1.10912250e-06, 9.20648047e-07, -1.14843279e-06, -1.07424749e-06,
316  9.58269576e-07, 1.11630634e-06, -6.55978270e-07, -9.28498858e-07,
317  4.76804701e-07, 6.28233782e-07, -5.38190060e-07, -4.50811695e-07,
318  7.21907128e-07, 5.13721134e-07, -7.91898561e-07, -6.98775094e-07,
319  6.29976309e-07, 7.69949518e-07, -3.53680414e-07, -6.09083564e-07,
320  1.98436165e-07, 3.33739875e-07, -2.81850268e-07, -1.79363102e-07,
321  4.85847782e-07, 2.63576622e-07, -5.74511016e-07, -4.68319242e-07,
322  4.29779921e-07, 5.57683825e-07, -1.69312228e-07, -4.13618087e-07,
323  2.86418681e-08, 1.53785781e-07, -1.25475908e-07, -1.37261001e-08,
324  3.41831748e-07, 1.11150848e-07, -4.41876227e-07, -3.28081137e-07,
325  3.07626614e-07, 4.28685999e-07, -5.68117301e-08, -2.94983550e-07,
326  -7.49686950e-08, 4.47029437e-08, -3.00533793e-08, 8.65554823e-08,
327  2.53951480e-07, 1.89772734e-08, -3.60943678e-07, -2.43375449e-07,
328  2.33093935e-07, 3.50856976e-07, 1.18260092e-08, -2.23484904e-07,
329  -1.38176528e-07, -2.09699988e-08, 2.81524357e-08, 1.46868453e-07,
330  2.00354017e-07, -3.64049742e-08, -3.11591993e-07, -1.92528478e-07,
331  1.87653644e-07, 3.04180734e-07, 5.36633725e-08, -1.80643006e-07,
332  -1.76695163e-07, -6.02879224e-08, 6.36137567e-08, 1.82948329e-07,
333  1.67709606e-07, -6.95097238e-08, -2.81543035e-07, -1.62157223e-07,
334  1.59995687e-07, 2.76320615e-07, 7.91191565e-08, -1.55089035e-07,
335  -2.00122815e-07, -8.37247404e-08, 8.51732622e-08, 2.04441857e-07,
336  1.47871244e-07, -8.92194564e-08, -2.63290375e-07, -1.44085068e-07,
337  1.43203505e-07, 2.59751664e-07, 9.45666953e-08, -1.39899410e-07,
338  -2.14332489e-07, -9.76492812e-08, 9.82430663e-08, 2.17206296e-07,
339  1.35851448e-07, -1.00919826e-07, -2.52237817e-07, -1.33361007e-07,
340  1.33041511e-07, 2.49923380e-07, 1.03909218e-07, -1.30892598e-07};
341 
342 
343 extern const double mey512[512]=
344 { 7.10493973e-01, 4.50105050e-01, -3.38583052e-03, -1.49893593e-01,
345  3.38163200e-03, 8.97669622e-02, -3.37437293e-03, -6.39386112e-02,
346  3.36395551e-03, 4.95432853e-02, -3.35052556e-03, -4.03450647e-02,
347  3.33435448e-03, 3.39457751e-02, -3.31560344e-03, -2.92265276e-02,
348  3.29420529e-03, 2.55950296e-02, -3.26998220e-03, -2.27082463e-02,
349  3.24288110e-03, 2.03537102e-02, -3.21309123e-03, -1.83929344e-02,
350  3.18092593e-03, 1.67319834e-02, -3.14658664e-03, -1.53049395e-02,
351  3.11004467e-03, 1.40640442e-02, -3.07115860e-03, -1.29736823e-02,
352  3.02990966e-03, 1.20067463e-02, -2.98651908e-03, -1.11424056e-02,
353  2.94132983e-03, 1.03645576e-02, -2.89457042e-03, -9.66055605e-03,
354  2.84623667e-03, 9.02015914e-03, -2.79620900e-03, -8.43482325e-03,
355  2.74448774e-03, 7.89738929e-03, -2.69131044e-03, -7.40199712e-03,
356  2.63703363e-03, 6.94396076e-03, -2.58189657e-03, -6.51946507e-03,
357  2.52590301e-03, 6.12518664e-03, -2.46893857e-03, -5.75806376e-03,
358  2.41100600e-03, 5.41532586e-03, -2.35234263e-03, -5.09465878e-03,
359  2.29330211e-03, 4.79426613e-03, -2.23411832e-03, -4.51270574e-03,
360  2.17478716e-03, 4.24861679e-03, -2.11518405e-03, -4.00057150e-03,
361  2.05529934e-03, 3.76716783e-03, -1.99535580e-03, -3.54724443e-03,
362  1.93569056e-03, 3.33998117e-03, -1.87651914e-03, -3.14476677e-03,
363  1.81781735e-03, 2.96095098e-03, -1.75943900e-03, -2.78771644e-03,
364  1.70135136e-03, 2.62418785e-03, -1.64375292e-03, -2.46966028e-03,
365  1.58695542e-03, 2.32371079e-03, -1.53114803e-03, -2.18607515e-03,
366  1.47627945e-03, 2.05640760e-03, -1.42217568e-03, -1.93415902e-03,
367  1.36877577e-03, 1.81869136e-03, -1.31624961e-03, -1.70951048e-03,
368  1.26488015e-03, 1.60638137e-03, -1.21482769e-03, -1.50920820e-03,
369  1.16601218e-03, 1.41779664e-03, -1.11823104e-03, -1.33173436e-03,
370  1.07139507e-03, 1.25050747e-03, -1.02564633e-03, -1.17373488e-03,
371  9.81240476e-04, 1.10128507e-03, -9.38311185e-04, -1.03315730e-03,
372  8.96752444e-04, 9.69244962e-04, -8.56336572e-04, -9.09217036e-04,
373  8.16950091e-04, 8.52635196e-04, -7.78711755e-04, -7.99188920e-04,
374  7.41854891e-04, 7.48812799e-04, -7.06491882e-04, -7.01568186e-04,
375  6.72496507e-04, 6.57407050e-04, -6.39621979e-04, -6.16053707e-04,
376  6.07736844e-04, 5.77122285e-04, -5.76943027e-04, -5.40352082e-04,
377  5.47458177e-04, 5.05725097e-04, -5.19380158e-04, -4.73347879e-04,
378  4.92569378e-04, 4.43215546e-04, -4.66766820e-04, -4.15093655e-04,
379  4.41829924e-04, 3.88635799e-04, -4.17850610e-04, -3.63619059e-04,
380  3.95037601e-04, 3.40061635e-04, -3.73480883e-04, -3.18104773e-04,
381  3.53034002e-04, 2.97776851e-04, -3.33432052e-04, -2.78875312e-04,
382  3.14527524e-04, 2.61084311e-04, -2.96408287e-04, -2.44210221e-04,
383  2.79279866e-04, 2.28299299e-04, -2.63229856e-04, -2.13519662e-04,
384  2.48110176e-04, 1.99925402e-04, -2.33655007e-04, -1.87338559e-04,
385  2.19716596e-04, 1.75466793e-04, -2.06383186e-04, -1.64138929e-04,
386  1.93861253e-04, 1.53422646e-04, -1.82239870e-04, -1.43506489e-04,
387  1.71372913e-04, 1.34464008e-04, -1.60996961e-04, -1.26135762e-04,
388  1.50967048e-04, 1.18247020e-04, -1.41374558e-04, -1.10643329e-04,
389  1.32429417e-04, 1.03408242e-04, -1.24224420e-04, -9.67453457e-05,
390  1.16617400e-04, 9.07424406e-05, -1.09349093e-04, -8.52535646e-05,
391  1.02278849e-04, 8.00167265e-05, -9.55025085e-05, -7.48895008e-05,
392  8.92345538e-05, 6.99667838e-05, -8.35724159e-05, -6.54628522e-05,
393  7.83786137e-05, 6.14755686e-05, -7.33985944e-05, -5.78684340e-05,
394  6.84964282e-05, 5.43883476e-05, -6.37726560e-05, -5.09012284e-05,
395  5.94464334e-05, 4.75097984e-05, -5.56198142e-05, -4.44356672e-05,
396  5.21598775e-05, 4.17835628e-05, -4.88165543e-05, -3.94234083e-05,
397  4.54583110e-05, 3.71081051e-05, -4.21899877e-05, -3.47091780e-05,
398  3.92349341e-05, 3.23345806e-05, -3.66992861e-05, -3.02108027e-05,
399  3.44540868e-05, 2.84471212e-05, -3.22531073e-05, -2.69176971e-05,
400  2.99685273e-05, 2.53793769e-05, -2.77087707e-05, -2.37073561e-05,
401  2.57006403e-05, 2.20130032e-05, -2.40535937e-05, -2.05259829e-05,
402  2.26418630e-05, 1.93585240e-05, -2.12222765e-05, -1.83875304e-05,
403  1.96699394e-05, 1.73723962e-05, -1.80960721e-05, -1.61906832e-05,
404  1.67301470e-05, 1.49559571e-05, -1.56841670e-05, -1.38999239e-05,
405  1.48347873e-05, 1.31367082e-05, -1.39411401e-05, -1.25449747e-05,
406  1.28805182e-05, 1.18857527e-05, -1.17662176e-05, -1.10381227e-05,
407  1.08296779e-05, 1.01170613e-05, -1.01847652e-05, -9.35558678e-06,
408  9.70989693e-06, 8.86904380e-06, -9.16586978e-06, -8.53723207e-06,
409  8.43154781e-06, 8.12223683e-06, -7.62170889e-06, -7.50412092e-06,
410  6.96918995e-06, 6.79877622e-06, -6.58917319e-06, -6.24007395e-06,
411  6.36131473e-06, 5.94415493e-06, -6.04757568e-06, -5.79156150e-06,
412  5.52791386e-06, 5.54507187e-06, -4.91813425e-06, -5.08539585e-06,
413  4.45203824e-06, 4.52903010e-06, -4.24571307e-06, -4.11041193e-06,
414  4.17966331e-06, 3.94621241e-06, -4.01664463e-06, -3.91747061e-06,
415  3.63735844e-06, 3.78742988e-06, -3.15830574e-06, -3.43723581e-06,
416  2.81393946e-06, 2.98379398e-06, -2.72095664e-06, -2.66192601e-06,
417  2.76043189e-06, 2.58866444e-06, -2.69565313e-06, -2.64538810e-06,
418  2.40781851e-06, 2.59565951e-06, -2.01389286e-06, -2.32092414e-06,
419  1.74876250e-06, 1.93836981e-06, -1.72952932e-06, -1.68308433e-06,
420  1.83764570e-06, 1.67235142e-06, -1.83675144e-06, -1.78778666e-06,
421  1.60837239e-06, 1.79317530e-06, -1.26977880e-06, -1.57017243e-06,
422  1.05614207e-06, 1.23616353e-06, -1.08483019e-06, -1.02642323e-06,
423  1.23754344e-06, 1.05841191e-06, -1.27815217e-06, -1.21391182e-06,
424  1.08839649e-06, 1.25686512e-06, -7.85746162e-07, -1.06907453e-06,
425  6.05558756e-07, 7.68064793e-07, -6.65376007e-07, -5.89242374e-07,
426  8.47059978e-07, 6.50192461e-07, -9.14631182e-07, -8.32815083e-07,
427  7.49969032e-07, 9.01163121e-07, -4.70672877e-07, -7.37143320e-07,
428  3.12221293e-07, 4.58378507e-07, -3.92269027e-07, -3.00368021e-07,
429  5.92783323e-07, 3.80785029e-07, -6.77882136e-07, -5.81612514e-07,
430  5.29535101e-07, 6.66981139e-07, -2.65425466e-07, -5.18870518e-07,
431  1.21110269e-07, 2.54972098e-07, -2.14317680e-07, -1.10850236e-07,
432  4.27083269e-07, 2.04239738e-07, -5.23588141e-07, -4.17181631e-07,
433  3.85860226e-07, 5.13860803e-07, -1.31636960e-07, -3.76307458e-07,
434  -3.47381479e-09, 1.22260587e-07, -9.83037647e-08, 1.26703641e-08,
435  3.19049659e-07, 8.92921966e-08, -4.22985860e-07, -3.10229518e-07,
436  2.92177900e-07, 4.14363873e-07, -4.43981636e-08, -2.83760147e-07,
437  -8.47125812e-08, 3.61898674e-08, -2.26525947e-08, 9.27066003e-08,
438  2.48602585e-07, 1.48777684e-08, -3.57385893e-07, -2.41051851e-07,
439  2.31092271e-07, 3.50063432e-07, 1.24833882e-08, -2.24000884e-07,
440  -1.37678980e-07, -1.93422935e-08, 2.66672827e-08, 1.44304767e-07,
441  2.02679507e-07, -3.30594651e-08, -3.14627088e-07, -1.96521276e-07,
442  1.91280784e-07, 3.08702428e-07, 4.95500446e-08, -1.85588018e-07,
443  -1.72189419e-07, -5.50138378e-08, 5.87967660e-08, 1.77427745e-07,
444  1.72768179e-07, -6.38130678e-08, -2.86782308e-07, -1.67970564e-07,
445  1.65360928e-07, 2.82199598e-07, 7.36773437e-08, -1.60988361e-07,
446  -1.94647503e-07, -7.78454510e-08, 7.97001447e-08, 1.98617113e-07,
447  1.53313254e-07, -8.34768690e-08, -2.68676869e-07, -1.49724187e-07,
448  1.48512364e-07, 2.65270067e-07, 8.93557497e-08, -1.45281733e-07,
449  -2.09236571e-07, -9.24169614e-08, 9.32747632e-08, 2.12135158e-07,
450  1.40683752e-07, -9.60169615e-08, -2.56927920e-07, -1.38092286e-07,
451  1.37583411e-07, 2.54481563e-07, 9.95215647e-08, -1.35276034e-07,
452  -2.18692260e-07, -1.01696569e-07, 1.02069320e-07, 2.20741387e-07,
453  1.32505104e-07, -1.03998389e-07, -2.49323071e-07, -1.30690960e-07,
454  1.30512721e-07, 2.47618852e-07, 1.06095460e-07, -1.28913016e-07,
455  -2.24804194e-07, -1.07596454e-07, 1.07751334e-07, 2.26212111e-07,
456  1.27223603e-07, -1.09071092e-07, -2.44414725e-07, -1.25987801e-07,
457  1.25951644e-07, 2.43258848e-07, 1.10333813e-07, -1.24871296e-07,
458  -2.28742719e-07, -1.11343388e-07, 1.11410960e-07, 2.29686092e-07,
459  1.23823810e-07, -1.12291985e-07, -2.41257033e-07, -1.23001982e-07,
460  1.23019136e-07, 2.40491414e-07, 1.13057321e-07, -1.22306368e-07,
461  -2.31272322e-07, -1.13720961e-07, 1.13760301e-07, 2.31890395e-07,
462  1.21642470e-07, -1.14335688e-07, -2.39232315e-07, -1.21107556e-07,
463  1.21140049e-07, 2.38735804e-07, 1.14801440e-07, -1.20679489e-07,
464  -2.32891432e-07, -1.15228879e-07, 1.15263272e-07, 2.33288441e-07,
465  1.20247808e-07, -1.15631902e-07, -2.37938703e-07, -1.19906138e-07,
466  1.19940329e-07, 2.37622682e-07, 1.15914360e-07, -1.19648238e-07,
467  -2.33924171e-07, -1.16184642e-07, 1.16221580e-07, 2.34174666e-07,
468  1.19359023e-07, -1.16453708e-07, -2.37114915e-07, -1.19144431e-07,
469  1.19176927e-07, 2.36917106e-07, 1.16622092e-07, -1.18994732e-07,
470  -2.34580653e-07, -1.16790263e-07, 1.16830525e-07, 2.34736310e-07,
471  1.18794590e-07, -1.16974611e-07, -2.36592207e-07, -1.18661688e-07};
472 
473 ClassImp(Meyer<double>)
474 
475 // constructors
476 
477 template<class DataType_t> Meyer<DataType_t>::
478 Meyer() :
479 WaveDWT<DataType_t>()
480 {
481  pLForward = NULL;
482  pLInverse = NULL;
483  pHForward = NULL;
484  pHInverse = NULL;
485  setFilter();
486 }
487 
488 template<class DataType_t> Meyer<DataType_t>::
489 Meyer(const Wavelet &w) :
490 WaveDWT<DataType_t>(w)
491 {
492  pLForward = NULL;
493  pLInverse = NULL;
494  pHForward = NULL;
495  pHInverse = NULL;
496  setFilter();
497 }
498 
499 template<class DataType_t> Meyer<DataType_t>::
500 Meyer(const Meyer<DataType_t> &w) :
501 WaveDWT<DataType_t>(w)
502 {
503  pLForward = NULL;
504  pLInverse = NULL;
505  pHForward = NULL;
506  pHInverse = NULL;
507  setFilter();
508 }
509 
510 template<class DataType_t> Meyer<DataType_t>::
511 Meyer(int m, int tree, enum BORDER border) :
512 WaveDWT<DataType_t>(m,m,tree,border)
513 {
514  pLForward = NULL;
515  pLInverse = NULL;
516  pHForward = NULL;
517  pHInverse = NULL;
518  setFilter();
519 }
520 
521 // destructor
522 template<class DataType_t>
523 Meyer<DataType_t>::~Meyer()
524 {
525  if(pLForward) free(pLForward);
526  if(pLInverse) free(pLInverse);
527  if(pHForward) free(pHForward);
528  if(pHInverse) free(pHInverse);
529 }
530 
531 // clone
532 template<class DataType_t>
533 Meyer<DataType_t>* Meyer<DataType_t>::Clone() const
534 {
535  return new Meyer<DataType_t>(*this);
536 }
537 
538 template<class DataType_t>
539 int Meyer<DataType_t>::getMaxLevel()
540 {
541  int i;
542  int maxLevel = 0;
543  int n = int(this->nWWS);
544  int m = this->m_L;
545  int N[10] = {512,384,256,192,128,96,64,48,32,16};
546 
547  if(this->m_TreeType == 2) { // adaptive wavelet
548  for(i=0; i<10; i++) {if(N[i] == this->m_L/2) break;} // find initial wavelet
549 
550  while((n>=2*m) && !(n&1)){
551  n /= 2;
552  if(i < 10) m = N[i++]*2;
553  else m = 32;
554  maxLevel++;
555  }
556  return maxLevel;
557  }
558 
559  else return Wavelet::getMaxLevel(n);
560 }
561 
562 //: calculate wavelet filter
563 //!param: taper function order n
564 //!param: beta(n,n) - value of Euler's beta function
565 //!param: integration step
566 template<class DataType_t>
567 double Meyer<DataType_t>::filter(int n, double B, double s)
568 {
569  int i;
570  int k = this->m_H;
571  double pi3 = PI/3;
572  double pi2 = PI/2;
573  double x = s;
574  double y = 0;
575  double z = 0;
576  double v = 0;
577 
578  if(n<2) n = 2;
579  B = pow(B,1./(n-1));
580 
581  double* H = (double*)malloc(k*sizeof(double));
582  for(i=0; i<k; i++) H[i] = 0.;
583 
584  while(x<1.){
585  v += s*pow(x*(1-x)/B,n-1);
586  y = cos(v*pi2);
587  z = pi3*(1+x);
588  for(i=0; i<k; i++) H[i] += y*cos(i*z)*s;
589  x += s;
590 // if(int(x/0.1)*0.1==x) printf(" %6.5f \n",x);
591  }
592 
593  v = y = 0.;
594  printf("\n");
595  for(i=0; i<k; i++){
596  H[i] += i ? sin(pi3*i)/pi3/i : 1;
597  H[i] *= sqrt(2.)/3.;
598  v += 2*H[i]*H[i];
599  y += 2*H[i];
600  printf(" %16.8e,",H[i]);
601  if(((i+1)/4)*4==i+1) printf("\n");
602  }
603  printf("\n normalization: %9.6f, orthogonality: %9.6f \n",y-H[0],v-H[0]*H[0]);
604  free(H);
605  return v-H[0]*H[0];
606 }
607 
608 
609 template<class DataType_t>
610 void Meyer<DataType_t>::setFilter()
611 {
612  const double* pF;
613  this->m_H = (this->m_H>>1)<<1;
614  int n = this->m_H/2;
615  switch(n)
616  {
617  case 512: pF = mey512; break;
618  case 384: pF = mey384; break;
619  case 256: pF = mey256; break;
620  case 192: pF = mey128; break;
621  case 128: pF = mey064; break;
622  case 96: pF = mey032; break;
623  case 64: pF = mey016; break;
624  case 48: pF = mey008; break;
625  case 32: pF = mey004; break;
626  case 16: pF = mey; break;
627  default: pF = mey;
628  this->m_H = this->m_L = 32;
629  break;
630  }
631 
632  if(!pLForward) {
633  pLInverse = (double*)malloc(this->m_H*sizeof(double));
634  pLForward = (double*)malloc(this->m_H*sizeof(double));
635  pHInverse = (double*)malloc(this->m_H*sizeof(double));
636  pHForward = (double*)malloc(this->m_H*sizeof(double));
637  }
638  else {
639  pLInverse = (double*)realloc(pLInverse,this->m_H*sizeof(double));
640  pLForward = (double*)realloc(pLForward,this->m_H*sizeof(double));
641  pHInverse = (double*)realloc(pHInverse,this->m_H*sizeof(double));
642  pHForward = (double*)realloc(pHForward,this->m_H*sizeof(double));
643  }
644 
645  double* temp = (double*)malloc(this->m_H*sizeof(double));
646 
647  temp[0] = 0.;
648  for(int i=0; i<this->m_H/2; i++){
649  temp[i+this->m_H/2] = pF[i];
650  temp[this->m_H/2-i] = pF[i];
651  }
652 
653 // LP filter for db3: h0 h1 h2 h3 h4 h5
654 // HP filter for db3: h5 -h4 h3 -h2 h1 -h0
655 // iLP filter for db3: h4 h1 h2 h3 h0 h5
656 // iHP filter for db3: h5 -h0 h3 -h2 h1 -h4
657 
658  for(int i=0; i<this->m_H; i+=2){
659 
660  pLForward[i] = temp[i];
661  pLForward[i+1] = temp[i+1];
662  pHForward[i] = temp[this->m_H-1-i];
663  pHForward[i+1] = -temp[this->m_H-2-i];
664 
665  pLInverse[i] = temp[this->m_H-2-i];
666  pLInverse[i+1] = temp[i+1];
667  pHInverse[i] = temp[this->m_H-1-i];
668  pHInverse[i+1] = -temp[i];
669 
670  }
671 
672  this->m_WaveType = MEYER;
673  free(temp);
674 }
675 
676 // forward function does one step of forward transformation.
677 // <level> input parameter is the level to be reconstructed
678 // <layer> input parameter is the layer to be reconstructed.
679 template<class DataType_t>
680 void Meyer<DataType_t>::forward(int level,int layer)
681 {
682  if(this->m_TreeType == 2 && !layer){ // adaptive wavelet
683  int i;
684  int N[10] = {512,384,256,192,128,96,64,48,32,16};
685  for(i=0; i<10; i++) {if(N[i] == this->m_L/2) break;} // find initial wavelet
686  if(i+level > 9) { this->m_H = 32; }
687  else { this->m_H = N[i+level]*2; }
688  setFilter();
689  }
690  this->forwardFWT(level, layer, pLForward, pHForward);
691 }
692 
693 // inverse function does one step of inverse transformation.
694 // <level> input parameter is the level to be reconstructed
695 // <layer> input parameter is the layer to be reconstructed.
696 template<class DataType_t>
697 void Meyer<DataType_t>::inverse(int level,int layer)
698 {
699  if(this->m_TreeType == 2 && !layer) { // adaptive wavelet
700  int i;
701  int N[10] = {512,384,256,192,128,96,64,48,32,16};
702  for(i=0; i<10; i++) {if(N[i] == this->m_L/2) break;} // find initial wavelet
703  if(i+level > 9) { this->m_H = 32; }
704  else { this->m_H = N[i+level]*2; }
705  setFilter();
706  }
707  this->inverseFWT(level, layer, pLInverse, pHInverse);
708 }
709 
710 // instantiations
711 
712 #define CLASS_INSTANTIATION(class_) template class Meyer< class_ >;
713 
714 CLASS_INSTANTIATION(float)
715 CLASS_INSTANTIATION(double)
716 //CLASS_INSTANTIATION(std::complex<float>)
717 //CLASS_INSTANTIATION(std::complex<double>)
718 
719 #undef CLASS_INSTANTIATION
720 
721 //} // end namespace wat
722 //} // end namespace datacondAPI
723 
tree
TTree * tree
Definition: TimeSortTree.C:20
Meyer::pLInverse
double * pLInverse
Definition: Meyer.hh:25
Meyer::setFilter
void setFilter()
Definition: Meyer.cc:610
printf
printf("total live time: non-zero lags = %10.1f \n", liveTot)
mey512
const double mey512[512]
v
WSeries< float > v[nIFO]
Definition: cwb_net.C:62
B
#define B
n
int n
Definition: cwb_net.C:10
z
wavearray< double > z
Definition: Test10.C:32
mey384
const double mey384[384]
BORDER
BORDER
Definition: Wavelet.hh:18
mey064
const double mey064[128]
mey016
const double mey016[64]
x
cout<< endl;cout<< "ts size = "<< ts.size()<< " ts rate = "<< ts.rate()<< endl;tf.Forward(ts, wdm);int levels=tf.getLevel();cout<< "tf size = "<< tf.size()<< endl;double dF=tf.resolution();double dT=1./(2 *dF);cout<< "rate(hz) : "<< RATE<< "\t layers : "<< nLAYERS<< "\t dF(hz) : "<< dF<< "\t dT(ms) : "<< dT *1000.<< endl;int itime=TIME_PIXEL_INDEX;int ifreq=FREQ_PIXEL_INDEX;int index=(levels+1)*itime+ifreq;double time=itime *dT;double freq=(ifreq >0)?ifreq *dF:dF/4;cout<< endl;cout<< "PIXEL TIME = "<< time<< " sec "<< endl;cout<< "PIXEL FREQ = "<< freq<< " Hz "<< endl;cout<< endl;wavearray< double > x
Definition: BaseFunctionWDM.C:51
mey004
const double mey004[32]
Meyer::Clone
virtual Meyer * Clone() const
return: Wavelet* - duplicate of *this, allocated on heap
Definition: Meyer.cc:533
m
int m
Definition: cwb_net.C:10
i
i drho i
Definition: cwb_epparameters.C:85
Meyer::filter
double filter(int, double, double=1.e-6)
param: taper function order n param: beta(n,n) - value of Euler's beta function param: integration st...
Definition: Meyer.cc:567
N
#define N
PI
#define PI
Definition: watfun.hh:14
w
wavearray< double > w
Definition: Test1.C:27
mey032
const double mey032[96]
Meyer::inverse
virtual void inverse(int level, int layer)
Definition: Meyer.cc:697
mey256
const double mey256[256]
Meyer::~Meyer
virtual ~Meyer()
Definition: Meyer.cc:523
Meyer::Meyer
Meyer()
Definition: Meyer.cc:478
setFilter
r setFilter(10)
int
i() int(T_cor *100))
Meyer::forward
virtual void forward(int level, int layer)
Definition: Meyer.cc:680
MEYER
Definition: Wavelet.hh:30
k
int k
Definition: cwb_report_prod_2.C:150
s
s s
Definition: cwb_net.C:137
Meyer::pHInverse
double * pHInverse
Definition: Meyer.hh:29
Meyer
Definition: Meyer.hh:18
Meyer::pHForward
double * pHForward
Definition: Meyer.hh:27
CLASS_INSTANTIATION
#define CLASS_INSTANTIATION(class_)
Definition: Meyer.cc:712
mey008
const double mey008[48]
Meyer::getMaxLevel
int getMaxLevel()
Definition: Meyer.cc:539
Wavelet::getMaxLevel
virtual int getMaxLevel(int)
Definition: Wavelet.hh:116
y
wavearray< double > y
Definition: Test10.C:31
Meyer::pLForward
double * pLForward
Definition: Meyer.hh:23
mey128
const double mey128[192]
mey
const double mey[16]