vertex_one_ring_barycenter.hpp

Go to the documentation of this file.

// Copyright (c) 2012-2019 University of Lyon and CNRS (France).
// All rights reserved.
//
// This file is part of MEPP2; you can redistribute it and/or modify
// it under the terms of the GNU Lesser General Public License as
// published by the Free Software Foundation; either version 3 of
// the License, or (at your option) any later version.
//
// This file is provided AS IS with NO WARRANTY OF ANY KIND, INCLUDING THE
// WARRANTY OF DESIGN, MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE.
#pragma once
 
#include <boost/graph/graph_traits.hpp>
#include <boost/graph/properties.hpp>
#include <vector>
 
#include "FEVV/Wrappings/Geometry_traits.h"
#include "FEVV/Operators/Generic/Manifold/k_ring.hpp"
 
 
namespace FEVV {
namespace Operators {
 
template< typename FaceGraph,
          typename PointMap,
          typename GeometryTraits = FEVV::Geometry_traits< FaceGraph > >
typename boost::property_traits< PointMap >::value_type
vertex_one_ring_barycenter(
    typename boost::graph_traits< FaceGraph >::vertex_descriptor v,
    const FaceGraph &g,
    const PointMap &pm,
    const typename GeometryTraits::Scalar smoothing_factor,
    const GeometryTraits &gt)
{
  typedef typename boost::property_traits< PointMap >::value_type
      Point; // this Point type cannot be const because the center Point must be
             // modifiable
  //typedef typename boost::property_traits< PointMap >::reference Reference;
 
  std::vector< typename boost::graph_traits< FaceGraph >::vertex_descriptor >
      qv;
  Operators::extract_1_ring_not_including_v< FaceGraph, GeometryTraits >(v, g, qv);
 
  Point center(0, 0, 0);
  unsigned int cpt = 0;
  typename std::vector< typename boost::graph_traits<
      const FaceGraph >::vertex_descriptor >::const_iterator it(qv.begin()),
      ite(qv.end());
  for(; it != ite; ++it, ++cpt)
  {
    Point tmp = get(pm, *it); // works for all mesh structures
    center = Point(gt.get_x(tmp) + gt.get_x(center),
                   gt.get_y(tmp) + gt.get_y(center),
                   gt.get_z(tmp) + gt.get_z(center));
  }
  if(cpt > 0)
  {
    return gt.add_pv(get(pm, v),
                    gt.scalar_mult(gt.sub_p(Point(gt.get_x(center) / cpt,
                                                  gt.get_y(center) / cpt,
                                                  gt.get_z(center) / cpt),
                                            get(pm, v)),
                                   smoothing_factor));
  }
  else
    return get(pm, v);
}
 
} // namespace Operators
} // namespace FEVV