polytope.hpp

Go to the documentation of this file.

 
#ifndef MONTECARLO_1_POLYTOPE_H
#define MONTECARLO_1_POLYTOPE_H
 
#include "integration_domain.hpp"
#include "../geometry.hpp"
#include <vector>
#include <stdexcept>
 
namespace mc::domains {
 
template <size_t dim>
class PolyTope : public IntegrationDomain<dim> {
public:
 
    PolyTope(const std::vector<geom::Point<dim>>&   vertices,
             const std::vector<std::array<double, dim>>&  norms,
             const std::vector<double>&      offs);
 
    geom::Bounds<dim> getBounds() const override;
 
    double getBoxVolume() const override;
 
    bool isInside(const geom::Point<dim> &p) const override;
 
private:
    std::vector<mc::geom::Point<dim>> vec;
    std::vector<std::array<double, dim>> normals;
    std::vector<double> offsets;
};
 
} // namespace mc::domains
 
#include "polytope.tpp"
 
#endif //MONTECARLO_1_POLYTOPE_H