/**
 * Copyright (c) 2011-2018 by Andrew Mustun. All rights reserved.
 * 
 * This file is part of the QCAD project.
 *
 * QCAD 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 3 of the License, or
 * (at your option) any later version.
 *
 * QCAD is distributed in the hope that it will be useful,
 * but WITHOUT ANY WARRANTY; without even the implied warranty of
 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
 * GNU General Public License for more details.
 *
 * You should have received a copy of the GNU General Public License
 * along with QCAD.
 */

#ifndef RSHAPE_H
#define RSHAPE_H

#include "../core_global.h"

#include <QList>
#include <QSharedPointer>

#include "RMath.h"
#include "RShapeProxy.h"
#include "RVector.h"

class RArc;
class RBox;
class RCircle;
class REllipse;
class RExplodable;
class RLine;
class RPolyline;
class RSpline;
class RTriangle;

#ifndef RDEFAULT_TOLERANCE_1E_MIN4
#define RDEFAULT_TOLERANCE_1E_MIN4 1.0e-4
#endif

class RShapeTransformation {
public:
    virtual ~RShapeTransformation() {}
    virtual RVector transform(const RVector& v) const = 0;
};

class RShapeTransformationScale : public RShapeTransformation {
public:
    RShapeTransformationScale(const RVector& factors, const RVector& center) : factors(factors), center(center) {}
    virtual ~RShapeTransformationScale() {}

    virtual RVector transform(const RVector& v) const {
        return v.getScaled(factors, center);
    }

    RVector factors;
    RVector center;
};

/**
 * Interface for geometrical shape classes.
 *
 * \ingroup math
 * \scriptable
 * \sharedPointerSupport
 */
class QCADCORE_EXPORT RShape {
public:
    enum Type {
        Unknown = -1,
        Point = 0,
        Line = 1,
        Arc = 2,
        Circle = 3,
        Ellipse = 4,
        Polyline = 5,
        Spline = 6,
        Triangle = 7,
        XLine = 8,
        Ray = 9
    };

public:
    RShape();
    virtual ~RShape();

    virtual bool isValid() const {
        return true;
    }

    virtual RShape::Type getShapeType() const {
        return Unknown;
    }

    static bool isPointShape(const RShape& s) {
        return s.getShapeType()==RShape::Point;
    }
    static bool isLineShape(const RShape& s) {
        return s.getShapeType()==RShape::Line;
    }
    static bool isArcShape(const RShape& s) {
        return s.getShapeType()==RShape::Arc;
    }
    static bool isCircleShape(const RShape& s) {
        return s.getShapeType()==RShape::Circle;
    }
    static bool isEllipseShape(const RShape& s) {
        return s.getShapeType()==RShape::Ellipse;
    }
    static bool isFullEllipseShape(const RShape& s);
    static bool isPolylineShape(const RShape& s) {
        return s.getShapeType()==RShape::Polyline;
    }
    static bool isSplineShape(const RShape& s) {
        return s.getShapeType()==RShape::Spline;
    }
    static bool isTriangleShape(const RShape& s) {
        return s.getShapeType()==RShape::Triangle;
    }
    static bool isXLineShape(const RShape& s) {
        return s.getShapeType()==RShape::XLine;
    }
    static bool isRayShape(const RShape& s) {
        return s.getShapeType()==RShape::Ray;
    }

    virtual QSharedPointer<RShape> clone() const = 0;

    // virtual QSharedPointer<RShape> cloneToSharedPointer() const {
    //     return QSharedPointer<RShape>(clone());
    // }

    /**
     * Reimplement and return true if this shape relies on interpolation
     * for geometrical operations such as finding intersection points. E.g.
     * splines.
     */
    virtual bool isInterpolated() const {
        return false;
    }

    virtual RVector getClosestPointOnShape(const RVector& p,
            bool limited = true, double strictRange = RMAXDOUBLE) const;

    /**
     * \return The bounding box of this shape.
     */
    virtual RBox getBoundingBox() const = 0;

