bspline.c

This library is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY without ev

  0.0746639850740189951912512662623,
  -0.0313329787073628846871956180962,
  -0.3011591259228349991008967259990,
  0.0264992409453454699696117210896,
  0.9516421218971785225243297231697,
  -0.9516421218971785225243297231697,
  -0.0264992409453454699696117210896,
  0.3011591259228349991008967259990,
  0.0313329787073628846871956180962,
  -0.0746639850740189951912512662623,
  -0.0168317654213106405344439270765,
  0.0090632583037826525954698068873,
  0.0030210861012608841984899356291
};

static const double h1_309[20] = { -0.0006797443727836989446602355165,
  0.0020392331183510968339807065496,
  0.0050603192196119810324706421788,
  -0.0206189126411055346546938106687,
  -0.0141127879301758447558029850103,
  0.0991347824942321571990197448581,
  0.0123001362694193142367090236328,
  -0.3201919683607785695513833204624,
  0.0020500227115698857061181706055,
  0.9421257006782067372990864259380,
  0.9421257006782067372990864259380,
  0.0020500227115698857061181706055,
  -0.3201919683607785695513833204624,
  0.0123001362694193142367090236328,
  0.0991347824942321571990197448581,
  -0.0141127879301758447558029850103,
  -0.0206189126411055346546938106687,
  0.0050603192196119810324706421788,
  0.0020392331183510968339807065496,
  -0.0006797443727836989446602355165
};

static const double g2_309[20] = { 0.0006797443727836989446602355165,
  0.0020392331183510968339807065496,
  -0.0050603192196119810324706421788,
  -0.0206189126411055346546938106687,
  0.0141127879301758447558029850103,
  0.0991347824942321571990197448581,
  -0.0123001362694193142367090236328,
  -0.3201919683607785695513833204624,
  -0.0020500227115698857061181706055,
  0.9421257006782067372990864259380,
  -0.9421257006782067372990864259380,
  0.0020500227115698857061181706055,
  0.3201919683607785695513833204624,
  0.0123001362694193142367090236328,
  -0.0991347824942321571990197448581,
  -0.0141127879301758447558029850103,
  0.0206189126411055346546938106687,
  0.0050603192196119810324706421788,
  -0.0020392331183510968339807065496,
  -0.0006797443727836989446602355165
};

static const double h2_3[20] = { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
  0.1767766952966368811002110905262,
  0.5303300858899106433006332715786,
  0.5303300858899106433006332715786,
  0.1767766952966368811002110905262,
  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0
};

static const double g1_3[20] = { 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
  -0.1767766952966368811002110905262,
  0.5303300858899106433006332715786,
  -0.5303300858899106433006332715786,
  0.1767766952966368811002110905262,
  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0
};

static int
bspline_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 103:
      *nc = 6;
      *h1 = h1_103;
      *g1 = &g1_1[2];
      *h2 = &h2_1[2];
      *g2 = g2_103;
      break;

    case 105:
      *nc = 10;
      *h1 = h1_105;
      *g1 = g1_1;
      *h2 = h2_1;
      *g2 = g2_105;
      break;

    case 202:
      *nc = 6;
      *h1 = h1_202;
      *g1 = &g1_2[6];
      *h2 = &h2_2[6];
      *g2 = g2_202;
      break;

    case 204:
      *nc = 10;
      *h1 = h1_204;
      *g1 = &g1_2[4];
      *h2 = &h2_2[4];
      *g2 = g2_204;
      break;

    case 206:
      *nc = 14;
      *h1 = h1_206;
      *g1 = &g1_2[2];
      *h2 = &h2_2[2];
      *g2 = g2_206;
      break;

    case 208:
      *nc = 18;
      *h1 = h1_208;
      *g1 = g1_2;
      *h2 = h2_2;
      *g2 = g2_208;
      break;

    case 301:
      *nc = 4;
      *h1 = h1_301;
      *g1 = &g1_3[8];
      *h2 = &h2_3[8];
      *g2 = g2_301;
      break;

    case 303:
      *nc = 8;
      *h1 = h1_303;
      *g1 = &g1_3[6];
      *h2 = &h2_3[6];
      *g2 = g2_303;
      break;

    case 305:
      *nc = 12;
      *h1 = h1_305;
      *g1 = &g1_3[4];
      *h2 = &h2_3[4];
      *g2 = g2_305;
      break;

    case 307:
      *nc = 16;
      *h1 = h1_307;
      *g1 = &g1_3[2];
      *h2 = &h2_3[2];
      *g2 = g2_307;
      break;

    case 309:
      *nc = 20;
      *h1 = h1_309;
      *g1 = g1_3;
      *h2 = h2_3;
      *g2 = g2_309;
      break;

    default:
      return GSL_FAILURE;
    }

  *offset = 0;

  return GSL_SUCCESS;
}

static int
bspline_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 103:
      *nc = 6;
      *h1 = h1_103;
      *g1 = &g1_1[2];
      *h2 = &h2_1[2];
      *g2 = g2_103;
      break;

    case 105:
      *nc = 10;
      *h1 = h1_105;
      *g1 = g1_1;
      *h2 = h2_1;
      *g2 = g2_105;
      break;

    case 202:
      *nc = 6;
      *h1 = h1_202;
      *g1 = &g1_2[6];
      *h2 = &h2_2[6];
      *g2 = g2_202;
      break;

    case 204:
      *nc = 10;
      *h1 = h1_204;
      *g1 = &g1_2[4];
      *h2 = &h2_2[4];
      *g2 = g2_204;
      break;

    case 206:
      *nc = 14;
      *h1 = h1_206;
      *g1 = &g1_2[2];
      *h2 = &h2_2[2];
      *g2 = g2_206;
      break;

    case 208:
      *nc = 18;
      *h1 = h1_208;
      *g1 = g1_2;
      *h2 = h2_2;
      *g2 = g2_208;
      break;

    case 301:
      *nc = 4;
      *h1 = h1_301;
      *g1 = &g1_3[8];
      *h2 = &h2_3[8];
      *g2 = g2_301;
      break;

    case 303:
      *nc = 8;
      *h1 = h1_303;
      *g1 = &g1_3[6];
      *h2 = &h2_3[6];
      *g2 = g2_303;
      break;

    case 305:
      *nc = 12;
      *h1 = h1_305;
      *g1 = &g1_3[4];
      *h2 = &h2_3[4];
      *g2 = g2_305;
      break;

    case 307:
      *nc = 16;
      *h1 = h1_307;
      *g1 = &g1_3[2];
      *h2 = &h2_3[2];
      *g2 = g2_307;
      break;

    case 309:
      *nc = 20;
      *h1 = h1_309;
      *g1 = g1_3;
      *h2 = h2_3;
      *g2 = g2_309;
      break;

    default:
      return GSL_FAILURE;
    }

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

  return GSL_SUCCESS;
}

static const gsl_wavelet_type bspline_type = {
  "bspline",
  &bspline_init
};

static const gsl_wavelet_type bspline_centered_type = { 
  "bspline-centered",
  &bspline_centered_init
};

const gsl_wavelet_type *gsl_wavelet_bspline = &bspline_type;
const gsl_wavelet_type *gsl_wavelet_bspline_centered = &bspline_centered_type;