libs/qm-dsp/dsp/wavelet/Wavelet.cpp

   1 /* -*- c-basic-offset: 4 indent-tabs-mode: nil -*-  vi:set ts=8 sts=4 sw=4: */
   2
   3 /*
   4     QM DSP Library
   5
   6     Centre for Digital Music, Queen Mary, University of London.
   7     This file copyright 2009 Thomas Wilmering.
   8
   9     This program is free software; you can redistribute it and/or
  10     modify it under the terms of the GNU General Public License as
  11     published by the Free Software Foundation; either version 2 of the
  12     License, or (at your option) any later version.  See the file
  13     COPYING included with this distribution for more information.
  14 */
  15
  16 #ifdef COMPILER_MSVC
  17 #pragma warning(disable:4305)
  18 #endif
  19
  20 #include "Wavelet.h"
  21
  22 #include <cassert>
  23
  24 std::string
  25 Wavelet::getWaveletName(Type wavelet)
  26 {
  27     switch (wavelet) {
  28         case Haar: return "Haar";
  29         case Daubechies_2: return "Daubechies 2";
  30         case Daubechies_3: return "Daubechies 3";
  31         case Daubechies_4: return "Daubechies 4";
  32         case Daubechies_5: return "Daubechies 5";
  33         case Daubechies_6: return "Daubechies 6";
  34         case Daubechies_7: return "Daubechies 7";
  35         case Daubechies_8: return "Daubechies 8";
  36         case Daubechies_9: return "Daubechies 9";
  37         case Daubechies_10: return "Daubechies 10";
  38         case Daubechies_20: return "Daubechies 20";
  39         case Daubechies_40: return "Daubechies 40";
  40         case Symlet_2: return "Symlet 2";
  41         case Symlet_3: return "Symlet 3";
  42         case Symlet_4: return "Symlet 4";
  43         case Symlet_5: return "Symlet 5";
  44         case Symlet_6: return "Symlet 6";
  45         case Symlet_7: return "Symlet 7";
  46         case Symlet_8: return "Symlet 8";
  47         case Symlet_9: return "Symlet 9";
  48         case Symlet_10: return "Symlet 10";
  49         case Symlet_20: return "Symlet 20";
  50         case Symlet_30: return "Symlet 30";
  51         case Coiflet_1: return "Coiflet 1";
  52         case Coiflet_2: return "Coiflet 2";
  53         case Coiflet_3: return "Coiflet 3";
  54         case Coiflet_4: return "Coiflet 4";
  55         case Coiflet_5: return "Coiflet 5";
  56         case Biorthogonal_1_3: return "Biorthogonal 1.3";
  57         case Biorthogonal_1_5: return "Biorthogonal 1.5";
  58         case Biorthogonal_2_2: return "Biorthogonal 2.2";
  59         case Biorthogonal_2_4: return "Biorthogonal 2.4";
  60         case Biorthogonal_2_6: return "Biorthogonal 2.6";
  61         case Biorthogonal_2_8: return "Biorthogonal 2.8";
  62         case Biorthogonal_3_1: return "Biorthogonal 3.1";
  63         case Biorthogonal_3_3: return "Biorthogonal 3.3";
  64         case Biorthogonal_3_5: return "Biorthogonal 3.5";
  65         case Biorthogonal_3_7: return "Biorthogonal 3.7";
  66         case Biorthogonal_3_9: return "Biorthogonal 3.9";
  67         case Biorthogonal_4_4: return "Biorthogonal 4.4";
  68         case Biorthogonal_5_5: return "Biorthogonal 5.5";
  69         case Biorthogonal_6_8: return "Biorthogonal 6.8";
  70         case Meyer: return "Meyer";
  71     }
  72
  73     return "(unknown)";
  74 }
  75
  76 void
  77 Wavelet::createDecompositionFilters(Type wavelet,
  78                                     std::vector<float> &lpd,
  79                                     std::vector<float> &hpd)
  80 {
  81     lpd.clear();
  82     hpd.clear();
  83
  84     unsigned int flength = 0;
  85
  86     switch (wavelet) {
  87
  88     case Haar:
  89         lpd.push_back(0.70710678118655);
  90         lpd.push_back(0.70710678118655);
  91         hpd.push_back(-0.70710678118655);
  92         hpd.push_back(0.70710678118655);
  93         flength = 2;
  94         break;
  95
  96     case Daubechies_2:
  97         lpd.push_back(-0.12940952255092);
  98         lpd.push_back(0.22414386804186);
  99         lpd.push_back(0.83651630373747);
 100         lpd.push_back(0.48296291314469);
 101         hpd.push_back(-0.48296291314469);
 102         hpd.push_back(0.83651630373747);
 103         hpd.push_back(-0.22414386804186);
 104         hpd.push_back(-0.12940952255092);
 105         flength = 4;
 106         break;
 107
 108     case Daubechies_3:
 109         lpd.push_back(0.03522629188210);
 110         lpd.push_back(-0.08544127388224);
 111         lpd.push_back(-0.13501102001039);
 112         lpd.push_back(0.45987750211933);
 113         lpd.push_back(0.80689150931334);
 114         lpd.push_back(0.33267055295096);
 115         hpd.push_back(-0.33267055295096);
 116         hpd.push_back(0.80689150931334);
 117         hpd.push_back(-0.45987750211933);
 118         hpd.push_back(-0.13501102001039);
 119         hpd.push_back(0.08544127388224);
 120         hpd.push_back(0.03522629188210);
 121         flength = 6;
 122         break;
 123
 124     case Daubechies_4:
 125         lpd.push_back(-0.01059740178500);
 126         lpd.push_back(0.03288301166698);
 127         lpd.push_back(0.03084138183599);
 128         lpd.push_back(-0.18703481171888);
 129         lpd.push_back(-0.02798376941698);
 130         lpd.push_back(0.63088076792959);
 131         lpd.push_back(0.71484657055254);
 132         lpd.push_back(0.23037781330886);
 133         hpd.push_back(-0.23037781330886);
 134         hpd.push_back(0.71484657055254);
 135         hpd.push_back(-0.63088076792959);
 136         hpd.push_back(-0.02798376941698);
 137         hpd.push_back(0.18703481171888);
 138         hpd.push_back(0.03084138183599);
 139         hpd.push_back(-0.03288301166698);
 140         hpd.push_back(-0.01059740178500);
 141         flength = 8;
 142         break;
 143
 144     case Daubechies_5:
 145         lpd.push_back(0.00333572528500);
 146         lpd.push_back(-0.01258075199902);
 147         lpd.push_back(-0.00624149021301);
 148         lpd.push_back(0.07757149384007);
 149         lpd.push_back(-0.03224486958503);
 150         lpd.push_back(-0.24229488706619);
 151         lpd.push_back(0.13842814590110);
 152         lpd.push_back(0.72430852843857);
 153         lpd.push_back(0.60382926979747);
 154         lpd.push_back(0.16010239797413);
 155         hpd.push_back(-0.16010239797413);
 156         hpd.push_back(0.60382926979747);
 157         hpd.push_back(-0.72430852843857);
 158         hpd.push_back(0.13842814590110);
 159         hpd.push_back(0.24229488706619);
 160         hpd.push_back(-0.03224486958503);
 161         hpd.push_back(-0.07757149384007);
 162         hpd.push_back(-0.00624149021301);
 163         hpd.push_back(0.01258075199902);
 164         hpd.push_back(0.00333572528500);
 165         flength = 10;
 166         break;
 167
 168     case Daubechies_6:
 169         lpd.push_back(-0.00107730108500);
 170         lpd.push_back(0.00477725751101);
 171         lpd.push_back(0.00055384220099);
 172         lpd.push_back(-0.03158203931803);
 173         lpd.push_back(0.02752286553002);
 174         lpd.push_back(0.09750160558708);
 175         lpd.push_back(-0.12976686756710);
 176         lpd.push_back(-0.22626469396517);
 177         lpd.push_back(0.31525035170924);
 178         lpd.push_back(0.75113390802158);
 179         lpd.push_back(0.49462389039839);
 180         lpd.push_back(0.11154074335008);
 181         hpd.push_back(-0.11154074335008);
 182         hpd.push_back(0.49462389039839);
 183         hpd.push_back(-0.75113390802158);
 184         hpd.push_back(0.31525035170924);
 185         hpd.push_back(0.22626469396517);
 186         hpd.push_back(-0.12976686756710);
 187         hpd.push_back(-0.09750160558708);
 188         hpd.push_back(0.02752286553002);
 189         hpd.push_back(0.03158203931803);
 190         hpd.push_back(0.00055384220099);
 191         hpd.push_back(-0.00477725751101);
 192         hpd.push_back(-0.00107730108500);
 193         flength = 12;
 194         break;
 195
 196     case Daubechies_7:
 197         lpd.push_back(0.00035371380000);
 198         lpd.push_back(-0.00180164070400);
 199         lpd.push_back(0.00042957797300);
 200         lpd.push_back(0.01255099855601);
 201         lpd.push_back(-0.01657454163102);
 202         lpd.push_back(-0.03802993693503);
 203         lpd.push_back(0.08061260915107);
 204         lpd.push_back(0.07130921926705);
 205         lpd.push_back(-0.22403618499417);
 206         lpd.push_back(-0.14390600392911);
 207         lpd.push_back(0.46978228740536);
 208         lpd.push_back(0.72913209084656);
 209         lpd.push_back(0.39653931948231);
 210         lpd.push_back(0.07785205408506);
 211         hpd.push_back(-0.07785205408506);
 212         hpd.push_back(0.39653931948231);
 213         hpd.push_back(-0.72913209084656);
 214         hpd.push_back(0.46978228740536);
 215         hpd.push_back(0.14390600392911);
 216         hpd.push_back(-0.22403618499417);
 217         hpd.push_back(-0.07130921926705);
 218         hpd.push_back(0.08061260915107);
 219         hpd.push_back(0.03802993693503);
 220         hpd.push_back(-0.01657454163102);
 221         hpd.push_back(-0.01255099855601);
 222         hpd.push_back(0.00042957797300);
 223         hpd.push_back(0.00180164070400);
 224         hpd.push_back(0.00035371380000);
 225         flength = 14;
 226         break;
 227
 228     case Daubechies_8:
 229         lpd.push_back(-0.00011747678400);
 230         lpd.push_back(0.00067544940600);
 231         lpd.push_back(-0.00039174037300);
 232         lpd.push_back(-0.00487035299301);
 233         lpd.push_back(0.00874609404702);
 234         lpd.push_back(0.01398102791702);
 235         lpd.push_back(-0.04408825393106);
 236         lpd.push_back(-0.01736930100202);
 237         lpd.push_back(0.12874742662019);
 238         lpd.push_back(0.00047248457400);
 239         lpd.push_back(-0.28401554296243);
 240         lpd.push_back(-0.01582910525602);
 241         lpd.push_back(0.58535468365487);
 242         lpd.push_back(0.67563073629801);
 243         lpd.push_back(0.31287159091447);
 244         lpd.push_back(0.05441584224308);
 245         hpd.push_back(-0.05441584224308);
 246         hpd.push_back(0.31287159091447);
 247         hpd.push_back(-0.67563073629801);
 248         hpd.push_back(0.58535468365487);
 249         hpd.push_back(0.01582910525602);
 250         hpd.push_back(-0.28401554296243);
 251         hpd.push_back(-0.00047248457400);
 252         hpd.push_back(0.12874742662019);
 253         hpd.push_back(0.01736930100202);
 254         hpd.push_back(-0.04408825393106);
 255         hpd.push_back(-0.01398102791702);
 256         hpd.push_back(0.00874609404702);
 257         hpd.push_back(0.00487035299301);
 258         hpd.push_back(-0.00039174037300);
 259         hpd.push_back(-0.00067544940600);
 260         hpd.push_back(-0.00011747678400);
 261         flength = 16;
 262         break;
 263
 264     case Daubechies_9:
 265         lpd.push_back(0.00003934732000);
 266         lpd.push_back(-0.00025196318900);
 267         lpd.push_back(0.00023038576400);
 268         lpd.push_back(0.00184764688296);
 269         lpd.push_back(-0.00428150368190);
 270         lpd.push_back(-0.00472320475789);
 271         lpd.push_back(0.02236166212352);
 272         lpd.push_back(0.00025094711499);
 273         lpd.push_back(-0.06763282905952);
 274         lpd.push_back(0.03072568147832);
 275         lpd.push_back(0.14854074933476);
 276         lpd.push_back(-0.09684078322088);
 277         lpd.push_back(-0.29327378327259);
 278         lpd.push_back(0.13319738582209);
 279         lpd.push_back(0.65728807803664);
 280         lpd.push_back(0.60482312367678);
 281         lpd.push_back(0.24383467463767);
 282         lpd.push_back(0.03807794736317);
 283         hpd.push_back(-0.03807794736317);
 284         hpd.push_back(0.24383467463767);
 285         hpd.push_back(-0.60482312367678);
 286         hpd.push_back(0.65728807803664);
 287         hpd.push_back(-0.13319738582209);
 288         hpd.push_back(-0.29327378327259);
 289         hpd.push_back(0.09684078322088);
 290         hpd.push_back(0.14854074933476);
 291         hpd.push_back(-0.03072568147832);
 292         hpd.push_back(-0.06763282905952);
 293         hpd.push_back(-0.00025094711499);
 294         hpd.push_back(0.02236166212352);
 295         hpd.push_back(0.00472320475789);
 296         hpd.push_back(-0.00428150368190);
 297         hpd.push_back(-0.00184764688296);
 298         hpd.push_back(0.00023038576400);
 299         hpd.push_back(0.00025196318900);
 300         hpd.push_back(0.00003934732000);
 301         flength = 18;
 302         break;
 303
 304     case Daubechies_10:
 305         lpd.push_back(-0.00001326420300);
 306         lpd.push_back(0.00009358867000);
 307         lpd.push_back(-0.00011646685499);
 308         lpd.push_back(-0.00068585669500);
 309         lpd.push_back(0.00199240529499);
 310         lpd.push_back(0.00139535174699);
 311         lpd.push_back(-0.01073317548298);
 312         lpd.push_back(0.00360655356699);
 313         lpd.push_back(0.03321267405893);
 314         lpd.push_back(-0.02945753682195);
 315         lpd.push_back(-0.07139414716586);
 316         lpd.push_back(0.09305736460381);
 317         lpd.push_back(0.12736934033574);
 318         lpd.push_back(-0.19594627437660);
 319         lpd.push_back(-0.24984642432649);
 320         lpd.push_back(0.28117234366043);
 321         lpd.push_back(0.68845903945259);
 322         lpd.push_back(0.52720118893092);
 323         lpd.push_back(0.18817680007762);
 324         lpd.push_back(0.02667005790095);
 325         hpd.push_back(-0.02667005790095);
 326         hpd.push_back(0.18817680007762);
 327         hpd.push_back(-0.52720118893092);
 328         hpd.push_back(0.68845903945259);
 329         hpd.push_back(-0.28117234366043);
 330         hpd.push_back(-0.24984642432649);
 331         hpd.push_back(0.19594627437660);
 332         hpd.push_back(0.12736934033574);
 333         hpd.push_back(-0.09305736460381);
 334         hpd.push_back(-0.07139414716586);
 335         hpd.push_back(0.02945753682195);
 336         hpd.push_back(0.03321267405893);
 337         hpd.push_back(-0.00360655356699);
 338         hpd.push_back(-0.01073317548298);
 339         hpd.push_back(-0.00139535174699);
 340         hpd.push_back(0.00199240529499);
 341         hpd.push_back(0.00068585669500);
 342         hpd.push_back(-0.00011646685499);
 343         hpd.push_back(-0.00009358867000);
 344         hpd.push_back(-0.00001326420300);
 345         flength = 20;
 346         break;
 347
 348     case Daubechies_20:
 349         lpd.push_back(-0.00000000029988);
 350         lpd.push_back(0.00000000405613);
 351         lpd.push_back(-0.00000001814843);
 352         lpd.push_back(0.00000000020143);
 353         lpd.push_back(0.00000026339242);
 354         lpd.push_back(-0.00000068470796);
 355         lpd.push_back(-0.00000101199401);
 356         lpd.push_back(0.00000724124829);
 357         lpd.push_back(-0.00000437614386);
 358         lpd.push_back(-0.00003710586183);
 359         lpd.push_back(0.00006774280828);
 360         lpd.push_back(0.00010153288973);
 361         lpd.push_back(-0.00038510474870);
 362         lpd.push_back(-0.00005349759845);
 363         lpd.push_back(0.00139255961930);
 364         lpd.push_back(-0.00083156217284);
 365         lpd.push_back(-0.00358149425960);
 366         lpd.push_back(0.00442054238705);
 367         lpd.push_back(0.00672162730228);
 368         lpd.push_back(-0.01381052613705);
 369         lpd.push_back(-0.00878932492387);
 370         lpd.push_back(0.03229429953057);
 371         lpd.push_back(0.00587468181179);
 372         lpd.push_back(-0.06172289962416);
 373         lpd.push_back(0.00563224685793);
 374         lpd.push_back(0.10229171917503);
 375         lpd.push_back(-0.02471682733721);
 376         lpd.push_back(-0.15545875070532);
 377         lpd.push_back(0.03985024645854);
 378         lpd.push_back(0.22829105081947);
 379         lpd.push_back(-0.01672708830868);
 380         lpd.push_back(-0.32678680043239);
 381         lpd.push_back(-0.13921208801080);
 382         lpd.push_back(0.36150229873767);
 383         lpd.push_back(0.61049323893578);
 384         lpd.push_back(0.47269618530872);
 385         lpd.push_back(0.21994211355038);
 386         lpd.push_back(0.06342378045879);
 387         lpd.push_back(0.01054939462490);
 388         lpd.push_back(0.00077995361366);
 389         hpd.push_back(-0.00077995361366);
 390         hpd.push_back(0.01054939462490);
 391         hpd.push_back(-0.06342378045879);
 392         hpd.push_back(0.21994211355038);
 393         hpd.push_back(-0.47269618530872);
 394         hpd.push_back(0.61049323893578);
 395         hpd.push_back(-0.36150229873767);
 396         hpd.push_back(-0.13921208801080);
 397         hpd.push_back(0.32678680043239);
 398         hpd.push_back(-0.01672708830868);
 399         hpd.push_back(-0.22829105081947);
 400         hpd.push_back(0.03985024645854);
 401         hpd.push_back(0.15545875070532);
 402         hpd.push_back(-0.02471682733721);
 403         hpd.push_back(-0.10229171917503);
 404         hpd.push_back(0.00563224685793);
 405         hpd.push_back(0.06172289962416);
 406         hpd.push_back(0.00587468181179);
 407         hpd.push_back(-0.03229429953057);
 408         hpd.push_back(-0.00878932492387);
 409         hpd.push_back(0.01381052613705);
 410         hpd.push_back(0.00672162730228);
 411         hpd.push_back(-0.00442054238705);
 412         hpd.push_back(-0.00358149425960);
 413         hpd.push_back(0.00083156217284);
 414         hpd.push_back(0.00139255961930);
 415         hpd.push_back(0.00005349759845);
 416         hpd.push_back(-0.00038510474870);
 417         hpd.push_back(-0.00010153288973);
 418         hpd.push_back(0.00006774280828);
 419         hpd.push_back(0.00003710586183);
 420         hpd.push_back(-0.00000437614386);
 421         hpd.push_back(-0.00000724124829);
 422         hpd.push_back(-0.00000101199401);
 423         hpd.push_back(0.00000068470796);
 424         hpd.push_back(0.00000026339242);
 425         hpd.push_back(-0.00000000020143);
 426         hpd.push_back(-0.00000001814843);
 427         hpd.push_back(-0.00000000405613);
 428         hpd.push_back(-0.00000000029988);
 429         flength = 40;
 430         break;
 431
 432     case Daubechies_40:
 433         lpd.push_back(-0.00000000000000);
 434         lpd.push_back(0.00000000000000);
 435         lpd.push_back(-0.00000000000000);
 436         lpd.push_back(0.00000000000000);
 437         lpd.push_back(-0.00000000000000);
 438         lpd.push_back(-0.00000000000001);
 439         lpd.push_back(0.00000000000003);
 440         lpd.push_back(-0.00000000000001);
 441         lpd.push_back(-0.00000000000037);
 442         lpd.push_back(0.00000000000121);
 443         lpd.push_back(0.00000000000108);
 444         lpd.push_back(-0.00000000001441);
 445         lpd.push_back(0.00000000001995);
 446         lpd.push_back(0.00000000008134);
 447         lpd.push_back(-0.00000000029288);
 448         lpd.push_back(-0.00000000009963);
 449         lpd.push_back(0.00000000205938);
 450         lpd.push_back(-0.00000000227332);
 451         lpd.push_back(-0.00000000873967);
 452         lpd.push_back(0.00000002288390);
 453         lpd.push_back(0.00000001741059);
 454         lpd.push_back(-0.00000012745008);
 455         lpd.push_back(0.00000005361186);
 456         lpd.push_back(0.00000048834273);
 457         lpd.push_back(-0.00000066863337);
 458         lpd.push_back(-0.00000128043283);
 459         lpd.push_back(0.00000350778202);
 460         lpd.push_back(0.00000165493877);
 461         lpd.push_back(-0.00001288842174);
 462         lpd.push_back(0.00000406557792);
 463         lpd.push_back(0.00003629399945);
 464         lpd.push_back(-0.00003566632371);
 465         lpd.push_back(-0.00007878802614);
 466         lpd.push_back(0.00014251943335);
 467         lpd.push_back(0.00012109867291);
 468         lpd.push_back(-0.00041818330353);
 469         lpd.push_back(-0.00007088554491);
 470         lpd.push_back(0.00099285797789);
 471         lpd.push_back(-0.00031193527356);
 472         lpd.push_back(-0.00197480748041);
 473         lpd.push_back(0.00148388984367);
 474         lpd.push_back(0.00332132500594);
 475         lpd.push_back(-0.00412431015737);
 476         lpd.push_back(-0.00467322831674);
 477         lpd.push_back(0.00901902126682);
 478         lpd.push_back(0.00524573780794);
 479         lpd.push_back(-0.01682072694998);
 480         lpd.push_back(-0.00388701203892);
 481         lpd.push_back(0.02781129406516);
 482         lpd.push_back(-0.00059370667303);
 483         lpd.push_back(-0.04179364874023);
 484         lpd.push_back(0.00895082555981);
 485         lpd.push_back(0.05834125954654);
 486         lpd.push_back(-0.02094317414735);
 487         lpd.push_back(-0.07729538082012);
 488         lpd.push_back(0.03492832601473);
 489         lpd.push_back(0.09941959711973);
 490         lpd.push_back(-0.04741061635274);
 491         lpd.push_back(-0.12673146574763);
 492         lpd.push_back(0.05216571236437);
 493         lpd.push_back(0.16185493481917);
 494         lpd.push_back(-0.03814703831047);
 495         lpd.push_back(-0.20452405393126);
 496         lpd.push_back(-0.01395065907035);
 497         lpd.push_back(0.24017584511066);
 498         lpd.push_back(0.13119716231525);
 499         lpd.push_back(-0.21027451192673);
 500         lpd.push_back(-0.31275608334174);
 501         lpd.push_back(-0.02817037418845);
 502         lpd.push_back(0.35209361261033);
 503         lpd.push_back(0.51902874766301);
 504         lpd.push_back(0.43815806016110);
 505         lpd.push_back(0.25977619424035);
 506         lpd.push_back(0.11476477757177);
 507         lpd.push_back(0.03842788992792);
 508         lpd.push_back(0.00969841552091);
 509         lpd.push_back(0.00179808920458);
 510         lpd.push_back(0.00023209364055);
 511         lpd.push_back(0.00001869221611);
 512         lpd.push_back(0.00000070900603);
 513         hpd.push_back(-0.00000070900603);
 514         hpd.push_back(0.00001869221611);
 515         hpd.push_back(-0.00023209364055);
 516         hpd.push_back(0.00179808920458);
 517         hpd.push_back(-0.00969841552091);
 518         hpd.push_back(0.03842788992792);
 519         hpd.push_back(-0.11476477757177);
 520         hpd.push_back(0.25977619424035);
 521         hpd.push_back(-0.43815806016110);
 522         hpd.push_back(0.51902874766301);
 523         hpd.push_back(-0.35209361261033);
 524         hpd.push_back(-0.02817037418845);
 525         hpd.push_back(0.31275608334174);
 526         hpd.push_back(-0.21027451192673);
 527         hpd.push_back(-0.13119716231525);
 528         hpd.push_back(0.24017584511066);
 529         hpd.push_back(0.01395065907035);
 530         hpd.push_back(-0.20452405393126);
 531         hpd.push_back(0.03814703831047);
 532         hpd.push_back(0.16185493481917);
 533         hpd.push_back(-0.05216571236437);
 534         hpd.push_back(-0.12673146574763);
 535         hpd.push_back(0.04741061635274);
 536         hpd.push_back(0.09941959711973);
 537         hpd.push_back(-0.03492832601473);
 538         hpd.push_back(-0.07729538082012);
 539         hpd.push_back(0.02094317414735);
 540         hpd.push_back(0.05834125954654);
 541         hpd.push_back(-0.00895082555981);
 542         hpd.push_back(-0.04179364874023);
 543         hpd.push_back(0.00059370667303);
 544         hpd.push_back(0.02781129406516);
 545         hpd.push_back(0.00388701203892);
 546         hpd.push_back(-0.01682072694998);
 547         hpd.push_back(-0.00524573780794);
 548         hpd.push_back(0.00901902126682);
 549         hpd.push_back(0.00467322831674);
 550         hpd.push_back(-0.00412431015737);
 551         hpd.push_back(-0.00332132500594);
 552         hpd.push_back(0.00148388984367);
 553         hpd.push_back(0.00197480748041);
 554         hpd.push_back(-0.00031193527356);
 555         hpd.push_back(-0.00099285797789);
 556         hpd.push_back(-0.00007088554491);
 557         hpd.push_back(0.00041818330353);
 558         hpd.push_back(0.00012109867291);
 559         hpd.push_back(-0.00014251943335);
 560         hpd.push_back(-0.00007878802614);
 561         hpd.push_back(0.00003566632371);
 562         hpd.push_back(0.00003629399945);
 563         hpd.push_back(-0.00000406557792);
 564         hpd.push_back(-0.00001288842174);
 565         hpd.push_back(-0.00000165493877);
 566         hpd.push_back(0.00000350778202);
 567         hpd.push_back(0.00000128043283);
 568         hpd.push_back(-0.00000066863337);
 569         hpd.push_back(-0.00000048834273);
 570         hpd.push_back(0.00000005361186);
 571         hpd.push_back(0.00000012745008);
 572         hpd.push_back(0.00000001741059);
 573         hpd.push_back(-0.00000002288390);
 574         hpd.push_back(-0.00000000873967);
 575         hpd.push_back(0.00000000227332);
 576         hpd.push_back(0.00000000205938);
 577         hpd.push_back(0.00000000009963);
 578         hpd.push_back(-0.00000000029288);
 579         hpd.push_back(-0.00000000008134);
 580         hpd.push_back(0.00000000001995);
 581         hpd.push_back(0.00000000001441);
 582         hpd.push_back(0.00000000000108);
 583         hpd.push_back(-0.00000000000121);
 584         hpd.push_back(-0.00000000000037);
 585         hpd.push_back(0.00000000000001);
 586         hpd.push_back(0.00000000000003);
 587         hpd.push_back(0.00000000000001);
 588         hpd.push_back(-0.00000000000000);
 589         hpd.push_back(-0.00000000000000);
 590         hpd.push_back(-0.00000000000000);
 591         hpd.push_back(-0.00000000000000);
 592         hpd.push_back(-0.00000000000000);
 593         flength = 80;
 594         break;
 595
 596     case Symlet_2:
 597         lpd.push_back(-0.12940952255092);
 598         lpd.push_back(0.22414386804186);
 599         lpd.push_back(0.83651630373747);
 600         lpd.push_back(0.48296291314469);
 601         hpd.push_back(-0.48296291314469);
 602         hpd.push_back(0.83651630373747);
 603         hpd.push_back(-0.22414386804186);
 604         hpd.push_back(-0.12940952255092);
 605         flength = 4;
 606         break;
 607
 608     case Symlet_3:
 609         lpd.push_back(0.03522629188210);
 610         lpd.push_back(-0.08544127388224);
 611         lpd.push_back(-0.13501102001039);
 612         lpd.push_back(0.45987750211933);
 613         lpd.push_back(0.80689150931334);
 614         lpd.push_back(0.33267055295096);
 615         hpd.push_back(-0.33267055295096);
 616         hpd.push_back(0.80689150931334);
 617         hpd.push_back(-0.45987750211933);
 618         hpd.push_back(-0.13501102001039);
 619         hpd.push_back(0.08544127388224);
 620         hpd.push_back(0.03522629188210);
 621         flength = 6;
 622         break;
 623
 624     case Symlet_4:
 625         lpd.push_back(-0.07576571478927);
 626         lpd.push_back(-0.02963552764600);
 627         lpd.push_back(0.49761866763202);
 628         lpd.push_back(0.80373875180592);
 629         lpd.push_back(0.29785779560528);
 630         lpd.push_back(-0.09921954357685);
 631         lpd.push_back(-0.01260396726204);
 632         lpd.push_back(0.03222310060404);
 633         hpd.push_back(-0.03222310060404);
 634         hpd.push_back(-0.01260396726204);
 635         hpd.push_back(0.09921954357685);
 636         hpd.push_back(0.29785779560528);
 637         hpd.push_back(-0.80373875180592);
 638         hpd.push_back(0.49761866763202);
 639         hpd.push_back(0.02963552764600);
 640         hpd.push_back(-0.07576571478927);
 641         flength = 8;
 642         break;
 643
 644     case Symlet_5:
 645         lpd.push_back(0.02733306834508);
 646         lpd.push_back(0.02951949092577);
 647         lpd.push_back(-0.03913424930238);
 648         lpd.push_back(0.19939753397739);
 649         lpd.push_back(0.72340769040242);
 650         lpd.push_back(0.63397896345821);
 651         lpd.push_back(0.01660210576452);
 652         lpd.push_back(-0.17532808990845);
 653         lpd.push_back(-0.02110183402476);
 654         lpd.push_back(0.01953888273529);
 655         hpd.push_back(-0.01953888273529);
 656         hpd.push_back(-0.02110183402476);
 657         hpd.push_back(0.17532808990845);
 658         hpd.push_back(0.01660210576452);
 659         hpd.push_back(-0.63397896345821);
 660         hpd.push_back(0.72340769040242);
 661         hpd.push_back(-0.19939753397739);
 662         hpd.push_back(-0.03913424930238);
 663         hpd.push_back(-0.02951949092577);
 664         hpd.push_back(0.02733306834508);
 665         flength = 10;
 666         break;
 667
 668     case Symlet_6:
 669         lpd.push_back(0.01540410932703);
 670         lpd.push_back(0.00349071208422);
 671         lpd.push_back(-0.11799011114819);
 672         lpd.push_back(-0.04831174258563);
 673         lpd.push_back(0.49105594192675);
 674         lpd.push_back(0.78764114103019);
 675         lpd.push_back(0.33792942172762);
 676         lpd.push_back(-0.07263752278646);
 677         lpd.push_back(-0.02106029251230);
 678         lpd.push_back(0.04472490177067);
 679         lpd.push_back(0.00176771186424);
 680         lpd.push_back(-0.00780070832503);
 681         hpd.push_back(0.00780070832503);
 682         hpd.push_back(0.00176771186424);
 683         hpd.push_back(-0.04472490177067);
 684         hpd.push_back(-0.02106029251230);
 685         hpd.push_back(0.07263752278646);
 686         hpd.push_back(0.33792942172762);
 687         hpd.push_back(-0.78764114103019);
 688         hpd.push_back(0.49105594192675);
 689         hpd.push_back(0.04831174258563);
 690         hpd.push_back(-0.11799011114819);
 691         hpd.push_back(-0.00349071208422);
 692         hpd.push_back(0.01540410932703);
 693         flength = 12;
 694         break;
 695
 696     case Symlet_7:
 697         lpd.push_back(0.00268181456826);
 698         lpd.push_back(-0.00104738488868);
 699         lpd.push_back(-0.01263630340325);
 700         lpd.push_back(0.03051551316596);
 701         lpd.push_back(0.06789269350137);
 702         lpd.push_back(-0.04955283493713);
 703         lpd.push_back(0.01744125508686);
 704         lpd.push_back(0.53610191709176);
 705         lpd.push_back(0.76776431700316);
 706         lpd.push_back(0.28862963175151);
 707         lpd.push_back(-0.14004724044296);
 708         lpd.push_back(-0.10780823770382);
 709         lpd.push_back(0.00401024487153);
 710         lpd.push_back(0.01026817670851);
 711         hpd.push_back(-0.01026817670851);
 712         hpd.push_back(0.00401024487153);
 713         hpd.push_back(0.10780823770382);
 714         hpd.push_back(-0.14004724044296);
 715         hpd.push_back(-0.28862963175151);
 716         hpd.push_back(0.76776431700316);
 717         hpd.push_back(-0.53610191709176);
 718         hpd.push_back(0.01744125508686);
 719         hpd.push_back(0.04955283493713);
 720         hpd.push_back(0.06789269350137);
 721         hpd.push_back(-0.03051551316596);
 722         hpd.push_back(-0.01263630340325);
 723         hpd.push_back(0.00104738488868);
 724         hpd.push_back(0.00268181456826);
 725         flength = 14;
 726         break;
 727
 728     case Symlet_8:
 729         lpd.push_back(-0.00338241595101);
 730         lpd.push_back(-0.00054213233179);
 731         lpd.push_back(0.03169508781149);
 732         lpd.push_back(0.00760748732492);
 733         lpd.push_back(-0.14329423835081);
 734         lpd.push_back(-0.06127335906766);
 735         lpd.push_back(0.48135965125837);
 736         lpd.push_back(0.77718575170052);
 737         lpd.push_back(0.36444189483533);
 738         lpd.push_back(-0.05194583810771);
 739         lpd.push_back(-0.02721902991706);
 740         lpd.push_back(0.04913717967361);
 741         lpd.push_back(0.00380875201389);
 742         lpd.push_back(-0.01495225833705);
 743         lpd.push_back(-0.00030292051472);
 744         lpd.push_back(0.00188995033276);
 745         hpd.push_back(-0.00188995033276);
 746         hpd.push_back(-0.00030292051472);
 747         hpd.push_back(0.01495225833705);
 748         hpd.push_back(0.00380875201389);
 749         hpd.push_back(-0.04913717967361);
 750         hpd.push_back(-0.02721902991706);
 751         hpd.push_back(0.05194583810771);
 752         hpd.push_back(0.36444189483533);
 753         hpd.push_back(-0.77718575170052);
 754         hpd.push_back(0.48135965125837);
 755         hpd.push_back(0.06127335906766);
 756         hpd.push_back(-0.14329423835081);
 757         hpd.push_back(-0.00760748732492);
 758         hpd.push_back(0.03169508781149);
 759         hpd.push_back(0.00054213233179);
 760         hpd.push_back(-0.00338241595101);
 761         flength = 16;
 762         break;
 763
 764     case Symlet_9:
 765         lpd.push_back(0.00140091552591);
 766         lpd.push_back(0.00061978088899);
 767         lpd.push_back(-0.01327196778182);
 768         lpd.push_back(-0.01152821020768);
 769         lpd.push_back(0.03022487885828);
 770         lpd.push_back(0.00058346274612);
 771         lpd.push_back(-0.05456895843083);
 772         lpd.push_back(0.23876091460730);
 773         lpd.push_back(0.71789708276441);
 774         lpd.push_back(0.61733844914094);
 775         lpd.push_back(0.03527248803527);
 776         lpd.push_back(-0.19155083129728);
 777         lpd.push_back(-0.01823377077940);
 778         lpd.push_back(0.06207778930289);
 779         lpd.push_back(0.00885926749340);
 780         lpd.push_back(-0.01026406402763);
 781         lpd.push_back(-0.00047315449868);
 782         lpd.push_back(0.00106949003291);
 783         hpd.push_back(-0.00106949003291);
 784         hpd.push_back(-0.00047315449868);
 785         hpd.push_back(0.01026406402763);
 786         hpd.push_back(0.00885926749340);
 787         hpd.push_back(-0.06207778930289);
 788         hpd.push_back(-0.01823377077940);
 789         hpd.push_back(0.19155083129728);
 790         hpd.push_back(0.03527248803527);
 791         hpd.push_back(-0.61733844914094);
 792         hpd.push_back(0.71789708276441);
 793         hpd.push_back(-0.23876091460730);
 794         hpd.push_back(-0.05456895843083);
 795         hpd.push_back(-0.00058346274612);
 796         hpd.push_back(0.03022487885828);
 797         hpd.push_back(0.01152821020768);
 798         hpd.push_back(-0.01327196778182);
 799         hpd.push_back(-0.00061978088899);
 800         hpd.push_back(0.00140091552591);
 801         flength = 18;
 802         break;
 803
 804     case Symlet_10:
 805         lpd.push_back(0.00077015980911);
 806         lpd.push_back(0.00009563267072);
 807         lpd.push_back(-0.00864129927702);
 808         lpd.push_back(-0.00146538258130);
 809         lpd.push_back(0.04592723923109);
 810         lpd.push_back(0.01160989390371);
 811         lpd.push_back(-0.15949427888491);
 812         lpd.push_back(-0.07088053578323);
 813         lpd.push_back(0.47169066693845);
 814         lpd.push_back(0.76951003702110);
 815         lpd.push_back(0.38382676106707);
 816         lpd.push_back(-0.03553674047383);
 817         lpd.push_back(-0.03199005688243);
 818         lpd.push_back(0.04999497207737);
 819         lpd.push_back(0.00576491203358);
 820         lpd.push_back(-0.02035493981231);
 821         lpd.push_back(-0.00080435893202);
 822         lpd.push_back(0.00459317358531);
 823         lpd.push_back(0.00005703608362);
 824         lpd.push_back(-0.00045932942100);
 825         hpd.push_back(0.00045932942100);
 826         hpd.push_back(0.00005703608362);
 827         hpd.push_back(-0.00459317358531);
 828         hpd.push_back(-0.00080435893202);
 829         hpd.push_back(0.02035493981231);
 830         hpd.push_back(0.00576491203358);
 831         hpd.push_back(-0.04999497207737);
 832         hpd.push_back(-0.03199005688243);
 833         hpd.push_back(0.03553674047383);
 834         hpd.push_back(0.38382676106707);
 835         hpd.push_back(-0.76951003702110);
 836         hpd.push_back(0.47169066693845);
 837         hpd.push_back(0.07088053578323);
 838         hpd.push_back(-0.15949427888491);
 839         hpd.push_back(-0.01160989390371);
 840         hpd.push_back(0.04592723923109);
 841         hpd.push_back(0.00146538258130);
 842         hpd.push_back(-0.00864129927702);
 843         hpd.push_back(-0.00009563267072);
 844         hpd.push_back(0.00077015980911);
 845         flength = 20;
 846         break;
 847
 848     case Symlet_20:
 849         lpd.push_back(0.00000036955375);
 850         lpd.push_back(-0.00000019015676);
 851         lpd.push_back(-0.00000791936141);
 852         lpd.push_back(0.00000302566606);
 853         lpd.push_back(0.00007992967836);
 854         lpd.push_back(-0.00001928412301);
 855         lpd.push_back(-0.00049473109157);
 856         lpd.push_back(0.00007215991188);
 857         lpd.push_back(0.00208899470819);
 858         lpd.push_back(-0.00030526283181);
 859         lpd.push_back(-0.00660658579912);
 860         lpd.push_back(0.00142308735944);
 861         lpd.push_back(0.01700404902335);
 862         lpd.push_back(-0.00331385738375);
 863         lpd.push_back(-0.03162943714501);
 864         lpd.push_back(0.00812322835637);
 865         lpd.push_back(0.02557934951027);
 866         lpd.push_back(-0.07899434492693);
 867         lpd.push_back(-0.02981936887758);
 868         lpd.push_back(0.40583144435233);
 869         lpd.push_back(0.75116272842520);
 870         lpd.push_back(0.47199147510012);
 871         lpd.push_back(-0.05108834292497);
 872         lpd.push_back(-0.16057829841831);
 873         lpd.push_back(0.03625095165278);
 874         lpd.push_back(0.08891966802764);
 875         lpd.push_back(-0.00684370196580);
 876         lpd.push_back(-0.03537333675714);
 877         lpd.push_back(0.00193859706711);
 878         lpd.push_back(0.01215704094879);
 879         lpd.push_back(-0.00061112638583);
 880         lpd.push_back(-0.00347164780287);
 881         lpd.push_back(0.00012544091723);
 882         lpd.push_back(0.00074761085979);
 883         lpd.push_back(-0.00002661555034);
 884         lpd.push_back(-0.00011739133516);
 885         lpd.push_back(0.00000452542221);
 886         lpd.push_back(0.00001228725278);
 887         lpd.push_back(-0.00000032567026);
 888         lpd.push_back(-0.00000063291290);
 889         hpd.push_back(0.00000063291290);
 890         hpd.push_back(-0.00000032567026);
 891         hpd.push_back(-0.00001228725278);
 892         hpd.push_back(0.00000452542221);
 893         hpd.push_back(0.00011739133516);
 894         hpd.push_back(-0.00002661555034);
 895         hpd.push_back(-0.00074761085979);
 896         hpd.push_back(0.00012544091723);
 897         hpd.push_back(0.00347164780287);
 898         hpd.push_back(-0.00061112638583);
 899         hpd.push_back(-0.01215704094879);
 900         hpd.push_back(0.00193859706711);
 901         hpd.push_back(0.03537333675714);
 902         hpd.push_back(-0.00684370196580);
 903         hpd.push_back(-0.08891966802764);
 904         hpd.push_back(0.03625095165278);
 905         hpd.push_back(0.16057829841831);
 906         hpd.push_back(-0.05108834292497);
 907         hpd.push_back(-0.47199147510012);
 908         hpd.push_back(0.75116272842520);
 909         hpd.push_back(-0.40583144435233);
 910         hpd.push_back(-0.02981936887758);
 911         hpd.push_back(0.07899434492693);
 912         hpd.push_back(0.02557934951027);
 913         hpd.push_back(-0.00812322835637);
 914         hpd.push_back(-0.03162943714501);
 915         hpd.push_back(0.00331385738375);
 916         hpd.push_back(0.01700404902335);
 917         hpd.push_back(-0.00142308735944);
 918         hpd.push_back(-0.00660658579912);
 919         hpd.push_back(0.00030526283181);
 920         hpd.push_back(0.00208899470819);
 921         hpd.push_back(-0.00007215991188);
 922         hpd.push_back(-0.00049473109157);
 923         hpd.push_back(0.00001928412301);
 924         hpd.push_back(0.00007992967836);
 925         hpd.push_back(-0.00000302566606);
 926         hpd.push_back(-0.00000791936141);
 927         hpd.push_back(0.00000019015676);
 928         hpd.push_back(0.00000036955375);
 929         flength = 40;
 930         break;
 931
 932     case Symlet_30:
 933         lpd.push_back(-0.00000000032641);
 934         lpd.push_back(0.00000000014934);
 935         lpd.push_back(0.00000001014810);
 936         lpd.push_back(-0.00000000380480);
 937         lpd.push_back(-0.00000014986093);
 938         lpd.push_back(0.00000005465704);
 939         lpd.push_back(0.00000142300535);
 940         lpd.push_back(-0.00000055432777);
 941         lpd.push_back(-0.00000989667666);
 942         lpd.push_back(0.00000402586608);
 943         lpd.push_back(0.00005384705071);
 944         lpd.push_back(-0.00002120203482);
 945         lpd.push_back(-0.00023724381672);
 946         lpd.push_back(0.00008209438737);
 947         lpd.push_back(0.00086145325795);
 948         lpd.push_back(-0.00023156681831);
 949         lpd.push_back(-0.00258964866287);
 950         lpd.push_back(0.00046642571155);
 951         lpd.push_back(0.00645416706920);
 952         lpd.push_back(-0.00051919282891);
 953         lpd.push_back(-0.01284118837655);
 954         lpd.push_back(0.00037212345706);
 955         lpd.push_back(0.01829725025850);
 956         lpd.push_back(-0.00795200356952);
 957         lpd.push_back(-0.02521794435506);
 958         lpd.push_back(0.02965858822312);
 959         lpd.push_back(0.02927100294045);
 960         lpd.push_back(-0.09630227559219);
 961         lpd.push_back(-0.03849155121162);
 962         lpd.push_back(0.40072250790936);
 963         lpd.push_back(0.74269376814242);
 964         lpd.push_back(0.48469805553706);
 965         lpd.push_back(-0.02843778080778);
 966         lpd.push_back(-0.16204679936904);
 967         lpd.push_back(0.02619852204894);
 968         lpd.push_back(0.08910160832786);
 969         lpd.push_back(-0.01523693234758);
 970         lpd.push_back(-0.04889669606591);
 971         lpd.push_back(0.00760939409040);
 972         lpd.push_back(0.02556854558808);
 973         lpd.push_back(-0.00162932098036);
 974         lpd.push_back(-0.01038876672659);
 975         lpd.push_back(0.00054379809563);
 976         lpd.push_back(0.00371380689803);
 977         lpd.push_back(-0.00027628401612);
 978         lpd.push_back(-0.00119692132405);
 979         lpd.push_back(0.00011735865251);
 980         lpd.push_back(0.00034385220669);
 981         lpd.push_back(-0.00003293175202);
 982         lpd.push_back(-0.00008272690387);
 983         lpd.push_back(0.00000647382532);
 984         lpd.push_back(0.00001589700426);
 985         lpd.push_back(-0.00000089828423);
 986         lpd.push_back(-0.00000232354967);
 987         lpd.push_back(0.00000008995011);
 988         lpd.push_back(0.00000024412960);
 989         lpd.push_back(-0.00000000612778);
 990         lpd.push_back(-0.00000001650488);
 991         lpd.push_back(0.00000000025363);
 992         lpd.push_back(0.00000000055439);
 993         hpd.push_back(-0.00000000055439);
 994         hpd.push_back(0.00000000025363);
 995         hpd.push_back(0.00000001650488);
 996         hpd.push_back(-0.00000000612778);
 997         hpd.push_back(-0.00000024412960);
 998         hpd.push_back(0.00000008995011);
 999         hpd.push_back(0.00000232354967);
1000         hpd.push_back(-0.00000089828423);
1001         hpd.push_back(-0.00001589700426);
1002         hpd.push_back(0.00000647382532);
1003         hpd.push_back(0.00008272690387);
1004         hpd.push_back(-0.00003293175202);
1005         hpd.push_back(-0.00034385220669);
1006         hpd.push_back(0.00011735865251);
1007         hpd.push_back(0.00119692132405);
1008         hpd.push_back(-0.00027628401612);
1009         hpd.push_back(-0.00371380689803);
1010         hpd.push_back(0.00054379809563);
1011         hpd.push_back(0.01038876672659);
1012         hpd.push_back(-0.00162932098036);
1013         hpd.push_back(-0.02556854558808);
1014         hpd.push_back(0.00760939409040);
1015         hpd.push_back(0.04889669606591);
1016         hpd.push_back(-0.01523693234758);
1017         hpd.push_back(-0.08910160832786);
1018         hpd.push_back(0.02619852204894);
1019         hpd.push_back(0.16204679936904);
1020         hpd.push_back(-0.02843778080778);
1021         hpd.push_back(-0.48469805553706);
1022         hpd.push_back(0.74269376814242);
1023         hpd.push_back(-0.40072250790936);
1024         hpd.push_back(-0.03849155121162);
1025         hpd.push_back(0.09630227559219);
1026         hpd.push_back(0.02927100294045);
1027         hpd.push_back(-0.02965858822312);
1028         hpd.push_back(-0.02521794435506);
1029         hpd.push_back(0.00795200356952);
1030         hpd.push_back(0.01829725025850);
1031         hpd.push_back(-0.00037212345706);
1032         hpd.push_back(-0.01284118837655);
1033         hpd.push_back(0.00051919282891);
1034         hpd.push_back(0.00645416706920);
1035         hpd.push_back(-0.00046642571155);
1036         hpd.push_back(-0.00258964866287);
1037         hpd.push_back(0.00023156681831);
1038         hpd.push_back(0.00086145325795);
1039         hpd.push_back(-0.00008209438737);
1040         hpd.push_back(-0.00023724381672);
1041         hpd.push_back(0.00002120203482);
1042         hpd.push_back(0.00005384705071);
1043         hpd.push_back(-0.00000402586608);
1044         hpd.push_back(-0.00000989667666);
1045         hpd.push_back(0.00000055432777);
1046         hpd.push_back(0.00000142300535);
1047         hpd.push_back(-0.00000005465704);
1048         hpd.push_back(-0.00000014986093);
1049         hpd.push_back(0.00000000380480);
1050         hpd.push_back(0.00000001014810);
1051         hpd.push_back(-0.00000000014934);
1052         hpd.push_back(-0.00000000032641);
1053         flength = 60;
1054         break;
1055
1056     case Coiflet_1:
1057         lpd.push_back(-0.01565572813546);
1058         lpd.push_back(-0.07273261951285);
1059         lpd.push_back(0.38486484686420);
1060         lpd.push_back(0.85257202021226);
1061         lpd.push_back(0.33789766245781);
1062         lpd.push_back(-0.07273261951285);
1063         hpd.push_back(0.07273261951285);
1064         hpd.push_back(0.33789766245781);
1065         hpd.push_back(-0.85257202021226);
1066         hpd.push_back(0.38486484686420);
1067         hpd.push_back(0.07273261951285);
1068         hpd.push_back(-0.01565572813546);
1069         flength = 6;
1070         break;
1071
1072     case Coiflet_2:
1073         lpd.push_back(-0.00072054944536);
1074         lpd.push_back(-0.00182320887070);
1075         lpd.push_back(0.00561143481939);
1076         lpd.push_back(0.02368017194633);
1077         lpd.push_back(-0.05943441864646);
1078         lpd.push_back(-0.07648859907831);
1079         lpd.push_back(0.41700518442169);
1080         lpd.push_back(0.81272363544554);
1081         lpd.push_back(0.38611006682116);
1082         lpd.push_back(-0.06737255472196);
1083         lpd.push_back(-0.04146493678176);
1084         lpd.push_back(0.01638733646352);
1085         hpd.push_back(-0.01638733646352);
1086         hpd.push_back(-0.04146493678176);
1087         hpd.push_back(0.06737255472196);
1088         hpd.push_back(0.38611006682116);
1089         hpd.push_back(-0.81272363544554);
1090         hpd.push_back(0.41700518442169);
1091         hpd.push_back(0.07648859907831);
1092         hpd.push_back(-0.05943441864646);
1093         hpd.push_back(-0.02368017194633);
1094         hpd.push_back(0.00561143481939);
1095         hpd.push_back(0.00182320887070);
1096         hpd.push_back(-0.00072054944536);
1097         flength = 12;
1098         break;
1099
1100     case Coiflet_3:
1101         lpd.push_back(-0.00003459977284);
1102         lpd.push_back(-0.00007098330314);
1103         lpd.push_back(0.00046621696011);
1104         lpd.push_back(0.00111751877089);
1105         lpd.push_back(-0.00257451768875);
1106         lpd.push_back(-0.00900797613666);
1107         lpd.push_back(0.01588054486362);
1108         lpd.push_back(0.03455502757306);
1109         lpd.push_back(-0.08230192710689);
1110         lpd.push_back(-0.07179982161931);
1111         lpd.push_back(0.42848347637762);
1112         lpd.push_back(0.79377722262562);
1113         lpd.push_back(0.40517690240962);
1114         lpd.push_back(-0.06112339000267);
1115         lpd.push_back(-0.06577191128186);
1116         lpd.push_back(0.02345269614184);
1117         lpd.push_back(0.00778259642733);
1118         lpd.push_back(-0.00379351286449);
1119         hpd.push_back(0.00379351286449);
1120         hpd.push_back(0.00778259642733);
1121         hpd.push_back(-0.02345269614184);
1122         hpd.push_back(-0.06577191128186);
1123         hpd.push_back(0.06112339000267);
1124         hpd.push_back(0.40517690240962);
1125         hpd.push_back(-0.79377722262562);
1126         hpd.push_back(0.42848347637762);
1127         hpd.push_back(0.07179982161931);
1128         hpd.push_back(-0.08230192710689);
1129         hpd.push_back(-0.03455502757306);
1130         hpd.push_back(0.01588054486362);
1131         hpd.push_back(0.00900797613666);
1132         hpd.push_back(-0.00257451768875);
1133         hpd.push_back(-0.00111751877089);
1134         hpd.push_back(0.00046621696011);
1135         hpd.push_back(0.00007098330314);
1136         hpd.push_back(-0.00003459977284);
1137         flength = 18;
1138         break;
1139
1140     case Coiflet_4:
1141         lpd.push_back(-0.00000178498500);
1142         lpd.push_back(-0.00000325968024);
1143         lpd.push_back(0.00003122987587);
1144         lpd.push_back(0.00006233903446);
1145         lpd.push_back(-0.00025997455249);
1146         lpd.push_back(-0.00058902075624);
1147         lpd.push_back(0.00126656192930);
1148         lpd.push_back(0.00375143615728);
1149         lpd.push_back(-0.00565828668661);
1150         lpd.push_back(-0.01521173152795);
1151         lpd.push_back(0.02508226184486);
1152         lpd.push_back(0.03933442712334);
1153         lpd.push_back(-0.09622044203399);
1154         lpd.push_back(-0.06662747426343);
1155         lpd.push_back(0.43438605649147);
1156         lpd.push_back(0.78223893092050);
1157         lpd.push_back(0.41530840703043);
1158         lpd.push_back(-0.05607731331675);
1159         lpd.push_back(-0.08126669968088);
1160         lpd.push_back(0.02668230015605);
1161         lpd.push_back(0.01606894396478);
1162         lpd.push_back(-0.00734616632764);
1163         lpd.push_back(-0.00162949201260);
1164         lpd.push_back(0.00089231366858);
1165         hpd.push_back(-0.00089231366858);
1166         hpd.push_back(-0.00162949201260);
1167         hpd.push_back(0.00734616632764);
1168         hpd.push_back(0.01606894396478);
1169         hpd.push_back(-0.02668230015605);
1170         hpd.push_back(-0.08126669968088);
1171         hpd.push_back(0.05607731331675);
1172         hpd.push_back(0.41530840703043);
1173         hpd.push_back(-0.78223893092050);
1174         hpd.push_back(0.43438605649147);
1175         hpd.push_back(0.06662747426343);
1176         hpd.push_back(-0.09622044203399);
1177         hpd.push_back(-0.03933442712334);
1178         hpd.push_back(0.02508226184486);
1179         hpd.push_back(0.01521173152795);
1180         hpd.push_back(-0.00565828668661);
1181         hpd.push_back(-0.00375143615728);
1182         hpd.push_back(0.00126656192930);
1183         hpd.push_back(0.00058902075624);
1184         hpd.push_back(-0.00025997455249);
1185         hpd.push_back(-0.00006233903446);
1186         hpd.push_back(0.00003122987587);
1187         hpd.push_back(0.00000325968024);
1188         hpd.push_back(-0.00000178498500);
1189         flength = 24;
1190         break;
1191
1192     case Coiflet_5:
1193         lpd.push_back(-0.00000009517657);
1194         lpd.push_back(-0.00000016744289);
1195         lpd.push_back(0.00000206376185);
1196         lpd.push_back(0.00000373465518);
1197         lpd.push_back(-0.00002131502681);
1198         lpd.push_back(-0.00004134043227);
1199         lpd.push_back(0.00014054114970);
1200         lpd.push_back(0.00030225958181);
1201         lpd.push_back(-0.00063813134305);
1202         lpd.push_back(-0.00166286370201);
1203         lpd.push_back(0.00243337321266);
1204         lpd.push_back(0.00676418544805);
1205         lpd.push_back(-0.00916423116248);
1206         lpd.push_back(-0.01976177894257);
1207         lpd.push_back(0.03268357426711);
1208         lpd.push_back(0.04128920875018);
1209         lpd.push_back(-0.10557420870334);
1210         lpd.push_back(-0.06203596396290);
1211         lpd.push_back(0.43799162617184);
1212         lpd.push_back(0.77428960365296);
1213         lpd.push_back(0.42156620669085);
1214         lpd.push_back(-0.05204316317624);
1215         lpd.push_back(-0.09192001055970);
1216         lpd.push_back(0.02816802897094);
1217         lpd.push_back(0.02340815678584);
1218         lpd.push_back(-0.01013111751985);
1219         lpd.push_back(-0.00415935878139);
1220         lpd.push_back(0.00217823635811);
1221         lpd.push_back(0.00035858968790);
1222         lpd.push_back(-0.00021208083980);
1223         hpd.push_back(0.00021208083980);
1224         hpd.push_back(0.00035858968790);
1225         hpd.push_back(-0.00217823635811);
1226         hpd.push_back(-0.00415935878139);
1227         hpd.push_back(0.01013111751985);
1228         hpd.push_back(0.02340815678584);
1229         hpd.push_back(-0.02816802897094);
1230         hpd.push_back(-0.09192001055970);
1231         hpd.push_back(0.05204316317624);
1232         hpd.push_back(0.42156620669085);
1233         hpd.push_back(-0.77428960365296);
1234         hpd.push_back(0.43799162617184);
1235         hpd.push_back(0.06203596396290);
1236         hpd.push_back(-0.10557420870334);
1237         hpd.push_back(-0.04128920875018);
1238         hpd.push_back(0.03268357426711);
1239         hpd.push_back(0.01976177894257);
1240         hpd.push_back(-0.00916423116248);
1241         hpd.push_back(-0.00676418544805);
1242         hpd.push_back(0.00243337321266);
1243         hpd.push_back(0.00166286370201);
1244         hpd.push_back(-0.00063813134305);
1245         hpd.push_back(-0.00030225958181);
1246         hpd.push_back(0.00014054114970);
1247         hpd.push_back(0.00004134043227);
1248         hpd.push_back(-0.00002131502681);
1249         hpd.push_back(-0.00000373465518);
1250         hpd.push_back(0.00000206376185);
1251         hpd.push_back(0.00000016744289);
1252         hpd.push_back(-0.00000009517657);
1253         flength = 30;
1254         break;
1255
1256     case Biorthogonal_1_3:
1257         lpd.push_back(-0.08838834764832);
1258         lpd.push_back(0.08838834764832);
1259         lpd.push_back(0.70710678118655);
1260         lpd.push_back(0.70710678118655);
1261         lpd.push_back(0.08838834764832);
1262         lpd.push_back(-0.08838834764832);
1263         hpd.push_back(-0.00000000000000);
1264         hpd.push_back(0.00000000000000);
1265         hpd.push_back(-0.70710678118655);
1266         hpd.push_back(0.70710678118655);
1267         hpd.push_back(-0.00000000000000);
1268         hpd.push_back(0.00000000000000);
1269         flength = 6;
1270         break;
1271
1272     case Biorthogonal_1_5:
1273         lpd.push_back(0.01657281518406);
1274         lpd.push_back(-0.01657281518406);
1275         lpd.push_back(-0.12153397801644);
1276         lpd.push_back(0.12153397801644);
1277         lpd.push_back(0.70710678118655);
1278         lpd.push_back(0.70710678118655);
1279         lpd.push_back(0.12153397801644);
1280         lpd.push_back(-0.12153397801644);
1281         lpd.push_back(-0.01657281518406);
1282         lpd.push_back(0.01657281518406);
1283         hpd.push_back(-0.00000000000000);
1284         hpd.push_back(0.00000000000000);
1285         hpd.push_back(-0.00000000000000);
1286         hpd.push_back(0.00000000000000);
1287         hpd.push_back(-0.70710678118655);
1288         hpd.push_back(0.70710678118655);
1289         hpd.push_back(-0.00000000000000);
1290         hpd.push_back(0.00000000000000);
1291         hpd.push_back(-0.00000000000000);
1292         hpd.push_back(0.00000000000000);
1293         flength = 10;
1294         break;
1295
1296     case Biorthogonal_2_2:
1297         lpd.push_back(0.00000000000000);
1298         lpd.push_back(-0.17677669529664);
1299         lpd.push_back(0.35355339059327);
1300         lpd.push_back(1.06066017177982);
1301         lpd.push_back(0.35355339059327);
1302         lpd.push_back(-0.17677669529664);
1303         hpd.push_back(-0.00000000000000);
1304         hpd.push_back(0.35355339059327);
1305         hpd.push_back(-0.70710678118655);
1306         hpd.push_back(0.35355339059327);
1307         hpd.push_back(-0.00000000000000);
1308         hpd.push_back(0.00000000000000);
1309         flength = 6;
1310         break;
1311
1312     case Biorthogonal_2_4:
1313         lpd.push_back(0.00000000000000);
1314         lpd.push_back(0.03314563036812);
1315         lpd.push_back(-0.06629126073624);
1316         lpd.push_back(-0.17677669529664);
1317         lpd.push_back(0.41984465132951);
1318         lpd.push_back(0.99436891104358);
1319         lpd.push_back(0.41984465132951);
1320         lpd.push_back(-0.17677669529664);
1321         lpd.push_back(-0.06629126073624);
1322         lpd.push_back(0.03314563036812);
1323         hpd.push_back(-0.00000000000000);
1324         hpd.push_back(0.00000000000000);
1325         hpd.push_back(-0.00000000000000);
1326         hpd.push_back(0.35355339059327);
1327         hpd.push_back(-0.70710678118655);
1328         hpd.push_back(0.35355339059327);
1329         hpd.push_back(-0.00000000000000);
1330         hpd.push_back(0.00000000000000);
1331         hpd.push_back(-0.00000000000000);
1332         hpd.push_back(0.00000000000000);
1333         flength = 10;
1334         break;
1335
1336
1337     case Biorthogonal_2_6:
1338         lpd.push_back(0.00000000000000);
1339         lpd.push_back(-0.00690533966002);
1340         lpd.push_back(0.01381067932005);
1341         lpd.push_back(0.04695630968817);
1342         lpd.push_back(-0.10772329869639);
1343         lpd.push_back(-0.16987135563661);
1344         lpd.push_back(0.44746600996961);
1345         lpd.push_back(0.96674755240348);
1346         lpd.push_back(0.44746600996961);
1347         lpd.push_back(-0.16987135563661);
1348         lpd.push_back(-0.10772329869639);
1349         lpd.push_back(0.04695630968817);
1350         lpd.push_back(0.01381067932005);
1351         lpd.push_back(-0.00690533966002);
1352         hpd.push_back(-0.00000000000000);
1353         hpd.push_back(0.00000000000000);
1354         hpd.push_back(-0.00000000000000);
1355         hpd.push_back(0.00000000000000);
1356         hpd.push_back(-0.00000000000000);
1357         hpd.push_back(0.35355339059327);
1358         hpd.push_back(-0.70710678118655);
1359         hpd.push_back(0.35355339059327);
1360         hpd.push_back(-0.00000000000000);
1361         hpd.push_back(0.00000000000000);
1362         hpd.push_back(-0.00000000000000);
1363         hpd.push_back(0.00000000000000);
1364         hpd.push_back(-0.00000000000000);
1365         hpd.push_back(0.00000000000000);
1366         flength = 14;
1367         break;
1368
1369     case Biorthogonal_2_8:
1370         lpd.push_back(0.00000000000000);
1371         lpd.push_back(0.00151054305063);
1372         lpd.push_back(-0.00302108610126);
1373         lpd.push_back(-0.01294751186255);
1374         lpd.push_back(0.02891610982635);
1375         lpd.push_back(0.05299848189069);
1376         lpd.push_back(-0.13491307360774);
1377         lpd.push_back(-0.16382918343409);
1378         lpd.push_back(0.46257144047592);
1379         lpd.push_back(0.95164212189718);
1380         lpd.push_back(0.46257144047592);
1381         lpd.push_back(-0.16382918343409);
1382         lpd.push_back(-0.13491307360774);
1383         lpd.push_back(0.05299848189069);
1384         lpd.push_back(0.02891610982635);
1385         lpd.push_back(-0.01294751186255);
1386         lpd.push_back(-0.00302108610126);
1387         lpd.push_back(0.00151054305063);
1388         hpd.push_back(-0.00000000000000);
1389         hpd.push_back(0.00000000000000);
1390         hpd.push_back(-0.00000000000000);
1391         hpd.push_back(0.00000000000000);
1392         hpd.push_back(-0.00000000000000);
1393         hpd.push_back(0.00000000000000);
1394         hpd.push_back(-0.00000000000000);
1395         hpd.push_back(0.35355339059327);
1396         hpd.push_back(-0.70710678118655);
1397         hpd.push_back(0.35355339059327);
1398         hpd.push_back(-0.00000000000000);
1399         hpd.push_back(0.00000000000000);
1400         hpd.push_back(-0.00000000000000);
1401         hpd.push_back(0.00000000000000);
1402         hpd.push_back(-0.00000000000000);
1403         hpd.push_back(0.00000000000000);
1404         hpd.push_back(-0.00000000000000);
1405         hpd.push_back(0.00000000000000);
1406         flength = 18;
1407         break;
1408
1409     case Biorthogonal_3_1:
1410         lpd.push_back(-0.35355339059327);
1411         lpd.push_back(1.06066017177982);
1412         lpd.push_back(1.06066017177982);
1413         lpd.push_back(-0.35355339059327);
1414         hpd.push_back(-0.17677669529664);
1415         hpd.push_back(0.53033008588991);
1416         hpd.push_back(-0.53033008588991);
1417         hpd.push_back(0.17677669529664);
1418         flength = 4;
1419         break;
1420
1421     case Biorthogonal_3_3:
1422         lpd.push_back(0.06629126073624);
1423         lpd.push_back(-0.19887378220872);
1424         lpd.push_back(-0.15467960838456);
1425         lpd.push_back(0.99436891104358);
1426         lpd.push_back(0.99436891104358);
1427         lpd.push_back(-0.15467960838456);
1428         lpd.push_back(-0.19887378220872);
1429         lpd.push_back(0.06629126073624);
1430         hpd.push_back(-0.00000000000000);
1431         hpd.push_back(0.00000000000000);
1432         hpd.push_back(-0.17677669529664);
1433         hpd.push_back(0.53033008588991);
1434         hpd.push_back(-0.53033008588991);
1435         hpd.push_back(0.17677669529664);
1436         hpd.push_back(-0.00000000000000);
1437         hpd.push_back(0.00000000000000);
1438         flength = 8;
1439         break;
1440
1441     case Biorthogonal_3_5:
1442         lpd.push_back(-0.01381067932005);
1443         lpd.push_back(0.04143203796015);
1444         lpd.push_back(0.05248058141619);
1445         lpd.push_back(-0.26792717880897);
1446         lpd.push_back(-0.07181553246426);
1447         lpd.push_back(0.96674755240348);
1448         lpd.push_back(0.96674755240348);
1449         lpd.push_back(-0.07181553246426);
1450         lpd.push_back(-0.26792717880897);
1451         lpd.push_back(0.05248058141619);
1452         lpd.push_back(0.04143203796015);
1453         lpd.push_back(-0.01381067932005);
1454         hpd.push_back(-0.00000000000000);
1455         hpd.push_back(0.00000000000000);
1456         hpd.push_back(-0.00000000000000);
1457         hpd.push_back(0.00000000000000);
1458         hpd.push_back(-0.17677669529664);
1459         hpd.push_back(0.53033008588991);
1460         hpd.push_back(-0.53033008588991);
1461         hpd.push_back(0.17677669529664);
1462         hpd.push_back(-0.00000000000000);
1463         hpd.push_back(0.00000000000000);
1464         hpd.push_back(-0.00000000000000);
1465         hpd.push_back(0.00000000000000);
1466         flength = 12;
1467         break;
1468
1469     case Biorthogonal_3_7:
1470         lpd.push_back(0.00302108610126);
1471         lpd.push_back(-0.00906325830378);
1472         lpd.push_back(-0.01683176542131);
1473         lpd.push_back(0.07466398507402);
1474         lpd.push_back(0.03133297870736);
1475         lpd.push_back(-0.30115912592284);
1476         lpd.push_back(-0.02649924094535);
1477         lpd.push_back(0.95164212189718);
1478         lpd.push_back(0.95164212189718);
1479         lpd.push_back(-0.02649924094535);
1480         lpd.push_back(-0.30115912592284);
1481         lpd.push_back(0.03133297870736);
1482         lpd.push_back(0.07466398507402);
1483         lpd.push_back(-0.01683176542131);
1484         lpd.push_back(-0.00906325830378);
1485         lpd.push_back(0.00302108610126);
1486         hpd.push_back(-0.00000000000000);
1487         hpd.push_back(0.00000000000000);
1488         hpd.push_back(-0.00000000000000);
1489         hpd.push_back(0.00000000000000);
1490         hpd.push_back(-0.00000000000000);
1491         hpd.push_back(0.00000000000000);
1492         hpd.push_back(-0.17677669529664);
1493         hpd.push_back(0.53033008588991);
1494         hpd.push_back(-0.53033008588991);
1495         hpd.push_back(0.17677669529664);
1496         hpd.push_back(-0.00000000000000);
1497         hpd.push_back(0.00000000000000);
1498         hpd.push_back(-0.00000000000000);
1499         hpd.push_back(0.00000000000000);
1500         hpd.push_back(-0.00000000000000);
1501         hpd.push_back(0.00000000000000);
1502         flength = 16;
1503         break;
1504
1505     case Biorthogonal_3_9:
1506         lpd.push_back(-0.00067974437278);
1507         lpd.push_back(0.00203923311835);
1508         lpd.push_back(0.00506031921961);
1509         lpd.push_back(-0.02061891264111);
1510         lpd.push_back(-0.01411278793018);
1511         lpd.push_back(0.09913478249423);
1512         lpd.push_back(0.01230013626942);
1513         lpd.push_back(-0.32019196836078);
1514         lpd.push_back(0.00205002271157);
1515         lpd.push_back(0.94212570067821);
1516         lpd.push_back(0.94212570067821);
1517         lpd.push_back(0.00205002271157);
1518         lpd.push_back(-0.32019196836078);
1519         lpd.push_back(0.01230013626942);
1520         lpd.push_back(0.09913478249423);
1521         lpd.push_back(-0.01411278793018);
1522         lpd.push_back(-0.02061891264111);
1523         lpd.push_back(0.00506031921961);
1524         lpd.push_back(0.00203923311835);
1525         lpd.push_back(-0.00067974437278);
1526         hpd.push_back(-0.00000000000000);
1527         hpd.push_back(0.00000000000000);
1528         hpd.push_back(-0.00000000000000);
1529         hpd.push_back(0.00000000000000);
1530         hpd.push_back(-0.00000000000000);
1531         hpd.push_back(0.00000000000000);
1532         hpd.push_back(-0.00000000000000);
1533         hpd.push_back(0.00000000000000);
1534         hpd.push_back(-0.17677669529664);
1535         hpd.push_back(0.53033008588991);
1536         hpd.push_back(-0.53033008588991);
1537         hpd.push_back(0.17677669529664);
1538         hpd.push_back(-0.00000000000000);
1539         hpd.push_back(0.00000000000000);
1540         hpd.push_back(-0.00000000000000);
1541         hpd.push_back(0.00000000000000);
1542         hpd.push_back(-0.00000000000000);
1543         hpd.push_back(0.00000000000000);
1544         hpd.push_back(-0.00000000000000);
1545         hpd.push_back(0.00000000000000);
1546         flength = 20;
1547         break;
1548
1549     case Biorthogonal_4_4:
1550         lpd.push_back(0.00000000000000);
1551         lpd.push_back(0.03782845550726);
1552         lpd.push_back(-0.02384946501956);
1553         lpd.push_back(-0.11062440441844);
1554         lpd.push_back(0.37740285561283);
1555         lpd.push_back(0.85269867900889);
1556         lpd.push_back(0.37740285561283);
1557         lpd.push_back(-0.11062440441844);
1558         lpd.push_back(-0.02384946501956);
1559         lpd.push_back(0.03782845550726);
1560         hpd.push_back(-0.00000000000000);
1561         hpd.push_back(-0.06453888262870);
1562         hpd.push_back(0.04068941760916);
1563         hpd.push_back(0.41809227322162);
1564         hpd.push_back(-0.78848561640558);
1565         hpd.push_back(0.41809227322162);
1566         hpd.push_back(0.04068941760916);
1567         hpd.push_back(-0.06453888262870);
1568         hpd.push_back(-0.00000000000000);
1569         hpd.push_back(0.00000000000000);
1570         flength = 10;
1571         break;
1572
1573     case Biorthogonal_5_5:
1574         lpd.push_back(0.00000000000000);
1575         lpd.push_back(0.00000000000000);
1576         lpd.push_back(0.03968708834741);
1577         lpd.push_back(0.00794810863724);
1578         lpd.push_back(-0.05446378846824);
1579         lpd.push_back(0.34560528195603);
1580         lpd.push_back(0.73666018142821);
1581         lpd.push_back(0.34560528195603);
1582         lpd.push_back(-0.05446378846824);
1583         lpd.push_back(0.00794810863724);
1584         lpd.push_back(0.03968708834741);
1585         lpd.push_back(0.00000000000000);
1586         hpd.push_back(-0.01345670945912);
1587         hpd.push_back(-0.00269496688011);
1588         hpd.push_back(0.13670658466433);
1589         hpd.push_back(-0.09350469740094);
1590         hpd.push_back(-0.47680326579848);
1591         hpd.push_back(0.89950610974865);
1592         hpd.push_back(-0.47680326579848);
1593         hpd.push_back(-0.09350469740094);
1594         hpd.push_back(0.13670658466433);
1595         hpd.push_back(-0.00269496688011);
1596         hpd.push_back(-0.01345670945912);
1597         hpd.push_back(0.00000000000000);
1598         flength = 12;
1599         break;
1600
1601     case Biorthogonal_6_8:
1602         lpd.push_back(0.00000000000000);
1603         lpd.push_back(0.00190883173648);
1604         lpd.push_back(-0.00191428612909);
1605         lpd.push_back(-0.01699063986760);
1606         lpd.push_back(0.01193456527973);
1607         lpd.push_back(0.04973290349094);
1608         lpd.push_back(-0.07726317316720);
1609         lpd.push_back(-0.09405920349574);
1610         lpd.push_back(0.42079628460983);
1611         lpd.push_back(0.82592299745840);
1612         lpd.push_back(0.42079628460983);
1613         lpd.push_back(-0.09405920349574);
1614         lpd.push_back(-0.07726317316720);
1615         lpd.push_back(0.04973290349094);
1616         lpd.push_back(0.01193456527973);
1617         lpd.push_back(-0.01699063986760);
1618         lpd.push_back(-0.00191428612909);
1619         lpd.push_back(0.00190883173648);
1620         hpd.push_back(0.00000000000000);
1621         hpd.push_back(-0.00000000000000);
1622         hpd.push_back(0.00000000000000);
1623         hpd.push_back(0.01442628250562);
1624         hpd.push_back(-0.01446750489679);
1625         hpd.push_back(-0.07872200106263);
1626         hpd.push_back(0.04036797903034);
1627         hpd.push_back(0.41784910915027);
1628         hpd.push_back(-0.75890772945365);
1629         hpd.push_back(0.41784910915027);
1630         hpd.push_back(0.04036797903034);
1631         hpd.push_back(-0.07872200106263);
1632         hpd.push_back(-0.01446750489679);
1633         hpd.push_back(0.01442628250562);
1634         hpd.push_back(0.00000000000000);
1635         hpd.push_back(-0.00000000000000);
1636         hpd.push_back(0.00000000000000);
1637         hpd.push_back(-0.00000000000000);
1638         flength = 18;
1639         break;
1640
1641     case Meyer:
1642         lpd.push_back(0.00000000000000);
1643         lpd.push_back(-0.00000150974086);
1644         lpd.push_back(0.00000127876676);
1645         lpd.push_back(0.00000044958556);
1646         lpd.push_back(-0.00000209656887);
1647         lpd.push_back(0.00000172322355);
1648         lpd.push_back(0.00000069808228);
1649         lpd.push_back(-0.00000287940803);
1650         lpd.push_back(0.00000238314839);
1651         lpd.push_back(0.00000098251560);
1652         lpd.push_back(-0.00000421778919);
1653         lpd.push_back(0.00000335350154);
1654         lpd.push_back(0.00000167472186);
1655         lpd.push_back(-0.00000603450134);
1656         lpd.push_back(0.00000483755580);
1657         lpd.push_back(0.00000240228802);
1658         lpd.push_back(-0.00000955630985);
1659         lpd.push_back(0.00000721652769);
1660         lpd.push_back(0.00000484907830);
1661         lpd.push_back(-0.00001420692858);
1662         lpd.push_back(0.00001050391427);
1663         lpd.push_back(0.00000618758030);
1664         lpd.push_back(-0.00002443800585);
1665         lpd.push_back(0.00002010638769);
1666         lpd.push_back(0.00001499352360);
1667         lpd.push_back(-0.00004642876428);
1668         lpd.push_back(0.00003234131191);
1669         lpd.push_back(0.00003740966576);
1670         lpd.push_back(-0.00010277900508);
1671         lpd.push_back(0.00002446195684);
1672         lpd.push_back(0.00014971351539);
1673         lpd.push_back(-0.00007559287026);
1674         lpd.push_back(-0.00013991314822);
1675         lpd.push_back(-0.00009351289388);
1676         lpd.push_back(0.00016118981973);
1677         lpd.push_back(0.00085950021376);
1678         lpd.push_back(-0.00057818579527);
1679         lpd.push_back(-0.00270216873394);
1680         lpd.push_back(0.00219477533646);
1681         lpd.push_back(0.00604551059646);
1682         lpd.push_back(-0.00638672861855);
1683         lpd.push_back(-0.01104464190054);
1684         lpd.push_back(0.01525091315859);
1685         lpd.push_back(0.01740388821018);
1686         lpd.push_back(-0.03209406335451);
1687         lpd.push_back(-0.02432178395952);
1688         lpd.push_back(0.06366730088447);
1689         lpd.push_back(0.03062124394342);
1690         lpd.push_back(-0.13269661535886);
1691         lpd.push_back(-0.03504828739060);
1692         lpd.push_back(0.44409503076653);
1693         lpd.push_back(0.74375100490379);
1694         lpd.push_back(0.44409503076653);
1695         lpd.push_back(-0.03504828739060);
1696         lpd.push_back(-0.13269661535886);
1697         lpd.push_back(0.03062124394342);
1698         lpd.push_back(0.06366730088447);
1699         lpd.push_back(-0.02432178395952);
1700         lpd.push_back(-0.03209406335451);
1701         lpd.push_back(0.01740388821018);
1702         lpd.push_back(0.01525091315859);
1703         lpd.push_back(-0.01104464190054);
1704         lpd.push_back(-0.00638672861855);
1705         lpd.push_back(0.00604551059646);
1706         lpd.push_back(0.00219477533646);
1707         lpd.push_back(-0.00270216873394);
1708         lpd.push_back(-0.00057818579527);
1709         lpd.push_back(0.00085950021376);
1710         lpd.push_back(0.00016118981973);
1711         lpd.push_back(-0.00009351289388);
1712         lpd.push_back(-0.00013991314822);
1713         lpd.push_back(-0.00007559287026);
1714         lpd.push_back(0.00014971351539);
1715         lpd.push_back(0.00002446195684);
1716         lpd.push_back(-0.00010277900508);
1717         lpd.push_back(0.00003740966576);
1718         lpd.push_back(0.00003234131191);
1719         lpd.push_back(-0.00004642876428);
1720         lpd.push_back(0.00001499352360);
1721         lpd.push_back(0.00002010638769);
1722         lpd.push_back(-0.00002443800585);
1723         lpd.push_back(0.00000618758030);
1724         lpd.push_back(0.00001050391427);
1725         lpd.push_back(-0.00001420692858);
1726         lpd.push_back(0.00000484907830);
1727         lpd.push_back(0.00000721652769);
1728         lpd.push_back(-0.00000955630985);
1729         lpd.push_back(0.00000240228802);
1730         lpd.push_back(0.00000483755580);
1731         lpd.push_back(-0.00000603450134);
1732         lpd.push_back(0.00000167472186);
1733         lpd.push_back(0.00000335350154);
1734         lpd.push_back(-0.00000421778919);
1735         lpd.push_back(0.00000098251560);
1736         lpd.push_back(0.00000238314839);
1737         lpd.push_back(-0.00000287940803);
1738         lpd.push_back(0.00000069808228);
1739         lpd.push_back(0.00000172322355);
1740         lpd.push_back(-0.00000209656887);
1741         lpd.push_back(0.00000044958556);
1742         lpd.push_back(0.00000127876676);
1743         lpd.push_back(-0.00000150974086);
1744         hpd.push_back(0.00000150974086);
1745         hpd.push_back(0.00000127876676);
1746         hpd.push_back(-0.00000044958556);
1747         hpd.push_back(-0.00000209656887);
1748         hpd.push_back(-0.00000172322355);
1749         hpd.push_back(0.00000069808228);
1750         hpd.push_back(0.00000287940803);
1751         hpd.push_back(0.00000238314839);
1752         hpd.push_back(-0.00000098251560);
1753         hpd.push_back(-0.00000421778919);
1754         hpd.push_back(-0.00000335350154);
1755         hpd.push_back(0.00000167472186);
1756         hpd.push_back(0.00000603450134);
1757         hpd.push_back(0.00000483755580);
1758         hpd.push_back(-0.00000240228802);
1759         hpd.push_back(-0.00000955630985);
1760         hpd.push_back(-0.00000721652769);
1761         hpd.push_back(0.00000484907830);
1762         hpd.push_back(0.00001420692858);
1763         hpd.push_back(0.00001050391427);
1764         hpd.push_back(-0.00000618758030);
1765         hpd.push_back(-0.00002443800585);
1766         hpd.push_back(-0.00002010638769);
1767         hpd.push_back(0.00001499352360);
1768         hpd.push_back(0.00004642876428);
1769         hpd.push_back(0.00003234131191);
1770         hpd.push_back(-0.00003740966576);
1771         hpd.push_back(-0.00010277900508);
1772         hpd.push_back(-0.00002446195684);
1773         hpd.push_back(0.00014971351539);
1774         hpd.push_back(0.00007559287026);
1775         hpd.push_back(-0.00013991314822);
1776         hpd.push_back(0.00009351289388);
1777         hpd.push_back(0.00016118981973);
1778         hpd.push_back(-0.00085950021376);
1779         hpd.push_back(-0.00057818579527);
1780         hpd.push_back(0.00270216873394);
1781         hpd.push_back(0.00219477533646);
1782         hpd.push_back(-0.00604551059646);
1783         hpd.push_back(-0.00638672861855);
1784         hpd.push_back(0.01104464190054);
1785         hpd.push_back(0.01525091315859);
1786         hpd.push_back(-0.01740388821018);
1787         hpd.push_back(-0.03209406335451);
1788         hpd.push_back(0.02432178395952);
1789         hpd.push_back(0.06366730088447);
1790         hpd.push_back(-0.03062124394342);
1791         hpd.push_back(-0.13269661535886);
1792         hpd.push_back(0.03504828739060);
1793         hpd.push_back(0.44409503076653);
1794         hpd.push_back(-0.74375100490379);
1795         hpd.push_back(0.44409503076653);
1796         hpd.push_back(0.03504828739060);
1797         hpd.push_back(-0.13269661535886);
1798         hpd.push_back(-0.03062124394342);
1799         hpd.push_back(0.06366730088447);
1800         hpd.push_back(0.02432178395952);
1801         hpd.push_back(-0.03209406335451);
1802         hpd.push_back(-0.01740388821018);
1803         hpd.push_back(0.01525091315859);
1804         hpd.push_back(0.01104464190054);
1805         hpd.push_back(-0.00638672861855);
1806         hpd.push_back(-0.00604551059646);
1807         hpd.push_back(0.00219477533646);
1808         hpd.push_back(0.00270216873394);
1809         hpd.push_back(-0.00057818579527);
1810         hpd.push_back(-0.00085950021376);
1811         hpd.push_back(0.00016118981973);
1812         hpd.push_back(0.00009351289388);
1813         hpd.push_back(-0.00013991314822);
1814         hpd.push_back(0.00007559287026);
1815         hpd.push_back(0.00014971351539);
1816         hpd.push_back(-0.00002446195684);
1817         hpd.push_back(-0.00010277900508);
1818         hpd.push_back(-0.00003740966576);
1819         hpd.push_back(0.00003234131191);
1820         hpd.push_back(0.00004642876428);
1821         hpd.push_back(0.00001499352360);
1822         hpd.push_back(-0.00002010638769);
1823         hpd.push_back(-0.00002443800585);
1824         hpd.push_back(-0.00000618758030);
1825         hpd.push_back(0.00001050391427);
1826         hpd.push_back(0.00001420692858);
1827         hpd.push_back(0.00000484907830);
1828         hpd.push_back(-0.00000721652769);
1829         hpd.push_back(-0.00000955630985);
1830         hpd.push_back(-0.00000240228802);
1831         hpd.push_back(0.00000483755580);
1832         hpd.push_back(0.00000603450134);
1833         hpd.push_back(0.00000167472186);
1834         hpd.push_back(-0.00000335350154);
1835         hpd.push_back(-0.00000421778919);
1836         hpd.push_back(-0.00000098251560);
1837         hpd.push_back(0.00000238314839);
1838         hpd.push_back(0.00000287940803);
1839         hpd.push_back(0.00000069808228);
1840         hpd.push_back(-0.00000172322355);
1841         hpd.push_back(-0.00000209656887);
1842         hpd.push_back(-0.00000044958556);
1843         hpd.push_back(0.00000127876676);
1844         hpd.push_back(0.00000150974086);
1845         hpd.push_back(0.00000000000000);
1846         flength = 102;
1847         break;
1848     }
1849
1850     assert(flength == lpd.size());
1851     assert(flength == hpd.size());
1852 }
1853