qhullwrapper.h

#ifndef QHULLWRAPPER_H_PROTECTOR
#define QHULLWRAPPER_H_PROTECTOR

#include <iostream>
#include <vector>
//#include <conio.h>

extern "C"
{

#include <stdio.h>
#include <stdlib.h>
#include "qhull/qhull.h"
#include "qhull/mem.h"
#include "qhull/qset.h"
#include "qhull/geom.h"
#include "qhull/merge.h"
#include "qhull/poly.h"
#include "qhull/io.h"
#include "qhull/stat.h"

}

#include "element.h"
#include "afltypes.h"
#include "util/aflutil.h"
#include "ddf/ddf.h"

#include <set>
#include <vector>

namespace afl 
{
    template <size_t D>
    class Locator
    {
        public:
            Locator( size_t n, const ddf::Node* min, const ddf::Node* max ) { throw exc::Exception("Locator is only implemented from 2D to 4D? This is likely a bug in AFL"); }

            size_t locate( const ddf::Node& x ) const { throw exc::Exception("Locator is only implemented from 2D to 4D? This is likely a bug in AFL"); }

            size_t locate( const Point& x ) const { throw exc::Exception("Locator is only implemented from 2D to 4D? This is likely a bug in AFL"); }

            std::vector<size_t> surround( size_t x, int l ) const;

            size_t getSide() { return _n; }

        private:
            size_t _n;
            const ddf::Node* _min;
            const ddf::Node* _max;
    };

    template <size_t D>
    void triangulate( const ddf::Property* xs, std::vector<Element<D> >& elems,
            std::vector<std::set<size_t> >& cells );

    template<size_t D>
    Element<D> findElement( const ddf::Property* xs, const Point& p );

}

#include "qhullwrapper.tcc"

#endif

---