|
Packit |
67cb25 |
/* wavelet/bspline.c
|
|
Packit |
67cb25 |
*
|
|
Packit |
67cb25 |
* Copyright (C) 2004 Ivo Alxneit
|
|
Packit |
67cb25 |
*
|
|
Packit |
67cb25 |
* This program is free software; you can redistribute it and/or modify
|
|
Packit |
67cb25 |
* it under the terms of the GNU General Public License as published by
|
|
Packit |
67cb25 |
* the Free Software Foundation; either version 3 of the License, or (at
|
|
Packit |
67cb25 |
* your option) any later version.
|
|
Packit |
67cb25 |
*
|
|
Packit |
67cb25 |
* This program is distributed in the hope that it will be useful, but
|
|
Packit |
67cb25 |
* WITHOUT ANY WARRANTY; without even the implied warranty of
|
|
Packit |
67cb25 |
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
|
|
Packit |
67cb25 |
* General Public License for more details.
|
|
Packit |
67cb25 |
*
|
|
Packit |
67cb25 |
* You should have received a copy of the GNU General Public License
|
|
Packit |
67cb25 |
* along with this program; if not, write to the Free Software
|
|
Packit |
67cb25 |
* Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301, USA.
|
|
Packit |
67cb25 |
*/
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
/* Coefficients are from A. Cohen, I. Daubechies, and J.-C. Feauveau;
|
|
Packit |
67cb25 |
* "Biorthogonal Bases of Compactly Supported Wavelets", Communications
|
|
Packit |
67cb25 |
* on Pure and Applied Mathematics, 45 (1992) 485--560 (table 6.1).
|
|
Packit |
67cb25 |
*
|
|
Packit |
67cb25 |
* Note the following errors in table 1:
|
|
Packit |
67cb25 |
*
|
|
Packit |
67cb25 |
* N = 2, N~ = 4, m0~
|
|
Packit |
67cb25 |
* the second term in z^-1 (45/64 z^-1) should be left out.
|
|
Packit |
67cb25 |
*
|
|
Packit |
67cb25 |
* N = 3, N~ = 7, m0~
|
|
Packit |
67cb25 |
* the term 336z^-3 should read 363z^-3.
|
|
Packit |
67cb25 |
*/
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
#include <config.h>
|
|
Packit |
67cb25 |
#include <gsl/gsl_errno.h>
|
|
Packit |
67cb25 |
#include <gsl/gsl_math.h>
|
|
Packit |
67cb25 |
#include <gsl/gsl_wavelet.h>
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
static const double h1_103[6] = { -0.0883883476483184405501055452631,
|
|
Packit |
67cb25 |
0.0883883476483184405501055452631,
|
|
Packit |
67cb25 |
M_SQRT1_2,
|
|
Packit |
67cb25 |
M_SQRT1_2,
|
|
Packit |
67cb25 |
0.0883883476483184405501055452631,
|
|
Packit |
67cb25 |
-0.0883883476483184405501055452631
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
static const double g2_103[6] = { -0.0883883476483184405501055452631,
|
|
Packit |
67cb25 |
-0.0883883476483184405501055452631,
|
|
Packit |
67cb25 |
M_SQRT1_2,
|
|
Packit |
67cb25 |
-(M_SQRT1_2),
|
|
Packit |
67cb25 |
0.0883883476483184405501055452631,
|
|
Packit |
67cb25 |
0.0883883476483184405501055452631
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
static const double h1_105[10] = { 0.0165728151840597076031447897368,
|
|
Packit |
67cb25 |
-0.0165728151840597076031447897368,
|
|
Packit |
67cb25 |
-0.1215339780164378557563951247368,
|
|
Packit |
67cb25 |
0.1215339780164378557563951247368,
|
|
Packit |
67cb25 |
M_SQRT1_2,
|
|
Packit |
67cb25 |
M_SQRT1_2,
|
|
Packit |
67cb25 |
0.1215339780164378557563951247368,
|
|
Packit |
67cb25 |
-0.1215339780164378557563951247368,
|
|
Packit |
67cb25 |
-0.0165728151840597076031447897368,
|
|
Packit |
67cb25 |
0.0165728151840597076031447897368
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
static const double g2_105[10] = { 0.0165728151840597076031447897368,
|
|
Packit |
67cb25 |
0.0165728151840597076031447897368,
|
|
Packit |
67cb25 |
-0.1215339780164378557563951247368,
|
|
Packit |
67cb25 |
-0.1215339780164378557563951247368,
|
|
Packit |
67cb25 |
M_SQRT1_2,
|
|
Packit |
67cb25 |
-(M_SQRT1_2),
|
|
Packit |
67cb25 |
0.1215339780164378557563951247368,
|
|
Packit |
67cb25 |
0.1215339780164378557563951247368,
|
|
Packit |
67cb25 |
-0.0165728151840597076031447897368,
|
|
Packit |
67cb25 |
-0.0165728151840597076031447897368
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
static const double g1_1[10] = { 0.0, 0.0, 0.0, 0.0,
|
|
Packit |
67cb25 |
M_SQRT1_2,
|
|
Packit |
67cb25 |
-(M_SQRT1_2),
|
|
Packit |
67cb25 |
0.0, 0.0, 0.0, 0.0
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
static const double h2_1[10] = { 0.0, 0.0, 0.0, 0.0,
|
|
Packit |
67cb25 |
M_SQRT1_2,
|
|
Packit |
67cb25 |
M_SQRT1_2,
|
|
Packit |
67cb25 |
0.0, 0.0, 0.0, 0.0
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
static const double h1_202[6] = { -0.1767766952966368811002110905262,
|
|
Packit |
67cb25 |
0.3535533905932737622004221810524,
|
|
Packit |
67cb25 |
1.0606601717798212866012665431573,
|
|
Packit |
67cb25 |
0.3535533905932737622004221810524,
|
|
Packit |
67cb25 |
-0.1767766952966368811002110905262,
|
|
Packit |
67cb25 |
0.0
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
static const double g2_202[6] = { 0.0,
|
|
Packit |
67cb25 |
-0.1767766952966368811002110905262,
|
|
Packit |
67cb25 |
-0.3535533905932737622004221810524,
|
|
Packit |
67cb25 |
1.0606601717798212866012665431573,
|
|
Packit |
67cb25 |
-0.3535533905932737622004221810524,
|
|
Packit |
67cb25 |
-0.1767766952966368811002110905262
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
static const double h1_204[10] = { 0.0331456303681194152062895794737,
|
|
Packit |
67cb25 |
-0.0662912607362388304125791589473,
|
|
Packit |
67cb25 |
-0.1767766952966368811002110905262,
|
|
Packit |
67cb25 |
0.4198446513295125926130013399998,
|
|
Packit |
67cb25 |
0.9943689110435824561886873842099,
|
|
Packit |
67cb25 |
0.4198446513295125926130013399998,
|
|
Packit |
67cb25 |
-0.1767766952966368811002110905262,
|
|
Packit |
67cb25 |
-0.0662912607362388304125791589473,
|
|
Packit |
67cb25 |
0.0331456303681194152062895794737,
|
|
Packit |
67cb25 |
0.0
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
static const double g2_204[10] = { 0.0,
|
|
Packit |
67cb25 |
0.0331456303681194152062895794737,
|
|
Packit |
67cb25 |
0.0662912607362388304125791589473,
|
|
Packit |
67cb25 |
-0.1767766952966368811002110905262,
|
|
Packit |
67cb25 |
-0.4198446513295125926130013399998,
|
|
Packit |
67cb25 |
0.9943689110435824561886873842099,
|
|
Packit |
67cb25 |
-0.4198446513295125926130013399998,
|
|
Packit |
67cb25 |
-0.1767766952966368811002110905262,
|
|
Packit |
67cb25 |
0.0662912607362388304125791589473,
|
|
Packit |
67cb25 |
0.0331456303681194152062895794737
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
static const double h1_206[14] = { -0.0069053396600248781679769957237,
|
|
Packit |
67cb25 |
0.0138106793200497563359539914474,
|
|
Packit |
67cb25 |
0.0469563096881691715422435709210,
|
|
Packit |
67cb25 |
-0.1077232986963880994204411332894,
|
|
Packit |
67cb25 |
-0.1698713556366120029322340948025,
|
|
Packit |
67cb25 |
0.4474660099696121052849093228945,
|
|
Packit |
67cb25 |
0.9667475524034829435167794013152,
|
|
Packit |
67cb25 |
0.4474660099696121052849093228945,
|
|
Packit |
67cb25 |
-0.1698713556366120029322340948025,
|
|
Packit |
67cb25 |
-0.1077232986963880994204411332894,
|
|
Packit |
67cb25 |
0.0469563096881691715422435709210,
|
|
Packit |
67cb25 |
0.0138106793200497563359539914474,
|
|
Packit |
67cb25 |
-0.0069053396600248781679769957237,
|
|
Packit |
67cb25 |
0.0
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
static const double g2_206[14] = { 0.0,
|
|
Packit |
67cb25 |
-0.0069053396600248781679769957237,
|
|
Packit |
67cb25 |
-0.0138106793200497563359539914474,
|
|
Packit |
67cb25 |
0.0469563096881691715422435709210,
|
|
Packit |
67cb25 |
0.1077232986963880994204411332894,
|
|
Packit |
67cb25 |
-0.1698713556366120029322340948025,
|
|
Packit |
67cb25 |
-0.4474660099696121052849093228945,
|
|
Packit |
67cb25 |
0.9667475524034829435167794013152,
|
|
Packit |
67cb25 |
-0.4474660099696121052849093228945,
|
|
Packit |
67cb25 |
-0.1698713556366120029322340948025,
|
|
Packit |
67cb25 |
0.1077232986963880994204411332894,
|
|
Packit |
67cb25 |
0.0469563096881691715422435709210,
|
|
Packit |
67cb25 |
-0.0138106793200497563359539914474,
|
|
Packit |
67cb25 |
-0.0069053396600248781679769957237,
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
static const double h1_208[18] = { 0.0015105430506304420992449678146,
|
|
Packit |
67cb25 |
-0.0030210861012608841984899356291,
|
|
Packit |
67cb25 |
-0.0129475118625466465649568669819,
|
|
Packit |
67cb25 |
0.0289161098263541773284036695929,
|
|
Packit |
67cb25 |
0.0529984818906909399392234421792,
|
|
Packit |
67cb25 |
-0.1349130736077360572068505539514,
|
|
Packit |
67cb25 |
-0.1638291834340902345352542235443,
|
|
Packit |
67cb25 |
0.4625714404759165262773590010400,
|
|
Packit |
67cb25 |
0.9516421218971785225243297231697,
|
|
Packit |
67cb25 |
0.4625714404759165262773590010400,
|
|
Packit |
67cb25 |
-0.1638291834340902345352542235443,
|
|
Packit |
67cb25 |
-0.1349130736077360572068505539514,
|
|
Packit |
67cb25 |
0.0529984818906909399392234421792,
|
|
Packit |
67cb25 |
0.0289161098263541773284036695929,
|
|
Packit |
67cb25 |
-0.0129475118625466465649568669819,
|
|
Packit |
67cb25 |
-0.0030210861012608841984899356291,
|
|
Packit |
67cb25 |
0.0015105430506304420992449678146,
|
|
Packit |
67cb25 |
0.0
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
static const double g2_208[18] = { 0.0,
|
|
Packit |
67cb25 |
0.0015105430506304420992449678146,
|
|
Packit |
67cb25 |
0.0030210861012608841984899356291,
|
|
Packit |
67cb25 |
-0.0129475118625466465649568669819,
|
|
Packit |
67cb25 |
-0.0289161098263541773284036695929,
|
|
Packit |
67cb25 |
0.0529984818906909399392234421792,
|
|
Packit |
67cb25 |
0.1349130736077360572068505539514,
|
|
Packit |
67cb25 |
-0.1638291834340902345352542235443,
|
|
Packit |
67cb25 |
-0.4625714404759165262773590010400,
|
|
Packit |
67cb25 |
0.9516421218971785225243297231697,
|
|
Packit |
67cb25 |
-0.4625714404759165262773590010400,
|
|
Packit |
67cb25 |
-0.1638291834340902345352542235443,
|
|
Packit |
67cb25 |
0.1349130736077360572068505539514,
|
|
Packit |
67cb25 |
0.0529984818906909399392234421792,
|
|
Packit |
67cb25 |
-0.0289161098263541773284036695929,
|
|
Packit |
67cb25 |
-0.0129475118625466465649568669819,
|
|
Packit |
67cb25 |
0.0030210861012608841984899356291,
|
|
Packit |
67cb25 |
0.0015105430506304420992449678146,
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
static const double h2_2[18] = { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
|
|
Packit |
67cb25 |
0.3535533905932737622004221810524,
|
|
Packit |
67cb25 |
0.7071067811865475244008443621048,
|
|
Packit |
67cb25 |
0.3535533905932737622004221810524,
|
|
Packit |
67cb25 |
0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
static const double g1_2[18] = { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
|
|
Packit |
67cb25 |
-0.3535533905932737622004221810524,
|
|
Packit |
67cb25 |
0.7071067811865475244008443621048,
|
|
Packit |
67cb25 |
-0.3535533905932737622004221810524,
|
|
Packit |
67cb25 |
0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
static const double h1_301[4] = { -0.3535533905932737622004221810524,
|
|
Packit |
67cb25 |
1.0606601717798212866012665431573,
|
|
Packit |
67cb25 |
1.0606601717798212866012665431573,
|
|
Packit |
67cb25 |
-0.3535533905932737622004221810524
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
static const double g2_301[4] = { 0.3535533905932737622004221810524,
|
|
Packit |
67cb25 |
1.0606601717798212866012665431573,
|
|
Packit |
67cb25 |
-1.0606601717798212866012665431573,
|
|
Packit |
67cb25 |
-0.3535533905932737622004221810524
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
static const double h1_303[8] = { 0.0662912607362388304125791589473,
|
|
Packit |
67cb25 |
-0.1988737822087164912377374768420,
|
|
Packit |
67cb25 |
-0.1546796083845572709626847042104,
|
|
Packit |
67cb25 |
0.9943689110435824561886873842099,
|
|
Packit |
67cb25 |
0.9943689110435824561886873842099,
|
|
Packit |
67cb25 |
-0.1546796083845572709626847042104,
|
|
Packit |
67cb25 |
-0.1988737822087164912377374768420,
|
|
Packit |
67cb25 |
0.0662912607362388304125791589473
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
static const double g2_303[8] = { -0.0662912607362388304125791589473,
|
|
Packit |
67cb25 |
-0.1988737822087164912377374768420,
|
|
Packit |
67cb25 |
0.1546796083845572709626847042104,
|
|
Packit |
67cb25 |
0.9943689110435824561886873842099,
|
|
Packit |
67cb25 |
-0.9943689110435824561886873842099,
|
|
Packit |
67cb25 |
-0.1546796083845572709626847042104,
|
|
Packit |
67cb25 |
0.1988737822087164912377374768420,
|
|
Packit |
67cb25 |
0.0662912607362388304125791589473
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
static const double h1_305[12] = { -0.0138106793200497563359539914474,
|
|
Packit |
67cb25 |
0.0414320379601492690078619743421,
|
|
Packit |
67cb25 |
0.0524805814161890740766251675000,
|
|
Packit |
67cb25 |
-0.2679271788089652729175074340788,
|
|
Packit |
67cb25 |
-0.0718155324642587329469607555263,
|
|
Packit |
67cb25 |
0.9667475524034829435167794013152,
|
|
Packit |
67cb25 |
0.9667475524034829435167794013152,
|
|
Packit |
67cb25 |
-0.0718155324642587329469607555263,
|
|
Packit |
67cb25 |
-0.2679271788089652729175074340788,
|
|
Packit |
67cb25 |
0.0524805814161890740766251675000,
|
|
Packit |
67cb25 |
0.0414320379601492690078619743421,
|
|
Packit |
67cb25 |
-0.0138106793200497563359539914474
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
static const double g2_305[12] = { 0.0138106793200497563359539914474,
|
|
Packit |
67cb25 |
0.0414320379601492690078619743421,
|
|
Packit |
67cb25 |
-0.0524805814161890740766251675000,
|
|
Packit |
67cb25 |
-0.2679271788089652729175074340788,
|
|
Packit |
67cb25 |
0.0718155324642587329469607555263,
|
|
Packit |
67cb25 |
0.9667475524034829435167794013152,
|
|
Packit |
67cb25 |
-0.9667475524034829435167794013152,
|
|
Packit |
67cb25 |
-0.0718155324642587329469607555263,
|
|
Packit |
67cb25 |
0.2679271788089652729175074340788,
|
|
Packit |
67cb25 |
0.0524805814161890740766251675000,
|
|
Packit |
67cb25 |
-0.0414320379601492690078619743421,
|
|
Packit |
67cb25 |
-0.0138106793200497563359539914474
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
static const double h1_307[16] = { 0.0030210861012608841984899356291,
|
|
Packit |
67cb25 |
-0.0090632583037826525954698068873,
|
|
Packit |
67cb25 |
-0.0168317654213106405344439270765,
|
|
Packit |
67cb25 |
0.0746639850740189951912512662623,
|
|
Packit |
67cb25 |
0.0313329787073628846871956180962,
|
|
Packit |
67cb25 |
-0.3011591259228349991008967259990,
|
|
Packit |
67cb25 |
-0.0264992409453454699696117210896,
|
|
Packit |
67cb25 |
0.9516421218971785225243297231697,
|
|
Packit |
67cb25 |
0.9516421218971785225243297231697,
|
|
Packit |
67cb25 |
-0.0264992409453454699696117210896,
|
|
Packit |
67cb25 |
-0.3011591259228349991008967259990,
|
|
Packit |
67cb25 |
0.0313329787073628846871956180962,
|
|
Packit |
67cb25 |
0.0746639850740189951912512662623,
|
|
Packit |
67cb25 |
-0.0168317654213106405344439270765,
|
|
Packit |
67cb25 |
-0.0090632583037826525954698068873,
|
|
Packit |
67cb25 |
0.0030210861012608841984899356291
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
static const double g2_307[16] = { -0.0030210861012608841984899356291,
|
|
Packit |
67cb25 |
-0.0090632583037826525954698068873,
|
|
Packit |
67cb25 |
0.0168317654213106405344439270765,
|
|
Packit |
67cb25 |
0.0746639850740189951912512662623,
|
|
Packit |
67cb25 |
-0.0313329787073628846871956180962,
|
|
Packit |
67cb25 |
-0.3011591259228349991008967259990,
|
|
Packit |
67cb25 |
0.0264992409453454699696117210896,
|
|
Packit |
67cb25 |
0.9516421218971785225243297231697,
|
|
Packit |
67cb25 |
-0.9516421218971785225243297231697,
|
|
Packit |
67cb25 |
-0.0264992409453454699696117210896,
|
|
Packit |
67cb25 |
0.3011591259228349991008967259990,
|
|
Packit |
67cb25 |
0.0313329787073628846871956180962,
|
|
Packit |
67cb25 |
-0.0746639850740189951912512662623,
|
|
Packit |
67cb25 |
-0.0168317654213106405344439270765,
|
|
Packit |
67cb25 |
0.0090632583037826525954698068873,
|
|
Packit |
67cb25 |
0.0030210861012608841984899356291
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
static const double h1_309[20] = { -0.0006797443727836989446602355165,
|
|
Packit |
67cb25 |
0.0020392331183510968339807065496,
|
|
Packit |
67cb25 |
0.0050603192196119810324706421788,
|
|
Packit |
67cb25 |
-0.0206189126411055346546938106687,
|
|
Packit |
67cb25 |
-0.0141127879301758447558029850103,
|
|
Packit |
67cb25 |
0.0991347824942321571990197448581,
|
|
Packit |
67cb25 |
0.0123001362694193142367090236328,
|
|
Packit |
67cb25 |
-0.3201919683607785695513833204624,
|
|
Packit |
67cb25 |
0.0020500227115698857061181706055,
|
|
Packit |
67cb25 |
0.9421257006782067372990864259380,
|
|
Packit |
67cb25 |
0.9421257006782067372990864259380,
|
|
Packit |
67cb25 |
0.0020500227115698857061181706055,
|
|
Packit |
67cb25 |
-0.3201919683607785695513833204624,
|
|
Packit |
67cb25 |
0.0123001362694193142367090236328,
|
|
Packit |
67cb25 |
0.0991347824942321571990197448581,
|
|
Packit |
67cb25 |
-0.0141127879301758447558029850103,
|
|
Packit |
67cb25 |
-0.0206189126411055346546938106687,
|
|
Packit |
67cb25 |
0.0050603192196119810324706421788,
|
|
Packit |
67cb25 |
0.0020392331183510968339807065496,
|
|
Packit |
67cb25 |
-0.0006797443727836989446602355165
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
static const double g2_309[20] = { 0.0006797443727836989446602355165,
|
|
Packit |
67cb25 |
0.0020392331183510968339807065496,
|
|
Packit |
67cb25 |
-0.0050603192196119810324706421788,
|
|
Packit |
67cb25 |
-0.0206189126411055346546938106687,
|
|
Packit |
67cb25 |
0.0141127879301758447558029850103,
|
|
Packit |
67cb25 |
0.0991347824942321571990197448581,
|
|
Packit |
67cb25 |
-0.0123001362694193142367090236328,
|
|
Packit |
67cb25 |
-0.3201919683607785695513833204624,
|
|
Packit |
67cb25 |
-0.0020500227115698857061181706055,
|
|
Packit |
67cb25 |
0.9421257006782067372990864259380,
|
|
Packit |
67cb25 |
-0.9421257006782067372990864259380,
|
|
Packit |
67cb25 |
0.0020500227115698857061181706055,
|
|
Packit |
67cb25 |
0.3201919683607785695513833204624,
|
|
Packit |
67cb25 |
0.0123001362694193142367090236328,
|
|
Packit |
67cb25 |
-0.0991347824942321571990197448581,
|
|
Packit |
67cb25 |
-0.0141127879301758447558029850103,
|
|
Packit |
67cb25 |
0.0206189126411055346546938106687,
|
|
Packit |
67cb25 |
0.0050603192196119810324706421788,
|
|
Packit |
67cb25 |
-0.0020392331183510968339807065496,
|
|
Packit |
67cb25 |
-0.0006797443727836989446602355165
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
static const double h2_3[20] = { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
|
|
Packit |
67cb25 |
0.1767766952966368811002110905262,
|
|
Packit |
67cb25 |
0.5303300858899106433006332715786,
|
|
Packit |
67cb25 |
0.5303300858899106433006332715786,
|
|
Packit |
67cb25 |
0.1767766952966368811002110905262,
|
|
Packit |
67cb25 |
0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
static const double g1_3[20] = { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
|
|
Packit |
67cb25 |
-0.1767766952966368811002110905262,
|
|
Packit |
67cb25 |
0.5303300858899106433006332715786,
|
|
Packit |
67cb25 |
-0.5303300858899106433006332715786,
|
|
Packit |
67cb25 |
0.1767766952966368811002110905262,
|
|
Packit |
67cb25 |
0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
static int
|
|
Packit |
67cb25 |
bspline_init (const double **h1, const double **g1,
|
|
Packit |
67cb25 |
const double **h2, const double **g2, size_t * nc,
|
|
Packit |
67cb25 |
size_t * offset, size_t member)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
switch (member)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
case 103:
|
|
Packit |
67cb25 |
*nc = 6;
|
|
Packit |
67cb25 |
*h1 = h1_103;
|
|
Packit |
67cb25 |
*g1 = &g1_1[2];
|
|
Packit |
67cb25 |
*h2 = &h2_1[2];
|
|
Packit |
67cb25 |
*g2 = g2_103;
|
|
Packit |
67cb25 |
break;
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
case 105:
|
|
Packit |
67cb25 |
*nc = 10;
|
|
Packit |
67cb25 |
*h1 = h1_105;
|
|
Packit |
67cb25 |
*g1 = g1_1;
|
|
Packit |
67cb25 |
*h2 = h2_1;
|
|
Packit |
67cb25 |
*g2 = g2_105;
|
|
Packit |
67cb25 |
break;
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
case 202:
|
|
Packit |
67cb25 |
*nc = 6;
|
|
Packit |
67cb25 |
*h1 = h1_202;
|
|
Packit |
67cb25 |
*g1 = &g1_2[6];
|
|
Packit |
67cb25 |
*h2 = &h2_2[6];
|
|
Packit |
67cb25 |
*g2 = g2_202;
|
|
Packit |
67cb25 |
break;
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
case 204:
|
|
Packit |
67cb25 |
*nc = 10;
|
|
Packit |
67cb25 |
*h1 = h1_204;
|
|
Packit |
67cb25 |
*g1 = &g1_2[4];
|
|
Packit |
67cb25 |
*h2 = &h2_2[4];
|
|
Packit |
67cb25 |
*g2 = g2_204;
|
|
Packit |
67cb25 |
break;
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
case 206:
|
|
Packit |
67cb25 |
*nc = 14;
|
|
Packit |
67cb25 |
*h1 = h1_206;
|
|
Packit |
67cb25 |
*g1 = &g1_2[2];
|
|
Packit |
67cb25 |
*h2 = &h2_2[2];
|
|
Packit |
67cb25 |
*g2 = g2_206;
|
|
Packit |
67cb25 |
break;
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
case 208:
|
|
Packit |
67cb25 |
*nc = 18;
|
|
Packit |
67cb25 |
*h1 = h1_208;
|
|
Packit |
67cb25 |
*g1 = g1_2;
|
|
Packit |
67cb25 |
*h2 = h2_2;
|
|
Packit |
67cb25 |
*g2 = g2_208;
|
|
Packit |
67cb25 |
break;
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
case 301:
|
|
Packit |
67cb25 |
*nc = 4;
|
|
Packit |
67cb25 |
*h1 = h1_301;
|
|
Packit |
67cb25 |
*g1 = &g1_3[8];
|
|
Packit |
67cb25 |
*h2 = &h2_3[8];
|
|
Packit |
67cb25 |
*g2 = g2_301;
|
|
Packit |
67cb25 |
break;
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
case 303:
|
|
Packit |
67cb25 |
*nc = 8;
|
|
Packit |
67cb25 |
*h1 = h1_303;
|
|
Packit |
67cb25 |
*g1 = &g1_3[6];
|
|
Packit |
67cb25 |
*h2 = &h2_3[6];
|
|
Packit |
67cb25 |
*g2 = g2_303;
|
|
Packit |
67cb25 |
break;
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
case 305:
|
|
Packit |
67cb25 |
*nc = 12;
|
|
Packit |
67cb25 |
*h1 = h1_305;
|
|
Packit |
67cb25 |
*g1 = &g1_3[4];
|
|
Packit |
67cb25 |
*h2 = &h2_3[4];
|
|
Packit |
67cb25 |
*g2 = g2_305;
|
|
Packit |
67cb25 |
break;
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
case 307:
|
|
Packit |
67cb25 |
*nc = 16;
|
|
Packit |
67cb25 |
*h1 = h1_307;
|
|
Packit |
67cb25 |
*g1 = &g1_3[2];
|
|
Packit |
67cb25 |
*h2 = &h2_3[2];
|
|
Packit |
67cb25 |
*g2 = g2_307;
|
|
Packit |
67cb25 |
break;
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
case 309:
|
|
Packit |
67cb25 |
*nc = 20;
|
|
Packit |
67cb25 |
*h1 = h1_309;
|
|
Packit |
67cb25 |
*g1 = g1_3;
|
|
Packit |
67cb25 |
*h2 = h2_3;
|
|
Packit |
67cb25 |
*g2 = g2_309;
|
|
Packit |
67cb25 |
break;
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
default:
|
|
Packit |
67cb25 |
return GSL_FAILURE;
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
*offset = 0;
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
return GSL_SUCCESS;
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
static int
|
|
Packit |
67cb25 |
bspline_centered_init (const double **h1, const double **g1,
|
|
Packit |
67cb25 |
const double **h2, const double **g2, size_t * nc,
|
|
Packit |
67cb25 |
size_t * offset, size_t member)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
switch (member)
|
|
Packit |
67cb25 |
{
|
|
Packit |
67cb25 |
case 103:
|
|
Packit |
67cb25 |
*nc = 6;
|
|
Packit |
67cb25 |
*h1 = h1_103;
|
|
Packit |
67cb25 |
*g1 = &g1_1[2];
|
|
Packit |
67cb25 |
*h2 = &h2_1[2];
|
|
Packit |
67cb25 |
*g2 = g2_103;
|
|
Packit |
67cb25 |
break;
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
case 105:
|
|
Packit |
67cb25 |
*nc = 10;
|
|
Packit |
67cb25 |
*h1 = h1_105;
|
|
Packit |
67cb25 |
*g1 = g1_1;
|
|
Packit |
67cb25 |
*h2 = h2_1;
|
|
Packit |
67cb25 |
*g2 = g2_105;
|
|
Packit |
67cb25 |
break;
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
case 202:
|
|
Packit |
67cb25 |
*nc = 6;
|
|
Packit |
67cb25 |
*h1 = h1_202;
|
|
Packit |
67cb25 |
*g1 = &g1_2[6];
|
|
Packit |
67cb25 |
*h2 = &h2_2[6];
|
|
Packit |
67cb25 |
*g2 = g2_202;
|
|
Packit |
67cb25 |
break;
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
case 204:
|
|
Packit |
67cb25 |
*nc = 10;
|
|
Packit |
67cb25 |
*h1 = h1_204;
|
|
Packit |
67cb25 |
*g1 = &g1_2[4];
|
|
Packit |
67cb25 |
*h2 = &h2_2[4];
|
|
Packit |
67cb25 |
*g2 = g2_204;
|
|
Packit |
67cb25 |
break;
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
case 206:
|
|
Packit |
67cb25 |
*nc = 14;
|
|
Packit |
67cb25 |
*h1 = h1_206;
|
|
Packit |
67cb25 |
*g1 = &g1_2[2];
|
|
Packit |
67cb25 |
*h2 = &h2_2[2];
|
|
Packit |
67cb25 |
*g2 = g2_206;
|
|
Packit |
67cb25 |
break;
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
case 208:
|
|
Packit |
67cb25 |
*nc = 18;
|
|
Packit |
67cb25 |
*h1 = h1_208;
|
|
Packit |
67cb25 |
*g1 = g1_2;
|
|
Packit |
67cb25 |
*h2 = h2_2;
|
|
Packit |
67cb25 |
*g2 = g2_208;
|
|
Packit |
67cb25 |
break;
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
case 301:
|
|
Packit |
67cb25 |
*nc = 4;
|
|
Packit |
67cb25 |
*h1 = h1_301;
|
|
Packit |
67cb25 |
*g1 = &g1_3[8];
|
|
Packit |
67cb25 |
*h2 = &h2_3[8];
|
|
Packit |
67cb25 |
*g2 = g2_301;
|
|
Packit |
67cb25 |
break;
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
case 303:
|
|
Packit |
67cb25 |
*nc = 8;
|
|
Packit |
67cb25 |
*h1 = h1_303;
|
|
Packit |
67cb25 |
*g1 = &g1_3[6];
|
|
Packit |
67cb25 |
*h2 = &h2_3[6];
|
|
Packit |
67cb25 |
*g2 = g2_303;
|
|
Packit |
67cb25 |
break;
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
case 305:
|
|
Packit |
67cb25 |
*nc = 12;
|
|
Packit |
67cb25 |
*h1 = h1_305;
|
|
Packit |
67cb25 |
*g1 = &g1_3[4];
|
|
Packit |
67cb25 |
*h2 = &h2_3[4];
|
|
Packit |
67cb25 |
*g2 = g2_305;
|
|
Packit |
67cb25 |
break;
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
case 307:
|
|
Packit |
67cb25 |
*nc = 16;
|
|
Packit |
67cb25 |
*h1 = h1_307;
|
|
Packit |
67cb25 |
*g1 = &g1_3[2];
|
|
Packit |
67cb25 |
*h2 = &h2_3[2];
|
|
Packit |
67cb25 |
*g2 = g2_307;
|
|
Packit |
67cb25 |
break;
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
case 309:
|
|
Packit |
67cb25 |
*nc = 20;
|
|
Packit |
67cb25 |
*h1 = h1_309;
|
|
Packit |
67cb25 |
*g1 = g1_3;
|
|
Packit |
67cb25 |
*h2 = h2_3;
|
|
Packit |
67cb25 |
*g2 = g2_309;
|
|
Packit |
67cb25 |
break;
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
default:
|
|
Packit |
67cb25 |
return GSL_FAILURE;
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
*offset = ((*nc) >> 1);
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
return GSL_SUCCESS;
|
|
Packit |
67cb25 |
}
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
static const gsl_wavelet_type bspline_type = {
|
|
Packit |
67cb25 |
"bspline",
|
|
Packit |
67cb25 |
&bspline_init
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
static const gsl_wavelet_type bspline_centered_type = {
|
|
Packit |
67cb25 |
"bspline-centered",
|
|
Packit |
67cb25 |
&bspline_centered_init
|
|
Packit |
67cb25 |
};
|
|
Packit |
67cb25 |
|
|
Packit |
67cb25 |
const gsl_wavelet_type *gsl_wavelet_bspline = &bspline_type;
|
|
Packit |
67cb25 |
const gsl_wavelet_type *gsl_wavelet_bspline_centered = &bspline_centered_type;
|