splineitem.h

/***************************************************************************
                          splineitem.h  -  description
                             -------------------
    begin                : Wed Feb 27 2002
    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 SPLINEITEM_H
#define SPLINEITEM_H
 
#include "funitem.h"
 
class Q_DECL_EXPORT SplineItem : public ScaledItem
{
public:
  enum {
    Type = UserType + Spline
  };
  SplineItem(int degree = 3, int derivOrder = 0, double xLow = 0,
             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);
  SplineItem(const SplineItem& a);
  SplineItem(const KConfigGroup* plo, unsigned color = 0,
             Symbol symbol = SolidLine, QGraphicsItem* parent = nullptr);
  virtual ~SplineItem();
  virtual QRectF boundingRect() const override;
  int calcTable(bool logX) const;
  int calcTable(const double* x, double* values, int n = 1, int derivOrder = 0,
                 double xLow = 0) const;
  void copySpline(const SplineItem* a);
  inline int degree() const { return m_degree; };
  inline double deltaX() const { return m_deltaX; };
  inline int derivOrder() const { return m_derivOrder; };
  void exponent(int& expX, int& expY, double& fx, double& fy, bool logX) const;
  void exportTable(QTextStream& stream, char separator, char format,
                   int precision) const;
  int fit(const ArrayItem* ai, double xMin, double xMax,
          double* chiSquare = nullptr, int* numKnots = nullptr,
          double smooth = 0, FunItem::DataErrorStruct* err = nullptr,
          double* avgErr = nullptr);
  int fit(const double* x, const double* y, const double* errors, int n,
          double xMin, double xMax, double* chiSquare = nullptr,
          int* numKnots = nullptr, double smooth = 0,
          FunItem::DataErrorStruct* err = nullptr, double* avgErr = nullptr);
  int interpolation(const double* x, const double *y, int n, double xMin,
                    double xMax);
  void limits(double& xMin, double& xMax, double& yMin, double& yMax,
              bool logX) const;
  inline int numKnots() const { return m_numKnots; };
  virtual void paint(QPainter* painter, const QStyleOptionGraphicsItem* option,
                     QWidget* widget = nullptr) override;
  int roots(QVector<double>& zero, int derivOrder = 0,
            double offset = 0) const;
  void setDashOffset();
  inline void setDegree(int degree) { m_degree = degree; };
  inline void setDeltaX(double deltaX) { m_deltaX = deltaX; };
  inline void setDerivOrder(int derivOrder) { m_derivOrder = derivOrder; };
  inline void setXLow(double xLow) { m_xLow = xLow; };
  inline void setXMax(double xMax) { m_xMax = xMax; };
  inline void setXMin(double xMin) { m_xMin = xMin; };
  virtual QPainterPath shape() const override;
  static int sortArrays(const double* x0, const double* y0,
                        const double* errors0, int& n,
                        FunItem::DataErrorStruct* err, double* x, double* y,
                        double* errors);
  virtual int type() const override { return Type; };
  void writePlo(KConfigGroup* plo, char format, int precision) const;
  inline double xLow() const { return m_xLow; };
  inline double xMax() const { return m_xMax; };
  inline double xMin() const { return m_xMin; };
  double* xValues() const;
  double* yValues() const;
 
protected:
  void paintCond(QPainter* painter) const;
  mutable bool logXo;
  int m_degree, m_numKnots, m_derivOrder;
  mutable int degreeo, numKnotso, derivOrdero;
  double m_xMin, m_xMax, m_deltaX, m_xLow;
  mutable double xMino, xMaxo, deltaXo, xLowo;
  mutable QVector<double> m_xValues, m_yValues, t, c;
};
 
#endif