/* +------------------------------------------------------------------------+
   |                     Mobile Robot Programming Toolkit (MRPT)            |
   |                          https://www.mrpt.org/                         |
   |                                                                        |
   | Copyright (c) 2005-2022, Individual contributors, see AUTHORS file     |
   | See: https://www.mrpt.org/Authors - All rights reserved.               |
   | Released under BSD License. See: https://www.mrpt.org/License          |
   +------------------------------------------------------------------------+ */
/////////////////////////////////////////////////////////////////////////////
// Name:            mathplot.cpp
// Purpose:         Framework for plotting in wxWindows
// Original Author: David Schalig
// Maintainer:      Davide Rondini
// Contributors:    Jose Luis Blanco, Val Greene
// Created:         21/07/2003
// Last edit:       22/02/2009
// Copyright:       (c) David Schalig, Davide Rondini
// Licence:         wxWindows licence
/////////////////////////////////////////////////////////////////////////////

#pragma once

#include <memory>

// JL: This is VERY ugly, but ask MS why we cannot export a DLL class with STL
// members !!
#if defined(_MSC_VER)
#pragma warning(push)
#pragma warning(disable : 4251)
#endif

/** @file mathplot.h */
/** @mainpage wxMathPlot
	wxMathPlot is a framework for mathematical graph plotting in wxWindows.

	The framework is designed for convenience and ease of use.

	@section screenshots Screenshots
	<a href="http://wxmathplot.sourceforge.net/screenshot.shtml" >Go to the
   screenshots page.</a>

	@section overview Overview
	The heart of wxMathPlot is mpWindow, which is a 2D canvas for plot layers.
	mpWindow can be embedded as subwindow in a wxPane, a wxFrame, or any other
   wxWindow. mpWindow provides a zoomable and moveable view of the layers. The
   current view can be controlled with the mouse, the scrollbars, and a context
   menu.

	Plot layers are implementations of the abstract base class mpLayer. Those
   can be function plots, scale rulers, or any other vector data visualisation.
   wxMathPlot provides two mpLayer implementations for plotting horizontal and
   vertical rulers: mpScaleX and mpScaleY. For convenient function plotting a
   series of classes derived from mpLayer are provided, like mpFX, mpProfile,
   mpLegend and so on. These base classes already come with plot code, user's
   own functions can be implemented by overriding just one member for retrieving
   a function value.

	mpWindow has built-in support for mouse-based pan and zoom through intuitive
   combinations of buttons and the mouse wheel. It also incorporates an optional
   double buffering mechanism to avoid flicker. Plots can be easily sent to
   printer evices or exported in bitmap formats like PNG, BMP or JPEG.

	@section coding Coding conventions
	wxMathPlot sticks to wxWindow's coding conventions. All entities defined by
   wxMathPlot have the prefix <i>mp</i>.

	@section author Author and license
	wxMathPlot is published under the terms of the wxWindow license.<br>
	The original author is David Schalig <mrhill@users.sourceforge.net>.<br>
	From June 2007 the project is maintained by Davide Rondini
   <cdron77@users.sourceforge.net>.<br> Authors can be contacted via the
   wxMathPlot's homepage at https://sourceforge.net/projects/wxmathplot<br>
	Contributors:<br>
	Jose Luis Blanco, Val Greene.<br>
*/

// this definition uses windows dll to export function.
// WXDLLIMPEXP_MATHPLOT definition definition changed to WXDLLIMPEXP_MATHPLOT
// mathplot_EXPORTS will be defined by cmake
//#ifdef mathplot_EXPORTS
// #define WXDLLIMPEXP_MATHPLOT WXEXPORT
// #define WXDLLIMPEXP_DATA_MATHPLOT(type) WXEXPORT type
//#else // not making DLL
// #define WXDLLIMPEXP_MATHPLOT
// #define WXDLLIMPEXP_DATA_MATHPLOT(type) type
//#endif

// Hack for MRPT: Link as part of mrpt-gui itself.
#define WXDLLIMPEXP_MATHPLOT

#if defined(__GNUG__) && !defined(__clang__)
#pragma interface "mathplot.h"
#endif

#include <vector>

// #include <wx/wx.h>
#include <wx/dcmemory.h>
#include <wx/defs.h>
#include <wx/dynarray.h>
#include <wx/event.h>
#include <wx/image.h>
#include <wx/menu.h>
#include <wx/pen.h>
#include <wx/print.h>
#include <wx/scrolwin.h>
#include <wx/string.h>

#include <deque>

// For memory leak debug
#ifdef _WINDOWS
#ifdef _DEBUG
#include <crtdbg.h>
#define DEBUG_NEW new (_NORMAL_BLOCK, __FILE__, __LINE__)
#else
#define DEBUG_NEW new
#endif	// _DEBUG
#endif	// _WINDOWS

// Separation for axes when set close to border
#define X_BORDER_SEPARATION 40
#define Y_BORDER_SEPARATION 60

//-----------------------------------------------------------------------------
// classes
//-----------------------------------------------------------------------------

class WXDLLIMPEXP_MATHPLOT mpLayer;
class WXDLLIMPEXP_MATHPLOT mpFX;
class WXDLLIMPEXP_MATHPLOT mpFY;
class WXDLLIMPEXP_MATHPLOT mpFXY;
class WXDLLIMPEXP_MATHPLOT mpFXYVector;
class WXDLLIMPEXP_MATHPLOT mpScaleX;
class WXDLLIMPEXP_MATHPLOT mpScaleY;
class WXDLLIMPEXP_MATHPLOT mpWindow;
class WXDLLIMPEXP_MATHPLOT mpText;
class WXDLLIMPEXP_MATHPLOT mpPrintout;

/** Command IDs used by mpWindow */
enum
{
	/** Fit view to match bounding box of all layers */
	mpID_FIT = 2000,
	/** Zoom into view at clickposition / window center */
	mpID_ZOOM_IN,
	/** Zoom out */
	mpID_ZOOM_OUT,
	/** Center view on click position */
	mpID_CENTER,
	/** Lock x/y scaling aspect */
	mpID_LOCKASPECT,
	/** Shows information about the mouse commands */
	mpID_HELP_MOUSE,
	/** JL: Prints the graph */
	mpID_PRINT
};

//-----------------------------------------------------------------------------
// mpLayer
//-----------------------------------------------------------------------------

typedef enum __mp_Layer_Type
{
	/** Layer type undefined */
	mpLAYER_UNDEF,
	/** Axis type layer */
	mpLAYER_AXIS,
	/** Plot type layer */
	mpLAYER_PLOT,
	/** Info box type layer */
	mpLAYER_INFO,
	/** Bitmap type layer */
	mpLAYER_BITMAP
} mpLayerType;

/** Plot layer, abstract base class.
	Any number of mpLayer implementations can be attached to mpWindow.
	Examples for mpLayer implementations are function graphs, or scale rulers.

	For convenience mpLayer defines a name, a font (wxFont), a pen (wxPen),
	and a continuity property (bool) as class members.
	The default values at constructor are the default font, a black pen, and
	 continuity set to false (draw separate points).
	These may or may not be used by implementations.
*/
class WXDLLIMPEXP_MATHPLOT mpLayer : public wxObject
{
   public:
	mpLayer();

	~mpLayer() override = default;

	/** Check whether this layer has a bounding box.
		The default implementation returns \a TRUE. Override and return
		FALSE if your mpLayer implementation should be ignored by the
	   calculation of the global bounding box for all layers in a mpWindow.
		@retval TRUE Has bounding box
		@retval FALSE Has not bounding box
	*/
	virtual bool HasBBox() { return TRUE; }

	/** Check whether the layer is an info box.
		The default implementation returns \a FALSE. It is overrided to \a TRUE
	   for mpInfoLayer class and its derivative. It is necessary to define mouse
	   actions behaviour over info boxes.
		@return whether the layer is an info boxes
		@sa mpInfoLayer::IsInfo */
	virtual bool IsInfo() { return false; };

	/** Get inclusive left border of bounding box.
		@return Value
	*/
	virtual double GetMinX() { return -1.0; }

	/** Get inclusive right border of bounding box.
		@return Value
	*/
	virtual double GetMaxX() { return 1.0; }

	/** Get inclusive bottom border of bounding box.
		@return Value
	*/
	virtual double GetMinY() { return -1.0; }

	/** Get inclusive top border of bounding box.
		@return Value
	*/
	virtual double GetMaxY() { return 1.0; }

