compute_mesh_bounding_box.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 "FEVV/Wrappings/Geometry_traits.h"
 
namespace FEVV {
namespace Filters {
template< typename PropertyGraph,
          typename PointMap,
          typename GeometryTraits = FEVV::Geometry_traits< PropertyGraph > >
void
compute_mesh_bounding_box(
    const PropertyGraph &g,
    const PointMap &pm,
    typename boost::property_traits< PointMap >::value_type &min_aabb,
    typename boost::property_traits< PointMap >::value_type &max_aabb,
    const GeometryTraits &gt)
{
  typedef typename boost::property_traits< PointMap >::value_type Point;
 
  //typedef boost::graph_traits< PropertyGraph > GraphTraits;
  auto iterator_pair = vertices(g); // vertices() returns a vertex_iterator pair
  auto vi = iterator_pair.first;
  auto vi_end = iterator_pair.second;
  if(vi != vi_end)
  {
    Point p = get(pm, *vi);
    min_aabb = Point(gt.get_x(p), gt.get_y(p), gt.get_z(p));
    max_aabb = Point(gt.get_x(p), gt.get_y(p), gt.get_z(p));
    ++vi;
    for(; vi != vi_end; ++vi)
    {
      Point p = get(pm, *vi);
      min_aabb = Point(std::min(gt.get_x(p), gt.get_x(min_aabb)),
                       std::min(gt.get_y(p), gt.get_y(min_aabb)),
                       std::min(gt.get_z(p), gt.get_z(min_aabb)));
 
      max_aabb = Point(std::max(gt.get_x(p), gt.get_x(max_aabb)),
                       std::max(gt.get_y(p), gt.get_y(max_aabb)),
                       std::max(gt.get_z(p), gt.get_z(max_aabb)));
    }
  }
}
 
} // namespace Filters
} // namespace FEVV