Imported Upstream version 1.72.0
[platform/upstream/boost.git] / libs / multiprecision / test / ublas_interop / test22.cpp
index 297bb06..ce3d803 100644 (file)
 //  GeNeSys mbH & Co. KG in producing this work.
 //
 
+#if defined(__GNUC__) && (__GNUC__ >= 9)
+#pragma GCC diagnostic ignored "-Wdeprecated-copy"
+#endif
+
 #include "test2.hpp"
 
-template<class V, class M, int N>
-void test_blas_2<V, M, N>::test () {
-    {
-        V v1 (N), v2 (N);
-        M m (N, N);
+template <class V, class M, int N>
+void test_blas_2<V, M, N>::test()
+{
+   {
+      V v1(N), v2(N);
+      M m(N, N);
 
-        // _t_mv
-        initialize_vector (v1);
-        initialize_matrix (m);
-        ublas::blas_2::tmv (v1, m);
-        std::cout << "tmv (v1, m) = " << v1 << std::endl;
-        initialize_vector (v1);
-        initialize_matrix (m);
-        ublas::blas_2::tmv (v1, ublas::trans (m));
-        std::cout << "tmv (v1, trans (m)) = " << v1 << std::endl;
+      // _t_mv
+      initialize_vector(v1);
+      initialize_matrix(m);
+      ublas::blas_2::tmv(v1, m);
+      std::cout << "tmv (v1, m) = " << v1 << std::endl;
+      initialize_vector(v1);
+      initialize_matrix(m);
+      ublas::blas_2::tmv(v1, ublas::trans(m));
+      std::cout << "tmv (v1, trans (m)) = " << v1 << std::endl;
 #ifdef USE_STD_COMPLEX
-        initialize_vector (v1);
-        initialize_matrix (m);
-        ublas::blas_2::tmv (v1, ublas::herm (m));
-        std::cout << "tmv (v1, herm (m)) = " << v1 << std::endl;
+      initialize_vector(v1);
+      initialize_matrix(m);
+      ublas::blas_2::tmv(v1, ublas::herm(m));
+      std::cout << "tmv (v1, herm (m)) = " << v1 << std::endl;
 #endif
 
-        // _t_sv
-        initialize_vector (v1);
-        initialize_vector (v2);
-        initialize_matrix (m, ublas::lower_tag ());
-        ublas::blas_2::tsv (v1, m, ublas::lower_tag ());
-        std::cout << "tsv (v1, m) = " << v1 << " " << ublas::prod (m, v1) - v2 << std::endl;
-        initialize_vector (v1);
-        initialize_vector (v2);
-        initialize_matrix (m, ublas::upper_tag ());
-        ublas::blas_2::tsv (v1, ublas::trans (m), ublas::lower_tag ());
-        std::cout << "tsv (v1, trans (m)) = " << v1 << " " << ublas::prod (ublas::trans (m), v1) - v2 << std::endl;
+      // _t_sv
+      initialize_vector(v1);
+      initialize_vector(v2);
+      initialize_matrix(m, ublas::lower_tag());
+      ublas::blas_2::tsv(v1, m, ublas::lower_tag());
+      std::cout << "tsv (v1, m) = " << v1 << " " << ublas::prod(m, v1) - v2 << std::endl;
+      initialize_vector(v1);
+      initialize_vector(v2);
+      initialize_matrix(m, ublas::upper_tag());
+      ublas::blas_2::tsv(v1, ublas::trans(m), ublas::lower_tag());
+      std::cout << "tsv (v1, trans (m)) = " << v1 << " " << ublas::prod(ublas::trans(m), v1) - v2 << std::endl;
 #ifdef USE_STD_COMPLEX
-        initialize_vector (v1);
-        initialize_vector (v2);
-        initialize_matrix (m, ublas::upper_tag ());
-        ublas::blas_2::tsv (v1, ublas::herm (m), ublas::lower_tag ());
-        std::cout << "tsv (v1, herm (m)) = " << v1 << " " << ublas::prod (ublas::herm (m), v1) - v2 << std::endl;
+      initialize_vector(v1);
+      initialize_vector(v2);
+      initialize_matrix(m, ublas::upper_tag());
+      ublas::blas_2::tsv(v1, ublas::herm(m), ublas::lower_tag());
+      std::cout << "tsv (v1, herm (m)) = " << v1 << " " << ublas::prod(ublas::herm(m), v1) - v2 << std::endl;
 #endif
-        initialize_vector (v1);
-        initialize_vector (v2);
-        initialize_matrix (m, ublas::upper_tag ());
-        ublas::blas_2::tsv (v1, m, ublas::upper_tag ());
-        std::cout << "tsv (v1, m) = " << v1 << " " << ublas::prod (m, v1) - v2 << std::endl;
-        initialize_vector (v1);
-        initialize_vector (v2);
-        initialize_matrix (m, ublas::lower_tag ());
-        ublas::blas_2::tsv (v1, ublas::trans (m), ublas::upper_tag ());
-        std::cout << "tsv (v1, trans (m)) = " << v1 << " " << ublas::prod (ublas::trans (m), v1) - v2 << std::endl;
+      initialize_vector(v1);
+      initialize_vector(v2);
+      initialize_matrix(m, ublas::upper_tag());
+      ublas::blas_2::tsv(v1, m, ublas::upper_tag());
+      std::cout << "tsv (v1, m) = " << v1 << " " << ublas::prod(m, v1) - v2 << std::endl;
+      initialize_vector(v1);
+      initialize_vector(v2);
+      initialize_matrix(m, ublas::lower_tag());
+      ublas::blas_2::tsv(v1, ublas::trans(m), ublas::upper_tag());
+      std::cout << "tsv (v1, trans (m)) = " << v1 << " " << ublas::prod(ublas::trans(m), v1) - v2 << std::endl;
 #ifdef USE_STD_COMPLEX
-        initialize_vector (v1);
-        initialize_vector (v2);
-        initialize_matrix (m, ublas::lower_tag ());
-        ublas::blas_2::tsv (v1, ublas::herm (m), ublas::upper_tag ());
-        std::cout << "tsv (v1, herm (m)) = " << v1 << " " << ublas::prod (ublas::herm (m), v1) - v2 << std::endl;
+      initialize_vector(v1);
+      initialize_vector(v2);
+      initialize_matrix(m, ublas::lower_tag());
+      ublas::blas_2::tsv(v1, ublas::herm(m), ublas::upper_tag());
+      std::cout << "tsv (v1, herm (m)) = " << v1 << " " << ublas::prod(ublas::herm(m), v1) - v2 << std::endl;
 #endif
 
-        // _g_mv
-        // _s_mv
-        // _h_mv
-        initialize_vector (v1);
-        initialize_vector (v2);
-        initialize_matrix (m);
-        ublas::blas_2::gmv (v1, value_type (1), value_type (1), m, v2);
-        std::cout << "gmv (v1, 1, 1, m, v2) = " << v1 << std::endl;
-        ublas::blas_2::gmv (v1, value_type (1), value_type (1), ublas::trans (m), v2);
-        std::cout << "gmv (v1, 1, 1, trans (m), v2) = " << v1 << std::endl;
+      // _g_mv
+      // _s_mv
+      // _h_mv
+      initialize_vector(v1);
+      initialize_vector(v2);
+      initialize_matrix(m);
+      ublas::blas_2::gmv(v1, value_type(1), value_type(1), m, v2);
+      std::cout << "gmv (v1, 1, 1, m, v2) = " << v1 << std::endl;
+      ublas::blas_2::gmv(v1, value_type(1), value_type(1), ublas::trans(m), v2);
+      std::cout << "gmv (v1, 1, 1, trans (m), v2) = " << v1 << std::endl;
 #ifdef USE_STD_COMPLEX
-        ublas::blas_2::gmv (v1, value_type (1), value_type (1), ublas::herm (m), v2);
-        std::cout << "gmv (v1, 1, 1, herm (m), v2) = " << v1 << std::endl;
+      ublas::blas_2::gmv(v1, value_type(1), value_type(1), ublas::herm(m), v2);
+      std::cout << "gmv (v1, 1, 1, herm (m), v2) = " << v1 << std::endl;
 #endif
 
-        // _g_r
-        // _g_ru
-        // _g_rc
-        initialize_vector (v1);
-        initialize_vector (v2);
-        initialize_matrix (m);
-        ublas::blas_2::gr (m, value_type (1), v1, v2);
-        std::cout << "gr (m, 1, v1, v2) = " << m << std::endl;
-        ublas::blas_2::gr (m, value_type (1), v1, ublas::conj (v2));
-        std::cout << "gr (m, 1, v1, conj (v2)) = " << m << std::endl;
+      // _g_r
+      // _g_ru
+      // _g_rc
+      initialize_vector(v1);
+      initialize_vector(v2);
+      initialize_matrix(m);
+      ublas::blas_2::gr(m, value_type(1), v1, v2);
+      std::cout << "gr (m, 1, v1, v2) = " << m << std::endl;
+      ublas::blas_2::gr(m, value_type(1), v1, ublas::conj(v2));
+      std::cout << "gr (m, 1, v1, conj (v2)) = " << m << std::endl;
 
-        // _s_r
-        initialize_vector (v1);
-        initialize_matrix (m);
-        ublas::blas_2::sr (m, value_type (1), v1);
-        std::cout << "sr (m, 1, v1) = " << m << std::endl;
+      // _s_r
+      initialize_vector(v1);
+      initialize_matrix(m);
+      ublas::blas_2::sr(m, value_type(1), v1);
+      std::cout << "sr (m, 1, v1) = " << m << std::endl;
 
 #ifdef USE_STD_COMPLEX
-        // _h_r
-        initialize_vector (v1);
-        initialize_matrix (m);
-        ublas::blas_2::hr (m, value_type (1), v1);
-        std::cout << "hr (m, 1, v1) = " << m << std::endl;
+      // _h_r
+      initialize_vector(v1);
+      initialize_matrix(m);
+      ublas::blas_2::hr(m, value_type(1), v1);
+      std::cout << "hr (m, 1, v1) = " << m << std::endl;
 #endif
 
-        // _s_r2
-        initialize_vector (v1);
-        initialize_vector (v2);
-        initialize_matrix (m);
-        ublas::blas_2::sr2 (m, value_type (1), v1, v2);
-        std::cout << "sr2 (m, 1, v1, v2) = " << m << std::endl;
+      // _s_r2
+      initialize_vector(v1);
+      initialize_vector(v2);
+      initialize_matrix(m);
+      ublas::blas_2::sr2(m, value_type(1), v1, v2);
+      std::cout << "sr2 (m, 1, v1, v2) = " << m << std::endl;
 
 #ifdef USE_STD_COMPLEX
-        // _h_r2
-        initialize_vector (v1);
-        initialize_vector (v2);
-        initialize_matrix (m);
-        ublas::blas_2::hr2 (m, value_type (1), v1, v2);
-        std::cout << "hr2 (m, 1, v1, v2) = " << m << std::endl;
+      // _h_r2
+      initialize_vector(v1);
+      initialize_vector(v2);
+      initialize_matrix(m);
+      ublas::blas_2::hr2(m, value_type(1), v1, v2);
+      std::cout << "hr2 (m, 1, v1, v2) = " << m << std::endl;
 #endif
-    }
+   }
 }
 
 #ifdef USE_FLOAT