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_DETACHMENT_HH 00027 # define MLN_MAKE_DETACHMENT_HH 00028 00032 00033 # include <mln/core/image/complex_image.hh> 00034 # include <mln/make/cell.hh> 00035 # include <mln/topo/is_facet.hh> 00036 00037 namespace mln 00038 { 00039 00040 namespace make 00041 { 00042 00057 template <unsigned D, typename G, typename V> 00058 p_set< complex_psite<D, G> > 00059 detachment(const complex_psite<D, G>& f, 00060 const complex_image<D, G, V>& ima); 00061 00062 00063 # ifndef MLN_INCLUDE_ONLY 00064 00065 template <unsigned D, typename G, typename V> 00066 inline 00067 p_set< complex_psite<D, G> > 00068 detachment(const complex_psite<D, G>& f, 00069 const complex_image<D, G, V>& ima) 00070 { 00071 mln_precondition(topo::is_facet(f)); 00072 mlc_equal(V, bool)::check(); 00073 00074 typedef complex_psite<D, G> psite; 00075 typedef p_set<psite> faces_t; 00076 00077 faces_t f_hat = make::cell(f); 00078 // Initialize DETACH_F to F_HAT. 00079 faces_t detach_f = f_hat; 00080 00081 typedef complex_lower_higher_neighborhood<D, G> adj_nbh_t; 00082 adj_nbh_t adj_nbh; 00083 mln_piter(faces_t) g(f_hat); 00084 mln_niter(adj_nbh_t) n(adj_nbh, g); 00085 for_all(g) 00086 for_all(n) 00087 if (ima(n) && !f_hat.has(n)) 00088 { 00089 detach_f.remove(g); 00090 break; 00091 } 00092 return detach_f; 00093 } 00094 00095 # endif // MLN_INCLUDE_ONLY 00096 00097 } // end of namespace mln::make 00098 00099 } // end of namespace mln 00100 00101 #endif // ! MLN_MAKE_DETACHMENT_HH