rtrsdetectgt.cpp

Go to the documentation of this file.

// rtrsdetectgt.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 "rtrsdetectgt_.h"
#include "drtrsdetectgt_.h"
#include "irtrsdetectgt_.h"
#include "artrsdetectgt_.h"
#include "dartrsdetectgt_.h"
#include "iartrsdetectgt_.h"
#endif


void rtrsdetectgt( 
          const int& Msize,
          const int& mmax,
          const int& mrastep,
          const float* mracoef,
          const float* mracoefrs,
          const float* thresh,
                int* hit
) 
{ 
#ifdef STORM_USEF
   rtrsdetectgt_( 
     &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 rtrsdetectgt( 
          const int& Msize,
          const int& mmax,
          const int& mrastep,
          const double* mracoef,
          const double* mracoefrs,
          const double* thresh,
                int* hit
) 
{ 
#ifdef STORM_USEF
   drtrsdetectgt_( 
     &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 rtrsdetectgt( 
          const int& Msize,
          const int& mmax,
          const int& mrastep,
          const int* mracoef,
          const int* mracoefrs,
          const int* thresh,
                int* hit
) 
{ 
#ifdef STORM_USEF
   irtrsdetectgt_( 
     &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 rtrsdetectgt(
          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
   artrsdetectgt_(
     &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 rtrsdetectgt(
          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
   dartrsdetectgt_(
     &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 rtrsdetectgt(
          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
   iartrsdetectgt_(
     &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
}

 

---