EQ2EMu/EQ2/source/depends/boost_1_72_0/libs/multiprecision/test/ublas_interop/test21.cpp

//
//  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, int N>
void test_blas_1<V, N>::test()
{
   {
      value_type t;
      real_type  n;
      V          v1(N), v2(N);

      // _asum
      initialize_vector(v1);
      n = ublas::blas_1::asum(v1);
      std::cout << "asum (v1) = " << n << std::endl;

      // _amax
      initialize_vector(v1);
      n = ublas::blas_1::amax(v1);
      std::cout << "amax (v1) = " << n << std::endl;

      // _nrm2
      initialize_vector(v1);
      n = ublas::blas_1::nrm2(v1);
      std::cout << "nrm2 (v1) = " << n << std::endl;

      // _dot
      // _dotu
      // _dotc
      initialize_vector(v1);
      initialize_vector(v2);
      t = ublas::blas_1::dot(v1, v2);
      std::cout << "dot (v1, v2) = " << t << std::endl;
      t = ublas::blas_1::dot(ublas::conj(v1), v2);
      std::cout << "dot (conj (v1), v2) = " << t << std::endl;

      // _copy
      initialize_vector(v2);
      ublas::blas_1::copy(v1, v2);
      std::cout << "copy (v1, v2) = " << v1 << std::endl;

      // _swap
      initialize_vector(v1);
      initialize_vector(v2);
      ublas::blas_1::swap(v1, v2);
      std::cout << "swap (v1, v2) = " << v1 << " " << v2 << std::endl;

      // _scal
      // csscal
      // zdscal
      initialize_vector(v1);
      ublas::blas_1::scal(v1, value_type(1));
      std::cout << "scal (v1, 1) = " << v1 << std::endl;

      // _axpy
      initialize_vector(v1);
      initialize_vector(v2);
      ublas::blas_1::axpy(v1, value_type(1), v2);
      std::cout << "axpy (v1, 1, v2) = " << v1 << std::endl;

      // _rot
      initialize_vector(v1);
      initialize_vector(v2);
      ublas::blas_1::rot(value_type(1), v1, value_type(1), v2);
      std::cout << "rot (1, v1, 1, v2) = " << v1 << " " << v2 << std::endl;
   }
}

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

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

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

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