Milena (Olena)  User documentation 2.0a Id
 All Classes Namespaces Functions Variables Typedefs Enumerator Groups Pages
detachment.hh
1 // Copyright (C) 2009 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 MLN_MAKE_DETACHMENT_HH
27 # define MLN_MAKE_DETACHMENT_HH
28 
32 
33 # include <mln/core/image/complex_image.hh>
34 # include <mln/make/cell.hh>
35 # include <mln/topo/is_facet.hh>
36 
37 namespace mln
38 {
39 
40  namespace make
41  {
42 
57  template <unsigned D, typename G, typename V>
58  p_set< complex_psite<D, G> >
59  detachment(const complex_psite<D, G>& f,
60  const complex_image<D, G, V>& ima);
61 
62 
63 # ifndef MLN_INCLUDE_ONLY
64 
65  template <unsigned D, typename G, typename V>
66  inline
67  p_set< complex_psite<D, G> >
69  const complex_image<D, G, V>& ima)
70  {
71  mln_precondition(topo::is_facet(f));
72  mlc_equal(V, bool)::check();
73 
74  typedef complex_psite<D, G> psite;
75  typedef p_set<psite> faces_t;
76 
77  faces_t f_hat = make::cell(f);
78  // Initialize DETACH_F to F_HAT.
79  faces_t detach_f = f_hat;
80 
81  typedef complex_lower_higher_neighborhood<D, G> adj_nbh_t;
82  adj_nbh_t adj_nbh;
83  mln_piter(faces_t) g(f_hat);
84  mln_niter(adj_nbh_t) n(adj_nbh, g);
85  for_all(g)
86  for_all(n)
87  if (ima(n) && !f_hat.has(n))
88  {
89  detach_f.remove(g);
90  break;
91  }
92  return detach_f;
93  }
94 
95 # endif // MLN_INCLUDE_ONLY
96 
97  } // end of namespace mln::make
98 
99 } // end of namespace mln
100 
101 #endif // ! MLN_MAKE_DETACHMENT_HH