rtmrartwt.cpp

Go to the documentation of this file.

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

#ifdef STORM_USEF
#include "rtmrartwt_.h"
#include "drtmrartwt_.h"
#include "irtmrartwt_.h"
#include "artmrartwt_.h"
#include "dartmrartwt_.h"
#include "iartmrartwt_.h"
#endif

#include "stormdef.h"
#define INVSQRT2 STORMDEF_MATHCONST_INVSQRT2

void rtmrartwt( 
          const int& Msize,
          const int& mmax,
                int& mrastep,
                float* mracoef,
          const float& datavalue
) { 
#ifdef STORM_USEF
   rtmrartwt_( &Msize, &mmax, &mrastep, mracoef, &datavalue );
#else
   mrastep++;
   int mcycle = 0;
   while ( ( mrastep & (1<<mcycle) ) && ( mcycle < mmax ) ) { mcycle++; }
   float tmp=datavalue;
   int m = 1;
   float tmpold;
   while ( m <= mcycle ) {
      tmpold = mracoef[ 2*(m-1) ];
      mracoef[ 2*(m-1) ] = tmp;
      mracoef[ 2*m +1 ] = (tmpold-tmp) * INVSQRT2;
      tmp = (tmpold+tmp) * INVSQRT2;
      m++;
   }
   mracoef[ 2*(m-1) ] = tmp;
#endif
} 


void rtmrartwt( 
          const int& Msize,
          const int& mmax,
                int& mrastep,
                double* mracoef,
          const double& datavalue
) {
#ifdef STORM_USEF 
   drtmrartwt_( &Msize, &mmax, &mrastep, mracoef, &datavalue );
#else
   mrastep++;
   int mcycle = 0;
   while ( (mrastep & (1<<mcycle) ) && ( mcycle < mmax ) ) { mcycle++; }
   double tmp = datavalue;
   int m = 1;
   double tmpold,tmps,tmpw;
   while ( m <= mcycle ) {
      tmpold = mracoef[ 2*(m-1) ];
      mracoef[ 2*(m-1) ] = tmp;
      mracoef[ 2*m + 1 ] = (tmpold-tmp) * INVSQRT2;
      tmp = (tmpold+tmp) * INVSQRT2;
      m++;
   }
   mracoef[ 2*(m-1) ] = tmp;
#endif
}


void rtmrartwt( 
          const int& Msize,
          const int& mmax,
                int& mrastep,
                int* mracoef,
          const int& datavalue
) {
#ifdef STORM_USEF 
   irtmrartwt_( &Msize, &mmax, &mrastep, mracoef, &datavalue );
#else
   mrastep++;
   int mcycle = 0;
   while ( (mrastep & (1<<mcycle) ) && ( mcycle < mmax ) ) { mcycle++; }
   int tmp = datavalue;
   int m = 1;
   int tmpold,tmps,tmpw;
   while ( m <= mcycle ) {
      tmpold = mracoef[ 2*(m-1) ];
      mracoef[ 2*(m-1) ] = tmp;
      mracoef[ 2*m + 1 ] = tmpold-tmp;
      tmp = tmpold + tmp;
      if ( (tmp >= 0) || ((tmp & 0x00000001) != 0) ) {
         tmp = tmp/2;
      } else {
         tmp = tmp/2 -1;
      }
      m++;
   }
   mracoef[ 2*(m-1) ] = tmp;
#endif
}


void rtmrartwt(
          const int& Asize,
          const int& Msize,
          const int& amax,
          const int& mmax,
                int& mrastep,
                float* mracoef,
          const float* datavalue
) {
#ifdef STORM_USEF
   artmrartwt_( &Asize, &Msize, &amax, &mmax, &mrastep, mracoef, datavalue);
#else
   mrastep++;
   int mcycle = 0;
   while ( ( mrastep & (1<<mcycle) ) && ( mcycle < mmax ) ) { mcycle++; }
   int m;
   float tmp,tmpold,tmps,tmpw;
   for(int a = 0; a <= amax; a++ ) {
      tmp = datavalue[a];
      m = 1;
      while ( m <= mcycle ) {
         tmpold = mracoef[ a*2*Msize + 2*(m-1) ];
         mracoef[ a*2*Msize + 2*(m-1) ] = tmp;
         mracoef[ a*2*Msize + 2*m + 1 ] = (tmpold-tmp) * INVSQRT2;
         tmp = (tmpold+tmp) * INVSQRT2;
         m++;
      }
      mracoef[ a*2*Msize + 2*(m-1) ] = tmp;
   }
#endif
}


void rtmrartwt(
          const int& Asize,
          const int& Msize,
          const int& amax,
          const int& mmax,
                int& mrastep,
                double* mracoef,
          const double* datavalue
) {
#ifdef STORM_USEF
   dartmrartwt_( &Asize, &Msize, &amax, &mmax, &mrastep, mracoef, datavalue );
#else
   mrastep++;
   int mcycle = 0;
   while ( ( mrastep & (1<<mcycle) ) && ( mcycle < mmax ) ) { mcycle++; }
   int m;
   double tmp,tmpold,tmps,tmpw;
   for(int a = 0; a <= amax; a++) {
      tmp = datavalue[a];
      m = 1;
      while ( m <= mcycle ) {
         tmpold = mracoef[ a*2*Msize + 2*(m-1) ];
         mracoef[ a*2*Msize + 2*(m-1) ] = tmp;
         mracoef[ a*2*Msize + 2*m + 1 ] = (tmpold-tmp) * INVSQRT2;
         tmp = (tmpold+tmp) * INVSQRT2;
         m++;
      }
      mracoef[ a*2*Msize + 2*(m-1) ] = tmp;
   }
#endif
}

 
void rtmrartwt(
          const int& Asize,
          const int& Msize,
          const int& amax,
          const int& mmax,
                int& mrastep,
                int* mracoef,
          const int* datavalue
) {
#ifdef STORM_USEF
   iartmrartwt_( &Asize, &Msize, &amax, &mmax, &mrastep, mracoef, datavalue );
#else
   mrastep++;
   int mcycle = 0;
   while ( ( mrastep & (1<<mcycle) ) && ( mcycle < mmax ) ) { mcycle++; }
   int m;
   int tmp,tmpold;
   for(int a = 0; a <= amax; a++) {
      tmp = datavalue[a];
      m = 1;
      while ( m <= mcycle ) {
         tmpold = mracoef[ a*2*Msize + 2*(m-1) ];
         mracoef[ a*2*Msize + 2*(m-1) ] = tmp;
         mracoef[ a*2*Msize + 2*m + 1 ] = tmpold-tmp;
         tmp = tmpold + tmp;
         if ( (tmp >= 0) || ((tmp & 0x00000001) != 0) ) {
            tmp = tmp/2;
         } else {
            tmp = tmp/2-1;
         }
         m++;
      }
      mracoef[ a*2*Msize + 2*(m-1) ] = tmp;
   }
#endif
}
 

---