constant.c

大师写的二代小波经典之作

/*
 *  -*- Mode: ANSI C -*-
 *  $Id: constant.c,v 1.5 1996/11/18 18:49:11 fernande Exp $
 *  $Source: /sgi.acct/sweldens/cvs/liftpack/Util/constant.c,v $
 *  Author: Peter Schroder
 *  Modified: Gabriel Fernandez
 *
 *  This file defines the values of the constants declared in constants.h.
 */
/* do not edit anything above this line */

/* FLWT header files */
#include <flwtdef.h>
#include <constant.h>

/* smallest value of error throughout the program */
#ifdef DOUBLE
Flt m_eps_f = (Flt)6E-8;
#else
Flt m_eps_f = (Flt)6E-16;
#endif

/* inverse powers of two */
Flt twomm[32] =
{
  (Flt)1./(Flt)( 1 << 0 ),
  (Flt)1./(Flt)( 1 << 1 ),
  (Flt)1./(Flt)( 1 << 2 ),
  (Flt)1./(Flt)( 1 << 3 ),
  (Flt)1./(Flt)( 1 << 4 ),
  (Flt)1./(Flt)( 1 << 5 ),
  (Flt)1./(Flt)( 1 << 6 ),
  (Flt)1./(Flt)( 1 << 7 ),
  (Flt)1./(Flt)( 1 << 8 ),
  (Flt)1./(Flt)( 1 << 9 ),
  (Flt)1./(Flt)( 1 << 10 ),
  (Flt)1./(Flt)( 1 << 11 ),
  (Flt)1./(Flt)( 1 << 12 ),
  (Flt)1./(Flt)( 1 << 13 ),
  (Flt)1./(Flt)( 1 << 14 ),
  (Flt)1./(Flt)( 1 << 15 ),
  (Flt)1./(Flt)( 1 << 16 ),
  (Flt)1./(Flt)( 1 << 17 ),
  (Flt)1./(Flt)( 1 << 18 ),
  (Flt)1./(Flt)( 1 << 19 ),
  (Flt)1./(Flt)( 1 << 20 ),
  (Flt)1./(Flt)( 1 << 21 ),
  (Flt)1./(Flt)( 1 << 22 ),
  (Flt)1./(Flt)( 1 << 23 ),
  (Flt)1./(Flt)( 1 << 24 ),
  (Flt)1./(Flt)( 1 << 25 ),
  (Flt)1./(Flt)( 1 << 26 ),
  (Flt)1./(Flt)( 1 << 27 ),
  (Flt)1./(Flt)( 1 << 28 ),
  (Flt)1./(Flt)( 1 << 29 ),
  (Flt)1./(Flt)( 1 << 30 ),
  (Flt)1./(Flt)( 1 << 31 ) };

/* inverse powers of sqrt(2) */
Flt twomhalfm[64] =
{
  (Flt)1.,
  (Flt)0.70710678118654752,
  (Flt)0.5,
  (Flt)0.35355339059327376,
  (Flt)0.25,
  (Flt)0.17677669529663688,
  (Flt)0.125,
  (Flt)0.08838834764831844,
  (Flt)0.0625,
  (Flt)0.04419417382415922,
  (Flt)0.03125,
  (Flt)0.02209708691207961,
  (Flt)0.015625,
  (Flt)0.011048543456039805,
  (Flt)0.0078125,
  (Flt)0.005524271728019903,
  (Flt)0.00390625,
  (Flt)0.002762135864009951,
  (Flt)0.001953125,
  (Flt)0.001381067932004976,
  (Flt)0.0009765625,
  (Flt)0.0006905339660024878,
  (Flt)0.00048828125,
  (Flt)0.0003452669830012439,
  (Flt)0.000244140625,
  (Flt)0.000172633491500622,
  (Flt)0.0001220703125,
  (Flt)0.00008631674575031098,
  (Flt)0.00006103515625,
  (Flt)0.00004315837287515549,
  (Flt)0.000030517578125,
  (Flt)0.00002157918643757774,
  (Flt)0.0000152587890625,
  (Flt)0.00001078959321878887,
  (Flt)7.62939453125e-6,
  (Flt)5.394796609394436e-6,
  (Flt)3.814697265625e-6,
  (Flt)2.697398304697218e-6,
  (Flt)1.9073486328125e-6,
  (Flt)1.348699152348609e-6,
  (Flt)9.5367431640625e-7,
  (Flt)6.74349576174305e-7,
  (Flt)4.76837158203125e-7,
  (Flt)3.371747880871523e-7,
  (Flt)2.384185791015625e-7,
  (Flt)1.685873940435761e-7,
  (Flt)1.1920928955078125e-7,
  (Flt)8.42936970217881e-8,
  (Flt)5.960464477539062e-8,
  (Flt)4.214684851089403e-8,
  (Flt)2.9802322387695312e-8,
  (Flt)2.107342425544702e-8,
  (Flt)1.4901161193847656e-8,
  (Flt)1.053671212772351e-8,
  (Flt)7.450580596923828e-9,
  (Flt)5.268356063861754e-9,
  (Flt)3.725290298461914e-9,
  (Flt)2.634178031930877e-9,
  (Flt)1.862645149230957e-9,
  (Flt)1.317089015965438e-9,
  (Flt)9.313225746154785e-10,
  (Flt)6.58544507982719e-10,
  (Flt)4.656612873077393e-10,
  (Flt)3.292722539913596e-10
};

