00001 // Copyright (C) 2007, 2008, 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 #include <cmath> 00027 #include <cassert> 00028 00029 #include <tests/norm/common.hh> 00030 00031 #include <mln/algebra/vec.hh> 00032 #include <mln/norm/l2.hh> 00033 00034 00035 namespace test 00036 { 00037 template <typename V, typename S> 00038 void 00039 check_l2_norm_and_distance(const V& vec1, const V& vec2, const S& ref_val) 00040 { 00041 // Pointer on mln::norm::l2. 00042 typedef mln_sum_(int) (*l2_t)(const V&); 00043 l2_t l2 = mln::norm::l2; 00044 00045 test::check_norm(l2, vec1, vec2); 00046 00047 // Pointer on mln::norm::l2_distance. 00048 typedef mln_sum_(int) (*l2_distance_t)(const V&, const V&); 00049 l2_distance_t l2_distance = mln::norm::l2_distance; 00050 00051 test::check_distance(l2_distance, vec1, vec2, ref_val); 00052 } 00053 } 00054 00055 00056 int main() 00057 { 00058 typedef mln::algebra::vec<3, int> vec_t; 00059 00060 // Reference value. 00061 float d = std::sqrt(float((4 - 2) * (4 - 2) + 00062 (1 + 2) * (1 + 2) + 00063 (0 - 3) * (0 - 3))); 00064 00065 // Tests using mln::algebra::vec. 00066 vec_t t, u; 00067 t.set(2, -2, 3); 00068 u.set(4, 1, 0); 00069 test::check_l2_norm_and_distance(t, u, d); 00070 00071 // Tests using plain arrays. 00072 int v[] = {2, -2, 3}; 00073 int w[] = {4, 1, 0}; 00074 test::check_l2_norm_and_distance(v, w, d); 00075 }