NURBSShape_fwd.hpp

Go to the documentation of this file.

#ifndef MEDUSA_BITS_DOMAINS_NURBSSHAPE_FWD_HPP_
#define MEDUSA_BITS_DOMAINS_NURBSSHAPE_FWD_HPP_
 
#include <medusa/Config.hpp>
#include <medusa/bits/domains/DomainShape_fwd.hpp>
#include <medusa/bits/types/Vec.hpp>
#include <medusa/bits/types/Range.hpp>
#include <medusa/bits/domains/NURBSPatch_fwd.hpp>
#include <medusa/bits/utils/randutils.hpp>
 
namespace mm {
 
namespace nurbs_shape_internal {
template <typename vec_t, typename param_vec_t>
struct NURBSShapeHelper {};
}  // namespace nurbs_shape_internal
 
template <typename vec_t, typename param_vec_t>
class NURBSShape : public DomainShape<vec_t> {
    friend struct nurbs_shape_internal::NURBSShapeHelper<vec_t, param_vec_t>;
 
  public:
    typedef typename vec_t::scalar_t scalar_t;  
    enum {  dim = vec_t::dim, proj_dim = dim + 1, param_dim = param_vec_t::dim};
    typedef Vec<scalar_t, proj_dim> proj_vec_t;  
 
  private:
    Range<NURBSPatch<vec_t, param_vec_t>> patches;  
    int max_points = 5000000;  
    int seed_;  
    int n_samples = 15;  
    scalar_t zeta = 1 - 1e-10;  
    scalar_t epsilon = 0;  
 
  public:
    NURBSShape(const Range<NURBSPatch<vec_t, param_vec_t>>& patches_in);
 
    NURBSShape(Range<NURBSPatch<vec_t, param_vec_t>>&& patches_in);
 
    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 {
        assert_msg(false, "This function is not available for this shape. A"
                          "bounding box of its DomainDiscretization should be taken.");
        return {};
    }
 
    DomainDiscretization<vec_t> discretizeBoundaryWithDensity(
            const std::function<scalar_t(vec_t)>&, int) const override;
 
    std::ostream& print(std::ostream& ostream) const override {
        return ostream << "NURBS shape made of " << patches.size() << " patches.";
    }
 
    NURBSShape* clone() const override {
        return new NURBSShape<vec_t, param_vec_t>(*this);
    }
 
    NURBSShape& maxPoints(int max_points) { this->max_points = max_points; return *this; }
 
    NURBSShape& seed(int seed) { seed_ = seed; return *this; }
 
    NURBSShape& proximityTolerance(scalar_t zeta);
 
    NURBSShape& numSamples(int n_samples) { this->n_samples = n_samples; return *this; }
 
    NURBSShape& boundaryProximity(scalar_t epsilon);
};
 
 
}  // namespace mm
 
#endif  // MEDUSA_BITS_DOMAINS_NURBSSHAPE_FWD_HPP_