/* ordinary powers of 2 */
Flt twom[32] =
{
  (Flt)(1 << 0),
  (Flt)(1 << 1),
  (Flt)(1 << 2),
  (Flt)(1 << 3),
  (Flt)(1 << 4),
  (Flt)(1 << 5),
  (Flt)(1 << 6),
  (Flt)(1 << 7),
  (Flt)(1 << 8),
  (Flt)(1 << 9),
  (Flt)(1 << 10),
  (Flt)(1 << 11),
  (Flt)(1 << 12),
  (Flt)(1 << 13),
  (Flt)(1 << 14),
  (Flt)(1 << 15),
  (Flt)(1 << 16),
  (Flt)(1 << 17),
  (Flt)(1 << 18),
  (Flt)(1 << 19),
  (Flt)(1 << 20),
  (Flt)(1 << 21),
  (Flt)(1 << 22),
  (Flt)(1 << 23),
  (Flt)(1 << 24),
  (Flt)(1 << 25),
  (Flt)(1 << 26),
  (Flt)(1 << 27),
  (Flt)(1 << 28),
  (Flt)(1 << 29),
  (Flt)(1 << 30),
  (Flt)(1 << 31) };

/* ordinary powers of sqrt(2) */
Flt twohalfm[64] =
{
  (Flt)1.,
  (Flt)1.4142135623730950488,
  (Flt)2.,
  (Flt)2.8284271247461900976,
  (Flt)4.,
  (Flt)5.656854249492380195,
  (Flt)8.,
  (Flt)11.31370849898476039,
  (Flt)16.,
  (Flt)22.627416997969520781,
  (Flt)32.,
  (Flt)45.25483399593904156,
  (Flt)64.,
  (Flt)90.50966799187808312,
  (Flt)128.,
  (Flt)181.01933598375616625,
  (Flt)256.,
  (Flt)362.0386719675123325,
  (Flt)512.,
  (Flt)724.077343935024665,
  (Flt)1024.,
  (Flt)1448.15468787004933,
  (Flt)2048.,
  (Flt)2896.30937574009866,
  (Flt)4096.,
  (Flt)5792.61875148019732,
  (Flt)8192.,
  (Flt)11585.23750296039464,
  (Flt)16384.,
  (Flt)23170.47500592078928,
  (Flt)32768.,
  (Flt)46340.95001184157856,
  (Flt)65536.,
  (Flt)92681.9000236831571,
  (Flt)131072.,
  (Flt)185363.8000473663142,
  (Flt)262144.,
  (Flt)370727.6000947326285,
  (Flt)524288.,
  (Flt)741455.200189465257,
  (Flt)1.048576e6,
  (Flt)1.482910400378930514e6,
  (Flt)2.097152e6,
  (Flt)2.965820800757861028e6,
  (Flt)4.194304e6,
  (Flt)5.931641601515722056e6,
  (Flt)8.388608e6,
  (Flt)1.186328320303144411e7,
  (Flt)1.6777216e7,
  (Flt)2.372656640606288822e7,
  (Flt)3.3554432e7,
  (Flt)4.745313281212577644e7,
  (Flt)6.7108864e7,
  (Flt)9.49062656242515529e7,
  (Flt)1.34217728e8,
  (Flt)1.898125312485031058e8,
  (Flt)2.68435456e8,
  (Flt)3.796250624970062116e8,
  (Flt)5.36870912e8,
  (Flt)7.59250124994012423e8,
  (Flt)1.073741824e9,
  (Flt)1.518500249988024846e9,
  (Flt)2.147483648e9,
  (Flt)3.037000499976049692e9
};

/* integer powers of two */
int itwom[31] = {
    1 << 0,
    1 << 1,
    1 << 2,
    1 << 3,
    1 << 4,
    1 << 5,
    1 << 6,
    1 << 7,
    1 << 8,
    1 << 9,
    1 << 10,
    1 << 11,
    1 << 12,
    1 << 13,
    1 << 14,
    1 << 15,
    1 << 16,
    1 << 17,
    1 << 18,
    1 << 19,
    1 << 20,
    1 << 21,
    1 << 22,
    1 << 23,
    1 << 24,
    1 << 25,
    1 << 26,
    1 << 27,
    1 << 28,
    1 << 29,
    1 << 30,
};