codac_ConnectedSubset.h

Go to the documentation of this file.

 
#ifndef __CODAC_CONNECTEDSUBSET_H__
#define __CODAC_CONNECTEDSUBSET_H__
 
#include <functional>
#include "codac_IntervalMatrix.h"
#include "codac_Vector.h"
#include "codac_Set.h"
 
namespace codac
{
  class Paving;
 
  class ConnectedSubset : public Set
  {
    public:
 
 
      ConnectedSubset(const std::vector<const Paving*>& v_subset_items);
 
      ~ConnectedSubset();
 
      bool is_strictly_included_in_paving() const;
 
      bool contains(const Vector& p) const;
 
      const Paving* get_paving() const;
 
      const std::vector<const Paving*>& get_items() const; // items of type k
 
      std::vector<IntervalVector> get_boxes() const; // items of type k
 
      std::vector<IntervalVector> get_boundary(SetValue value_boundary = SetValue::UNKNOWN, SetValue value_out = SetValue::OUT) const; // items of type k-1
 
      static const std::vector<IntervalVector> get_boxed_hulls(const std::vector<ConnectedSubset> v_subsets);
 
 
      bool zero_proven(const std::function<IntervalVector(const IntervalVector&)>& f);
 
      int zeros_number(const std::function<IntervalVector(const IntervalVector&)>& f, const std::function<IntervalMatrix(const IntervalVector&)>& Jf, float precision);
 
 
    protected:
 
 
      int topological_degree(const std::function<IntervalVector(const IntervalVector&)>& f);
 
      bool all_positive_signs(const std::vector<int>& v_s) const;
 
      int orientation(const IntervalVector& b, const IntervalVector& parent_coface, int orientation) const;
 
      int compute_local_degree(const std::function<IntervalVector(const IntervalVector&)>& f, const IntervalVector& b, const IntervalVector& common_cocoface) const;
      
      std::vector<int> sign_vector(const std::function<IntervalVector(const IntervalVector&)>& f, const IntervalVector& b, const IntervalVector& common_cocoface) const;
      
      std::vector<IntervalVector> get_cofaces(const IntervalVector& b) const;
 
      std::vector<IntervalVector> get_cofaces(const IntervalVector& b, const IntervalVector& common_cocoface) const;
 
      int box_dimension(const IntervalVector& b) const;
 
      bool non_singular_jacobian(const std::function<IntervalMatrix(const IntervalVector&)>& Jf, float precision);
 
 
    protected:
 
      std::vector<const Paving*> m_v_subset_items; 
  };
}
 
#endif