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 <iostream> 00028 #include <mln/fun/x2x/translation.hh> 00029 #include <mln/fun/i2v/all_to.hh> 00030 00031 00032 00033 int main() 00034 { 00035 using namespace mln; 00036 00037 float 00038 a = 2.3, 00039 b = 0, 00040 c = 2.9; 00041 00042 algebra::vec<3,float> vec1 = make::vec(a, b, c); 00043 fun::x2x::translation<3,float> tr1(all_to(1.6)); 00044 00045 std::cout << vec1 << std::endl; 00046 std::cout << tr1(vec1) << std::endl; 00047 std::cout << tr1.inv()(vec1) << std::endl; 00048 mln_assertion(tr1(tr1.inv()(vec1)) == vec1); 00049 mln_assertion(fabs(tr1(vec1)[0] - 3.9) <= 0.125 && 00050 fabs(tr1(vec1)[1] - 1.6) <= 0.125 && 00051 fabs(tr1(vec1)[2] - 4.5) <= 0.125); 00052 }