Rescale.h

//-*-C++-*-
/***************************************************************************
 *
 *   Copyright (C) 2006-2024 by Willem van Straten
 *   Licensed under the Academic Free License version 2.1
 *
 ***************************************************************************/
 
// dspsr/Signal/General/dsp/Rescale.h
 
#ifndef __dsp_Rescale_h
#define __dsp_Rescale_h
 
#include "dsp/Transformation.h"
#include "dsp/TimeSeries.h"
 
#include <vector>
 
namespace dsp
{
  class Rescale : public Transformation<TimeSeries,TimeSeries>
  {
 
  public:
 
    Rescale ();
 
    ~Rescale () = default;
 
    void prepare ();
 
    void transformation ();
 
    void set_interval_seconds (double seconds);
 
    void set_interval_samples (uint64_t samples);
 
    void set_exact (bool exact);
 
    void set_constant (bool constant_offset_scale);
 
    void set_decay (float decay_constant);
 
    void set_output_after_interval (bool output_after_interval);
 
    void set_output_time_total (bool output_time_total);
 
    MJD get_update_epoch () const;
 
    const float* get_offset (unsigned ipol) const;
 
    const float* get_scale (unsigned ipol) const;
 
    const double* get_mean (unsigned ipol) const;
 
    const double* get_variance (unsigned ipol) const;
 
    uint64_t get_nsample () const;
 
    const float* get_time (unsigned ipol) const;
 
  private:
 
    std::vector< std::vector<double> > freq_total;
 
    std::vector< std::vector<double> > freq_totalsq;
 
    std::vector< std::vector<float> > time_total;
 
    std::vector< std::vector<float> > scale;
 
    std::vector< std::vector<float> > offset;
 
    std::vector< std::vector<float> > decay_offset;
 
    bool exact{false};
 
    bool output_time_total{false};
 
    bool output_after_interval{false};
 
    double interval_seconds{0.0};
 
    uint64_t interval_samples{0};
 
    float decay_constant{1e4};
 
    bool do_decay{false};
 
    uint64_t nsample{0};
 
    uint64_t isample{0};
 
    MJD update_epoch{};
 
    bool constant_offset_scale{false};
 
    bool first_integration{false};
 
    void init ();
 
    void compute_various (bool first_call = false, bool first_integration = false);
  };
 
} // namespace dsp
 
#endif // __dsp_Rescale_h