1 //
2 //  Copyright (c) 2000-2002
3 //  Joerg Walter, Mathias Koch
4 //
5 //  Distributed under the Boost Software License, Version 1.0. (See
6 //  accompanying file LICENSE_1_0.txt or copy at
7 //  http://www.boost.org/LICENSE_1_0.txt)
8 //
9 //  The authors gratefully acknowledge the support of
10 //  GeNeSys mbH & Co. KG in producing this work.
11 //
12 
13 #include "test2.hpp"
14 
15 template<class V, class M, std::size_t N>
test()16 void test_blas_2<V, M, N>::test () {
17     {
18         V v1 (N), v2 (N);
19         M m (N, N);
20 
21         // _t_mv
22         initialize_vector (v1);
23         initialize_matrix (m);
24         ublas::blas_2::tmv (v1, m);
25         std::cout << "tmv (v1, m) = " << v1 << std::endl;
26         initialize_vector (v1);
27         initialize_matrix (m);
28         ublas::blas_2::tmv (v1, ublas::trans (m));
29         std::cout << "tmv (v1, trans (m)) = " << v1 << std::endl;
30 #ifdef USE_STD_COMPLEX
31         initialize_vector (v1);
32         initialize_matrix (m);
33         ublas::blas_2::tmv (v1, ublas::herm (m));
34         std::cout << "tmv (v1, herm (m)) = " << v1 << std::endl;
35 #endif
36 
37         // _t_sv
38         initialize_vector (v1);
39         initialize_vector (v2);
40         initialize_matrix (m, ublas::lower_tag ());
41         ublas::blas_2::tsv (v1, m, ublas::lower_tag ());
42         std::cout << "tsv (v1, m) = " << v1 << " " << ublas::prod (m, v1) - v2 << std::endl;
43         initialize_vector (v1);
44         initialize_vector (v2);
45         initialize_matrix (m, ublas::upper_tag ());
46         ublas::blas_2::tsv (v1, ublas::trans (m), ublas::lower_tag ());
47         std::cout << "tsv (v1, trans (m)) = " << v1 << " " << ublas::prod (ublas::trans (m), v1) - v2 << std::endl;
48 #ifdef USE_STD_COMPLEX
49         initialize_vector (v1);
50         initialize_vector (v2);
51         initialize_matrix (m, ublas::upper_tag ());
52         ublas::blas_2::tsv (v1, ublas::herm (m), ublas::lower_tag ());
53         std::cout << "tsv (v1, herm (m)) = " << v1 << " " << ublas::prod (ublas::herm (m), v1) - v2 << std::endl;
54 #endif
55         initialize_vector (v1);
56         initialize_vector (v2);
57         initialize_matrix (m, ublas::upper_tag ());
58         ublas::blas_2::tsv (v1, m, ublas::upper_tag ());
59         std::cout << "tsv (v1, m) = " << v1 << " " << ublas::prod (m, v1) - v2 << std::endl;
60         initialize_vector (v1);
61         initialize_vector (v2);
62         initialize_matrix (m, ublas::lower_tag ());
63         ublas::blas_2::tsv (v1, ublas::trans (m), ublas::upper_tag ());
64         std::cout << "tsv (v1, trans (m)) = " << v1 << " " << ublas::prod (ublas::trans (m), v1) - v2 << std::endl;
65 #ifdef USE_STD_COMPLEX
66         initialize_vector (v1);
67         initialize_vector (v2);
68         initialize_matrix (m, ublas::lower_tag ());
69         ublas::blas_2::tsv (v1, ublas::herm (m), ublas::upper_tag ());
70         std::cout << "tsv (v1, herm (m)) = " << v1 << " " << ublas::prod (ublas::herm (m), v1) - v2 << std::endl;
71 #endif
72 
73         // _g_mv
74         // _s_mv
75         // _h_mv
76         initialize_vector (v1);
77         initialize_vector (v2);
78         initialize_matrix (m);
79         ublas::blas_2::gmv (v1, value_type (1), value_type (1), m, v2);
80         std::cout << "gmv (v1, 1, 1, m, v2) = " << v1 << std::endl;
81         ublas::blas_2::gmv (v1, value_type (1), value_type (1), ublas::trans (m), v2);
82         std::cout << "gmv (v1, 1, 1, trans (m), v2) = " << v1 << std::endl;
83 #ifdef USE_STD_COMPLEX
84         ublas::blas_2::gmv (v1, value_type (1), value_type (1), ublas::herm (m), v2);
85         std::cout << "gmv (v1, 1, 1, herm (m), v2) = " << v1 << std::endl;
86 #endif
87 
88         // _g_r
89         // _g_ru
90         // _g_rc
91         initialize_vector (v1);
92         initialize_vector (v2);
93         initialize_matrix (m);
94         ublas::blas_2::gr (m, value_type (1), v1, v2);
95         std::cout << "gr (m, 1, v1, v2) = " << m << std::endl;
96         ublas::blas_2::gr (m, value_type (1), v1, ublas::conj (v2));
97         std::cout << "gr (m, 1, v1, conj (v2)) = " << m << std::endl;
98 
99         // _s_r
100         initialize_vector (v1);
101         initialize_matrix (m);
102         ublas::blas_2::sr (m, value_type (1), v1);
103         std::cout << "sr (m, 1, v1) = " << m << std::endl;
104 
105 #ifdef USE_STD_COMPLEX
106         // _h_r
107         initialize_vector (v1);
108         initialize_matrix (m);
109         ublas::blas_2::hr (m, value_type (1), v1);
110         std::cout << "hr (m, 1, v1) = " << m << std::endl;
111 #endif
112 
113         // _s_r2
114         initialize_vector (v1);
115         initialize_vector (v2);
116         initialize_matrix (m);
117         ublas::blas_2::sr2 (m, value_type (1), v1, v2);
118         std::cout << "sr2 (m, 1, v1, v2) = " << m << std::endl;
119 
120 #ifdef USE_STD_COMPLEX
121         // _h_r2
122         initialize_vector (v1);
123         initialize_vector (v2);
124         initialize_matrix (m);
125         ublas::blas_2::hr2 (m, value_type (1), v1, v2);
126         std::cout << "hr2 (m, 1, v1, v2) = " << m << std::endl;
127 #endif
128     }
129 }
130 
131 #ifdef USE_FLOAT
132 template struct test_blas_2<ublas::vector<float>, ublas::matrix<float>, 3>;
133 #endif
134 
135 #ifdef USE_DOUBLE
136 template struct test_blas_2<ublas::vector<double>, ublas::matrix<double>, 3>;
137 #endif
138 
139 #ifdef USE_STD_COMPLEX
140 #ifdef USE_FLOAT
141 template struct test_blas_2<ublas::vector<std::complex<float> >, ublas::matrix<std::complex<float> >, 3>;
142 #endif
143 
144 #ifdef USE_DOUBLE
145 template struct test_blas_2<ublas::vector<std::complex<double> >, ublas::matrix<std::complex<double> >, 3>;
146 #endif
147 #endif
148