rtrsdetectle.cpp

Go to the documentation of this file.

// rtrsdetectle.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 "rtrsdetectle_.h"
#include "drtrsdetectle_.h"
#include "irtrsdetectle_.h"
#include "artrsdetectle_.h"
#include "dartrsdetectle_.h"
#include "iartrsdetectle_.h"
#endif


void rtrsdetectle( 
          const int& Msize,
          const int& mmax,
          const int& mrastep,
          const float* mracoef,
          const float* mracoefrs,
          const float* thresh,
                int* hit
) 
{ 
#ifdef STORM_USEF
   rtrsdetectle_( 
     &Msize, 
     &mmax, 
     &mrastep, mracoef, mracoefrs,
     thresh, hit
   );
#else
   for( int m = 0; m <= mmax; m++ ) {
      if ( abs( mracoef[ 2*m + 1 ] + mracoefrs[ 2*m + 1 ] ) <= thresh[m] ) { 
         hit[m] = 1; 
      }
      else { 
         hit[m] = 0; 
      }
   }
#endif
}


void rtrsdetectle( 
          const int& Msize,
          const int& mmax,
          const int& mrastep,
          const double* mracoef,
          const double* mracoefrs,
          const double* thresh,
                int* hit
) 
{ 
#ifdef STORM_USEF
   drtrsdetectle_( 
     &Msize, 
     &mmax, 
     &mrastep, mracoef, mracoefrs,
     thresh, hit
   );
#else
   for( int m = 0; m <= mmax; m++ ) {
      if ( abs( mracoef[ 2*m + 1 ] + mracoefrs[ 2*m + 1 ]) <= thresh[m] ) { 
         hit[m] = 1; 
      }
      else { 
         hit[m] = 0; 
      }
   }
#endif
} 


void rtrsdetectle( 
          const int& Msize,
          const int& mmax,
          const int& mrastep,
          const int* mracoef,
          const int* mracoefrs,
          const int* thresh,
                int* hit
) 
{ 
#ifdef STORM_USEF
   irtrsdetectle_( 
     &Msize, 
     &mmax, 
     &mrastep, mracoef, mracoefrs,
     thresh, hit
   );
#else
   for( int m = 0; m <= mmax; m++ ) {
      if ( abs( mracoef[ 2*m + 1 ] + mracoefrs[ 2*m + 1 ]) <= thresh[m] ) { 
         hit[m] = 1; 
      }
      else { 
         hit[m] = 0; 
      }
   }
#endif
} 


void rtrsdetectle(
          const int& Asize,
          const int& Msize,
          const int& amax,
          const int& mmax,
          const int& mrastep,
          const float* mracoef,
          const float* mracoefrs,
          const float* thresh,
                int* hit
)
{
#ifdef STORM_USEF
   artrsdetectle_(
     &Asize, &Msize,
     &amax, &mmax,
     &mrastep, mracoef, mracoefrs,
     thresh, hit
   );
#else
   for( int m = 0; m <= mmax; m++ ) {   
      hit[m] = 0;
      for( int a = 0; a <= amax; a++ ) {
         if ( abs( mracoef[ Msize*2*a + 2*m + 1 ]
           +mracoefrs[ Msize*2*a + 2*m + 1 ]) <= thresh[m] ) { hit[m]++; }
      }
   }
#endif
}


void rtrsdetectle(
          const int& Asize,
          const int& Msize,
          const int& amax,
          const int& mmax,
          const int& mrastep,
          const double* mracoef,
          const double* mracoefrs,
          const double* thresh,
                int* hit
)
{
#ifdef STORM_USEF
   dartrsdetectle_(
     &Asize, &Msize,
     &amax, &mmax,
     &mrastep, mracoef, mracoefrs,
     thresh, hit
   );
#else
   for( int m = 0; m <= mmax; m++ ) {   
      hit[m] = 0;
      for( int a = 0; a <= amax; a++ ) {
         if ( abs( mracoef[ Msize*2*a + 2*m + 1 ]
           + mracoefrs[ Msize*2*a + 2*m + 1 ] ) <= thresh[m] ) { hit[m]++; }
      }
   }
#endif
}

void rtrsdetectle(
          const int& Asize,
          const int& Msize,
          const int& amax,
          const int& mmax,
          const int& mrastep,
          const int* mracoef,
          const int* mracoefrs,
          const int* thresh,
                int* hit
)
{
#ifdef STORM_USEF
   iartrsdetectle_(
     &Asize, &Msize,
     &amax, &mmax,
     &mrastep, mracoef, mracoefrs,
     thresh, hit
   );
#else
   for( int m = 0; m <= mmax; m++ ) {   
      hit[m] = 0;
      for( int a = 0; a <= amax; a++ ) {
         if ( abs( mracoef[ Msize*2*a + 2*m + 1 ]
           + mracoefrs[ Msize*2*a + 2*m + 1 ] ) <= thresh[m] ) { hit[m]++; }
      }
   }
#endif
}

---