Milena (Olena)  User documentation 2.0a Id
 All Classes Namespaces Functions Variables Typedefs Enumerator Groups Pages
registration/internal/rms.hh
1 // Copyright (C) 2008, 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_REGISTRATION_INTERNAL_RMS_HH
27 # define MLN_REGISTRATION_INTERNAL_RMS_HH
28 
29 # include <mln/norm/l2.hh>
30 # include <mln/core/site_set/p_array.hh>
31 
32 namespace mln
33 {
34 
35  namespace registration
36  {
37 
38  namespace internal
39  {
40 
41  template <typename P, typename M, typename T>
42  float rms(const p_array<P>& a1,
43  M& map,
44  const unsigned length,
45  const T& q1,
46  const T& q2);
47 
48 # ifndef MLN_INCLUDE_ONLY
49 
50 
51  template <typename P, typename M, typename T>
52  float rms(const p_array<P>& a,
53  M& map,
54  const T& q1,
55  const T& q2)
56  {
57  float f = 0.f;
58  for (unsigned i = 0; i < a.nsites(); ++i)
59  {
60  algebra::vec<P::dim,float> a2f = convert::to< algebra::vec<P::dim,float> >(a[i]);
61  algebra::vec<P::dim,float> a1f =
62  convert::to< algebra::vec<P::dim,float> > \
63  (map(q1(convert::to< algebra::vec<P::dim,float> >(a[i]))));
64 
65  f += norm::l2(a1f - q2(a2f));
66  }
67  return f / a.nsites();
68  }
69 
70 
71 # endif // ! MLN_INCLUDE_ONLY
72 
73  } // end of nqmespace mln::registration::internal
74 
75  } // end of namespace mln::registration
76 
77 } // end of namespace mln
78 
79 
80 #endif // ! MLN_REGISTRATION_INTERNAL_RMS_HH