	/** Plot given view of layer to the given device context.
		An implementation of this function has to transform layer coordinates to
		wxDC coordinates based on the view parameters retrievable from the
	mpWindow passed in \a w. Note that the public methods of mpWindow: x2p,y2p
	and p2x,p2y are already provided which transform layer coordinates to DC
	pixel coordinates, and <b>user code should rely on them</b> for portability
	and future changes to be applied transparently, instead of implementing the
	following formulas manually.

	The passed device context \a dc has its coordinate origin set to the
	top-left corner of the visible area (the default). The coordinate
	orientation is as shown in the following picture: <pre> (wxDC origin 0,0)
			   x-------------> ascending X ----------------+
			   |                                           |
			   |                                           |
			   V ascending Y                               |
			   |                                           |
			   |                                           |
			   |                                           |
			   +-------------------------------------------+  <-- right-bottom
	corner of the mpWindow visible area.
		</pre>
		Note that Y ascends in downward direction, whereas the usual vertical
	orientation for mathematical plots is vice versa. Thus Y-orientation will be
	swapped usually, when transforming between wxDC and mpLayer coordinates.
	This change of coordinates is taken into account in the methods
	p2x,p2y,x2p,y2p.

		<b> Rules for transformation between mpLayer and wxDC coordinates </b>
		@code
		dc_X = (layer_X - mpWindow::GetPosX()) * mpWindow::GetScaleX()
		dc_Y = (mpWindow::GetPosY() - layer_Y) * mpWindow::GetScaleY() //
	swapping Y-orientation

		layer_X = (dc_X / mpWindow::GetScaleX()) + mpWindow::GetPosX() // scale
	guaranteed to be not 0 layer_Y = mpWindow::GetPosY() - (dc_Y /
	mpWindow::GetScaleY()) // swapping Y-orientation
		@endcode

		@param dc Device context to plot to.
		@param w  View to plot. The visible area can be retrieved from this
	object.
	@sa mpWindow::p2x,mpWindow::p2y,mpWindow::x2p,mpWindow::y2p
	*/
	virtual void Plot(wxDC& dc, mpWindow& w) = 0;

	/** Get layer name.
		@return Name
	*/
	wxString GetName() const { return m_name; }

	/** Get font set for this layer.
		@return Font
	*/
	const wxFont& GetFont() const { return m_font; }

	/** Get pen set for this layer.
		@return Pen
	*/
	const wxPen& GetPen() const { return m_pen; }

	/** Set the 'continuity' property of the layer (true:draws a continuous
	 * line, false:draws separate points).
	 * @sa GetContinuity
	 */
	void SetContinuity(bool continuity) { m_continuous = continuity; }

	/** Gets the 'continuity' property of the layer.
	 * @sa SetContinuity
	 */
	bool GetContinuity() const { return m_continuous; }

	/** Shows or hides the text label with the name of the layer (default is
	 * visible).
	 */
	void ShowName(bool show) { m_showName = show; };

	/** Set layer name
		@param name Name, will be copied to internal class member
	*/
	void SetName(wxString name) { m_name = name; }

	/** Set layer font
		@param font Font, will be copied to internal class member
	*/
	void SetFont(wxFont& font) { m_font = font; }

	/** Set layer pen
		@param pen Pen, will be copied to internal class member
	*/
	void SetPen(wxPen pen) { m_pen = pen; }

	/** Set Draw mode: inside or outside margins. Default is outside, which
	   allows the layer to draw up to the mpWindow border.
		@param drawModeOutside The draw mode to be set */
	void SetDrawOutsideMargins(bool drawModeOutside)
	{
		m_drawOutsideMargins = drawModeOutside;
	};

	/** Get Draw mode: inside or outside margins.
		@return The draw mode */
	bool GetDrawOutsideMargins() { return m_drawOutsideMargins; };

	/** Get a small square bitmap filled with the colour of the pen used in the
	   layer. Useful to create legends or similar reference to the layers.
		@param side side length in pixels
		@return a wxBitmap filled with layer's colour */
	wxBitmap GetColourSquare(int side = 16);

	/** Get layer type: a Layer can be of different types: plot lines, axis,
	   info boxes, etc, this method returns the right value.
		@return An integer indicating layer type */
	mpLayerType GetLayerType() { return m_type; };

	/** Checks whether the layer is visible or not.
		@return \a true if visible */
	bool IsVisible() { return m_visible; };

	/** Sets layer visibility.
		@param show visibility bool. */
	void SetVisible(bool show) { m_visible = show; };

	/** Get brush set for this layer.
		@return brush. */
	const wxBrush& GetBrush() const { return m_brush; };

	/** Set layer brush
		@param brush brush, will be copied to internal class member	*/
	void SetBrush(wxBrush brush) { m_brush = brush; };

   protected:
	/** Layer's font */
	wxFont m_font;
	/** Layer's pen */
	wxPen m_pen;
	/** Layer's brush */
	wxBrush m_brush;
	/** Layer's name */
	wxString m_name;
	/** Specify if the layer will be plotted as a continuous line or a set of
	 * points. */
	bool m_continuous;
	/** States whether the name of the layer must be shown (default is true). */
	bool m_showName;
	/** select if the layer should draw only inside margins or over all DC */
	bool m_drawOutsideMargins;
	/** Define layer type, which is assigned by constructor */
	mpLayerType m_type{mpLAYER_UNDEF};
	/** Toggles layer visibility */
	bool m_visible;
	DECLARE_DYNAMIC_CLASS(mpLayer)
};

//-----------------------------------------------------------------------------
// mpInfoLayer
//-----------------------------------------------------------------------------

/** @class mpInfoLayer
	@brief Base class to create small rectangular info boxes
	mpInfoLayer is the base class to create a small rectangular info box in
   transparent overlay over plot layers. It is used to implement objects like
   legends.
*/
class WXDLLIMPEXP_MATHPLOT mpInfoLayer : public mpLayer
{
   public:
	/** Default constructor. */
	mpInfoLayer();

	/** Complete constructor.
		@param rect Sets the initial size rectangle of the layer.
		@param brush pointer to a fill brush. Default is transparent */
	mpInfoLayer(wxRect rect, const wxBrush* brush = wxTRANSPARENT_BRUSH);

	/** Destructor */
	~mpInfoLayer() override;

	/** Updates the content of the info box. Should be overidden by derived
	   classes. Update may behave in different ways according to the type of
	   event which called it.
		@param w parent mpWindow from which to obtain informations
		@param event The event which called the update. */
	virtual void UpdateInfo(mpWindow& w, wxEvent& event);

	/** mpInfoLayer has not bounding box. @sa mpLayer::HasBBox
		@return always \a FALSE */
	bool HasBBox() override { return false; };

	/** Plot method. Can be overidden by derived classes.
		@param dc the device content where to plot
		@param w the window to plot
		@sa mpLayer::Plot */
	void Plot(wxDC& dc, mpWindow& w) override;

	/** Specifies that this is an Info box layer.
		@return always \a TRUE
		@sa mpLayer::IsInfo */
	bool IsInfo() override { return true; };

	/** Checks whether a point is inside the info box rectangle.
		@param point The point to be checked
		@return \a true if the point is inside the bounding box */
	virtual bool Inside(wxPoint& point);

	/** Moves the layer rectangle of given pixel deltas.
		@param delta The wxPoint container for delta coordinates along x and y.
	   Units are in pixels. */
	virtual void Move(wxPoint delta);

	/** Updates the rectangle reference point. Used by internal methods of
	 * mpWindow to correctly move mpInfoLayers. */
	virtual void UpdateReference();

	/** Returns the position of the upper left corner of the box (in pixels)
		@return The rectangle position */
	wxPoint GetPosition() const;

	/** Returns the size of the box (in pixels)
		@return The rectangle size */
	wxSize GetSize();

	/** Returns the current rectangle coordinates.
		@return The info layer rectangle */
	const wxRect& GetRectangle() { return m_dim; };

   protected:
	/** The bounding rectangle of the box. It may be resized dynamically by the
	 * Plot method. */
	wxRect m_dim;
	/** Holds the reference point for movements */
	wxPoint m_reference;
	/** The brush to be used for the background */
	wxBrush m_brush;
	/** Holds the mpWindow size. Used to rescale position when window is
	 * resized. */
	int m_winX, m_winY;

	DECLARE_DYNAMIC_CLASS(mpInfoLayer)
};

/** @class mpInfoCoords
	@brief Implements an overlay box which shows the mouse coordinates in plot
   units. When an mpInfoCoords layer is activated, when mouse is moved over the
   mpWindow, its coordinates (in mpWindow units, not pixels) are continuously
   reported inside the layer box. */
class WXDLLIMPEXP_MATHPLOT mpInfoCoords : public mpInfoLayer
{
   public:
	/** Default constructor */
	mpInfoCoords();
	/** Complete constructor, setting initial rectangle and background brush.
		@param rect The initial bounding rectangle.
		@param brush The wxBrush to be used for box background: default is
	   transparent */
	mpInfoCoords(wxRect rect, const wxBrush* brush = wxTRANSPARENT_BRUSH);

	/** Default destructor */
	~mpInfoCoords() override;

	/** Updates the content of the info box. It is used to update coordinates.
		@param w parent mpWindow from which to obtain information
		@param event The event which called the update. */
	void UpdateInfo(mpWindow& w, wxEvent& event) override;

