funitem.h

/***************************************************************************
                          funitem.h  -  description
                             -------------------
 
    This file is a part of kpl - a program for graphical presentation of
    data sets and functions.
 
    begin                : Sun Aug 29 1999
    copyright            : (C) 2024 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 FUNITEM_H
#define FUNITEM_H
 
#include <QVector>
#include <QUrl>
#include "scaleditem.h"
 
class QLibrary;
class QTextStream;
class ArrayItem;
class LMFit;
 
class Q_DECL_EXPORT FunItem : public ScaledItem
{
public:
  static const int NumParMax = 20;
  enum {
    Type = UserType + Function
  };
  struct DataErrorStruct {
    enum ErrorModelArgument {
      ColX, 
      ColY  
    };
    DataErrorStruct() : errModArg(ColY) {};
    bool fitErrCol;
    ErrorModelArgument errModArg;
    double parErrMod[NumParMax];
    QString errModName;
    QString errModPath;
  };
  FunItem();
  FunItem(const FunItem& a);
  FunItem(const KConfigGroup* plo, const QUrl& urlPlo, unsigned color = 0,
          Symbol symbol = SolidLine, QGraphicsItem* parent = nullptr);
  FunItem(const QUrl& url, const QString& name, unsigned color = 0,
          Symbol symbol = SolidLine, double fx = 1, double fy = 1,
          QGraphicsItem* parent = nullptr, double xMin = 0, double xMax = 0,
          double deltaX = 0, bool visible = true,
          FillStyle fillStyle = NoBrush, double xOffset = 0,
          double yOffset = 0, double relSize = 1);
  virtual ~FunItem();
  virtual QRectF boundingRect() const override;
  static void calcTable(const double* x, double* values,
                        double (*func)(double, const double*),
                        const double* par, int n);
  int calcTable() const;
  int calcTable(bool logX) const;
  inline double deltaT() const { return m_deltaT; };
  static bool errorStructChanged(const DataErrorStruct& e1,
                                 const DataErrorStruct& e2);
  void exponent(int& expX, int& expY, double& fx, double& fy, bool logX) const;
  void exportTable(QTextStream& stream, char separator, char format,
                   int precision) const;
  static int fit(const QList<const ArrayItem*>* aList, QList<FunItem*>* fList,
                 const bool* bFit, double* par,
                 QVector<double>* resVector = nullptr,
                 QVector<double>* sigma = nullptr,
                 FunItem::DataErrorStruct* err = nullptr, bool nonLin = true,
                 double tol = 1e-14, int itMax = 100, LMFit* lm = nullptr,
                 double* chiSquare = nullptr, double* corr = nullptr,
                 double* parErr = nullptr, double* avgErr = nullptr,
                 QObject* parent = nullptr);
  double (*funcY)(double, const double*);
  inline const QLibrary* libY() const { return m_libY; };
  void limits(double& xMin, double& xMax, double& yMin, double& yMax,
              bool logX) const;
  inline QString nameY() const { return m_nameY; };
  virtual void paint(QPainter* painter, const QStyleOptionGraphicsItem* option,
                     QWidget* widget = nullptr) override;
  static bool parChanged(const double* p1, const double* p2);
  static QString parToString(const double* par, char format, int precision);
  inline double* parY() { return m_parY; };
  void setDashOffset();
  inline void setDeltaT(double deltaT) { m_deltaT = deltaT; };
  bool setFuncY(const QUrl& url, const QString& name);
  virtual void setPar(int i, double value, bool yFunc = true);
  inline void setTMax(double xMax) { m_tMax = xMax; };
  inline void setTMin(double xMin) { m_tMin = xMin; };
  virtual QPainterPath shape() const override;
  static void stringToPar(const QString& s, double* par);
  inline double tMax() const { return m_tMax; };
  inline double tMin() const { return m_tMin; };
  const double* tValues() const;
  virtual int type() const override { return Type; };
  inline const QUrl& urlY() const { return m_urlY; };
  void writePlo(KConfigGroup* plo, const QUrl& urlPlo, bool absolute,
                char format, int precision) const;
  double* yValues() const;
 
protected:
  void init();
  void paintCond(QPainter* painter) const;
  void unload();
  mutable bool logXo;
  double m_tMin, m_tMax, m_deltaT;
  mutable double tMino, tMaxo, deltaTo;
  mutable double (*funcYo)(double, const double*);
  double m_parY[NumParMax];
  mutable double parYo[NumParMax];
  mutable QVector<double> m_tValues, m_yValues;
  QString m_nameY;
  QUrl m_urlY;
  QLibrary* m_libY;
};
 
#endif