hifwt.cpp

Go to the documentation of this file.

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

#ifdef STORM_USEF
#include "hifwt_.h"
#include "dhifwt_.h"
#include "ihifwt_.h"
#include "hifwtws_.h"
#include "dhifwtws_.h"
#include "ihifwtws_.h"
#endif

#include "stormdef.h"
#define INVSQRT2 STORMDEF_MATHCONST_INVSQRT2

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

#ifdef STORM_USEF

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

#else

   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++ ) {
         sdata[2*i] = ( xx[i] + xx[size2+i] ) * INVSQRT2;
         sdata[2*i+1] = ( xx[i] - xx[size2+i] ) * INVSQRT2;
      }
      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 hifwt( 
          const int& Sizef,
          const int& ns,
          const int& size,
          const double* wdata,
                double* sdata
) { 

#ifdef STORM_USEF

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

#else

   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++ ) {
         sdata[2*i] = ( xx[i] + xx[size2+i] ) * INVSQRT2;
         sdata[2*i+1] = ( xx[i] - xx[size2+i] ) * INVSQRT2;
      }
      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 hifwt( 
          const int& Sizef,
          const int& ns,
          const int& size,
          const int* wdata,
                int* sdata
) { 

#ifdef STORM_USEF

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

#else

   int size2;
   int size1 = ns;
   int *xx = new int[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++ ) {
         sdata[2*i] = ( xx[i] + xx[size2+i] ) / 2;
         sdata[2*i+1] = xx[i] - xx[size2+i];
      }
      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 hifwt( 
          const int& Sizef,
          const int& ns,
          const int& size,
          const float* wdata,
                float* sdata,
                float* ws
) { 

#ifdef STORM_USEF

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

#else

   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++ ) {
         sdata[2*i] = ( ws[i] + ws[size2+i] ) * INVSQRT2;
         sdata[2*i+1] = ( ws[i] - ws[size2+i] ) * INVSQRT2;
      }
      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 hifwt( 
          const int& Sizef,
          const int& ns,
          const int& size,
          const double* wdata,
                double* sdata,
                double* ws
) { 

#ifdef STORM_USEF

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

#else

   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++ ) {
         sdata[2*i] = ( ws[i] + ws[size2+i] ) * INVSQRT2;
         sdata[2*i+1] = ( ws[i] - ws[size2+i] ) * INVSQRT2;
      }
      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 hifwt( 
          const int& Sizef,
          const int& ns,
          const int& size,
          const int* wdata,
                int* sdata,
                int* ws
) { 

#ifdef STORM_USEF

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

#else

   int size2;
   int size1 = ns;

   if ( size < ns ) return 1;

   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++ ) {
         sdata[2*i] = ( ws[i] + ws[size2+i] ) / 2;
         sdata[2*i+1] = ws[i] - ws[size2+i];
      }
      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 INVSQRT2

---