politopy.h

Go to the documentation of this file.

// politopix allows to make computations on polytopes such as finding vertices, intersecting, Minkowski sums, ...
//     Copyright (C) 2021 : Delos Vincent
//
//     This program 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.
//
//     This program 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 this program.  If not, see <http://www.gnu.org/licenses/>.
//
// I2M (UMR CNRS 5295 / University of Bordeaux)
 
#include <boost/python.hpp>
#include <boost/python/suite/indexing/vector_indexing_suite.hpp>
 
using namespace boost::python;
 
 
// For boost python
static polito_EXPORT int computeMinkowskiSumOfPolytopesPy(const boost::python::list& listOfPol, boost::shared_ptr<Polytope_Rn>& C) {
  std::vector< boost::shared_ptr<Polytope_Rn> > arrayOfPol;
  for (unsigned int i=0; i<boost::python::len(listOfPol); ++i)
    arrayOfPol.push_back(boost::python::extract< boost::shared_ptr<Polytope_Rn> >(listOfPol[i]));
  politopixAPI::computeMinkowskiSumOfPolytopes(arrayOfPol, C);
}
 
int (*computeMinkowskiSumOfPolytopes1)(const boost::shared_ptr<Polytope_Rn>&, const boost::shared_ptr<Polytope_Rn>&, boost::shared_ptr<Polytope_Rn>&) = &politopixAPI::computeMinkowskiSumOfPolytopes;
int (*computeMinkowskiSumOfPolytopes2)(const std::vector< boost::shared_ptr<Polytope_Rn> >&, boost::shared_ptr<Polytope_Rn>&) = &politopixAPI::computeMinkowskiSumOfPolytopes;
int (*computeMinkowskiSumOfPolytopes3)(const boost::python::list&, boost::shared_ptr<Polytope_Rn>&) = &computeMinkowskiSumOfPolytopesPy;
int (*computeIntersection)(const boost::shared_ptr<Polytope_Rn>&, const boost::shared_ptr<Polytope_Rn>&, boost::shared_ptr<Polytope_Rn>&) = &politopixAPI::computeIntersection;
int (*computeDoubleDescription2)(boost::shared_ptr<Polytope_Rn>&, double bb_size) = &politopixAPI::computeDoubleDescription;
int (*computeDoubleDescription1)(boost::shared_ptr<Polytope_Rn>&) = &politopixAPI::computeDoubleDescription;
 
 
BOOST_PYTHON_MODULE(libpolito) {
  Py_Initialize();
 
  class_<Generator_Rn, boost::shared_ptr<Generator_Rn> >("generator", init<>())
    .def("print", &Generator_Rn::print)
    .def("setCoord", &Generator_Rn::setCoordinate)
    .def("getCoord", &Generator_Rn::getCoordinate)
    ;
  
  class_<HalfSpace_Rn, boost::shared_ptr<HalfSpace_Rn> >("halfspace", init<>())
    .def("print", &HalfSpace_Rn::print)
    .def("setCoef", &HalfSpace_Rn::setCoefficient)
    .def("getCoef", &HalfSpace_Rn::getCoefficient)
    .def("setConstant", &HalfSpace_Rn::setConstant)
    .def("getConstant", &HalfSpace_Rn::getConstant)
    ;
  
  class_<PolyhedralCone_Rn, boost::shared_ptr<PolyhedralCone_Rn> >("polyhedralcone", init<>())
    .def("getGen", &PolyhedralCone_Rn::getGenerator, return_value_policy<copy_const_reference>())
    .def("nbGen", &PolyhedralCone_Rn::numberOfGenerators)    
    .def("getHS", &PolyhedralCone_Rn::getHalfSpace, return_value_policy<copy_const_reference>())
    .def("nbHS",  &PolyhedralCone_Rn::numberOfHalfSpaces)
    ;
 
  class_<Polytope_Rn, boost::shared_ptr<Polytope_Rn>, bases<PolyhedralCone_Rn> >("polytope", init<>())
    .def("create", &Polytope_Rn::create).staticmethod("create")
    ;
 
  class_<politopixAPI>("politopix")
    .def("loadPolytope", &politopixAPI::loadPolytope).staticmethod("loadPolytope")
    .def("savePolytope", &politopixAPI::savePolytope).staticmethod("savePolytope")
    .def("volume", &politopixAPI::computeVolume).staticmethod("volume")
    .def("check", &politopixAPI::checkTopologyAndGeometry).staticmethod("check")
    .def("cube", &politopixAPI::makeCube).staticmethod("cube")    
    .def("sum", computeMinkowskiSumOfPolytopes1)
    .def("sum", computeMinkowskiSumOfPolytopes2)
    .def("sum", computeMinkowskiSumOfPolytopes3).staticmethod("sum")
    .def("doubledesc", computeDoubleDescription1)
    .def("doubledesc", computeDoubleDescription2).staticmethod("doubledesc")
    .def("intersect", computeIntersection).staticmethod("intersect")
    ;
 
  class_<Rn>("Rn")
    .def("getDim", &Rn::getDimension).staticmethod("getDim")
    .def("setDim", &Rn::setDimension).staticmethod("setDim")
    .def("getTol", &Rn::getTolerance).staticmethod("getTol")
    .def("setTol", &Rn::setTolerance).staticmethod("setTol")
    ;
 
}