	/** Plot method.
		@param dc the device content where to plot
		@param w the window to plot
		@sa mpLayer::Plot */
	void Plot(wxDC& dc, mpWindow& w) override;

   protected:
	/** string holding the coordinates to be drawn. */
	wxString m_content;
};

/** @class mpInfoLegend
	@brief Implements the legend to be added to the plot
	This layer allows you to add a legend to describe the plots in the window.
   The legend uses the layer name as a label, and displays only layers of type
   mpLAYER_PLOT. */
class WXDLLIMPEXP_MATHPLOT mpInfoLegend : public mpInfoLayer
{
   public:
	/** Default constructor */
	mpInfoLegend();

	/** Complete constructor, setting initial rectangle and background brush.
		@param rect The initial bounding rectangle.
		@param brush The wxBrush to be used for box background: default is
	   transparent
		@sa mpInfoLayer::mpInfoLayer */
	mpInfoLegend(wxRect rect, const wxBrush* brush = wxTRANSPARENT_BRUSH);

	/**  Default destructor */
	~mpInfoLegend() override;

	/** Updates the content of the info box. Unused in this class.
		@param w parent mpWindow from which to obtain information
		@param event The event which called the update. */
	void UpdateInfo(mpWindow& w, wxEvent& event) override;

	/** Plot method.
		@param dc the device content where to plot
		@param w the window to plot
		@sa mpLayer::Plot */
	void Plot(wxDC& dc, mpWindow& w) override;

   protected:
};

//-----------------------------------------------------------------------------
// mpLayer implementations - functions
//-----------------------------------------------------------------------------

/** @name Label alignment constants
@{*/

/** @internal */
#define mpALIGNMASK 0x03
/** Aligns label to the right. For use with mpFX. */
#define mpALIGN_RIGHT 0x00
/** Aligns label to the center. For use with mpFX and mpFY. */
#define mpALIGN_CENTER 0x01
/** Aligns label to the left. For use with mpFX. */
#define mpALIGN_LEFT 0x02
/** Aligns label to the top. For use with mpFY. */
#define mpALIGN_TOP mpALIGN_RIGHT
/** Aligns label to the bottom. For use with mpFY. */
#define mpALIGN_BOTTOM mpALIGN_LEFT
/** Aligns X axis to bottom border. For mpScaleX */
#define mpALIGN_BORDER_BOTTOM 0x04
/** Aligns X axis to top border. For mpScaleX */
#define mpALIGN_BORDER_TOP 0x05
/** Set label for X axis in normal mode */
#define mpX_NORMAL 0x00
/** Set label for X axis in time mode: the value is represented as
 * minutes:seconds.milliseconds if time is less than 2 minutes,
 * hours:minutes:seconds otherwise. */
#define mpX_TIME 0x01
/** Set label for X axis in hours mode: the value is always represented as
 * hours:minutes:seconds. */
#define mpX_HOURS 0x02
/** Set label for X axis in date mode: the value is always represented as
 * yyyy-mm-dd. */
#define mpX_DATE 0x03
/** Set label for X axis in datetime mode: the value is always represented as
 * yyyy-mm-ddThh:mm:ss. */
#define mpX_DATETIME 0x04
/** Aligns Y axis to left border. For mpScaleY */
#define mpALIGN_BORDER_LEFT mpALIGN_BORDER_BOTTOM
/** Aligns Y axis to right border. For mpScaleY */
#define mpALIGN_BORDER_RIGHT mpALIGN_BORDER_TOP
/** Aligns label to north-east. For use with mpFXY. */
#define mpALIGN_NE 0x00
/** Aligns label to north-west. For use with mpFXY. */
#define mpALIGN_NW 0x01
/** Aligns label to south-west. For use with mpFXY. */
#define mpALIGN_SW 0x02
/** Aligns label to south-east. For use with mpFXY. */
#define mpALIGN_SE 0x03

/*@}*/

/** @name mpLayer implementations - functions
@{*/

/** Abstract base class providing plot and labeling functionality for functions
   F:X->Y. Override mpFX::GetY to implement a function. Optionally implement a
   constructor and pass a name (label) and a label alignment to the constructor
   mpFX::mpFX. If the layer name is empty, no label will be plotted.
*/
class WXDLLIMPEXP_MATHPLOT mpFX : public mpLayer
{
   public:
	/** @param name  Label
		@param flags Label alignment, pass one of #mpALIGN_RIGHT,
	   #mpALIGN_CENTER, #mpALIGN_LEFT.
	*/
	mpFX(wxString name = wxEmptyString, int flags = mpALIGN_RIGHT);

	/** Get function value for argument.
		Override this function in your implementation.
		@param x Argument
		@return Function value
	*/
	virtual double GetY(double x) = 0;

	/** Layer plot handler.
		This implementation will plot the function in the visible area and
		put a label according to the aligment specified.
	*/
	void Plot(wxDC& dc, mpWindow& w) override;

   protected:
	/** Holds label alignment */
	int m_flags = 0;

	DECLARE_DYNAMIC_CLASS(mpFX)
};

/** Abstract base class providing plot and labeling functionality for functions
   F:Y->X. Override mpFY::GetX to implement a function. Optionally implement a
   constructor and pass a name (label) and a label alignment to the constructor
   mpFY::mpFY. If the layer name is empty, no label will be plotted.
*/
class WXDLLIMPEXP_MATHPLOT mpFY : public mpLayer
{
   public:
	/** @param name  Label
		@param flags Label alignment, pass one of #mpALIGN_BOTTOM,
	   #mpALIGN_CENTER, #mpALIGN_TOP.
	*/
	mpFY(wxString name = wxEmptyString, int flags = mpALIGN_TOP);

	/** Get function value for argument.
		Override this function in your implementation.
		@param y Argument
		@return Function value
	*/
	virtual double GetX(double y) = 0;

	/** Layer plot handler.
		This implementation will plot the function in the visible area and
		put a label according to the aligment specified.
	*/
	void Plot(wxDC& dc, mpWindow& w) override;

   protected:
	/** Holds label alignment */
	int m_flags = 0;

	DECLARE_DYNAMIC_CLASS(mpFY)
};

/** Abstract base class providing plot and labeling functionality for a locus
   plot F:N->X,Y. Locus argument N is assumed to be in range 0 .. MAX_N, and
   implicitly derived by enumerating all locus values. Override mpFXY::Rewind
   and mpFXY::GetNextXY to implement a locus. Optionally implement a constructor
   and pass a name (label) and a label alignment to the constructor
   mpFXY::mpFXY. If the layer name is empty, no label will be plotted.
*/
class WXDLLIMPEXP_MATHPLOT mpFXY : public mpLayer
{
   public:
	/** @param name  Label
		@param flags Label alignment, pass one of #mpALIGN_NE, #mpALIGN_NW,
	   #mpALIGN_SW, #mpALIGN_SE.
	*/
	mpFXY(wxString name = wxEmptyString, int flags = mpALIGN_NE);

	/** Rewind value enumeration with mpFXY::GetNextXY.
		Override this function in your implementation.
	*/
	virtual void Rewind() = 0;

	/** Get locus value for next N.
		Override this function in your implementation.
		@param x Returns X value
		@param y Returns Y value
	*/
	virtual bool GetNextXY(double& x, double& y) = 0;

	/** Layer plot handler.
		This implementation will plot the locus in the visible area and
		put a label according to the alignment specified.
	*/
	void Plot(wxDC& dc, mpWindow& w) override;

   protected:
	/** Holds label alignment */
	int m_flags = 0;

	// Data to calculate label positioning
	wxCoord maxDrawX = 0, minDrawX = 0, maxDrawY = 0, minDrawY = 0;
	// int drawnPoints;

	/** Update label positioning data
		@param xnew New x coordinate
		@param ynew New y coordinate
	*/
	void UpdateViewBoundary(wxCoord xnew, wxCoord ynew);

	DECLARE_DYNAMIC_CLASS(mpFXY)
};

/** Abstract base class providing plot and labeling functionality for functions
   F:Y->X. Override mpProfile::GetX to implement a function. This class is
   similar to mpFY, but the Plot method is different. The plot is in fact
   represented by lines instead of points, which gives best rendering of
   rapidly-varying functions, and in general, data which are not so close one to
   another. Optionally implement a constructor and pass a name (label) and a
   label alignment to the constructor mpProfile::mpProfile. If the layer name is
   empty, no label will be plotted.
*/
class WXDLLIMPEXP_MATHPLOT mpProfile : public mpLayer
{
   public:
	/** @param name  Label
		@param flags Label alignment, pass one of #mpALIGN_BOTTOM,
	   #mpALIGN_CENTER, #mpALIGN_TOP.
	*/
	mpProfile(wxString name = wxEmptyString, int flags = mpALIGN_TOP);

