Milena (Olena)
User documentation 2.0a Id
|
00001 // Copyright (C) 2009 EPITA Research and Development Laboratory (LRDE) 00002 // 00003 // This file is part of Olena. 00004 // 00005 // Olena is free software: you can redistribute it and/or modify it under 00006 // the terms of the GNU General Public License as published by the Free 00007 // Software Foundation, version 2 of the License. 00008 // 00009 // Olena is distributed in the hope that it will be useful, 00010 // but WITHOUT ANY WARRANTY; without even the implied warranty of 00011 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU 00012 // General Public License for more details. 00013 // 00014 // You should have received a copy of the GNU General Public License 00015 // along with Olena. If not, see <http://www.gnu.org/licenses/>. 00016 // 00017 // As a special exception, you may use this file as part of a free 00018 // software project without restriction. Specifically, if other files 00019 // instantiate templates or use macros or inline functions from this 00020 // file, or you compile this file and link it with other files to produce 00021 // an executable, this file does not by itself cause the resulting 00022 // executable to be covered by the GNU General Public License. This 00023 // exception does not however invalidate any other reasons why the 00024 // executable file might be covered by the GNU General Public License. 00025 00026 #ifndef MLN_MAKE_CELL_HH 00027 # define MLN_MAKE_CELL_HH 00028 00031 00032 # include <mln/core/site_set/p_set.hh> 00033 # include <mln/core/site_set/complex_psite.hh> 00034 # include <mln/core/image/complex_neighborhoods.hh> 00035 # include <mln/core/image/complex_neighborhood_piter.hh> 00036 00037 # include <mln/topo/is_facet.hh> 00038 00039 namespace mln 00040 { 00041 00042 namespace make 00043 { 00044 00053 template <unsigned D, typename G> 00054 p_set< complex_psite<D, G> > 00055 cell(const complex_psite<D, G>& f); 00056 00057 00058 # ifndef MLN_INCLUDE_ONLY 00059 00060 template <unsigned D, typename G> 00061 inline 00062 p_set< complex_psite<D, G> > 00063 cell(const complex_psite<D, G>& f) 00064 { 00065 mln_precondition(topo::is_facet(f)); 00066 00067 typedef complex_psite<D, G> psite; 00068 typedef p_set<psite> faces_t; 00069 00070 // Compute the cell F^HAT. 00071 faces_t f_hat; 00072 /* FIXME: We need a cell-iterator here 00073 (see https://trac.lrde.org/olena/ticket/162). */ 00074 typedef complex_m_face_neighborhood<D, G> m_faces_nbh_t; 00075 m_faces_nbh_t m_faces_nbh; 00076 mln_niter(m_faces_nbh_t) g(m_faces_nbh, f); 00077 for (unsigned m = 0; m < f.n(); ++m) 00078 { 00079 g.iter().set_m(m); 00080 for_all(g) 00081 { 00082 f_hat.insert(g); 00083 } 00084 } 00085 f_hat.insert(f); 00086 return f_hat; 00087 } 00088 00089 # endif // MLN_INCLUDE_ONLY 00090 00091 } // end of namespace mln::make 00092 00093 } // end of namespace mln 00094 00095 #endif // ! MLN_MAKE_CELL_HH