    /**
     * Flattens this shape to the X/Y plane
     */
    virtual void to2D() {
        setZ(0.0);
    }

    virtual void setZ(double z) = 0;

    /**
     * \return The total length of this shape if applicable or NaN if not.
     */
    virtual double getLength() const = 0;

    virtual bool equals(const RShape& other, double tolerance=RS::PointTolerance) const;
    virtual QList<RVector> getVectorProperties() const { return QList<RVector>(); }
    virtual QList<int> getIntProperties() const { return QList<int>(); }
    virtual QList<double> getDoubleProperties() const { return QList<double>(); }
    virtual QList<bool> getBoolProperties() const { return QList<bool>(); }

    /**
     * \return Shortest vector from this shape to the given point.
     *
     * \param point the given point
     *
     * \param limited If false, the shape is treated as an unlimited
     *      shape that extends in a logical way. For example an
     *      unlimited line for a line, full circle for an arc or
     *      a plane for a triangle.
     */
    virtual RVector getVectorTo(const RVector& point,
            bool limited = true, double strictRange = RMAXDOUBLE) const = 0;

    virtual double getDistanceTo(const RVector& point, bool limited = true, double strictRange = RMAXDOUBLE) const;
    virtual double getMaxDistanceTo(const QList<RVector>& points, bool limited = true, double strictRange = RMAXDOUBLE) const;
    virtual bool isOnShape(const RVector& point,
                           bool limited = true,
                           double tolerance = RDEFAULT_TOLERANCE_1E_MIN4) const;
    virtual QList<RVector> filterOnShape(const QList<RVector>& pointList,
        bool limited = true,
        double tolerance = RDEFAULT_TOLERANCE_1E_MIN4) const;
    virtual RVector getVectorFromEndpointTo(const RVector& point) const;

    /**
     * \return The end point(s) of this shape.
     */
    virtual QList<RVector> getEndPoints() const = 0;

    /**
     * \return The middle point(s) of this shape.
     */
    virtual QList<RVector> getMiddlePoints() const = 0;

    /**
     * \return The center point(s) of this shape.
     */
    virtual QList<RVector> getCenterPoints() const = 0;

    /**
     * \return The reference point(s) of this shape.
     */
    virtual QList<RVector> getArcReferencePoints() const {
        return QList<RVector>();
    }

    virtual RVector getPointOnShape() const;

    virtual QList<RVector> getPointCloud(double segmentLength) const = 0;

    /**
     * \return All points on this shape with the given distance to an endpoint.
     */
    virtual QList<RVector> getPointsWithDistanceToEnd(
        double distance, int from = RS::FromAny) const = 0;

    virtual RVector getPointWithDistanceToStart(double distance) const {
        QList<RVector> res = getPointsWithDistanceToEnd(distance, RS::FromStart|RS::AlongPolyline);
        if (res.isEmpty()) {
            return RVector::invalid;
        }
        return res[0];
    }

    virtual RVector getPointWithDistanceToEnd(double distance) const {
        QList<RVector> res = getPointsWithDistanceToEnd(distance, RS::FromEnd|RS::AlongPolyline);
        if (res.isEmpty()) {
            return RVector::invalid;
        }
        return res[0];
    }

    /**
     * \return Angle on the entity at the given distance from the start point.
     */
    virtual double getAngleAt(double distance, RS::From from = RS::FromStart) const {
        Q_UNUSED(distance)
        Q_UNUSED(from)
        return RNANDOUBLE;
    }
    virtual double getAngleAtPoint(const RVector& pos) const {
        double d = getDistanceFromStart(pos);
        return getAngleAt(d);
    }

    virtual RVector getPointAtPercent(double p) const;
    virtual double getAngleAtPercent(double p) const;

    virtual bool intersectsWith(const RShape& other, 
        bool limited = true) const;

    QList<RVector> getIntersectionPoints(const RShape& other,
            bool limited = true, bool same = false, bool force = false) const;

