rtmrartjiwt.cpp

Go to the documentation of this file.

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

#ifdef STORM_USEF
#include "rtmrartjiwt_.h"
#include "drtmrartjiwt_.h"
#include "irtmrartjiwt_.h"
#include "artmrartjiwt_.h"
#include "dartmrartjiwt_.h"
#include "iartmrartjiwt_.h"
#endif

#ifdef USE_STD
#include <iostream>
using namespace std;
#else
#include <iostream.h>
#endif

#include "stormdef.h"
#define INVSQRT2 STORMDEF_MATHCONST_INVSQRT2

void rtmrartjiwt( 
          const int& Msize,
          const int& mmax,
                int& mrastep,
                float* mracoef,
                float& data
) { 
#ifdef STORM_USEF
   rtmrartjiwt_( &Msize, &mmax, &mrastep, mracoef, &data );
#else
   ++mrastep;
   for(int m = (mmax-1); m >= 0; --m ) {
      if ( (mrastep & (1<<m)) != 0 ) {
         mracoef[2*m] = ( mracoef[2*(m+1)] - mracoef[2*(m+1)+1] ) * INVSQRT2;
      } else {
         mracoef[2*m] = ( mracoef[2*(m+1)] + mracoef[2*(m+1)+1] ) * INVSQRT2;
      }
   }
   data = mracoef[0];
#endif
}


void rtmrartjiwt( 
          const int& Msize,
          const int& mmax,
                int& mrastep,
                double* mracoef,
                double& data
) { 
#ifdef STORM_USEF
   drtmrartjiwt_( &Msize, &mmax, &mrastep, mracoef, &data );
#else
   ++mrastep;
   for(int m = (mmax-1); m >= 0; --m ) {
      if ( (mrastep & (1<<m)) != 0 ) {
         mracoef[2*m] = ( mracoef[2*(m+1)] - mracoef[2*(m+1)+1] ) * INVSQRT2;
      } else {
         mracoef[2*m] = ( mracoef[2*(m+1)] + mracoef[2*(m+1)+1] ) * INVSQRT2;
      }
   }
   data = mracoef[0];
#endif
} 


void rtmrartjiwt( 
          const int& Msize,
          const int& mmax,
                int& mrastep,
                int* mracoef,
                int& data
) { 
#ifdef STORM_USEF
   irtmrartjiwt_( &Msize, &mmax, &mrastep, mracoef, &data );
#else
   int tmp;
   ++mrastep;
   for(int m = (mmax-1); m >= 0; --m ) {
      if ( (mrastep & (1<<m)) != 0 ) { //right
         tmp = mracoef[2*(m+1)+1];
         if ( (tmp >= 0) || ((tmp & 0x00000001) == 0) ) {
            mracoef[2*m] = mracoef[2*(m+1)] - tmp/2;
         } else {
            mracoef[2*m] = mracoef[2*(m+1)] - (tmp/2-1);
         }
      } else { //left
         tmp = mracoef[2*(m+1)+1] + 1;
         if ( (tmp >= 0) || ((tmp & 0x00000001) == 0) ) {
            mracoef[2*m] = mracoef[2*(m+1)] + tmp/2;
         } else {
            mracoef[2*m] = mracoef[2*(m+1)] + (tmp/2-1);
         }
      }
   }
   data = mracoef[0];
#endif
} 


void rtmrartjiwt( 
          const int& Asize,
          const int& Msize,
          const int& amax,
          const int& mmax,
                int& mrastep,
                float* mracoef,
                float* data
) { 
#ifdef STORM_USEF
   artmrartjiwt_( &Asize,&Msize, &amax,&mmax, &mrastep, mracoef, data );
#else
   ++mrastep;
   for( int a =0; a <= amax; a++ ) {
      for(int m = (mmax-1); m >= 0; --m ) {
         if ( (mrastep & (1<<m)) != 0 ) {
            mracoef[ a*2*Msize + 2*m] = ( mracoef[ a*2*Msize + 2*(m+1)] 
              - mracoef[ a*2*Msize + 2*(m+1)+1] ) * INVSQRT2;
         } else {
            mracoef[ a*2*Msize + 2*m] = ( mracoef[ a*2*Msize + 2*(m+1)] 
              + mracoef[ a*2*Msize + 2*(m+1)+1] ) * INVSQRT2;
         }
      }
      data[a] = mracoef[a*2*Msize];
   }
#endif
} 


void rtmrartjiwt( 
          const int& Asize,
          const int& Msize,
          const int& amax,
          const int& mmax,
                int& mrastep,
                double* mracoef,
                double* data
) { 
#ifdef STORM_USEF
   dartmrartjiwt_( &Asize,&Msize, &amax,&mmax, &mrastep, mracoef, data );
#else
   ++mrastep;
   for( int a =0; a <= amax; a++ ) {
      for(int m = (mmax-1); m >= 0; --m ) {
         if ( (mrastep & (1<<m)) != 0 ) {
            mracoef[ a*2*Msize + 2*m] = ( mracoef[ a*2*Msize + 2*(m+1)] 
              - mracoef[ a*2*Msize + 2*(m+1)+1] ) * INVSQRT2;
         } else {
            mracoef[ a*2*Msize + 2*m] = ( mracoef[ a*2*Msize + 2*(m+1)] 
              + mracoef[ a*2*Msize + 2*(m+1)+1] ) * INVSQRT2;
         }
      }
      data[a] = mracoef[a*2*Msize];
   }
#endif
} 


void rtmrartjiwt( 
          const int& Asize,
          const int& Msize,
          const int& amax,
          const int& mmax,
                int& mrastep,
                int* mracoef,
                int* data
) { 
#ifdef STORM_USEF
   iartmrartjiwt_( &Asize,&Msize, &amax,&mmax, &mrastep, mracoef, data );
#else
   ++mrastep;
   for( int a =0; a <= amax; a++ ) {
      for(int m = (mmax-1); m >= 0; --m ) {
         if ( (mrastep & (1<<m)) != 0 ) { //right
            int tmp = mracoef[a*2*Msize+2*(m+1)+1];
            if ( (tmp >= 0) || ((tmp & 0x00000001) == 0) ) {
               mracoef[a*2*Msize+2*m]=mracoef[a*2*Msize+2*(m+1)] - tmp/2;
            } else {
               mracoef[a*2*Msize+2*m]=mracoef[a*2*Msize+2*(m+1)] - (tmp/2-1);
            }
         } else { //left
            int tmp = mracoef[a*2*Msize+2*(m+1)+1] + 1;
            if ( (tmp >= 0) || ((tmp & 0x00000001) == 0) ) {
               mracoef[a*2*Msize+2*m]=mracoef[a*2*Msize+2*(m+1)] + tmp/2;
            } else {
               mracoef[a*2*Msize+2*m]=mracoef[a*2*Msize+2*(m+1)] + (tmp/2-1);
            }
         }
      }
      data[a] = mracoef[a*2*Msize];
   }
#endif
} 

---