---

qscontour.h

/***************************************************************************
                                qscontour.h
                             -------------------                                         
    begin                : 01-January-2000
    copyright            : (C) 2000 by Kamil Dobkowski                         
    email                : kamildobk@poczta.onet.pl                                     
 ***************************************************************************/

/***************************************************************************
 *                                                                         *
 *   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 QSCONTOUR_H
#define QSCONTOUR_H

#ifdef HAVE_CONFIG_H
#include <config.h>
#endif
#include"qsaxes.h"
#include"qsplot.h"
#include<math.h>


class QSContour : public QSPlot2D
  {
   Q_OBJECT
    Q_PROPERTY( bool contourFills READ contourFills WRITE setContourFills )
    Q_PROPERTY( bool contourLines READ contourLines WRITE setContourLines )
    Q_PROPERTY( bool contourLabels READ contourLabels WRITE setContourLabels )
    Q_PROPERTY( double labelSpacing READ labelSpacing WRITE setLabelSpacing )

   public:
    QSContour(QSAxes* parent, const char * name=0);
    virtual ~QSContour();

    void setContourFills( bool visible );
    void setContourLines( bool visible );
    void setContourLabels( bool visible );
    void setLabelSpacing( double labelSpacing );

    bool contourFills() const { return m_fills; }
    bool contourLines() const { return m_lines; }
    bool contourLabels() const { return m_labels; }
    double labelSpacing() const { return m_label_spacing; }

    virtual QString posInfo( QSPt2f& pos );
    virtual bool isClicked( const QSPt2f& pos );
    virtual QSPt2f legendItemSize( QSDrv *drv );
    virtual void drawLegendItem( const QSPt2f& pos, QSDrv *drv );

    virtual void loadStateFromStream( QDataStream& stream, QSObjectFactory *factory );
    virtual void saveStateToStream( QDataStream& stream, QSObjectFactory *factory );
    enum PointElement { PointMark = 0 };
    enum LineElement { Grid = 0 };

    protected:

    static const double EPS_VALUE;

    virtual void allocRuntimeData();
    virtual void freeRuntimeData();
    virtual bool start();
    virtual bool prepare() { return true; }
    virtual bool step();
    virtual void end();
    virtual int triangles() = 0;
    virtual void getTriangle( int number, QSPt3f pts[3], double level = -10.0 ) = 0;
    virtual bool getNeighbouringTriangle( int *number, int *edge ) = 0;
    virtual int triangleAtPos( const QSPt2f& pos ) = 0;
    virtual QString vertexInfo( int t, int v, QSPt2f& pos ) = 0;
    virtual void drawGrid() = 0;
    virtual void drawPoints() = 0;
    virtual bool rectangleGrid() const { return false; }
        virtual void getRectangle( int triangleNumber, QSPt3f pts[4] ) {}


    private:
    bool m_labels;
    bool m_fills;
    bool m_lines;
    double m_label_spacing;
    struct contour_runtime_data;
    struct contour_levels_data;
    struct contour_runtime_data *d;

    void init_triangle_buffer();
    void call_prepare();
    void start_drawing_fills();
    void drawing_fills();
    void start_drawing_grid();
    void drawing_grid();
    void start_drawing_lines();
    void drawing_lines();
    void start_drawing_labels();
    void drawing_labels();
    void start_drawing_points();
    void drawing_points();

    void calculate_levels_data( contour_levels_data *data );
    void free_levels_data( contour_levels_data *data );
    void start_contour( const QSPt3f& p1, const QSPt3f& p2, int level );
    void draw_contour(  const QSPt3f& p1, const QSPt3f& p2, int level );
    void end_contour( int level );
    void draw_line_buffer( int level );

    void cut_polygon( double level, const QSPt3f *pts, int npts,
                          QSPt3f *above_pts, int *above_npts,
                          QSPt3f *under_pts, int *under_npts,
                          QSPt3f *cross_line = NULL, int *cross_npts = NULL );
    inline double distance( const QSPt2f& p1, const QSPt2f& p2 );
    inline int find_level_greater_than( double value );
    void flush_rectangle_buffer();
    bool draw_rectangle( const QSPt3f pts[4] );
   };


inline double QSContour::distance( const QSPt2f& p1, const QSPt2f& p2 )
  {
   return sqrt((p1.x-p2.x)*(p1.x-p2.x)+(p1.y-p2.y)*(p1.y-p2.y));
  }

//------------------------------------------------------------------------------------------------------------------//


class QSGriddedContour : public QSContour
  {
    Q_OBJECT
   public:
    enum Channels {
                XVector   = 0,
                YVector   = 1,     // Column vector, monotone. Size 1 x h+1
                Data      = 2
                };
    QSGriddedContour(QSAxes* parent, const char * name=0);
    virtual ~QSGriddedContour();
    virtual ColumnType columnType( int channel, int column ) const;
    virtual QString channelVariable( int channel ) const;
    virtual void loadStateFromStream( QDataStream& stream, QSObjectFactory *factory );
    virtual void saveStateToStream( QDataStream& stream, QSObjectFactory *factory );

  protected:
    virtual bool getAxisRange( QSAxis *axis, double& min, double& max );
    virtual void dataChanged( int channel = -1 );
    virtual void allocRuntimeData();
    virtual void freeRuntimeData();
    virtual int triangles();
    virtual void getTriangle( int number, QSPt3f pts[3], double level = -10.0 );
    virtual bool getNeighbouringTriangle( int *number, int *edge );
    virtual int triangleAtPos( const QSPt2f& pos );
    virtual QString vertexInfo( int t, int v, QSPt2f& pos );
    virtual void drawGrid();
    virtual void drawPoints();
    virtual bool rectangleGrid() const { return true; }
        virtual void getRectangle( int triangleNumber, QSPt3f pts[4] );

   private:
    double m_dmin;
    double m_dmax;
    bool m_evalid;
    void calculate_data_range();
    double xvector( int col );
    double yvector( int row );
    struct gridded_contour_runtime_data;
    gridded_contour_runtime_data *d;
    void load_mesh_into_buffer( int triangleNumber );
  };

//------------------------------------------------------------------------------------------------------------------//


class QSNonGriddedContour : public QSContour
 {  
    Q_OBJECT
  public:
    enum Channels {
                VXCoord   = 0,
        VYCoord   = 1,
        VZCoord   = 2,
        Triangles = 3     // size 3 x number of triangles
                };
    QSNonGriddedContour(QSAxes* parent, const char * name=0);
    virtual ~QSNonGriddedContour();
        virtual ColumnType columnType( int channel, int column ) const;
    virtual QString channelVariable( int channel ) const;
    virtual void loadStateFromStream( QDataStream& stream, QSObjectFactory *factory );
    virtual void saveStateToStream( QDataStream& stream, QSObjectFactory *factory );

  protected:
    virtual bool getAxisRange( QSAxis *axis, double& min, double& max );
    virtual void dataChanged( int channel = -1 );
    virtual bool prepare();
    virtual void allocRuntimeData();
    virtual void freeRuntimeData();
    virtual int triangles();
    virtual void getTriangle( int number, QSPt3f pts[3], double level = -10.0 );
    virtual bool getNeighbouringTriangle( int *number, int *edge );
    virtual int triangleAtPos( const QSPt2f& pos );
    virtual QString vertexInfo( int t, int v, QSPt2f& pos );
    virtual void drawGrid();
    virtual void drawPoints();

  private:
    bool m_evalid;
    double m_xmax;
    double m_ymax;
    double m_zmax;
    double m_xmin;
    double m_ymin;
    double m_zmin;
    int *m_neighbours[3];
    struct non_gridded_contour_runtime_data;    
    non_gridded_contour_runtime_data *d;
    QSPt3f point( int index );
    void search_for_neighbours();
    void make_visible_triangles_list();
 };



#endif

---