28 #include <mln/fun/x2x/translation.hh>
29 #include <mln/fun/i2v/all_to.hh>
42 algebra::vec<3,float> vec1 =
make::vec(a, b, c);
45 std::cout << vec1 << std::endl;
46 std::cout << tr1(vec1) << std::endl;
47 std::cout << tr1.inv()(vec1) << std::endl;
48 mln_assertion(tr1(tr1.inv()(vec1)) == vec1);
49 mln_assertion(fabs(tr1(vec1)[0] - 3.9) <= 0.125 &&
50 fabs(tr1(vec1)[1] - 1.6) <= 0.125 &&
51 fabs(tr1(vec1)[2] - 4.5) <= 0.125);