lmfit.h

/***************************************************************************
                          lmfit.h  -  description
                             -------------------
    begin                : Sat May 22 1999
    copyright            : (C) 2017 by Werner Stille
    email                : kpl@werner-stille.de
 ***************************************************************************/

/***************************************************************************
 *                                                                         *
 *   This program is free software; you can redistribute it and/or modify  *
 *   it under the terms of the GNU General Public License as published by  *
 *   the Free Software Foundation; either version 2 of the License, or     *
 *   (at your option) any later version.                                   *
 *                                                                         *
 ***************************************************************************/

#ifndef LMDIF_H
#define LMDIF_H

#include <QObject>

class QString;
class ArrayItem;
class FunItem;

class Q_DECL_EXPORT LMFit : public QObject
{
  Q_OBJECT
public:
  LMFit(const QList<const ArrayItem*>* aList, QList<FunItem*>* fList,
        const bool* bFit, double* par, QObject* parent = 0,
        const double* sigma = 0);
  ~LMFit();
  static double enorm(int n, const double x[]);
  static double igamc(double a, double x);
  bool linFit(int n, int nPar, double par[], double resVector[], double fJac[],
              int iPvt[]);
  void lmdif(int n, int nPar, double par[], double resVector[], double tol,
             int itMax, int* info, double fJac[], int iPvt[]);
  static int minv(double a[], double x[], int n);
  inline void setSigma(const double* sigma) { m_sigma = sigma; };
  inline void setUserBreak(bool userBreak = true) { m_userBreak = userBreak; };
  inline bool userBreak() const { return m_userBreak; };

protected:
  void fcn(int n, int np, const double p[], double fvec[], int* iflag);
  void fdjac2(int m, int n, double x[], const double fvec[], double fjac[],
              int ldfjac, int* iflag, double epsfcn, double wa[]);
  void lmdif(int m, int np, double px[], double fvec[], double ftol,
             double ptol, double gtol, int maxfev, double epsfcn,
             double diag[], int mode, double factor, int nprint,
             int* info, int* nfev, double fjac[], int ldfjac,
             int ipvt[], double qtf[], double wa1[], double wa2[],
             double wa3[], double wa4[]);
  bool m_userBreak;
  const bool* m_bFit;
  double* m_par;
  const double* m_sigma;
  const QList<const ArrayItem*>* m_aList;
  QList<FunItem*>* m_fList;

signals:
  void updateMessage(const QString& message);
};

#endif