src/include/Hpolyhedron.h

/*
 * Hpolyhedron.h
 *
 *  Created on: Nov 20, 2012
 *      Author: romain
 */

#ifndef HPOLYHEDRON_H_
#define HPOLYHEDRON_H_

#include "Set.h"
#include "utils.h"
#include "Box.h"
#include "Simplex.h"

namespace nltool {

class Hpolyhedron: public Set {
public:
        Hpolyhedron();

        Hpolyhedron(int dim, Template& mat, Vector& vec);

        Hpolyhedron(int dim, Template& mat, Vector& vec, Box& box);

        Hpolyhedron(Simplex& simplex);

        Hpolyhedron(Box& box);

        virtual ~Hpolyhedron();

        bool includes(Hpolyhedron& poly);

        Number optimize(const Vector & direction, Vector & optPoint,
                        bool max = true) {
                return 0;
        }

        bool includes(Vector& point);

        std::vector<Hpolyhedron> split();

        int clean();

        Number support(Vector& a);

        /*
         * ##########################################
         * GETTTERS and SETTERS
         * ##########################################
         */

        Vector getCoefficients() const;
        Template getTemplates() const;

        void setCoefficients(Vector coefficients);
        void setTemplates(Template templates);

        Template getAllTemplates();
        Vector getAllCoefficients();

        Vertices getVertices();

        void inflate(Number inflat);

        const Box& getBbox() const {
                return bbox;
        }

        void setBbox(const Box& bbox) {
                this->bbox = bbox;
        }

private:
        Template templates;
        Vector coefficients;
        Box bbox; //Optional information
        Vertices vertices; //Optional information
};

inline std::ostream& operator<<(std::ostream& out, Hpolyhedron poly) {
        for (unsigned int i = 0; i < poly.getTemplates().size(); i++) {
                for (unsigned int j = 0; j < poly.getDimension() - 1; j++) {
                        out << poly.getTemplates()[i][j] << "x" << j << " + ";
                }
                out << poly.getTemplates()[i][poly.getDimension() - 1] << "x"
                                << poly.getDimension() - 1 << " <= "
                                << poly.getCoefficients()[i] << std::endl;
        }
        out << std::endl;
        return out;
}

} // namespace nltool

#endif /* HPOLYHEDRON_H_ */