23       V v1(N), v2(N);  in test()  local
27       initialize_vector(v1);  in test()
29       ublas::blas_2::tmv(v1, m);  in test()
30       std::cout << "tmv (v1, m) = " << v1 << std::endl;  in test()
31       initialize_vector(v1);  in test()
33       ublas::blas_2::tmv(v1, ublas::trans(m));  in test()
34       std::cout << "tmv (v1, trans (m)) = " << v1 << std::endl;  in test()
36       initialize_vector(v1);  in test()
38       ublas::blas_2::tmv(v1, ublas::herm(m));  in test()
39       std::cout << "tmv (v1, herm (m)) = " << v1 << std::endl;  in test()
43       initialize_vector(v1);  in test()
46       ublas::blas_2::tsv(v1, m, ublas::lower_tag());  in test()
47       std::cout << "tsv (v1, m) = " << v1 << " " << ublas::prod(m, v1) - v2 << std::endl;  in test()
48       initialize_vector(v1);  in test()
51       ublas::blas_2::tsv(v1, ublas::trans(m), ublas::lower_tag());  in test()
52 …std::cout << "tsv (v1, trans (m)) = " << v1 << " " << ublas::prod(ublas::trans(m), v1) - v2 << std…  in test()
54       initialize_vector(v1);  in test()
57       ublas::blas_2::tsv(v1, ublas::herm(m), ublas::lower_tag());  in test()
58 …std::cout << "tsv (v1, herm (m)) = " << v1 << " " << ublas::prod(ublas::herm(m), v1) - v2 << std::…  in test()
60       initialize_vector(v1);  in test()
63       ublas::blas_2::tsv(v1, m, ublas::upper_tag());  in test()
64       std::cout << "tsv (v1, m) = " << v1 << " " << ublas::prod(m, v1) - v2 << std::endl;  in test()
65       initialize_vector(v1);  in test()
68       ublas::blas_2::tsv(v1, ublas::trans(m), ublas::upper_tag());  in test()
69 …std::cout << "tsv (v1, trans (m)) = " << v1 << " " << ublas::prod(ublas::trans(m), v1) - v2 << std…  in test()
71       initialize_vector(v1);  in test()
74       ublas::blas_2::tsv(v1, ublas::herm(m), ublas::upper_tag());  in test()
75 …std::cout << "tsv (v1, herm (m)) = " << v1 << " " << ublas::prod(ublas::herm(m), v1) - v2 << std::…  in test()
81       initialize_vector(v1);  in test()
84       ublas::blas_2::gmv(v1, value_type(1), value_type(1), m, v2);  in test()
85       std::cout << "gmv (v1, 1, 1, m, v2) = " << v1 << std::endl;  in test()
86       ublas::blas_2::gmv(v1, value_type(1), value_type(1), ublas::trans(m), v2);  in test()
87       std::cout << "gmv (v1, 1, 1, trans (m), v2) = " << v1 << std::endl;  in test()
89       ublas::blas_2::gmv(v1, value_type(1), value_type(1), ublas::herm(m), v2);  in test()
90       std::cout << "gmv (v1, 1, 1, herm (m), v2) = " << v1 << std::endl;  in test()
96       initialize_vector(v1);  in test()
99       ublas::blas_2::gr(m, value_type(1), v1, v2);  in test()
101       ublas::blas_2::gr(m, value_type(1), v1, ublas::conj(v2));  in test()
105       initialize_vector(v1);  in test()
107       ublas::blas_2::sr(m, value_type(1), v1);  in test()
112       initialize_vector(v1);  in test()
114       ublas::blas_2::hr(m, value_type(1), v1);  in test()
119       initialize_vector(v1);  in test()
122       ublas::blas_2::sr2(m, value_type(1), v1, v2);  in test()
127       initialize_vector(v1);  in test()
130       ublas::blas_2::hr2(m, value_type(1), v1, v2);  in test()