diff options
Diffstat (limited to 'src/qhull/src/libqhullcpp/Qhull.h')
| -rw-r--r-- | src/qhull/src/libqhullcpp/Qhull.h | 132 |
1 files changed, 132 insertions, 0 deletions
diff --git a/src/qhull/src/libqhullcpp/Qhull.h b/src/qhull/src/libqhullcpp/Qhull.h new file mode 100644 index 000000000..e3ad9d8e1 --- /dev/null +++ b/src/qhull/src/libqhullcpp/Qhull.h @@ -0,0 +1,132 @@ +/**************************************************************************** +** +** Copyright (c) 2008-2015 C.B. Barber. All rights reserved. +** $Id: //main/2015/qhull/src/libqhullcpp/Qhull.h#3 $$Change: 2066 $ +** $DateTime: 2016/01/18 19:29:17 $$Author: bbarber $ +** +****************************************************************************/ + +#ifndef QHULLCPP_H +#define QHULLCPP_H + +#include "libqhullcpp/QhullPoint.h" +#include "libqhullcpp/QhullVertex.h" +#include "libqhullcpp/QhullFacet.h" + +namespace orgQhull { + +/*** + Compile qhullcpp and libqhull with the same compiler. setjmp() and longjmp() must be the same. + + #define QHULL_NO_STL + Do not supply conversions to STL + Coordinates.h requires <vector>. It could be rewritten for another vector class such as QList + #define QHULL_USES_QT + Supply conversions to QT + qhulltest requires QT. It is defined in RoadTest.h + + #define QHULL_ASSERT + Defined by QhullError.h + It invokes assert() +*/ + +#//!\name Used here + class QhullFacetList; + class QhullPoints; + class QhullQh; + class RboxPoints; + +#//!\name Defined here + class Qhull; + +//! Interface to Qhull from C++ +class Qhull { + +private: +#//!\name Members and friends + QhullQh * qh_qh; //! qhT for this instance + Coordinates origin_point; //! origin for qh_qh->hull_dim. Set by runQhull() + bool run_called; //! True at start of runQhull. Errors if call again. + Coordinates feasible_point; //! feasible point for half-space intersection (alternative to qh.feasible_string for qh.feasible_point) + +public: +#//!\name Constructors + Qhull(); //!< call runQhull() next + Qhull(const RboxPoints &rboxPoints, const char *qhullCommand2); + Qhull(const char *inputComment2, int pointDimension, int pointCount, const realT *pointCoordinates, const char *qhullCommand2); + ~Qhull() throw(); +private: //! Disable copy constructor and assignment. Qhull owns QhullQh. + Qhull(const Qhull &); + Qhull & operator=(const Qhull &); + +private: + void allocateQhullQh(); + +public: + +#//!\name GetSet + void checkIfQhullInitialized(); + int dimension() const { return qh_qh->input_dim; } //!< Dimension of input and result + void disableOutputStream() { qh_qh->disableOutputStream(); } + void enableOutputStream() { qh_qh->enableOutputStream(); } + countT facetCount() const { return qh_qh->num_facets; } + Coordinates feasiblePoint() const; + int hullDimension() const { return qh_qh->hull_dim; } //!< Dimension of the computed hull + bool hasOutputStream() const { return qh_qh->hasOutputStream(); } + bool initialized() const { return (qh_qh->hull_dim>0); } + const char * inputComment() const { return qh_qh->rbox_command; } + QhullPoint inputOrigin(); + //! non-const due to QhullPoint + QhullPoint origin() { QHULL_ASSERT(initialized()); return QhullPoint(qh_qh, origin_point.data()); } + QhullQh * qh() const { return qh_qh; }; + const char * qhullCommand() const { return qh_qh->qhull_command; } + const char * rboxCommand() const { return qh_qh->rbox_command; } + int rotateRandom() const { return qh_qh->ROTATErandom; } //!< Return QRn for repeating QR0 runs + void setFeasiblePoint(const Coordinates &c) { feasible_point= c; } //!< Sets qh.feasible_point via initializeFeasiblePoint + countT vertexCount() const { return qh_qh->num_vertices; } + +#//!\name Delegated to QhullQh + double angleEpsilon() const { return qh_qh->angleEpsilon(); } //!< Epsilon for hyperplane angle equality + void appendQhullMessage(const std::string &s) { qh_qh->appendQhullMessage(s); } + void clearQhullMessage() { qh_qh->clearQhullMessage(); } + double distanceEpsilon() const { return qh_qh->distanceEpsilon(); } //!< Epsilon for distance to hyperplane + double factorEpsilon() const { return qh_qh->factorEpsilon(); } //!< Factor for angleEpsilon and distanceEpsilon + std::string qhullMessage() const { return qh_qh->qhullMessage(); } + bool hasQhullMessage() const { return qh_qh->hasQhullMessage(); } + int qhullStatus() const { return qh_qh->qhullStatus(); } + void setErrorStream(std::ostream *os) { qh_qh->setErrorStream(os); } + void setFactorEpsilon(double a) { qh_qh->setFactorEpsilon(a); } + void setOutputStream(std::ostream *os) { qh_qh->setOutputStream(os); } + +#//!\name ForEach + QhullFacet beginFacet() const { return QhullFacet(qh_qh, qh_qh->facet_list); } + QhullVertex beginVertex() const { return QhullVertex(qh_qh, qh_qh->vertex_list); } + void defineVertexNeighborFacets(); //!< Automatically called if merging facets or Voronoi diagram + QhullFacet endFacet() const { return QhullFacet(qh_qh, qh_qh->facet_tail); } + QhullVertex endVertex() const { return QhullVertex(qh_qh, qh_qh->vertex_tail); } + QhullFacetList facetList() const; + QhullFacet firstFacet() const { return beginFacet(); } + QhullVertex firstVertex() const { return beginVertex(); } + QhullPoints points() const; + QhullPointSet otherPoints() const; + //! Same as points().coordinates() + coordT * pointCoordinateBegin() const { return qh_qh->first_point; } + coordT * pointCoordinateEnd() const { return qh_qh->first_point + qh_qh->num_points*qh_qh->hull_dim; } + QhullVertexList vertexList() const; + +#//!\name Methods + double area(); + void outputQhull(); + void outputQhull(const char * outputflags); + void runQhull(const RboxPoints &rboxPoints, const char *qhullCommand2); + void runQhull(const char *inputComment2, int pointDimension, int pointCount, const realT *pointCoordinates, const char *qhullCommand2); + double volume(); + +#//!\name Helpers +private: + void initializeFeasiblePoint(int hulldim); +};//Qhull + +}//namespace orgQhull + +#endif // QHULLCPP_H |