//
//  Copyright (c) 2000-2002
//  Joerg Walter, Mathias Koch
//
//  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)
//
//  The authors gratefully acknowledge the support of
//  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);

      // _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;
#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;
#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;
#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;
#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;
#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;
#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;
#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;

      // _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;
#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;

#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;
#endif
   }
}

#ifdef USE_FLOAT
template struct test_blas_2<ublas::vector<mp_test_type>, ublas::matrix<mp_test_type>, 3>;
#endif

#ifdef USE_DOUBLE
template struct test_blas_2<ublas::vector<double>, ublas::matrix<double>, 3>;
#endif

#ifdef USE_STD_COMPLEX
#ifdef USE_FLOAT
template struct test_blas_2<ublas::vector<std::complex<mp_test_type> >, ublas::matrix<std::complex<mp_test_type> >, 3>;
#endif

#ifdef USE_DOUBLE
template struct test_blas_2<ublas::vector<std::complex<double> >, ublas::matrix<std::complex<double> >, 3>;
#endif
#endif
