MHO_IonosphericFringeFitterOpenMP.hh

Go to the documentation of this file.

#ifndef MHO_IonosphericFringeFitterOpenMP_HH__
#define MHO_IonosphericFringeFitterOpenMP_HH__
 
#include <memory>
#include <vector>
 
#include "MHO_BasicFringeFitter.hh"
#include "MHO_IonosphericPhaseCorrection.hh"
#include "MHO_Tokenizer.hh"
 
#ifdef _OPENMP
    #include <omp.h>
#endif
 
namespace hops
{
 
struct IonLoopResult
{
        double famp = 0.0;
        double ion_diff = 0.0;
};
 
class MHO_IonosphericFringeFitterOpenMP: public MHO_BasicFringeFitter
{
 
    public:
        MHO_IonosphericFringeFitterOpenMP(MHO_FringeData* data);
        virtual ~MHO_IonosphericFringeFitterOpenMP();
 
        virtual void Run() override;
 
        void Finalize() override;
 
        virtual void Accept(MHO_FringeFitterVisitor* visitor) override { visitor->Visit(this); };
 
    protected:
        int rjc_ion_search();
 
        void sort_tecs(int nion, std::vector< std::vector< double > >& dtec);
 
        int ion_search_smooth();
 
        void smoother(double* f, double* g, double* tec_step, int* npts);
 
        double calculate_approx_snr();
 
        double fInitialSBWin[2]; //save the initial SBD window
        int ion_npts;
 
        //per-thread infrastructure for parallel ion search
        int fNIonThreads;
        visibility_type fVisRef; //clean reference copy of vis_data (before any ion correction)
 
        std::vector< visibility_type > fPerThreadVis;                  //per-thread working vis copies
        std::vector< visibility_type > fPerThreadSBD;                  //per-thread SBD output buffers
        std::vector< std::unique_ptr< MHO_NormFX > > fPerThreadNormFX; //per-thread norm fx operator
 
        std::vector< std::unique_ptr< MHO_MBDelaySearch > > fPerThreadMBDSearch;
        std::vector< std::unique_ptr< MHO_InterpolateFringePeakOptimized > > fPerThreadPeakInterp;
        std::vector< std::unique_ptr< MHO_IonosphericPhaseCorrection > > fPerThreadIono;
 
        std::vector< IonLoopResult > fIonLoopResults; //one slot per ionloop index
 
    private:
        void initialize_ion_threads();
};
 
} // namespace hops
 
#endif