1 /** test move semantics - run with and without BOOST_UBLAS_MOVE_SEMANTICS defined */
2
3 // Copyright Nasos Iliopoulos, Gunter Winkler 2009.
4 // Distributed under the Boost Software License, Version 1.0.
5 // (See accompanying file LICENSE_1_0.txt or copy at
6 // http://www.boost.org/LICENSE_1_0.txt)
7
8 #define BOOST_UBLAS_MOVE_SEMANTICS
9 #include <boost/numeric/ublas/vector.hpp>
10 #include <boost/numeric/ublas/matrix.hpp>
11 #include <boost/numeric/ublas/io.hpp>
12
13 namespace ublas= boost::numeric::ublas;
14 std::vector<double> a;
15
doubleit(ublas::vector<double> m)16 ublas::vector<double> doubleit(ublas::vector<double> m)
17 {
18 ublas::vector<double> r;
19 r=2.0*m;
20 std::cout << "Temporary pointer: " << &r[0] << std::endl;
21 return r;
22 }
23 template <class T,size_t N>
doubleit(ublas::bounded_vector<T,N> m)24 ublas::bounded_vector<T,N > doubleit(ublas::bounded_vector<T, N> m)
25 {
26 ublas::bounded_vector<T,N> r;
27 r=2.0*m;
28 std::cout << "Temporary pointer: " << &r[0] << std::endl;
29 return r;
30 }
31
32 template <class T,size_t N>
doubleit(ublas::c_vector<T,N> m)33 ublas::c_vector<T,N > doubleit(ublas::c_vector<T, N> m)
34 {
35 ublas::c_vector<T,N> r;
36 r=2.0*m;
37 std::cout << "Temporary pointer: " << &r[0] << std::endl;
38 return r;
39 }
40
doubleit(ublas::matrix<double> m)41 ublas::matrix<double> doubleit(ublas::matrix<double> m)
42 {
43 ublas::matrix<double> r;
44 r=2.0*m;
45 std::cout << "Temporary pointer r: " << &r(0,0) << std::endl;
46 return r;
47 }
48 template <class T,size_t N, size_t M>
doubleit(ublas::bounded_matrix<T,N,M> m)49 ublas::bounded_matrix<T,N, M > doubleit(ublas::bounded_matrix<T, N, M> m)
50 {
51 ublas::bounded_matrix<T,N, M> r;
52 r=2.0*m;
53 std::cout << "Temporary pointer: " << &(r(0,0)) << std::endl;
54 return r;
55 }
56
57 template <class T,size_t N, size_t M>
doubleit(ublas::c_matrix<T,N,M> m)58 ublas::c_matrix<T,N, M > doubleit(ublas::c_matrix<T, N, M> m)
59 {
60 ublas::c_matrix<T,N, M> r;
61 r=2.0*m;
62 std::cout << "Temporary pointer: " << &(r(0,0)) << std::endl;
63 return r;
64 }
65
66
test1()67 void test1()
68 {
69 std::cout << "vector<double> --------------------------------------------------------------------" << std::endl;
70 ublas::vector<double> a(ublas::scalar_vector<double>(2,2.0));
71 a = doubleit(a);
72 std::cout << "Pointer (must be equal to temp. pointer if move semantics are enabled) : " << &a[0] << std::endl;
73
74 std::cout << a << std::endl;
75
76 std::cout << "bounded_vector<double,2> --------------------------------------------------------------------" << std::endl;
77 ublas::bounded_vector<double,2> b(ublas::scalar_vector<double>(2,2.0));
78 ublas::bounded_vector<double,2> c;
79 noalias(c)=doubleit(b);
80 std::cout << "Pointer (bounded_vector swaps by copy this should be different than temp. pointer) : " << &c[0] << std::endl;
81 c(1)=0.0;
82 std::cout << b << std::endl;
83 std::cout << c << std::endl;
84
85 std::cout << "c_vector<double,2> --------------------------------------------------------------------" << std::endl;
86 ublas::c_vector<double,2> e=ublas::scalar_vector<double>(2,2.0);
87 ublas::c_vector<double,2> f;
88 f=doubleit(e);
89 std::cout << "Pointer (c_vector swaps by copy this should be different than temp. pointer) : " << &f[0] << std::endl;
90 f(1)=0;
91 std::cout << e << std::endl;
92 std::cout << f << std::endl;
93
94 }
95
test2()96 void test2() {
97 std::cout << "matrix<double> --------------------------------------------------------------------" << std::endl;
98 ublas::matrix<double> a(ublas::scalar_matrix<double>(2, 3, 2.0));
99 a = doubleit(a);
100 std::cout << "Pointer (must be equal to temp. pointer if move semantics are enabled) : " << &(a(0,0)) << std::endl;
101 std::cout << a << std::endl;
102
103 std::cout << "bounded_matrix<double,2, 3> --------------------------------------------------------------------" << std::endl;
104 ublas::bounded_matrix<double,2, 3> b(ublas::scalar_matrix<double>(2,3, 2.0));
105 ublas::bounded_matrix<double,2, 3> c;
106 noalias(c)=doubleit(b);
107 std::cout << "Pointer (bounded_matrix swaps by copy this should be different than temp. pointer) : " << &(c(0,0)) << std::endl;
108 c(1,1)=0.0;
109 std::cout << b << std::endl;
110 std::cout << c << std::endl;
111
112 std::cout << "c_matrix<double,2 ,3> --------------------------------------------------------------------" << std::endl;
113 ublas::c_matrix<double,2, 3> e=ublas::scalar_matrix<double>(2,3, 2.0);
114 ublas::c_matrix<double,2, 3> f;
115 f=doubleit(e);
116 std::cout << "Pointer (c_matrix swaps by copy this should be different than temp. pointer) : " << &(f(0,0)) << std::endl;
117 f(1,1)=0;
118 std::cout << e << std::endl;
119 std::cout << f << std::endl;
120 }
121
main()122 int main(){
123 test1();
124 test2();
125 return 0;
126 }
127
128