Milena (Olena)  User documentation 2.0a Id
 All Classes Namespaces Functions Variables Typedefs Enumerator Groups Pages
diagonal2d.hh
1 // Copyright (C) 2007, 2008, 2009, 2010 EPITA Research and Development
2 // Laboratory (LRDE)
3 //
4 // This file is part of Olena.
5 //
6 // Olena is free software: you can redistribute it and/or modify it under
7 // the terms of the GNU General Public License as published by the Free
8 // Software Foundation, version 2 of the License.
9 //
10 // Olena is distributed in the hope that it will be useful,
11 // but WITHOUT ANY WARRANTY; without even the implied warranty of
12 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
13 // General Public License for more details.
14 //
15 // You should have received a copy of the GNU General Public License
16 // along with Olena. If not, see <http://www.gnu.org/licenses/>.
17 //
18 // As a special exception, you may use this file as part of a free
19 // software project without restriction. Specifically, if other files
20 // instantiate templates or use macros or inline functions from this
21 // file, or you compile this file and link it with other files to produce
22 // an executable, this file does not by itself cause the resulting
23 // executable to be covered by the GNU General Public License. This
24 // exception does not however invalidate any other reasons why the
25 // executable file might be covered by the GNU General Public License.
26 
27 #ifndef MLN_CANVAS_BROWSING_DIAGONAL2D_HH
28 # define MLN_CANVAS_BROWSING_DIAGONAL2D_HH
29 
33 
34 # include <mln/core/concept/browsing.hh>
35 # include <mln/core/concept/image.hh>
36 
37 namespace mln
38 {
39 
40  namespace canvas
41  {
42 
43  namespace browsing
44  {
45 
47 
82  struct diagonal2d_t : public Browsing< diagonal2d_t >
83  {
84  // This default constructor is needed for compilation with gcc
85  // 4.6.0, gcc 4.6.1 and Clang.
86  diagonal2d_t();
87 
88  template <typename F>
89  void operator()(F& f) const;
90  };
91 
92  extern const diagonal2d_t diagonal2d;
93 
94 # ifndef MLN_INCLUDE_ONLY
95 
96 # ifndef MLN_WO_GLOBAL_VARS
97 
98  const diagonal2d_t diagonal2d;
99 
100 # endif // ! MLN_WO_GLOBAL_VARS
101 
102  inline
103  diagonal2d_t::diagonal2d_t()
104  {
105  }
106 
107 
108  template <typename F>
109  inline
110  void
111  diagonal2d_t::operator()(F& f) const
112  {
113  trace::entering("canvas::browsing::diagonal2d");
114 
115  typedef typename F::I I;
116  typedef mln_deduce(I, psite, delta) dpsite;
117  typedef mln_psite(I) psite;
118 
119  // Directions
120  dpsite dp_first(1, 0);
121  dpsite dp_second(0, 1);
122  dpsite dp_diag(-1, 1);
123 
124  mln_deduce(I, psite, delta) diag;
125 
126  psite
127  pmin = f.input.domain().pmin(),
128  pmax = f.input.domain().pmax();
129 
130  f.p = pmin;
131 
132  f.init();
133 
134  while (f.input.domain().has(f.p))
135  {
136  psite start_diag = f.p;
137  // Browse one diag.
138  f.init_diag();
139  while (f.input.domain().has(f.p))
140  {
141  f.next();
142  f.p = f.p + dp_diag;
143  }
144 
145  // Goto next diag start.
146  f.p = start_diag + dp_first;
147  }
148  f.p = f.p - dp_first;
149 
150  f.p = f.p + dp_second;
151  while (f.input.domain().has(f.p))
152  {
153  psite start_diag = f.p;
154 
155  // Browse one diag.
156  f.init_diag();
157  while (f.input.domain().has(f.p))
158  {
159  f.next();
160  f.p = f.p + dp_diag;
161  }
162  f.p = f.p - dp_diag;
163 
164  // Goto next diag start.
165  f.p = start_diag + dp_second;
166  }
167 
168  f.final();
169  trace::exiting("canvas::browsing::diagonal2d");
170  }
171 
172 # endif // ! MLN_INCLUDE_ONLY
173 
174  } // end of namespace mln::canvas::browsing
175 
176  } // end of namespace mln::canvas
177 
178 } // end of namespace mln
179 
180 #endif // ! MLN_CANVAS_BROWSING_DIAGONAL2D_HH