	/** Get function value for argument.
		Override this function in your implementation.
		@param x Argument
		@return Function value
	*/
	virtual double GetY(double x) = 0;

	/** Layer plot handler.
		This implementation will plot the function in the visible area and
		put a label according to the aligment specified.
	*/
	void Plot(wxDC& dc, mpWindow& w) override;

   protected:
	/** Holds label alignment */
	int m_flags = 0;

	DECLARE_DYNAMIC_CLASS(mpProfile)
};

/*@}*/

//-----------------------------------------------------------------------------
// mpLayer implementations - furniture (scales, ...)
//-----------------------------------------------------------------------------

/** @name mpLayer implementations - furniture (scales, ...)
@{*/

/** Plot layer implementing a x-scale ruler.
	The ruler is fixed at Y=0 in the coordinate system. A label is plotted at
	the bottom-right hand of the ruler. The scale numbering automatically
	adjusts to view and zoom factor.
*/
class WXDLLIMPEXP_MATHPLOT mpScaleX : public mpLayer
{
   public:
	/** Full constructor.
		@param name Label to plot by the ruler
		@param flags Set the position of the scale with respect to the window.
		@param ticks Select ticks or grid. Give TRUE (default) for drawing axis
	   ticks, FALSE for drawing the grid.
		@param type mpX_NORMAL for normal labels, mpX_TIME for time axis in
	   hours, minutes, seconds. */
	mpScaleX(
		wxString name = wxT("X"), int flags = mpALIGN_CENTER, bool ticks = true,
		unsigned int type = mpX_NORMAL);

	/** Layer plot handler.
		This implementation will plot the ruler adjusted to the visible area. */
	void Plot(wxDC& dc, mpWindow& w) override;

	/** Check whether this layer has a bounding box.
		This implementation returns \a FALSE thus making the ruler invisible
		to the plot layer bounding box calculation by mpWindow. */
	bool HasBBox() override { return FALSE; }

	/** Set X axis alignment.
		@param align alignment (choose between mpALIGN_BORDER_BOTTOM,
	   mpALIGN_BOTTOM, mpALIGN_CENTER, mpALIGN_TOP, mpALIGN_BORDER_TOP */
	void SetAlign(int align) { m_flags = align; };

	/** Set X axis ticks or grid
		@param ticks TRUE to plot axis ticks, FALSE to plot grid. */
	void SetTicks(bool ticks) { m_ticks = ticks; };

	/** Get X axis ticks or grid
		@return TRUE if plot is drawing axis ticks, FALSE if the grid is active.
	 */
	bool GetTicks() { return m_ticks; };

	/** Get X axis label view mode.
		@return mpX_NORMAL for normal labels, mpX_TIME for time axis in hours,
	   minutes, seconds. */
	unsigned int GetLabelMode() { return m_labelType; };

	/** Set X axis label view mode.
		@param mode mpX_NORMAL for normal labels, mpX_TIME for time axis in
	   hours, minutes, seconds. */
	void SetLabelMode(unsigned int mode) { m_labelType = mode; };

	/** Set X axis Label format (used for mpX_NORMAL draw mode).
		@param format The format string */
	void SetLabelFormat(const wxString& format) { m_labelFormat = format; };

	/** Get X axis Label format (used for mpX_NORMAL draw mode).
	@return The format string */
	const wxString& SetLabelFormat() { return m_labelFormat; };

   protected:
	/** Flag for axis alignment */
	int m_flags = 0;
	/** Flag to toggle between ticks or grid */
	bool m_ticks;
	/** Select labels mode: mpX_NORMAL for normal labels, mpX_TIME for time axis
	 * in hours, minutes, seconds */
	unsigned int m_labelType;
	/** Format string used to print labels */
	wxString m_labelFormat;

	DECLARE_DYNAMIC_CLASS(mpScaleX)
};

/** Plot layer implementing a y-scale ruler.
	If align is set to mpALIGN_CENTER, the ruler is fixed at X=0 in the
   coordinate system. If the align is set to mpALIGN_TOP or mpALIGN_BOTTOM, the
   axis is always drawn respectively at top or bottom of the window. A label is
   plotted at the top-right hand of the ruler. The scale numbering automatically
	adjusts to view and zoom factor.
*/
class WXDLLIMPEXP_MATHPLOT mpScaleY : public mpLayer
{
   public:
	/** @param name Label to plot by the ruler
		@param flags Set position of the scale respect to the window.
		@param ticks Select ticks or grid. Give TRUE (default) for drawing axis
	   ticks, FALSE for drawing the grid */
	mpScaleY(
		wxString name = wxT("Y"), int flags = mpALIGN_CENTER,
		bool ticks = true);

	/** Layer plot handler.
		This implementation will plot the ruler adjusted to the visible area.
	*/
	void Plot(wxDC& dc, mpWindow& w) override;

	/** Check whether this layer has a bounding box.
		This implementation returns \a FALSE thus making the ruler invisible
		to the plot layer bounding box calculation by mpWindow.
	*/
	bool HasBBox() override { return FALSE; }

	/** Set Y axis alignment.
		@param align alignment (choose between mpALIGN_BORDER_LEFT,
	   mpALIGN_LEFT, mpALIGN_CENTER, mpALIGN_RIGHT, mpALIGN_BORDER_RIGHT) */
	void SetAlign(int align) { m_flags = align; };

	/** Set Y axis ticks or grid
		@param ticks TRUE to plot axis ticks, FALSE to plot grid. */
	void SetTicks(bool ticks) { m_ticks = ticks; };

	/** Get Y axis ticks or grid
		@return TRUE if plot is drawing axis ticks, FALSE if the grid is active.
	 */
	bool GetTicks() { return m_ticks; };

	/** Set Y axis Label format.
	@param format The format string */
	void SetLabelFormat(const wxString& format) { m_labelFormat = format; };

	/** Get Y axis Label format.
	@return The format string */
	const wxString& SetLabelFormat() { return m_labelFormat; };

   protected:
	/** Flag for axis alignment */
	int m_flags = 0;
	/** Flag to toggle between ticks or grid */
	bool m_ticks;
	/** Format string used to print labels */
	wxString m_labelFormat;

	DECLARE_DYNAMIC_CLASS(mpScaleY)
};

//-----------------------------------------------------------------------------
// mpWindow
//-----------------------------------------------------------------------------

/** @name Constants defining mouse modes for mpWindow
@{*/

/** Mouse panning drags the view. Mouse mode for mpWindow. */
#define mpMOUSEMODE_DRAG 0
/** Mouse panning creates a zoom box. Mouse mode for mpWindow. */
#define mpMOUSEMODE_ZOOMBOX 1

/*@}*/
/** Define the type for the list of layers inside mpWindow */
// WX_DECLARE_HASH_MAP( int, mpLayer*, wxIntegerHash, wxIntegerEqual,
// wxLayerList );
typedef std::deque<mpLayer*> wxLayerList;

/** Canvas for plotting mpLayer implementations.

	This class defines a zoomable and moveable 2D plot canvas. Any number
	of mpLayer implementations (scale rulers, function plots, ...) can be
	attached using mpWindow::AddLayer.

	The canvas window provides a context menu with actions for navigating the
   view. The context menu can be retrieved with mpWindow::GetPopupMenu, e.g. for
   extending it externally.

	Since wxMathPlot version 0.03, the mpWindow incorporates the following
   features:
		- Mouse based pan/zoom (Default=enabled): Can be set with
   EnableMousePanZoom.

	The mouse commands can be visualized by the user through the popup menu, and
   are:
		- Mouse Move+CTRL: Pan (Move)
		- Mouse Wheel: Vertical scroll
		- Mouse Wheel+SHIFT: Horizontal scroll
		- Mouse Wheel UP+CTRL: Zoom in
		- Mouse Wheel DOWN+CTRL: Zoom out

*/
class WXDLLIMPEXP_MATHPLOT mpWindow : public wxWindow
{
   public:
	mpWindow() = default;
	mpWindow(
		wxWindow* parent, wxWindowID id, const wxPoint& pos = wxDefaultPosition,
		const wxSize& size = wxDefaultSize, long flags = 0);
	~mpWindow() override;

	/** Get reference to context menu of the plot canvas.
		@return Pointer to menu. The menu can be modified.
	*/
	wxMenu* GetPopupMenu() { return &m_popmenu; }

	/** Add a plot layer to the canvas.
		@param layer Pointer to layer. The mpLayer object will get under control
	   of mpWindow, i.e. it will be delete'd on mpWindow destruction
		@param refreshDisplay States whether to refresh the display (UpdateAll)
	   after adding the layer.
		@retval TRUE Success
		@retval FALSE Failure due to out of memory.
	*/
	bool AddLayer(mpLayer* layer, bool refreshDisplay = true);

