AdaptiveSmooth.h

//-*-C++-*-
/***************************************************************************
 *
 *   Copyright (C) 2008 by Paul Demorest
 *   Licensed under the Academic Free License version 2.1
 *
 ***************************************************************************/
 
#ifndef __Pulsar_AdaptiveSmooth_h
#define __Pulsar_AdaptiveSmooth_h
 
#include "Pulsar/Transformation.h"
 
namespace Pulsar {
 
  class Profile;
 
  class AdaptiveSmooth : public Transformation<Profile> {
 
  public:
 
    AdaptiveSmooth ();
 
    virtual ~AdaptiveSmooth ();
 
    float get_turns () const;
 
    unsigned get_max_harm() const;
 
    void set_hold (bool);
 
    bool get_hold() const;
 
    void transform(Profile *);
 
    enum Method { Wiener, Sinc_MSE, Sinc_IC, VonMises };
 
    void set_method(Method m);
 
  protected:
 
    std::vector<float> filter;
 
    int max_harm;
 
    bool hold;
 
    float noise_band;
 
    Method method;
 
  private:
 
    void compute(const float *fprof, int nbin);
 
    void compute_wiener(const float *pspec, double sigma2, int nh);
 
    void compute_sinc_mse(const float *pspec, double sigma2, int nh);
 
    void compute_sinc_ic(const float *pspec, double sigma2, int nh);
 
    void compute_lpf(const float *pspec, double sigma2, int nh,
        float (*filter_func)(float,float));
 
    void compute_vonmises_lpf(const float *pspec, double sigma2, int nh);
 
  };
 
}
 
#endif