rtrsmrawavdrs.cpp

Go to the documentation of this file.

// rtrsmrawavdrs.cpp

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

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

#ifdef STORM_USEF
#include "rtrsmrawavdrs_.h"
#include "drtrsmrawavdrs_.h"
#include "irtrsmrawavdrs_.h"
#include "artrsmrawavdrs1_.h"
#include "dartrsmrawavdrs1_.h"
#include "iartrsmrawavdrs1_.h"
#endif

#include "stormdef.h"
#define INVSQRT2 STORMDEF_MATHCONST_INVSQRT2

void rtrsmrawavdrs( 
          const int& Msize,
          const int& mmax,
          const int& mrastep,
          const float* mracoef,
          const float* mracoefrs,
          const int& m,
                float& wavdrs
) 
{ 
#ifdef STORM_USEF
   rtrsmrawavdrs_( 
     &Msize, 
     &mmax, 
     &mrastep, mracoef, mracoefrs,
     &m, &wavdrs
   );
#else
//   if ( ( m >= 0 ) && ( m <= mmax ) ) {
      wavdrs = mracoefrs[ 2*m +1 ] * pow(INVSQRT2,m);
//   } else {
//       wavdrs = 1.0/0;
//   }
#endif
}


void rtrsmrawavdrs( 
          const int& Msize,
          const int& mmax,
          const int& mrastep,
          const double* mracoef,
          const double* mracoefrs,
          const int& m,
                double& wavdrs
) 
{ 
#ifdef STORM_USEF
   drtrsmrawavdrs_( 
     &Msize, 
     &mmax, 
     &mrastep, mracoef, mracoefrs,
     &m, &wavdrs
   );
#else
//   if ( ( m >= 0 ) && ( m <= mmax ) ) {
      wavdrs = mracoefrs[ 2*m + 1 ] * pow(INVSQRT2,m);
//   } else {
//      wavdrs = 1.0/0;
//   }
#endif
} 


void rtrsmrawavdrs( 
          const int& Msize,
          const int& mmax,
          const int& mrastep,
          const int* mracoef,
          const int* mracoefrs,
          const int& m,
                int& wavdrs
) 
{ 
#ifdef STORM_USEF
   irtrsmrawavdrs_( 
     &Msize, 
     &mmax, 
     &mrastep, mracoef, mracoefrs,
     &m, &wavdrs
   );
#else
//   if ( ( m >= 0 ) && ( m <= mmax ) ) {
      int tmp = mracoefrs[ 2*m + 1 ];
      if (tmp >= 0) {
         wavdrs = tmp >> m;
      } else {
         if ((tmp & 0x00000001) != 0) {
            wavdrs = - ( (-tmp) >> m);
         } else {
            wavdrs = - ( (-tmp) >> m) - 1;
         }
      }
//   } else {
//      wavdrs = 1.0/0;
//   }
#endif
} 


void rtrsmrawavdrs(
          const int& Asize,
          const int& Msize,
          const int& amax,
          const int& mmax,
          const int& mrastep,
          const float* mracoef,
          const float* mracoefrs,
          const int& a,
          const int& m,
                float& wavdrs
)
{
#ifdef STORM_USEF
   artrsmrawavdrs1_(
     &Asize, &Msize,
     &amax, &mmax,
     &mrastep, mracoef, mracoefrs,
     &a, &m, &wavdrs
   );
#else
//   if ( ( a >= 0 ) && ( a <= amax ) && ( m >= 0 ) && ( m <= mmax ) ) {
      wavdrs = mracoefrs[ Msize*2*a + 2*m + 1 ] * pow(INVSQRT2,m);
//   } else {
//      wavdrs = 1.0/0;
//   }
#endif
}


void rtrsmrawavdrs(
          const int& Asize,
          const int& Msize,
          const int& amax,
          const int& mmax,
          const int& mrastep,
          const double* mracoef,
          const double* mracoefrs,
          const int& a,
          const int& m,
                double& wavdrs
)
{
#ifdef STORM_USEF
   dartrsmrawavdrs1_(
     &Asize, &Msize,
     &amax, &mmax,
     &mrastep, mracoef, mracoefrs,
     &a, &m, &wavdrs
   );
#else
//   if ( ( a >= 0 ) && ( a <= amax ) && ( m >= 0 ) && ( m <= mmax ) ) {
      wavdrs = mracoefrs[ Msize*2*a + 2*m + 1 ] * pow(INVSQRT2,m);
//   } else {
//      wavdrs = 1.0/0;
//   }
#endif
}


void rtrsmrawavdrs(
          const int& Asize,
          const int& Msize,
          const int& amax,
          const int& mmax,
          const int& mrastep,
          const int* mracoef,
          const int* mracoefrs,
          const int& a,
          const int& m,
                int& wavdrs
)
{
#ifdef STORM_USEF
   iartrsmrawavdrs1_(
     &Asize, &Msize,
     &amax, &mmax,
     &mrastep, mracoef, mracoefrs,
     &a, &m, &wavdrs
   );
#else
//   if ( ( a >= 0 ) && ( a <= amax ) && ( m >= 0 ) && ( m <= mmax ) ) {
      int tmp = mracoefrs[ Msize*2*a + 2*m + 1 ];
      if (tmp >= 0) {
         wavdrs = tmp >> m;
      } else {
         if ((tmp & 0x00000001) != 0) {
            wavdrs = - ( (-tmp) >> m);
         } else {
            wavdrs = - ( (-tmp) >> m) - 1;
         }
      }
//   } else {
//      wavdrs = 1/0;
//   }
#endif
}

---