PolygonShape_fwd.hpp

Go to the documentation of this file.

#ifndef MEDUSA_BITS_DOMAINS_POLYGONSHAPE_FWD_HPP_
#define MEDUSA_BITS_DOMAINS_POLYGONSHAPE_FWD_HPP_
 
#include <medusa/Config.hpp>
#include "DomainShape_fwd.hpp"
 
namespace mm {
 
template <typename vec_t>
class PolygonShape : public DomainShape<vec_t> {
    static_assert(vec_t::dim == 2, "Only available in 2 dimensions.");
    using base_t = DomainShape<vec_t>;  
 
    std::vector<vec_t> points_;  
    std::vector<vec_t> points_with_margin_;
    using DomainShape<vec_t>::margin_;
 
  public:
    using typename base_t::scalar_t;
    using typename base_t::vector_t;
    using base_t::dim;
    using base_t::discretizeBoundaryWithDensity;
    using base_t::discretizeBoundaryWithStep;
 
    explicit PolygonShape(const std::vector<vec_t>& points);
 
    const std::vector<vec_t>& points() const { return points_; }
 
    void setMargin(scalar_t margin) override;
 
    bool contains(const vec_t& point) const override;
 
    std::pair<vec_t, vec_t> bbox() const override;
 
    DomainDiscretization<vec_t> discretizeBoundaryWithStep(scalar_t step, int type) const override;
 
     DomainDiscretization<vec_t> discretizeBoundaryWithDensity(
            const std::function<scalar_t(vec_t)>& dr, int type) const override;
 
    PolygonShape<vec_t>* clone() const override { return new PolygonShape<vec_t>(*this); }
 
    std::ostream& print(std::ostream& os) const override {
        return os << "Polygon shape with " << points_.size() << " points.";
    }
 
  private:
    static inline scalar_t isLeft(const vec_t& P0, const vec_t& P1, const vec_t& P2) {
        return (P1[0] - P0[0]) * (P2[1] - P0[1]) - (P2[0] - P0[0]) * (P1[1] - P0[1]);
    }
 
    void extendByMargin();
};
 
}  // namespace mm
 
#endif  // MEDUSA_BITS_DOMAINS_POLYGONSHAPE_FWD_HPP_