$extrastylesheet
Olena  User documentation 2.1
An Image Processing Platform
 All Classes Namespaces Functions Variables Typedefs Enumerations Enumerator Friends Groups Pages
box_intersection.hh
1 // Copyright (C) 2011 EPITA Research and Development Laboratory (LRDE)
2 //
3 // This file is part of Olena.
4 //
5 // Olena is free software: you can redistribute it and/or modify it under
6 // the terms of the GNU General Public License as published by the Free
7 // Software Foundation, version 2 of the License.
8 //
9 // Olena is distributed in the hope that it will be useful,
10 // but WITHOUT ANY WARRANTY; without even the implied warranty of
11 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
12 // General Public License for more details.
13 //
14 // You should have received a copy of the GNU General Public License
15 // along with Olena. If not, see <http://www.gnu.org/licenses/>.
16 //
17 // As a special exception, you may use this file as part of a free
18 // software project without restriction. Specifically, if other files
19 // instantiate templates or use macros or inline functions from this
20 // file, or you compile this file and link it with other files to produce
21 // an executable, this file does not by itself cause the resulting
22 // executable to be covered by the GNU General Public License. This
23 // exception does not however invalidate any other reasons why the
24 // executable file might be covered by the GNU General Public License.
25 
26 #ifndef SCRIBO_UTIL_BOX_INTERSECTION_HH
27 # define SCRIBO_UTIL_BOX_INTERSECTION_HH
28 
32 
33 
34 #include <mln/core/site_set/box.hh>
35 #include <mln/math/min.hh>
36 #include <mln/math/max.hh>
37 
38 namespace scribo
39 {
40 
41  namespace util
42  {
43  using namespace mln;
44 
50  template <typename P>
51  box<P>
52  box_intersection(const box<P>& lhs, const box<P>& rhs);
53 
54 
55 # ifndef MLN_INCLUDE_ONLY
56 
57  template <typename P>
58  box<P>
59  box_intersection(const box<P>& lhs, const box<P>& rhs)
60  {
61  mln_trace("scribo::util::box_intersection");
62 
63  P pmin, pmax;
64  for (unsigned dim = 0; dim < P::dim; ++dim)
65  {
66  pmin[dim] = math::max(lhs.pmin()[dim], rhs.pmin()[dim]);
67  pmax[dim] = math::min(lhs.pmax()[dim], rhs.pmax()[dim]);
68 
69  if (pmin[dim] > pmax[dim]) // No intersection.
70  return box<P>();
71  }
72 
73  box<P> output(pmin, pmax);
74 
75  return output;
76  }
77 
78 # endif // ! MLN_INCLUDE_ONLY
79 
80  } // end of namespace scribo::util
81 
82 } // end of namespace scribo
83 
84 #endif // ! SCRIBO_UTIL_BOX_INTERSECTION_HH