pd4fwt.cpp

Go to the documentation of this file.

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

#ifdef STORM_USEF
#include "pd4fwt_.h"
#include "dpd4fwt_.h"
#include "pd4fwtws_.h"
#include "dpd4fwtws_.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 pd4fwt( 
          const int& Sizef,
          const int& ns,
          const int& ks,
          const int& size,
          const float* sdata,
                float* wdata
) { 

#ifdef STORM_USEF

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

#else

   int j;
   int j1,j2,j3;
   int size2;
   int size1 = size;
   float *xx = new float[Sizef];

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

   for( int i = 0; i < size1; i++ ) { wdata[i] = sdata[i]; }
   if ( size1 == ns ) return 0;
   while ( size1 > ns ) {
      size2 = size1 / 2;
      for( int i = 0; i < size2; i++ ) {
         j = 2 * i;
         j1 = (j+1)%size1;
         j2 = (j+2)%size1;
         j3 = (j+3)%size1;
         xx[i] = H0*wdata[j] + H1*wdata[j+1] 
           + H2*wdata[j2] + H3*wdata[j3];
         xx[size2+i] = H3*wdata[j] - H2*wdata[j1] 
           + H1*wdata[j2] + H0*wdata[j3];
      }
      for( int i = 0; i < size1; i++ ) { wdata[i] = xx[i]; }
      size1 = size2;
   }
   delete [] xx;

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

   return 0;

#endif

} 



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

#ifdef STORM_USEF

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

#else

   int j;
   int j1,j2,j3;
   int size2;
   int size1 = size;
   double *xx = new double[Sizef];

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

   for( int i = 0; i < size1; i++ ) { wdata[i] = sdata[i]; }
   if ( size1 == ns ) return 0;
   while ( size1 > ns ) {
      size2 = size1 / 2;
      for( int i = 0; i < size2; i++ ) {
         j = 2 * i;
         j1 = (j+1)%size1;
         j2 = (j+2)%size1;
         j3 = (j+3)%size1;
         xx[i] = H0*wdata[j] + H1*wdata[j+1] 
           + H2*wdata[j2] + H3*wdata[j3];
         xx[size2+i] = H3*wdata[j] - H2*wdata[j1] 
           + H1*wdata[j2] + H0*wdata[j3];
      }
      for( int i = 0; i < size1; i++ ) { wdata[i] = xx[i]; }
      size1 = size2;
   }
   delete [] xx;

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

   return 0;

#endif

} 



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

#ifdef STORM_USEF

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

#else

   int j;
   int j1,j2,j3;
   int size2;
   int size1 = size;
   if ( size1 < ns ) return 1;
   for( int i = 0; i < size1; i++ ) { wdata[i] = sdata[i]; }

#ifdef STORM_FWT_CCHECKUSAGE
   if ( size1 == ns ) return 0;
#endif

   while ( size1 > ns ) {
      size2 = size1 / 2;
      for( int i = 0; i < size2; i++ ) {
         j = 2 * i;
         j1 = (j+1)%size1;
         j2 = (j+2)%size1;
         j3 = (j+3)%size1;
         ws[i] = H0*wdata[j] + H1*wdata[j+1] 
           + H2*wdata[j2] + H3*wdata[j3];
         ws[size2+i] = H3*wdata[j] - H2*wdata[j1] 
           + H1*wdata[j2] + H0*wdata[j3];
      }
      for( int i = 0; i < size1; i++ ) { wdata[i] = ws[i]; }
      size1 = size2;
   }

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

   return 0;

#endif

} 



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

#ifdef STORM_USEF

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

#else

   int j;
   int j1,j2,j3;
   int size2;
   int size1 = size;

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

   for( int i = 0; i < size1; i++ ) { wdata[i] = sdata[i]; }
   if ( size1 == ns ) return 0;
   while ( size1 > ns ) {
      size2 = size1 / 2;
      for( int i = 0; i < size2; i++ ) {
         j = 2 * i;
         j1 = (j+1)%size1;
         j2 = (j+2)%size1;
         j3 = (j+3)%size1;
         ws[i] = H0*wdata[j] + H1*wdata[j+1] 
           + H2*wdata[j2] + H3*wdata[j3];
         ws[size2+i] = H3*wdata[j] - H2*wdata[j1] 
           + H1*wdata[j2] + H0*wdata[j3];
      }
      for( int i = 0; i < size1; i++ ) { wdata[i] = ws[i]; }
      size1 = size2;
   }

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

   return 0;

#endif

} 

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

---