    virtual QList<RVector> getSelfIntersectionPoints(double tolerance=RS::PointTolerance) const {
        Q_UNUSED(tolerance)
        return QList<RVector>();
    }

    virtual bool isDirected() const {
        return false;
    }

    virtual double getDirection1() const {
        return RNANDOUBLE;
    }
    virtual double getDirection2() const {
        return RNANDOUBLE;
    }

    virtual RS::Side getSideOfPoint(const RVector& point) const {
        Q_UNUSED(point)
        return RS::NoSide;
    }

    virtual RVector getStartPoint() const {
        return RVector::invalid;
    }
    virtual RVector getEndPoint() const {
        return RVector::invalid;
    }
    virtual RVector getMiddlePoint() const {
        return RVector::invalid;
    }

    virtual bool reverse() {
        return false;
    }

    /**
     * \param trimPoint Trim start point or end point to this coordinate.
     * \param clickPoint Point that was clicked by user.
     * \param extend True if the shape is to be extended.
     * Used for polylines to determine that the first segment has to be trimmed even if another segment was clicked.
     */
    virtual bool trimStartPoint(const RVector& trimPoint, const RVector& clickPoint = RVector::invalid, bool extend = false) {
        Q_UNUSED(trimPoint)
        Q_UNUSED(clickPoint)
        Q_UNUSED(extend)

        return false;
    }

    virtual bool trimStartPoint(double trimDist) {
        RVector p = getPointWithDistanceToStart(trimDist);
        return trimStartPoint(p);
    }

    /**
     * \param trimPoint Trim start point or end point to this coordinate.
     * \param clickPoint Point that was clicked by user.
     * \param extend True if the shape is to be extended.
     * Used for polylines to determine that the last segment has to be trimmed even if another segment was clicked.
     */
    virtual bool trimEndPoint(const RVector& trimPoint, const RVector& clickPoint = RVector::invalid, bool extend = false) {
        Q_UNUSED(trimPoint)
        Q_UNUSED(clickPoint)
        Q_UNUSED(extend)

        return false;
    }

    virtual bool trimEndPoint(double trimDist) {
        RVector p = getPointWithDistanceToStart(trimDist);
        return trimEndPoint(p);
    }

    /**
     * \param trimPoint Trim start point or end point to this coordinate.
     * \param clickPoint Point that was clicked by user.
     */
    virtual RS::Ending getTrimEnd(const RVector& trimPoint, const RVector& clickPoint) {
        Q_UNUSED(trimPoint)
        Q_UNUSED(clickPoint)

        return RS::EndingNone;
    }

    /**
     * \return Distance of given point to start point of shape along shape.
     * \param p Point (assumed to be on shape).
     */
    virtual double getDistanceFromStart(const RVector& p) const {
        Q_UNUSED(p)
        return RMAXDOUBLE;
    }

    /**
     * \return All possible distances of given point to start point of shape along shape.
     * \param p Point (assumed to be on shape).
     */
    virtual QList<double> getDistancesFromStart(const RVector& p) const {
        return QList<double>() << getDistanceFromStart(p);
    }

    /**
     * \return The intersection point(s) between this shape and the given
     *      other shape.
     */
    static QList<RVector> getIntersectionPoints(const RShape& shape1,
            const RShape& shape2, bool limited = true, bool same = false, bool force = false);

    virtual bool move(const RVector& offset) = 0;
    virtual bool rotate(double rotation, const RVector& center = RDEFAULT_RVECTOR) = 0;
    virtual bool scale(double scaleFactor, const RVector& center = RVector());
    virtual bool scale(const RVector& scaleFactors, const RVector& center = RVector()) = 0;
    virtual bool mirror(const RLine& axis) = 0;
    virtual bool flipHorizontal();
    virtual bool flipVertical();
    virtual bool stretch(const RBox& area, const RVector& offset);
    virtual bool stretch(const RPolyline& area, const RVector& offset);

    virtual QSharedPointer<RShape> getTransformed(const QTransform& transform) const = 0;

