DCFmp_Rtmra_double_array.h

Go to the documentation of this file.

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

template<>
void DCFmp< Rtmra<double,array> >::encode( 
  Rtmra<double,array>& x, const char mcode
) 
{
   double tmp;
   bytepack bp;

   int n = 0;
   if (mcode_to_sw(mcode)) { 
      n = ( (x.mrastep()) & ((1<<x.mmax())-1) ) >> mcode_to_m(mcode); 
   } 

   clear();
   set( mcode, n );

   for( int a = 0; a <= x.amax(); ++a ) {
      if ( mcode_to_sw(mcode) == 0 ) { tmp = x.scal( a, mcode_to_m(mcode) ); }
      else { tmp = x.wav( a, mcode_to_m(mcode) ); }
      if ( fabs(tmp) >= _threshold ) {
         push_back( bp( a, _bplen_a ) );
         push_back( bp( tmp, _bplen_coef ) );
      }
   }

};

template<>
void DCFmp< Rtmra<double,array> >::replace( Rtmra<double,array>& x ) 
{
   int a=0;

   if ( !null() && !empty() ) {
      double tmp;
      bytepack bp;
      int next=0;
      int a1;
      for(int i=0; i < ((*_data).size() / ( _bplen_a + _bplen_coef )); ++i ) {
         bp( 0, _bplen_a );
         for( int j = 0; j < bp.size(); ++j ) { bp[j] = (*_data)[next++]; }
         a1=int(bp);

         bp( 0, _bplen_coef );
         for( int j = 0; j < bp.size(); ++j) { bp[j] = (*_data)[next++]; }
         tmp = double(bp);

         while ( a < a1 ) { 
            if ( sw() ) { x.wav( a, m(), 0 ); }
            else { x.scal( a, m(), 0 ); }
            ++a; 
         }
         if ( sw() ) { x.wav( a, m(), tmp ); }
         else { x.scal( a, m(), tmp ); }
         ++a;
      }
   }

   while ( a <= x.amax() ) {
      if ( sw() ) { x.wav( a, m(), 0 ); }
      else { x.scal( a, m(), 0 ); }
      ++a; 
   }
};

---