rtrsdetectlt.cpp

Go to the documentation of this file.

// rtrsdetectlt.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 "rtrsdetectlt_.h"
#include "drtrsdetectlt_.h"
#include "irtrsdetectlt_.h"
#include "artrsdetectlt_.h"
#include "dartrsdetectlt_.h"
#include "iartrsdetectlt_.h"
#endif


void rtrsdetectlt( 
          const int& Msize,
          const int& mmax,
          const int& mrastep,
          const float* mracoef,
          const float* mracoefrs,
          const float* thresh,
                int* hit
) 
{ 
#ifdef STORM_USEF
   rtrsdetectlt_( 
     &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 rtrsdetectlt( 
          const int& Msize,
          const int& mmax,
          const int& mrastep,
          const double* mracoef,
          const double* mracoefrs,
          const double* thresh,
                int* hit
) 
{ 
#ifdef STORM_USEF
   drtrsdetectlt_( 
     &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 rtrsdetectlt( 
          const int& Msize,
          const int& mmax,
          const int& mrastep,
          const int* mracoef,
          const int* mracoefrs,
          const int* thresh,
                int* hit
) 
{ 
#ifdef STORM_USEF
   irtrsdetectlt_( 
     &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 rtrsdetectlt(
          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
   artrsdetectlt_(
     &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 rtrsdetectlt(
          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
   dartrsdetectlt_(
     &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 rtrsdetectlt(
          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
   iartrsdetectlt_(
     &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
}

 

---