STLShape_fwd.hpp

Go to the documentation of this file.

#ifndef MEDUSA_BITS_DOMAINS_STLSHAPE_FWD_HPP_
#define MEDUSA_BITS_DOMAINS_STLSHAPE_FWD_HPP_
 
#include "DomainShape_fwd.hpp"
#include <medusa/bits/io/STL_fwd.hpp>
#include <utility>
 
namespace mm {
 
template <typename vec_t>
class STLShape : public DomainShape<vec_t> {
  public:
    using typename DomainShape<vec_t>::scalar_t;
    using DomainShape<vec_t>::dim;
 
  private:
    static_assert(dim == 3, "STL files are for 3D models only.");
    std::vector<vec_t> vertices_;  
    std::vector<std::array<int, 3>> faces_;  
    std::vector<vec_t> normals_;  
    std::pair<vec_t, vec_t> bbox_;  
 
  public:
    STLShape(const std::vector<STL::Triangle>& stl);
 
    bool contains(const vec_t& /* point */) const override { return true; }
 
    bool hasContains() const override { return false; }
 
    std::pair<vec_t, vec_t> bbox() const override { return bbox_; }
 
    using DomainShape<vec_t>::discretizeBoundaryWithDensity;
    DomainDiscretization<vec_t> discretizeBoundaryWithDensity(
            const std::function<scalar_t(vec_t)>&, int) const override;
 
    std::ostream& print(std::ostream& ostream) const override {
        return ostream << "STL shape with " << vertices_.size() << " vertices and "
                       << faces_.size() << " faces.";
    }
 
    STLShape* clone() const override {
        return new STLShape<vec_t>(*this);
    }
 
  private:
    static vec_t v(const STL::Point& p) {
        return vec_t(p.x, p.y, p.z);
    }
};
 
}  // namespace mm
 
#endif  // MEDUSA_BITS_DOMAINS_STLSHAPE_FWD_HPP_