Milena (Olena)  User documentation 2.0a Id
 All Classes Namespaces Functions Variables Typedefs Enumerator Groups Pages
cell.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_CELL_HH
27 # define MLN_MAKE_CELL_HH
28 
31 
32 # include <mln/core/site_set/p_set.hh>
33 # include <mln/core/site_set/complex_psite.hh>
34 # include <mln/core/image/complex_neighborhoods.hh>
35 # include <mln/core/image/complex_neighborhood_piter.hh>
36 
37 # include <mln/topo/is_facet.hh>
38 
39 namespace mln
40 {
41 
42  namespace make
43  {
44 
53  template <unsigned D, typename G>
54  p_set< complex_psite<D, G> >
55  cell(const complex_psite<D, G>& f);
56 
57 
58 # ifndef MLN_INCLUDE_ONLY
59 
60  template <unsigned D, typename G>
61  inline
62  p_set< complex_psite<D, G> >
64  {
65  mln_precondition(topo::is_facet(f));
66 
67  typedef complex_psite<D, G> psite;
68  typedef p_set<psite> faces_t;
69 
70  // Compute the cell F^HAT.
71  faces_t f_hat;
72  /* FIXME: We need a cell-iterator here
73  (see https://trac.lrde.org/olena/ticket/162). */
74  typedef complex_m_face_neighborhood<D, G> m_faces_nbh_t;
75  m_faces_nbh_t m_faces_nbh;
76  mln_niter(m_faces_nbh_t) g(m_faces_nbh, f);
77  for (unsigned m = 0; m < f.n(); ++m)
78  {
79  g.iter().set_m(m);
80  for_all(g)
81  {
82  f_hat.insert(g);
83  }
84  }
85  f_hat.insert(f);
86  return f_hat;
87  }
88 
89 # endif // MLN_INCLUDE_ONLY
90 
91  } // end of namespace mln::make
92 
93 } // end of namespace mln
94 
95 #endif // ! MLN_MAKE_CELL_HH