	/** Remove a plot layer from the canvas.
		@param layer Pointer to layer. The mpLayer object will be destructed
	   using delete.
		@param alsoDeleteObject If set to true, the mpLayer object will be also
	   "deleted", not just removed from the internal list.
		@param refreshDisplay States whether to refresh the display (UpdateAll)
	   after removing the layer.
		@return true if layer is deleted correctly

		N.B. Only the layer reference in the mpWindow is deleted, the layer
	   object still exists!
	*/
	bool DelLayer(
		mpLayer* layer, bool alsoDeleteObject = false,
		bool refreshDisplay = true);

	/** Remove all layers from the plot.
		@param alsoDeleteObject If set to true, the mpLayer objects will be also
	   "deleted", not just removed from the internal list.
		@param refreshDisplay States whether to refresh the display (UpdateAll)
	   after removing the layers.
	*/
	void DelAllLayers(bool alsoDeleteObject, bool refreshDisplay = true);

	/*! Get the layer in list position indicated.
		N.B. You <i>must</i> know the index of the layer inside the list!
		@param position position of the layer in the layers list
		@return pointer to mpLayer
	*/
	mpLayer* GetLayer(int position);

	/*! Get the layer by its name (case sensitive).
		@param name The name of the layer to retrieve
		@return A pointer to the mpLayer object, or nullptr if not found.
	*/
	mpLayer* GetLayerByName(const wxString& name);

	/** Get current view's X scale.
		See @ref mpLayer::Plot "rules for coordinate transformation"
		@return Scale
	*/
	double GetXscl() { return m_scaleX; }
	double GetScaleX() const
	{
		return m_scaleX;
	};	// Schaling's method: maybe another method esists with the same name

	/** Get current view's Y scale.
		See @ref mpLayer::Plot "rules for coordinate transformation"
		@return Scale
	*/
	double GetYscl() const { return m_scaleY; }
	double GetScaleY() const
	{
		return m_scaleY;
	}  // Schaling's method: maybe another method exists with the same name

	/** Get current view's X position.
		See @ref mpLayer::Plot "rules for coordinate transformation"
		@return X Position in layer coordinate system, that corresponds to the
	   center point of the view.
	*/
	double GetXpos() const { return m_posX; }
	double GetPosX() const { return m_posX; }

	/** Get current view's Y position.
		See @ref mpLayer::Plot "rules for coordinate transformation"
		@return Y Position in layer coordinate system, that corresponds to the
	   center point of the view.
	*/
	double GetYpos() const { return m_posY; }
	double GetPosY() const { return m_posY; }

	/** Get current view's X dimension in device context units.
		Usually this is equal to wxDC::GetSize, but it might differ thus mpLayer
		implementations should rely on the value returned by the function.
		See @ref mpLayer::Plot "rules for coordinate transformation"
		@return X dimension.
	*/
	int GetScrX() const { return m_scrX; }
	int GetXScreen() const { return m_scrX; }

	/** Get current view's Y dimension in device context units.
		Usually this is equal to wxDC::GetSize, but it might differ thus mpLayer
		implementations should rely on the value returned by the function.
		See @ref mpLayer::Plot "rules for coordinate transformation"
		@return Y dimension.
	*/
	int GetScrY() const { return m_scrY; }
	int GetYScreen() const { return m_scrY; }

	/** Set current view's X scale and refresh display.
		@param scaleX New scale, must not be 0.
	*/
	void SetScaleX(double scaleX);

	/** Set current view's Y scale and refresh display.
		@param scaleY New scale, must not be 0.
	*/
	void SetScaleY(double scaleY)
	{
		if (scaleY != 0) m_scaleY = scaleY;
		UpdateAll();
	}

	/** Set current view's X position and refresh display.
		@param posX New position that corresponds to the center point of the
	   view.
	*/
	void SetPosX(double posX)
	{
		m_posX = posX;
		UpdateAll();
	}

	/** Set current view's Y position and refresh display.
		@param posY New position that corresponds to the center point of the
	   view.
	*/
	void SetPosY(double posY)
	{
		m_posY = posY;
		UpdateAll();
	}

	/** Set current view's X and Y position and refresh display.
		@param posX New position that corresponds to the center point of the
	   view.
		@param posY New position that corresponds to the center point of the
	   view.
	*/
	void SetPos(double posX, double posY)
	{
		m_posX = posX;
		m_posY = posY;
		UpdateAll();
	}

	/** Set current view's dimensions in device context units.
		Needed by plotting functions. It doesn't refresh display.
		@param scrX New position that corresponds to the center point of the
	   view.
		@param scrY New position that corresponds to the center point of the
	   view.
	*/
	void SetScr(int scrX, int scrY)
	{
		m_scrX = scrX;
		m_scrY = scrY;
	}

	/** Converts mpWindow (screen) pixel coordinates into graph (floating point)
	 * coordinates, using current mpWindow position and scale.
	 * @sa p2y,x2p,y2p */
	//     double p2x(wxCoord pixelCoordX, bool drawOutside = true ); // {
	//     return m_posX + pixelCoordX/m_scaleX; }
	inline double p2x(wxCoord pixelCoordX)
	{
		return m_posX + pixelCoordX / m_scaleX;
	}

	/** Converts mpWindow (screen) pixel coordinates into graph (floating point)
	 * coordinates, using current mpWindow position and scale.
	 * @sa p2x,x2p,y2p */
	//     double p2y(wxCoord pixelCoordY, bool drawOutside = true ); //{ return
	//     m_posY - pixelCoordY/m_scaleY; }
	inline double p2y(wxCoord pixelCoordY)
	{
		return m_posY - pixelCoordY / m_scaleY;
	}

	/** Converts graph (floating point) coordinates into mpWindow (screen) pixel
	 * coordinates, using current mpWindow position and scale.
	 * @sa p2x,p2y,y2p */
	//     wxCoord x2p(double x, bool drawOutside = true); // { return (wxCoord)
	//     ( (x-m_posX) * m_scaleX); }
	inline wxCoord x2p(double x) { return (wxCoord)((x - m_posX) * m_scaleX); }

	/** Converts graph (floating point) coordinates into mpWindow (screen) pixel
	 * coordinates, using current mpWindow position and scale.
	 * @sa p2x,p2y,x2p */
	//     wxCoord y2p(double y, bool drawOutside = true); // { return (wxCoord)
	//     ( (m_posY-y) * m_scaleY); }
	inline wxCoord y2p(double y) { return (wxCoord)((m_posY - y) * m_scaleY); }

	/** Enable/disable the feature of pan/zoom with the mouse (default=enabled)
	 */
	void EnableMousePanZoom(bool enabled) { m_enableMouseNavigation = enabled; }

	/** Enable or disable X/Y scale aspect locking for the view.
		@note Explicit calls to mpWindow::SetScaleX and mpWindow::SetScaleY will
	   set an unlocked aspect, but any other action changing the view scale will
			  lock the aspect again.
	*/
	void LockAspect(bool enable = TRUE);

	/** Checks whether the X/Y scale aspect is locked.
		@retval TRUE Locked
		@retval FALSE Unlocked
	*/
	inline bool IsAspectLocked() { return m_lockaspect; }

	/** Set view to fit global bounding box of all plot layers and refresh
	   display. Scale and position will be set to show all attached mpLayers.
		The X/Y scale aspect lock is taken into account.
	*/
	void Fit() override;

	/** Set view to fit a given bounding box and refresh display.
		The X/Y scale aspect lock is taken into account.
	If provided, the parameters printSizeX and printSizeY are taken as the DC
	size, and the pixel scales are computed accordingly. Also, in this case the
	passed borders are not saved as the "desired borders", since this use will
	be invoked only when printing.
	*/
	void Fit(
		double xMin, double xMax, double yMin, double yMax,
		wxCoord* printSizeX = nullptr, wxCoord* printSizeY = nullptr);

	/** Zoom into current view and refresh display
	 * @param centerPoint The point (pixel coordinates) that will stay in the
	 * same position on the screen after the zoom (by default, the center of the
	 * mpWindow).
	 */
	void ZoomIn(const wxPoint& centerPoint = wxDefaultPosition);

	/** Zoom out current view and refresh display
	 * @param centerPoint The point (pixel coordinates) that will stay in the
	 * same position on the screen after the zoom (by default, the center of the
	 * mpWindow).
	 */
	void ZoomOut(const wxPoint& centerPoint = wxDefaultPosition);

	/** Zoom in current view along X and refresh display */
	void ZoomInX();
	/** Zoom out current view along X and refresh display */
	void ZoomOutX();
	/** Zoom in current view along Y and refresh display */
	void ZoomInY();
	/** Zoom out current view along Y and refresh display */
	void ZoomOutY();

	/** Zoom view fitting given coordinates to the window (p0 and p1 do not need
	 * to be in any specific order) */
	void ZoomRect(wxPoint p0, wxPoint p1);

	/** Refresh display */
	void UpdateAll();

	// Added methods by Davide Rondini

