blas_2.cxx

Multivac 的Level set包

  template <class Prop0, class Allocator0,
	    class Allocator1>
  void
  Mlt(const enum CBLAS_TRANSPOSE TransA,
      const enum CBLAS_DIAG DiagA,
      const Matrix<complex<float>, Prop0, RowLoTriang, Allocator0>& A,
      Vector<complex<float>, Vect_Full, Allocator1>& X)
  {
    cblas_ctrmv(CblasRowMajor, CblasLower, TransA, DiagA,
		A.GetN(), reinterpret_cast<const void*>(A.GetData()),
		A.GetM(), reinterpret_cast<void*>(X.GetData()), 1);
  }


  template <class Prop0, class Allocator0,
	    class Allocator1>
  void
  Mlt(const enum CBLAS_TRANSPOSE TransA,
      const enum CBLAS_DIAG DiagA,
      const Matrix<complex<double>, Prop0, RowLoTriang, Allocator0>& A,
      Vector<complex<double>, Vect_Full, Allocator1>& X)
  {
    cblas_ztrmv(CblasRowMajor, CblasLower, TransA, DiagA,
		A.GetN(), reinterpret_cast<const void*>(A.GetData()),
		A.GetM(), reinterpret_cast<void*>(X.GetData()), 1);
  }


  /*** ColUpTriangPacked, NoTrans and NonUnit ***/


  template <class Prop0, class Allocator0,
	    class Allocator1>
  void Mlt(const Matrix<float, Prop0, ColUpTriangPacked, Allocator0>& A,
	   Vector<float, Vect_Full, Allocator1>& X)
  {
    cblas_stpmv(CblasColMajor, CblasUpper, CblasNoTrans, CblasNonUnit,
		A.GetN(), A.GetData(), X.GetData(), 1);
  }


  template <class Prop0, class Allocator0,
	    class Allocator1>
  void Mlt(const Matrix<double, Prop0, ColUpTriangPacked, Allocator0>& A,
	   Vector<double, Vect_Full, Allocator1>& X)
  {
    cblas_dtpmv(CblasColMajor, CblasUpper, CblasNoTrans, CblasNonUnit,
		A.GetN(), A.GetData(), X.GetData(), 1);
  }


  template <class Prop0, class Allocator0,
	    class Allocator1>
  void
  Mlt(const Matrix<complex<float>, Prop0, ColUpTriangPacked, Allocator0>& A,
      Vector<complex<float>, Vect_Full, Allocator1>& X)
  {
    cblas_ctpmv(CblasColMajor, CblasUpper, CblasNoTrans, CblasNonUnit,
		A.GetN(), reinterpret_cast<const void*>(A.GetData()),
		reinterpret_cast<void*>(X.GetData()), 1);
  }


  template <class Prop0, class Allocator0,
	    class Allocator1>
  void
  Mlt(const Matrix<complex<double>, Prop0, ColUpTriangPacked, Allocator0>& A,
      Vector<complex<double>, Vect_Full, Allocator1>& X)
  {
    cblas_ztpmv(CblasColMajor, CblasUpper, CblasNoTrans, CblasNonUnit,
		A.GetN(), reinterpret_cast<const void*>(A.GetData()),
		reinterpret_cast<void*>(X.GetData()), 1);
  }


  /*** ColUpTriangPacked ***/


  template <class Prop0, class Allocator0,
	    class Allocator1>
  void Mlt(const enum CBLAS_TRANSPOSE TransA,
	   const enum CBLAS_DIAG DiagA,
	   const Matrix<float, Prop0, ColUpTriangPacked, Allocator0>& A,
	   Vector<float, Vect_Full, Allocator1>& X)
  {
    cblas_stpmv(CblasColMajor, CblasUpper, TransA, DiagA,
		A.GetN(), A.GetData(), X.GetData(), 1);
  }


  template <class Prop0, class Allocator0,
	    class Allocator1>
  void Mlt(const enum CBLAS_TRANSPOSE TransA,
	   const enum CBLAS_DIAG DiagA,
	   const Matrix<double, Prop0, ColUpTriangPacked, Allocator0>& A,
	   Vector<double, Vect_Full, Allocator1>& X)
  {
    cblas_dtpmv(CblasColMajor, CblasUpper, TransA, DiagA,
		A.GetN(), A.GetData(), X.GetData(), 1);
  }


  template <class Prop0, class Allocator0,
	    class Allocator1>
  void
  Mlt(const enum CBLAS_TRANSPOSE TransA,
      const enum CBLAS_DIAG DiagA,
      const Matrix<complex<float>, Prop0, ColUpTriangPacked, Allocator0>& A,
      Vector<complex<float>, Vect_Full, Allocator1>& X)
  {
    cblas_ctpmv(CblasColMajor, CblasUpper, TransA, DiagA,
		A.GetN(), reinterpret_cast<const void*>(A.GetData()),
		reinterpret_cast<void*>(X.GetData()), 1);
  }


  template <class Prop0, class Allocator0,
	    class Allocator1>
  void
  Mlt(const enum CBLAS_TRANSPOSE TransA,
      const enum CBLAS_DIAG DiagA,
      const Matrix<complex<double>, Prop0, ColUpTriangPacked, Allocator0>& A,
      Vector<complex<double>, Vect_Full, Allocator1>& X)
  {
    cblas_ztpmv(CblasColMajor, CblasUpper, TransA, DiagA,
		A.GetN(), reinterpret_cast<const void*>(A.GetData()),
		reinterpret_cast<void*>(X.GetData()), 1);
  }


  /*** ColLoTriangPacked, NoTrans and NonUnit ***/


  template <class Prop0, class Allocator0,
	    class Allocator1>
  void Mlt(const Matrix<float, Prop0, ColLoTriangPacked, Allocator0>& A,
	   Vector<float, Vect_Full, Allocator1>& X)
  {
    cblas_stpmv(CblasColMajor, CblasLower, CblasNoTrans, CblasNonUnit,
		A.GetN(), A.GetData(), X.GetData(), 1);
  }


  template <class Prop0, class Allocator0,
	    class Allocator1>
  void Mlt(const Matrix<double, Prop0, ColLoTriangPacked, Allocator0>& A,
	   Vector<double, Vect_Full, Allocator1>& X)
  {
    cblas_dtpmv(CblasColMajor, CblasLower, CblasNoTrans, CblasNonUnit,
		A.GetN(), A.GetData(), X.GetData(), 1);
  }


  template <class Prop0, class Allocator0,
	    class Allocator1>
  void
  Mlt(const Matrix<complex<float>, Prop0, ColLoTriangPacked, Allocator0>& A,
      Vector<complex<float>, Vect_Full, Allocator1>& X)
  {
    cblas_ctpmv(CblasColMajor, CblasLower, CblasNoTrans, CblasNonUnit,
		A.GetN(), reinterpret_cast<const void*>(A.GetData()),
		reinterpret_cast<void*>(X.GetData()), 1);
  }


  template <class Prop0, class Allocator0,
	    class Allocator1>
  void
  Mlt(const Matrix<complex<double>, Prop0, ColLoTriangPacked, Allocator0>& A,
      Vector<complex<double>, Vect_Full, Allocator1>& X)
  {
    cblas_ztpmv(CblasColMajor, CblasLower, CblasNoTrans, CblasNonUnit,
		A.GetN(), reinterpret_cast<const void*>(A.GetData()),
		reinterpret_cast<void*>(X.GetData()), 1);
  }


  /*** ColLoTriangPacked ***/


  template <class Prop0, class Allocator0,
	    class Allocator1>
  void Mlt(const enum CBLAS_TRANSPOSE TransA,
	   const enum CBLAS_DIAG DiagA,
	   const Matrix<float, Prop0, ColLoTriangPacked, Allocator0>& A,
	   Vector<float, Vect_Full, Allocator1>& X)
  {
    cblas_stpmv(CblasColMajor, CblasLower, TransA, DiagA,
		A.GetN(), A.GetData(), X.GetData(), 1);
  }


  template <class Prop0, class Allocator0,
	    class Allocator1>
  void Mlt(const enum CBLAS_TRANSPOSE TransA,
	   const enum CBLAS_DIAG DiagA,
	   const Matrix<double, Prop0, ColLoTriangPacked, Allocator0>& A,
	   Vector<double, Vect_Full, Allocator1>& X)
  {
    cblas_dtpmv(CblasColMajor, CblasLower, TransA, DiagA,
		A.GetN(), A.GetData(), X.GetData(), 1);
  }


  template <class Prop0, class Allocator0,
	    class Allocator1>
  void
  Mlt(const enum CBLAS_TRANSPOSE TransA,
      const enum CBLAS_DIAG DiagA,
      const Matrix<complex<float>, Prop0, ColLoTriangPacked, Allocator0>& A,
      Vector<complex<float>, Vect_Full, Allocator1>& X)
  {
    cblas_ctpmv(CblasColMajor, CblasLower, TransA, DiagA,
		A.GetN(), reinterpret_cast<const void*>(A.GetData()),
		reinterpret_cast<void*>(X.GetData()), 1);
  }


  template <class Prop0, class Allocator0,
	    class Allocator1>
  void
  Mlt(const enum CBLAS_TRANSPOSE TransA,
      const enum CBLAS_DIAG DiagA,
      const Matrix<complex<double>, Prop0, ColLoTriangPacked, Allocator0>& A,
      Vector<complex<double>, Vect_Full, Allocator1>& X)
  {
    cblas_ztpmv(CblasColMajor, CblasLower, TransA, DiagA,
		A.GetN(), reinterpret_cast<const void*>(A.GetData()),
		reinterpret_cast<void*>(X.GetData()), 1);
  }


  /*** RowUpTriangPacked, NoTrans and NonUnit ***/


  template <class Prop0, class Allocator0,
	    class Allocator1>
  void Mlt(const Matrix<float, Prop0, RowUpTriangPacked, Allocator0>& A,
	   Vector<float, Vect_Full, Allocator1>& X)
  {
    cblas_stpmv(CblasRowMajor, CblasUpper, CblasNoTrans, CblasNonUnit,
		A.GetN(), A.GetData(), X.GetData(), 1);
  }


  template <class Prop0, class Allocator0,
	    class Allocator1>
  void Mlt(const Matrix<double, Prop0, RowUpTriangPacked, Allocator0>& A,
	   Vector<double, Vect_Full, Allocator1>& X)
  {
    cblas_dtpmv(CblasRowMajor, CblasUpper, CblasNoTrans, CblasNonUnit,
		A.GetN(), A.GetData(), X.GetData(), 1);
  }


  template <class Prop0, class Allocator0,
	    class Allocator1>
  void
  Mlt(const Matrix<complex<float>, Prop0, RowUpTriangPacked, Allocator0>& A,
      Vector<complex<float>, Vect_Full, Allocator1>& X)
  {
    cblas_ctpmv(CblasRowMajor, CblasUpper, CblasNoTrans, CblasNonUnit,
		A.GetN(), reinterpret_cast<const void*>(A.GetData()),
		reinterpret_cast<void*>(X.GetData()), 1);
  }


  template <class Prop0, class Allocator0,
	    class Allocator1>
  void
  Mlt(const Matrix<complex<double>, Prop0, RowUpTriangPacked, Allocator0>& A,
      Vector<complex<double>, Vect_Full, Allocator1>& X)
  {
    cblas_ztpmv(CblasRowMajor, CblasUpper, CblasNoTrans, CblasNonUnit,
		A.GetN(), reinterpret_cast<const void*>(A.GetData()),
		reinterpret_cast<void*>(X.GetData()), 1);
  }


  /*** RowUpTriangPacked ***/


  template <class Prop0, class Allocator0,
	    class Allocator1>
  void Mlt(const enum CBLAS_TRANSPOSE TransA,
	   const enum CBLAS_DIAG DiagA,
	   const Matrix<float, Prop0, RowUpTriangPacked, Allocator0>& A,
	   Vector<float, Vect_Full, Allocator1>& X)
  {
    cblas_stpmv(CblasRowMajor, CblasUpper, TransA, DiagA,
		A.GetN(), A.GetData(), X.GetData(), 1);
  }


  template <class Prop0, class Allocator0,
	    class Allocator1>
  void Mlt(const enum CBLAS_TRANSPOSE TransA,
	   const enum CBLAS_DIAG DiagA,
	   const Matrix<double, Prop0, RowUpTriangPacked, Allocator0>& A,
	   Vector<double, Vect_Full, Allocator1>& X)
  {
    cblas_dtpmv(CblasRowMajor, CblasUpper, TransA, DiagA,
		A.GetN(), A.GetData(), X.GetData(), 1);
  }


  template <class Prop0, class Allocator0,
	    class Allocator1>
  void
  Mlt(const enum CBLAS_TRANSPOSE TransA,
      const enum CBLAS_DIAG DiagA,
      const Matrix<complex<float>, Prop0, RowUpTriangPacked, Allocator0>& A,
      Vector<complex<float>, Vect_Full, Allocator1>& X)
  {
    cblas_ctpmv(CblasRowMajor, CblasUpper, TransA, DiagA,
		A.GetN(), reinterpret_cast<const void*>(A.GetData()),
		reinterpret_cast<void*>(X.GetData()), 1);
  }


  template <class Prop0, class Allocator0,
	    class Allocator1>
  void
  Mlt(const enum CBLAS_TRANSPOSE TransA,
      const enum CBLAS_DIAG DiagA,
      const Matrix<complex<double>, Prop0, RowUpTriangPacked, Allocator0>& A,
      Vector<complex<double>, Vect_Full, Allocator1>& X)
  {
    cblas_ztpmv(CblasRowMajor, CblasUpper, TransA, DiagA,
		A.GetN(), reinterpret_cast<const void*>(A.GetData()),
		reinterpret_cast<void*>(X.GetData()), 1);
  }


  /*** RowLoTriangPacked, NoTrans and NonUnit ***/


  template <class Prop0, class Allocator0,
	    class Allocator1>
  void Mlt(const Matrix<float, Prop0, RowLoTriangPacked, Allocator0>& A,
	   Vector<float, Vect_Full, Allocator1>& X)
  {
    cblas_stpmv(CblasRowMajor, CblasLower, CblasNoTrans, CblasNonUnit,
		A.GetN(), A.GetData(), X.GetData(), 1);
  }


  template <class Prop0, class Allocator0,
	    class Allocator1>
  void Mlt(const Matrix<double, Prop0, RowLoTriangPacked, Allocator0>& A,
	   Vector<double, Vect_Full, Allocator1>& X)
  {
    cblas_dtpmv(CblasRowMajor, CblasLower, CblasNoTrans, CblasNonUnit,
		A.GetN(), A.GetData(), X.GetData(), 1);
  }


  template <class Prop0, class Allocator0,
	    class Allocator1>
  void
  Mlt(const Matrix<complex<float>, Prop0, RowLoTriangPacked, Allocator0>& A,
      Vector<complex<float>, Vect_Full, Allocator1>& X)
  {
    cblas_ctpmv(CblasRowMajor, CblasLower, CblasNoTrans, CblasNonUnit,
		A.GetN(), reinterpret_cast<const void*>(A.GetData()),
		reinterpret_cast<void*>(X.GetData()), 1);
  }


  template <class Prop0, class Allocator0,
	    class Allocator1>
  void
  Mlt(const Matrix<complex<double>, Prop0, RowLoTriangPacked, Allocator0>& A,
      Vector<complex<double>, Vect_Full, Allocator1>& X)
  {
    cblas_ztpmv(CblasRowMajor, CblasLower, CblasNoTrans, CblasNonUnit,
		A.GetN(), reinterpret_cast<const void*>(A.GetData()),
		reinterpret_cast<void*>(X.GetData()), 1);
  }


  /*** RowLoTriangPacked ***/


  template <class Prop0, class Allocator0,
	    class Allocator1>
  void Mlt(const enum CBLAS_TRANSPOSE TransA,
	   const enum CBLAS_DIAG DiagA,
	   const Matrix<float, Prop0, RowLoTriangPacked, Allocator0>& A,
	   Vector<float, Vect_Full, Allocator1>& X)
  {
    cblas_stpmv(CblasRowMajor, CblasLower, TransA, DiagA,
		A.GetN(), A.GetData(), X.GetData(), 1);
  }


  template <class Prop0, class Allocator0,
	    class Allocator1>
  void Mlt(const enum CBLAS_TRANSPOSE TransA,
	   const enum CBLAS_DIAG DiagA,
	   const Matrix<double, Prop0, RowLoTriangPacked, Allocator0>& A,
	   Vector<double, Vect_Full, Allocator1>& X)
  {
    cblas_dtpmv(CblasRowMajor, CblasLower, TransA, DiagA,
		A.GetN(), A.GetData(), X.GetData(), 1);
  }


  template <class Prop0, class Allocator0,
	    class Allocator1>
  void
  Mlt(const enum CBLAS_TRANSPOSE TransA,
      const enum CBLAS_DIAG DiagA,
      const Matrix<complex<float>, Prop0, RowLoTriangPacked, Allocator0>& A,
      Vector<complex<float>, Vect_Full, Allocator1>& X)
  {