    static QList<RPolyline> getPolylines(const QList<QSharedPointer<RShape> >& shapes);
    static QList<QSharedPointer<RShape> > getOrderedShapes(const QList<QSharedPointer<RShape> >& shapes);
    static bool order(QList<QList<QSharedPointer<RShape> > >& boundary);

    static QList<RVector> getIntersectionPointsLL(const RLine& line1,
            const RLine& line2, bool limited = true) {
        return getIntersectionPointsLL(line1, line2, limited, limited);
    }
    static QList<RVector> getIntersectionPointsLL(const RLine& line1,
            const RLine& line2, bool limited1, bool limited2);
    static QList<RVector> getIntersectionPointsLA(const RLine& line1,
           const RArc& arc2, bool limited = true) {
        return getIntersectionPointsLA(line1, arc2, limited, limited);
    }
    static QList<RVector> getIntersectionPointsLA(const RLine& line1,
            const RArc& arc2, bool limited1, bool limited2);
    static QList<RVector> getIntersectionPointsLC(const RLine& line1,
            const RCircle& circle2, bool limited = true);
    static QList<RVector> getIntersectionPointsLE(const RLine& line1,
            const REllipse& ellipse2, bool limited = true) {
        return getIntersectionPointsLE(line1, ellipse2, limited, limited);
    }
    static QList<RVector> getIntersectionPointsLE(const RLine& line1,
            const REllipse& ellipse2, bool limited1, bool limited2);
    static QList<RVector> getIntersectionPointsLT(const RLine& line1,
            const RTriangle& triangle2, bool limited = true) {
        return getIntersectionPointsLT(line1, triangle2, limited, limited);
    }
    static QList<RVector> getIntersectionPointsLT(const RLine& line1,
            const RTriangle& triangle2, bool limited1, bool limited2);
    static QList<RVector> getIntersectionPointsLS(const RLine& line1,
            const RSpline& spline2, bool limited = true);
    static QList<RVector> getIntersectionPointsLX(const RLine& line1,
            const RExplodable& explodable2, bool limited = true);

    static QList<RVector> getIntersectionPointsAA(const RArc& arc1,
            const RArc& arc2, bool limited = true);
    static QList<RVector> getIntersectionPointsAC(const RArc& arc1,
            const RCircle& circle2, bool limited = true);
    static QList<RVector> getIntersectionPointsAE(const RArc& arc1,
            const REllipse& ellipse2, bool limited = true);
    static QList<RVector> getIntersectionPointsAT(const RArc& arc1,
            const RTriangle& triangle2, bool limited = true);
    static QList<RVector> getIntersectionPointsAS(const RArc& arc1,
            const RSpline& spline2, bool limited = true);
    static QList<RVector> getIntersectionPointsAX(const RArc& arc1,
            const RExplodable& explodable2, bool limited = true);

    static QList<RVector> getIntersectionPointsCC(const RCircle& circle1,
            const RCircle& circle2);
    static QList<RVector> getIntersectionPointsCE(const RCircle& circle1,
            const REllipse& ellipse2);
    static QList<RVector> getIntersectionPointsCS(const RCircle& circle1,
            const RSpline& spline2, bool limited = true);
    static QList<RVector> getIntersectionPointsCX(const RCircle& circle1,
            const RExplodable& explodable2, bool limited = true);

    static QList<RVector> getIntersectionPointsEE(const REllipse& ellipse1,
            const REllipse& ellipse2);
    static QList<RVector> getIntersectionPointsEE(const REllipse& ellipse1,
            const REllipse& ellipse2, bool limited);
    static QList<RVector> getIntersectionPointsES(const REllipse& ellipse1,
            const RSpline& spline2, bool limited = true);
    static QList<RVector> getIntersectionPointsEX(const REllipse& ellipse1,
            const RExplodable& explodable2, bool limited = true);

    static QList<RVector> getIntersectionPointsSX(const RSpline& spline1,
            const RExplodable& explodable2, bool limited);