	/** Counts the number of plot layers, excluding axes or text: this is to
	   count only the layers which have a bounding box. \return The number of
	   profiles plotted.
	*/
	unsigned int CountLayers();

	/** Counts the number of plot layers, whether or not they have a bounding
	   box. \return The number of layers in the mpWindow. */
	size_t CountAllLayers() { return m_layers.size(); };

	/** Draws the mpWindow on a page for printing
		\param print the mpPrintout where to print the graph */
	// void PrintGraph(mpPrintout *print);

	void ShowPrintDialog()
	{
		wxCommandEvent dum;
		OnPrintMenu(dum);
	}

	/** Returns the left-border layer coordinate that the user wants the
	 * mpWindow to show (it may be not exactly the actual shown coordinate in
	 * the case of locked aspect ratio).
	 * @sa Fit
	 */
	double GetDesiredXmin() { return m_desiredXmin; }

	/** Returns the right-border layer coordinate that the user wants the
	 * mpWindow to show (it may be not exactly the actual shown coordinate in
	 * the case of locked aspect ratio).
	 * @sa Fit
	 */
	double GetDesiredXmax() { return m_desiredXmax; }

	/** Returns the bottom-border layer coordinate that the user wants the
	 * mpWindow to show (it may be not exactly the actual shown coordinate in
	 * the case of locked aspect ratio).
	 * @sa Fit
	 */
	double GetDesiredYmin() { return m_desiredYmin; }

	/** Returns the top layer-border coordinate that the user wants the mpWindow
	 * to show (it may be not exactly the actual shown coordinate in the case of
	 * locked aspect ratio).
	 * @sa Fit
	 */
	double GetDesiredYmax() { return m_desiredYmax; }

	/** Returns the bounding box coordinates
		@param bbox Pointer to a 6-element double array where to store bounding
	   box coordinates. */
	void GetBoundingBox(double* bbox);

	/** Enable/disable scrollbars
	  @param status Set to true to show scrollbars */
	void SetMPScrollbars(bool status);

	/** Get scrollbars status.
	  @return true if scrollbars are visible */
	bool GetMPScrollbars() { return m_enableScrollBars; };

	/** Draw the window on a wxBitmap, then save it to a file.
	  @param filename File name where to save the screenshot
	  @param type image type to be saved: see wxImage output file types for
	  flags
	  @param imageSize Set a size for the output image. Default is the same as
	  the screen size
	  @param fit Decide whether to fit the plot into the size*/
	bool SaveScreenshot(
		const wxString& filename, int type = wxBITMAP_TYPE_BMP,
		wxSize imageSize = wxDefaultSize, bool fit = false);

	/** This value sets the zoom steps whenever the user clicks "Zoom in/out" or
	 * performs zoom with the mouse wheel. It must be a number above unity. This
	 * number is used for zoom in, and its inverse for zoom out. Set to 1.5 by
	 * default. */
	static double zoomIncrementalFactor;

	/** Set window margins, creating a blank area where some kinds of layers
	   cannot draw. This is useful for example to draw axes outside the area
	   where the plots are drawn.
		@param top Top border
		@param right Right border
		@param bottom Bottom border
		@param left Left border */
	void SetMargins(int top, int right, int bottom, int left);

	/** Set the top margin. @param top Top Margin */
	void SetMarginTop(int top) { m_marginTop = top; };
	/** Set the right margin. @param right Right Margin */
	void SetMarginRight(int right) { m_marginRight = right; };
	/** Set the bottom margin. @param bottom Bottom Margin */
	void SetMarginBottom(int bottom) { m_marginBottom = bottom; };
	/** Set the left margin. @param left Left Margin */
	void SetMarginLeft(int left) { m_marginLeft = left; };

	/** Get the top margin. @param top Top Margin */
	int GetMarginTop() { return m_marginTop; };
	/** Get the right margin. @param right Right Margin */
	int GetMarginRight() { return m_marginRight; };
	/** Get the bottom margin. @param bottom Bottom Margin */
	int GetMarginBottom() { return m_marginBottom; };
	/** Get the left margin. @param left Left Margin */
	int GetMarginLeft() { return m_marginLeft; };

	/** Sets whether to show coordinate tooltip when mouse passes over the plot.
	 * \param value true for enable, false for disable */
	// void EnableCoordTooltip(bool value = true);
	/** Gets coordinate tooltip status. \return true for enable, false for
	 * disable */
	// bool GetCoordTooltip() { return m_coordTooltip; };

	/** Check if a given point is inside the area of a mpInfoLayer and
	   eventually returns its pointer.
		@param point The position to be checked
		@return If an info layer is found, returns its pointer, nullptr
	   otherwise */
	mpInfoLayer* IsInsideInfoLayer(wxPoint& point);

	/** Sets the visibility of a layer by its name.
		@param name The layer name to set visibility
		@param viewable the view status to be set */
	void SetLayerVisible(const wxString& name, bool viewable);

	/** Check whether a layer with given name is visible
		@param name The layer name
		@return layer visibility status */
	bool IsLayerVisible(const wxString& name);

	/** Sets the visibility of a layer by its position in layer list.
		@param position The layer position in layer list
		@param viewable the view status to be set */
	void SetLayerVisible(const unsigned int position, bool viewable);

	/** Check whether the layer at given position is visible
		@param position The layer position in layer list
		@return layer visibility status */
	bool IsLayerVisible(const unsigned int position);

	/** Set Color theme. Provide colours to set a new colour theme.
		@param bgColour Background colour
		@param drawColour The colour used to draw all elements in foreground,
	   axes excluded
		@param axesColour The colour used to draw axes (but not their labels) */
	void SetColourTheme(
		const wxColour& bgColour, const wxColour& drawColour,
		const wxColour& axesColour);

	/** Get axes draw colour
		@return reference to axis colour used in theme */
	const wxColour& GetAxesColour() { return m_axColour; };

	/** Recalculate global layer bounding box, and save it in m_minX,...
	 * \return true if there is any valid BBox information.
	 */
	virtual bool UpdateBBox();

   protected:
	/** Paint handler, will plot all attached layers */
	void OnPaint(wxPaintEvent& event);
	/** Size handler, will update scroll bar sizes */
	void OnSize(wxSizeEvent& event);
	// void OnScroll2       (wxScrollWinEvent &event); //!< Scroll handler, will
	// move canvas
	/** Mouse handler, will show context menu */
	void OnShowPopupMenu(wxMouseEvent& event);
	/** Mouse handler, for detecting when the user drags with the right button
	 * or just "clicks" for the menu */
	void OnMouseRightDown(wxMouseEvent& event);
	/** Context menu handler */
	void OnCenter(wxCommandEvent& event);
	/** Context menu handler */
	void OnFit(wxCommandEvent& event);
	/** Context menu handler */
	void OnZoomIn(wxCommandEvent& event);
	/** Context menu handler */
	void OnZoomOut(wxCommandEvent& event);
	/** Context menu handler */
	void OnLockAspect(wxCommandEvent& event);
	/** Context menu handler */
	void OnMouseHelp(wxCommandEvent& event);
	/** Context menu handler */
	void OnPrintMenu(wxCommandEvent& event);
	/** Mouse handler for the wheel */
	void OnMouseWheel(wxMouseEvent& event);
	/** Mouse handler for mouse motion (for pan) */
	void OnMouseMove(wxMouseEvent& event);
	/** Mouse left click (for rect zoom) */
	void OnMouseLeftDown(wxMouseEvent& event);
	/** Mouse left click (for rect zoom) */
	void OnMouseLeftRelease(wxMouseEvent& event);
	/** Scroll thumb on scroll bar moving */
	void OnScrollThumbTrack(wxScrollWinEvent& event);
	/** Scroll page up */
	void OnScrollPageUp(wxScrollWinEvent& event);
	/** Scroll page down */
	void OnScrollPageDown(wxScrollWinEvent& event);
	/** Scroll line up */
	void OnScrollLineUp(wxScrollWinEvent& event);
	/** Scroll line down */
	void OnScrollLineDown(wxScrollWinEvent& event);
	/** Scroll to top */
	void OnScrollTop(wxScrollWinEvent& event);
	/** Scroll to bottom */
	void OnScrollBottom(wxScrollWinEvent& event);

	void DoScrollCalc(const int position, const int orientation);

	void DoZoomInXCalc(const int staticXpixel);
	void DoZoomInYCalc(const int staticYpixel);
	void DoZoomOutXCalc(const int staticXpixel);
	void DoZoomOutYCalc(const int staticYpixel);

	// wxList m_layers;    //!< List of attached plot layers
	/** List of attached plot layers */
	wxLayerList m_layers;
	/** Canvas' context menu */
	wxMenu m_popmenu;
	/** Scale aspect is locked or not */
	bool m_lockaspect;
	// bool   m_coordTooltip; //!< Selects whether to show coordinate tooltip
	/** Background Colour */
	wxColour m_bgColour;
	/** Foreground Colour */
	wxColour m_fgColour;
	/** Axes Colour */
	wxColour m_axColour;

