ParallelIIEstimator.h

 
#include <DGtal/kernel/domains/DomainSplitter.h>
#include <DGtal/topology/helpers/Surfaces.h>
 
#ifndef DGTAL_WITH_OPENMP
#error You need to have activated OpenMP (DGTAL_WITH_OPENMP) to include this file.
#endif
 
// We require openmp for this estimator, even though it could run without
#include <omp.h>
 
namespace DGtal
{
  template<class TEstimator, typename TSplitter>
  class ParallelIIEstimator
  {
  public:
    using Estimator = TEstimator;
    using Splitter = TSplitter;
    using Domain = typename TEstimator::Domain;
    using Scalar = typename TEstimator::Scalar;
 
    using KSpace = typename TEstimator::KSpace;
    using PointPredicate = typename TEstimator::PointPredicate;
    using Surfel = typename KSpace::Surfel;
    using SurfelSet = typename KSpace::SurfelSet;
    using EstimatorQuantity = typename TEstimator::Quantity;
    using Quantity = EstimatorQuantity;
 
    // Building the surface
    using Boundary = LightImplicitDigitalSurface<KSpace, PointPredicate>;
    using Surface = DigitalSurface<Boundary>;
    using Visitor = DepthFirstVisitor<Surface>;
    using VisitorRange = GraphVisitorRange<Visitor>;
 

    template<typename... Args>
    ParallelIIEstimator(Splitter splitter, int32_t nbThread, Args&&... args);

    void clear();

    Scalar h() const;

    void attach(ConstAlias<KSpace> K,
                ConstAlias<PointPredicate> aPointPredicate);

    void setParams(double dRadius);
 

     bool isValid() const;

    template<typename ItA, typename ItB>
    void init(double h_, ItA ite, ItB itb);

    template<typename It>
    Quantity eval(It it);

    template<typename It, typename Oit>
    Oit eval(It itb, It ite, Oit result);

    void selfDisplay ( std::ostream & out ) const;
 
  private:
    std::vector<Estimator> myEstimators;
    Splitter mySplitter;
 
    CountedConstPtrOrConstPtr<PointPredicate> myPointPredicate;
    CountedConstPtrOrConstPtr<KSpace> myKSpace;
 
    double myH;
    double myRadius;
  };
}
 
#include "ParallelIIEstimator.ih"