src/parsec/disk-image/parsec/parsec-benchmark/pkgs/libs/gsl/src/wavelet/daubechies.c - public/gem5-resources - Git at Google

 /* wavelet/daubechies.c
  *
  * Copyright (C) 2004 Ivo Alxneit
  *
  * This program is free software; you can redistribute it and/or modify
  * it under the terms of the GNU General Public License as published by
  * the Free Software Foundation; either version 2 of the License, or (at
  * your option) any later version.
  *
  * This program is distributed in the hope that it will be useful, but
  * WITHOUT ANY WARRANTY; without even the implied warranty of
  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  * General Public License for more details.
  *
  * You should have received a copy of the GNU General Public License
  * along with this program; if not, write to the Free Software
  * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301, USA.
  */

 /*
  * Coefficients for Daubechies wavelets of extremal phase are from
  * I. Daubechies, "Orthonormal Bases of Compactly Supported Wavelets",
  * Communications on Pure and Applied Mathematics, 41 (1988) 909--996
  * (table 1).
  * Additional digits have been obtained using the Mathematica package
  * Daubechies.m by Tong Chen & Meng Xu available at
  * http://www.cwp.mines.edu/wavelets/.
  */

 #include <config.h>
 #include <gsl/gsl_errno.h>
 #include <gsl/gsl_wavelet.h>

 static const double h_4[4] = { 0.48296291314453414337487159986,
   0.83651630373780790557529378092,
   0.22414386804201338102597276224,
   -0.12940952255126038117444941881
 };

 static const double g_4[4] = { -0.12940952255126038117444941881,
   -0.22414386804201338102597276224,
   0.83651630373780790557529378092,
   -0.48296291314453414337487159986
 };

 static const double h_6[6] = { 0.33267055295008261599851158914,
   0.80689150931109257649449360409,
   0.45987750211849157009515194215,
   -0.13501102001025458869638990670,
   -0.08544127388202666169281916918,
   0.03522629188570953660274066472
 };

 static const double g_6[6] = { 0.03522629188570953660274066472,
   0.08544127388202666169281916918,
   -0.13501102001025458869638990670,
   -0.45987750211849157009515194215,
   0.80689150931109257649449360409,
   -0.33267055295008261599851158914
 };

 static const double h_8[8] = { 0.23037781330889650086329118304,
   0.71484657055291564708992195527,
   0.63088076792985890788171633830,
   -0.02798376941685985421141374718,
   -0.18703481171909308407957067279,
   0.03084138183556076362721936253,
   0.03288301166688519973540751355,
   -0.01059740178506903210488320852
 };

 static const double g_8[8] = { -0.01059740178506903210488320852,
   -0.03288301166688519973540751355,
   0.03084138183556076362721936253,
   0.18703481171909308407957067279,
   -0.02798376941685985421141374718,
   -0.63088076792985890788171633830,
   0.71484657055291564708992195527,
   -0.23037781330889650086329118304
 };

 static const double h_10[10] = { 0.16010239797419291448072374802,
   0.60382926979718967054011930653,
   0.72430852843777292772807124410,
   0.13842814590132073150539714634,
   -0.24229488706638203186257137947,
   -0.03224486958463837464847975506,
   0.07757149384004571352313048939,
   -0.00624149021279827427419051911,
   -0.01258075199908199946850973993,
   0.00333572528547377127799818342
 };

 static const double g_10[10] = { 0.00333572528547377127799818342,
   0.01258075199908199946850973993,
   -0.00624149021279827427419051911,
   -0.07757149384004571352313048939,
   -0.03224486958463837464847975506,
   0.24229488706638203186257137947,
   0.13842814590132073150539714634,
   -0.72430852843777292772807124410,
   0.60382926979718967054011930653,
   -0.16010239797419291448072374802
 };

 static const double h_12[12] = { 0.11154074335010946362132391724,
   0.49462389039845308567720417688,
   0.75113390802109535067893449844,
   0.31525035170919762908598965481,
   -0.22626469396543982007631450066,
   -0.12976686756726193556228960588,
   0.09750160558732304910234355254,
   0.02752286553030572862554083950,
   -0.03158203931748602956507908070,
   0.00055384220116149613925191840,
   0.00477725751094551063963597525,
   -0.00107730108530847956485262161
 };

 static const double g_12[12] = { -0.00107730108530847956485262161,
   -0.00477725751094551063963597525,
   0.00055384220116149613925191840,
   0.03158203931748602956507908070,
   0.02752286553030572862554083950,
   -0.09750160558732304910234355254,
   -0.12976686756726193556228960588,
   0.22626469396543982007631450066,
   0.31525035170919762908598965481,
   -0.75113390802109535067893449844,
   0.49462389039845308567720417688,
   -0.11154074335010946362132391724
 };

 static const double h_14[14] = { 0.07785205408500917901996352196,
   0.39653931948191730653900039094,
   0.72913209084623511991694307034,
   0.46978228740519312247159116097,
   -0.14390600392856497540506836221,
   -0.22403618499387498263814042023,
   0.07130921926683026475087657050,
   0.08061260915108307191292248036,
   -0.03802993693501441357959206160,
   -0.01657454163066688065410767489,
   0.01255099855609984061298988603,
   0.00042957797292136652113212912,
   -0.00180164070404749091526826291,
   0.00035371379997452024844629584
 };

 static const double g_14[14] = { 0.00035371379997452024844629584,
   0.00180164070404749091526826291,
   0.00042957797292136652113212912,
   -0.01255099855609984061298988603,
   -0.01657454163066688065410767489,
   0.03802993693501441357959206160,
   0.08061260915108307191292248036,
   -0.07130921926683026475087657050,
   -0.22403618499387498263814042023,
   0.14390600392856497540506836221,
   0.46978228740519312247159116097,
   -0.72913209084623511991694307034,
   0.39653931948191730653900039094,
   -0.07785205408500917901996352196
 };

 static const double h_16[16] = { 0.05441584224310400995500940520,
   0.31287159091429997065916237551,
   0.67563073629728980680780076705,
   0.58535468365420671277126552005,
   -0.01582910525634930566738054788,
   -0.28401554296154692651620313237,
   0.00047248457391328277036059001,
   0.12874742662047845885702928751,
   -0.01736930100180754616961614887,
   -0.04408825393079475150676372324,
   0.01398102791739828164872293057,
   0.00874609404740577671638274325,
   -0.00487035299345157431042218156,
   -0.00039174037337694704629808036,
   0.00067544940645056936636954757,
   -0.00011747678412476953373062823
 };

 static const double g_16[16] = { -0.00011747678412476953373062823,
   -0.00067544940645056936636954757,
   -0.00039174037337694704629808036,
   0.00487035299345157431042218156,
   0.00874609404740577671638274325,
   -0.01398102791739828164872293057,
   -0.04408825393079475150676372324,
   0.01736930100180754616961614887,
   0.12874742662047845885702928751,
   -0.00047248457391328277036059001,
   -0.28401554296154692651620313237,
   0.01582910525634930566738054788,
   0.58535468365420671277126552005,
   -0.67563073629728980680780076705,
   0.31287159091429997065916237551,
   -0.05441584224310400995500940520
 };

 static const double h_18[18] = { 0.03807794736387834658869765888,
   0.24383467461259035373204158165,
   0.60482312369011111190307686743,
   0.65728807805130053807821263905,
   0.13319738582500757619095494590,
   -0.29327378327917490880640319524,
   -0.09684078322297646051350813354,
   0.14854074933810638013507271751,
   0.03072568147933337921231740072,
   -0.06763282906132997367564227483,
   0.00025094711483145195758718975,
   0.02236166212367909720537378270,
   -0.00472320475775139727792570785,
   -0.00428150368246342983449679500,
   0.00184764688305622647661912949,
   0.00023038576352319596720521639,
   -0.00025196318894271013697498868,
   0.00003934732031627159948068988
 };

 static const double g_18[18] = { 0.00003934732031627159948068988,
   0.00025196318894271013697498868,
   0.00023038576352319596720521639,
   -0.00184764688305622647661912949,
   -0.00428150368246342983449679500,
   0.00472320475775139727792570785,
   0.02236166212367909720537378270,
   -0.00025094711483145195758718975,
   -0.06763282906132997367564227483,
   -0.03072568147933337921231740072,
   0.14854074933810638013507271751,
   0.09684078322297646051350813354,
   -0.29327378327917490880640319524,
   -0.13319738582500757619095494590,
   0.65728807805130053807821263905,
   -0.60482312369011111190307686743,
   0.24383467461259035373204158165,
   -0.03807794736387834658869765888
 };

 static const double h_20[20] = { 0.02667005790055555358661744877,
   0.18817680007769148902089297368,
   0.52720118893172558648174482796,
   0.68845903945360356574187178255,
   0.28117234366057746074872699845,
   -0.24984642432731537941610189792,
   -0.19594627437737704350429925432,
   0.12736934033579326008267723320,
   0.09305736460357235116035228984,
   -0.07139414716639708714533609308,
   -0.02945753682187581285828323760,
   0.03321267405934100173976365318,
   0.00360655356695616965542329142,
   -0.01073317548333057504431811411,
   0.00139535174705290116578931845,
   0.00199240529518505611715874224,
   -0.00068585669495971162656137098,
   -0.00011646685512928545095148097,
   0.00009358867032006959133405013,
   -0.00001326420289452124481243668
 };

 static const double g_20[20] = { -0.00001326420289452124481243668,
   -0.00009358867032006959133405013,
   -0.00011646685512928545095148097,
   0.00068585669495971162656137098,
   0.00199240529518505611715874224,
   -0.00139535174705290116578931845,
   -0.01073317548333057504431811411,
   -0.00360655356695616965542329142,
   0.03321267405934100173976365318,
   0.02945753682187581285828323760,
   -0.07139414716639708714533609308,
   -0.09305736460357235116035228984,
   0.12736934033579326008267723320,
   0.19594627437737704350429925432,
   -0.24984642432731537941610189792,
   -0.28117234366057746074872699845,
   0.68845903945360356574187178255,
   -0.52720118893172558648174482796,
   0.18817680007769148902089297368,
   -0.02667005790055555358661744877
 };

 static int
 daubechies_init (const double **h1, const double **g1, const double **h2,
                  const double **g2, size_t * nc, size_t * offset,
                  size_t member)
 {
   switch (member)
     {
     case 4:
       *h1 = h_4;
       *g1 = g_4;
       *h2 = h_4;
       *g2 = g_4;
       break;

     case 6:
       *h1 = h_6;
       *g1 = g_6;
       *h2 = h_6;
       *g2 = g_6;
       break;

     case 8:
       *h1 = h_8;
       *g1 = g_8;
       *h2 = h_8;
       *g2 = g_8;
       break;

     case 10:
       *h1 = h_10;
       *g1 = g_10;
       *h2 = h_10;
       *g2 = g_10;
       break;

     case 12:
       *h1 = h_12;
       *g1 = g_12;
       *h2 = h_12;
       *g2 = g_12;
       break;

     case 14:
       *h1 = h_14;
       *g1 = g_14;
       *h2 = h_14;
       *g2 = g_14;
       break;

     case 16:
       *h1 = h_16;
       *g1 = g_16;
       *h2 = h_16;
       *g2 = g_16;
       break;

     case 18:
       *h1 = h_18;
       *g1 = g_18;
       *h2 = h_18;
       *g2 = g_18;
       break;

     case 20:
       *h1 = h_20;
       *g1 = g_20;
       *h2 = h_20;
       *g2 = g_20;
       break;

     default:
       return GSL_FAILURE;
     }

   *nc = member;
   *offset = 0;

   return GSL_SUCCESS;
 }

 static int
 daubechies_centered_init (const double **h1, const double **g1,
                           const double **h2, const double **g2, size_t * nc,
                           size_t * offset, size_t member)
 {
   switch (member)
     {
     case 4:
       *h1 = h_4;
       *g1 = g_4;
       *h2 = h_4;
       *g2 = g_4;
       break;

     case 6:
       *h1 = h_6;
       *g1 = g_6;
       *h2 = h_6;
       *g2 = g_6;
       break;

     case 8:
       *h1 = h_8;
       *g1 = g_8;
       *h2 = h_8;
       *g2 = g_8;
       break;

     case 10:
       *h1 = h_10;
       *g1 = g_10;
       *h2 = h_10;
       *g2 = g_10;
       break;

     case 12:
       *h1 = h_12;
       *g1 = g_12;
       *h2 = h_12;
       *g2 = g_12;
       break;

     case 14:
       *h1 = h_14;
       *g1 = g_14;
       *h2 = h_14;
       *g2 = g_14;
       break;

     case 16:
       *h1 = h_16;
       *g1 = g_16;
       *h2 = h_16;
       *g2 = g_16;
       break;

     case 18:
       *h1 = h_18;
       *g1 = g_18;
       *h2 = h_18;
       *g2 = g_18;
       break;

     case 20:
       *h1 = h_20;
       *g1 = g_20;
       *h2 = h_20;
       *g2 = g_20;
       break;

     default:
       return GSL_FAILURE;
     }

   *nc = member;
   *offset = (member >> 1);

   return GSL_SUCCESS;
 }

 static const gsl_wavelet_type daubechies_type = {
   "daubechies",
   &daubechies_init
 };

 static const gsl_wavelet_type daubechies_centered_type = {
   "daubechies-centered",
   &daubechies_centered_init
 };

 const gsl_wavelet_type *gsl_wavelet_daubechies = &daubechies_type;
 const gsl_wavelet_type *gsl_wavelet_daubechies_centered =
   &daubechies_centered_type;
	/* wavelet/daubechies.c
	*
	* Copyright (C) 2004 Ivo Alxneit
	*
	* This program is free software; you can redistribute it and/or modify
	* it under the terms of the GNU General Public License as published by
	* the Free Software Foundation; either version 2 of the License, or (at
	* your option) any later version.
	*
	* This program is distributed in the hope that it will be useful, but
	* WITHOUT ANY WARRANTY; without even the implied warranty of
	* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
	* General Public License for more details.
	*
	* You should have received a copy of the GNU General Public License
	* along with this program; if not, write to the Free Software
	* Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301, USA.
	*/

	/*
	* Coefficients for Daubechies wavelets of extremal phase are from
	* I. Daubechies, "Orthonormal Bases of Compactly Supported Wavelets",
	* Communications on Pure and Applied Mathematics, 41 (1988) 909--996
	* (table 1).
	* Additional digits have been obtained using the Mathematica package
	* Daubechies.m by Tong Chen & Meng Xu available at
	* http://www.cwp.mines.edu/wavelets/.
	*/

	#include <config.h>
	#include <gsl/gsl_errno.h>
	#include <gsl/gsl_wavelet.h>

	static const double h_4[4] = { 0.48296291314453414337487159986,
	0.83651630373780790557529378092,
	0.22414386804201338102597276224,
	-0.12940952255126038117444941881
	};

	static const double g_4[4] = { -0.12940952255126038117444941881,
	-0.22414386804201338102597276224,
	0.83651630373780790557529378092,
	-0.48296291314453414337487159986
	};

	static const double h_6[6] = { 0.33267055295008261599851158914,
	0.80689150931109257649449360409,
	0.45987750211849157009515194215,
	-0.13501102001025458869638990670,
	-0.08544127388202666169281916918,
	0.03522629188570953660274066472
	};

	static const double g_6[6] = { 0.03522629188570953660274066472,
	0.08544127388202666169281916918,
	-0.13501102001025458869638990670,
	-0.45987750211849157009515194215,
	0.80689150931109257649449360409,
	-0.33267055295008261599851158914
	};

	static const double h_8[8] = { 0.23037781330889650086329118304,
	0.71484657055291564708992195527,
	0.63088076792985890788171633830,
	-0.02798376941685985421141374718,
	-0.18703481171909308407957067279,
	0.03084138183556076362721936253,
	0.03288301166688519973540751355,
	-0.01059740178506903210488320852
	};

	static const double g_8[8] = { -0.01059740178506903210488320852,
	-0.03288301166688519973540751355,
	0.03084138183556076362721936253,
	0.18703481171909308407957067279,
	-0.02798376941685985421141374718,
	-0.63088076792985890788171633830,
	0.71484657055291564708992195527,
	-0.23037781330889650086329118304
	};

	static const double h_10[10] = { 0.16010239797419291448072374802,
	0.60382926979718967054011930653,
	0.72430852843777292772807124410,
	0.13842814590132073150539714634,
	-0.24229488706638203186257137947,
	-0.03224486958463837464847975506,
	0.07757149384004571352313048939,
	-0.00624149021279827427419051911,
	-0.01258075199908199946850973993,
	0.00333572528547377127799818342
	};

	static const double g_10[10] = { 0.00333572528547377127799818342,
	0.01258075199908199946850973993,
	-0.00624149021279827427419051911,
	-0.07757149384004571352313048939,
	-0.03224486958463837464847975506,
	0.24229488706638203186257137947,
	0.13842814590132073150539714634,
	-0.72430852843777292772807124410,
	0.60382926979718967054011930653,
	-0.16010239797419291448072374802
	};

	static const double h_12[12] = { 0.11154074335010946362132391724,
	0.49462389039845308567720417688,
	0.75113390802109535067893449844,
	0.31525035170919762908598965481,
	-0.22626469396543982007631450066,
	-0.12976686756726193556228960588,
	0.09750160558732304910234355254,
	0.02752286553030572862554083950,
	-0.03158203931748602956507908070,
	0.00055384220116149613925191840,
	0.00477725751094551063963597525,
	-0.00107730108530847956485262161
	};

	static const double g_12[12] = { -0.00107730108530847956485262161,
	-0.00477725751094551063963597525,
	0.00055384220116149613925191840,
	0.03158203931748602956507908070,
	0.02752286553030572862554083950,
	-0.09750160558732304910234355254,
	-0.12976686756726193556228960588,
	0.22626469396543982007631450066,
	0.31525035170919762908598965481,
	-0.75113390802109535067893449844,
	0.49462389039845308567720417688,
	-0.11154074335010946362132391724
	};

	static const double h_14[14] = { 0.07785205408500917901996352196,
	0.39653931948191730653900039094,
	0.72913209084623511991694307034,
	0.46978228740519312247159116097,
	-0.14390600392856497540506836221,
	-0.22403618499387498263814042023,
	0.07130921926683026475087657050,
	0.08061260915108307191292248036,
	-0.03802993693501441357959206160,
	-0.01657454163066688065410767489,
	0.01255099855609984061298988603,
	0.00042957797292136652113212912,
	-0.00180164070404749091526826291,
	0.00035371379997452024844629584
	};

	static const double g_14[14] = { 0.00035371379997452024844629584,
	0.00180164070404749091526826291,
	0.00042957797292136652113212912,
	-0.01255099855609984061298988603,
	-0.01657454163066688065410767489,
	0.03802993693501441357959206160,
	0.08061260915108307191292248036,
	-0.07130921926683026475087657050,
	-0.22403618499387498263814042023,
	0.14390600392856497540506836221,
	0.46978228740519312247159116097,
	-0.72913209084623511991694307034,
	0.39653931948191730653900039094,
	-0.07785205408500917901996352196
	};

	static const double h_16[16] = { 0.05441584224310400995500940520,
	0.31287159091429997065916237551,
	0.67563073629728980680780076705,
	0.58535468365420671277126552005,
	-0.01582910525634930566738054788,
	-0.28401554296154692651620313237,
	0.00047248457391328277036059001,
	0.12874742662047845885702928751,
	-0.01736930100180754616961614887,
	-0.04408825393079475150676372324,
	0.01398102791739828164872293057,
	0.00874609404740577671638274325,
	-0.00487035299345157431042218156,
	-0.00039174037337694704629808036,
	0.00067544940645056936636954757,
	-0.00011747678412476953373062823
	};

	static const double g_16[16] = { -0.00011747678412476953373062823,
	-0.00067544940645056936636954757,
	-0.00039174037337694704629808036,
	0.00487035299345157431042218156,
	0.00874609404740577671638274325,
	-0.01398102791739828164872293057,
	-0.04408825393079475150676372324,
	0.01736930100180754616961614887,
	0.12874742662047845885702928751,
	-0.00047248457391328277036059001,
	-0.28401554296154692651620313237,
	0.01582910525634930566738054788,
	0.58535468365420671277126552005,
	-0.67563073629728980680780076705,
	0.31287159091429997065916237551,
	-0.05441584224310400995500940520
	};

	static const double h_18[18] = { 0.03807794736387834658869765888,
	0.24383467461259035373204158165,
	0.60482312369011111190307686743,
	0.65728807805130053807821263905,
	0.13319738582500757619095494590,
	-0.29327378327917490880640319524,
	-0.09684078322297646051350813354,
	0.14854074933810638013507271751,
	0.03072568147933337921231740072,
	-0.06763282906132997367564227483,
	0.00025094711483145195758718975,
	0.02236166212367909720537378270,
	-0.00472320475775139727792570785,
	-0.00428150368246342983449679500,
	0.00184764688305622647661912949,
	0.00023038576352319596720521639,
	-0.00025196318894271013697498868,
	0.00003934732031627159948068988
	};

	static const double g_18[18] = { 0.00003934732031627159948068988,
	0.00025196318894271013697498868,
	0.00023038576352319596720521639,
	-0.00184764688305622647661912949,
	-0.00428150368246342983449679500,
	0.00472320475775139727792570785,
	0.02236166212367909720537378270,
	-0.00025094711483145195758718975,
	-0.06763282906132997367564227483,
	-0.03072568147933337921231740072,
	0.14854074933810638013507271751,
	0.09684078322297646051350813354,
	-0.29327378327917490880640319524,
	-0.13319738582500757619095494590,
	0.65728807805130053807821263905,
	-0.60482312369011111190307686743,
	0.24383467461259035373204158165,
	-0.03807794736387834658869765888
	};

	static const double h_20[20] = { 0.02667005790055555358661744877,
	0.18817680007769148902089297368,
	0.52720118893172558648174482796,
	0.68845903945360356574187178255,
	0.28117234366057746074872699845,
	-0.24984642432731537941610189792,
	-0.19594627437737704350429925432,
	0.12736934033579326008267723320,
	0.09305736460357235116035228984,
	-0.07139414716639708714533609308,
	-0.02945753682187581285828323760,
	0.03321267405934100173976365318,
	0.00360655356695616965542329142,
	-0.01073317548333057504431811411,
	0.00139535174705290116578931845,
	0.00199240529518505611715874224,
	-0.00068585669495971162656137098,
	-0.00011646685512928545095148097,
	0.00009358867032006959133405013,
	-0.00001326420289452124481243668
	};

	static const double g_20[20] = { -0.00001326420289452124481243668,
	-0.00009358867032006959133405013,
	-0.00011646685512928545095148097,
	0.00068585669495971162656137098,
	0.00199240529518505611715874224,
	-0.00139535174705290116578931845,
	-0.01073317548333057504431811411,
	-0.00360655356695616965542329142,
	0.03321267405934100173976365318,
	0.02945753682187581285828323760,
	-0.07139414716639708714533609308,
	-0.09305736460357235116035228984,
	0.12736934033579326008267723320,
	0.19594627437737704350429925432,
	-0.24984642432731537941610189792,
	-0.28117234366057746074872699845,
	0.68845903945360356574187178255,
	-0.52720118893172558648174482796,
	0.18817680007769148902089297368,
	-0.02667005790055555358661744877
	};

	static int
	daubechies_init (const double h1, const double g1, const double **h2,
	const double *g2, size_t nc, size_t * offset,
	size_t member)
	{
	switch (member)
	{
	case 4:
	*h1 = h_4;
	*g1 = g_4;
	*h2 = h_4;
	*g2 = g_4;
	break;

	case 6:
	*h1 = h_6;
	*g1 = g_6;
	*h2 = h_6;
	*g2 = g_6;
	break;

	case 8:
	*h1 = h_8;
	*g1 = g_8;
	*h2 = h_8;
	*g2 = g_8;
	break;

	case 10:
	*h1 = h_10;
	*g1 = g_10;
	*h2 = h_10;
	*g2 = g_10;
	break;

	case 12:
	*h1 = h_12;
	*g1 = g_12;
	*h2 = h_12;
	*g2 = g_12;
	break;

	case 14:
	*h1 = h_14;
	*g1 = g_14;
	*h2 = h_14;
	*g2 = g_14;
	break;

	case 16:
	*h1 = h_16;
	*g1 = g_16;
	*h2 = h_16;
	*g2 = g_16;
	break;

	case 18:
	*h1 = h_18;
	*g1 = g_18;
	*h2 = h_18;
	*g2 = g_18;
	break;

	case 20:
	*h1 = h_20;
	*g1 = g_20;
	*h2 = h_20;
	*g2 = g_20;
	break;

	default:
	return GSL_FAILURE;
	}

	*nc = member;
	*offset = 0;

	return GSL_SUCCESS;
	}

	static int
	daubechies_centered_init (const double h1, const double g1,
	const double h2, const double g2, size_t * nc,
	size_t * offset, size_t member)
	{
	switch (member)
	{
	case 4:
	*h1 = h_4;
	*g1 = g_4;
	*h2 = h_4;
	*g2 = g_4;
	break;

	case 6:
	*h1 = h_6;
	*g1 = g_6;
	*h2 = h_6;
	*g2 = g_6;
	break;

	case 8:
	*h1 = h_8;
	*g1 = g_8;
	*h2 = h_8;
	*g2 = g_8;
	break;

	case 10:
	*h1 = h_10;
	*g1 = g_10;
	*h2 = h_10;
	*g2 = g_10;
	break;

	case 12:
	*h1 = h_12;
	*g1 = g_12;
	*h2 = h_12;
	*g2 = g_12;
	break;

	case 14:
	*h1 = h_14;
	*g1 = g_14;
	*h2 = h_14;
	*g2 = g_14;
	break;

	case 16:
	*h1 = h_16;
	*g1 = g_16;
	*h2 = h_16;
	*g2 = g_16;
	break;

	case 18:
	*h1 = h_18;
	*g1 = g_18;
	*h2 = h_18;
	*g2 = g_18;
	break;

	case 20:
	*h1 = h_20;
	*g1 = g_20;
	*h2 = h_20;
	*g2 = g_20;
	break;

	default:
	return GSL_FAILURE;
	}

	*nc = member;
	*offset = (member >> 1);

	return GSL_SUCCESS;
	}

	static const gsl_wavelet_type daubechies_type = {
	"daubechies",
	&daubechies_init
	};

	static const gsl_wavelet_type daubechies_centered_type = {
	"daubechies-centered",
	&daubechies_centered_init
	};

	const gsl_wavelet_type *gsl_wavelet_daubechies = &daubechies_type;
	const gsl_wavelet_type *gsl_wavelet_daubechies_centered =
	&daubechies_centered_type;