Skeleton_blocker_geometric_complex.h

 /*    This file is part of the Gudhi Library. The Gudhi library
  *    (Geometric Understanding in Higher Dimensions) is a generic C++
  *    library for computational topology.
  *
  *    Author(s):       David Salinas
  *
  *    Copyright (C) 2014  INRIA Sophia Antipolis-Mediterranee (France)
  *
  *    This program is free software: you can redistribute it and/or modify
  *    it under the terms of the GNU General Public License as published by
  *    the Free Software Foundation, either version 3 of the License, or
  *    (at your option) any later version.
  *
  *    This program is distributed in the hope that it will be useful,
  *    but WITHOUT ANY WARRANTY; without even the implied warranty of
  *    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
  *    GNU General Public License for more details.
  *
  *    You should have received a copy of the GNU General Public License
  *    along with this program.  If not, see <http://www.gnu.org/licenses/>.
  */
#ifndef SKELETON_BLOCKER_GEOMETRIC_COMPLEX_H_
#define SKELETON_BLOCKER_GEOMETRIC_COMPLEX_H_


#include "gudhi/Utils.h"
#include "gudhi/Skeleton_blocker_simplifiable_complex.h"
#include "gudhi/Skeleton_blocker/Skeleton_blocker_sub_complex.h"



namespace Gudhi{

namespace skbl {

template<typename SkeletonBlockerGeometricDS>
class Skeleton_blocker_geometric_complex : public Skeleton_blocker_simplifiable_complex<SkeletonBlockerGeometricDS>
{
public:
    typedef typename SkeletonBlockerGeometricDS::GT GT;

    typedef Skeleton_blocker_simplifiable_complex<SkeletonBlockerGeometricDS> SimplifiableSkeletonblocker;

    typedef typename SimplifiableSkeletonblocker::Vertex_handle Vertex_handle;
    typedef typename SimplifiableSkeletonblocker::Root_vertex_handle Root_vertex_handle;
    typedef typename SimplifiableSkeletonblocker::Edge_handle Edge_handle;
    typedef typename SimplifiableSkeletonblocker::Simplex_handle Simplex_handle;


    typedef typename SimplifiableSkeletonblocker::Graph_vertex Graph_vertex;

    typedef typename SkeletonBlockerGeometricDS::Point Point;


    Vertex_handle add_vertex(const Point& point){
        Vertex_handle ad = SimplifiableSkeletonblocker::add_vertex();
        (*this)[ad].point() = point;
        return ad;
    }


    const Point& point(Vertex_handle v) const{
        assert(this->contains_vertex(v));
        return (*this)[v].point();
    }

    Point& point(Vertex_handle v) {
        assert(this->contains_vertex(v));
        return (*this)[v].point();
    }

    const Point& point(Root_vertex_handle global_v) const{
        Vertex_handle local_v ( (*this)[global_v]) ;
        assert(this->contains_vertex(local_v));
        return (*this)[local_v].point();
    }

    Point& point(Root_vertex_handle global_v) {
        Vertex_handle local_v ( (*this)[global_v]) ;
        assert(this->contains_vertex(local_v));
        return (*this)[local_v].point();
    }


    typedef Skeleton_blocker_link_complex<Skeleton_blocker_geometric_complex> Geometric_link;

    Geometric_link link(const Simplex_handle& simplex) const{
        Geometric_link link(*this,simplex);
        //we now add the point info
        add_points_to_link(link);
        return link;
    }

    Geometric_link link(Edge_handle edge) const{
        Geometric_link link(*this,edge);
        //we now add the point info
        add_points_to_link(link);
        return link;
    }
private:
    void add_points_to_link(Geometric_link& link) const{
        for(Vertex_handle v : link.vertex_range()){
            Root_vertex_handle v_root(link.get_id(v));
            link.point(v) = (*this).point(v_root);
        }
    }
};

}

}  // namespace GUDHI

#endif /* SKELETON_BLOCKER_GEOMETRIC_COMPLEX_H_ */