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