	/** Global layer bounding box, left border incl. */
	double m_minX;
	/** Global layer bounding box, right border incl. */
	double m_maxX;
	/** Global layer bounding box, bottom border incl. */
	double m_minY;
	/** Global layer bounding box, top border incl. */
	double m_maxY;
	/** Current view's X scale */
	double m_scaleX;
	/** Current view's Y scale */
	double m_scaleY;
	/** Current view's X position */
	double m_posX;
	/** Current view's Y position */
	double m_posY;
	/** Current view's X dimension */
	int m_scrX;
	/** Current view's Y dimension */
	int m_scrY;
	/** Last mouse click X position, for centering and zooming the view */
	int m_clickedX = 0;
	/** Last mouse click Y position, for centering and zooming the view */
	int m_clickedY = 0;

	/** These are updated in Fit() only, and may be different from the real
	 * borders (layer coordinates) only if lock aspect ratio is true.
	 */
	double m_desiredXmin, m_desiredXmax, m_desiredYmin, m_desiredYmax;

	int m_marginTop, m_marginRight, m_marginBottom, m_marginLeft;

	/** For pan/zoom with the mouse. */
	bool m_enableMouseNavigation;
	bool m_mouseMovedAfterRightClick;
	/** For the right button "drag" feature */
	long m_mouseRClick_X = 0, m_mouseRClick_Y = 0;
	/** Starting coords for rectangular zoom selection */
	int m_mouseLClick_X = 0, m_mouseLClick_Y = 0;
	bool m_enableScrollBars;
	int m_scrollX = 0, m_scrollY = 0;
	/** For moving info layers over the window area */
	mpInfoLayer* m_movingInfoLayer;

	bool m_drawDottedSelectedWindow = false;
	wxPoint m_dottedWindowSize{0, 0};

	DECLARE_DYNAMIC_CLASS(mpWindow)
	DECLARE_EVENT_TABLE()
};

//-----------------------------------------------------------------------------
// mpFXYVector - provided by Jose Luis Blanco
//-----------------------------------------------------------------------------

/** A class providing graphs functionality for a 2D plot (either continuous or a
   set of points), from vectors of data. This class can be used directly, the
   user does not need to derive any new class. Simply pass the data as two
   vectors with the same length containing the X and Y coordinates to the method
   SetData.

	 To generate a graph with a set of points, call
	 \code
	 layerVar->SetContinuity(false)
	 \endcode

	 or

	 \code
	 layerVar->SetContinuity(true)
	 \endcode

	 to render the sequence of coordinates as a continuous line.

	 (Added: Jose Luis Blanco, AGO-2007)
*/
class WXDLLIMPEXP_MATHPLOT mpFXYVector : public mpFXY
{
   public:
	/** @param name  Label
		@param flags Label alignment, pass one of #mpALIGN_NE, #mpALIGN_NW,
	   #mpALIGN_SW, #mpALIGN_SE.
	*/
	mpFXYVector(wxString name = wxEmptyString, int flags = mpALIGN_NE);

	/** Changes the internal data: the set of points to draw.
		Both vectors MUST be of the same length. This method DOES NOT refresh
	  the mpWindow; do it manually.
	  * @sa Clear
	*/
	void SetData(const std::vector<double>& xs, const std::vector<double>& ys);

	/** Changes the internal data: the set of points to draw.
		Both vectors MUST be of the same length. This method DOES NOT refresh
	  the mpWindow; do it manually.
	  * @sa Clear
	*/
	void SetData(const std::vector<float>& xs, const std::vector<float>& ys);

	/** Clears all the data, leaving the layer empty.
	 * @sa SetData
	 */
	void Clear();

	/** Returns the number of data points currently hold in X & Y.
	 * @sa SetData
	 */
	size_t GetDataLength() const { return m_xs.size(); }

	/** Append a new data point (x,y)
	 * @sa SetData
	 */
	void AppendDataPoint(float x, float y);

   protected:
	/** The internal copy of the set of data to draw.
	 */
	std::vector<double> m_xs, m_ys;

	/** The internal counter for the "GetNextXY" interface
	 */
	size_t m_index;

	/** Loaded at SetData
	 */
	double m_minX, m_maxX, m_minY, m_maxY;

	/** Rewind value enumeration with mpFXY::GetNextXY.
		Overridden in this implementation.
	*/
	void Rewind() override;

	/** Get locus value for next N.
		Overridden in this implementation.
		@param x Returns X value
		@param y Returns Y value
	*/
	bool GetNextXY(double& x, double& y) override;

   public:
	/** Returns the actual minimum X data (loaded in SetData).
	 */
	double GetMinX() override { return m_minX; }

	/** Returns the actual minimum Y data (loaded in SetData).
	 */
	double GetMinY() override { return m_minY; }

	/** Returns the actual maximum X data (loaded in SetData).
	 */
	double GetMaxX() override { return m_maxX; }

	/** Returns the actual maximum Y data (loaded in SetData).
	 */
	double GetMaxY() override { return m_maxY; }

   protected:
	/** Holds label alignment */
	int m_flags = 0;

	DECLARE_DYNAMIC_CLASS(mpFXYVector)
};

//-----------------------------------------------------------------------------
// mpText - provided by Val Greene
//-----------------------------------------------------------------------------

/** Plot layer implementing a text string.
The text is plotted using a percentage system 0-100%, so the actual
coordinates for the location are not required, and the text stays
on the plot reguardless of the other layers location and scaling
factors.
*/
class WXDLLIMPEXP_MATHPLOT mpText : public mpLayer
{
   public:
	/** @param name text to be drawn in the plot
		@param offsetx holds offset for the X location in percentage (0-100)
		@param offsety holds offset for the Y location in percentage (0-100) */
	mpText(wxString name = wxT("Title"), int offsetx = 5, int offsety = 50);

	/** Text Layer plot handler.
		This implementation will plot text adjusted to the visible area. */
	void Plot(wxDC& dc, mpWindow& w) override;

	/** mpText should not be used for scaling decisions. */
	bool HasBBox() override { return FALSE; }

   protected:
	/** Holds offset for X in percentage */
	int m_offsetx;
	/** Holds offset for Y in percentage */
	int m_offsety;

	DECLARE_DYNAMIC_CLASS(mpText)
};

//-----------------------------------------------------------------------------
// mpPrintout - provided by Davide Rondini
//-----------------------------------------------------------------------------

/** Printout class used by mpWindow to draw in the objects to be printed.
	The object itself can then used by the default wxWidgets printing system
	to print mppWindow objects.
*/
class WXDLLIMPEXP_MATHPLOT mpPrintout : public wxPrintout
{
   public:
	mpPrintout(
		mpWindow* drawWindow,
		const wxChar* title = _T("wxMathPlot print output"));
	~mpPrintout() override = default;

	void SetDrawState(bool drawState) { drawn = drawState; };
	bool OnPrintPage(int page) override;
	bool HasPage(int page) override;

   private:
	bool drawn;
	mpWindow* plotWindow;
};

//-----------------------------------------------------------------------------
// mpMovableObject  - provided by Jose Luis Blanco
//-----------------------------------------------------------------------------
/** This virtual class represents objects that can be moved to an arbitrary 2D
 * location+rotation. The current transformation is set through
 * SetCoordinateBase. To ease the implementation of descendent classes,
 * mpMovableObject will be in charge of Bounding Box computation and layer
 * rendering, assuming that the object updates its shape in m_shape_xs &
 * m_shape_ys.
 */
class WXDLLIMPEXP_MATHPLOT mpMovableObject : public mpLayer
{
   public:
	/** Default constructor (sets location and rotation to (0,0,0))
	 */
	mpMovableObject() : m_shape_xs(0), m_shape_ys(0) { m_type = mpLAYER_PLOT; }

	~mpMovableObject() override = default;

	/** Get the current coordinate transformation.
	 */
	void GetCoordinateBase(double& x, double& y, double& phi) const
	{
		x = m_reference_x;
		y = m_reference_y;
		phi = m_reference_phi;
	}

	/** Set the coordinate transformation (phi in radians, 0 means no rotation).
	 */
	void SetCoordinateBase(double x, double y, double phi = 0)
	{
		m_reference_x = x;
		m_reference_y = y;
		m_reference_phi = phi;
		m_flags = mpALIGN_NE;
		ShapeUpdated();
	}

	bool HasBBox() override { return m_trans_shape_xs.size() != 0; }

	/** Get inclusive left border of bounding box.
	 */
	double GetMinX() override { return m_bbox_min_x; }

	/** Get inclusive right border of bounding box.
	 */
	double GetMaxX() override { return m_bbox_max_x; }

	/** Get inclusive bottom border of bounding box.
	 */
	double GetMinY() override { return m_bbox_min_y; }

