Zone2.h

// Copyright (C) Geom Software e.U, Bernhard Kornberger, Graz/Austria
//
// This file is part of the Fade2D library. The student license is free
// of charge and covers personal non-commercial research. Licensees
// holding a commercial license may use this file in accordance with
// the Commercial License Agreement.
//
// This software is provided AS IS with NO WARRANTY OF ANY KIND,
// INCLUDING THE WARRANTY OF DESIGN, MERCHANTABILITY AND FITNESS
// FOR A PARTICULAR PURPOSE.
//
// Please contact the author if any conditions of this licensing are
// not clear to you.
//
// Author: Bernhard Kornberger, bkorn (at) geom.at
// http://www.geom.at

#pragma once

#include "common.h"
#include "Bbox2.h"
#include "Edge2.h"
#include "Segment2.h"
#include "UserPredicates.h"

#if GEOM_PSEUDO3D==GEOM_TRUE
    namespace GEOM_FADE25D {
#elif GEOM_PSEUDO3D==GEOM_FALSE
    namespace GEOM_FADE2D {
#else
    #error GEOM_PSEUDO3D is not defined
#endif


class Dt2; // Fwd
class ConstraintGraph2; // Fwd
class Triangle2; // Fwd
class Point2; // Fwd
class Visualizer2; // Fwd
class CLASS_DECLSPEC Zone2
{
public:
    Zone2(Dt2* pDt_,ZoneLocation zoneLoc_);
    Zone2(Dt2* pDt_,ZoneLocation zoneLoc_,ConstraintGraph2* pConstraintGraph_);
    Zone2(Dt2* pDt_,ZoneLocation zoneLoc_,const std::vector<ConstraintGraph2*>& vConstraintGraphs_);

    Zone2(Dt2* pDt_,const std::vector<ConstraintGraph2*>& vConstraintGraphs_,ZoneLocation zoneLoc_,std::vector<Point2>& vStartPoints);

    ZoneLocation getZoneLocation() const;

    Zone2* convertToBoundedZone();

    void show(const std::string& postscriptFilename,bool bShowFull,bool bWithConstraints) const;

    void show(Visualizer2* pVisualizer,bool bShowFull,bool bWithConstraints) const;

#if GEOM_PSEUDO3D==GEOM_TRUE

    void showGeomview(const std::string& filename,const std::string& color) const;
#endif

#if GEOM_PSEUDO3D==GEOM_TRUE


void analyzeAngles();
void slopeValleyRidgeOptimization();

    void optimizeValleysAndRidges(double tolerance2D,double lowerThreshold25D);
#endif

void unifyGrid(double tolerance);

    bool assignDt2(Dt2* pDt_);
    void getTriangles(std::vector<Triangle2*>& vTriangles_) const;

    void getVertices(std::vector<Point2*>& vVertices_) const;


    ConstraintGraph2* getConstraintGraph() const;

    size_t getNumberOfTriangles() const;


    void getConstraintGraphs(std::vector<ConstraintGraph2*>& vConstraintGraphs_) const;

    Dt2* getDelaunayTriangulation() const;

    size_t numberOfConstraintGraphs() const;

    void debug(std::string name="");

    Bbox2 getBoundingBox() const;

    void getBoundaryEdges(std::vector<Edge2>& vEdges) const;

    void getBoundarySegments(std::vector<Segment2>& vSegments) const;

    double getArea2D() const;


#if GEOM_PSEUDO3D==GEOM_TRUE

    double getArea25D() const;
#endif

    void getBorderEdges(std::vector<Edge2>& vBorderEdgesOut) const;


protected:
    Zone2(const Zone2&);
    void getTriangles_RESULT(std::vector<Triangle2*>& vTriangles) const;
    void initWorkspace(bool bInside,std::vector<Triangle2*>& vWorkspace) const;
    void bfsFromWorkspace(std::vector<Triangle2*>& vWorkspace,std::vector<Triangle2*>& vTriangles) const;

    Zone2* ctbz_treatCC(std::vector<Triangle2*>& vOneCC);
    void optimize(std::vector<Triangle2*>& vT,std::vector<Triangle2*>& vChangedT);
    // Data
    Dt2* pDt;
    ZoneLocation zoneLoc;
    //Zone2* pInputZone0;
    //Zone2* pInputZone1;
    CLASS_DECLSPEC
    friend Zone2* zoneUnion(Zone2* pZone0,Zone2* pZone1);
    CLASS_DECLSPEC
    friend Zone2* zoneIntersection(Zone2* pZone0,Zone2* pZone1);
    CLASS_DECLSPEC
    friend Zone2* zoneDifference(Zone2* pZone0,Zone2* pZone1);
    CLASS_DECLSPEC
    friend Zone2* zoneSymmetricDifference(Zone2* pZone0,Zone2* pZone1);
    CLASS_DECLSPEC
    friend Zone2* peelOffIf(Zone2* pZone, UserPredicateT* pPredicate,bool bVerbose);

private:
#ifdef _WIN32
#pragma warning(push)
#pragma warning(disable:4251)
#endif
    std::vector<ConstraintGraph2*> vConstraintGraphs;
    std::vector<Point2> vStartPoints;
    std::vector<Zone2*> vInputZones;
#ifdef _WIN32
#pragma warning(pop)
#endif

};

// Free functions

CLASS_DECLSPEC
Zone2* zoneUnion(Zone2* pZone0,Zone2* pZone1);
CLASS_DECLSPEC
Zone2* zoneIntersection(Zone2* pZone0,Zone2* pZone1);
CLASS_DECLSPEC
Zone2* zoneDifference(Zone2* pZone0,Zone2* pZone1);
CLASS_DECLSPEC
Zone2* zoneSymmetricDifference(Zone2* pZone0,Zone2* pZone1);





} // (namespace)