prd_service.hpp

/* prd_service.hpp */
/* $Id: prd_service.hpp 21566 2012-04-18 12:25:05Z sasha $ */
#ifndef __prd_service_hpp
#define __prd_service_hpp

#include <mix/Types.h>       // for "Fortran" macros

#include <mth/lev.h>
#include <mth/fft.h>

#include <s2proc/prd_service_f.hpp>

#define MIN_ACF_ZERO_SAMPLE_VALUE 1.e-20

#include <stdio.h>
#include <mix/wenv.h>
#include <stdarg.h>
#include <vector>

typedef struct {
  unsigned int uTeTa;
  float fAlfa;
  unsigned int opLen;
  unsigned int fOpDomain; /* 0 - frequency domain, 1 - time domain */
} OperatorParams;

bool getTrACF(std::vector<float>* wBody           /* current trace buffer */, 
          float* pfACfunction                /* current trace ACF function (OUTPUT) */);

bool PRDconvStep(float* pfACF                    /* operator spectrum */, 
         std::vector<float> *curTrace     /* current trace */, 
         int uSamplesFFT,
         float* pfTrace);

int TrFindGain(std::vector<float> *curTrace, float& Gain);

bool PredOpForm(unsigned int uCount,  
        OperatorParams *opPar, 
        float* pfACF, 
        std::vector<double> &operLevOut);

bool PredErrOpForm(const unsigned int uCount, 
           OperatorParams *opPar, 
           float* pfACF, 
           std::vector<double> &oper);

bool PredErrOpSpecForm(const unsigned int uCount,  
               OperatorParams *opPar           /* operator parameters */, 
               float* pfACF                   /* operator spectrum (filtered) (OUTPUT) */, 
               unsigned int uCurrentWinSize);

#endif /* prd_service.hpp */