	/** Get inclusive top border of bounding box.
	 */
	double GetMaxY() override { return m_bbox_max_y; }

	void Plot(wxDC& dc, mpWindow& w) override;

	/** Set label axis alignment.
	 *  @param align alignment (choose between mpALIGN_NE, mpALIGN_NW,
	 * mpALIGN_SW, mpALIGN_SE
	 */
	void SetAlign(int align) { m_flags = align; };

   protected:
	/** Holds label alignment */
	int m_flags = 0;

	/** The coordinates of the object (orientation "phi" is in radians).
	 */
	double m_reference_x{0}, m_reference_y{0}, m_reference_phi{0};

	/** A method for 2D translation and rotation, using the current
	 * transformation stored in m_reference_x,m_reference_y,m_reference_phi.
	 */
	void TranslatePoint(double x, double y, double& out_x, double& out_y);

	/** This contains the object points, in local coordinates (to be transformed
	 * by the current transformation).
	 */
	std::vector<double> m_shape_xs, m_shape_ys;

	/** The buffer for the translated & rotated points (to avoid recomputing
	 * them with each mpWindow refresh).
	 *
	 */
	std::vector<double> m_trans_shape_xs, m_trans_shape_ys;

	/** The precomputed bounding box:
	 * @sa ShapeUpdated
	 */
	double m_bbox_min_x, m_bbox_max_x, m_bbox_min_y, m_bbox_max_y;

	/** Must be called by the descendent class after updating the shape
	 * (m_shape_xs/ys), or when the transformation changes. This method updates
	 * the buffers m_trans_shape_xs/ys, and the precomputed bounding box.
	 */
	void ShapeUpdated();
};

//-----------------------------------------------------------------------------
// mpCovarianceEllipse  - provided by Jose Luis Blanco
//-----------------------------------------------------------------------------
/** A 2D ellipse, described by a 2x2 covariance matrix.
 *  The relation between the multivariate Gaussian confidence interval and
 *   the "quantiles" in this class is:
 *     - 1 : 68.27% confidence interval
 *     - 2 : 95.45%
 *     - 3 : 99.73%
 *     - 4 : 99.994%
 * For example, see
 * http://en.wikipedia.org/wiki/Normal_distribution#Standard_deviation_and_confidence_intervals
 *
 * The ellipse will be always centered at the origin. Use
 * mpMovableObject::SetCoordinateBase to move it.
 */
class WXDLLIMPEXP_MATHPLOT mpCovarianceEllipse : public mpMovableObject
{
   public:
	/** Default constructor.
	 * Initializes to a unity diagonal covariance matrix, a 95% confidence
	 * interval (2 sigmas), 32 segments, and a continuous plot
	 * (m_continuous=true).
	 */
	mpCovarianceEllipse(
		double cov_00 = 1, double cov_11 = 1, double cov_01 = 0,
		double quantiles = 2, int segments = 32,
		const wxString& layerName = wxT(""))
		: m_cov_00(cov_00),
		  m_cov_11(cov_11),
		  m_cov_01(cov_01),
		  m_quantiles(quantiles),
		  m_segments(segments)
	{
		m_continuous = true;
		m_name = layerName;
		RecalculateShape();
		m_type = mpLAYER_PLOT;
	}

	~mpCovarianceEllipse() override = default;

	double GetQuantiles() const { return m_quantiles; }

	/** Set how many "quantiles" to draw, that is, the confidence interval of
	 * the ellipse (see above).
	 */
	void SetQuantiles(double q)
	{
		m_quantiles = q;
		RecalculateShape();
	}

	void SetSegments(int segments) { m_segments = segments; }
	int GetSegments() const { return m_segments; }

	/** Returns the elements of the current covariance matrix:
	 */
	void GetCovarianceMatrix(
		double& cov_00, double& cov_01, double& cov_11) const
	{
		cov_00 = m_cov_00;
		cov_01 = m_cov_01;
		cov_11 = m_cov_11;
	}

	/** Changes the covariance matrix:
	 */
	void SetCovarianceMatrix(double cov_00, double cov_01, double cov_11)
	{
		m_cov_00 = cov_00;
		m_cov_01 = cov_01;
		m_cov_11 = cov_11;
		RecalculateShape();
	}

   protected:
	/** The elements of the matrix (only 3 since cov(0,1)=cov(1,0) in any
	 * positive definite matrix).
	 */
	double m_cov_00, m_cov_11, m_cov_01;
	double m_quantiles;

	/** The number of line segments that build up the ellipse.
	 */
	int m_segments;

	/** Called to update the m_shape_xs, m_shape_ys vectors, whenever a
	 * parameter changes.
	 */
	void RecalculateShape();
};

//-----------------------------------------------------------------------------
// mpPolygon - provided by Jose Luis Blanco
//-----------------------------------------------------------------------------
/** An arbitrary polygon, descendant of mpMovableObject.
 *  Use "setPoints" to set the list of N points. This class also can draw
 * non-closed polygons by passing the appropriate parameters to "setPoints". To
 * draw a point-cloud, call "SetContinuity(false)".
 */
class WXDLLIMPEXP_MATHPLOT mpPolygon : public mpMovableObject
{
   public:
	/** Default constructor.
	 */
	mpPolygon(const wxString& layerName = wxT(""))
	{
		m_continuous = true;
		m_name = layerName;
	}

	~mpPolygon() override = default;

	/** Set the points in the polygon.
	 * @param points_xs  The X coordinates of the points.
	 * @param points_ys  The Y coordinates of the points.
	 * @param closedShape If set to true, an additional segment will be added
	 * from the last to the first point.
	 */
	void setPoints(
		const std::vector<double>& points_xs,
		const std::vector<double>& points_ys, bool closedShape = true);

	/** Set the points in the polygon.
	 * @param points_xs  The X coordinates of the points.
	 * @param points_ys  The Y coordinates of the points.
	 * @param closedShape If set to true, an additional segment will be added
	 * from the last to the first point.
	 */
	void setPoints(
		const std::vector<float>& points_xs,
		const std::vector<float>& points_ys, bool closedShape = true);
};

//-----------------------------------------------------------------------------
// mpMovableObject  - provided by Jose Luis Blanco
//-----------------------------------------------------------------------------
/** This virtual class represents objects that can be moved to an arbitrary 2D
 * location+rotation. The current transformation is set through
 * SetCoordinateBase. To ease the implementation of descendent classes,
 * mpMovableObject will be in charge of Bounding Box computation and layer
 * render, assuming that the object updates its shape in m_shape_xs &
 * m_shape_ys.
 */
class WXDLLIMPEXP_MATHPLOT mpBitmapLayer : public mpLayer
{
   public:
	/** Default constructor.
	 */
	mpBitmapLayer()
	{
		m_min_x = m_max_x = m_min_y = m_max_y = 0;
		m_validImg = false;
		m_type = mpLAYER_BITMAP;
	}

	~mpBitmapLayer() override = default;

	/** Returns a copy of the current bitmap assigned to the layer.
	 */
	void GetBitmapCopy(wxImage& outBmp) const;

	/** Change the bitmap associated with the layer (to update the screen,
	 * refresh the mpWindow).
	 *  @param inBmp The bitmap to associate. A copy is made, thus it can be
	 * released after calling this.
	 *  @param x The left corner X coordinate (in plot units).
	 *  @param y The top corner Y coordinate (in plot units).
	 *  @param lx The width in plot units.
	 *  @param ly The height in plot units.
	 */
	void SetBitmap(
		const wxImage& inBmp, double x, double y, double lx, double ly);

	bool HasBBox() override { return true; }

	/** Get inclusive left border of bounding box.
	 */
	double GetMinX() override { return m_min_x; }

	/** Get inclusive right border of bounding box.
	 */
	double GetMaxX() override { return m_max_x; }

	/** Get inclusive bottom border of bounding box.
	 */
	double GetMinY() override { return m_min_y; }

	/** Get inclusive top border of bounding box.
	 */
	double GetMaxY() override { return m_max_y; }

	void Plot(wxDC& dc, mpWindow& w) override;

	/** Set label axis alignment.
	 *  @param align alignment (choose between mpALIGN_NE, mpALIGN_NW,
	 * mpALIGN_SW, mpALIGN_SE
	 */
	void SetAlign(int align) { m_flags = align; };

   protected:
	/** Holds label alignment */
	int m_flags = 0;

	/** The internal copy of the Bitmap:
	 */
	wxImage m_bitmap;
	wxBitmap m_scaledBitmap;
	wxCoord m_scaledBitmap_offset_x, m_scaledBitmap_offset_y;

	bool m_validImg;

	/** The shape of the bitmap:
	 */
	double m_min_x, m_max_x, m_min_y, m_max_y;
};

/*@}*/

#if defined(_MSC_VER)
#pragma warning(pop)
#endif
