|
|
Boost-Commit : |
Subject: [Boost-commit] svn:boost r79908 - in sandbox/numeric_bindings/boost/numeric/bindings/blas: level1 level2 level3
From: rutger_at_[hidden]
Date: 2012-08-07 09:40:40
Author: rutger
Date: 2012-08-07 09:40:31 EDT (Tue, 07 Aug 2012)
New Revision: 79908
URL: http://svn.boost.org/trac/boost/changeset/79908
Log:
Update of BLAS bindings (result of updated translator/generator)
Text files modified:
sandbox/numeric_bindings/boost/numeric/bindings/blas/level1/asum.hpp | 2
sandbox/numeric_bindings/boost/numeric/bindings/blas/level1/axpy.hpp | 2
sandbox/numeric_bindings/boost/numeric/bindings/blas/level1/copy.hpp | 2
sandbox/numeric_bindings/boost/numeric/bindings/blas/level1/dotc.hpp | 2
sandbox/numeric_bindings/boost/numeric/bindings/blas/level1/iamax.hpp | 2
sandbox/numeric_bindings/boost/numeric/bindings/blas/level1/nrm2.hpp | 2
sandbox/numeric_bindings/boost/numeric/bindings/blas/level1/rot.hpp | 2
sandbox/numeric_bindings/boost/numeric/bindings/blas/level1/rotg.hpp | 4
sandbox/numeric_bindings/boost/numeric/bindings/blas/level1/scal.hpp | 2
sandbox/numeric_bindings/boost/numeric/bindings/blas/level1/swap.hpp | 2
sandbox/numeric_bindings/boost/numeric/bindings/blas/level2/gbmv.hpp | 78 +++++++++++++++++----------------
sandbox/numeric_bindings/boost/numeric/bindings/blas/level2/gemv.hpp | 64 ++++++++++++---------------
sandbox/numeric_bindings/boost/numeric/bindings/blas/level2/ger.hpp | 40 ++++++++--------
sandbox/numeric_bindings/boost/numeric/bindings/blas/level2/gerc.hpp | 30 ++++++------
sandbox/numeric_bindings/boost/numeric/bindings/blas/level2/geru.hpp | 30 ++++++------
sandbox/numeric_bindings/boost/numeric/bindings/blas/level2/hbmv.hpp | 70 +++++++++++++----------------
sandbox/numeric_bindings/boost/numeric/bindings/blas/level2/hemv.hpp | 46 ++++++++++---------
sandbox/numeric_bindings/boost/numeric/bindings/blas/level2/her.hpp | 50 +++++++++-----------
sandbox/numeric_bindings/boost/numeric/bindings/blas/level2/her2.hpp | 46 ++++++++++---------
sandbox/numeric_bindings/boost/numeric/bindings/blas/level2/hpmv.hpp | 46 ++++++++++---------
sandbox/numeric_bindings/boost/numeric/bindings/blas/level2/hpr.hpp | 44 +++++++++---------
sandbox/numeric_bindings/boost/numeric/bindings/blas/level2/hpr2.hpp | 44 +++++++++---------
sandbox/numeric_bindings/boost/numeric/bindings/blas/level2/sbmv.hpp | 35 ++++++--------
sandbox/numeric_bindings/boost/numeric/bindings/blas/level2/spmv.hpp | 31 ++++++------
sandbox/numeric_bindings/boost/numeric/bindings/blas/level2/spr.hpp | 22 ++++----
sandbox/numeric_bindings/boost/numeric/bindings/blas/level2/spr2.hpp | 30 ++++++------
sandbox/numeric_bindings/boost/numeric/bindings/blas/level2/symv.hpp | 31 ++++++------
sandbox/numeric_bindings/boost/numeric/bindings/blas/level2/syr.hpp | 24 ++++------
sandbox/numeric_bindings/boost/numeric/bindings/blas/level2/syr2.hpp | 31 ++++++------
sandbox/numeric_bindings/boost/numeric/bindings/blas/level2/tbmv.hpp | 80 ++++++++++++++++------------------
sandbox/numeric_bindings/boost/numeric/bindings/blas/level2/tbsv.hpp | 80 ++++++++++++++++------------------
sandbox/numeric_bindings/boost/numeric/bindings/blas/level2/tpmv.hpp | 69 ++++++++++++++---------------
sandbox/numeric_bindings/boost/numeric/bindings/blas/level2/tpsv.hpp | 69 ++++++++++++++---------------
sandbox/numeric_bindings/boost/numeric/bindings/blas/level2/trmv.hpp | 73 ++++++++++++++++---------------
sandbox/numeric_bindings/boost/numeric/bindings/blas/level2/trsv.hpp | 67 ++++++++++++++--------------
sandbox/numeric_bindings/boost/numeric/bindings/blas/level3/gemm.hpp | 72 +++++++++++++++---------------
sandbox/numeric_bindings/boost/numeric/bindings/blas/level3/hemm.hpp | 93 +++++++++++++++++++--------------------
sandbox/numeric_bindings/boost/numeric/bindings/blas/level3/her2k.hpp | 59 ++++++++++++------------
sandbox/numeric_bindings/boost/numeric/bindings/blas/level3/herk.hpp | 81 +++++++++++++++++-----------------
sandbox/numeric_bindings/boost/numeric/bindings/blas/level3/symm.hpp | 93 +++++++++++++++++++--------------------
sandbox/numeric_bindings/boost/numeric/bindings/blas/level3/syr2k.hpp | 59 ++++++++++++------------
sandbox/numeric_bindings/boost/numeric/bindings/blas/level3/syrk.hpp | 74 +++++++++++++++---------------
sandbox/numeric_bindings/boost/numeric/bindings/blas/level3/trmm.hpp | 92 +++++++++++++++++++--------------------
sandbox/numeric_bindings/boost/numeric/bindings/blas/level3/trsm.hpp | 88 ++++++++++++++++++-------------------
44 files changed, 967 insertions(+), 996 deletions(-)
Modified: sandbox/numeric_bindings/boost/numeric/bindings/blas/level1/asum.hpp
==============================================================================
--- sandbox/numeric_bindings/boost/numeric/bindings/blas/level1/asum.hpp (original)
+++ sandbox/numeric_bindings/boost/numeric/bindings/blas/level1/asum.hpp 2012-08-07 09:40:31 EDT (Tue, 07 Aug 2012)
@@ -129,7 +129,7 @@
//
inline double asum( const int n, const std::complex<double>* x,
const int incx ) {
- return // NOT FOUND();
+ return cublasDzasum( n, x, incx );
}
#else
Modified: sandbox/numeric_bindings/boost/numeric/bindings/blas/level1/axpy.hpp
==============================================================================
--- sandbox/numeric_bindings/boost/numeric/bindings/blas/level1/axpy.hpp (original)
+++ sandbox/numeric_bindings/boost/numeric/bindings/blas/level1/axpy.hpp 2012-08-07 09:40:31 EDT (Tue, 07 Aug 2012)
@@ -137,7 +137,7 @@
inline void axpy( const int n, const std::complex<double> a,
const std::complex<double>* x, const int incx,
std::complex<double>* y, const int incy ) {
- // NOT FOUND();
+ cublasZaxpy( n, a, x, incx, y, incy );
}
#else
Modified: sandbox/numeric_bindings/boost/numeric/bindings/blas/level1/copy.hpp
==============================================================================
--- sandbox/numeric_bindings/boost/numeric/bindings/blas/level1/copy.hpp (original)
+++ sandbox/numeric_bindings/boost/numeric/bindings/blas/level1/copy.hpp 2012-08-07 09:40:31 EDT (Tue, 07 Aug 2012)
@@ -133,7 +133,7 @@
//
inline void copy( const int n, const std::complex<double>* x, const int incx,
std::complex<double>* y, const int incy ) {
- // NOT FOUND();
+ cublasZcopy( n, x, incx, y, incy );
}
#else
Modified: sandbox/numeric_bindings/boost/numeric/bindings/blas/level1/dotc.hpp
==============================================================================
--- sandbox/numeric_bindings/boost/numeric/bindings/blas/level1/dotc.hpp (original)
+++ sandbox/numeric_bindings/boost/numeric/bindings/blas/level1/dotc.hpp 2012-08-07 09:40:31 EDT (Tue, 07 Aug 2012)
@@ -97,7 +97,7 @@
//
inline std::complex<double> dotc( const int n, const std::complex<double>* x,
const int incx, const std::complex<double>* y, const int incy ) {
- return // NOT FOUND();
+ return cublasZdotc( n, x, incx, y, incy );
}
#else
Modified: sandbox/numeric_bindings/boost/numeric/bindings/blas/level1/iamax.hpp
==============================================================================
--- sandbox/numeric_bindings/boost/numeric/bindings/blas/level1/iamax.hpp (original)
+++ sandbox/numeric_bindings/boost/numeric/bindings/blas/level1/iamax.hpp 2012-08-07 09:40:31 EDT (Tue, 07 Aug 2012)
@@ -131,7 +131,7 @@
//
inline std::ptrdiff_t iamax( const int n, const std::complex<double>* x,
const int incx ) {
- return // NOT FOUND();
+ return cublasIzamax( n, x, incx );
}
#else
Modified: sandbox/numeric_bindings/boost/numeric/bindings/blas/level1/nrm2.hpp
==============================================================================
--- sandbox/numeric_bindings/boost/numeric/bindings/blas/level1/nrm2.hpp (original)
+++ sandbox/numeric_bindings/boost/numeric/bindings/blas/level1/nrm2.hpp 2012-08-07 09:40:31 EDT (Tue, 07 Aug 2012)
@@ -129,7 +129,7 @@
//
inline double nrm2( const int n, const std::complex<double>* x,
const int incx ) {
- return // NOT FOUND();
+ return cublasDznrm2( n, x, incx );
}
#else
Modified: sandbox/numeric_bindings/boost/numeric/bindings/blas/level1/rot.hpp
==============================================================================
--- sandbox/numeric_bindings/boost/numeric/bindings/blas/level1/rot.hpp (original)
+++ sandbox/numeric_bindings/boost/numeric/bindings/blas/level1/rot.hpp 2012-08-07 09:40:31 EDT (Tue, 07 Aug 2012)
@@ -137,7 +137,7 @@
inline void rot( const int n, std::complex<double>* x, const int incx,
std::complex<double>* y, const int incy, const double c,
const double s ) {
- // NOT FOUND();
+ cublasZdrot( n, x, incx, y, incy, c, s );
}
#else
Modified: sandbox/numeric_bindings/boost/numeric/bindings/blas/level1/rotg.hpp
==============================================================================
--- sandbox/numeric_bindings/boost/numeric/bindings/blas/level1/rotg.hpp (original)
+++ sandbox/numeric_bindings/boost/numeric/bindings/blas/level1/rotg.hpp 2012-08-07 09:40:31 EDT (Tue, 07 Aug 2012)
@@ -119,7 +119,7 @@
//
inline void rotg( std::complex<float>& a, std::complex<float>& b, float& c,
std::complex<float>& s ) {
- // NOT FOUND();
+ cublasCrotg( &a, b, &c, &s );
}
//
@@ -129,7 +129,7 @@
//
inline void rotg( std::complex<double>& a, std::complex<double>& b, double& c,
std::complex<double>& s ) {
- // NOT FOUND();
+ cublasZrotg( &a, b, &c, &s );
}
#else
Modified: sandbox/numeric_bindings/boost/numeric/bindings/blas/level1/scal.hpp
==============================================================================
--- sandbox/numeric_bindings/boost/numeric/bindings/blas/level1/scal.hpp (original)
+++ sandbox/numeric_bindings/boost/numeric/bindings/blas/level1/scal.hpp 2012-08-07 09:40:31 EDT (Tue, 07 Aug 2012)
@@ -149,7 +149,7 @@
//
inline void scal( const int n, const double a, std::complex<double>* x,
const int incx ) {
- // NOT FOUND();
+ cublasZdscal( n, a, x, incx );
}
//
Modified: sandbox/numeric_bindings/boost/numeric/bindings/blas/level1/swap.hpp
==============================================================================
--- sandbox/numeric_bindings/boost/numeric/bindings/blas/level1/swap.hpp (original)
+++ sandbox/numeric_bindings/boost/numeric/bindings/blas/level1/swap.hpp 2012-08-07 09:40:31 EDT (Tue, 07 Aug 2012)
@@ -133,7 +133,7 @@
//
inline void swap( const int n, std::complex<double>* x, const int incx,
std::complex<double>* y, const int incy ) {
- // NOT FOUND();
+ cublasZswap( n, x, incx, y, incy );
}
#else
Modified: sandbox/numeric_bindings/boost/numeric/bindings/blas/level2/gbmv.hpp
==============================================================================
--- sandbox/numeric_bindings/boost/numeric/bindings/blas/level2/gbmv.hpp (original)
+++ sandbox/numeric_bindings/boost/numeric/bindings/blas/level2/gbmv.hpp 2012-08-07 09:40:31 EDT (Tue, 07 Aug 2012)
@@ -64,10 +64,10 @@
// * float value-type.
//
template< typename Order, typename Trans >
-inline void gbmv( const Order order, const Trans trans, const int m,
- const int n, const int kl, const int ku, const float alpha,
- const float* a, const int lda, const float* x, const int incx,
- const float beta, float* y, const int incy ) {
+inline void gbmv( const Order, const Trans, const int m, const int n,
+ const int kl, const int ku, const float alpha, const float* a,
+ const int lda, const float* x, const int incx, const float beta,
+ float* y, const int incy ) {
cblas_sgbmv( cblas_option< Order >::value, cblas_option< Trans >::value,
m, n, kl, ku, alpha, a, lda, x, incx, beta, y, incy );
}
@@ -78,10 +78,10 @@
// * double value-type.
//
template< typename Order, typename Trans >
-inline void gbmv( const Order order, const Trans trans, const int m,
- const int n, const int kl, const int ku, const double alpha,
- const double* a, const int lda, const double* x, const int incx,
- const double beta, double* y, const int incy ) {
+inline void gbmv( const Order, const Trans, const int m, const int n,
+ const int kl, const int ku, const double alpha, const double* a,
+ const int lda, const double* x, const int incx, const double beta,
+ double* y, const int incy ) {
cblas_dgbmv( cblas_option< Order >::value, cblas_option< Trans >::value,
m, n, kl, ku, alpha, a, lda, x, incx, beta, y, incy );
}
@@ -92,10 +92,10 @@
// * complex<float> value-type.
//
template< typename Order, typename Trans >
-inline void gbmv( const Order order, const Trans trans, const int m,
- const int n, const int kl, const int ku,
- const std::complex<float> alpha, const std::complex<float>* a,
- const int lda, const std::complex<float>* x, const int incx,
+inline void gbmv( const Order, const Trans, const int m, const int n,
+ const int kl, const int ku, const std::complex<float> alpha,
+ const std::complex<float>* a, const int lda,
+ const std::complex<float>* x, const int incx,
const std::complex<float> beta, std::complex<float>* y,
const int incy ) {
cblas_cgbmv( cblas_option< Order >::value, cblas_option< Trans >::value,
@@ -108,10 +108,10 @@
// * complex<double> value-type.
//
template< typename Order, typename Trans >
-inline void gbmv( const Order order, const Trans trans, const int m,
- const int n, const int kl, const int ku,
- const std::complex<double> alpha, const std::complex<double>* a,
- const int lda, const std::complex<double>* x, const int incx,
+inline void gbmv( const Order, const Trans, const int m, const int n,
+ const int kl, const int ku, const std::complex<double> alpha,
+ const std::complex<double>* a, const int lda,
+ const std::complex<double>* x, const int incx,
const std::complex<double> beta, std::complex<double>* y,
const int incy ) {
cblas_zgbmv( cblas_option< Order >::value, cblas_option< Trans >::value,
@@ -125,10 +125,10 @@
// * float value-type.
//
template< typename Order, typename Trans >
-inline void gbmv( const Order order, const Trans trans, const int m,
- const int n, const int kl, const int ku, const float alpha,
- const float* a, const int lda, const float* x, const int incx,
- const float beta, float* y, const int incy ) {
+inline void gbmv( const Order, const Trans, const int m, const int n,
+ const int kl, const int ku, const float alpha, const float* a,
+ const int lda, const float* x, const int incx, const float beta,
+ float* y, const int incy ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
cublasSgbmv( blas_option< Trans >::value, m, n, kl, ku, alpha, a, lda, x,
incx, beta, y, incy );
@@ -140,12 +140,13 @@
// * double value-type.
//
template< typename Order, typename Trans >
-inline void gbmv( const Order order, const Trans trans, const int m,
- const int n, const int kl, const int ku, const double alpha,
- const double* a, const int lda, const double* x, const int incx,
- const double beta, double* y, const int incy ) {
+inline void gbmv( const Order, const Trans, const int m, const int n,
+ const int kl, const int ku, const double alpha, const double* a,
+ const int lda, const double* x, const int incx, const double beta,
+ double* y, const int incy ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
- // NOT FOUND();
+ cublasDgbmv( blas_option< Trans >::value, m, n, kl, ku, alpha, a, lda, x,
+ incx, beta, y, incy );
}
//
@@ -154,10 +155,10 @@
// * complex<float> value-type.
//
template< typename Order, typename Trans >
-inline void gbmv( const Order order, const Trans trans, const int m,
- const int n, const int kl, const int ku,
- const std::complex<float> alpha, const std::complex<float>* a,
- const int lda, const std::complex<float>* x, const int incx,
+inline void gbmv( const Order, const Trans, const int m, const int n,
+ const int kl, const int ku, const std::complex<float> alpha,
+ const std::complex<float>* a, const int lda,
+ const std::complex<float>* x, const int incx,
const std::complex<float> beta, std::complex<float>* y,
const int incy ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
@@ -171,14 +172,15 @@
// * complex<double> value-type.
//
template< typename Order, typename Trans >
-inline void gbmv( const Order order, const Trans trans, const int m,
- const int n, const int kl, const int ku,
- const std::complex<double> alpha, const std::complex<double>* a,
- const int lda, const std::complex<double>* x, const int incx,
+inline void gbmv( const Order, const Trans, const int m, const int n,
+ const int kl, const int ku, const std::complex<double> alpha,
+ const std::complex<double>* a, const int lda,
+ const std::complex<double>* x, const int incx,
const std::complex<double> beta, std::complex<double>* y,
const int incy ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
- // NOT FOUND();
+ cublasZgbmv( blas_option< Trans >::value, m, n, kl, ku, alpha, a, lda, x,
+ incx, beta, y, incy );
}
#else
@@ -188,7 +190,7 @@
// * float value-type.
//
template< typename Order, typename Trans >
-inline void gbmv( const Order order, const Trans trans, const fortran_int_t m,
+inline void gbmv( const Order, const Trans, const fortran_int_t m,
const fortran_int_t n, const fortran_int_t kl, const fortran_int_t ku,
const float alpha, const float* a, const fortran_int_t lda,
const float* x, const fortran_int_t incx, const float beta, float* y,
@@ -204,7 +206,7 @@
// * double value-type.
//
template< typename Order, typename Trans >
-inline void gbmv( const Order order, const Trans trans, const fortran_int_t m,
+inline void gbmv( const Order, const Trans, const fortran_int_t m,
const fortran_int_t n, const fortran_int_t kl, const fortran_int_t ku,
const double alpha, const double* a, const fortran_int_t lda,
const double* x, const fortran_int_t incx, const double beta,
@@ -220,7 +222,7 @@
// * complex<float> value-type.
//
template< typename Order, typename Trans >
-inline void gbmv( const Order order, const Trans trans, const fortran_int_t m,
+inline void gbmv( const Order, const Trans, const fortran_int_t m,
const fortran_int_t n, const fortran_int_t kl, const fortran_int_t ku,
const std::complex<float> alpha, const std::complex<float>* a,
const fortran_int_t lda, const std::complex<float>* x,
@@ -237,7 +239,7 @@
// * complex<double> value-type.
//
template< typename Order, typename Trans >
-inline void gbmv( const Order order, const Trans trans, const fortran_int_t m,
+inline void gbmv( const Order, const Trans, const fortran_int_t m,
const fortran_int_t n, const fortran_int_t kl, const fortran_int_t ku,
const std::complex<double> alpha, const std::complex<double>* a,
const fortran_int_t lda, const std::complex<double>* x,
Modified: sandbox/numeric_bindings/boost/numeric/bindings/blas/level2/gemv.hpp
==============================================================================
--- sandbox/numeric_bindings/boost/numeric/bindings/blas/level2/gemv.hpp (original)
+++ sandbox/numeric_bindings/boost/numeric/bindings/blas/level2/gemv.hpp 2012-08-07 09:40:31 EDT (Tue, 07 Aug 2012)
@@ -63,10 +63,9 @@
// * float value-type.
//
template< typename Order, typename Trans >
-inline void gemv( const Order order, const Trans trans, const int m,
- const int n, const float alpha, const float* a, const int lda,
- const float* x, const int incx, const float beta, float* y,
- const int incy ) {
+inline void gemv( const Order, const Trans, const int m, const int n,
+ const float alpha, const float* a, const int lda, const float* x,
+ const int incx, const float beta, float* y, const int incy ) {
cblas_sgemv( cblas_option< Order >::value, cblas_option< Trans >::value,
m, n, alpha, a, lda, x, incx, beta, y, incy );
}
@@ -77,10 +76,9 @@
// * double value-type.
//
template< typename Order, typename Trans >
-inline void gemv( const Order order, const Trans trans, const int m,
- const int n, const double alpha, const double* a, const int lda,
- const double* x, const int incx, const double beta, double* y,
- const int incy ) {
+inline void gemv( const Order, const Trans, const int m, const int n,
+ const double alpha, const double* a, const int lda, const double* x,
+ const int incx, const double beta, double* y, const int incy ) {
cblas_dgemv( cblas_option< Order >::value, cblas_option< Trans >::value,
m, n, alpha, a, lda, x, incx, beta, y, incy );
}
@@ -91,10 +89,9 @@
// * complex<float> value-type.
//
template< typename Order, typename Trans >
-inline void gemv( const Order order, const Trans trans, const int m,
- const int n, const std::complex<float> alpha,
- const std::complex<float>* a, const int lda,
- const std::complex<float>* x, const int incx,
+inline void gemv( const Order, const Trans, const int m, const int n,
+ const std::complex<float> alpha, const std::complex<float>* a,
+ const int lda, const std::complex<float>* x, const int incx,
const std::complex<float> beta, std::complex<float>* y,
const int incy ) {
cblas_cgemv( cblas_option< Order >::value, cblas_option< Trans >::value,
@@ -107,10 +104,9 @@
// * complex<double> value-type.
//
template< typename Order, typename Trans >
-inline void gemv( const Order order, const Trans trans, const int m,
- const int n, const std::complex<double> alpha,
- const std::complex<double>* a, const int lda,
- const std::complex<double>* x, const int incx,
+inline void gemv( const Order, const Trans, const int m, const int n,
+ const std::complex<double> alpha, const std::complex<double>* a,
+ const int lda, const std::complex<double>* x, const int incx,
const std::complex<double> beta, std::complex<double>* y,
const int incy ) {
cblas_zgemv( cblas_option< Order >::value, cblas_option< Trans >::value,
@@ -124,10 +120,9 @@
// * float value-type.
//
template< typename Order, typename Trans >
-inline void gemv( const Order order, const Trans trans, const int m,
- const int n, const float alpha, const float* a, const int lda,
- const float* x, const int incx, const float beta, float* y,
- const int incy ) {
+inline void gemv( const Order, const Trans, const int m, const int n,
+ const float alpha, const float* a, const int lda, const float* x,
+ const int incx, const float beta, float* y, const int incy ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
cublasSgemv( blas_option< Trans >::value, m, n, alpha, a, lda, x, incx,
beta, y, incy );
@@ -139,10 +134,9 @@
// * double value-type.
//
template< typename Order, typename Trans >
-inline void gemv( const Order order, const Trans trans, const int m,
- const int n, const double alpha, const double* a, const int lda,
- const double* x, const int incx, const double beta, double* y,
- const int incy ) {
+inline void gemv( const Order, const Trans, const int m, const int n,
+ const double alpha, const double* a, const int lda, const double* x,
+ const int incx, const double beta, double* y, const int incy ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
cublasDgemv( blas_option< Trans >::value, m, n, alpha, a, lda, x, incx,
beta, y, incy );
@@ -154,10 +148,9 @@
// * complex<float> value-type.
//
template< typename Order, typename Trans >
-inline void gemv( const Order order, const Trans trans, const int m,
- const int n, const std::complex<float> alpha,
- const std::complex<float>* a, const int lda,
- const std::complex<float>* x, const int incx,
+inline void gemv( const Order, const Trans, const int m, const int n,
+ const std::complex<float> alpha, const std::complex<float>* a,
+ const int lda, const std::complex<float>* x, const int incx,
const std::complex<float> beta, std::complex<float>* y,
const int incy ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
@@ -171,10 +164,9 @@
// * complex<double> value-type.
//
template< typename Order, typename Trans >
-inline void gemv( const Order order, const Trans trans, const int m,
- const int n, const std::complex<double> alpha,
- const std::complex<double>* a, const int lda,
- const std::complex<double>* x, const int incx,
+inline void gemv( const Order, const Trans, const int m, const int n,
+ const std::complex<double> alpha, const std::complex<double>* a,
+ const int lda, const std::complex<double>* x, const int incx,
const std::complex<double> beta, std::complex<double>* y,
const int incy ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
@@ -189,7 +181,7 @@
// * float value-type.
//
template< typename Order, typename Trans >
-inline void gemv( const Order order, const Trans trans, const fortran_int_t m,
+inline void gemv( const Order, const Trans, const fortran_int_t m,
const fortran_int_t n, const float alpha, const float* a,
const fortran_int_t lda, const float* x, const fortran_int_t incx,
const float beta, float* y, const fortran_int_t incy ) {
@@ -204,7 +196,7 @@
// * double value-type.
//
template< typename Order, typename Trans >
-inline void gemv( const Order order, const Trans trans, const fortran_int_t m,
+inline void gemv( const Order, const Trans, const fortran_int_t m,
const fortran_int_t n, const double alpha, const double* a,
const fortran_int_t lda, const double* x, const fortran_int_t incx,
const double beta, double* y, const fortran_int_t incy ) {
@@ -219,7 +211,7 @@
// * complex<float> value-type.
//
template< typename Order, typename Trans >
-inline void gemv( const Order order, const Trans trans, const fortran_int_t m,
+inline void gemv( const Order, const Trans, const fortran_int_t m,
const fortran_int_t n, const std::complex<float> alpha,
const std::complex<float>* a, const fortran_int_t lda,
const std::complex<float>* x, const fortran_int_t incx,
@@ -236,7 +228,7 @@
// * complex<double> value-type.
//
template< typename Order, typename Trans >
-inline void gemv( const Order order, const Trans trans, const fortran_int_t m,
+inline void gemv( const Order, const Trans, const fortran_int_t m,
const fortran_int_t n, const std::complex<double> alpha,
const std::complex<double>* a, const fortran_int_t lda,
const std::complex<double>* x, const fortran_int_t incx,
Modified: sandbox/numeric_bindings/boost/numeric/bindings/blas/level2/ger.hpp
==============================================================================
--- sandbox/numeric_bindings/boost/numeric/bindings/blas/level2/ger.hpp (original)
+++ sandbox/numeric_bindings/boost/numeric/bindings/blas/level2/ger.hpp 2012-08-07 09:40:31 EDT (Tue, 07 Aug 2012)
@@ -62,9 +62,9 @@
// * float value-type.
//
template< typename Order >
-inline void ger( const Order order, const int m, const int n,
- const float alpha, const float* x, const int incx, const float* y,
- const int incy, float* a, const int lda ) {
+inline void ger( const Order, const int m, const int n, const float alpha,
+ const float* x, const int incx, const float* y, const int incy,
+ float* a, const int lda ) {
cblas_sger( cblas_option< Order >::value, m, n, alpha, x, incx, y, incy,
a, lda );
}
@@ -75,9 +75,9 @@
// * double value-type.
//
template< typename Order >
-inline void ger( const Order order, const int m, const int n,
- const double alpha, const double* x, const int incx, const double* y,
- const int incy, double* a, const int lda ) {
+inline void ger( const Order, const int m, const int n, const double alpha,
+ const double* x, const int incx, const double* y, const int incy,
+ double* a, const int lda ) {
cblas_dger( cblas_option< Order >::value, m, n, alpha, x, incx, y, incy,
a, lda );
}
@@ -89,9 +89,9 @@
// * float value-type.
//
template< typename Order >
-inline void ger( const Order order, const int m, const int n,
- const float alpha, const float* x, const int incx, const float* y,
- const int incy, float* a, const int lda ) {
+inline void ger( const Order, const int m, const int n, const float alpha,
+ const float* x, const int incx, const float* y, const int incy,
+ float* a, const int lda ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
cublasSger( m, n, alpha, x, incx, y, incy, a, lda );
}
@@ -102,9 +102,9 @@
// * double value-type.
//
template< typename Order >
-inline void ger( const Order order, const int m, const int n,
- const double alpha, const double* x, const int incx, const double* y,
- const int incy, double* a, const int lda ) {
+inline void ger( const Order, const int m, const int n, const double alpha,
+ const double* x, const int incx, const double* y, const int incy,
+ double* a, const int lda ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
cublasDger( m, n, alpha, x, incx, y, incy, a, lda );
}
@@ -116,10 +116,10 @@
// * float value-type.
//
template< typename Order >
-inline void ger( const Order order, const fortran_int_t m,
- const fortran_int_t n, const float alpha, const float* x,
- const fortran_int_t incx, const float* y, const fortran_int_t incy,
- float* a, const fortran_int_t lda ) {
+inline void ger( const Order, const fortran_int_t m, const fortran_int_t n,
+ const float alpha, const float* x, const fortran_int_t incx,
+ const float* y, const fortran_int_t incy, float* a,
+ const fortran_int_t lda ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
BLAS_SGER( &m, &n, &alpha, x, &incx, y, &incy, a, &lda );
}
@@ -130,10 +130,10 @@
// * double value-type.
//
template< typename Order >
-inline void ger( const Order order, const fortran_int_t m,
- const fortran_int_t n, const double alpha, const double* x,
- const fortran_int_t incx, const double* y, const fortran_int_t incy,
- double* a, const fortran_int_t lda ) {
+inline void ger( const Order, const fortran_int_t m, const fortran_int_t n,
+ const double alpha, const double* x, const fortran_int_t incx,
+ const double* y, const fortran_int_t incy, double* a,
+ const fortran_int_t lda ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
BLAS_DGER( &m, &n, &alpha, x, &incx, y, &incy, a, &lda );
}
Modified: sandbox/numeric_bindings/boost/numeric/bindings/blas/level2/gerc.hpp
==============================================================================
--- sandbox/numeric_bindings/boost/numeric/bindings/blas/level2/gerc.hpp (original)
+++ sandbox/numeric_bindings/boost/numeric/bindings/blas/level2/gerc.hpp 2012-08-07 09:40:31 EDT (Tue, 07 Aug 2012)
@@ -62,7 +62,7 @@
// * complex<float> value-type.
//
template< typename Order >
-inline void gerc( const Order order, const int m, const int n,
+inline void gerc( const Order, const int m, const int n,
const std::complex<float> alpha, const std::complex<float>* x,
const int incx, const std::complex<float>* y, const int incy,
std::complex<float>* a, const int lda ) {
@@ -76,7 +76,7 @@
// * complex<double> value-type.
//
template< typename Order >
-inline void gerc( const Order order, const int m, const int n,
+inline void gerc( const Order, const int m, const int n,
const std::complex<double> alpha, const std::complex<double>* x,
const int incx, const std::complex<double>* y, const int incy,
std::complex<double>* a, const int lda ) {
@@ -91,7 +91,7 @@
// * complex<float> value-type.
//
template< typename Order >
-inline void gerc( const Order order, const int m, const int n,
+inline void gerc( const Order, const int m, const int n,
const std::complex<float> alpha, const std::complex<float>* x,
const int incx, const std::complex<float>* y, const int incy,
std::complex<float>* a, const int lda ) {
@@ -105,12 +105,12 @@
// * complex<double> value-type.
//
template< typename Order >
-inline void gerc( const Order order, const int m, const int n,
+inline void gerc( const Order, const int m, const int n,
const std::complex<double> alpha, const std::complex<double>* x,
const int incx, const std::complex<double>* y, const int incy,
std::complex<double>* a, const int lda ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
- // NOT FOUND();
+ cublasZgerc( m, n, alpha, x, incx, y, incy, a, lda );
}
#else
@@ -120,11 +120,11 @@
// * complex<float> value-type.
//
template< typename Order >
-inline void gerc( const Order order, const fortran_int_t m,
- const fortran_int_t n, const std::complex<float> alpha,
- const std::complex<float>* x, const fortran_int_t incx,
- const std::complex<float>* y, const fortran_int_t incy,
- std::complex<float>* a, const fortran_int_t lda ) {
+inline void gerc( const Order, const fortran_int_t m, const fortran_int_t n,
+ const std::complex<float> alpha, const std::complex<float>* x,
+ const fortran_int_t incx, const std::complex<float>* y,
+ const fortran_int_t incy, std::complex<float>* a,
+ const fortran_int_t lda ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
BLAS_CGERC( &m, &n, &alpha, x, &incx, y, &incy, a, &lda );
}
@@ -135,11 +135,11 @@
// * complex<double> value-type.
//
template< typename Order >
-inline void gerc( const Order order, const fortran_int_t m,
- const fortran_int_t n, const std::complex<double> alpha,
- const std::complex<double>* x, const fortran_int_t incx,
- const std::complex<double>* y, const fortran_int_t incy,
- std::complex<double>* a, const fortran_int_t lda ) {
+inline void gerc( const Order, const fortran_int_t m, const fortran_int_t n,
+ const std::complex<double> alpha, const std::complex<double>* x,
+ const fortran_int_t incx, const std::complex<double>* y,
+ const fortran_int_t incy, std::complex<double>* a,
+ const fortran_int_t lda ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
BLAS_ZGERC( &m, &n, &alpha, x, &incx, y, &incy, a, &lda );
}
Modified: sandbox/numeric_bindings/boost/numeric/bindings/blas/level2/geru.hpp
==============================================================================
--- sandbox/numeric_bindings/boost/numeric/bindings/blas/level2/geru.hpp (original)
+++ sandbox/numeric_bindings/boost/numeric/bindings/blas/level2/geru.hpp 2012-08-07 09:40:31 EDT (Tue, 07 Aug 2012)
@@ -62,7 +62,7 @@
// * complex<float> value-type.
//
template< typename Order >
-inline void geru( const Order order, const int m, const int n,
+inline void geru( const Order, const int m, const int n,
const std::complex<float> alpha, const std::complex<float>* x,
const int incx, const std::complex<float>* y, const int incy,
std::complex<float>* a, const int lda ) {
@@ -76,7 +76,7 @@
// * complex<double> value-type.
//
template< typename Order >
-inline void geru( const Order order, const int m, const int n,
+inline void geru( const Order, const int m, const int n,
const std::complex<double> alpha, const std::complex<double>* x,
const int incx, const std::complex<double>* y, const int incy,
std::complex<double>* a, const int lda ) {
@@ -91,7 +91,7 @@
// * complex<float> value-type.
//
template< typename Order >
-inline void geru( const Order order, const int m, const int n,
+inline void geru( const Order, const int m, const int n,
const std::complex<float> alpha, const std::complex<float>* x,
const int incx, const std::complex<float>* y, const int incy,
std::complex<float>* a, const int lda ) {
@@ -105,12 +105,12 @@
// * complex<double> value-type.
//
template< typename Order >
-inline void geru( const Order order, const int m, const int n,
+inline void geru( const Order, const int m, const int n,
const std::complex<double> alpha, const std::complex<double>* x,
const int incx, const std::complex<double>* y, const int incy,
std::complex<double>* a, const int lda ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
- // NOT FOUND();
+ cublasZgeru( m, n, alpha, x, incx, y, incy, a, lda );
}
#else
@@ -120,11 +120,11 @@
// * complex<float> value-type.
//
template< typename Order >
-inline void geru( const Order order, const fortran_int_t m,
- const fortran_int_t n, const std::complex<float> alpha,
- const std::complex<float>* x, const fortran_int_t incx,
- const std::complex<float>* y, const fortran_int_t incy,
- std::complex<float>* a, const fortran_int_t lda ) {
+inline void geru( const Order, const fortran_int_t m, const fortran_int_t n,
+ const std::complex<float> alpha, const std::complex<float>* x,
+ const fortran_int_t incx, const std::complex<float>* y,
+ const fortran_int_t incy, std::complex<float>* a,
+ const fortran_int_t lda ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
BLAS_CGERU( &m, &n, &alpha, x, &incx, y, &incy, a, &lda );
}
@@ -135,11 +135,11 @@
// * complex<double> value-type.
//
template< typename Order >
-inline void geru( const Order order, const fortran_int_t m,
- const fortran_int_t n, const std::complex<double> alpha,
- const std::complex<double>* x, const fortran_int_t incx,
- const std::complex<double>* y, const fortran_int_t incy,
- std::complex<double>* a, const fortran_int_t lda ) {
+inline void geru( const Order, const fortran_int_t m, const fortran_int_t n,
+ const std::complex<double> alpha, const std::complex<double>* x,
+ const fortran_int_t incx, const std::complex<double>* y,
+ const fortran_int_t incy, std::complex<double>* a,
+ const fortran_int_t lda ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
BLAS_ZGERU( &m, &n, &alpha, x, &incx, y, &incy, a, &lda );
}
Modified: sandbox/numeric_bindings/boost/numeric/bindings/blas/level2/hbmv.hpp
==============================================================================
--- sandbox/numeric_bindings/boost/numeric/bindings/blas/level2/hbmv.hpp (original)
+++ sandbox/numeric_bindings/boost/numeric/bindings/blas/level2/hbmv.hpp 2012-08-07 09:40:31 EDT (Tue, 07 Aug 2012)
@@ -64,10 +64,9 @@
// * float value-type.
//
template< typename Order, typename UpLo >
-inline void hbmv( const Order order, const UpLo uplo, const int n,
- const int k, const float alpha, const float* a, const int lda,
- const float* x, const int incx, const float beta, float* y,
- const int incy ) {
+inline void hbmv( const Order, const UpLo, const int n, const int k,
+ const float alpha, const float* a, const int lda, const float* x,
+ const int incx, const float beta, float* y, const int incy ) {
cblas_ssbmv( cblas_option< Order >::value, cblas_option< UpLo >::value, n,
k, alpha, a, lda, x, incx, beta, y, incy );
}
@@ -78,10 +77,9 @@
// * double value-type.
//
template< typename Order, typename UpLo >
-inline void hbmv( const Order order, const UpLo uplo, const int n,
- const int k, const double alpha, const double* a, const int lda,
- const double* x, const int incx, const double beta, double* y,
- const int incy ) {
+inline void hbmv( const Order, const UpLo, const int n, const int k,
+ const double alpha, const double* a, const int lda, const double* x,
+ const int incx, const double beta, double* y, const int incy ) {
cblas_dsbmv( cblas_option< Order >::value, cblas_option< UpLo >::value, n,
k, alpha, a, lda, x, incx, beta, y, incy );
}
@@ -92,10 +90,9 @@
// * complex<float> value-type.
//
template< typename Order, typename UpLo >
-inline void hbmv( const Order order, const UpLo uplo, const int n,
- const int k, const std::complex<float> alpha,
- const std::complex<float>* a, const int lda,
- const std::complex<float>* x, const int incx,
+inline void hbmv( const Order, const UpLo, const int n, const int k,
+ const std::complex<float> alpha, const std::complex<float>* a,
+ const int lda, const std::complex<float>* x, const int incx,
const std::complex<float> beta, std::complex<float>* y,
const int incy ) {
cblas_chbmv( cblas_option< Order >::value, cblas_option< UpLo >::value, n,
@@ -108,10 +105,9 @@
// * complex<double> value-type.
//
template< typename Order, typename UpLo >
-inline void hbmv( const Order order, const UpLo uplo, const int n,
- const int k, const std::complex<double> alpha,
- const std::complex<double>* a, const int lda,
- const std::complex<double>* x, const int incx,
+inline void hbmv( const Order, const UpLo, const int n, const int k,
+ const std::complex<double> alpha, const std::complex<double>* a,
+ const int lda, const std::complex<double>* x, const int incx,
const std::complex<double> beta, std::complex<double>* y,
const int incy ) {
cblas_zhbmv( cblas_option< Order >::value, cblas_option< UpLo >::value, n,
@@ -125,10 +121,9 @@
// * float value-type.
//
template< typename Order, typename UpLo >
-inline void hbmv( const Order order, const UpLo uplo, const int n,
- const int k, const float alpha, const float* a, const int lda,
- const float* x, const int incx, const float beta, float* y,
- const int incy ) {
+inline void hbmv( const Order, const UpLo, const int n, const int k,
+ const float alpha, const float* a, const int lda, const float* x,
+ const int incx, const float beta, float* y, const int incy ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
cublasSsbmv( blas_option< UpLo >::value, n, k, alpha, a, lda, x, incx,
beta, y, incy );
@@ -140,12 +135,12 @@
// * double value-type.
//
template< typename Order, typename UpLo >
-inline void hbmv( const Order order, const UpLo uplo, const int n,
- const int k, const double alpha, const double* a, const int lda,
- const double* x, const int incx, const double beta, double* y,
- const int incy ) {
+inline void hbmv( const Order, const UpLo, const int n, const int k,
+ const double alpha, const double* a, const int lda, const double* x,
+ const int incx, const double beta, double* y, const int incy ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
- // NOT FOUND();
+ cublasDsbmv( blas_option< UpLo >::value, n, k, alpha, a, lda, x, incx,
+ beta, y, incy );
}
//
@@ -154,10 +149,9 @@
// * complex<float> value-type.
//
template< typename Order, typename UpLo >
-inline void hbmv( const Order order, const UpLo uplo, const int n,
- const int k, const std::complex<float> alpha,
- const std::complex<float>* a, const int lda,
- const std::complex<float>* x, const int incx,
+inline void hbmv( const Order, const UpLo, const int n, const int k,
+ const std::complex<float> alpha, const std::complex<float>* a,
+ const int lda, const std::complex<float>* x, const int incx,
const std::complex<float> beta, std::complex<float>* y,
const int incy ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
@@ -171,14 +165,14 @@
// * complex<double> value-type.
//
template< typename Order, typename UpLo >
-inline void hbmv( const Order order, const UpLo uplo, const int n,
- const int k, const std::complex<double> alpha,
- const std::complex<double>* a, const int lda,
- const std::complex<double>* x, const int incx,
+inline void hbmv( const Order, const UpLo, const int n, const int k,
+ const std::complex<double> alpha, const std::complex<double>* a,
+ const int lda, const std::complex<double>* x, const int incx,
const std::complex<double> beta, std::complex<double>* y,
const int incy ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
- // NOT FOUND();
+ cublasZhbmv( blas_option< UpLo >::value, n, k, alpha, a, lda, x, incx,
+ beta, y, incy );
}
#else
@@ -188,7 +182,7 @@
// * float value-type.
//
template< typename Order, typename UpLo >
-inline void hbmv( const Order order, const UpLo uplo, const fortran_int_t n,
+inline void hbmv( const Order, const UpLo, const fortran_int_t n,
const fortran_int_t k, const float alpha, const float* a,
const fortran_int_t lda, const float* x, const fortran_int_t incx,
const float beta, float* y, const fortran_int_t incy ) {
@@ -203,7 +197,7 @@
// * double value-type.
//
template< typename Order, typename UpLo >
-inline void hbmv( const Order order, const UpLo uplo, const fortran_int_t n,
+inline void hbmv( const Order, const UpLo, const fortran_int_t n,
const fortran_int_t k, const double alpha, const double* a,
const fortran_int_t lda, const double* x, const fortran_int_t incx,
const double beta, double* y, const fortran_int_t incy ) {
@@ -218,7 +212,7 @@
// * complex<float> value-type.
//
template< typename Order, typename UpLo >
-inline void hbmv( const Order order, const UpLo uplo, const fortran_int_t n,
+inline void hbmv( const Order, const UpLo, const fortran_int_t n,
const fortran_int_t k, const std::complex<float> alpha,
const std::complex<float>* a, const fortran_int_t lda,
const std::complex<float>* x, const fortran_int_t incx,
@@ -235,7 +229,7 @@
// * complex<double> value-type.
//
template< typename Order, typename UpLo >
-inline void hbmv( const Order order, const UpLo uplo, const fortran_int_t n,
+inline void hbmv( const Order, const UpLo, const fortran_int_t n,
const fortran_int_t k, const std::complex<double> alpha,
const std::complex<double>* a, const fortran_int_t lda,
const std::complex<double>* x, const fortran_int_t incx,
Modified: sandbox/numeric_bindings/boost/numeric/bindings/blas/level2/hemv.hpp
==============================================================================
--- sandbox/numeric_bindings/boost/numeric/bindings/blas/level2/hemv.hpp (original)
+++ sandbox/numeric_bindings/boost/numeric/bindings/blas/level2/hemv.hpp 2012-08-07 09:40:31 EDT (Tue, 07 Aug 2012)
@@ -63,9 +63,9 @@
// * float value-type.
//
template< typename Order, typename UpLo >
-inline void hemv( const Order order, const UpLo uplo, const int n,
- const float alpha, const float* a, const int lda, const float* x,
- const int incx, const float beta, float* y, const int incy ) {
+inline void hemv( const Order, const UpLo, const int n, const float alpha,
+ const float* a, const int lda, const float* x, const int incx,
+ const float beta, float* y, const int incy ) {
cblas_ssymv( cblas_option< Order >::value, cblas_option< UpLo >::value, n,
alpha, a, lda, x, incx, beta, y, incy );
}
@@ -76,9 +76,9 @@
// * double value-type.
//
template< typename Order, typename UpLo >
-inline void hemv( const Order order, const UpLo uplo, const int n,
- const double alpha, const double* a, const int lda, const double* x,
- const int incx, const double beta, double* y, const int incy ) {
+inline void hemv( const Order, const UpLo, const int n, const double alpha,
+ const double* a, const int lda, const double* x, const int incx,
+ const double beta, double* y, const int incy ) {
cblas_dsymv( cblas_option< Order >::value, cblas_option< UpLo >::value, n,
alpha, a, lda, x, incx, beta, y, incy );
}
@@ -89,7 +89,7 @@
// * complex<float> value-type.
//
template< typename Order, typename UpLo >
-inline void hemv( const Order order, const UpLo uplo, const int n,
+inline void hemv( const Order, const UpLo, const int n,
const std::complex<float> alpha, const std::complex<float>* a,
const int lda, const std::complex<float>* x, const int incx,
const std::complex<float> beta, std::complex<float>* y,
@@ -104,7 +104,7 @@
// * complex<double> value-type.
//
template< typename Order, typename UpLo >
-inline void hemv( const Order order, const UpLo uplo, const int n,
+inline void hemv( const Order, const UpLo, const int n,
const std::complex<double> alpha, const std::complex<double>* a,
const int lda, const std::complex<double>* x, const int incx,
const std::complex<double> beta, std::complex<double>* y,
@@ -120,9 +120,9 @@
// * float value-type.
//
template< typename Order, typename UpLo >
-inline void hemv( const Order order, const UpLo uplo, const int n,
- const float alpha, const float* a, const int lda, const float* x,
- const int incx, const float beta, float* y, const int incy ) {
+inline void hemv( const Order, const UpLo, const int n, const float alpha,
+ const float* a, const int lda, const float* x, const int incx,
+ const float beta, float* y, const int incy ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
cublasSsymv( blas_option< UpLo >::value, n, alpha, a, lda, x, incx, beta,
y, incy );
@@ -134,11 +134,12 @@
// * double value-type.
//
template< typename Order, typename UpLo >
-inline void hemv( const Order order, const UpLo uplo, const int n,
- const double alpha, const double* a, const int lda, const double* x,
- const int incx, const double beta, double* y, const int incy ) {
+inline void hemv( const Order, const UpLo, const int n, const double alpha,
+ const double* a, const int lda, const double* x, const int incx,
+ const double beta, double* y, const int incy ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
- // NOT FOUND();
+ cublasDsymv( blas_option< UpLo >::value, n, alpha, a, lda, x, incx, beta,
+ y, incy );
}
//
@@ -147,7 +148,7 @@
// * complex<float> value-type.
//
template< typename Order, typename UpLo >
-inline void hemv( const Order order, const UpLo uplo, const int n,
+inline void hemv( const Order, const UpLo, const int n,
const std::complex<float> alpha, const std::complex<float>* a,
const int lda, const std::complex<float>* x, const int incx,
const std::complex<float> beta, std::complex<float>* y,
@@ -163,13 +164,14 @@
// * complex<double> value-type.
//
template< typename Order, typename UpLo >
-inline void hemv( const Order order, const UpLo uplo, const int n,
+inline void hemv( const Order, const UpLo, const int n,
const std::complex<double> alpha, const std::complex<double>* a,
const int lda, const std::complex<double>* x, const int incx,
const std::complex<double> beta, std::complex<double>* y,
const int incy ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
- // NOT FOUND();
+ cublasZhemv( blas_option< UpLo >::value, n, alpha, a, lda, x, incx, beta,
+ y, incy );
}
#else
@@ -179,7 +181,7 @@
// * float value-type.
//
template< typename Order, typename UpLo >
-inline void hemv( const Order order, const UpLo uplo, const fortran_int_t n,
+inline void hemv( const Order, const UpLo, const fortran_int_t n,
const float alpha, const float* a, const fortran_int_t lda,
const float* x, const fortran_int_t incx, const float beta, float* y,
const fortran_int_t incy ) {
@@ -194,7 +196,7 @@
// * double value-type.
//
template< typename Order, typename UpLo >
-inline void hemv( const Order order, const UpLo uplo, const fortran_int_t n,
+inline void hemv( const Order, const UpLo, const fortran_int_t n,
const double alpha, const double* a, const fortran_int_t lda,
const double* x, const fortran_int_t incx, const double beta,
double* y, const fortran_int_t incy ) {
@@ -209,7 +211,7 @@
// * complex<float> value-type.
//
template< typename Order, typename UpLo >
-inline void hemv( const Order order, const UpLo uplo, const fortran_int_t n,
+inline void hemv( const Order, const UpLo, const fortran_int_t n,
const std::complex<float> alpha, const std::complex<float>* a,
const fortran_int_t lda, const std::complex<float>* x,
const fortran_int_t incx, const std::complex<float> beta,
@@ -225,7 +227,7 @@
// * complex<double> value-type.
//
template< typename Order, typename UpLo >
-inline void hemv( const Order order, const UpLo uplo, const fortran_int_t n,
+inline void hemv( const Order, const UpLo, const fortran_int_t n,
const std::complex<double> alpha, const std::complex<double>* a,
const fortran_int_t lda, const std::complex<double>* x,
const fortran_int_t incx, const std::complex<double> beta,
Modified: sandbox/numeric_bindings/boost/numeric/bindings/blas/level2/her.hpp
==============================================================================
--- sandbox/numeric_bindings/boost/numeric/bindings/blas/level2/her.hpp (original)
+++ sandbox/numeric_bindings/boost/numeric/bindings/blas/level2/her.hpp 2012-08-07 09:40:31 EDT (Tue, 07 Aug 2012)
@@ -63,9 +63,8 @@
// * float value-type.
//
template< typename Order, typename UpLo >
-inline void her( const Order order, const UpLo uplo, const int n,
- const float alpha, const float* x, const int incx, float* a,
- const int lda ) {
+inline void her( const Order, const UpLo, const int n, const float alpha,
+ const float* x, const int incx, float* a, const int lda ) {
cblas_ssyr( cblas_option< Order >::value, cblas_option< UpLo >::value, n,
alpha, x, incx, a, lda );
}
@@ -76,9 +75,8 @@
// * double value-type.
//
template< typename Order, typename UpLo >
-inline void her( const Order order, const UpLo uplo, const int n,
- const double alpha, const double* x, const int incx, double* a,
- const int lda ) {
+inline void her( const Order, const UpLo, const int n, const double alpha,
+ const double* x, const int incx, double* a, const int lda ) {
cblas_dsyr( cblas_option< Order >::value, cblas_option< UpLo >::value, n,
alpha, x, incx, a, lda );
}
@@ -89,9 +87,9 @@
// * complex<float> value-type.
//
template< typename Order, typename UpLo >
-inline void her( const Order order, const UpLo uplo, const int n,
- const float alpha, const std::complex<float>* x, const int incx,
- std::complex<float>* a, const int lda ) {
+inline void her( const Order, const UpLo, const int n, const float alpha,
+ const std::complex<float>* x, const int incx, std::complex<float>* a,
+ const int lda ) {
cblas_cher( cblas_option< Order >::value, cblas_option< UpLo >::value, n,
alpha, x, incx, a, lda );
}
@@ -102,8 +100,8 @@
// * complex<double> value-type.
//
template< typename Order, typename UpLo >
-inline void her( const Order order, const UpLo uplo, const int n,
- const double alpha, const std::complex<double>* x, const int incx,
+inline void her( const Order, const UpLo, const int n, const double alpha,
+ const std::complex<double>* x, const int incx,
std::complex<double>* a, const int lda ) {
cblas_zher( cblas_option< Order >::value, cblas_option< UpLo >::value, n,
alpha, x, incx, a, lda );
@@ -116,9 +114,8 @@
// * float value-type.
//
template< typename Order, typename UpLo >
-inline void her( const Order order, const UpLo uplo, const int n,
- const float alpha, const float* x, const int incx, float* a,
- const int lda ) {
+inline void her( const Order, const UpLo, const int n, const float alpha,
+ const float* x, const int incx, float* a, const int lda ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
cublasSsyr( blas_option< UpLo >::value, n, alpha, x, incx, a, lda );
}
@@ -129,9 +126,8 @@
// * double value-type.
//
template< typename Order, typename UpLo >
-inline void her( const Order order, const UpLo uplo, const int n,
- const double alpha, const double* x, const int incx, double* a,
- const int lda ) {
+inline void her( const Order, const UpLo, const int n, const double alpha,
+ const double* x, const int incx, double* a, const int lda ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
cublasDsyr( blas_option< UpLo >::value, n, alpha, x, incx, a, lda );
}
@@ -142,9 +138,9 @@
// * complex<float> value-type.
//
template< typename Order, typename UpLo >
-inline void her( const Order order, const UpLo uplo, const int n,
- const float alpha, const std::complex<float>* x, const int incx,
- std::complex<float>* a, const int lda ) {
+inline void her( const Order, const UpLo, const int n, const float alpha,
+ const std::complex<float>* x, const int incx, std::complex<float>* a,
+ const int lda ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
cublasCher( blas_option< UpLo >::value, n, alpha, x, incx, a, lda );
}
@@ -155,11 +151,11 @@
// * complex<double> value-type.
//
template< typename Order, typename UpLo >
-inline void her( const Order order, const UpLo uplo, const int n,
- const double alpha, const std::complex<double>* x, const int incx,
+inline void her( const Order, const UpLo, const int n, const double alpha,
+ const std::complex<double>* x, const int incx,
std::complex<double>* a, const int lda ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
- // NOT FOUND();
+ cublasZher( blas_option< UpLo >::value, n, alpha, x, incx, a, lda );
}
#else
@@ -169,7 +165,7 @@
// * float value-type.
//
template< typename Order, typename UpLo >
-inline void her( const Order order, const UpLo uplo, const fortran_int_t n,
+inline void her( const Order, const UpLo, const fortran_int_t n,
const float alpha, const float* x, const fortran_int_t incx, float* a,
const fortran_int_t lda ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
@@ -182,7 +178,7 @@
// * double value-type.
//
template< typename Order, typename UpLo >
-inline void her( const Order order, const UpLo uplo, const fortran_int_t n,
+inline void her( const Order, const UpLo, const fortran_int_t n,
const double alpha, const double* x, const fortran_int_t incx,
double* a, const fortran_int_t lda ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
@@ -195,7 +191,7 @@
// * complex<float> value-type.
//
template< typename Order, typename UpLo >
-inline void her( const Order order, const UpLo uplo, const fortran_int_t n,
+inline void her( const Order, const UpLo, const fortran_int_t n,
const float alpha, const std::complex<float>* x,
const fortran_int_t incx, std::complex<float>* a,
const fortran_int_t lda ) {
@@ -209,7 +205,7 @@
// * complex<double> value-type.
//
template< typename Order, typename UpLo >
-inline void her( const Order order, const UpLo uplo, const fortran_int_t n,
+inline void her( const Order, const UpLo, const fortran_int_t n,
const double alpha, const std::complex<double>* x,
const fortran_int_t incx, std::complex<double>* a,
const fortran_int_t lda ) {
Modified: sandbox/numeric_bindings/boost/numeric/bindings/blas/level2/her2.hpp
==============================================================================
--- sandbox/numeric_bindings/boost/numeric/bindings/blas/level2/her2.hpp (original)
+++ sandbox/numeric_bindings/boost/numeric/bindings/blas/level2/her2.hpp 2012-08-07 09:40:31 EDT (Tue, 07 Aug 2012)
@@ -63,9 +63,9 @@
// * float value-type.
//
template< typename Order, typename UpLo >
-inline void her2( const Order order, const UpLo uplo, const int n,
- const float alpha, const float* x, const int incx, const float* y,
- const int incy, float* a, const int lda ) {
+inline void her2( const Order, const UpLo, const int n, const float alpha,
+ const float* x, const int incx, const float* y, const int incy,
+ float* a, const int lda ) {
cblas_ssyr2( cblas_option< Order >::value, cblas_option< UpLo >::value, n,
alpha, x, incx, y, incy, a, lda );
}
@@ -76,9 +76,9 @@
// * double value-type.
//
template< typename Order, typename UpLo >
-inline void her2( const Order order, const UpLo uplo, const int n,
- const double alpha, const double* x, const int incx, const double* y,
- const int incy, double* a, const int lda ) {
+inline void her2( const Order, const UpLo, const int n, const double alpha,
+ const double* x, const int incx, const double* y, const int incy,
+ double* a, const int lda ) {
cblas_dsyr2( cblas_option< Order >::value, cblas_option< UpLo >::value, n,
alpha, x, incx, y, incy, a, lda );
}
@@ -89,7 +89,7 @@
// * complex<float> value-type.
//
template< typename Order, typename UpLo >
-inline void her2( const Order order, const UpLo uplo, const int n,
+inline void her2( const Order, const UpLo, const int n,
const std::complex<float> alpha, const std::complex<float>* x,
const int incx, const std::complex<float>* y, const int incy,
std::complex<float>* a, const int lda ) {
@@ -103,7 +103,7 @@
// * complex<double> value-type.
//
template< typename Order, typename UpLo >
-inline void her2( const Order order, const UpLo uplo, const int n,
+inline void her2( const Order, const UpLo, const int n,
const std::complex<double> alpha, const std::complex<double>* x,
const int incx, const std::complex<double>* y, const int incy,
std::complex<double>* a, const int lda ) {
@@ -118,9 +118,9 @@
// * float value-type.
//
template< typename Order, typename UpLo >
-inline void her2( const Order order, const UpLo uplo, const int n,
- const float alpha, const float* x, const int incx, const float* y,
- const int incy, float* a, const int lda ) {
+inline void her2( const Order, const UpLo, const int n, const float alpha,
+ const float* x, const int incx, const float* y, const int incy,
+ float* a, const int lda ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
cublasSsyr2( blas_option< UpLo >::value, n, alpha, x, incx, y, incy, a,
lda );
@@ -132,11 +132,12 @@
// * double value-type.
//
template< typename Order, typename UpLo >
-inline void her2( const Order order, const UpLo uplo, const int n,
- const double alpha, const double* x, const int incx, const double* y,
- const int incy, double* a, const int lda ) {
+inline void her2( const Order, const UpLo, const int n, const double alpha,
+ const double* x, const int incx, const double* y, const int incy,
+ double* a, const int lda ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
- // NOT FOUND();
+ cublasDsyr2( blas_option< UpLo >::value, n, alpha, x, incx, y, incy, a,
+ lda );
}
//
@@ -145,7 +146,7 @@
// * complex<float> value-type.
//
template< typename Order, typename UpLo >
-inline void her2( const Order order, const UpLo uplo, const int n,
+inline void her2( const Order, const UpLo, const int n,
const std::complex<float> alpha, const std::complex<float>* x,
const int incx, const std::complex<float>* y, const int incy,
std::complex<float>* a, const int lda ) {
@@ -160,12 +161,13 @@
// * complex<double> value-type.
//
template< typename Order, typename UpLo >
-inline void her2( const Order order, const UpLo uplo, const int n,
+inline void her2( const Order, const UpLo, const int n,
const std::complex<double> alpha, const std::complex<double>* x,
const int incx, const std::complex<double>* y, const int incy,
std::complex<double>* a, const int lda ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
- // NOT FOUND();
+ cublasZher2( blas_option< UpLo >::value, n, alpha, x, incx, y, incy, a,
+ lda );
}
#else
@@ -175,7 +177,7 @@
// * float value-type.
//
template< typename Order, typename UpLo >
-inline void her2( const Order order, const UpLo uplo, const fortran_int_t n,
+inline void her2( const Order, const UpLo, const fortran_int_t n,
const float alpha, const float* x, const fortran_int_t incx,
const float* y, const fortran_int_t incy, float* a,
const fortran_int_t lda ) {
@@ -190,7 +192,7 @@
// * double value-type.
//
template< typename Order, typename UpLo >
-inline void her2( const Order order, const UpLo uplo, const fortran_int_t n,
+inline void her2( const Order, const UpLo, const fortran_int_t n,
const double alpha, const double* x, const fortran_int_t incx,
const double* y, const fortran_int_t incy, double* a,
const fortran_int_t lda ) {
@@ -205,7 +207,7 @@
// * complex<float> value-type.
//
template< typename Order, typename UpLo >
-inline void her2( const Order order, const UpLo uplo, const fortran_int_t n,
+inline void her2( const Order, const UpLo, const fortran_int_t n,
const std::complex<float> alpha, const std::complex<float>* x,
const fortran_int_t incx, const std::complex<float>* y,
const fortran_int_t incy, std::complex<float>* a,
@@ -221,7 +223,7 @@
// * complex<double> value-type.
//
template< typename Order, typename UpLo >
-inline void her2( const Order order, const UpLo uplo, const fortran_int_t n,
+inline void her2( const Order, const UpLo, const fortran_int_t n,
const std::complex<double> alpha, const std::complex<double>* x,
const fortran_int_t incx, const std::complex<double>* y,
const fortran_int_t incy, std::complex<double>* a,
Modified: sandbox/numeric_bindings/boost/numeric/bindings/blas/level2/hpmv.hpp
==============================================================================
--- sandbox/numeric_bindings/boost/numeric/bindings/blas/level2/hpmv.hpp (original)
+++ sandbox/numeric_bindings/boost/numeric/bindings/blas/level2/hpmv.hpp 2012-08-07 09:40:31 EDT (Tue, 07 Aug 2012)
@@ -63,9 +63,9 @@
// * float value-type.
//
template< typename Order, typename UpLo >
-inline void hpmv( const Order order, const UpLo uplo, const int n,
- const float alpha, const float* ap, const float* x, const int incx,
- const float beta, float* y, const int incy ) {
+inline void hpmv( const Order, const UpLo, const int n, const float alpha,
+ const float* ap, const float* x, const int incx, const float beta,
+ float* y, const int incy ) {
cblas_sspmv( cblas_option< Order >::value, cblas_option< UpLo >::value, n,
alpha, ap, x, incx, beta, y, incy );
}
@@ -76,9 +76,9 @@
// * double value-type.
//
template< typename Order, typename UpLo >
-inline void hpmv( const Order order, const UpLo uplo, const int n,
- const double alpha, const double* ap, const double* x, const int incx,
- const double beta, double* y, const int incy ) {
+inline void hpmv( const Order, const UpLo, const int n, const double alpha,
+ const double* ap, const double* x, const int incx, const double beta,
+ double* y, const int incy ) {
cblas_dspmv( cblas_option< Order >::value, cblas_option< UpLo >::value, n,
alpha, ap, x, incx, beta, y, incy );
}
@@ -89,7 +89,7 @@
// * complex<float> value-type.
//
template< typename Order, typename UpLo >
-inline void hpmv( const Order order, const UpLo uplo, const int n,
+inline void hpmv( const Order, const UpLo, const int n,
const std::complex<float> alpha, const std::complex<float>* ap,
const std::complex<float>* x, const int incx,
const std::complex<float> beta, std::complex<float>* y,
@@ -104,7 +104,7 @@
// * complex<double> value-type.
//
template< typename Order, typename UpLo >
-inline void hpmv( const Order order, const UpLo uplo, const int n,
+inline void hpmv( const Order, const UpLo, const int n,
const std::complex<double> alpha, const std::complex<double>* ap,
const std::complex<double>* x, const int incx,
const std::complex<double> beta, std::complex<double>* y,
@@ -120,9 +120,9 @@
// * float value-type.
//
template< typename Order, typename UpLo >
-inline void hpmv( const Order order, const UpLo uplo, const int n,
- const float alpha, const float* ap, const float* x, const int incx,
- const float beta, float* y, const int incy ) {
+inline void hpmv( const Order, const UpLo, const int n, const float alpha,
+ const float* ap, const float* x, const int incx, const float beta,
+ float* y, const int incy ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
cublasSspmv( blas_option< UpLo >::value, n, alpha, ap, x, incx, beta, y,
incy );
@@ -134,11 +134,12 @@
// * double value-type.
//
template< typename Order, typename UpLo >
-inline void hpmv( const Order order, const UpLo uplo, const int n,
- const double alpha, const double* ap, const double* x, const int incx,
- const double beta, double* y, const int incy ) {
+inline void hpmv( const Order, const UpLo, const int n, const double alpha,
+ const double* ap, const double* x, const int incx, const double beta,
+ double* y, const int incy ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
- // NOT FOUND();
+ cublasDspmv( blas_option< UpLo >::value, n, alpha, ap, x, incx, beta, y,
+ incy );
}
//
@@ -147,7 +148,7 @@
// * complex<float> value-type.
//
template< typename Order, typename UpLo >
-inline void hpmv( const Order order, const UpLo uplo, const int n,
+inline void hpmv( const Order, const UpLo, const int n,
const std::complex<float> alpha, const std::complex<float>* ap,
const std::complex<float>* x, const int incx,
const std::complex<float> beta, std::complex<float>* y,
@@ -163,13 +164,14 @@
// * complex<double> value-type.
//
template< typename Order, typename UpLo >
-inline void hpmv( const Order order, const UpLo uplo, const int n,
+inline void hpmv( const Order, const UpLo, const int n,
const std::complex<double> alpha, const std::complex<double>* ap,
const std::complex<double>* x, const int incx,
const std::complex<double> beta, std::complex<double>* y,
const int incy ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
- // NOT FOUND();
+ cublasZhpmv( blas_option< UpLo >::value, n, alpha, ap, x, incx, beta, y,
+ incy );
}
#else
@@ -179,7 +181,7 @@
// * float value-type.
//
template< typename Order, typename UpLo >
-inline void hpmv( const Order order, const UpLo uplo, const fortran_int_t n,
+inline void hpmv( const Order, const UpLo, const fortran_int_t n,
const float alpha, const float* ap, const float* x,
const fortran_int_t incx, const float beta, float* y,
const fortran_int_t incy ) {
@@ -194,7 +196,7 @@
// * double value-type.
//
template< typename Order, typename UpLo >
-inline void hpmv( const Order order, const UpLo uplo, const fortran_int_t n,
+inline void hpmv( const Order, const UpLo, const fortran_int_t n,
const double alpha, const double* ap, const double* x,
const fortran_int_t incx, const double beta, double* y,
const fortran_int_t incy ) {
@@ -209,7 +211,7 @@
// * complex<float> value-type.
//
template< typename Order, typename UpLo >
-inline void hpmv( const Order order, const UpLo uplo, const fortran_int_t n,
+inline void hpmv( const Order, const UpLo, const fortran_int_t n,
const std::complex<float> alpha, const std::complex<float>* ap,
const std::complex<float>* x, const fortran_int_t incx,
const std::complex<float> beta, std::complex<float>* y,
@@ -225,7 +227,7 @@
// * complex<double> value-type.
//
template< typename Order, typename UpLo >
-inline void hpmv( const Order order, const UpLo uplo, const fortran_int_t n,
+inline void hpmv( const Order, const UpLo, const fortran_int_t n,
const std::complex<double> alpha, const std::complex<double>* ap,
const std::complex<double>* x, const fortran_int_t incx,
const std::complex<double> beta, std::complex<double>* y,
Modified: sandbox/numeric_bindings/boost/numeric/bindings/blas/level2/hpr.hpp
==============================================================================
--- sandbox/numeric_bindings/boost/numeric/bindings/blas/level2/hpr.hpp (original)
+++ sandbox/numeric_bindings/boost/numeric/bindings/blas/level2/hpr.hpp 2012-08-07 09:40:31 EDT (Tue, 07 Aug 2012)
@@ -63,8 +63,8 @@
// * float value-type.
//
template< typename Order, typename UpLo >
-inline void hpr( const Order order, const UpLo uplo, const int n,
- const float alpha, const float* x, const int incx, float* ap ) {
+inline void hpr( const Order, const UpLo, const int n, const float alpha,
+ const float* x, const int incx, float* ap ) {
cblas_sspr( cblas_option< Order >::value, cblas_option< UpLo >::value, n,
alpha, x, incx, ap );
}
@@ -75,8 +75,8 @@
// * double value-type.
//
template< typename Order, typename UpLo >
-inline void hpr( const Order order, const UpLo uplo, const int n,
- const double alpha, const double* x, const int incx, double* ap ) {
+inline void hpr( const Order, const UpLo, const int n, const double alpha,
+ const double* x, const int incx, double* ap ) {
cblas_dspr( cblas_option< Order >::value, cblas_option< UpLo >::value, n,
alpha, x, incx, ap );
}
@@ -87,8 +87,8 @@
// * complex<float> value-type.
//
template< typename Order, typename UpLo >
-inline void hpr( const Order order, const UpLo uplo, const int n,
- const float alpha, const std::complex<float>* x, const int incx,
+inline void hpr( const Order, const UpLo, const int n, const float alpha,
+ const std::complex<float>* x, const int incx,
std::complex<float>* ap ) {
cblas_chpr( cblas_option< Order >::value, cblas_option< UpLo >::value, n,
alpha, x, incx, ap );
@@ -100,8 +100,8 @@
// * complex<double> value-type.
//
template< typename Order, typename UpLo >
-inline void hpr( const Order order, const UpLo uplo, const int n,
- const double alpha, const std::complex<double>* x, const int incx,
+inline void hpr( const Order, const UpLo, const int n, const double alpha,
+ const std::complex<double>* x, const int incx,
std::complex<double>* ap ) {
cblas_zhpr( cblas_option< Order >::value, cblas_option< UpLo >::value, n,
alpha, x, incx, ap );
@@ -114,8 +114,8 @@
// * float value-type.
//
template< typename Order, typename UpLo >
-inline void hpr( const Order order, const UpLo uplo, const int n,
- const float alpha, const float* x, const int incx, float* ap ) {
+inline void hpr( const Order, const UpLo, const int n, const float alpha,
+ const float* x, const int incx, float* ap ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
cublasSspr( blas_option< UpLo >::value, n, alpha, x, incx, ap );
}
@@ -126,10 +126,10 @@
// * double value-type.
//
template< typename Order, typename UpLo >
-inline void hpr( const Order order, const UpLo uplo, const int n,
- const double alpha, const double* x, const int incx, double* ap ) {
+inline void hpr( const Order, const UpLo, const int n, const double alpha,
+ const double* x, const int incx, double* ap ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
- // NOT FOUND();
+ cublasDspr( blas_option< UpLo >::value, n, alpha, x, incx, ap );
}
//
@@ -138,8 +138,8 @@
// * complex<float> value-type.
//
template< typename Order, typename UpLo >
-inline void hpr( const Order order, const UpLo uplo, const int n,
- const float alpha, const std::complex<float>* x, const int incx,
+inline void hpr( const Order, const UpLo, const int n, const float alpha,
+ const std::complex<float>* x, const int incx,
std::complex<float>* ap ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
cublasChpr( blas_option< UpLo >::value, n, alpha, x, incx, ap );
@@ -151,11 +151,11 @@
// * complex<double> value-type.
//
template< typename Order, typename UpLo >
-inline void hpr( const Order order, const UpLo uplo, const int n,
- const double alpha, const std::complex<double>* x, const int incx,
+inline void hpr( const Order, const UpLo, const int n, const double alpha,
+ const std::complex<double>* x, const int incx,
std::complex<double>* ap ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
- // NOT FOUND();
+ cublasZhpr( blas_option< UpLo >::value, n, alpha, x, incx, ap );
}
#else
@@ -165,7 +165,7 @@
// * float value-type.
//
template< typename Order, typename UpLo >
-inline void hpr( const Order order, const UpLo uplo, const fortran_int_t n,
+inline void hpr( const Order, const UpLo, const fortran_int_t n,
const float alpha, const float* x, const fortran_int_t incx,
float* ap ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
@@ -178,7 +178,7 @@
// * double value-type.
//
template< typename Order, typename UpLo >
-inline void hpr( const Order order, const UpLo uplo, const fortran_int_t n,
+inline void hpr( const Order, const UpLo, const fortran_int_t n,
const double alpha, const double* x, const fortran_int_t incx,
double* ap ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
@@ -191,7 +191,7 @@
// * complex<float> value-type.
//
template< typename Order, typename UpLo >
-inline void hpr( const Order order, const UpLo uplo, const fortran_int_t n,
+inline void hpr( const Order, const UpLo, const fortran_int_t n,
const float alpha, const std::complex<float>* x,
const fortran_int_t incx, std::complex<float>* ap ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
@@ -204,7 +204,7 @@
// * complex<double> value-type.
//
template< typename Order, typename UpLo >
-inline void hpr( const Order order, const UpLo uplo, const fortran_int_t n,
+inline void hpr( const Order, const UpLo, const fortran_int_t n,
const double alpha, const std::complex<double>* x,
const fortran_int_t incx, std::complex<double>* ap ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
Modified: sandbox/numeric_bindings/boost/numeric/bindings/blas/level2/hpr2.hpp
==============================================================================
--- sandbox/numeric_bindings/boost/numeric/bindings/blas/level2/hpr2.hpp (original)
+++ sandbox/numeric_bindings/boost/numeric/bindings/blas/level2/hpr2.hpp 2012-08-07 09:40:31 EDT (Tue, 07 Aug 2012)
@@ -63,9 +63,9 @@
// * float value-type.
//
template< typename Order, typename UpLo >
-inline void hpr2( const Order order, const UpLo uplo, const int n,
- const float alpha, const float* x, const int incx, const float* y,
- const int incy, float* ap ) {
+inline void hpr2( const Order, const UpLo, const int n, const float alpha,
+ const float* x, const int incx, const float* y, const int incy,
+ float* ap ) {
cblas_sspr2( cblas_option< Order >::value, cblas_option< UpLo >::value, n,
alpha, x, incx, y, incy, ap );
}
@@ -76,9 +76,9 @@
// * double value-type.
//
template< typename Order, typename UpLo >
-inline void hpr2( const Order order, const UpLo uplo, const int n,
- const double alpha, const double* x, const int incx, const double* y,
- const int incy, double* ap ) {
+inline void hpr2( const Order, const UpLo, const int n, const double alpha,
+ const double* x, const int incx, const double* y, const int incy,
+ double* ap ) {
cblas_dspr2( cblas_option< Order >::value, cblas_option< UpLo >::value, n,
alpha, x, incx, y, incy, ap );
}
@@ -89,7 +89,7 @@
// * complex<float> value-type.
//
template< typename Order, typename UpLo >
-inline void hpr2( const Order order, const UpLo uplo, const int n,
+inline void hpr2( const Order, const UpLo, const int n,
const std::complex<float> alpha, const std::complex<float>* x,
const int incx, const std::complex<float>* y, const int incy,
std::complex<float>* ap ) {
@@ -103,7 +103,7 @@
// * complex<double> value-type.
//
template< typename Order, typename UpLo >
-inline void hpr2( const Order order, const UpLo uplo, const int n,
+inline void hpr2( const Order, const UpLo, const int n,
const std::complex<double> alpha, const std::complex<double>* x,
const int incx, const std::complex<double>* y, const int incy,
std::complex<double>* ap ) {
@@ -118,9 +118,9 @@
// * float value-type.
//
template< typename Order, typename UpLo >
-inline void hpr2( const Order order, const UpLo uplo, const int n,
- const float alpha, const float* x, const int incx, const float* y,
- const int incy, float* ap ) {
+inline void hpr2( const Order, const UpLo, const int n, const float alpha,
+ const float* x, const int incx, const float* y, const int incy,
+ float* ap ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
cublasSspr2( blas_option< UpLo >::value, n, alpha, x, incx, y, incy, ap );
}
@@ -131,11 +131,11 @@
// * double value-type.
//
template< typename Order, typename UpLo >
-inline void hpr2( const Order order, const UpLo uplo, const int n,
- const double alpha, const double* x, const int incx, const double* y,
- const int incy, double* ap ) {
+inline void hpr2( const Order, const UpLo, const int n, const double alpha,
+ const double* x, const int incx, const double* y, const int incy,
+ double* ap ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
- // NOT FOUND();
+ cublasDspr2( blas_option< UpLo >::value, n, alpha, x, incx, y, incy, ap );
}
//
@@ -144,7 +144,7 @@
// * complex<float> value-type.
//
template< typename Order, typename UpLo >
-inline void hpr2( const Order order, const UpLo uplo, const int n,
+inline void hpr2( const Order, const UpLo, const int n,
const std::complex<float> alpha, const std::complex<float>* x,
const int incx, const std::complex<float>* y, const int incy,
std::complex<float>* ap ) {
@@ -158,12 +158,12 @@
// * complex<double> value-type.
//
template< typename Order, typename UpLo >
-inline void hpr2( const Order order, const UpLo uplo, const int n,
+inline void hpr2( const Order, const UpLo, const int n,
const std::complex<double> alpha, const std::complex<double>* x,
const int incx, const std::complex<double>* y, const int incy,
std::complex<double>* ap ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
- // NOT FOUND();
+ cublasZhpr2( blas_option< UpLo >::value, n, alpha, x, incx, y, incy, ap );
}
#else
@@ -173,7 +173,7 @@
// * float value-type.
//
template< typename Order, typename UpLo >
-inline void hpr2( const Order order, const UpLo uplo, const fortran_int_t n,
+inline void hpr2( const Order, const UpLo, const fortran_int_t n,
const float alpha, const float* x, const fortran_int_t incx,
const float* y, const fortran_int_t incy, float* ap ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
@@ -187,7 +187,7 @@
// * double value-type.
//
template< typename Order, typename UpLo >
-inline void hpr2( const Order order, const UpLo uplo, const fortran_int_t n,
+inline void hpr2( const Order, const UpLo, const fortran_int_t n,
const double alpha, const double* x, const fortran_int_t incx,
const double* y, const fortran_int_t incy, double* ap ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
@@ -201,7 +201,7 @@
// * complex<float> value-type.
//
template< typename Order, typename UpLo >
-inline void hpr2( const Order order, const UpLo uplo, const fortran_int_t n,
+inline void hpr2( const Order, const UpLo, const fortran_int_t n,
const std::complex<float> alpha, const std::complex<float>* x,
const fortran_int_t incx, const std::complex<float>* y,
const fortran_int_t incy, std::complex<float>* ap ) {
@@ -216,7 +216,7 @@
// * complex<double> value-type.
//
template< typename Order, typename UpLo >
-inline void hpr2( const Order order, const UpLo uplo, const fortran_int_t n,
+inline void hpr2( const Order, const UpLo, const fortran_int_t n,
const std::complex<double> alpha, const std::complex<double>* x,
const fortran_int_t incx, const std::complex<double>* y,
const fortran_int_t incy, std::complex<double>* ap ) {
Modified: sandbox/numeric_bindings/boost/numeric/bindings/blas/level2/sbmv.hpp
==============================================================================
--- sandbox/numeric_bindings/boost/numeric/bindings/blas/level2/sbmv.hpp (original)
+++ sandbox/numeric_bindings/boost/numeric/bindings/blas/level2/sbmv.hpp 2012-08-07 09:40:31 EDT (Tue, 07 Aug 2012)
@@ -64,10 +64,9 @@
// * float value-type.
//
template< typename Order, typename UpLo >
-inline void sbmv( const Order order, const UpLo uplo, const int n,
- const int k, const float alpha, const float* a, const int lda,
- const float* x, const int incx, const float beta, float* y,
- const int incy ) {
+inline void sbmv( const Order, const UpLo, const int n, const int k,
+ const float alpha, const float* a, const int lda, const float* x,
+ const int incx, const float beta, float* y, const int incy ) {
cblas_ssbmv( cblas_option< Order >::value, cblas_option< UpLo >::value, n,
k, alpha, a, lda, x, incx, beta, y, incy );
}
@@ -78,10 +77,9 @@
// * double value-type.
//
template< typename Order, typename UpLo >
-inline void sbmv( const Order order, const UpLo uplo, const int n,
- const int k, const double alpha, const double* a, const int lda,
- const double* x, const int incx, const double beta, double* y,
- const int incy ) {
+inline void sbmv( const Order, const UpLo, const int n, const int k,
+ const double alpha, const double* a, const int lda, const double* x,
+ const int incx, const double beta, double* y, const int incy ) {
cblas_dsbmv( cblas_option< Order >::value, cblas_option< UpLo >::value, n,
k, alpha, a, lda, x, incx, beta, y, incy );
}
@@ -93,10 +91,9 @@
// * float value-type.
//
template< typename Order, typename UpLo >
-inline void sbmv( const Order order, const UpLo uplo, const int n,
- const int k, const float alpha, const float* a, const int lda,
- const float* x, const int incx, const float beta, float* y,
- const int incy ) {
+inline void sbmv( const Order, const UpLo, const int n, const int k,
+ const float alpha, const float* a, const int lda, const float* x,
+ const int incx, const float beta, float* y, const int incy ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
cublasSsbmv( blas_option< UpLo >::value, n, k, alpha, a, lda, x, incx,
beta, y, incy );
@@ -108,12 +105,12 @@
// * double value-type.
//
template< typename Order, typename UpLo >
-inline void sbmv( const Order order, const UpLo uplo, const int n,
- const int k, const double alpha, const double* a, const int lda,
- const double* x, const int incx, const double beta, double* y,
- const int incy ) {
+inline void sbmv( const Order, const UpLo, const int n, const int k,
+ const double alpha, const double* a, const int lda, const double* x,
+ const int incx, const double beta, double* y, const int incy ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
- // NOT FOUND();
+ cublasDsbmv( blas_option< UpLo >::value, n, k, alpha, a, lda, x, incx,
+ beta, y, incy );
}
#else
@@ -123,7 +120,7 @@
// * float value-type.
//
template< typename Order, typename UpLo >
-inline void sbmv( const Order order, const UpLo uplo, const fortran_int_t n,
+inline void sbmv( const Order, const UpLo, const fortran_int_t n,
const fortran_int_t k, const float alpha, const float* a,
const fortran_int_t lda, const float* x, const fortran_int_t incx,
const float beta, float* y, const fortran_int_t incy ) {
@@ -138,7 +135,7 @@
// * double value-type.
//
template< typename Order, typename UpLo >
-inline void sbmv( const Order order, const UpLo uplo, const fortran_int_t n,
+inline void sbmv( const Order, const UpLo, const fortran_int_t n,
const fortran_int_t k, const double alpha, const double* a,
const fortran_int_t lda, const double* x, const fortran_int_t incx,
const double beta, double* y, const fortran_int_t incy ) {
Modified: sandbox/numeric_bindings/boost/numeric/bindings/blas/level2/spmv.hpp
==============================================================================
--- sandbox/numeric_bindings/boost/numeric/bindings/blas/level2/spmv.hpp (original)
+++ sandbox/numeric_bindings/boost/numeric/bindings/blas/level2/spmv.hpp 2012-08-07 09:40:31 EDT (Tue, 07 Aug 2012)
@@ -63,9 +63,9 @@
// * float value-type.
//
template< typename Order, typename UpLo >
-inline void spmv( const Order order, const UpLo uplo, const int n,
- const float alpha, const float* ap, const float* x, const int incx,
- const float beta, float* y, const int incy ) {
+inline void spmv( const Order, const UpLo, const int n, const float alpha,
+ const float* ap, const float* x, const int incx, const float beta,
+ float* y, const int incy ) {
cblas_sspmv( cblas_option< Order >::value, cblas_option< UpLo >::value, n,
alpha, ap, x, incx, beta, y, incy );
}
@@ -76,9 +76,9 @@
// * double value-type.
//
template< typename Order, typename UpLo >
-inline void spmv( const Order order, const UpLo uplo, const int n,
- const double alpha, const double* ap, const double* x, const int incx,
- const double beta, double* y, const int incy ) {
+inline void spmv( const Order, const UpLo, const int n, const double alpha,
+ const double* ap, const double* x, const int incx, const double beta,
+ double* y, const int incy ) {
cblas_dspmv( cblas_option< Order >::value, cblas_option< UpLo >::value, n,
alpha, ap, x, incx, beta, y, incy );
}
@@ -90,9 +90,9 @@
// * float value-type.
//
template< typename Order, typename UpLo >
-inline void spmv( const Order order, const UpLo uplo, const int n,
- const float alpha, const float* ap, const float* x, const int incx,
- const float beta, float* y, const int incy ) {
+inline void spmv( const Order, const UpLo, const int n, const float alpha,
+ const float* ap, const float* x, const int incx, const float beta,
+ float* y, const int incy ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
cublasSspmv( blas_option< UpLo >::value, n, alpha, ap, x, incx, beta, y,
incy );
@@ -104,11 +104,12 @@
// * double value-type.
//
template< typename Order, typename UpLo >
-inline void spmv( const Order order, const UpLo uplo, const int n,
- const double alpha, const double* ap, const double* x, const int incx,
- const double beta, double* y, const int incy ) {
+inline void spmv( const Order, const UpLo, const int n, const double alpha,
+ const double* ap, const double* x, const int incx, const double beta,
+ double* y, const int incy ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
- // NOT FOUND();
+ cublasDspmv( blas_option< UpLo >::value, n, alpha, ap, x, incx, beta, y,
+ incy );
}
#else
@@ -118,7 +119,7 @@
// * float value-type.
//
template< typename Order, typename UpLo >
-inline void spmv( const Order order, const UpLo uplo, const fortran_int_t n,
+inline void spmv( const Order, const UpLo, const fortran_int_t n,
const float alpha, const float* ap, const float* x,
const fortran_int_t incx, const float beta, float* y,
const fortran_int_t incy ) {
@@ -133,7 +134,7 @@
// * double value-type.
//
template< typename Order, typename UpLo >
-inline void spmv( const Order order, const UpLo uplo, const fortran_int_t n,
+inline void spmv( const Order, const UpLo, const fortran_int_t n,
const double alpha, const double* ap, const double* x,
const fortran_int_t incx, const double beta, double* y,
const fortran_int_t incy ) {
Modified: sandbox/numeric_bindings/boost/numeric/bindings/blas/level2/spr.hpp
==============================================================================
--- sandbox/numeric_bindings/boost/numeric/bindings/blas/level2/spr.hpp (original)
+++ sandbox/numeric_bindings/boost/numeric/bindings/blas/level2/spr.hpp 2012-08-07 09:40:31 EDT (Tue, 07 Aug 2012)
@@ -63,8 +63,8 @@
// * float value-type.
//
template< typename Order, typename UpLo >
-inline void spr( const Order order, const UpLo uplo, const int n,
- const float alpha, const float* x, const int incx, float* ap ) {
+inline void spr( const Order, const UpLo, const int n, const float alpha,
+ const float* x, const int incx, float* ap ) {
cblas_sspr( cblas_option< Order >::value, cblas_option< UpLo >::value, n,
alpha, x, incx, ap );
}
@@ -75,8 +75,8 @@
// * double value-type.
//
template< typename Order, typename UpLo >
-inline void spr( const Order order, const UpLo uplo, const int n,
- const double alpha, const double* x, const int incx, double* ap ) {
+inline void spr( const Order, const UpLo, const int n, const double alpha,
+ const double* x, const int incx, double* ap ) {
cblas_dspr( cblas_option< Order >::value, cblas_option< UpLo >::value, n,
alpha, x, incx, ap );
}
@@ -88,8 +88,8 @@
// * float value-type.
//
template< typename Order, typename UpLo >
-inline void spr( const Order order, const UpLo uplo, const int n,
- const float alpha, const float* x, const int incx, float* ap ) {
+inline void spr( const Order, const UpLo, const int n, const float alpha,
+ const float* x, const int incx, float* ap ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
cublasSspr( blas_option< UpLo >::value, n, alpha, x, incx, ap );
}
@@ -100,10 +100,10 @@
// * double value-type.
//
template< typename Order, typename UpLo >
-inline void spr( const Order order, const UpLo uplo, const int n,
- const double alpha, const double* x, const int incx, double* ap ) {
+inline void spr( const Order, const UpLo, const int n, const double alpha,
+ const double* x, const int incx, double* ap ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
- // NOT FOUND();
+ cublasDspr( blas_option< UpLo >::value, n, alpha, x, incx, ap );
}
#else
@@ -113,7 +113,7 @@
// * float value-type.
//
template< typename Order, typename UpLo >
-inline void spr( const Order order, const UpLo uplo, const fortran_int_t n,
+inline void spr( const Order, const UpLo, const fortran_int_t n,
const float alpha, const float* x, const fortran_int_t incx,
float* ap ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
@@ -126,7 +126,7 @@
// * double value-type.
//
template< typename Order, typename UpLo >
-inline void spr( const Order order, const UpLo uplo, const fortran_int_t n,
+inline void spr( const Order, const UpLo, const fortran_int_t n,
const double alpha, const double* x, const fortran_int_t incx,
double* ap ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
Modified: sandbox/numeric_bindings/boost/numeric/bindings/blas/level2/spr2.hpp
==============================================================================
--- sandbox/numeric_bindings/boost/numeric/bindings/blas/level2/spr2.hpp (original)
+++ sandbox/numeric_bindings/boost/numeric/bindings/blas/level2/spr2.hpp 2012-08-07 09:40:31 EDT (Tue, 07 Aug 2012)
@@ -63,9 +63,9 @@
// * float value-type.
//
template< typename Order, typename UpLo >
-inline void spr2( const Order order, const UpLo uplo, const int n,
- const float alpha, const float* x, const int incx, const float* y,
- const int incy, float* ap ) {
+inline void spr2( const Order, const UpLo, const int n, const float alpha,
+ const float* x, const int incx, const float* y, const int incy,
+ float* ap ) {
cblas_sspr2( cblas_option< Order >::value, cblas_option< UpLo >::value, n,
alpha, x, incx, y, incy, ap );
}
@@ -76,9 +76,9 @@
// * double value-type.
//
template< typename Order, typename UpLo >
-inline void spr2( const Order order, const UpLo uplo, const int n,
- const double alpha, const double* x, const int incx, const double* y,
- const int incy, double* ap ) {
+inline void spr2( const Order, const UpLo, const int n, const double alpha,
+ const double* x, const int incx, const double* y, const int incy,
+ double* ap ) {
cblas_dspr2( cblas_option< Order >::value, cblas_option< UpLo >::value, n,
alpha, x, incx, y, incy, ap );
}
@@ -90,9 +90,9 @@
// * float value-type.
//
template< typename Order, typename UpLo >
-inline void spr2( const Order order, const UpLo uplo, const int n,
- const float alpha, const float* x, const int incx, const float* y,
- const int incy, float* ap ) {
+inline void spr2( const Order, const UpLo, const int n, const float alpha,
+ const float* x, const int incx, const float* y, const int incy,
+ float* ap ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
cublasSspr2( blas_option< UpLo >::value, n, alpha, x, incx, y, incy, ap );
}
@@ -103,11 +103,11 @@
// * double value-type.
//
template< typename Order, typename UpLo >
-inline void spr2( const Order order, const UpLo uplo, const int n,
- const double alpha, const double* x, const int incx, const double* y,
- const int incy, double* ap ) {
+inline void spr2( const Order, const UpLo, const int n, const double alpha,
+ const double* x, const int incx, const double* y, const int incy,
+ double* ap ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
- // NOT FOUND();
+ cublasDspr2( blas_option< UpLo >::value, n, alpha, x, incx, y, incy, ap );
}
#else
@@ -117,7 +117,7 @@
// * float value-type.
//
template< typename Order, typename UpLo >
-inline void spr2( const Order order, const UpLo uplo, const fortran_int_t n,
+inline void spr2( const Order, const UpLo, const fortran_int_t n,
const float alpha, const float* x, const fortran_int_t incx,
const float* y, const fortran_int_t incy, float* ap ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
@@ -131,7 +131,7 @@
// * double value-type.
//
template< typename Order, typename UpLo >
-inline void spr2( const Order order, const UpLo uplo, const fortran_int_t n,
+inline void spr2( const Order, const UpLo, const fortran_int_t n,
const double alpha, const double* x, const fortran_int_t incx,
const double* y, const fortran_int_t incy, double* ap ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
Modified: sandbox/numeric_bindings/boost/numeric/bindings/blas/level2/symv.hpp
==============================================================================
--- sandbox/numeric_bindings/boost/numeric/bindings/blas/level2/symv.hpp (original)
+++ sandbox/numeric_bindings/boost/numeric/bindings/blas/level2/symv.hpp 2012-08-07 09:40:31 EDT (Tue, 07 Aug 2012)
@@ -63,9 +63,9 @@
// * float value-type.
//
template< typename Order, typename UpLo >
-inline void symv( const Order order, const UpLo uplo, const int n,
- const float alpha, const float* a, const int lda, const float* x,
- const int incx, const float beta, float* y, const int incy ) {
+inline void symv( const Order, const UpLo, const int n, const float alpha,
+ const float* a, const int lda, const float* x, const int incx,
+ const float beta, float* y, const int incy ) {
cblas_ssymv( cblas_option< Order >::value, cblas_option< UpLo >::value, n,
alpha, a, lda, x, incx, beta, y, incy );
}
@@ -76,9 +76,9 @@
// * double value-type.
//
template< typename Order, typename UpLo >
-inline void symv( const Order order, const UpLo uplo, const int n,
- const double alpha, const double* a, const int lda, const double* x,
- const int incx, const double beta, double* y, const int incy ) {
+inline void symv( const Order, const UpLo, const int n, const double alpha,
+ const double* a, const int lda, const double* x, const int incx,
+ const double beta, double* y, const int incy ) {
cblas_dsymv( cblas_option< Order >::value, cblas_option< UpLo >::value, n,
alpha, a, lda, x, incx, beta, y, incy );
}
@@ -90,9 +90,9 @@
// * float value-type.
//
template< typename Order, typename UpLo >
-inline void symv( const Order order, const UpLo uplo, const int n,
- const float alpha, const float* a, const int lda, const float* x,
- const int incx, const float beta, float* y, const int incy ) {
+inline void symv( const Order, const UpLo, const int n, const float alpha,
+ const float* a, const int lda, const float* x, const int incx,
+ const float beta, float* y, const int incy ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
cublasSsymv( blas_option< UpLo >::value, n, alpha, a, lda, x, incx, beta,
y, incy );
@@ -104,11 +104,12 @@
// * double value-type.
//
template< typename Order, typename UpLo >
-inline void symv( const Order order, const UpLo uplo, const int n,
- const double alpha, const double* a, const int lda, const double* x,
- const int incx, const double beta, double* y, const int incy ) {
+inline void symv( const Order, const UpLo, const int n, const double alpha,
+ const double* a, const int lda, const double* x, const int incx,
+ const double beta, double* y, const int incy ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
- // NOT FOUND();
+ cublasDsymv( blas_option< UpLo >::value, n, alpha, a, lda, x, incx, beta,
+ y, incy );
}
#else
@@ -118,7 +119,7 @@
// * float value-type.
//
template< typename Order, typename UpLo >
-inline void symv( const Order order, const UpLo uplo, const fortran_int_t n,
+inline void symv( const Order, const UpLo, const fortran_int_t n,
const float alpha, const float* a, const fortran_int_t lda,
const float* x, const fortran_int_t incx, const float beta, float* y,
const fortran_int_t incy ) {
@@ -133,7 +134,7 @@
// * double value-type.
//
template< typename Order, typename UpLo >
-inline void symv( const Order order, const UpLo uplo, const fortran_int_t n,
+inline void symv( const Order, const UpLo, const fortran_int_t n,
const double alpha, const double* a, const fortran_int_t lda,
const double* x, const fortran_int_t incx, const double beta,
double* y, const fortran_int_t incy ) {
Modified: sandbox/numeric_bindings/boost/numeric/bindings/blas/level2/syr.hpp
==============================================================================
--- sandbox/numeric_bindings/boost/numeric/bindings/blas/level2/syr.hpp (original)
+++ sandbox/numeric_bindings/boost/numeric/bindings/blas/level2/syr.hpp 2012-08-07 09:40:31 EDT (Tue, 07 Aug 2012)
@@ -63,9 +63,8 @@
// * float value-type.
//
template< typename Order, typename UpLo >
-inline void syr( const Order order, const UpLo uplo, const int n,
- const float alpha, const float* x, const int incx, float* a,
- const int lda ) {
+inline void syr( const Order, const UpLo, const int n, const float alpha,
+ const float* x, const int incx, float* a, const int lda ) {
cblas_ssyr( cblas_option< Order >::value, cblas_option< UpLo >::value, n,
alpha, x, incx, a, lda );
}
@@ -76,9 +75,8 @@
// * double value-type.
//
template< typename Order, typename UpLo >
-inline void syr( const Order order, const UpLo uplo, const int n,
- const double alpha, const double* x, const int incx, double* a,
- const int lda ) {
+inline void syr( const Order, const UpLo, const int n, const double alpha,
+ const double* x, const int incx, double* a, const int lda ) {
cblas_dsyr( cblas_option< Order >::value, cblas_option< UpLo >::value, n,
alpha, x, incx, a, lda );
}
@@ -90,9 +88,8 @@
// * float value-type.
//
template< typename Order, typename UpLo >
-inline void syr( const Order order, const UpLo uplo, const int n,
- const float alpha, const float* x, const int incx, float* a,
- const int lda ) {
+inline void syr( const Order, const UpLo, const int n, const float alpha,
+ const float* x, const int incx, float* a, const int lda ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
cublasSsyr( blas_option< UpLo >::value, n, alpha, x, incx, a, lda );
}
@@ -103,9 +100,8 @@
// * double value-type.
//
template< typename Order, typename UpLo >
-inline void syr( const Order order, const UpLo uplo, const int n,
- const double alpha, const double* x, const int incx, double* a,
- const int lda ) {
+inline void syr( const Order, const UpLo, const int n, const double alpha,
+ const double* x, const int incx, double* a, const int lda ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
cublasDsyr( blas_option< UpLo >::value, n, alpha, x, incx, a, lda );
}
@@ -117,7 +113,7 @@
// * float value-type.
//
template< typename Order, typename UpLo >
-inline void syr( const Order order, const UpLo uplo, const fortran_int_t n,
+inline void syr( const Order, const UpLo, const fortran_int_t n,
const float alpha, const float* x, const fortran_int_t incx, float* a,
const fortran_int_t lda ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
@@ -130,7 +126,7 @@
// * double value-type.
//
template< typename Order, typename UpLo >
-inline void syr( const Order order, const UpLo uplo, const fortran_int_t n,
+inline void syr( const Order, const UpLo, const fortran_int_t n,
const double alpha, const double* x, const fortran_int_t incx,
double* a, const fortran_int_t lda ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
Modified: sandbox/numeric_bindings/boost/numeric/bindings/blas/level2/syr2.hpp
==============================================================================
--- sandbox/numeric_bindings/boost/numeric/bindings/blas/level2/syr2.hpp (original)
+++ sandbox/numeric_bindings/boost/numeric/bindings/blas/level2/syr2.hpp 2012-08-07 09:40:31 EDT (Tue, 07 Aug 2012)
@@ -63,9 +63,9 @@
// * float value-type.
//
template< typename Order, typename UpLo >
-inline void syr2( const Order order, const UpLo uplo, const int n,
- const float alpha, const float* x, const int incx, const float* y,
- const int incy, float* a, const int lda ) {
+inline void syr2( const Order, const UpLo, const int n, const float alpha,
+ const float* x, const int incx, const float* y, const int incy,
+ float* a, const int lda ) {
cblas_ssyr2( cblas_option< Order >::value, cblas_option< UpLo >::value, n,
alpha, x, incx, y, incy, a, lda );
}
@@ -76,9 +76,9 @@
// * double value-type.
//
template< typename Order, typename UpLo >
-inline void syr2( const Order order, const UpLo uplo, const int n,
- const double alpha, const double* x, const int incx, const double* y,
- const int incy, double* a, const int lda ) {
+inline void syr2( const Order, const UpLo, const int n, const double alpha,
+ const double* x, const int incx, const double* y, const int incy,
+ double* a, const int lda ) {
cblas_dsyr2( cblas_option< Order >::value, cblas_option< UpLo >::value, n,
alpha, x, incx, y, incy, a, lda );
}
@@ -90,9 +90,9 @@
// * float value-type.
//
template< typename Order, typename UpLo >
-inline void syr2( const Order order, const UpLo uplo, const int n,
- const float alpha, const float* x, const int incx, const float* y,
- const int incy, float* a, const int lda ) {
+inline void syr2( const Order, const UpLo, const int n, const float alpha,
+ const float* x, const int incx, const float* y, const int incy,
+ float* a, const int lda ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
cublasSsyr2( blas_option< UpLo >::value, n, alpha, x, incx, y, incy, a,
lda );
@@ -104,11 +104,12 @@
// * double value-type.
//
template< typename Order, typename UpLo >
-inline void syr2( const Order order, const UpLo uplo, const int n,
- const double alpha, const double* x, const int incx, const double* y,
- const int incy, double* a, const int lda ) {
+inline void syr2( const Order, const UpLo, const int n, const double alpha,
+ const double* x, const int incx, const double* y, const int incy,
+ double* a, const int lda ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
- // NOT FOUND();
+ cublasDsyr2( blas_option< UpLo >::value, n, alpha, x, incx, y, incy, a,
+ lda );
}
#else
@@ -118,7 +119,7 @@
// * float value-type.
//
template< typename Order, typename UpLo >
-inline void syr2( const Order order, const UpLo uplo, const fortran_int_t n,
+inline void syr2( const Order, const UpLo, const fortran_int_t n,
const float alpha, const float* x, const fortran_int_t incx,
const float* y, const fortran_int_t incy, float* a,
const fortran_int_t lda ) {
@@ -133,7 +134,7 @@
// * double value-type.
//
template< typename Order, typename UpLo >
-inline void syr2( const Order order, const UpLo uplo, const fortran_int_t n,
+inline void syr2( const Order, const UpLo, const fortran_int_t n,
const double alpha, const double* x, const fortran_int_t incx,
const double* y, const fortran_int_t incy, double* a,
const fortran_int_t lda ) {
Modified: sandbox/numeric_bindings/boost/numeric/bindings/blas/level2/tbmv.hpp
==============================================================================
--- sandbox/numeric_bindings/boost/numeric/bindings/blas/level2/tbmv.hpp (original)
+++ sandbox/numeric_bindings/boost/numeric/bindings/blas/level2/tbmv.hpp 2012-08-07 09:40:31 EDT (Tue, 07 Aug 2012)
@@ -66,9 +66,9 @@
// * float value-type.
//
template< typename Order, typename UpLo, typename Trans, typename Diag >
-inline void tbmv( const Order order, const UpLo uplo, const Trans trans,
- const Diag diag, const int n, const int k, const float* a,
- const int lda, float* x, const int incx ) {
+inline void tbmv( const Order, const UpLo, const Trans, const Diag,
+ const int n, const int k, const float* a, const int lda, float* x,
+ const int incx ) {
cblas_stbmv( cblas_option< Order >::value, cblas_option< UpLo >::value,
cblas_option< Trans >::value, cblas_option< Diag >::value, n, k,
a, lda, x, incx );
@@ -80,9 +80,9 @@
// * double value-type.
//
template< typename Order, typename UpLo, typename Trans, typename Diag >
-inline void tbmv( const Order order, const UpLo uplo, const Trans trans,
- const Diag diag, const int n, const int k, const double* a,
- const int lda, double* x, const int incx ) {
+inline void tbmv( const Order, const UpLo, const Trans, const Diag,
+ const int n, const int k, const double* a, const int lda, double* x,
+ const int incx ) {
cblas_dtbmv( cblas_option< Order >::value, cblas_option< UpLo >::value,
cblas_option< Trans >::value, cblas_option< Diag >::value, n, k,
a, lda, x, incx );
@@ -94,10 +94,9 @@
// * complex<float> value-type.
//
template< typename Order, typename UpLo, typename Trans, typename Diag >
-inline void tbmv( const Order order, const UpLo uplo, const Trans trans,
- const Diag diag, const int n, const int k,
- const std::complex<float>* a, const int lda, std::complex<float>* x,
- const int incx ) {
+inline void tbmv( const Order, const UpLo, const Trans, const Diag,
+ const int n, const int k, const std::complex<float>* a, const int lda,
+ std::complex<float>* x, const int incx ) {
cblas_ctbmv( cblas_option< Order >::value, cblas_option< UpLo >::value,
cblas_option< Trans >::value, cblas_option< Diag >::value, n, k,
a, lda, x, incx );
@@ -109,10 +108,9 @@
// * complex<double> value-type.
//
template< typename Order, typename UpLo, typename Trans, typename Diag >
-inline void tbmv( const Order order, const UpLo uplo, const Trans trans,
- const Diag diag, const int n, const int k,
- const std::complex<double>* a, const int lda, std::complex<double>* x,
- const int incx ) {
+inline void tbmv( const Order, const UpLo, const Trans, const Diag,
+ const int n, const int k, const std::complex<double>* a,
+ const int lda, std::complex<double>* x, const int incx ) {
cblas_ztbmv( cblas_option< Order >::value, cblas_option< UpLo >::value,
cblas_option< Trans >::value, cblas_option< Diag >::value, n, k,
a, lda, x, incx );
@@ -125,9 +123,9 @@
// * float value-type.
//
template< typename Order, typename UpLo, typename Trans, typename Diag >
-inline void tbmv( const Order order, const UpLo uplo, const Trans trans,
- const Diag diag, const int n, const int k, const float* a,
- const int lda, float* x, const int incx ) {
+inline void tbmv( const Order, const UpLo, const Trans, const Diag,
+ const int n, const int k, const float* a, const int lda, float* x,
+ const int incx ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
cublasStbmv( blas_option< UpLo >::value, blas_option< Trans >::value,
blas_option< Diag >::value, n, k, a, lda, x, incx );
@@ -139,11 +137,12 @@
// * double value-type.
//
template< typename Order, typename UpLo, typename Trans, typename Diag >
-inline void tbmv( const Order order, const UpLo uplo, const Trans trans,
- const Diag diag, const int n, const int k, const double* a,
- const int lda, double* x, const int incx ) {
+inline void tbmv( const Order, const UpLo, const Trans, const Diag,
+ const int n, const int k, const double* a, const int lda, double* x,
+ const int incx ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
- // NOT FOUND();
+ cublasDtbmv( blas_option< UpLo >::value, blas_option< Trans >::value,
+ blas_option< Diag >::value, n, k, a, lda, x, incx );
}
//
@@ -152,10 +151,9 @@
// * complex<float> value-type.
//
template< typename Order, typename UpLo, typename Trans, typename Diag >
-inline void tbmv( const Order order, const UpLo uplo, const Trans trans,
- const Diag diag, const int n, const int k,
- const std::complex<float>* a, const int lda, std::complex<float>* x,
- const int incx ) {
+inline void tbmv( const Order, const UpLo, const Trans, const Diag,
+ const int n, const int k, const std::complex<float>* a, const int lda,
+ std::complex<float>* x, const int incx ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
cublasCtbmv( blas_option< UpLo >::value, blas_option< Trans >::value,
blas_option< Diag >::value, n, k, a, lda, x, incx );
@@ -167,12 +165,12 @@
// * complex<double> value-type.
//
template< typename Order, typename UpLo, typename Trans, typename Diag >
-inline void tbmv( const Order order, const UpLo uplo, const Trans trans,
- const Diag diag, const int n, const int k,
- const std::complex<double>* a, const int lda, std::complex<double>* x,
- const int incx ) {
+inline void tbmv( const Order, const UpLo, const Trans, const Diag,
+ const int n, const int k, const std::complex<double>* a,
+ const int lda, std::complex<double>* x, const int incx ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
- // NOT FOUND();
+ cublasZtbmv( blas_option< UpLo >::value, blas_option< Trans >::value,
+ blas_option< Diag >::value, n, k, a, lda, x, incx );
}
#else
@@ -182,10 +180,9 @@
// * float value-type.
//
template< typename Order, typename UpLo, typename Trans, typename Diag >
-inline void tbmv( const Order order, const UpLo uplo, const Trans trans,
- const Diag diag, const fortran_int_t n, const fortran_int_t k,
- const float* a, const fortran_int_t lda, float* x,
- const fortran_int_t incx ) {
+inline void tbmv( const Order, const UpLo, const Trans, const Diag,
+ const fortran_int_t n, const fortran_int_t k, const float* a,
+ const fortran_int_t lda, float* x, const fortran_int_t incx ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
BLAS_STBMV( &blas_option< UpLo >::value, &blas_option< Trans >::value,
&blas_option< Diag >::value, &n, &k, a, &lda, x, &incx );
@@ -197,10 +194,9 @@
// * double value-type.
//
template< typename Order, typename UpLo, typename Trans, typename Diag >
-inline void tbmv( const Order order, const UpLo uplo, const Trans trans,
- const Diag diag, const fortran_int_t n, const fortran_int_t k,
- const double* a, const fortran_int_t lda, double* x,
- const fortran_int_t incx ) {
+inline void tbmv( const Order, const UpLo, const Trans, const Diag,
+ const fortran_int_t n, const fortran_int_t k, const double* a,
+ const fortran_int_t lda, double* x, const fortran_int_t incx ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
BLAS_DTBMV( &blas_option< UpLo >::value, &blas_option< Trans >::value,
&blas_option< Diag >::value, &n, &k, a, &lda, x, &incx );
@@ -212,8 +208,8 @@
// * complex<float> value-type.
//
template< typename Order, typename UpLo, typename Trans, typename Diag >
-inline void tbmv( const Order order, const UpLo uplo, const Trans trans,
- const Diag diag, const fortran_int_t n, const fortran_int_t k,
+inline void tbmv( const Order, const UpLo, const Trans, const Diag,
+ const fortran_int_t n, const fortran_int_t k,
const std::complex<float>* a, const fortran_int_t lda,
std::complex<float>* x, const fortran_int_t incx ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
@@ -227,8 +223,8 @@
// * complex<double> value-type.
//
template< typename Order, typename UpLo, typename Trans, typename Diag >
-inline void tbmv( const Order order, const UpLo uplo, const Trans trans,
- const Diag diag, const fortran_int_t n, const fortran_int_t k,
+inline void tbmv( const Order, const UpLo, const Trans, const Diag,
+ const fortran_int_t n, const fortran_int_t k,
const std::complex<double>* a, const fortran_int_t lda,
std::complex<double>* x, const fortran_int_t incx ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
Modified: sandbox/numeric_bindings/boost/numeric/bindings/blas/level2/tbsv.hpp
==============================================================================
--- sandbox/numeric_bindings/boost/numeric/bindings/blas/level2/tbsv.hpp (original)
+++ sandbox/numeric_bindings/boost/numeric/bindings/blas/level2/tbsv.hpp 2012-08-07 09:40:31 EDT (Tue, 07 Aug 2012)
@@ -66,9 +66,9 @@
// * float value-type.
//
template< typename Order, typename UpLo, typename Trans, typename Diag >
-inline void tbsv( const Order order, const UpLo uplo, const Trans trans,
- const Diag diag, const int n, const int k, const float* a,
- const int lda, float* x, const int incx ) {
+inline void tbsv( const Order, const UpLo, const Trans, const Diag,
+ const int n, const int k, const float* a, const int lda, float* x,
+ const int incx ) {
cblas_stbsv( cblas_option< Order >::value, cblas_option< UpLo >::value,
cblas_option< Trans >::value, cblas_option< Diag >::value, n, k,
a, lda, x, incx );
@@ -80,9 +80,9 @@
// * double value-type.
//
template< typename Order, typename UpLo, typename Trans, typename Diag >
-inline void tbsv( const Order order, const UpLo uplo, const Trans trans,
- const Diag diag, const int n, const int k, const double* a,
- const int lda, double* x, const int incx ) {
+inline void tbsv( const Order, const UpLo, const Trans, const Diag,
+ const int n, const int k, const double* a, const int lda, double* x,
+ const int incx ) {
cblas_dtbsv( cblas_option< Order >::value, cblas_option< UpLo >::value,
cblas_option< Trans >::value, cblas_option< Diag >::value, n, k,
a, lda, x, incx );
@@ -94,10 +94,9 @@
// * complex<float> value-type.
//
template< typename Order, typename UpLo, typename Trans, typename Diag >
-inline void tbsv( const Order order, const UpLo uplo, const Trans trans,
- const Diag diag, const int n, const int k,
- const std::complex<float>* a, const int lda, std::complex<float>* x,
- const int incx ) {
+inline void tbsv( const Order, const UpLo, const Trans, const Diag,
+ const int n, const int k, const std::complex<float>* a, const int lda,
+ std::complex<float>* x, const int incx ) {
cblas_ctbsv( cblas_option< Order >::value, cblas_option< UpLo >::value,
cblas_option< Trans >::value, cblas_option< Diag >::value, n, k,
a, lda, x, incx );
@@ -109,10 +108,9 @@
// * complex<double> value-type.
//
template< typename Order, typename UpLo, typename Trans, typename Diag >
-inline void tbsv( const Order order, const UpLo uplo, const Trans trans,
- const Diag diag, const int n, const int k,
- const std::complex<double>* a, const int lda, std::complex<double>* x,
- const int incx ) {
+inline void tbsv( const Order, const UpLo, const Trans, const Diag,
+ const int n, const int k, const std::complex<double>* a,
+ const int lda, std::complex<double>* x, const int incx ) {
cblas_ztbsv( cblas_option< Order >::value, cblas_option< UpLo >::value,
cblas_option< Trans >::value, cblas_option< Diag >::value, n, k,
a, lda, x, incx );
@@ -125,9 +123,9 @@
// * float value-type.
//
template< typename Order, typename UpLo, typename Trans, typename Diag >
-inline void tbsv( const Order order, const UpLo uplo, const Trans trans,
- const Diag diag, const int n, const int k, const float* a,
- const int lda, float* x, const int incx ) {
+inline void tbsv( const Order, const UpLo, const Trans, const Diag,
+ const int n, const int k, const float* a, const int lda, float* x,
+ const int incx ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
cublasStbsv( blas_option< UpLo >::value, blas_option< Trans >::value,
blas_option< Diag >::value, n, k, a, lda, x, incx );
@@ -139,11 +137,12 @@
// * double value-type.
//
template< typename Order, typename UpLo, typename Trans, typename Diag >
-inline void tbsv( const Order order, const UpLo uplo, const Trans trans,
- const Diag diag, const int n, const int k, const double* a,
- const int lda, double* x, const int incx ) {
+inline void tbsv( const Order, const UpLo, const Trans, const Diag,
+ const int n, const int k, const double* a, const int lda, double* x,
+ const int incx ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
- // NOT FOUND();
+ cublasDtbsv( blas_option< UpLo >::value, blas_option< Trans >::value,
+ blas_option< Diag >::value, n, k, a, lda, x, incx );
}
//
@@ -152,10 +151,9 @@
// * complex<float> value-type.
//
template< typename Order, typename UpLo, typename Trans, typename Diag >
-inline void tbsv( const Order order, const UpLo uplo, const Trans trans,
- const Diag diag, const int n, const int k,
- const std::complex<float>* a, const int lda, std::complex<float>* x,
- const int incx ) {
+inline void tbsv( const Order, const UpLo, const Trans, const Diag,
+ const int n, const int k, const std::complex<float>* a, const int lda,
+ std::complex<float>* x, const int incx ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
cublasCtbsv( blas_option< UpLo >::value, blas_option< Trans >::value,
blas_option< Diag >::value, n, k, a, lda, x, incx );
@@ -167,12 +165,12 @@
// * complex<double> value-type.
//
template< typename Order, typename UpLo, typename Trans, typename Diag >
-inline void tbsv( const Order order, const UpLo uplo, const Trans trans,
- const Diag diag, const int n, const int k,
- const std::complex<double>* a, const int lda, std::complex<double>* x,
- const int incx ) {
+inline void tbsv( const Order, const UpLo, const Trans, const Diag,
+ const int n, const int k, const std::complex<double>* a,
+ const int lda, std::complex<double>* x, const int incx ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
- // NOT FOUND();
+ cublasZtbsv( blas_option< UpLo >::value, blas_option< Trans >::value,
+ blas_option< Diag >::value, n, k, a, lda, x, incx );
}
#else
@@ -182,10 +180,9 @@
// * float value-type.
//
template< typename Order, typename UpLo, typename Trans, typename Diag >
-inline void tbsv( const Order order, const UpLo uplo, const Trans trans,
- const Diag diag, const fortran_int_t n, const fortran_int_t k,
- const float* a, const fortran_int_t lda, float* x,
- const fortran_int_t incx ) {
+inline void tbsv( const Order, const UpLo, const Trans, const Diag,
+ const fortran_int_t n, const fortran_int_t k, const float* a,
+ const fortran_int_t lda, float* x, const fortran_int_t incx ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
BLAS_STBSV( &blas_option< UpLo >::value, &blas_option< Trans >::value,
&blas_option< Diag >::value, &n, &k, a, &lda, x, &incx );
@@ -197,10 +194,9 @@
// * double value-type.
//
template< typename Order, typename UpLo, typename Trans, typename Diag >
-inline void tbsv( const Order order, const UpLo uplo, const Trans trans,
- const Diag diag, const fortran_int_t n, const fortran_int_t k,
- const double* a, const fortran_int_t lda, double* x,
- const fortran_int_t incx ) {
+inline void tbsv( const Order, const UpLo, const Trans, const Diag,
+ const fortran_int_t n, const fortran_int_t k, const double* a,
+ const fortran_int_t lda, double* x, const fortran_int_t incx ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
BLAS_DTBSV( &blas_option< UpLo >::value, &blas_option< Trans >::value,
&blas_option< Diag >::value, &n, &k, a, &lda, x, &incx );
@@ -212,8 +208,8 @@
// * complex<float> value-type.
//
template< typename Order, typename UpLo, typename Trans, typename Diag >
-inline void tbsv( const Order order, const UpLo uplo, const Trans trans,
- const Diag diag, const fortran_int_t n, const fortran_int_t k,
+inline void tbsv( const Order, const UpLo, const Trans, const Diag,
+ const fortran_int_t n, const fortran_int_t k,
const std::complex<float>* a, const fortran_int_t lda,
std::complex<float>* x, const fortran_int_t incx ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
@@ -227,8 +223,8 @@
// * complex<double> value-type.
//
template< typename Order, typename UpLo, typename Trans, typename Diag >
-inline void tbsv( const Order order, const UpLo uplo, const Trans trans,
- const Diag diag, const fortran_int_t n, const fortran_int_t k,
+inline void tbsv( const Order, const UpLo, const Trans, const Diag,
+ const fortran_int_t n, const fortran_int_t k,
const std::complex<double>* a, const fortran_int_t lda,
std::complex<double>* x, const fortran_int_t incx ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
Modified: sandbox/numeric_bindings/boost/numeric/bindings/blas/level2/tpmv.hpp
==============================================================================
--- sandbox/numeric_bindings/boost/numeric/bindings/blas/level2/tpmv.hpp (original)
+++ sandbox/numeric_bindings/boost/numeric/bindings/blas/level2/tpmv.hpp 2012-08-07 09:40:31 EDT (Tue, 07 Aug 2012)
@@ -65,9 +65,8 @@
// * float value-type.
//
template< typename Order, typename UpLo, typename Trans, typename Diag >
-inline void tpmv( const Order order, const UpLo uplo, const Trans trans,
- const Diag diag, const int n, const float* ap, float* x,
- const int incx ) {
+inline void tpmv( const Order, const UpLo, const Trans, const Diag,
+ const int n, const float* ap, float* x, const int incx ) {
cblas_stpmv( cblas_option< Order >::value, cblas_option< UpLo >::value,
cblas_option< Trans >::value, cblas_option< Diag >::value, n, ap,
x, incx );
@@ -79,9 +78,8 @@
// * double value-type.
//
template< typename Order, typename UpLo, typename Trans, typename Diag >
-inline void tpmv( const Order order, const UpLo uplo, const Trans trans,
- const Diag diag, const int n, const double* ap, double* x,
- const int incx ) {
+inline void tpmv( const Order, const UpLo, const Trans, const Diag,
+ const int n, const double* ap, double* x, const int incx ) {
cblas_dtpmv( cblas_option< Order >::value, cblas_option< UpLo >::value,
cblas_option< Trans >::value, cblas_option< Diag >::value, n, ap,
x, incx );
@@ -93,9 +91,9 @@
// * complex<float> value-type.
//
template< typename Order, typename UpLo, typename Trans, typename Diag >
-inline void tpmv( const Order order, const UpLo uplo, const Trans trans,
- const Diag diag, const int n, const std::complex<float>* ap,
- std::complex<float>* x, const int incx ) {
+inline void tpmv( const Order, const UpLo, const Trans, const Diag,
+ const int n, const std::complex<float>* ap, std::complex<float>* x,
+ const int incx ) {
cblas_ctpmv( cblas_option< Order >::value, cblas_option< UpLo >::value,
cblas_option< Trans >::value, cblas_option< Diag >::value, n, ap,
x, incx );
@@ -107,9 +105,9 @@
// * complex<double> value-type.
//
template< typename Order, typename UpLo, typename Trans, typename Diag >
-inline void tpmv( const Order order, const UpLo uplo, const Trans trans,
- const Diag diag, const int n, const std::complex<double>* ap,
- std::complex<double>* x, const int incx ) {
+inline void tpmv( const Order, const UpLo, const Trans, const Diag,
+ const int n, const std::complex<double>* ap, std::complex<double>* x,
+ const int incx ) {
cblas_ztpmv( cblas_option< Order >::value, cblas_option< UpLo >::value,
cblas_option< Trans >::value, cblas_option< Diag >::value, n, ap,
x, incx );
@@ -122,9 +120,8 @@
// * float value-type.
//
template< typename Order, typename UpLo, typename Trans, typename Diag >
-inline void tpmv( const Order order, const UpLo uplo, const Trans trans,
- const Diag diag, const int n, const float* ap, float* x,
- const int incx ) {
+inline void tpmv( const Order, const UpLo, const Trans, const Diag,
+ const int n, const float* ap, float* x, const int incx ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
cublasStpmv( blas_option< UpLo >::value, blas_option< Trans >::value,
blas_option< Diag >::value, n, ap, x, incx );
@@ -136,11 +133,11 @@
// * double value-type.
//
template< typename Order, typename UpLo, typename Trans, typename Diag >
-inline void tpmv( const Order order, const UpLo uplo, const Trans trans,
- const Diag diag, const int n, const double* ap, double* x,
- const int incx ) {
+inline void tpmv( const Order, const UpLo, const Trans, const Diag,
+ const int n, const double* ap, double* x, const int incx ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
- // NOT FOUND();
+ cublasDtpmv( blas_option< UpLo >::value, blas_option< Trans >::value,
+ blas_option< Diag >::value, n, ap, x, incx );
}
//
@@ -149,9 +146,9 @@
// * complex<float> value-type.
//
template< typename Order, typename UpLo, typename Trans, typename Diag >
-inline void tpmv( const Order order, const UpLo uplo, const Trans trans,
- const Diag diag, const int n, const std::complex<float>* ap,
- std::complex<float>* x, const int incx ) {
+inline void tpmv( const Order, const UpLo, const Trans, const Diag,
+ const int n, const std::complex<float>* ap, std::complex<float>* x,
+ const int incx ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
cublasCtpmv( blas_option< UpLo >::value, blas_option< Trans >::value,
blas_option< Diag >::value, n, ap, x, incx );
@@ -163,11 +160,12 @@
// * complex<double> value-type.
//
template< typename Order, typename UpLo, typename Trans, typename Diag >
-inline void tpmv( const Order order, const UpLo uplo, const Trans trans,
- const Diag diag, const int n, const std::complex<double>* ap,
- std::complex<double>* x, const int incx ) {
+inline void tpmv( const Order, const UpLo, const Trans, const Diag,
+ const int n, const std::complex<double>* ap, std::complex<double>* x,
+ const int incx ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
- // NOT FOUND();
+ cublasZtpmv( blas_option< UpLo >::value, blas_option< Trans >::value,
+ blas_option< Diag >::value, n, ap, x, incx );
}
#else
@@ -177,8 +175,8 @@
// * float value-type.
//
template< typename Order, typename UpLo, typename Trans, typename Diag >
-inline void tpmv( const Order order, const UpLo uplo, const Trans trans,
- const Diag diag, const fortran_int_t n, const float* ap, float* x,
+inline void tpmv( const Order, const UpLo, const Trans, const Diag,
+ const fortran_int_t n, const float* ap, float* x,
const fortran_int_t incx ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
BLAS_STPMV( &blas_option< UpLo >::value, &blas_option< Trans >::value,
@@ -191,8 +189,8 @@
// * double value-type.
//
template< typename Order, typename UpLo, typename Trans, typename Diag >
-inline void tpmv( const Order order, const UpLo uplo, const Trans trans,
- const Diag diag, const fortran_int_t n, const double* ap, double* x,
+inline void tpmv( const Order, const UpLo, const Trans, const Diag,
+ const fortran_int_t n, const double* ap, double* x,
const fortran_int_t incx ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
BLAS_DTPMV( &blas_option< UpLo >::value, &blas_option< Trans >::value,
@@ -205,8 +203,8 @@
// * complex<float> value-type.
//
template< typename Order, typename UpLo, typename Trans, typename Diag >
-inline void tpmv( const Order order, const UpLo uplo, const Trans trans,
- const Diag diag, const fortran_int_t n, const std::complex<float>* ap,
+inline void tpmv( const Order, const UpLo, const Trans, const Diag,
+ const fortran_int_t n, const std::complex<float>* ap,
std::complex<float>* x, const fortran_int_t incx ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
BLAS_CTPMV( &blas_option< UpLo >::value, &blas_option< Trans >::value,
@@ -219,10 +217,9 @@
// * complex<double> value-type.
//
template< typename Order, typename UpLo, typename Trans, typename Diag >
-inline void tpmv( const Order order, const UpLo uplo, const Trans trans,
- const Diag diag, const fortran_int_t n,
- const std::complex<double>* ap, std::complex<double>* x,
- const fortran_int_t incx ) {
+inline void tpmv( const Order, const UpLo, const Trans, const Diag,
+ const fortran_int_t n, const std::complex<double>* ap,
+ std::complex<double>* x, const fortran_int_t incx ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
BLAS_ZTPMV( &blas_option< UpLo >::value, &blas_option< Trans >::value,
&blas_option< Diag >::value, &n, ap, x, &incx );
Modified: sandbox/numeric_bindings/boost/numeric/bindings/blas/level2/tpsv.hpp
==============================================================================
--- sandbox/numeric_bindings/boost/numeric/bindings/blas/level2/tpsv.hpp (original)
+++ sandbox/numeric_bindings/boost/numeric/bindings/blas/level2/tpsv.hpp 2012-08-07 09:40:31 EDT (Tue, 07 Aug 2012)
@@ -65,9 +65,8 @@
// * float value-type.
//
template< typename Order, typename UpLo, typename Trans, typename Diag >
-inline void tpsv( const Order order, const UpLo uplo, const Trans trans,
- const Diag diag, const int n, const float* ap, float* x,
- const int incx ) {
+inline void tpsv( const Order, const UpLo, const Trans, const Diag,
+ const int n, const float* ap, float* x, const int incx ) {
cblas_stpsv( cblas_option< Order >::value, cblas_option< UpLo >::value,
cblas_option< Trans >::value, cblas_option< Diag >::value, n, ap,
x, incx );
@@ -79,9 +78,8 @@
// * double value-type.
//
template< typename Order, typename UpLo, typename Trans, typename Diag >
-inline void tpsv( const Order order, const UpLo uplo, const Trans trans,
- const Diag diag, const int n, const double* ap, double* x,
- const int incx ) {
+inline void tpsv( const Order, const UpLo, const Trans, const Diag,
+ const int n, const double* ap, double* x, const int incx ) {
cblas_dtpsv( cblas_option< Order >::value, cblas_option< UpLo >::value,
cblas_option< Trans >::value, cblas_option< Diag >::value, n, ap,
x, incx );
@@ -93,9 +91,9 @@
// * complex<float> value-type.
//
template< typename Order, typename UpLo, typename Trans, typename Diag >
-inline void tpsv( const Order order, const UpLo uplo, const Trans trans,
- const Diag diag, const int n, const std::complex<float>* ap,
- std::complex<float>* x, const int incx ) {
+inline void tpsv( const Order, const UpLo, const Trans, const Diag,
+ const int n, const std::complex<float>* ap, std::complex<float>* x,
+ const int incx ) {
cblas_ctpsv( cblas_option< Order >::value, cblas_option< UpLo >::value,
cblas_option< Trans >::value, cblas_option< Diag >::value, n, ap,
x, incx );
@@ -107,9 +105,9 @@
// * complex<double> value-type.
//
template< typename Order, typename UpLo, typename Trans, typename Diag >
-inline void tpsv( const Order order, const UpLo uplo, const Trans trans,
- const Diag diag, const int n, const std::complex<double>* ap,
- std::complex<double>* x, const int incx ) {
+inline void tpsv( const Order, const UpLo, const Trans, const Diag,
+ const int n, const std::complex<double>* ap, std::complex<double>* x,
+ const int incx ) {
cblas_ztpsv( cblas_option< Order >::value, cblas_option< UpLo >::value,
cblas_option< Trans >::value, cblas_option< Diag >::value, n, ap,
x, incx );
@@ -122,9 +120,8 @@
// * float value-type.
//
template< typename Order, typename UpLo, typename Trans, typename Diag >
-inline void tpsv( const Order order, const UpLo uplo, const Trans trans,
- const Diag diag, const int n, const float* ap, float* x,
- const int incx ) {
+inline void tpsv( const Order, const UpLo, const Trans, const Diag,
+ const int n, const float* ap, float* x, const int incx ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
cublasStpsv( blas_option< UpLo >::value, blas_option< Trans >::value,
blas_option< Diag >::value, n, ap, x, incx );
@@ -136,11 +133,11 @@
// * double value-type.
//
template< typename Order, typename UpLo, typename Trans, typename Diag >
-inline void tpsv( const Order order, const UpLo uplo, const Trans trans,
- const Diag diag, const int n, const double* ap, double* x,
- const int incx ) {
+inline void tpsv( const Order, const UpLo, const Trans, const Diag,
+ const int n, const double* ap, double* x, const int incx ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
- // NOT FOUND();
+ cublasDtpsv( blas_option< UpLo >::value, blas_option< Trans >::value,
+ blas_option< Diag >::value, n, ap, x, incx );
}
//
@@ -149,9 +146,9 @@
// * complex<float> value-type.
//
template< typename Order, typename UpLo, typename Trans, typename Diag >
-inline void tpsv( const Order order, const UpLo uplo, const Trans trans,
- const Diag diag, const int n, const std::complex<float>* ap,
- std::complex<float>* x, const int incx ) {
+inline void tpsv( const Order, const UpLo, const Trans, const Diag,
+ const int n, const std::complex<float>* ap, std::complex<float>* x,
+ const int incx ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
cublasCtpsv( blas_option< UpLo >::value, blas_option< Trans >::value,
blas_option< Diag >::value, n, ap, x, incx );
@@ -163,11 +160,12 @@
// * complex<double> value-type.
//
template< typename Order, typename UpLo, typename Trans, typename Diag >
-inline void tpsv( const Order order, const UpLo uplo, const Trans trans,
- const Diag diag, const int n, const std::complex<double>* ap,
- std::complex<double>* x, const int incx ) {
+inline void tpsv( const Order, const UpLo, const Trans, const Diag,
+ const int n, const std::complex<double>* ap, std::complex<double>* x,
+ const int incx ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
- // NOT FOUND();
+ cublasZtpsv( blas_option< UpLo >::value, blas_option< Trans >::value,
+ blas_option< Diag >::value, n, ap, x, incx );
}
#else
@@ -177,8 +175,8 @@
// * float value-type.
//
template< typename Order, typename UpLo, typename Trans, typename Diag >
-inline void tpsv( const Order order, const UpLo uplo, const Trans trans,
- const Diag diag, const fortran_int_t n, const float* ap, float* x,
+inline void tpsv( const Order, const UpLo, const Trans, const Diag,
+ const fortran_int_t n, const float* ap, float* x,
const fortran_int_t incx ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
BLAS_STPSV( &blas_option< UpLo >::value, &blas_option< Trans >::value,
@@ -191,8 +189,8 @@
// * double value-type.
//
template< typename Order, typename UpLo, typename Trans, typename Diag >
-inline void tpsv( const Order order, const UpLo uplo, const Trans trans,
- const Diag diag, const fortran_int_t n, const double* ap, double* x,
+inline void tpsv( const Order, const UpLo, const Trans, const Diag,
+ const fortran_int_t n, const double* ap, double* x,
const fortran_int_t incx ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
BLAS_DTPSV( &blas_option< UpLo >::value, &blas_option< Trans >::value,
@@ -205,8 +203,8 @@
// * complex<float> value-type.
//
template< typename Order, typename UpLo, typename Trans, typename Diag >
-inline void tpsv( const Order order, const UpLo uplo, const Trans trans,
- const Diag diag, const fortran_int_t n, const std::complex<float>* ap,
+inline void tpsv( const Order, const UpLo, const Trans, const Diag,
+ const fortran_int_t n, const std::complex<float>* ap,
std::complex<float>* x, const fortran_int_t incx ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
BLAS_CTPSV( &blas_option< UpLo >::value, &blas_option< Trans >::value,
@@ -219,10 +217,9 @@
// * complex<double> value-type.
//
template< typename Order, typename UpLo, typename Trans, typename Diag >
-inline void tpsv( const Order order, const UpLo uplo, const Trans trans,
- const Diag diag, const fortran_int_t n,
- const std::complex<double>* ap, std::complex<double>* x,
- const fortran_int_t incx ) {
+inline void tpsv( const Order, const UpLo, const Trans, const Diag,
+ const fortran_int_t n, const std::complex<double>* ap,
+ std::complex<double>* x, const fortran_int_t incx ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
BLAS_ZTPSV( &blas_option< UpLo >::value, &blas_option< Trans >::value,
&blas_option< Diag >::value, &n, ap, x, &incx );
Modified: sandbox/numeric_bindings/boost/numeric/bindings/blas/level2/trmv.hpp
==============================================================================
--- sandbox/numeric_bindings/boost/numeric/bindings/blas/level2/trmv.hpp (original)
+++ sandbox/numeric_bindings/boost/numeric/bindings/blas/level2/trmv.hpp 2012-08-07 09:40:31 EDT (Tue, 07 Aug 2012)
@@ -65,8 +65,8 @@
// * float value-type.
//
template< typename Order, typename UpLo, typename Trans, typename Diag >
-inline void trmv( const Order order, const UpLo uplo, const Trans trans,
- const Diag diag, const int n, const float* a, const int lda, float* x,
+inline void trmv( const Order, const UpLo, const Trans, const Diag,
+ const int n, const float* a, const int lda, float* x,
const int incx ) {
cblas_strmv( cblas_option< Order >::value, cblas_option< UpLo >::value,
cblas_option< Trans >::value, cblas_option< Diag >::value, n, a,
@@ -79,9 +79,9 @@
// * double value-type.
//
template< typename Order, typename UpLo, typename Trans, typename Diag >
-inline void trmv( const Order order, const UpLo uplo, const Trans trans,
- const Diag diag, const int n, const double* a, const int lda,
- double* x, const int incx ) {
+inline void trmv( const Order, const UpLo, const Trans, const Diag,
+ const int n, const double* a, const int lda, double* x,
+ const int incx ) {
cblas_dtrmv( cblas_option< Order >::value, cblas_option< UpLo >::value,
cblas_option< Trans >::value, cblas_option< Diag >::value, n, a,
lda, x, incx );
@@ -93,9 +93,9 @@
// * complex<float> value-type.
//
template< typename Order, typename UpLo, typename Trans, typename Diag >
-inline void trmv( const Order order, const UpLo uplo, const Trans trans,
- const Diag diag, const int n, const std::complex<float>* a,
- const int lda, std::complex<float>* x, const int incx ) {
+inline void trmv( const Order, const UpLo, const Trans, const Diag,
+ const int n, const std::complex<float>* a, const int lda,
+ std::complex<float>* x, const int incx ) {
cblas_ctrmv( cblas_option< Order >::value, cblas_option< UpLo >::value,
cblas_option< Trans >::value, cblas_option< Diag >::value, n, a,
lda, x, incx );
@@ -107,9 +107,9 @@
// * complex<double> value-type.
//
template< typename Order, typename UpLo, typename Trans, typename Diag >
-inline void trmv( const Order order, const UpLo uplo, const Trans trans,
- const Diag diag, const int n, const std::complex<double>* a,
- const int lda, std::complex<double>* x, const int incx ) {
+inline void trmv( const Order, const UpLo, const Trans, const Diag,
+ const int n, const std::complex<double>* a, const int lda,
+ std::complex<double>* x, const int incx ) {
cblas_ztrmv( cblas_option< Order >::value, cblas_option< UpLo >::value,
cblas_option< Trans >::value, cblas_option< Diag >::value, n, a,
lda, x, incx );
@@ -122,8 +122,8 @@
// * float value-type.
//
template< typename Order, typename UpLo, typename Trans, typename Diag >
-inline void trmv( const Order order, const UpLo uplo, const Trans trans,
- const Diag diag, const int n, const float* a, const int lda, float* x,
+inline void trmv( const Order, const UpLo, const Trans, const Diag,
+ const int n, const float* a, const int lda, float* x,
const int incx ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
cublasStrmv( blas_option< UpLo >::value, blas_option< Trans >::value,
@@ -136,11 +136,12 @@
// * double value-type.
//
template< typename Order, typename UpLo, typename Trans, typename Diag >
-inline void trmv( const Order order, const UpLo uplo, const Trans trans,
- const Diag diag, const int n, const double* a, const int lda,
- double* x, const int incx ) {
+inline void trmv( const Order, const UpLo, const Trans, const Diag,
+ const int n, const double* a, const int lda, double* x,
+ const int incx ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
- // NOT FOUND();
+ cublasDtrmv( blas_option< UpLo >::value, blas_option< Trans >::value,
+ blas_option< Diag >::value, n, a, lda, x, incx );
}
//
@@ -149,11 +150,12 @@
// * complex<float> value-type.
//
template< typename Order, typename UpLo, typename Trans, typename Diag >
-inline void trmv( const Order order, const UpLo uplo, const Trans trans,
- const Diag diag, const int n, const std::complex<float>* a,
- const int lda, std::complex<float>* x, const int incx ) {
+inline void trmv( const Order, const UpLo, const Trans, const Diag,
+ const int n, const std::complex<float>* a, const int lda,
+ std::complex<float>* x, const int incx ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
- // NOT FOUND();
+ cublasCtrmv( blas_option< UpLo >::value, blas_option< Trans >::value,
+ blas_option< Diag >::value, n, a, lda, x, incx );
}
//
@@ -162,11 +164,12 @@
// * complex<double> value-type.
//
template< typename Order, typename UpLo, typename Trans, typename Diag >
-inline void trmv( const Order order, const UpLo uplo, const Trans trans,
- const Diag diag, const int n, const std::complex<double>* a,
- const int lda, std::complex<double>* x, const int incx ) {
+inline void trmv( const Order, const UpLo, const Trans, const Diag,
+ const int n, const std::complex<double>* a, const int lda,
+ std::complex<double>* x, const int incx ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
- // NOT FOUND();
+ cublasZtrmv( blas_option< UpLo >::value, blas_option< Trans >::value,
+ blas_option< Diag >::value, n, a, lda, x, incx );
}
#else
@@ -176,9 +179,9 @@
// * float value-type.
//
template< typename Order, typename UpLo, typename Trans, typename Diag >
-inline void trmv( const Order order, const UpLo uplo, const Trans trans,
- const Diag diag, const fortran_int_t n, const float* a,
- const fortran_int_t lda, float* x, const fortran_int_t incx ) {
+inline void trmv( const Order, const UpLo, const Trans, const Diag,
+ const fortran_int_t n, const float* a, const fortran_int_t lda,
+ float* x, const fortran_int_t incx ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
BLAS_STRMV( &blas_option< UpLo >::value, &blas_option< Trans >::value,
&blas_option< Diag >::value, &n, a, &lda, x, &incx );
@@ -190,9 +193,9 @@
// * double value-type.
//
template< typename Order, typename UpLo, typename Trans, typename Diag >
-inline void trmv( const Order order, const UpLo uplo, const Trans trans,
- const Diag diag, const fortran_int_t n, const double* a,
- const fortran_int_t lda, double* x, const fortran_int_t incx ) {
+inline void trmv( const Order, const UpLo, const Trans, const Diag,
+ const fortran_int_t n, const double* a, const fortran_int_t lda,
+ double* x, const fortran_int_t incx ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
BLAS_DTRMV( &blas_option< UpLo >::value, &blas_option< Trans >::value,
&blas_option< Diag >::value, &n, a, &lda, x, &incx );
@@ -204,8 +207,8 @@
// * complex<float> value-type.
//
template< typename Order, typename UpLo, typename Trans, typename Diag >
-inline void trmv( const Order order, const UpLo uplo, const Trans trans,
- const Diag diag, const fortran_int_t n, const std::complex<float>* a,
+inline void trmv( const Order, const UpLo, const Trans, const Diag,
+ const fortran_int_t n, const std::complex<float>* a,
const fortran_int_t lda, std::complex<float>* x,
const fortran_int_t incx ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
@@ -219,8 +222,8 @@
// * complex<double> value-type.
//
template< typename Order, typename UpLo, typename Trans, typename Diag >
-inline void trmv( const Order order, const UpLo uplo, const Trans trans,
- const Diag diag, const fortran_int_t n, const std::complex<double>* a,
+inline void trmv( const Order, const UpLo, const Trans, const Diag,
+ const fortran_int_t n, const std::complex<double>* a,
const fortran_int_t lda, std::complex<double>* x,
const fortran_int_t incx ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
Modified: sandbox/numeric_bindings/boost/numeric/bindings/blas/level2/trsv.hpp
==============================================================================
--- sandbox/numeric_bindings/boost/numeric/bindings/blas/level2/trsv.hpp (original)
+++ sandbox/numeric_bindings/boost/numeric/bindings/blas/level2/trsv.hpp 2012-08-07 09:40:31 EDT (Tue, 07 Aug 2012)
@@ -65,8 +65,8 @@
// * float value-type.
//
template< typename Order, typename UpLo, typename Trans, typename Diag >
-inline void trsv( const Order order, const UpLo uplo, const Trans trans,
- const Diag diag, const int n, const float* a, const int lda, float* x,
+inline void trsv( const Order, const UpLo, const Trans, const Diag,
+ const int n, const float* a, const int lda, float* x,
const int incx ) {
cblas_strsv( cblas_option< Order >::value, cblas_option< UpLo >::value,
cblas_option< Trans >::value, cblas_option< Diag >::value, n, a,
@@ -79,9 +79,9 @@
// * double value-type.
//
template< typename Order, typename UpLo, typename Trans, typename Diag >
-inline void trsv( const Order order, const UpLo uplo, const Trans trans,
- const Diag diag, const int n, const double* a, const int lda,
- double* x, const int incx ) {
+inline void trsv( const Order, const UpLo, const Trans, const Diag,
+ const int n, const double* a, const int lda, double* x,
+ const int incx ) {
cblas_dtrsv( cblas_option< Order >::value, cblas_option< UpLo >::value,
cblas_option< Trans >::value, cblas_option< Diag >::value, n, a,
lda, x, incx );
@@ -93,9 +93,9 @@
// * complex<float> value-type.
//
template< typename Order, typename UpLo, typename Trans, typename Diag >
-inline void trsv( const Order order, const UpLo uplo, const Trans trans,
- const Diag diag, const int n, const std::complex<float>* a,
- const int lda, std::complex<float>* x, const int incx ) {
+inline void trsv( const Order, const UpLo, const Trans, const Diag,
+ const int n, const std::complex<float>* a, const int lda,
+ std::complex<float>* x, const int incx ) {
cblas_ctrsv( cblas_option< Order >::value, cblas_option< UpLo >::value,
cblas_option< Trans >::value, cblas_option< Diag >::value, n, a,
lda, x, incx );
@@ -107,9 +107,9 @@
// * complex<double> value-type.
//
template< typename Order, typename UpLo, typename Trans, typename Diag >
-inline void trsv( const Order order, const UpLo uplo, const Trans trans,
- const Diag diag, const int n, const std::complex<double>* a,
- const int lda, std::complex<double>* x, const int incx ) {
+inline void trsv( const Order, const UpLo, const Trans, const Diag,
+ const int n, const std::complex<double>* a, const int lda,
+ std::complex<double>* x, const int incx ) {
cblas_ztrsv( cblas_option< Order >::value, cblas_option< UpLo >::value,
cblas_option< Trans >::value, cblas_option< Diag >::value, n, a,
lda, x, incx );
@@ -122,8 +122,8 @@
// * float value-type.
//
template< typename Order, typename UpLo, typename Trans, typename Diag >
-inline void trsv( const Order order, const UpLo uplo, const Trans trans,
- const Diag diag, const int n, const float* a, const int lda, float* x,
+inline void trsv( const Order, const UpLo, const Trans, const Diag,
+ const int n, const float* a, const int lda, float* x,
const int incx ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
cublasStrsv( blas_option< UpLo >::value, blas_option< Trans >::value,
@@ -136,9 +136,9 @@
// * double value-type.
//
template< typename Order, typename UpLo, typename Trans, typename Diag >
-inline void trsv( const Order order, const UpLo uplo, const Trans trans,
- const Diag diag, const int n, const double* a, const int lda,
- double* x, const int incx ) {
+inline void trsv( const Order, const UpLo, const Trans, const Diag,
+ const int n, const double* a, const int lda, double* x,
+ const int incx ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
cublasDtrsv( blas_option< UpLo >::value, blas_option< Trans >::value,
blas_option< Diag >::value, n, a, lda, x, incx );
@@ -150,9 +150,9 @@
// * complex<float> value-type.
//
template< typename Order, typename UpLo, typename Trans, typename Diag >
-inline void trsv( const Order order, const UpLo uplo, const Trans trans,
- const Diag diag, const int n, const std::complex<float>* a,
- const int lda, std::complex<float>* x, const int incx ) {
+inline void trsv( const Order, const UpLo, const Trans, const Diag,
+ const int n, const std::complex<float>* a, const int lda,
+ std::complex<float>* x, const int incx ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
cublasCtrsv( blas_option< UpLo >::value, blas_option< Trans >::value,
blas_option< Diag >::value, n, a, lda, x, incx );
@@ -164,11 +164,12 @@
// * complex<double> value-type.
//
template< typename Order, typename UpLo, typename Trans, typename Diag >
-inline void trsv( const Order order, const UpLo uplo, const Trans trans,
- const Diag diag, const int n, const std::complex<double>* a,
- const int lda, std::complex<double>* x, const int incx ) {
+inline void trsv( const Order, const UpLo, const Trans, const Diag,
+ const int n, const std::complex<double>* a, const int lda,
+ std::complex<double>* x, const int incx ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
- // NOT FOUND();
+ cublasZtrsv( blas_option< UpLo >::value, blas_option< Trans >::value,
+ blas_option< Diag >::value, n, a, lda, x, incx );
}
#else
@@ -178,9 +179,9 @@
// * float value-type.
//
template< typename Order, typename UpLo, typename Trans, typename Diag >
-inline void trsv( const Order order, const UpLo uplo, const Trans trans,
- const Diag diag, const fortran_int_t n, const float* a,
- const fortran_int_t lda, float* x, const fortran_int_t incx ) {
+inline void trsv( const Order, const UpLo, const Trans, const Diag,
+ const fortran_int_t n, const float* a, const fortran_int_t lda,
+ float* x, const fortran_int_t incx ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
BLAS_STRSV( &blas_option< UpLo >::value, &blas_option< Trans >::value,
&blas_option< Diag >::value, &n, a, &lda, x, &incx );
@@ -192,9 +193,9 @@
// * double value-type.
//
template< typename Order, typename UpLo, typename Trans, typename Diag >
-inline void trsv( const Order order, const UpLo uplo, const Trans trans,
- const Diag diag, const fortran_int_t n, const double* a,
- const fortran_int_t lda, double* x, const fortran_int_t incx ) {
+inline void trsv( const Order, const UpLo, const Trans, const Diag,
+ const fortran_int_t n, const double* a, const fortran_int_t lda,
+ double* x, const fortran_int_t incx ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
BLAS_DTRSV( &blas_option< UpLo >::value, &blas_option< Trans >::value,
&blas_option< Diag >::value, &n, a, &lda, x, &incx );
@@ -206,8 +207,8 @@
// * complex<float> value-type.
//
template< typename Order, typename UpLo, typename Trans, typename Diag >
-inline void trsv( const Order order, const UpLo uplo, const Trans trans,
- const Diag diag, const fortran_int_t n, const std::complex<float>* a,
+inline void trsv( const Order, const UpLo, const Trans, const Diag,
+ const fortran_int_t n, const std::complex<float>* a,
const fortran_int_t lda, std::complex<float>* x,
const fortran_int_t incx ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
@@ -221,8 +222,8 @@
// * complex<double> value-type.
//
template< typename Order, typename UpLo, typename Trans, typename Diag >
-inline void trsv( const Order order, const UpLo uplo, const Trans trans,
- const Diag diag, const fortran_int_t n, const std::complex<double>* a,
+inline void trsv( const Order, const UpLo, const Trans, const Diag,
+ const fortran_int_t n, const std::complex<double>* a,
const fortran_int_t lda, std::complex<double>* x,
const fortran_int_t incx ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
Modified: sandbox/numeric_bindings/boost/numeric/bindings/blas/level3/gemm.hpp
==============================================================================
--- sandbox/numeric_bindings/boost/numeric/bindings/blas/level3/gemm.hpp (original)
+++ sandbox/numeric_bindings/boost/numeric/bindings/blas/level3/gemm.hpp 2012-08-07 09:40:31 EDT (Tue, 07 Aug 2012)
@@ -63,10 +63,10 @@
// * float value-type.
//
template< typename Order, typename TransA, typename TransB >
-inline void gemm( const Order order, const TransA transa, const TransB transb,
- const int m, const int n, const int k, const float alpha,
- const float* a, const int lda, const float* b, const int ldb,
- const float beta, float* c, const int ldc ) {
+inline void gemm( const Order, const TransA, const TransB, const int m,
+ const int n, const int k, const float alpha, const float* a,
+ const int lda, const float* b, const int ldb, const float beta,
+ float* c, const int ldc ) {
cblas_sgemm( cblas_option< Order >::value, cblas_option< TransA >::value,
cblas_option< TransB >::value, m, n, k, alpha, a, lda, b, ldb,
beta, c, ldc );
@@ -78,10 +78,10 @@
// * double value-type.
//
template< typename Order, typename TransA, typename TransB >
-inline void gemm( const Order order, const TransA transa, const TransB transb,
- const int m, const int n, const int k, const double alpha,
- const double* a, const int lda, const double* b, const int ldb,
- const double beta, double* c, const int ldc ) {
+inline void gemm( const Order, const TransA, const TransB, const int m,
+ const int n, const int k, const double alpha, const double* a,
+ const int lda, const double* b, const int ldb, const double beta,
+ double* c, const int ldc ) {
cblas_dgemm( cblas_option< Order >::value, cblas_option< TransA >::value,
cblas_option< TransB >::value, m, n, k, alpha, a, lda, b, ldb,
beta, c, ldc );
@@ -93,10 +93,10 @@
// * complex<float> value-type.
//
template< typename Order, typename TransA, typename TransB >
-inline void gemm( const Order order, const TransA transa, const TransB transb,
- const int m, const int n, const int k,
- const std::complex<float> alpha, const std::complex<float>* a,
- const int lda, const std::complex<float>* b, const int ldb,
+inline void gemm( const Order, const TransA, const TransB, const int m,
+ const int n, const int k, const std::complex<float> alpha,
+ const std::complex<float>* a, const int lda,
+ const std::complex<float>* b, const int ldb,
const std::complex<float> beta, std::complex<float>* c,
const int ldc ) {
cblas_cgemm( cblas_option< Order >::value, cblas_option< TransA >::value,
@@ -110,10 +110,10 @@
// * complex<double> value-type.
//
template< typename Order, typename TransA, typename TransB >
-inline void gemm( const Order order, const TransA transa, const TransB transb,
- const int m, const int n, const int k,
- const std::complex<double> alpha, const std::complex<double>* a,
- const int lda, const std::complex<double>* b, const int ldb,
+inline void gemm( const Order, const TransA, const TransB, const int m,
+ const int n, const int k, const std::complex<double> alpha,
+ const std::complex<double>* a, const int lda,
+ const std::complex<double>* b, const int ldb,
const std::complex<double> beta, std::complex<double>* c,
const int ldc ) {
cblas_zgemm( cblas_option< Order >::value, cblas_option< TransA >::value,
@@ -128,10 +128,10 @@
// * float value-type.
//
template< typename Order, typename TransA, typename TransB >
-inline void gemm( const Order order, const TransA transa, const TransB transb,
- const int m, const int n, const int k, const float alpha,
- const float* a, const int lda, const float* b, const int ldb,
- const float beta, float* c, const int ldc ) {
+inline void gemm( const Order, const TransA, const TransB, const int m,
+ const int n, const int k, const float alpha, const float* a,
+ const int lda, const float* b, const int ldb, const float beta,
+ float* c, const int ldc ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
cublasSgemm( blas_option< TransA >::value, blas_option< TransB >::value,
m, n, k, alpha, a, lda, b, ldb, beta, c, ldc );
@@ -143,10 +143,10 @@
// * double value-type.
//
template< typename Order, typename TransA, typename TransB >
-inline void gemm( const Order order, const TransA transa, const TransB transb,
- const int m, const int n, const int k, const double alpha,
- const double* a, const int lda, const double* b, const int ldb,
- const double beta, double* c, const int ldc ) {
+inline void gemm( const Order, const TransA, const TransB, const int m,
+ const int n, const int k, const double alpha, const double* a,
+ const int lda, const double* b, const int ldb, const double beta,
+ double* c, const int ldc ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
cublasDgemm( blas_option< TransA >::value, blas_option< TransB >::value,
m, n, k, alpha, a, lda, b, ldb, beta, c, ldc );
@@ -158,10 +158,10 @@
// * complex<float> value-type.
//
template< typename Order, typename TransA, typename TransB >
-inline void gemm( const Order order, const TransA transa, const TransB transb,
- const int m, const int n, const int k,
- const std::complex<float> alpha, const std::complex<float>* a,
- const int lda, const std::complex<float>* b, const int ldb,
+inline void gemm( const Order, const TransA, const TransB, const int m,
+ const int n, const int k, const std::complex<float> alpha,
+ const std::complex<float>* a, const int lda,
+ const std::complex<float>* b, const int ldb,
const std::complex<float> beta, std::complex<float>* c,
const int ldc ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
@@ -175,10 +175,10 @@
// * complex<double> value-type.
//
template< typename Order, typename TransA, typename TransB >
-inline void gemm( const Order order, const TransA transa, const TransB transb,
- const int m, const int n, const int k,
- const std::complex<double> alpha, const std::complex<double>* a,
- const int lda, const std::complex<double>* b, const int ldb,
+inline void gemm( const Order, const TransA, const TransB, const int m,
+ const int n, const int k, const std::complex<double> alpha,
+ const std::complex<double>* a, const int lda,
+ const std::complex<double>* b, const int ldb,
const std::complex<double> beta, std::complex<double>* c,
const int ldc ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
@@ -193,7 +193,7 @@
// * float value-type.
//
template< typename Order, typename TransA, typename TransB >
-inline void gemm( const Order order, const TransA transa, const TransB transb,
+inline void gemm( const Order, const TransA, const TransB,
const fortran_int_t m, const fortran_int_t n, const fortran_int_t k,
const float alpha, const float* a, const fortran_int_t lda,
const float* b, const fortran_int_t ldb, const float beta, float* c,
@@ -209,7 +209,7 @@
// * double value-type.
//
template< typename Order, typename TransA, typename TransB >
-inline void gemm( const Order order, const TransA transa, const TransB transb,
+inline void gemm( const Order, const TransA, const TransB,
const fortran_int_t m, const fortran_int_t n, const fortran_int_t k,
const double alpha, const double* a, const fortran_int_t lda,
const double* b, const fortran_int_t ldb, const double beta,
@@ -225,7 +225,7 @@
// * complex<float> value-type.
//
template< typename Order, typename TransA, typename TransB >
-inline void gemm( const Order order, const TransA transa, const TransB transb,
+inline void gemm( const Order, const TransA, const TransB,
const fortran_int_t m, const fortran_int_t n, const fortran_int_t k,
const std::complex<float> alpha, const std::complex<float>* a,
const fortran_int_t lda, const std::complex<float>* b,
@@ -242,7 +242,7 @@
// * complex<double> value-type.
//
template< typename Order, typename TransA, typename TransB >
-inline void gemm( const Order order, const TransA transa, const TransB transb,
+inline void gemm( const Order, const TransA, const TransB,
const fortran_int_t m, const fortran_int_t n, const fortran_int_t k,
const std::complex<double> alpha, const std::complex<double>* a,
const fortran_int_t lda, const std::complex<double>* b,
Modified: sandbox/numeric_bindings/boost/numeric/bindings/blas/level3/hemm.hpp
==============================================================================
--- sandbox/numeric_bindings/boost/numeric/bindings/blas/level3/hemm.hpp (original)
+++ sandbox/numeric_bindings/boost/numeric/bindings/blas/level3/hemm.hpp 2012-08-07 09:40:31 EDT (Tue, 07 Aug 2012)
@@ -63,10 +63,10 @@
// * float value-type.
//
template< typename Order, typename Side, typename UpLo >
-inline void hemm( const Order order, const Side side, const UpLo uplo,
- const int m, const int n, const float alpha, const float* a,
- const int lda, const float* b, const int ldb, const float beta,
- float* c, const int ldc ) {
+inline void hemm( const Order, const Side, const UpLo, const int m,
+ const int n, const float alpha, const float* a, const int lda,
+ const float* b, const int ldb, const float beta, float* c,
+ const int ldc ) {
cblas_ssymm( cblas_option< Order >::value, cblas_option< Side >::value,
cblas_option< UpLo >::value, m, n, alpha, a, lda, b, ldb, beta, c,
ldc );
@@ -78,10 +78,10 @@
// * double value-type.
//
template< typename Order, typename Side, typename UpLo >
-inline void hemm( const Order order, const Side side, const UpLo uplo,
- const int m, const int n, const double alpha, const double* a,
- const int lda, const double* b, const int ldb, const double beta,
- double* c, const int ldc ) {
+inline void hemm( const Order, const Side, const UpLo, const int m,
+ const int n, const double alpha, const double* a, const int lda,
+ const double* b, const int ldb, const double beta, double* c,
+ const int ldc ) {
cblas_dsymm( cblas_option< Order >::value, cblas_option< Side >::value,
cblas_option< UpLo >::value, m, n, alpha, a, lda, b, ldb, beta, c,
ldc );
@@ -93,8 +93,8 @@
// * complex<float> value-type.
//
template< typename Order, typename Side, typename UpLo >
-inline void hemm( const Order order, const Side side, const UpLo uplo,
- const int m, const int n, const std::complex<float> alpha,
+inline void hemm( const Order, const Side, const UpLo, const int m,
+ const int n, const std::complex<float> alpha,
const std::complex<float>* a, const int lda,
const std::complex<float>* b, const int ldb,
const std::complex<float> beta, std::complex<float>* c,
@@ -110,8 +110,8 @@
// * complex<double> value-type.
//
template< typename Order, typename Side, typename UpLo >
-inline void hemm( const Order order, const Side side, const UpLo uplo,
- const int m, const int n, const std::complex<double> alpha,
+inline void hemm( const Order, const Side, const UpLo, const int m,
+ const int n, const std::complex<double> alpha,
const std::complex<double>* a, const int lda,
const std::complex<double>* b, const int ldb,
const std::complex<double> beta, std::complex<double>* c,
@@ -128,10 +128,10 @@
// * float value-type.
//
template< typename Order, typename Side, typename UpLo >
-inline void hemm( const Order order, const Side side, const UpLo uplo,
- const int m, const int n, const float alpha, const float* a,
- const int lda, const float* b, const int ldb, const float beta,
- float* c, const int ldc ) {
+inline void hemm( const Order, const Side, const UpLo, const int m,
+ const int n, const float alpha, const float* a, const int lda,
+ const float* b, const int ldb, const float beta, float* c,
+ const int ldc ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
cublasSsymm( blas_option< Side >::value, blas_option< UpLo >::value, m, n,
alpha, a, lda, b, ldb, beta, c, ldc );
@@ -143,10 +143,10 @@
// * double value-type.
//
template< typename Order, typename Side, typename UpLo >
-inline void hemm( const Order order, const Side side, const UpLo uplo,
- const int m, const int n, const double alpha, const double* a,
- const int lda, const double* b, const int ldb, const double beta,
- double* c, const int ldc ) {
+inline void hemm( const Order, const Side, const UpLo, const int m,
+ const int n, const double alpha, const double* a, const int lda,
+ const double* b, const int ldb, const double beta, double* c,
+ const int ldc ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
cublasDsymm( blas_option< Side >::value, blas_option< UpLo >::value, m, n,
alpha, a, lda, b, ldb, beta, c, ldc );
@@ -158,8 +158,8 @@
// * complex<float> value-type.
//
template< typename Order, typename Side, typename UpLo >
-inline void hemm( const Order order, const Side side, const UpLo uplo,
- const int m, const int n, const std::complex<float> alpha,
+inline void hemm( const Order, const Side, const UpLo, const int m,
+ const int n, const std::complex<float> alpha,
const std::complex<float>* a, const int lda,
const std::complex<float>* b, const int ldb,
const std::complex<float> beta, std::complex<float>* c,
@@ -175,14 +175,15 @@
// * complex<double> value-type.
//
template< typename Order, typename Side, typename UpLo >
-inline void hemm( const Order order, const Side side, const UpLo uplo,
- const int m, const int n, const std::complex<double> alpha,
+inline void hemm( const Order, const Side, const UpLo, const int m,
+ const int n, const std::complex<double> alpha,
const std::complex<double>* a, const int lda,
const std::complex<double>* b, const int ldb,
const std::complex<double> beta, std::complex<double>* c,
const int ldc ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
- // NOT FOUND();
+ cublasZhemm( blas_option< Side >::value, blas_option< UpLo >::value, m, n,
+ alpha, a, lda, b, ldb, beta, c, ldc );
}
#else
@@ -192,11 +193,10 @@
// * float value-type.
//
template< typename Order, typename Side, typename UpLo >
-inline void hemm( const Order order, const Side side, const UpLo uplo,
- const fortran_int_t m, const fortran_int_t n, const float alpha,
- const float* a, const fortran_int_t lda, const float* b,
- const fortran_int_t ldb, const float beta, float* c,
- const fortran_int_t ldc ) {
+inline void hemm( const Order, const Side, const UpLo, const fortran_int_t m,
+ const fortran_int_t n, const float alpha, const float* a,
+ const fortran_int_t lda, const float* b, const fortran_int_t ldb,
+ const float beta, float* c, const fortran_int_t ldc ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
BLAS_SSYMM( &blas_option< Side >::value, &blas_option< UpLo >::value, &m,
&n, &alpha, a, &lda, b, &ldb, &beta, c, &ldc );
@@ -208,11 +208,10 @@
// * double value-type.
//
template< typename Order, typename Side, typename UpLo >
-inline void hemm( const Order order, const Side side, const UpLo uplo,
- const fortran_int_t m, const fortran_int_t n, const double alpha,
- const double* a, const fortran_int_t lda, const double* b,
- const fortran_int_t ldb, const double beta, double* c,
- const fortran_int_t ldc ) {
+inline void hemm( const Order, const Side, const UpLo, const fortran_int_t m,
+ const fortran_int_t n, const double alpha, const double* a,
+ const fortran_int_t lda, const double* b, const fortran_int_t ldb,
+ const double beta, double* c, const fortran_int_t ldc ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
BLAS_DSYMM( &blas_option< Side >::value, &blas_option< UpLo >::value, &m,
&n, &alpha, a, &lda, b, &ldb, &beta, c, &ldc );
@@ -224,12 +223,12 @@
// * complex<float> value-type.
//
template< typename Order, typename Side, typename UpLo >
-inline void hemm( const Order order, const Side side, const UpLo uplo,
- const fortran_int_t m, const fortran_int_t n,
- const std::complex<float> alpha, const std::complex<float>* a,
- const fortran_int_t lda, const std::complex<float>* b,
- const fortran_int_t ldb, const std::complex<float> beta,
- std::complex<float>* c, const fortran_int_t ldc ) {
+inline void hemm( const Order, const Side, const UpLo, const fortran_int_t m,
+ const fortran_int_t n, const std::complex<float> alpha,
+ const std::complex<float>* a, const fortran_int_t lda,
+ const std::complex<float>* b, const fortran_int_t ldb,
+ const std::complex<float> beta, std::complex<float>* c,
+ const fortran_int_t ldc ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
BLAS_CHEMM( &blas_option< Side >::value, &blas_option< UpLo >::value, &m,
&n, &alpha, a, &lda, b, &ldb, &beta, c, &ldc );
@@ -241,12 +240,12 @@
// * complex<double> value-type.
//
template< typename Order, typename Side, typename UpLo >
-inline void hemm( const Order order, const Side side, const UpLo uplo,
- const fortran_int_t m, const fortran_int_t n,
- const std::complex<double> alpha, const std::complex<double>* a,
- const fortran_int_t lda, const std::complex<double>* b,
- const fortran_int_t ldb, const std::complex<double> beta,
- std::complex<double>* c, const fortran_int_t ldc ) {
+inline void hemm( const Order, const Side, const UpLo, const fortran_int_t m,
+ const fortran_int_t n, const std::complex<double> alpha,
+ const std::complex<double>* a, const fortran_int_t lda,
+ const std::complex<double>* b, const fortran_int_t ldb,
+ const std::complex<double> beta, std::complex<double>* c,
+ const fortran_int_t ldc ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
BLAS_ZHEMM( &blas_option< Side >::value, &blas_option< UpLo >::value, &m,
&n, &alpha, a, &lda, b, &ldb, &beta, c, &ldc );
Modified: sandbox/numeric_bindings/boost/numeric/bindings/blas/level3/her2k.hpp
==============================================================================
--- sandbox/numeric_bindings/boost/numeric/bindings/blas/level3/her2k.hpp (original)
+++ sandbox/numeric_bindings/boost/numeric/bindings/blas/level3/her2k.hpp 2012-08-07 09:40:31 EDT (Tue, 07 Aug 2012)
@@ -64,10 +64,10 @@
// * float value-type.
//
template< typename Order, typename UpLo, typename Trans >
-inline void her2k( const Order order, const UpLo uplo, const Trans trans,
- const int n, const int k, const float alpha, const float* a,
- const int lda, const float* b, const int ldb, const float beta,
- float* c, const int ldc ) {
+inline void her2k( const Order, const UpLo, const Trans, const int n,
+ const int k, const float alpha, const float* a, const int lda,
+ const float* b, const int ldb, const float beta, float* c,
+ const int ldc ) {
cblas_ssyr2k( cblas_option< Order >::value, cblas_option< UpLo >::value,
cblas_option< Trans >::value, n, k, alpha, a, lda, b, ldb, beta,
c, ldc );
@@ -79,10 +79,10 @@
// * double value-type.
//
template< typename Order, typename UpLo, typename Trans >
-inline void her2k( const Order order, const UpLo uplo, const Trans trans,
- const int n, const int k, const double alpha, const double* a,
- const int lda, const double* b, const int ldb, const double beta,
- double* c, const int ldc ) {
+inline void her2k( const Order, const UpLo, const Trans, const int n,
+ const int k, const double alpha, const double* a, const int lda,
+ const double* b, const int ldb, const double beta, double* c,
+ const int ldc ) {
cblas_dsyr2k( cblas_option< Order >::value, cblas_option< UpLo >::value,
cblas_option< Trans >::value, n, k, alpha, a, lda, b, ldb, beta,
c, ldc );
@@ -94,8 +94,8 @@
// * complex<float> value-type.
//
template< typename Order, typename UpLo, typename Trans >
-inline void her2k( const Order order, const UpLo uplo, const Trans trans,
- const int n, const int k, const std::complex<float> alpha,
+inline void her2k( const Order, const UpLo, const Trans, const int n,
+ const int k, const std::complex<float> alpha,
const std::complex<float>* a, const int lda,
const std::complex<float>* b, const int ldb, const float beta,
std::complex<float>* c, const int ldc ) {
@@ -110,8 +110,8 @@
// * complex<double> value-type.
//
template< typename Order, typename UpLo, typename Trans >
-inline void her2k( const Order order, const UpLo uplo, const Trans trans,
- const int n, const int k, const std::complex<double> alpha,
+inline void her2k( const Order, const UpLo, const Trans, const int n,
+ const int k, const std::complex<double> alpha,
const std::complex<double>* a, const int lda,
const std::complex<double>* b, const int ldb, const double beta,
std::complex<double>* c, const int ldc ) {
@@ -127,10 +127,10 @@
// * float value-type.
//
template< typename Order, typename UpLo, typename Trans >
-inline void her2k( const Order order, const UpLo uplo, const Trans trans,
- const int n, const int k, const float alpha, const float* a,
- const int lda, const float* b, const int ldb, const float beta,
- float* c, const int ldc ) {
+inline void her2k( const Order, const UpLo, const Trans, const int n,
+ const int k, const float alpha, const float* a, const int lda,
+ const float* b, const int ldb, const float beta, float* c,
+ const int ldc ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
cublasSsyr2k( blas_option< UpLo >::value, blas_option< Trans >::value, n,
k, alpha, a, lda, b, ldb, beta, c, ldc );
@@ -142,10 +142,10 @@
// * double value-type.
//
template< typename Order, typename UpLo, typename Trans >
-inline void her2k( const Order order, const UpLo uplo, const Trans trans,
- const int n, const int k, const double alpha, const double* a,
- const int lda, const double* b, const int ldb, const double beta,
- double* c, const int ldc ) {
+inline void her2k( const Order, const UpLo, const Trans, const int n,
+ const int k, const double alpha, const double* a, const int lda,
+ const double* b, const int ldb, const double beta, double* c,
+ const int ldc ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
cublasDsyr2k( blas_option< UpLo >::value, blas_option< Trans >::value, n,
k, alpha, a, lda, b, ldb, beta, c, ldc );
@@ -157,8 +157,8 @@
// * complex<float> value-type.
//
template< typename Order, typename UpLo, typename Trans >
-inline void her2k( const Order order, const UpLo uplo, const Trans trans,
- const int n, const int k, const std::complex<float> alpha,
+inline void her2k( const Order, const UpLo, const Trans, const int n,
+ const int k, const std::complex<float> alpha,
const std::complex<float>* a, const int lda,
const std::complex<float>* b, const int ldb, const float beta,
std::complex<float>* c, const int ldc ) {
@@ -173,13 +173,14 @@
// * complex<double> value-type.
//
template< typename Order, typename UpLo, typename Trans >
-inline void her2k( const Order order, const UpLo uplo, const Trans trans,
- const int n, const int k, const std::complex<double> alpha,
+inline void her2k( const Order, const UpLo, const Trans, const int n,
+ const int k, const std::complex<double> alpha,
const std::complex<double>* a, const int lda,
const std::complex<double>* b, const int ldb, const double beta,
std::complex<double>* c, const int ldc ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
- // NOT FOUND();
+ cublasZher2k( blas_option< UpLo >::value, blas_option< Trans >::value, n,
+ k, alpha, a, lda, b, ldb, beta, c, ldc );
}
#else
@@ -189,7 +190,7 @@
// * float value-type.
//
template< typename Order, typename UpLo, typename Trans >
-inline void her2k( const Order order, const UpLo uplo, const Trans trans,
+inline void her2k( const Order, const UpLo, const Trans,
const fortran_int_t n, const fortran_int_t k, const float alpha,
const float* a, const fortran_int_t lda, const float* b,
const fortran_int_t ldb, const float beta, float* c,
@@ -205,7 +206,7 @@
// * double value-type.
//
template< typename Order, typename UpLo, typename Trans >
-inline void her2k( const Order order, const UpLo uplo, const Trans trans,
+inline void her2k( const Order, const UpLo, const Trans,
const fortran_int_t n, const fortran_int_t k, const double alpha,
const double* a, const fortran_int_t lda, const double* b,
const fortran_int_t ldb, const double beta, double* c,
@@ -221,7 +222,7 @@
// * complex<float> value-type.
//
template< typename Order, typename UpLo, typename Trans >
-inline void her2k( const Order order, const UpLo uplo, const Trans trans,
+inline void her2k( const Order, const UpLo, const Trans,
const fortran_int_t n, const fortran_int_t k,
const std::complex<float> alpha, const std::complex<float>* a,
const fortran_int_t lda, const std::complex<float>* b,
@@ -238,7 +239,7 @@
// * complex<double> value-type.
//
template< typename Order, typename UpLo, typename Trans >
-inline void her2k( const Order order, const UpLo uplo, const Trans trans,
+inline void her2k( const Order, const UpLo, const Trans,
const fortran_int_t n, const fortran_int_t k,
const std::complex<double> alpha, const std::complex<double>* a,
const fortran_int_t lda, const std::complex<double>* b,
Modified: sandbox/numeric_bindings/boost/numeric/bindings/blas/level3/herk.hpp
==============================================================================
--- sandbox/numeric_bindings/boost/numeric/bindings/blas/level3/herk.hpp (original)
+++ sandbox/numeric_bindings/boost/numeric/bindings/blas/level3/herk.hpp 2012-08-07 09:40:31 EDT (Tue, 07 Aug 2012)
@@ -64,9 +64,9 @@
// * float value-type.
//
template< typename Order, typename UpLo, typename Trans >
-inline void herk( const Order order, const UpLo uplo, const Trans trans,
- const int n, const int k, const float alpha, const float* a,
- const int lda, const float beta, float* c, const int ldc ) {
+inline void herk( const Order, const UpLo, const Trans, const int n,
+ const int k, const float alpha, const float* a, const int lda,
+ const float beta, float* c, const int ldc ) {
cblas_ssyrk( cblas_option< Order >::value, cblas_option< UpLo >::value,
cblas_option< Trans >::value, n, k, alpha, a, lda, beta, c, ldc );
}
@@ -77,9 +77,9 @@
// * double value-type.
//
template< typename Order, typename UpLo, typename Trans >
-inline void herk( const Order order, const UpLo uplo, const Trans trans,
- const int n, const int k, const double alpha, const double* a,
- const int lda, const double beta, double* c, const int ldc ) {
+inline void herk( const Order, const UpLo, const Trans, const int n,
+ const int k, const double alpha, const double* a, const int lda,
+ const double beta, double* c, const int ldc ) {
cblas_dsyrk( cblas_option< Order >::value, cblas_option< UpLo >::value,
cblas_option< Trans >::value, n, k, alpha, a, lda, beta, c, ldc );
}
@@ -90,10 +90,10 @@
// * complex<float> value-type.
//
template< typename Order, typename UpLo, typename Trans >
-inline void herk( const Order order, const UpLo uplo, const Trans trans,
- const int n, const int k, const float alpha,
- const std::complex<float>* a, const int lda, const float beta,
- std::complex<float>* c, const int ldc ) {
+inline void herk( const Order, const UpLo, const Trans, const int n,
+ const int k, const float alpha, const std::complex<float>* a,
+ const int lda, const float beta, std::complex<float>* c,
+ const int ldc ) {
cblas_cherk( cblas_option< Order >::value, cblas_option< UpLo >::value,
cblas_option< Trans >::value, n, k, alpha, a, lda, beta, c, ldc );
}
@@ -104,10 +104,10 @@
// * complex<double> value-type.
//
template< typename Order, typename UpLo, typename Trans >
-inline void herk( const Order order, const UpLo uplo, const Trans trans,
- const int n, const int k, const double alpha,
- const std::complex<double>* a, const int lda, const double beta,
- std::complex<double>* c, const int ldc ) {
+inline void herk( const Order, const UpLo, const Trans, const int n,
+ const int k, const double alpha, const std::complex<double>* a,
+ const int lda, const double beta, std::complex<double>* c,
+ const int ldc ) {
cblas_zherk( cblas_option< Order >::value, cblas_option< UpLo >::value,
cblas_option< Trans >::value, n, k, alpha, a, lda, beta, c, ldc );
}
@@ -119,9 +119,9 @@
// * float value-type.
//
template< typename Order, typename UpLo, typename Trans >
-inline void herk( const Order order, const UpLo uplo, const Trans trans,
- const int n, const int k, const float alpha, const float* a,
- const int lda, const float beta, float* c, const int ldc ) {
+inline void herk( const Order, const UpLo, const Trans, const int n,
+ const int k, const float alpha, const float* a, const int lda,
+ const float beta, float* c, const int ldc ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
cublasSsyrk( blas_option< UpLo >::value, blas_option< Trans >::value, n,
k, alpha, a, lda, beta, c, ldc );
@@ -133,9 +133,9 @@
// * double value-type.
//
template< typename Order, typename UpLo, typename Trans >
-inline void herk( const Order order, const UpLo uplo, const Trans trans,
- const int n, const int k, const double alpha, const double* a,
- const int lda, const double beta, double* c, const int ldc ) {
+inline void herk( const Order, const UpLo, const Trans, const int n,
+ const int k, const double alpha, const double* a, const int lda,
+ const double beta, double* c, const int ldc ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
cublasDsyrk( blas_option< UpLo >::value, blas_option< Trans >::value, n,
k, alpha, a, lda, beta, c, ldc );
@@ -147,10 +147,10 @@
// * complex<float> value-type.
//
template< typename Order, typename UpLo, typename Trans >
-inline void herk( const Order order, const UpLo uplo, const Trans trans,
- const int n, const int k, const float alpha,
- const std::complex<float>* a, const int lda, const float beta,
- std::complex<float>* c, const int ldc ) {
+inline void herk( const Order, const UpLo, const Trans, const int n,
+ const int k, const float alpha, const std::complex<float>* a,
+ const int lda, const float beta, std::complex<float>* c,
+ const int ldc ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
cublasCherk( blas_option< UpLo >::value, blas_option< Trans >::value, n,
k, alpha, a, lda, beta, c, ldc );
@@ -162,12 +162,13 @@
// * complex<double> value-type.
//
template< typename Order, typename UpLo, typename Trans >
-inline void herk( const Order order, const UpLo uplo, const Trans trans,
- const int n, const int k, const double alpha,
- const std::complex<double>* a, const int lda, const double beta,
- std::complex<double>* c, const int ldc ) {
+inline void herk( const Order, const UpLo, const Trans, const int n,
+ const int k, const double alpha, const std::complex<double>* a,
+ const int lda, const double beta, std::complex<double>* c,
+ const int ldc ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
- // NOT FOUND();
+ cublasZherk( blas_option< UpLo >::value, blas_option< Trans >::value, n,
+ k, alpha, a, lda, beta, c, ldc );
}
#else
@@ -177,9 +178,9 @@
// * float value-type.
//
template< typename Order, typename UpLo, typename Trans >
-inline void herk( const Order order, const UpLo uplo, const Trans trans,
- const fortran_int_t n, const fortran_int_t k, const float alpha,
- const float* a, const fortran_int_t lda, const float beta, float* c,
+inline void herk( const Order, const UpLo, const Trans, const fortran_int_t n,
+ const fortran_int_t k, const float alpha, const float* a,
+ const fortran_int_t lda, const float beta, float* c,
const fortran_int_t ldc ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
BLAS_SSYRK( &blas_option< UpLo >::value, &blas_option< Trans >::value, &n,
@@ -192,10 +193,10 @@
// * double value-type.
//
template< typename Order, typename UpLo, typename Trans >
-inline void herk( const Order order, const UpLo uplo, const Trans trans,
- const fortran_int_t n, const fortran_int_t k, const double alpha,
- const double* a, const fortran_int_t lda, const double beta,
- double* c, const fortran_int_t ldc ) {
+inline void herk( const Order, const UpLo, const Trans, const fortran_int_t n,
+ const fortran_int_t k, const double alpha, const double* a,
+ const fortran_int_t lda, const double beta, double* c,
+ const fortran_int_t ldc ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
BLAS_DSYRK( &blas_option< UpLo >::value, &blas_option< Trans >::value, &n,
&k, &alpha, a, &lda, &beta, c, &ldc );
@@ -207,8 +208,8 @@
// * complex<float> value-type.
//
template< typename Order, typename UpLo, typename Trans >
-inline void herk( const Order order, const UpLo uplo, const Trans trans,
- const fortran_int_t n, const fortran_int_t k, const float alpha,
+inline void herk( const Order, const UpLo, const Trans, const fortran_int_t n,
+ const fortran_int_t k, const float alpha,
const std::complex<float>* a, const fortran_int_t lda,
const float beta, std::complex<float>* c, const fortran_int_t ldc ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
@@ -222,8 +223,8 @@
// * complex<double> value-type.
//
template< typename Order, typename UpLo, typename Trans >
-inline void herk( const Order order, const UpLo uplo, const Trans trans,
- const fortran_int_t n, const fortran_int_t k, const double alpha,
+inline void herk( const Order, const UpLo, const Trans, const fortran_int_t n,
+ const fortran_int_t k, const double alpha,
const std::complex<double>* a, const fortran_int_t lda,
const double beta, std::complex<double>* c, const fortran_int_t ldc ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
Modified: sandbox/numeric_bindings/boost/numeric/bindings/blas/level3/symm.hpp
==============================================================================
--- sandbox/numeric_bindings/boost/numeric/bindings/blas/level3/symm.hpp (original)
+++ sandbox/numeric_bindings/boost/numeric/bindings/blas/level3/symm.hpp 2012-08-07 09:40:31 EDT (Tue, 07 Aug 2012)
@@ -63,10 +63,10 @@
// * float value-type.
//
template< typename Order, typename Side, typename UpLo >
-inline void symm( const Order order, const Side side, const UpLo uplo,
- const int m, const int n, const float alpha, const float* a,
- const int lda, const float* b, const int ldb, const float beta,
- float* c, const int ldc ) {
+inline void symm( const Order, const Side, const UpLo, const int m,
+ const int n, const float alpha, const float* a, const int lda,
+ const float* b, const int ldb, const float beta, float* c,
+ const int ldc ) {
cblas_ssymm( cblas_option< Order >::value, cblas_option< Side >::value,
cblas_option< UpLo >::value, m, n, alpha, a, lda, b, ldb, beta, c,
ldc );
@@ -78,10 +78,10 @@
// * double value-type.
//
template< typename Order, typename Side, typename UpLo >
-inline void symm( const Order order, const Side side, const UpLo uplo,
- const int m, const int n, const double alpha, const double* a,
- const int lda, const double* b, const int ldb, const double beta,
- double* c, const int ldc ) {
+inline void symm( const Order, const Side, const UpLo, const int m,
+ const int n, const double alpha, const double* a, const int lda,
+ const double* b, const int ldb, const double beta, double* c,
+ const int ldc ) {
cblas_dsymm( cblas_option< Order >::value, cblas_option< Side >::value,
cblas_option< UpLo >::value, m, n, alpha, a, lda, b, ldb, beta, c,
ldc );
@@ -93,8 +93,8 @@
// * complex<float> value-type.
//
template< typename Order, typename Side, typename UpLo >
-inline void symm( const Order order, const Side side, const UpLo uplo,
- const int m, const int n, const std::complex<float> alpha,
+inline void symm( const Order, const Side, const UpLo, const int m,
+ const int n, const std::complex<float> alpha,
const std::complex<float>* a, const int lda,
const std::complex<float>* b, const int ldb,
const std::complex<float> beta, std::complex<float>* c,
@@ -110,8 +110,8 @@
// * complex<double> value-type.
//
template< typename Order, typename Side, typename UpLo >
-inline void symm( const Order order, const Side side, const UpLo uplo,
- const int m, const int n, const std::complex<double> alpha,
+inline void symm( const Order, const Side, const UpLo, const int m,
+ const int n, const std::complex<double> alpha,
const std::complex<double>* a, const int lda,
const std::complex<double>* b, const int ldb,
const std::complex<double> beta, std::complex<double>* c,
@@ -128,10 +128,10 @@
// * float value-type.
//
template< typename Order, typename Side, typename UpLo >
-inline void symm( const Order order, const Side side, const UpLo uplo,
- const int m, const int n, const float alpha, const float* a,
- const int lda, const float* b, const int ldb, const float beta,
- float* c, const int ldc ) {
+inline void symm( const Order, const Side, const UpLo, const int m,
+ const int n, const float alpha, const float* a, const int lda,
+ const float* b, const int ldb, const float beta, float* c,
+ const int ldc ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
cublasSsymm( blas_option< Side >::value, blas_option< UpLo >::value, m, n,
alpha, a, lda, b, ldb, beta, c, ldc );
@@ -143,10 +143,10 @@
// * double value-type.
//
template< typename Order, typename Side, typename UpLo >
-inline void symm( const Order order, const Side side, const UpLo uplo,
- const int m, const int n, const double alpha, const double* a,
- const int lda, const double* b, const int ldb, const double beta,
- double* c, const int ldc ) {
+inline void symm( const Order, const Side, const UpLo, const int m,
+ const int n, const double alpha, const double* a, const int lda,
+ const double* b, const int ldb, const double beta, double* c,
+ const int ldc ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
cublasDsymm( blas_option< Side >::value, blas_option< UpLo >::value, m, n,
alpha, a, lda, b, ldb, beta, c, ldc );
@@ -158,8 +158,8 @@
// * complex<float> value-type.
//
template< typename Order, typename Side, typename UpLo >
-inline void symm( const Order order, const Side side, const UpLo uplo,
- const int m, const int n, const std::complex<float> alpha,
+inline void symm( const Order, const Side, const UpLo, const int m,
+ const int n, const std::complex<float> alpha,
const std::complex<float>* a, const int lda,
const std::complex<float>* b, const int ldb,
const std::complex<float> beta, std::complex<float>* c,
@@ -175,14 +175,15 @@
// * complex<double> value-type.
//
template< typename Order, typename Side, typename UpLo >
-inline void symm( const Order order, const Side side, const UpLo uplo,
- const int m, const int n, const std::complex<double> alpha,
+inline void symm( const Order, const Side, const UpLo, const int m,
+ const int n, const std::complex<double> alpha,
const std::complex<double>* a, const int lda,
const std::complex<double>* b, const int ldb,
const std::complex<double> beta, std::complex<double>* c,
const int ldc ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
- // NOT FOUND();
+ cublasZsymm( blas_option< Side >::value, blas_option< UpLo >::value, m, n,
+ alpha, a, lda, b, ldb, beta, c, ldc );
}
#else
@@ -192,11 +193,10 @@
// * float value-type.
//
template< typename Order, typename Side, typename UpLo >
-inline void symm( const Order order, const Side side, const UpLo uplo,
- const fortran_int_t m, const fortran_int_t n, const float alpha,
- const float* a, const fortran_int_t lda, const float* b,
- const fortran_int_t ldb, const float beta, float* c,
- const fortran_int_t ldc ) {
+inline void symm( const Order, const Side, const UpLo, const fortran_int_t m,
+ const fortran_int_t n, const float alpha, const float* a,
+ const fortran_int_t lda, const float* b, const fortran_int_t ldb,
+ const float beta, float* c, const fortran_int_t ldc ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
BLAS_SSYMM( &blas_option< Side >::value, &blas_option< UpLo >::value, &m,
&n, &alpha, a, &lda, b, &ldb, &beta, c, &ldc );
@@ -208,11 +208,10 @@
// * double value-type.
//
template< typename Order, typename Side, typename UpLo >
-inline void symm( const Order order, const Side side, const UpLo uplo,
- const fortran_int_t m, const fortran_int_t n, const double alpha,
- const double* a, const fortran_int_t lda, const double* b,
- const fortran_int_t ldb, const double beta, double* c,
- const fortran_int_t ldc ) {
+inline void symm( const Order, const Side, const UpLo, const fortran_int_t m,
+ const fortran_int_t n, const double alpha, const double* a,
+ const fortran_int_t lda, const double* b, const fortran_int_t ldb,
+ const double beta, double* c, const fortran_int_t ldc ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
BLAS_DSYMM( &blas_option< Side >::value, &blas_option< UpLo >::value, &m,
&n, &alpha, a, &lda, b, &ldb, &beta, c, &ldc );
@@ -224,12 +223,12 @@
// * complex<float> value-type.
//
template< typename Order, typename Side, typename UpLo >
-inline void symm( const Order order, const Side side, const UpLo uplo,
- const fortran_int_t m, const fortran_int_t n,
- const std::complex<float> alpha, const std::complex<float>* a,
- const fortran_int_t lda, const std::complex<float>* b,
- const fortran_int_t ldb, const std::complex<float> beta,
- std::complex<float>* c, const fortran_int_t ldc ) {
+inline void symm( const Order, const Side, const UpLo, const fortran_int_t m,
+ const fortran_int_t n, const std::complex<float> alpha,
+ const std::complex<float>* a, const fortran_int_t lda,
+ const std::complex<float>* b, const fortran_int_t ldb,
+ const std::complex<float> beta, std::complex<float>* c,
+ const fortran_int_t ldc ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
BLAS_CSYMM( &blas_option< Side >::value, &blas_option< UpLo >::value, &m,
&n, &alpha, a, &lda, b, &ldb, &beta, c, &ldc );
@@ -241,12 +240,12 @@
// * complex<double> value-type.
//
template< typename Order, typename Side, typename UpLo >
-inline void symm( const Order order, const Side side, const UpLo uplo,
- const fortran_int_t m, const fortran_int_t n,
- const std::complex<double> alpha, const std::complex<double>* a,
- const fortran_int_t lda, const std::complex<double>* b,
- const fortran_int_t ldb, const std::complex<double> beta,
- std::complex<double>* c, const fortran_int_t ldc ) {
+inline void symm( const Order, const Side, const UpLo, const fortran_int_t m,
+ const fortran_int_t n, const std::complex<double> alpha,
+ const std::complex<double>* a, const fortran_int_t lda,
+ const std::complex<double>* b, const fortran_int_t ldb,
+ const std::complex<double> beta, std::complex<double>* c,
+ const fortran_int_t ldc ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
BLAS_ZSYMM( &blas_option< Side >::value, &blas_option< UpLo >::value, &m,
&n, &alpha, a, &lda, b, &ldb, &beta, c, &ldc );
Modified: sandbox/numeric_bindings/boost/numeric/bindings/blas/level3/syr2k.hpp
==============================================================================
--- sandbox/numeric_bindings/boost/numeric/bindings/blas/level3/syr2k.hpp (original)
+++ sandbox/numeric_bindings/boost/numeric/bindings/blas/level3/syr2k.hpp 2012-08-07 09:40:31 EDT (Tue, 07 Aug 2012)
@@ -64,10 +64,10 @@
// * float value-type.
//
template< typename Order, typename UpLo, typename Trans >
-inline void syr2k( const Order order, const UpLo uplo, const Trans trans,
- const int n, const int k, const float alpha, const float* a,
- const int lda, const float* b, const int ldb, const float beta,
- float* c, const int ldc ) {
+inline void syr2k( const Order, const UpLo, const Trans, const int n,
+ const int k, const float alpha, const float* a, const int lda,
+ const float* b, const int ldb, const float beta, float* c,
+ const int ldc ) {
cblas_ssyr2k( cblas_option< Order >::value, cblas_option< UpLo >::value,
cblas_option< Trans >::value, n, k, alpha, a, lda, b, ldb, beta,
c, ldc );
@@ -79,10 +79,10 @@
// * double value-type.
//
template< typename Order, typename UpLo, typename Trans >
-inline void syr2k( const Order order, const UpLo uplo, const Trans trans,
- const int n, const int k, const double alpha, const double* a,
- const int lda, const double* b, const int ldb, const double beta,
- double* c, const int ldc ) {
+inline void syr2k( const Order, const UpLo, const Trans, const int n,
+ const int k, const double alpha, const double* a, const int lda,
+ const double* b, const int ldb, const double beta, double* c,
+ const int ldc ) {
cblas_dsyr2k( cblas_option< Order >::value, cblas_option< UpLo >::value,
cblas_option< Trans >::value, n, k, alpha, a, lda, b, ldb, beta,
c, ldc );
@@ -94,8 +94,8 @@
// * complex<float> value-type.
//
template< typename Order, typename UpLo, typename Trans >
-inline void syr2k( const Order order, const UpLo uplo, const Trans trans,
- const int n, const int k, const std::complex<float> alpha,
+inline void syr2k( const Order, const UpLo, const Trans, const int n,
+ const int k, const std::complex<float> alpha,
const std::complex<float>* a, const int lda,
const std::complex<float>* b, const int ldb,
const std::complex<float> beta, std::complex<float>* c,
@@ -111,8 +111,8 @@
// * complex<double> value-type.
//
template< typename Order, typename UpLo, typename Trans >
-inline void syr2k( const Order order, const UpLo uplo, const Trans trans,
- const int n, const int k, const std::complex<double> alpha,
+inline void syr2k( const Order, const UpLo, const Trans, const int n,
+ const int k, const std::complex<double> alpha,
const std::complex<double>* a, const int lda,
const std::complex<double>* b, const int ldb,
const std::complex<double> beta, std::complex<double>* c,
@@ -129,10 +129,10 @@
// * float value-type.
//
template< typename Order, typename UpLo, typename Trans >
-inline void syr2k( const Order order, const UpLo uplo, const Trans trans,
- const int n, const int k, const float alpha, const float* a,
- const int lda, const float* b, const int ldb, const float beta,
- float* c, const int ldc ) {
+inline void syr2k( const Order, const UpLo, const Trans, const int n,
+ const int k, const float alpha, const float* a, const int lda,
+ const float* b, const int ldb, const float beta, float* c,
+ const int ldc ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
cublasSsyr2k( blas_option< UpLo >::value, blas_option< Trans >::value, n,
k, alpha, a, lda, b, ldb, beta, c, ldc );
@@ -144,10 +144,10 @@
// * double value-type.
//
template< typename Order, typename UpLo, typename Trans >
-inline void syr2k( const Order order, const UpLo uplo, const Trans trans,
- const int n, const int k, const double alpha, const double* a,
- const int lda, const double* b, const int ldb, const double beta,
- double* c, const int ldc ) {
+inline void syr2k( const Order, const UpLo, const Trans, const int n,
+ const int k, const double alpha, const double* a, const int lda,
+ const double* b, const int ldb, const double beta, double* c,
+ const int ldc ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
cublasDsyr2k( blas_option< UpLo >::value, blas_option< Trans >::value, n,
k, alpha, a, lda, b, ldb, beta, c, ldc );
@@ -159,8 +159,8 @@
// * complex<float> value-type.
//
template< typename Order, typename UpLo, typename Trans >
-inline void syr2k( const Order order, const UpLo uplo, const Trans trans,
- const int n, const int k, const std::complex<float> alpha,
+inline void syr2k( const Order, const UpLo, const Trans, const int n,
+ const int k, const std::complex<float> alpha,
const std::complex<float>* a, const int lda,
const std::complex<float>* b, const int ldb,
const std::complex<float> beta, std::complex<float>* c,
@@ -176,14 +176,15 @@
// * complex<double> value-type.
//
template< typename Order, typename UpLo, typename Trans >
-inline void syr2k( const Order order, const UpLo uplo, const Trans trans,
- const int n, const int k, const std::complex<double> alpha,
+inline void syr2k( const Order, const UpLo, const Trans, const int n,
+ const int k, const std::complex<double> alpha,
const std::complex<double>* a, const int lda,
const std::complex<double>* b, const int ldb,
const std::complex<double> beta, std::complex<double>* c,
const int ldc ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
- // NOT FOUND();
+ cublasZsyr2k( blas_option< UpLo >::value, blas_option< Trans >::value, n,
+ k, alpha, a, lda, b, ldb, beta, c, ldc );
}
#else
@@ -193,7 +194,7 @@
// * float value-type.
//
template< typename Order, typename UpLo, typename Trans >
-inline void syr2k( const Order order, const UpLo uplo, const Trans trans,
+inline void syr2k( const Order, const UpLo, const Trans,
const fortran_int_t n, const fortran_int_t k, const float alpha,
const float* a, const fortran_int_t lda, const float* b,
const fortran_int_t ldb, const float beta, float* c,
@@ -209,7 +210,7 @@
// * double value-type.
//
template< typename Order, typename UpLo, typename Trans >
-inline void syr2k( const Order order, const UpLo uplo, const Trans trans,
+inline void syr2k( const Order, const UpLo, const Trans,
const fortran_int_t n, const fortran_int_t k, const double alpha,
const double* a, const fortran_int_t lda, const double* b,
const fortran_int_t ldb, const double beta, double* c,
@@ -225,7 +226,7 @@
// * complex<float> value-type.
//
template< typename Order, typename UpLo, typename Trans >
-inline void syr2k( const Order order, const UpLo uplo, const Trans trans,
+inline void syr2k( const Order, const UpLo, const Trans,
const fortran_int_t n, const fortran_int_t k,
const std::complex<float> alpha, const std::complex<float>* a,
const fortran_int_t lda, const std::complex<float>* b,
@@ -242,7 +243,7 @@
// * complex<double> value-type.
//
template< typename Order, typename UpLo, typename Trans >
-inline void syr2k( const Order order, const UpLo uplo, const Trans trans,
+inline void syr2k( const Order, const UpLo, const Trans,
const fortran_int_t n, const fortran_int_t k,
const std::complex<double> alpha, const std::complex<double>* a,
const fortran_int_t lda, const std::complex<double>* b,
Modified: sandbox/numeric_bindings/boost/numeric/bindings/blas/level3/syrk.hpp
==============================================================================
--- sandbox/numeric_bindings/boost/numeric/bindings/blas/level3/syrk.hpp (original)
+++ sandbox/numeric_bindings/boost/numeric/bindings/blas/level3/syrk.hpp 2012-08-07 09:40:31 EDT (Tue, 07 Aug 2012)
@@ -64,9 +64,9 @@
// * float value-type.
//
template< typename Order, typename UpLo, typename Trans >
-inline void syrk( const Order order, const UpLo uplo, const Trans trans,
- const int n, const int k, const float alpha, const float* a,
- const int lda, const float beta, float* c, const int ldc ) {
+inline void syrk( const Order, const UpLo, const Trans, const int n,
+ const int k, const float alpha, const float* a, const int lda,
+ const float beta, float* c, const int ldc ) {
cblas_ssyrk( cblas_option< Order >::value, cblas_option< UpLo >::value,
cblas_option< Trans >::value, n, k, alpha, a, lda, beta, c, ldc );
}
@@ -77,9 +77,9 @@
// * double value-type.
//
template< typename Order, typename UpLo, typename Trans >
-inline void syrk( const Order order, const UpLo uplo, const Trans trans,
- const int n, const int k, const double alpha, const double* a,
- const int lda, const double beta, double* c, const int ldc ) {
+inline void syrk( const Order, const UpLo, const Trans, const int n,
+ const int k, const double alpha, const double* a, const int lda,
+ const double beta, double* c, const int ldc ) {
cblas_dsyrk( cblas_option< Order >::value, cblas_option< UpLo >::value,
cblas_option< Trans >::value, n, k, alpha, a, lda, beta, c, ldc );
}
@@ -90,8 +90,8 @@
// * complex<float> value-type.
//
template< typename Order, typename UpLo, typename Trans >
-inline void syrk( const Order order, const UpLo uplo, const Trans trans,
- const int n, const int k, const std::complex<float> alpha,
+inline void syrk( const Order, const UpLo, const Trans, const int n,
+ const int k, const std::complex<float> alpha,
const std::complex<float>* a, const int lda,
const std::complex<float> beta, std::complex<float>* c,
const int ldc ) {
@@ -106,8 +106,8 @@
// * complex<double> value-type.
//
template< typename Order, typename UpLo, typename Trans >
-inline void syrk( const Order order, const UpLo uplo, const Trans trans,
- const int n, const int k, const std::complex<double> alpha,
+inline void syrk( const Order, const UpLo, const Trans, const int n,
+ const int k, const std::complex<double> alpha,
const std::complex<double>* a, const int lda,
const std::complex<double> beta, std::complex<double>* c,
const int ldc ) {
@@ -123,9 +123,9 @@
// * float value-type.
//
template< typename Order, typename UpLo, typename Trans >
-inline void syrk( const Order order, const UpLo uplo, const Trans trans,
- const int n, const int k, const float alpha, const float* a,
- const int lda, const float beta, float* c, const int ldc ) {
+inline void syrk( const Order, const UpLo, const Trans, const int n,
+ const int k, const float alpha, const float* a, const int lda,
+ const float beta, float* c, const int ldc ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
cublasSsyrk( blas_option< UpLo >::value, blas_option< Trans >::value, n,
k, alpha, a, lda, beta, c, ldc );
@@ -137,9 +137,9 @@
// * double value-type.
//
template< typename Order, typename UpLo, typename Trans >
-inline void syrk( const Order order, const UpLo uplo, const Trans trans,
- const int n, const int k, const double alpha, const double* a,
- const int lda, const double beta, double* c, const int ldc ) {
+inline void syrk( const Order, const UpLo, const Trans, const int n,
+ const int k, const double alpha, const double* a, const int lda,
+ const double beta, double* c, const int ldc ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
cublasDsyrk( blas_option< UpLo >::value, blas_option< Trans >::value, n,
k, alpha, a, lda, beta, c, ldc );
@@ -151,8 +151,8 @@
// * complex<float> value-type.
//
template< typename Order, typename UpLo, typename Trans >
-inline void syrk( const Order order, const UpLo uplo, const Trans trans,
- const int n, const int k, const std::complex<float> alpha,
+inline void syrk( const Order, const UpLo, const Trans, const int n,
+ const int k, const std::complex<float> alpha,
const std::complex<float>* a, const int lda,
const std::complex<float> beta, std::complex<float>* c,
const int ldc ) {
@@ -167,8 +167,8 @@
// * complex<double> value-type.
//
template< typename Order, typename UpLo, typename Trans >
-inline void syrk( const Order order, const UpLo uplo, const Trans trans,
- const int n, const int k, const std::complex<double> alpha,
+inline void syrk( const Order, const UpLo, const Trans, const int n,
+ const int k, const std::complex<double> alpha,
const std::complex<double>* a, const int lda,
const std::complex<double> beta, std::complex<double>* c,
const int ldc ) {
@@ -184,9 +184,9 @@
// * float value-type.
//
template< typename Order, typename UpLo, typename Trans >
-inline void syrk( const Order order, const UpLo uplo, const Trans trans,
- const fortran_int_t n, const fortran_int_t k, const float alpha,
- const float* a, const fortran_int_t lda, const float beta, float* c,
+inline void syrk( const Order, const UpLo, const Trans, const fortran_int_t n,
+ const fortran_int_t k, const float alpha, const float* a,
+ const fortran_int_t lda, const float beta, float* c,
const fortran_int_t ldc ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
BLAS_SSYRK( &blas_option< UpLo >::value, &blas_option< Trans >::value, &n,
@@ -199,10 +199,10 @@
// * double value-type.
//
template< typename Order, typename UpLo, typename Trans >
-inline void syrk( const Order order, const UpLo uplo, const Trans trans,
- const fortran_int_t n, const fortran_int_t k, const double alpha,
- const double* a, const fortran_int_t lda, const double beta,
- double* c, const fortran_int_t ldc ) {
+inline void syrk( const Order, const UpLo, const Trans, const fortran_int_t n,
+ const fortran_int_t k, const double alpha, const double* a,
+ const fortran_int_t lda, const double beta, double* c,
+ const fortran_int_t ldc ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
BLAS_DSYRK( &blas_option< UpLo >::value, &blas_option< Trans >::value, &n,
&k, &alpha, a, &lda, &beta, c, &ldc );
@@ -214,11 +214,11 @@
// * complex<float> value-type.
//
template< typename Order, typename UpLo, typename Trans >
-inline void syrk( const Order order, const UpLo uplo, const Trans trans,
- const fortran_int_t n, const fortran_int_t k,
- const std::complex<float> alpha, const std::complex<float>* a,
- const fortran_int_t lda, const std::complex<float> beta,
- std::complex<float>* c, const fortran_int_t ldc ) {
+inline void syrk( const Order, const UpLo, const Trans, const fortran_int_t n,
+ const fortran_int_t k, const std::complex<float> alpha,
+ const std::complex<float>* a, const fortran_int_t lda,
+ const std::complex<float> beta, std::complex<float>* c,
+ const fortran_int_t ldc ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
BLAS_CSYRK( &blas_option< UpLo >::value, &blas_option< Trans >::value, &n,
&k, &alpha, a, &lda, &beta, c, &ldc );
@@ -230,11 +230,11 @@
// * complex<double> value-type.
//
template< typename Order, typename UpLo, typename Trans >
-inline void syrk( const Order order, const UpLo uplo, const Trans trans,
- const fortran_int_t n, const fortran_int_t k,
- const std::complex<double> alpha, const std::complex<double>* a,
- const fortran_int_t lda, const std::complex<double> beta,
- std::complex<double>* c, const fortran_int_t ldc ) {
+inline void syrk( const Order, const UpLo, const Trans, const fortran_int_t n,
+ const fortran_int_t k, const std::complex<double> alpha,
+ const std::complex<double>* a, const fortran_int_t lda,
+ const std::complex<double> beta, std::complex<double>* c,
+ const fortran_int_t ldc ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
BLAS_ZSYRK( &blas_option< UpLo >::value, &blas_option< Trans >::value, &n,
&k, &alpha, a, &lda, &beta, c, &ldc );
Modified: sandbox/numeric_bindings/boost/numeric/bindings/blas/level3/trmm.hpp
==============================================================================
--- sandbox/numeric_bindings/boost/numeric/bindings/blas/level3/trmm.hpp (original)
+++ sandbox/numeric_bindings/boost/numeric/bindings/blas/level3/trmm.hpp 2012-08-07 09:40:31 EDT (Tue, 07 Aug 2012)
@@ -66,10 +66,9 @@
//
template< typename Order, typename Side, typename UpLo, typename TransA,
typename Diag >
-inline void trmm( const Order order, const Side side, const UpLo uplo,
- const TransA transa, const Diag diag, const int m, const int n,
- const float alpha, const float* a, const int lda, float* b,
- const int ldb ) {
+inline void trmm( const Order, const Side, const UpLo, const TransA,
+ const Diag, const int m, const int n, const float alpha,
+ const float* a, const int lda, float* b, const int ldb ) {
cblas_strmm( cblas_option< Order >::value, cblas_option< Side >::value,
cblas_option< UpLo >::value, cblas_option< TransA >::value,
cblas_option< Diag >::value, m, n, alpha, a, lda, b, ldb );
@@ -82,10 +81,9 @@
//
template< typename Order, typename Side, typename UpLo, typename TransA,
typename Diag >
-inline void trmm( const Order order, const Side side, const UpLo uplo,
- const TransA transa, const Diag diag, const int m, const int n,
- const double alpha, const double* a, const int lda, double* b,
- const int ldb ) {
+inline void trmm( const Order, const Side, const UpLo, const TransA,
+ const Diag, const int m, const int n, const double alpha,
+ const double* a, const int lda, double* b, const int ldb ) {
cblas_dtrmm( cblas_option< Order >::value, cblas_option< Side >::value,
cblas_option< UpLo >::value, cblas_option< TransA >::value,
cblas_option< Diag >::value, m, n, alpha, a, lda, b, ldb );
@@ -98,10 +96,10 @@
//
template< typename Order, typename Side, typename UpLo, typename TransA,
typename Diag >
-inline void trmm( const Order order, const Side side, const UpLo uplo,
- const TransA transa, const Diag diag, const int m, const int n,
- const std::complex<float> alpha, const std::complex<float>* a,
- const int lda, std::complex<float>* b, const int ldb ) {
+inline void trmm( const Order, const Side, const UpLo, const TransA,
+ const Diag, const int m, const int n, const std::complex<float> alpha,
+ const std::complex<float>* a, const int lda, std::complex<float>* b,
+ const int ldb ) {
cblas_ctrmm( cblas_option< Order >::value, cblas_option< Side >::value,
cblas_option< UpLo >::value, cblas_option< TransA >::value,
cblas_option< Diag >::value, m, n, &alpha, a, lda, b, ldb );
@@ -114,8 +112,8 @@
//
template< typename Order, typename Side, typename UpLo, typename TransA,
typename Diag >
-inline void trmm( const Order order, const Side side, const UpLo uplo,
- const TransA transa, const Diag diag, const int m, const int n,
+inline void trmm( const Order, const Side, const UpLo, const TransA,
+ const Diag, const int m, const int n,
const std::complex<double> alpha, const std::complex<double>* a,
const int lda, std::complex<double>* b, const int ldb ) {
cblas_ztrmm( cblas_option< Order >::value, cblas_option< Side >::value,
@@ -131,10 +129,9 @@
//
template< typename Order, typename Side, typename UpLo, typename TransA,
typename Diag >
-inline void trmm( const Order order, const Side side, const UpLo uplo,
- const TransA transa, const Diag diag, const int m, const int n,
- const float alpha, const float* a, const int lda, float* b,
- const int ldb ) {
+inline void trmm( const Order, const Side, const UpLo, const TransA,
+ const Diag, const int m, const int n, const float alpha,
+ const float* a, const int lda, float* b, const int ldb ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
cublasStrmm( blas_option< Side >::value, blas_option< UpLo >::value,
blas_option< TransA >::value, blas_option< Diag >::value, m, n,
@@ -148,10 +145,9 @@
//
template< typename Order, typename Side, typename UpLo, typename TransA,
typename Diag >
-inline void trmm( const Order order, const Side side, const UpLo uplo,
- const TransA transa, const Diag diag, const int m, const int n,
- const double alpha, const double* a, const int lda, double* b,
- const int ldb ) {
+inline void trmm( const Order, const Side, const UpLo, const TransA,
+ const Diag, const int m, const int n, const double alpha,
+ const double* a, const int lda, double* b, const int ldb ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
cublasDtrmm( blas_option< Side >::value, blas_option< UpLo >::value,
blas_option< TransA >::value, blas_option< Diag >::value, m, n,
@@ -165,10 +161,10 @@
//
template< typename Order, typename Side, typename UpLo, typename TransA,
typename Diag >
-inline void trmm( const Order order, const Side side, const UpLo uplo,
- const TransA transa, const Diag diag, const int m, const int n,
- const std::complex<float> alpha, const std::complex<float>* a,
- const int lda, std::complex<float>* b, const int ldb ) {
+inline void trmm( const Order, const Side, const UpLo, const TransA,
+ const Diag, const int m, const int n, const std::complex<float> alpha,
+ const std::complex<float>* a, const int lda, std::complex<float>* b,
+ const int ldb ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
cublasCtrmm( blas_option< Side >::value, blas_option< UpLo >::value,
blas_option< TransA >::value, blas_option< Diag >::value, m, n,
@@ -182,12 +178,14 @@
//
template< typename Order, typename Side, typename UpLo, typename TransA,
typename Diag >
-inline void trmm( const Order order, const Side side, const UpLo uplo,
- const TransA transa, const Diag diag, const int m, const int n,
+inline void trmm( const Order, const Side, const UpLo, const TransA,
+ const Diag, const int m, const int n,
const std::complex<double> alpha, const std::complex<double>* a,
const int lda, std::complex<double>* b, const int ldb ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
- // NOT FOUND();
+ cublasZtrmm( blas_option< Side >::value, blas_option< UpLo >::value,
+ blas_option< TransA >::value, blas_option< Diag >::value, m, n,
+ alpha, a, lda, b, ldb );
}
#else
@@ -198,10 +196,10 @@
//
template< typename Order, typename Side, typename UpLo, typename TransA,
typename Diag >
-inline void trmm( const Order order, const Side side, const UpLo uplo,
- const TransA transa, const Diag diag, const fortran_int_t m,
- const fortran_int_t n, const float alpha, const float* a,
- const fortran_int_t lda, float* b, const fortran_int_t ldb ) {
+inline void trmm( const Order, const Side, const UpLo, const TransA,
+ const Diag, const fortran_int_t m, const fortran_int_t n,
+ const float alpha, const float* a, const fortran_int_t lda, float* b,
+ const fortran_int_t ldb ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
BLAS_STRMM( &blas_option< Side >::value, &blas_option< UpLo >::value,
&blas_option< TransA >::value, &blas_option< Diag >::value, &m,
@@ -215,10 +213,10 @@
//
template< typename Order, typename Side, typename UpLo, typename TransA,
typename Diag >
-inline void trmm( const Order order, const Side side, const UpLo uplo,
- const TransA transa, const Diag diag, const fortran_int_t m,
- const fortran_int_t n, const double alpha, const double* a,
- const fortran_int_t lda, double* b, const fortran_int_t ldb ) {
+inline void trmm( const Order, const Side, const UpLo, const TransA,
+ const Diag, const fortran_int_t m, const fortran_int_t n,
+ const double alpha, const double* a, const fortran_int_t lda,
+ double* b, const fortran_int_t ldb ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
BLAS_DTRMM( &blas_option< Side >::value, &blas_option< UpLo >::value,
&blas_option< TransA >::value, &blas_option< Diag >::value, &m,
@@ -232,11 +230,11 @@
//
template< typename Order, typename Side, typename UpLo, typename TransA,
typename Diag >
-inline void trmm( const Order order, const Side side, const UpLo uplo,
- const TransA transa, const Diag diag, const fortran_int_t m,
- const fortran_int_t n, const std::complex<float> alpha,
- const std::complex<float>* a, const fortran_int_t lda,
- std::complex<float>* b, const fortran_int_t ldb ) {
+inline void trmm( const Order, const Side, const UpLo, const TransA,
+ const Diag, const fortran_int_t m, const fortran_int_t n,
+ const std::complex<float> alpha, const std::complex<float>* a,
+ const fortran_int_t lda, std::complex<float>* b,
+ const fortran_int_t ldb ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
BLAS_CTRMM( &blas_option< Side >::value, &blas_option< UpLo >::value,
&blas_option< TransA >::value, &blas_option< Diag >::value, &m,
@@ -250,11 +248,11 @@
//
template< typename Order, typename Side, typename UpLo, typename TransA,
typename Diag >
-inline void trmm( const Order order, const Side side, const UpLo uplo,
- const TransA transa, const Diag diag, const fortran_int_t m,
- const fortran_int_t n, const std::complex<double> alpha,
- const std::complex<double>* a, const fortran_int_t lda,
- std::complex<double>* b, const fortran_int_t ldb ) {
+inline void trmm( const Order, const Side, const UpLo, const TransA,
+ const Diag, const fortran_int_t m, const fortran_int_t n,
+ const std::complex<double> alpha, const std::complex<double>* a,
+ const fortran_int_t lda, std::complex<double>* b,
+ const fortran_int_t ldb ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
BLAS_ZTRMM( &blas_option< Side >::value, &blas_option< UpLo >::value,
&blas_option< TransA >::value, &blas_option< Diag >::value, &m,
Modified: sandbox/numeric_bindings/boost/numeric/bindings/blas/level3/trsm.hpp
==============================================================================
--- sandbox/numeric_bindings/boost/numeric/bindings/blas/level3/trsm.hpp (original)
+++ sandbox/numeric_bindings/boost/numeric/bindings/blas/level3/trsm.hpp 2012-08-07 09:40:31 EDT (Tue, 07 Aug 2012)
@@ -66,10 +66,9 @@
//
template< typename Order, typename Side, typename UpLo, typename TransA,
typename Diag >
-inline void trsm( const Order order, const Side side, const UpLo uplo,
- const TransA transa, const Diag diag, const int m, const int n,
- const float alpha, const float* a, const int lda, float* b,
- const int ldb ) {
+inline void trsm( const Order, const Side, const UpLo, const TransA,
+ const Diag, const int m, const int n, const float alpha,
+ const float* a, const int lda, float* b, const int ldb ) {
cblas_strsm( cblas_option< Order >::value, cblas_option< Side >::value,
cblas_option< UpLo >::value, cblas_option< TransA >::value,
cblas_option< Diag >::value, m, n, alpha, a, lda, b, ldb );
@@ -82,10 +81,9 @@
//
template< typename Order, typename Side, typename UpLo, typename TransA,
typename Diag >
-inline void trsm( const Order order, const Side side, const UpLo uplo,
- const TransA transa, const Diag diag, const int m, const int n,
- const double alpha, const double* a, const int lda, double* b,
- const int ldb ) {
+inline void trsm( const Order, const Side, const UpLo, const TransA,
+ const Diag, const int m, const int n, const double alpha,
+ const double* a, const int lda, double* b, const int ldb ) {
cblas_dtrsm( cblas_option< Order >::value, cblas_option< Side >::value,
cblas_option< UpLo >::value, cblas_option< TransA >::value,
cblas_option< Diag >::value, m, n, alpha, a, lda, b, ldb );
@@ -98,10 +96,10 @@
//
template< typename Order, typename Side, typename UpLo, typename TransA,
typename Diag >
-inline void trsm( const Order order, const Side side, const UpLo uplo,
- const TransA transa, const Diag diag, const int m, const int n,
- const std::complex<float> alpha, const std::complex<float>* a,
- const int lda, std::complex<float>* b, const int ldb ) {
+inline void trsm( const Order, const Side, const UpLo, const TransA,
+ const Diag, const int m, const int n, const std::complex<float> alpha,
+ const std::complex<float>* a, const int lda, std::complex<float>* b,
+ const int ldb ) {
cblas_ctrsm( cblas_option< Order >::value, cblas_option< Side >::value,
cblas_option< UpLo >::value, cblas_option< TransA >::value,
cblas_option< Diag >::value, m, n, &alpha, a, lda, b, ldb );
@@ -114,8 +112,8 @@
//
template< typename Order, typename Side, typename UpLo, typename TransA,
typename Diag >
-inline void trsm( const Order order, const Side side, const UpLo uplo,
- const TransA transa, const Diag diag, const int m, const int n,
+inline void trsm( const Order, const Side, const UpLo, const TransA,
+ const Diag, const int m, const int n,
const std::complex<double> alpha, const std::complex<double>* a,
const int lda, std::complex<double>* b, const int ldb ) {
cblas_ztrsm( cblas_option< Order >::value, cblas_option< Side >::value,
@@ -131,10 +129,9 @@
//
template< typename Order, typename Side, typename UpLo, typename TransA,
typename Diag >
-inline void trsm( const Order order, const Side side, const UpLo uplo,
- const TransA transa, const Diag diag, const int m, const int n,
- const float alpha, const float* a, const int lda, float* b,
- const int ldb ) {
+inline void trsm( const Order, const Side, const UpLo, const TransA,
+ const Diag, const int m, const int n, const float alpha,
+ const float* a, const int lda, float* b, const int ldb ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
cublasStrsm( blas_option< Side >::value, blas_option< UpLo >::value,
blas_option< TransA >::value, blas_option< Diag >::value, m, n,
@@ -148,10 +145,9 @@
//
template< typename Order, typename Side, typename UpLo, typename TransA,
typename Diag >
-inline void trsm( const Order order, const Side side, const UpLo uplo,
- const TransA transa, const Diag diag, const int m, const int n,
- const double alpha, const double* a, const int lda, double* b,
- const int ldb ) {
+inline void trsm( const Order, const Side, const UpLo, const TransA,
+ const Diag, const int m, const int n, const double alpha,
+ const double* a, const int lda, double* b, const int ldb ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
cublasDtrsm( blas_option< Side >::value, blas_option< UpLo >::value,
blas_option< TransA >::value, blas_option< Diag >::value, m, n,
@@ -165,10 +161,10 @@
//
template< typename Order, typename Side, typename UpLo, typename TransA,
typename Diag >
-inline void trsm( const Order order, const Side side, const UpLo uplo,
- const TransA transa, const Diag diag, const int m, const int n,
- const std::complex<float> alpha, const std::complex<float>* a,
- const int lda, std::complex<float>* b, const int ldb ) {
+inline void trsm( const Order, const Side, const UpLo, const TransA,
+ const Diag, const int m, const int n, const std::complex<float> alpha,
+ const std::complex<float>* a, const int lda, std::complex<float>* b,
+ const int ldb ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
cublasCtrsm( blas_option< Side >::value, blas_option< UpLo >::value,
blas_option< TransA >::value, blas_option< Diag >::value, m, n,
@@ -182,8 +178,8 @@
//
template< typename Order, typename Side, typename UpLo, typename TransA,
typename Diag >
-inline void trsm( const Order order, const Side side, const UpLo uplo,
- const TransA transa, const Diag diag, const int m, const int n,
+inline void trsm( const Order, const Side, const UpLo, const TransA,
+ const Diag, const int m, const int n,
const std::complex<double> alpha, const std::complex<double>* a,
const int lda, std::complex<double>* b, const int ldb ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
@@ -200,10 +196,10 @@
//
template< typename Order, typename Side, typename UpLo, typename TransA,
typename Diag >
-inline void trsm( const Order order, const Side side, const UpLo uplo,
- const TransA transa, const Diag diag, const fortran_int_t m,
- const fortran_int_t n, const float alpha, const float* a,
- const fortran_int_t lda, float* b, const fortran_int_t ldb ) {
+inline void trsm( const Order, const Side, const UpLo, const TransA,
+ const Diag, const fortran_int_t m, const fortran_int_t n,
+ const float alpha, const float* a, const fortran_int_t lda, float* b,
+ const fortran_int_t ldb ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
BLAS_STRSM( &blas_option< Side >::value, &blas_option< UpLo >::value,
&blas_option< TransA >::value, &blas_option< Diag >::value, &m,
@@ -217,10 +213,10 @@
//
template< typename Order, typename Side, typename UpLo, typename TransA,
typename Diag >
-inline void trsm( const Order order, const Side side, const UpLo uplo,
- const TransA transa, const Diag diag, const fortran_int_t m,
- const fortran_int_t n, const double alpha, const double* a,
- const fortran_int_t lda, double* b, const fortran_int_t ldb ) {
+inline void trsm( const Order, const Side, const UpLo, const TransA,
+ const Diag, const fortran_int_t m, const fortran_int_t n,
+ const double alpha, const double* a, const fortran_int_t lda,
+ double* b, const fortran_int_t ldb ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
BLAS_DTRSM( &blas_option< Side >::value, &blas_option< UpLo >::value,
&blas_option< TransA >::value, &blas_option< Diag >::value, &m,
@@ -234,11 +230,11 @@
//
template< typename Order, typename Side, typename UpLo, typename TransA,
typename Diag >
-inline void trsm( const Order order, const Side side, const UpLo uplo,
- const TransA transa, const Diag diag, const fortran_int_t m,
- const fortran_int_t n, const std::complex<float> alpha,
- const std::complex<float>* a, const fortran_int_t lda,
- std::complex<float>* b, const fortran_int_t ldb ) {
+inline void trsm( const Order, const Side, const UpLo, const TransA,
+ const Diag, const fortran_int_t m, const fortran_int_t n,
+ const std::complex<float> alpha, const std::complex<float>* a,
+ const fortran_int_t lda, std::complex<float>* b,
+ const fortran_int_t ldb ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
BLAS_CTRSM( &blas_option< Side >::value, &blas_option< UpLo >::value,
&blas_option< TransA >::value, &blas_option< Diag >::value, &m,
@@ -252,11 +248,11 @@
//
template< typename Order, typename Side, typename UpLo, typename TransA,
typename Diag >
-inline void trsm( const Order order, const Side side, const UpLo uplo,
- const TransA transa, const Diag diag, const fortran_int_t m,
- const fortran_int_t n, const std::complex<double> alpha,
- const std::complex<double>* a, const fortran_int_t lda,
- std::complex<double>* b, const fortran_int_t ldb ) {
+inline void trsm( const Order, const Side, const UpLo, const TransA,
+ const Diag, const fortran_int_t m, const fortran_int_t n,
+ const std::complex<double> alpha, const std::complex<double>* a,
+ const fortran_int_t lda, std::complex<double>* b,
+ const fortran_int_t ldb ) {
BOOST_STATIC_ASSERT( (is_same<Order, tag::column_major>::value) );
BLAS_ZTRSM( &blas_option< Side >::value, &blas_option< UpLo >::value,
&blas_option< TransA >::value, &blas_option< Diag >::value, &m,
Boost-Commit list run by bdawes at acm.org, david.abrahams at rcn.com, gregod at cs.rpi.edu, cpdaniel at pacbell.net, john at johnmaddock.co.uk