psc_accumulator_1st.hpp

/* psc_accumulator_1st.hpp */
/* $Id$ */
#ifndef __psc_accumulator_1st_hpp
#define __psc_accumulator_1st_hpp

#include <s2proc/psc_accumulator.hpp>

namespace psc
{
  class Accumulator1st : public Accumulator
  {
  public:
    
    Accumulator1st(const dba::Seismic &db, const Parameters &params);
    
    virtual ~Accumulator1st();
    
    bool operator()(const WaveField &second_order_row,
            WaveField &accumulated_row) const;
    
  protected:
    
    struct RowsTimeShift
    {
      int id1, id2, id3, id4;
      Trace::TimeT dt;
      Trace::AmplT w;
      bool used;
      
      static bool id1_comp(const RowsTimeShift &row1, const RowsTimeShift &row2)
      {return row1.id1 < row2.id1;}
    };
    
    Trace::Header m_def_hdr;
        
    bool accumulate(const int stride, 
            std::vector<RowsTimeShift> &rows_time_shifts,
            WaveField &accumulated_row,
            bool *all_are_used_flag) const;
    
    bool accumulate_row(const int row_id,
            const unsigned int stride,
            const Trace::TimeT time_shift,
            WaveField &accumulated_row) const;
    
    size_t get_best_row_id(const std::vector<RowsTimeShift> &rows_time_shifts,
               bool *all_are_used_flag) const;
  };
}; // namespace psc

#endif /* psc_accumulator_1st.hpp */