fitthreat.h

#ifndef FITTHREAT_H
#define FITTHREAT_H
#include <gsl/gsl_vector.h>
#include <gsl/gsl_complex.h>
#include <gsl/gsl_complex_math.h>
#include <gsl/gsl_eigen.h>
#include <gsl/gsl_sort_vector.h>
#include <gsl/gsl_math.h>
#include <gsl/gsl_matrix.h>
#include <gsl/gsl_linalg.h>
#include <gsl/gsl_blas.h>
#include <gsl/gsl_multifit_nlin.h>
#include <gsl/gsl_multimin.h>
#include <gsl/gsl_multiroots.h>
#include <gsl/gsl_errno.h>
#include <QProgressDialog>
#include <QMessageBox>
#include <QThread>
#include <QMutex>
#include <QTime>
#include "array.h"
#include "math.h"

typedef struct fitcontent{
    double (*chi2)(void*);
    double (*runchi2)(void*,int);
    void* instance;
    int cpara;
    int* IsChanged;
    double* param;
    int runi;
}fitcontent;
enum threattyp {fitseriell,fitglobal,errorseriell,errorglobal,notvalid};

class CorrelationEle{
    public:
        CorrelationEle(int ia, int ib, double vala,
                       double valb, double valab,
                       double xa, double xb);
        ~CorrelationEle();
        bool operator<(const CorrelationEle o2)const;
        int a,b;
        double va,vb;
        double correlation;
};


class fitthreat: public QThread{
Q_OBJECT
public:
    void run();
    void exit();
    fitthreat(double* para, int* changed, double* errors, int cpara,
              double (*f)(void*), void* instance, int MaxIter,
              double conv, bool forceSimplex);
    fitthreat(double* para, int runi, double* errors, int cpara,
              double (*f)(void*,int), void* instance, int MaxIter,
              double conv);
    fitthreat();
    ~fitthreat();
    fitthreat(double* para, int* changed, double *errors, int cpara,
              double (*chi2)(void*), void* instance, bool correlated);
    fitthreat(double* para, int runi, double *errors, int cpara,
              double (*chi2)(void*,int), void* instance);
public slots:
    void stop_calculations();

private:
    void fit_seriell();
    void fit_global();
    void fit_global_simplex();
    void errors_covariant();
    int correlation_possible(int xi, int yi);
    void Hesse_to_covariant_errors(gsl_matrix *H, double* errors,
                                   double* vals, int cpara);
    QString correlations(gsl_matrix *H, int cpara, double* vals);
    double fxx(double* v, double v0, double* ss, double chi0);
    double fx(double* v, double v0, double* ss, double chi0);
    double fxy(double* vx, double* vy,
               double v0x, double v0y,
               double fxx, double fyy,
               double ssx, double ssy, double chi0);
    double* para;
    int runi;
    double* errors;
    int* changed;
    int cpara;
    bool correlated;
    double (*f_global)(void*);
    double (*f_seriell)(void*,int);
    double f();
    double (*chi2)(void*);
    void* instance;
    int MaxIter;
    threattyp fitmodel;
    int proceed_calc;
    double conv;
    bool usesimplex;
signals:
    void refresh_progress(int ID, int citer,double chi2);
    void finished(int ID);
    void error_message(QString message);

};

#endif // FITTHREAT_H