pd4ifwt.cpp

Go to the documentation of this file.

// pd4ifwt.cpp
//////////////////////////////////////////////////////////////////////
//          (c) Copyright 1994-2002 Brown Deer Technology, LLC.
//                        All rights reserved.
//////////////////////////////////////////////////////////////////////

#ifdef STORM_USEF
#include "pd4ifwt_.h"
#include "dpd4ifwt_.h"
#include "pd4ifwtws_.h"
#include "dpd4ifwtws_.h"
#endif

#include "stormdef.h"
#define H0 STORMDEF_MATHCONST_D4H0
#define H1 STORMDEF_MATHCONST_D4H1
#define H2 STORMDEF_MATHCONST_D4H2
#define H3 STORMDEF_MATHCONST_D4H3

int pd4ifwt( 
          const int& Sizef,
          const int& ns,
          const int& ks,
          const int& size,
          const float* wdata,
                float* sdata
) { 

#ifdef STORM_USEF

   int iflag;
   pd4ifwt_( &Sizef, &ns, &ks, &size, wdata, sdata, &iflag );
   return iflag;

#else

   int i1;
   int size2;
   int size1 = ns;
   float* xx = new float[Sizef];

#ifdef STORM_FWT_CCHECKUSAGE
   if ( size < ns ) return 1;
#endif

   if ( size == ns ) {
      for( int i = 0; i < size; i++ ) { sdata[i] = wdata[i]; }
      return 0; 
   }
   for( int i = 0; i < size; i++ ) { xx[i] = wdata[i]; }
   while ( size1 < size ) {
      size2 = size1;
      size1 = 2 * size1;
      for( int i = 0; i < size2; i++ ) {
         i1 = (i-1+size2)%size2;
         sdata[2*i] = H0*xx[i] + H2*xx[i1] 
           + H3*xx[size2+i] + H1*xx[size2+i1];
         sdata[2*i+1] = H1*xx[i] + H3*xx[i1] 
           - H2*xx[size2+i] - H0*xx[size2+i1];
      }
      for( int i = 0; i < size1; i++ ) { xx[i] = sdata[i]; }
   }
   delete [] xx;

#ifdef STORM_FWT_CCHECKUSAGE
   if ( size1 != size ) return 2;
#endif

   return 0;

#endif

} 



int pd4ifwt( 
          const int& Sizef,
          const int& ns,
          const int& ks,
          const int& size,
          const double* wdata,
                double* sdata
) { 

#ifdef STORM_USEF

   int iflag;
   dpd4ifwt_( &Sizef, &ns, &ks, &size, wdata, sdata, &iflag );
   return iflag;

#else

   int i1;
   int size2;
   int size1 = ns;
   double *xx = new double[Sizef];

#ifdef STORM_FWT_CCHECKUSAGE
   if ( size < ns ) return 1;
#endif

   if ( size == ns ) {
      for( int i = 0; i < size; i++ ) { sdata[i] = wdata[i]; }
      return 0;
   }
   for( int i = 0; i < size; i++ ) { xx[i] = wdata[i]; }
   while ( size1 < size ) {
      size2 = size1; 
      size1 = 2 * size1;
      for( int i = 0; i < size2; i++ ) {
         i1 = (i-1+size2)%size2;
         sdata[2*i] = H0*xx[i] + H2*xx[i1] 
           + H3*xx[size2+i] + H1*xx[size2+i1];
         sdata[2*i+1] = H1*xx[i] + H3*xx[i1] 
           - H2*xx[size2+i] - H0*xx[size2+i1];
      }
      for( int i = 0; i < size1; i++ ) { xx[i] = sdata[i]; }
   }
   delete [] xx;

#ifdef STORM_FWT_CCHECKUSAGE
   if ( size1 != size ) return 2;
#endif

   return 0;

#endif

} 


int pd4ifwt( 
          const int& Sizef,
          const int& ns,
          const int& ks,
          const int& size,
          const float* wdata,
                float* sdata,
                float* ws
) { 

#ifdef STORM_USEF

   int iflag;
   pd4ifwtws_( &Sizef, &ns, &ks, &size, wdata, sdata, ws, &iflag );
   return iflag;

#else

   int i1;
   int size2;
   int size1 = ns;

#ifdef STORM_FWT_CCHECKUSAGE
   if ( size < ns ) return 1;
#endif

   if ( size == ns ) {
      for( int i = 0; i < size; i++ ) { sdata[i] = wdata[i]; }
      return 0; 
   }
   for( int i = 0; i < size; i++ ) { ws[i] = wdata[i]; }
   while ( size1 < size ) {
      size2 = size1;
      size1 = 2 * size1;
      for( int i = 0; i < size2; i++ ) {
         i1 = (i-1+size2)%size2;
         sdata[2*i] = H0*ws[i] + H2*ws[i1] 
           + H3*ws[size2+i] + H1*ws[size2+i1];
         sdata[2*i+1] = H1*ws[i] + H3*ws[i1] 
           - H2*ws[size2+i] - H0*ws[size2+i1];
      }
      for( int i = 0; i < size1; i++ ) { ws[i] = sdata[i]; }
   }

#ifdef STORM_FWT_CCHECKUSAGE
   if ( size1 != size ) return 2;
#endif

   return 0;

#endif

} 


int pd4ifwt( 
          const int& Sizef,
          const int& ns,
          const int& ks,
          const int& size,
          const double* wdata,
                double* sdata,
                double* ws
) { 

#ifdef STORM_USEF

   int iflag;
   dpd4ifwtws_( &Sizef, &ns, &ks, &size, wdata, sdata, ws, &iflag );
   return iflag;

#else

   int i1;
   int size2;
   int size1 = ns;

#ifdef STORM_FWT_CCHECKUSAGE
   if ( size < ns ) return 1;
#endif

   if ( size == ns ) {
      for( int i = 0; i < size; i++ ) { sdata[i] = wdata[i]; }
      return 0;
   }
   for( int i = 0; i < size; i++ ) { ws[i] = wdata[i]; }
   while ( size1 < size ) {
      size2 = size1; 
      size1 = 2 * size1;
      for( int i = 0; i < size2; i++ ) {
         i1 = (i-1+size2)%size2;
         sdata[2*i] = H0*ws[i] + H2*ws[i1] 
           + H3*ws[size2+i] + H1*ws[size2+i1];
         sdata[2*i+1] = H1*ws[i] + H3*ws[i1] 
           - H2*ws[size2+i] - H0*ws[size2+i1];
      }
      for( int i = 0; i < size1; i++ ) { ws[i] = sdata[i]; }
   }

#ifdef STORM_FWT_CCHECKUSAGE
   if ( size1 != size ) return 2;
#endif

   return 0;

#endif

} 


#undef H0
#undef H1
#undef H2
#undef H3

---