    static QList<RVector> getIntersectionPointsSS(const RSpline& spline1,
            const RSpline& spline2, bool limited = true, bool same = false, double tolerance=RS::PointTolerance);

    static QList<RVector> getIntersectionPointsXX(const RExplodable& explodable1,
            const RExplodable& explodable2, bool limited = true, bool same = false);

    static const RExplodable* castToExplodable(const RShape* shape);

    virtual QList<QSharedPointer<RShape> > getOffsetShapes(double distance, int number, RS::Side side, const RVector& position = RVector::invalid);

    static QList<QSharedPointer<RShape> > getOffsetLines(const RShape& shape, double distance, int number, RS::Side side, const RVector& position = RVector::invalid);
    static QList<QSharedPointer<RShape> > getOffsetArcs(const RShape& shape, double distance, int number, RS::Side side, const RVector& position = RVector::invalid);

    static QList<QSharedPointer<RShape> > getReversedShapeList(const QList<QSharedPointer<RShape> >& shapes);

    virtual QList<QSharedPointer<RShape> > splitAt(const QList<RVector>& points) const;

    static QList<QSharedPointer<RShape> > trim(
            const RShape& trimShape, const RVector& trimClickPos,
            const RShape& limitingShape, const RVector& limitingClickPos,
            bool trimBoth, bool samePolyline);

    static QList<QSharedPointer<RShape> > roundCorners(const QList<QSharedPointer<RShape> >& shapes, double radius);

    static QList<QSharedPointer<RShape> > roundShapes(
            const QSharedPointer<RShape> shape1, const RVector& clickPos1,
            const QSharedPointer<RShape> shape2, const RVector& clickPos2,
            bool trim, bool samePolyline, double radius, const RVector& pos);

    static QSharedPointer<RShape> xLineToRay(QSharedPointer<RShape> shape);
    static QSharedPointer<RShape> rayToLine(QSharedPointer<RShape> shape);

    static QSharedPointer<RShape> scaleArc(const RShape& shape, const RVector& scaleFactors, const RVector& center = RDEFAULT_RVECTOR) {
        RShapeTransformationScale t(scaleFactors, center);
        return transformArc(shape, t);
    }

    /**
     * \nonscriptable
     */
    static QSharedPointer<RShape> transformArc(const RShape& shape, RShapeTransformation& transformation);
    static QSharedPointer<RShape> ellipseToArcCircleEllipse(const REllipse& ellipse);

    static int getErrorCode() {
        return errorCode;
    }

    void dump() const;

    /**
     * \nonscriptable
     */
    friend QDebug operator<<(QDebug dbg, const RShape& s) {
        s.print(dbg);
        return dbg;
    }

    static bool hasProxy() {
        return shapeProxy!=NULL;
    }

    /**
     * \nonscriptable
     */
    static void setShapeProxy(RShapeProxy* p) {
        if (shapeProxy!=NULL) {
            delete shapeProxy;
        }
        shapeProxy = p;
    }

    /**
     * \nonscriptable
     */
    static RShapeProxy* getShapeProxy() {
        return shapeProxy;
    }

private:
    static double ellipse2tr(double x, double y, double AA, double BB,
                        double CC, double DD, double EE, double FF);

private:
    static double twopi;
    static double epsTolerance;

protected:
    virtual void print(QDebug dbg) const;
    static int errorCode;

private:
    static RShapeProxy* shapeProxy;
};

Q_DECLARE_METATYPE(RShape*)
Q_DECLARE_METATYPE(RShape::Type)
Q_DECLARE_METATYPE(const RShape*)
Q_DECLARE_METATYPE(QSharedPointer<RShape>)
Q_DECLARE_METATYPE(QSharedPointer<const RShape>)
Q_DECLARE_METATYPE(QSharedPointer<RShape>*)
Q_DECLARE_METATYPE(QList<QSharedPointer<RShape> >)
Q_DECLARE_METATYPE(QList<QSharedPointer<RShape> >*)
Q_DECLARE_METATYPE(QList<QSharedPointer<const RShape> >)

#endif
