ZeroBufferFit.h

#ifndef ZeroBufferFit_h
#define ZeroBufferFit_h
 
/* SPDX-License-Identifier: CC0-1.0 */
 
#include <string>
#include <vector>
 
#include "HiCalTypes.h"
#include "HiCalUtil.h"
#include "HiCalConf.h"
#include "NonLinearLSQ.h"
#include "Module.h"
#include "IException.h"
 
namespace Isis {
 
  class ZeroBufferFit : public NonLinearLSQ, public Module {
 
    public:
    ZeroBufferFit(const HiCalConf &conf);
    virtual ~ZeroBufferFit() { }
 
    void setBin(int bin) { _timet.setBin(bin); }
 
    void setLineTime(double ltime) { _timet.setLineTime(ltime); }
 
     inline int size() const { return (_data.dim()); }
 
    int nSize()  const { return (_b2.dim()); }
 
    int nParms() const { return (4); }
 
    void setabsErr(double absError) { _absErr = absError; }
    void setrelErr(double relError) { _relErr = relError; }
    double absErr() const { return (_absErr); }
    double relErr() const { return (_relErr); }
 
    HiVector Solve(const HiVector &d);
    NLVector guess();
    int checkIteration(const int Iter, const NLVector &fitcoefs,
                       const NLVector &uncerts, double cplxconj,
                       int Istatus);
 
    NLVector f_x(const NLVector &a);
    NLMatrix df_x(const NLVector &a);
 
    double Chisq() const { return (_chisq); }
    int    DoF()   const { return (nSize() - nParms()); }
 
    HiVector Yfit() const;
    HiVector Normalize(const HiVector &v);
 
  private:
     HiLineTimeEqn _timet;    //  This is the X data set
     HiVector      _data;     //  Typically will be the HiRISE buffer data
     HiVector      _b2;       //  Data buffer used in fitting
     double        _absErr;   //  Absolute error convergence test
     double        _relErr;   //  Relative error convergence test
     double        _maxLog;   //  Maximum log value to constrain
     int           _badLines; //  Exclude lines at end of buffer
     int           _sWidth;   //  Width of guestimate filter
     int           _sIters;   //  Filter interations
     bool          _skipFit;  //  Skip fitting and pass input through
     bool          _useLinFit; // Use linear fit on failure of LM, else Zf
     int           _minLines; //  Minimum number of lines to fit (Default: 100)
     HiVector      _cc;       //  Parameter of 2-D fit
     HiVector      _guess;    //  Initial guestimate of solutions
     HiVector      _coefs;    //  Coefficients of solution
     HiVector      _uncert;   //  Uncertanties
     double        _chisq;    //  ChiSq of NonLinear equation
 
     HiVector poly_fit(const HiVector &d, const double line0 = 0.0) const;
     virtual void printOn(std::ostream &o) const;
     int goodLines(const HiVector &d) const { return (d.dim() - _badLines); }
     bool gotGoodLines(const HiVector &d) const {
       return (goodLines(d) >= _minLines);
     }
  };
 
}     // namespace Isis
#endif