$include_dir="/home/hyper-archives/boost-commit/include"; include("$include_dir/msg-header.inc") ?>
Subject: [Boost-commit] svn:boost r56263 - in sandbox/numeric_bindings/boost/numeric/bindings/blas: level1 level2 level3
From: rutger_at_[hidden]
Date: 2009-09-17 07:18:47
Author: rutger
Date: 2009-09-17 07:18:43 EDT (Thu, 17 Sep 2009)
New Revision: 56263
URL: http://svn.boost.org/trac/boost/changeset/56263
Log:
Added support for detection of number of columns of op(A)
Text files modified: 
   sandbox/numeric_bindings/boost/numeric/bindings/blas/level1/asum.hpp  |     4 ++--                                    
   sandbox/numeric_bindings/boost/numeric/bindings/blas/level1/axpy.hpp  |     8 ++++----                                
   sandbox/numeric_bindings/boost/numeric/bindings/blas/level1/copy.hpp  |    16 ++++++++--------                        
   sandbox/numeric_bindings/boost/numeric/bindings/blas/level1/dot.hpp   |     8 ++++----                                
   sandbox/numeric_bindings/boost/numeric/bindings/blas/level1/dotc.hpp  |     8 ++++----                                
   sandbox/numeric_bindings/boost/numeric/bindings/blas/level1/dotu.hpp  |     8 ++++----                                
   sandbox/numeric_bindings/boost/numeric/bindings/blas/level1/drot.hpp  |     4 ++--                                    
   sandbox/numeric_bindings/boost/numeric/bindings/blas/level1/nrm2.hpp  |     4 ++--                                    
   sandbox/numeric_bindings/boost/numeric/bindings/blas/level1/rot.hpp   |     4 ++--                                    
   sandbox/numeric_bindings/boost/numeric/bindings/blas/level1/scal.hpp  |     8 ++++----                                
   sandbox/numeric_bindings/boost/numeric/bindings/blas/level1/sdot.hpp  |     4 ++--                                    
   sandbox/numeric_bindings/boost/numeric/bindings/blas/level1/srot.hpp  |     4 ++--                                    
   sandbox/numeric_bindings/boost/numeric/bindings/blas/level2/gbmv.hpp  |    12 ++++++------                            
   sandbox/numeric_bindings/boost/numeric/bindings/blas/level2/gemv.hpp  |    16 ++++++++--------                        
   sandbox/numeric_bindings/boost/numeric/bindings/blas/level2/ger.hpp   |     4 ++--                                    
   sandbox/numeric_bindings/boost/numeric/bindings/blas/level2/gerc.hpp  |     8 ++++----                                
   sandbox/numeric_bindings/boost/numeric/bindings/blas/level2/geru.hpp  |     8 ++++----                                
   sandbox/numeric_bindings/boost/numeric/bindings/blas/level2/hbmv.hpp  |     8 ++++----                                
   sandbox/numeric_bindings/boost/numeric/bindings/blas/level2/hemv.hpp  |     8 ++++----                                
   sandbox/numeric_bindings/boost/numeric/bindings/blas/level2/her.hpp   |     4 ++--                                    
   sandbox/numeric_bindings/boost/numeric/bindings/blas/level2/her2.hpp  |     8 ++++----                                
   sandbox/numeric_bindings/boost/numeric/bindings/blas/level2/hpmv.hpp  |     8 ++++----                                
   sandbox/numeric_bindings/boost/numeric/bindings/blas/level2/hpr.hpp   |     4 ++--                                    
   sandbox/numeric_bindings/boost/numeric/bindings/blas/level2/hpr2.hpp  |     8 ++++----                                
   sandbox/numeric_bindings/boost/numeric/bindings/blas/level2/sbmv.hpp  |     8 ++++----                                
   sandbox/numeric_bindings/boost/numeric/bindings/blas/level2/spmv.hpp  |     4 ++--                                    
   sandbox/numeric_bindings/boost/numeric/bindings/blas/level2/spr.hpp   |     4 ++--                                    
   sandbox/numeric_bindings/boost/numeric/bindings/blas/level2/spr2.hpp  |     4 ++--                                    
   sandbox/numeric_bindings/boost/numeric/bindings/blas/level2/symv.hpp  |     4 ++--                                    
   sandbox/numeric_bindings/boost/numeric/bindings/blas/level2/syr.hpp   |     4 ++--                                    
   sandbox/numeric_bindings/boost/numeric/bindings/blas/level2/syr2.hpp  |     4 ++--                                    
   sandbox/numeric_bindings/boost/numeric/bindings/blas/level2/tbmv.hpp  |     8 ++++----                                
   sandbox/numeric_bindings/boost/numeric/bindings/blas/level2/tbsv.hpp  |     8 ++++----                                
   sandbox/numeric_bindings/boost/numeric/bindings/blas/level2/tpmv.hpp  |     8 ++++----                                
   sandbox/numeric_bindings/boost/numeric/bindings/blas/level2/tpsv.hpp  |     8 ++++----                                
   sandbox/numeric_bindings/boost/numeric/bindings/blas/level2/trmv.hpp  |     8 ++++----                                
   sandbox/numeric_bindings/boost/numeric/bindings/blas/level2/trsv.hpp  |     8 ++++----                                
   sandbox/numeric_bindings/boost/numeric/bindings/blas/level3/gemm.hpp  |    33 +++++++++++++++++----------------       
   sandbox/numeric_bindings/boost/numeric/bindings/blas/level3/hemm.hpp  |     8 ++++----                                
   sandbox/numeric_bindings/boost/numeric/bindings/blas/level3/her2k.hpp |     8 ++++----                                
   sandbox/numeric_bindings/boost/numeric/bindings/blas/level3/herk.hpp  |     4 ++--                                    
   sandbox/numeric_bindings/boost/numeric/bindings/blas/level3/symm.hpp  |    16 ++++++++--------                        
   sandbox/numeric_bindings/boost/numeric/bindings/blas/level3/syr2k.hpp |    16 ++++++++--------                        
   sandbox/numeric_bindings/boost/numeric/bindings/blas/level3/syrk.hpp  |     8 ++++----                                
   sandbox/numeric_bindings/boost/numeric/bindings/blas/level3/trmm.hpp  |     8 ++++----                                
   sandbox/numeric_bindings/boost/numeric/bindings/blas/level3/trsm.hpp  |     8 ++++----                                
   46 files changed, 183 insertions(+), 182 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	2009-09-17 07:18:43 EDT (Thu, 17 Sep 2009)
@@ -29,11 +29,11 @@
 
 // overloaded functions to call blas
 namespace detail {
-    inline float asum( const integer_t n, float const* x,
+    inline float asum( const integer_t n, const float* x,
             const integer_t incx ) {
         return BLAS_SASUM( &n, x, &incx );
     }
-    inline double asum( const integer_t n, double const* x,
+    inline double asum( const integer_t n, const double* x,
             const integer_t incx ) {
         return BLAS_DASUM( &n, x, &incx );
     }
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	2009-09-17 07:18:43 EDT (Thu, 17 Sep 2009)
@@ -29,22 +29,22 @@
 
 // overloaded functions to call blas
 namespace detail {
-    inline void axpy( const integer_t n, const float a, float const* x,
+    inline void axpy( const integer_t n, const float a, const float* x,
             const integer_t incx, float* y, const integer_t incy ) {
         BLAS_SAXPY( &n, &a, x, &incx, y, &incy );
     }
-    inline void axpy( const integer_t n, const double a, double const* x,
+    inline void axpy( const integer_t n, const double a, const double* x,
             const integer_t incx, double* y, const integer_t incy ) {
         BLAS_DAXPY( &n, &a, x, &incx, y, &incy );
     }
     inline void axpy( const integer_t n, const traits::complex_f a,
-            traits::complex_f const* x, const integer_t incx,
+            const traits::complex_f* x, const integer_t incx,
             traits::complex_f* y, const integer_t incy ) {
         BLAS_CAXPY( &n, traits::complex_ptr(&a), traits::complex_ptr(x),
                 &incx, traits::complex_ptr(y), &incy );
     }
     inline void axpy( const integer_t n, const traits::complex_d a,
-            traits::complex_d const* x, const integer_t incx,
+            const traits::complex_d* x, const integer_t incx,
             traits::complex_d* y, const integer_t incy ) {
         BLAS_ZAXPY( &n, traits::complex_ptr(&a), traits::complex_ptr(x),
                 &incx, traits::complex_ptr(y), &incy );
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	2009-09-17 07:18:43 EDT (Thu, 17 Sep 2009)
@@ -29,22 +29,22 @@
 
 // overloaded functions to call blas
 namespace detail {
-    inline void copy( const integer_t n, float const* x, const integer_t incx,
-            float const* y, const integer_t incy ) {
+    inline void copy( const integer_t n, const float* x, const integer_t incx,
+            const float* y, const integer_t incy ) {
         BLAS_SCOPY( &n, x, &incx, y, &incy );
     }
-    inline void copy( const integer_t n, double const* x,
-            const integer_t incx, double const* y, const integer_t incy ) {
+    inline void copy( const integer_t n, const double* x,
+            const integer_t incx, const double* y, const integer_t incy ) {
         BLAS_DCOPY( &n, x, &incx, y, &incy );
     }
-    inline void copy( const integer_t n, traits::complex_f const* x,
-            const integer_t incx, traits::complex_f const* y,
+    inline void copy( const integer_t n, const traits::complex_f* x,
+            const integer_t incx, const traits::complex_f* y,
             const integer_t incy ) {
         BLAS_CCOPY( &n, traits::complex_ptr(x), &incx, traits::complex_ptr(y),
                 &incy );
     }
-    inline void copy( const integer_t n, traits::complex_d const* x,
-            const integer_t incx, traits::complex_d const* y,
+    inline void copy( const integer_t n, const traits::complex_d* x,
+            const integer_t incx, const traits::complex_d* y,
             const integer_t incy ) {
         BLAS_ZCOPY( &n, traits::complex_ptr(x), &incx, traits::complex_ptr(y),
                 &incy );
Modified: sandbox/numeric_bindings/boost/numeric/bindings/blas/level1/dot.hpp
==============================================================================
--- sandbox/numeric_bindings/boost/numeric/bindings/blas/level1/dot.hpp	(original)
+++ sandbox/numeric_bindings/boost/numeric/bindings/blas/level1/dot.hpp	2009-09-17 07:18:43 EDT (Thu, 17 Sep 2009)
@@ -29,12 +29,12 @@
 
 // overloaded functions to call blas
 namespace detail {
-    inline float dot( const integer_t n, float const* x, const integer_t incx,
-            float const* y, const integer_t incy ) {
+    inline float dot( const integer_t n, const float* x, const integer_t incx,
+            const float* y, const integer_t incy ) {
         return BLAS_SDOT( &n, x, &incx, y, &incy );
     }
-    inline double dot( const integer_t n, double const* x,
-            const integer_t incx, double const* y, const integer_t incy ) {
+    inline double dot( const integer_t n, const double* x,
+            const integer_t incx, const double* y, const integer_t incy ) {
         return BLAS_DDOT( &n, x, &incx, y, &incy );
     }
 }
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	2009-09-17 07:18:43 EDT (Thu, 17 Sep 2009)
@@ -29,14 +29,14 @@
 
 // overloaded functions to call blas
 namespace detail {
-    inline fcomplex_t dotc( const integer_t n, traits::complex_f const* x,
-            const integer_t incx, traits::complex_f const* y,
+    inline fcomplex_t dotc( const integer_t n, const traits::complex_f* x,
+            const integer_t incx, const traits::complex_f* y,
             const integer_t incy ) {
         return BLAS_CDOTC( &n, traits::complex_ptr(x), &incx,
                 traits::complex_ptr(y), &incy );
     }
-    inline dcomplex_t dotc( const integer_t n, traits::complex_d const* x,
-            const integer_t incx, traits::complex_d const* y,
+    inline dcomplex_t dotc( const integer_t n, const traits::complex_d* x,
+            const integer_t incx, const traits::complex_d* y,
             const integer_t incy ) {
         return BLAS_ZDOTC( &n, traits::complex_ptr(x), &incx,
                 traits::complex_ptr(y), &incy );
Modified: sandbox/numeric_bindings/boost/numeric/bindings/blas/level1/dotu.hpp
==============================================================================
--- sandbox/numeric_bindings/boost/numeric/bindings/blas/level1/dotu.hpp	(original)
+++ sandbox/numeric_bindings/boost/numeric/bindings/blas/level1/dotu.hpp	2009-09-17 07:18:43 EDT (Thu, 17 Sep 2009)
@@ -29,14 +29,14 @@
 
 // overloaded functions to call blas
 namespace detail {
-    inline fcomplex_t dotu( const integer_t n, traits::complex_f const* x,
-            const integer_t incx, traits::complex_f const* y,
+    inline fcomplex_t dotu( const integer_t n, const traits::complex_f* x,
+            const integer_t incx, const traits::complex_f* y,
             const integer_t incy ) {
         return BLAS_CDOTU( &n, traits::complex_ptr(x), &incx,
                 traits::complex_ptr(y), &incy );
     }
-    inline dcomplex_t dotu( const integer_t n, traits::complex_d const* x,
-            const integer_t incx, traits::complex_d const* y,
+    inline dcomplex_t dotu( const integer_t n, const traits::complex_d* x,
+            const integer_t incx, const traits::complex_d* y,
             const integer_t incy ) {
         return BLAS_ZDOTU( &n, traits::complex_ptr(x), &incx,
                 traits::complex_ptr(y), &incy );
Modified: sandbox/numeric_bindings/boost/numeric/bindings/blas/level1/drot.hpp
==============================================================================
--- sandbox/numeric_bindings/boost/numeric/bindings/blas/level1/drot.hpp	(original)
+++ sandbox/numeric_bindings/boost/numeric/bindings/blas/level1/drot.hpp	2009-09-17 07:18:43 EDT (Thu, 17 Sep 2009)
@@ -29,8 +29,8 @@
 
 // overloaded functions to call blas
 namespace detail {
-    inline void drot( const integer_t n, traits::complex_d const* cx,
-            const integer_t incx, traits::complex_d const* cy,
+    inline void drot( const integer_t n, const traits::complex_d* cx,
+            const integer_t incx, const traits::complex_d* cy,
             const integer_t incy, const double c, const double s ) {
         BLAS_ZDROT( &n, traits::complex_ptr(cx), &incx,
                 traits::complex_ptr(cy), &incy, &c, &s );
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	2009-09-17 07:18:43 EDT (Thu, 17 Sep 2009)
@@ -29,11 +29,11 @@
 
 // overloaded functions to call blas
 namespace detail {
-    inline float nrm2( const integer_t n, float const* x,
+    inline float nrm2( const integer_t n, const float* x,
             const integer_t incx ) {
         return BLAS_SNRM2( &n, x, &incx );
     }
-    inline double nrm2( const integer_t n, double const* x,
+    inline double nrm2( const integer_t n, const double* x,
             const integer_t incx ) {
         return BLAS_DNRM2( &n, x, &incx );
     }
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	2009-09-17 07:18:43 EDT (Thu, 17 Sep 2009)
@@ -29,11 +29,11 @@
 
 // overloaded functions to call blas
 namespace detail {
-    inline void rot( const integer_t n, float const* x, const integer_t incx,
+    inline void rot( const integer_t n, const float* x, const integer_t incx,
             float* y, const integer_t incy, const float c, const float s ) {
         BLAS_SROT( &n, x, &incx, y, &incy, &c, &s );
     }
-    inline void rot( const integer_t n, double const* x, const integer_t incx,
+    inline void rot( const integer_t n, const double* x, const integer_t incx,
             double* y, const integer_t incy, const double c, const double s ) {
         BLAS_DROT( &n, x, &incx, y, &incy, &c, &s );
     }
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	2009-09-17 07:18:43 EDT (Thu, 17 Sep 2009)
@@ -29,21 +29,21 @@
 
 // overloaded functions to call blas
 namespace detail {
-    inline void scal( const integer_t n, const float a, float const* x,
+    inline void scal( const integer_t n, const float a, const float* x,
             const integer_t incx ) {
         BLAS_SSCAL( &n, &a, x, &incx );
     }
-    inline void scal( const integer_t n, const double a, double const* x,
+    inline void scal( const integer_t n, const double a, const double* x,
             const integer_t incx ) {
         BLAS_DSCAL( &n, &a, x, &incx );
     }
     inline void scal( const integer_t n, const traits::complex_f a,
-            traits::complex_f const* x, const integer_t incx ) {
+            const traits::complex_f* x, const integer_t incx ) {
         BLAS_CSCAL( &n, traits::complex_ptr(&a), traits::complex_ptr(x),
                 &incx );
     }
     inline void scal( const integer_t n, const traits::complex_d a,
-            traits::complex_d const* x, const integer_t incx ) {
+            const traits::complex_d* x, const integer_t incx ) {
         BLAS_ZSCAL( &n, traits::complex_ptr(&a), traits::complex_ptr(x),
                 &incx );
     }
Modified: sandbox/numeric_bindings/boost/numeric/bindings/blas/level1/sdot.hpp
==============================================================================
--- sandbox/numeric_bindings/boost/numeric/bindings/blas/level1/sdot.hpp	(original)
+++ sandbox/numeric_bindings/boost/numeric/bindings/blas/level1/sdot.hpp	2009-09-17 07:18:43 EDT (Thu, 17 Sep 2009)
@@ -29,8 +29,8 @@
 
 // overloaded functions to call blas
 namespace detail {
-    inline double sdot( const integer_t n, float const* sx,
-            const integer_t incx, float const* sy, const integer_t incy ) {
+    inline double sdot( const integer_t n, const float* sx,
+            const integer_t incx, const float* sy, const integer_t incy ) {
         return BLAS_DSDOT( &n, sx, &incx, sy, &incy );
     }
 }
Modified: sandbox/numeric_bindings/boost/numeric/bindings/blas/level1/srot.hpp
==============================================================================
--- sandbox/numeric_bindings/boost/numeric/bindings/blas/level1/srot.hpp	(original)
+++ sandbox/numeric_bindings/boost/numeric/bindings/blas/level1/srot.hpp	2009-09-17 07:18:43 EDT (Thu, 17 Sep 2009)
@@ -29,8 +29,8 @@
 
 // overloaded functions to call blas
 namespace detail {
-    inline void srot( const integer_t n, traits::complex_f const* cx,
-            const integer_t incx, traits::complex_f const* cy,
+    inline void srot( const integer_t n, const traits::complex_f* cx,
+            const integer_t incx, const traits::complex_f* cy,
             const integer_t incy, const float c, const float s ) {
         BLAS_CSROT( &n, traits::complex_ptr(cx), &incx,
                 traits::complex_ptr(cy), &incy, &c, &s );
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	2009-09-17 07:18:43 EDT (Thu, 17 Sep 2009)
@@ -31,7 +31,7 @@
 namespace detail {
     inline void gbmv( const char trans, const integer_t m, const integer_t n,
             const integer_t kl, const integer_t ku, const float alpha,
-            float const* a, const integer_t lda, float const* x,
+            const float* a, const integer_t lda, const float* x,
             const integer_t incx, const float beta, float* y,
             const integer_t incy ) {
         BLAS_SGBMV( &trans, &m, &n, &kl, &ku, &alpha, a, &lda, x, &incx,
@@ -39,7 +39,7 @@
     }
     inline void gbmv( const char trans, const integer_t m, const integer_t n,
             const integer_t kl, const integer_t ku, const double alpha,
-            double const* a, const integer_t lda, double const* x,
+            const double* a, const integer_t lda, const double* x,
             const integer_t incx, const double beta, double* y,
             const integer_t incy ) {
         BLAS_DGBMV( &trans, &m, &n, &kl, &ku, &alpha, a, &lda, x, &incx,
@@ -47,8 +47,8 @@
     }
     inline void gbmv( const char trans, const integer_t m, const integer_t n,
             const integer_t kl, const integer_t ku,
-            const traits::complex_f alpha, traits::complex_f const* a,
-            const integer_t lda, traits::complex_f const* x,
+            const traits::complex_f alpha, const traits::complex_f* a,
+            const integer_t lda, const traits::complex_f* x,
             const integer_t incx, const traits::complex_f beta,
             traits::complex_f* y, const integer_t incy ) {
         BLAS_CGBMV( &trans, &m, &n, &kl, &ku, traits::complex_ptr(&alpha),
@@ -57,8 +57,8 @@
     }
     inline void gbmv( const char trans, const integer_t m, const integer_t n,
             const integer_t kl, const integer_t ku,
-            const traits::complex_d alpha, traits::complex_d const* a,
-            const integer_t lda, traits::complex_d const* x,
+            const traits::complex_d alpha, const traits::complex_d* a,
+            const integer_t lda, const traits::complex_d* x,
             const integer_t incx, const traits::complex_d beta,
             traits::complex_d* y, const integer_t incy ) {
         BLAS_ZGBMV( &trans, &m, &n, &kl, &ku, traits::complex_ptr(&alpha),
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	2009-09-17 07:18:43 EDT (Thu, 17 Sep 2009)
@@ -30,22 +30,22 @@
 // overloaded functions to call blas
 namespace detail {
     inline void gemv( const char trans, const integer_t m, const integer_t n,
-            const float alpha, float const* a, const integer_t lda,
-            float const* x, const integer_t incx, const float beta, float* y,
+            const float alpha, const float* a, const integer_t lda,
+            const float* x, const integer_t incx, const float beta, float* y,
             const integer_t incy ) {
         BLAS_SGEMV( &trans, &m, &n, &alpha, a, &lda, x, &incx, &beta, y,
                 &incy );
     }
     inline void gemv( const char trans, const integer_t m, const integer_t n,
-            const double alpha, double const* a, const integer_t lda,
-            double const* x, const integer_t incx, const double beta,
+            const double alpha, const double* a, const integer_t lda,
+            const double* x, const integer_t incx, const double beta,
             double* y, const integer_t incy ) {
         BLAS_DGEMV( &trans, &m, &n, &alpha, a, &lda, x, &incx, &beta, y,
                 &incy );
     }
     inline void gemv( const char trans, const integer_t m, const integer_t n,
-            const traits::complex_f alpha, traits::complex_f const* a,
-            const integer_t lda, traits::complex_f const* x,
+            const traits::complex_f alpha, const traits::complex_f* a,
+            const integer_t lda, const traits::complex_f* x,
             const integer_t incx, const traits::complex_f beta,
             traits::complex_f* y, const integer_t incy ) {
         BLAS_CGEMV( &trans, &m, &n, traits::complex_ptr(&alpha),
@@ -53,8 +53,8 @@
                 traits::complex_ptr(&beta), traits::complex_ptr(y), &incy );
     }
     inline void gemv( const char trans, const integer_t m, const integer_t n,
-            const traits::complex_d alpha, traits::complex_d const* a,
-            const integer_t lda, traits::complex_d const* x,
+            const traits::complex_d alpha, const traits::complex_d* a,
+            const integer_t lda, const traits::complex_d* x,
             const integer_t incx, const traits::complex_d beta,
             traits::complex_d* y, const integer_t incy ) {
         BLAS_ZGEMV( &trans, &m, &n, traits::complex_ptr(&alpha),
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	2009-09-17 07:18:43 EDT (Thu, 17 Sep 2009)
@@ -30,12 +30,12 @@
 // overloaded functions to call blas
 namespace detail {
     inline void ger( const integer_t m, const integer_t n, const float alpha,
-            float const* x, const integer_t incx, float const* y,
+            const float* x, const integer_t incx, const float* y,
             const integer_t incy, float* a, const integer_t lda ) {
         BLAS_SGER( &m, &n, &alpha, x, &incx, y, &incy, a, &lda );
     }
     inline void ger( const integer_t m, const integer_t n, const double alpha,
-            double const* x, const integer_t incx, double const* y,
+            const double* x, const integer_t incx, const double* y,
             const integer_t incy, double* a, const integer_t lda ) {
         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	2009-09-17 07:18:43 EDT (Thu, 17 Sep 2009)
@@ -30,16 +30,16 @@
 // overloaded functions to call blas
 namespace detail {
     inline void gerc( const integer_t m, const integer_t n,
-            const traits::complex_f alpha, traits::complex_f const* x,
-            const integer_t incx, traits::complex_f const* y,
+            const traits::complex_f alpha, const traits::complex_f* x,
+            const integer_t incx, const traits::complex_f* y,
             const integer_t incy, traits::complex_f* a, const integer_t lda ) {
         BLAS_CGERC( &m, &n, traits::complex_ptr(&alpha),
                 traits::complex_ptr(x), &incx, traits::complex_ptr(y), &incy,
                 traits::complex_ptr(a), &lda );
     }
     inline void gerc( const integer_t m, const integer_t n,
-            const traits::complex_d alpha, traits::complex_d const* x,
-            const integer_t incx, traits::complex_d const* y,
+            const traits::complex_d alpha, const traits::complex_d* x,
+            const integer_t incx, const traits::complex_d* y,
             const integer_t incy, traits::complex_d* a, const integer_t lda ) {
         BLAS_ZGERC( &m, &n, traits::complex_ptr(&alpha),
                 traits::complex_ptr(x), &incx, traits::complex_ptr(y), &incy,
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	2009-09-17 07:18:43 EDT (Thu, 17 Sep 2009)
@@ -30,16 +30,16 @@
 // overloaded functions to call blas
 namespace detail {
     inline void geru( const integer_t m, const integer_t n,
-            const traits::complex_f alpha, traits::complex_f const* x,
-            const integer_t incx, traits::complex_f const* y,
+            const traits::complex_f alpha, const traits::complex_f* x,
+            const integer_t incx, const traits::complex_f* y,
             const integer_t incy, traits::complex_f* a, const integer_t lda ) {
         BLAS_CGERU( &m, &n, traits::complex_ptr(&alpha),
                 traits::complex_ptr(x), &incx, traits::complex_ptr(y), &incy,
                 traits::complex_ptr(a), &lda );
     }
     inline void geru( const integer_t m, const integer_t n,
-            const traits::complex_d alpha, traits::complex_d const* x,
-            const integer_t incx, traits::complex_d const* y,
+            const traits::complex_d alpha, const traits::complex_d* x,
+            const integer_t incx, const traits::complex_d* y,
             const integer_t incy, traits::complex_d* a, const integer_t lda ) {
         BLAS_ZGERU( &m, &n, traits::complex_ptr(&alpha),
                 traits::complex_ptr(x), &incx, traits::complex_ptr(y), &incy,
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	2009-09-17 07:18:43 EDT (Thu, 17 Sep 2009)
@@ -30,8 +30,8 @@
 // overloaded functions to call blas
 namespace detail {
     inline void hbmv( const char uplo, const integer_t n, const integer_t k,
-            const traits::complex_f alpha, traits::complex_f const* a,
-            const integer_t lda, traits::complex_f const* x,
+            const traits::complex_f alpha, const traits::complex_f* a,
+            const integer_t lda, const traits::complex_f* x,
             const integer_t incx, const traits::complex_f beta,
             traits::complex_f* y, const integer_t incy ) {
         BLAS_CHBMV( &uplo, &n, &k, traits::complex_ptr(&alpha),
@@ -39,8 +39,8 @@
                 traits::complex_ptr(&beta), traits::complex_ptr(y), &incy );
     }
     inline void hbmv( const char uplo, const integer_t n, const integer_t k,
-            const traits::complex_d alpha, traits::complex_d const* a,
-            const integer_t lda, traits::complex_d const* x,
+            const traits::complex_d alpha, const traits::complex_d* a,
+            const integer_t lda, const traits::complex_d* x,
             const integer_t incx, const traits::complex_d beta,
             traits::complex_d* y, const integer_t incy ) {
         BLAS_ZHBMV( &uplo, &n, &k, traits::complex_ptr(&alpha),
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	2009-09-17 07:18:43 EDT (Thu, 17 Sep 2009)
@@ -30,8 +30,8 @@
 // overloaded functions to call blas
 namespace detail {
     inline void hemv( const char uplo, const integer_t n,
-            const traits::complex_f alpha, traits::complex_f const* a,
-            const integer_t lda, traits::complex_f const* x,
+            const traits::complex_f alpha, const traits::complex_f* a,
+            const integer_t lda, const traits::complex_f* x,
             const integer_t incx, const traits::complex_f beta,
             traits::complex_f* y, const integer_t incy ) {
         BLAS_CHEMV( &uplo, &n, traits::complex_ptr(&alpha),
@@ -39,8 +39,8 @@
                 traits::complex_ptr(&beta), traits::complex_ptr(y), &incy );
     }
     inline void hemv( const char uplo, const integer_t n,
-            const traits::complex_d alpha, traits::complex_d const* a,
-            const integer_t lda, traits::complex_d const* x,
+            const traits::complex_d alpha, const traits::complex_d* a,
+            const integer_t lda, const traits::complex_d* x,
             const integer_t incx, const traits::complex_d beta,
             traits::complex_d* y, const integer_t incy ) {
         BLAS_ZHEMV( &uplo, &n, traits::complex_ptr(&alpha),
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	2009-09-17 07:18:43 EDT (Thu, 17 Sep 2009)
@@ -30,13 +30,13 @@
 // overloaded functions to call blas
 namespace detail {
     inline void her( const char uplo, const integer_t n, const float alpha,
-            traits::complex_f const* x, const integer_t incx,
+            const traits::complex_f* x, const integer_t incx,
             traits::complex_f* a, const integer_t lda ) {
         BLAS_CHER( &uplo, &n, &alpha, traits::complex_ptr(x), &incx,
                 traits::complex_ptr(a), &lda );
     }
     inline void her( const char uplo, const integer_t n, const double alpha,
-            traits::complex_d const* x, const integer_t incx,
+            const traits::complex_d* x, const integer_t incx,
             traits::complex_d* a, const integer_t lda ) {
         BLAS_ZHER( &uplo, &n, &alpha, traits::complex_ptr(x), &incx,
                 traits::complex_ptr(a), &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	2009-09-17 07:18:43 EDT (Thu, 17 Sep 2009)
@@ -30,16 +30,16 @@
 // overloaded functions to call blas
 namespace detail {
     inline void her2( const char uplo, const integer_t n,
-            const traits::complex_f alpha, traits::complex_f const* x,
-            const integer_t incx, traits::complex_f const* y,
+            const traits::complex_f alpha, const traits::complex_f* x,
+            const integer_t incx, const traits::complex_f* y,
             const integer_t incy, traits::complex_f* a, const integer_t lda ) {
         BLAS_CHER2( &uplo, &n, traits::complex_ptr(&alpha),
                 traits::complex_ptr(x), &incx, traits::complex_ptr(y), &incy,
                 traits::complex_ptr(a), &lda );
     }
     inline void her2( const char uplo, const integer_t n,
-            const traits::complex_d alpha, traits::complex_d const* x,
-            const integer_t incx, traits::complex_d const* y,
+            const traits::complex_d alpha, const traits::complex_d* x,
+            const integer_t incx, const traits::complex_d* y,
             const integer_t incy, traits::complex_d* a, const integer_t lda ) {
         BLAS_ZHER2( &uplo, &n, traits::complex_ptr(&alpha),
                 traits::complex_ptr(x), &incx, traits::complex_ptr(y), &incy,
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	2009-09-17 07:18:43 EDT (Thu, 17 Sep 2009)
@@ -30,8 +30,8 @@
 // overloaded functions to call blas
 namespace detail {
     inline void hpmv( const char uplo, const integer_t n,
-            const traits::complex_f alpha, traits::complex_f const* ap,
-            traits::complex_f const* x, const integer_t incx,
+            const traits::complex_f alpha, const traits::complex_f* ap,
+            const traits::complex_f* x, const integer_t incx,
             const traits::complex_f beta, traits::complex_f* y,
             const integer_t incy ) {
         BLAS_CHPMV( &uplo, &n, traits::complex_ptr(&alpha),
@@ -39,8 +39,8 @@
                 traits::complex_ptr(&beta), traits::complex_ptr(y), &incy );
     }
     inline void hpmv( const char uplo, const integer_t n,
-            const traits::complex_d alpha, traits::complex_d const* ap,
-            traits::complex_d const* x, const integer_t incx,
+            const traits::complex_d alpha, const traits::complex_d* ap,
+            const traits::complex_d* x, const integer_t incx,
             const traits::complex_d beta, traits::complex_d* y,
             const integer_t incy ) {
         BLAS_ZHPMV( &uplo, &n, traits::complex_ptr(&alpha),
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	2009-09-17 07:18:43 EDT (Thu, 17 Sep 2009)
@@ -30,13 +30,13 @@
 // overloaded functions to call blas
 namespace detail {
     inline void hpr( const char uplo, const integer_t n, const float alpha,
-            traits::complex_f const* x, const integer_t incx,
+            const traits::complex_f* x, const integer_t incx,
             traits::complex_f* ap ) {
         BLAS_CHPR( &uplo, &n, &alpha, traits::complex_ptr(x), &incx,
                 traits::complex_ptr(ap) );
     }
     inline void hpr( const char uplo, const integer_t n, const double alpha,
-            traits::complex_d const* x, const integer_t incx,
+            const traits::complex_d* x, const integer_t incx,
             traits::complex_d* ap ) {
         BLAS_ZHPR( &uplo, &n, &alpha, traits::complex_ptr(x), &incx,
                 traits::complex_ptr(ap) );
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	2009-09-17 07:18:43 EDT (Thu, 17 Sep 2009)
@@ -30,16 +30,16 @@
 // overloaded functions to call blas
 namespace detail {
     inline void hpr2( const char uplo, const integer_t n,
-            const traits::complex_f alpha, traits::complex_f const* x,
-            const integer_t incx, traits::complex_f const* y,
+            const traits::complex_f alpha, const traits::complex_f* x,
+            const integer_t incx, const traits::complex_f* y,
             const integer_t incy, traits::complex_f* ap ) {
         BLAS_CHPR2( &uplo, &n, traits::complex_ptr(&alpha),
                 traits::complex_ptr(x), &incx, traits::complex_ptr(y), &incy,
                 traits::complex_ptr(ap) );
     }
     inline void hpr2( const char uplo, const integer_t n,
-            const traits::complex_d alpha, traits::complex_d const* x,
-            const integer_t incx, traits::complex_d const* y,
+            const traits::complex_d alpha, const traits::complex_d* x,
+            const integer_t incx, const traits::complex_d* y,
             const integer_t incy, traits::complex_d* ap ) {
         BLAS_ZHPR2( &uplo, &n, traits::complex_ptr(&alpha),
                 traits::complex_ptr(x), &incx, traits::complex_ptr(y), &incy,
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	2009-09-17 07:18:43 EDT (Thu, 17 Sep 2009)
@@ -30,15 +30,15 @@
 // overloaded functions to call blas
 namespace detail {
     inline void sbmv( const char uplo, const integer_t n, const integer_t k,
-            const float alpha, float const* a, const integer_t lda,
-            float const* x, const integer_t incx, const float beta, float* y,
+            const float alpha, const float* a, const integer_t lda,
+            const float* x, const integer_t incx, const float beta, float* y,
             const integer_t incy ) {
         BLAS_SSBMV( &uplo, &n, &k, &alpha, a, &lda, x, &incx, &beta, y,
                 &incy );
     }
     inline void sbmv( const char uplo, const integer_t n, const integer_t k,
-            const double alpha, double const* a, const integer_t lda,
-            double const* x, const integer_t incx, const double beta,
+            const double alpha, const double* a, const integer_t lda,
+            const double* x, const integer_t incx, const double beta,
             double* y, const integer_t incy ) {
         BLAS_DSBMV( &uplo, &n, &k, &alpha, a, &lda, x, &incx, &beta, y,
                 &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	2009-09-17 07:18:43 EDT (Thu, 17 Sep 2009)
@@ -30,12 +30,12 @@
 // overloaded functions to call blas
 namespace detail {
     inline void spmv( const char uplo, const integer_t n, const float alpha,
-            float const* ap, float const* x, const integer_t incx,
+            const float* ap, const float* x, const integer_t incx,
             const float beta, float* y, const integer_t incy ) {
         BLAS_SSPMV( &uplo, &n, &alpha, ap, x, &incx, &beta, y, &incy );
     }
     inline void spmv( const char uplo, const integer_t n, const double alpha,
-            double const* ap, double const* x, const integer_t incx,
+            const double* ap, const double* x, const integer_t incx,
             const double beta, double* y, const integer_t incy ) {
         BLAS_DSPMV( &uplo, &n, &alpha, ap, x, &incx, &beta, y, &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	2009-09-17 07:18:43 EDT (Thu, 17 Sep 2009)
@@ -30,11 +30,11 @@
 // overloaded functions to call blas
 namespace detail {
     inline void spr( const char uplo, const integer_t n, const float alpha,
-            float const* x, const integer_t incx, float* ap ) {
+            const float* x, const integer_t incx, float* ap ) {
         BLAS_SSPR( &uplo, &n, &alpha, x, &incx, ap );
     }
     inline void spr( const char uplo, const integer_t n, const double alpha,
-            double const* x, const integer_t incx, double* ap ) {
+            const double* x, const integer_t incx, double* ap ) {
         BLAS_DSPR( &uplo, &n, &alpha, x, &incx, ap );
     }
 }
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	2009-09-17 07:18:43 EDT (Thu, 17 Sep 2009)
@@ -30,12 +30,12 @@
 // overloaded functions to call blas
 namespace detail {
     inline void spr2( const char uplo, const integer_t n, const float alpha,
-            float const* x, const integer_t incx, float const* y,
+            const float* x, const integer_t incx, const float* y,
             const integer_t incy, float* ap ) {
         BLAS_SSPR2( &uplo, &n, &alpha, x, &incx, y, &incy, ap );
     }
     inline void spr2( const char uplo, const integer_t n, const double alpha,
-            double const* x, const integer_t incx, double const* y,
+            const double* x, const integer_t incx, const double* y,
             const integer_t incy, double* ap ) {
         BLAS_DSPR2( &uplo, &n, &alpha, x, &incx, y, &incy, ap );
     }
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	2009-09-17 07:18:43 EDT (Thu, 17 Sep 2009)
@@ -30,13 +30,13 @@
 // overloaded functions to call blas
 namespace detail {
     inline void symv( const char uplo, const integer_t n, const float alpha,
-            float const* a, const integer_t lda, float const* x,
+            const float* a, const integer_t lda, const float* x,
             const integer_t incx, const float beta, float* y,
             const integer_t incy ) {
         BLAS_SSYMV( &uplo, &n, &alpha, a, &lda, x, &incx, &beta, y, &incy );
     }
     inline void symv( const char uplo, const integer_t n, const double alpha,
-            double const* a, const integer_t lda, double const* x,
+            const double* a, const integer_t lda, const double* x,
             const integer_t incx, const double beta, double* y,
             const integer_t incy ) {
         BLAS_DSYMV( &uplo, &n, &alpha, a, &lda, x, &incx, &beta, y, &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	2009-09-17 07:18:43 EDT (Thu, 17 Sep 2009)
@@ -30,12 +30,12 @@
 // overloaded functions to call blas
 namespace detail {
     inline void syr( const char uplo, const integer_t n, const float alpha,
-            float const* x, const integer_t incx, float* a,
+            const float* x, const integer_t incx, float* a,
             const integer_t lda ) {
         BLAS_SSYR( &uplo, &n, &alpha, x, &incx, a, &lda );
     }
     inline void syr( const char uplo, const integer_t n, const double alpha,
-            double const* x, const integer_t incx, double* a,
+            const double* x, const integer_t incx, double* a,
             const integer_t lda ) {
         BLAS_DSYR( &uplo, &n, &alpha, x, &incx, a, &lda );
     }
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	2009-09-17 07:18:43 EDT (Thu, 17 Sep 2009)
@@ -30,12 +30,12 @@
 // overloaded functions to call blas
 namespace detail {
     inline void syr2( const char uplo, const integer_t n, const float alpha,
-            float const* x, const integer_t incx, float const* y,
+            const float* x, const integer_t incx, const float* y,
             const integer_t incy, float* a, const integer_t lda ) {
         BLAS_SSYR2( &uplo, &n, &alpha, x, &incx, y, &incy, a, &lda );
     }
     inline void syr2( const char uplo, const integer_t n, const double alpha,
-            double const* x, const integer_t incx, double const* y,
+            const double* x, const integer_t incx, const double* y,
             const integer_t incy, double* a, const integer_t lda ) {
         BLAS_DSYR2( &uplo, &n, &alpha, x, &incx, y, &incy, a, &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	2009-09-17 07:18:43 EDT (Thu, 17 Sep 2009)
@@ -30,23 +30,23 @@
 // overloaded functions to call blas
 namespace detail {
     inline void tbmv( const char uplo, const char trans, const char diag,
-            const integer_t n, const integer_t k, float const* a,
+            const integer_t n, const integer_t k, const float* a,
             const integer_t lda, float* x, const integer_t incx ) {
         BLAS_STBMV( &uplo, &trans, &diag, &n, &k, a, &lda, x, &incx );
     }
     inline void tbmv( const char uplo, const char trans, const char diag,
-            const integer_t n, const integer_t k, double const* a,
+            const integer_t n, const integer_t k, const double* a,
             const integer_t lda, double* x, const integer_t incx ) {
         BLAS_DTBMV( &uplo, &trans, &diag, &n, &k, a, &lda, x, &incx );
     }
     inline void tbmv( const char uplo, const char trans, const char diag,
-            const integer_t n, const integer_t k, traits::complex_f const* a,
+            const integer_t n, const integer_t k, const traits::complex_f* a,
             const integer_t lda, traits::complex_f* x, const integer_t incx ) {
         BLAS_CTBMV( &uplo, &trans, &diag, &n, &k, traits::complex_ptr(a),
                 &lda, traits::complex_ptr(x), &incx );
     }
     inline void tbmv( const char uplo, const char trans, const char diag,
-            const integer_t n, const integer_t k, traits::complex_d const* a,
+            const integer_t n, const integer_t k, const traits::complex_d* a,
             const integer_t lda, traits::complex_d* x, const integer_t incx ) {
         BLAS_ZTBMV( &uplo, &trans, &diag, &n, &k, traits::complex_ptr(a),
                 &lda, traits::complex_ptr(x), &incx );
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	2009-09-17 07:18:43 EDT (Thu, 17 Sep 2009)
@@ -30,23 +30,23 @@
 // overloaded functions to call blas
 namespace detail {
     inline void tbsv( const char uplo, const char trans, const char diag,
-            const integer_t n, const integer_t k, float const* a,
+            const integer_t n, const integer_t k, const float* a,
             const integer_t lda, float* x, const integer_t incx ) {
         BLAS_STBSV( &uplo, &trans, &diag, &n, &k, a, &lda, x, &incx );
     }
     inline void tbsv( const char uplo, const char trans, const char diag,
-            const integer_t n, const integer_t k, double const* a,
+            const integer_t n, const integer_t k, const double* a,
             const integer_t lda, double* x, const integer_t incx ) {
         BLAS_DTBSV( &uplo, &trans, &diag, &n, &k, a, &lda, x, &incx );
     }
     inline void tbsv( const char uplo, const char trans, const char diag,
-            const integer_t n, const integer_t k, traits::complex_f const* a,
+            const integer_t n, const integer_t k, const traits::complex_f* a,
             const integer_t lda, traits::complex_f* x, const integer_t incx ) {
         BLAS_CTBSV( &uplo, &trans, &diag, &n, &k, traits::complex_ptr(a),
                 &lda, traits::complex_ptr(x), &incx );
     }
     inline void tbsv( const char uplo, const char trans, const char diag,
-            const integer_t n, const integer_t k, traits::complex_d const* a,
+            const integer_t n, const integer_t k, const traits::complex_d* a,
             const integer_t lda, traits::complex_d* x, const integer_t incx ) {
         BLAS_ZTBSV( &uplo, &trans, &diag, &n, &k, traits::complex_ptr(a),
                 &lda, traits::complex_ptr(x), &incx );
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	2009-09-17 07:18:43 EDT (Thu, 17 Sep 2009)
@@ -30,23 +30,23 @@
 // overloaded functions to call blas
 namespace detail {
     inline void tpmv( const char uplo, const char trans, const char diag,
-            const integer_t n, float const* ap, float* x,
+            const integer_t n, const float* ap, float* x,
             const integer_t incx ) {
         BLAS_STPMV( &uplo, &trans, &diag, &n, ap, x, &incx );
     }
     inline void tpmv( const char uplo, const char trans, const char diag,
-            const integer_t n, double const* ap, double* x,
+            const integer_t n, const double* ap, double* x,
             const integer_t incx ) {
         BLAS_DTPMV( &uplo, &trans, &diag, &n, ap, x, &incx );
     }
     inline void tpmv( const char uplo, const char trans, const char diag,
-            const integer_t n, traits::complex_f const* ap,
+            const integer_t n, const traits::complex_f* ap,
             traits::complex_f* x, const integer_t incx ) {
         BLAS_CTPMV( &uplo, &trans, &diag, &n, traits::complex_ptr(ap),
                 traits::complex_ptr(x), &incx );
     }
     inline void tpmv( const char uplo, const char trans, const char diag,
-            const integer_t n, traits::complex_d const* ap,
+            const integer_t n, const traits::complex_d* ap,
             traits::complex_d* x, const integer_t incx ) {
         BLAS_ZTPMV( &uplo, &trans, &diag, &n, traits::complex_ptr(ap),
                 traits::complex_ptr(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	2009-09-17 07:18:43 EDT (Thu, 17 Sep 2009)
@@ -30,23 +30,23 @@
 // overloaded functions to call blas
 namespace detail {
     inline void tpsv( const char uplo, const char trans, const char diag,
-            const integer_t n, float const* ap, float* x,
+            const integer_t n, const float* ap, float* x,
             const integer_t incx ) {
         BLAS_STPSV( &uplo, &trans, &diag, &n, ap, x, &incx );
     }
     inline void tpsv( const char uplo, const char trans, const char diag,
-            const integer_t n, double const* ap, double* x,
+            const integer_t n, const double* ap, double* x,
             const integer_t incx ) {
         BLAS_DTPSV( &uplo, &trans, &diag, &n, ap, x, &incx );
     }
     inline void tpsv( const char uplo, const char trans, const char diag,
-            const integer_t n, traits::complex_f const* ap,
+            const integer_t n, const traits::complex_f* ap,
             traits::complex_f* x, const integer_t incx ) {
         BLAS_CTPSV( &uplo, &trans, &diag, &n, traits::complex_ptr(ap),
                 traits::complex_ptr(x), &incx );
     }
     inline void tpsv( const char uplo, const char trans, const char diag,
-            const integer_t n, traits::complex_d const* ap,
+            const integer_t n, const traits::complex_d* ap,
             traits::complex_d* x, const integer_t incx ) {
         BLAS_ZTPSV( &uplo, &trans, &diag, &n, traits::complex_ptr(ap),
                 traits::complex_ptr(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	2009-09-17 07:18:43 EDT (Thu, 17 Sep 2009)
@@ -30,23 +30,23 @@
 // overloaded functions to call blas
 namespace detail {
     inline void trmv( const char uplo, const char trans, const char diag,
-            const integer_t n, float const* a, const integer_t lda, float* x,
+            const integer_t n, const float* a, const integer_t lda, float* x,
             const integer_t incx ) {
         BLAS_STRMV( &uplo, &trans, &diag, &n, a, &lda, x, &incx );
     }
     inline void trmv( const char uplo, const char trans, const char diag,
-            const integer_t n, double const* a, const integer_t lda,
+            const integer_t n, const double* a, const integer_t lda,
             double* x, const integer_t incx ) {
         BLAS_DTRMV( &uplo, &trans, &diag, &n, a, &lda, x, &incx );
     }
     inline void trmv( const char uplo, const char trans, const char diag,
-            const integer_t n, traits::complex_f const* a,
+            const integer_t n, const traits::complex_f* a,
             const integer_t lda, traits::complex_f* x, const integer_t incx ) {
         BLAS_CTRMV( &uplo, &trans, &diag, &n, traits::complex_ptr(a), &lda,
                 traits::complex_ptr(x), &incx );
     }
     inline void trmv( const char uplo, const char trans, const char diag,
-            const integer_t n, traits::complex_d const* a,
+            const integer_t n, const traits::complex_d* a,
             const integer_t lda, traits::complex_d* x, const integer_t incx ) {
         BLAS_ZTRMV( &uplo, &trans, &diag, &n, traits::complex_ptr(a), &lda,
                 traits::complex_ptr(x), &incx );
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	2009-09-17 07:18:43 EDT (Thu, 17 Sep 2009)
@@ -30,23 +30,23 @@
 // overloaded functions to call blas
 namespace detail {
     inline void trsv( const char uplo, const char trans, const char diag,
-            const integer_t n, float const* a, const integer_t lda, float* x,
+            const integer_t n, const float* a, const integer_t lda, float* x,
             const integer_t incx ) {
         BLAS_STRSV( &uplo, &trans, &diag, &n, a, &lda, x, &incx );
     }
     inline void trsv( const char uplo, const char trans, const char diag,
-            const integer_t n, double const* a, const integer_t lda,
+            const integer_t n, const double* a, const integer_t lda,
             double* x, const integer_t incx ) {
         BLAS_DTRSV( &uplo, &trans, &diag, &n, a, &lda, x, &incx );
     }
     inline void trsv( const char uplo, const char trans, const char diag,
-            const integer_t n, traits::complex_f const* a,
+            const integer_t n, const traits::complex_f* a,
             const integer_t lda, traits::complex_f* x, const integer_t incx ) {
         BLAS_CTRSV( &uplo, &trans, &diag, &n, traits::complex_ptr(a), &lda,
                 traits::complex_ptr(x), &incx );
     }
     inline void trsv( const char uplo, const char trans, const char diag,
-            const integer_t n, traits::complex_d const* a,
+            const integer_t n, const traits::complex_d* a,
             const integer_t lda, traits::complex_d* x, const integer_t incx ) {
         BLAS_ZTRSV( &uplo, &trans, &diag, &n, traits::complex_ptr(a), &lda,
                 traits::complex_ptr(x), &incx );
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	2009-09-17 07:18:43 EDT (Thu, 17 Sep 2009)
@@ -31,7 +31,7 @@
 namespace detail {
     inline void gemm( const char transa, const char transb, const integer_t m,
             const integer_t n, const integer_t k, const float alpha,
-            float const* a, const integer_t lda, float const* b,
+            const float* a, const integer_t lda, const float* b,
             const integer_t ldb, const float beta, float* c,
             const integer_t ldc ) {
         BLAS_SGEMM( &transa, &transb, &m, &n, &k, &alpha, a, &lda, b, &ldb,
@@ -39,7 +39,7 @@
     }
     inline void gemm( const char transa, const char transb, const integer_t m,
             const integer_t n, const integer_t k, const double alpha,
-            double const* a, const integer_t lda, double const* b,
+            const double* a, const integer_t lda, const double* b,
             const integer_t ldb, const double beta, double* c,
             const integer_t ldc ) {
         BLAS_DGEMM( &transa, &transb, &m, &n, &k, &alpha, a, &lda, b, &ldb,
@@ -47,8 +47,8 @@
     }
     inline void gemm( const char transa, const char transb, const integer_t m,
             const integer_t n, const integer_t k,
-            const traits::complex_f alpha, traits::complex_f const* a,
-            const integer_t lda, traits::complex_f const* b,
+            const traits::complex_f alpha, const traits::complex_f* a,
+            const integer_t lda, const traits::complex_f* b,
             const integer_t ldb, const traits::complex_f beta,
             traits::complex_f* c, const integer_t ldc ) {
         BLAS_CGEMM( &transa, &transb, &m, &n, &k, traits::complex_ptr(&alpha),
@@ -57,8 +57,8 @@
     }
     inline void gemm( const char transa, const char transb, const integer_t m,
             const integer_t n, const integer_t k,
-            const traits::complex_d alpha, traits::complex_d const* a,
-            const integer_t lda, traits::complex_d const* b,
+            const traits::complex_d alpha, const traits::complex_d* a,
+            const integer_t lda, const traits::complex_d* b,
             const integer_t ldb, const traits::complex_d beta,
             traits::complex_d* c, const integer_t ldc ) {
         BLAS_ZGEMM( &transa, &transb, &m, &n, &k, traits::complex_ptr(&alpha),
@@ -78,19 +78,20 @@
     // templated specialization
     template< typename MatrixA, typename MatrixB, typename MatrixC >
     static return_type invoke( const char transa, const char transb,
-            const integer_t k, const value_type alpha, const MatrixA& a,
-            const MatrixB& b, const value_type beta, MatrixC& c ) {
+            const value_type alpha, const MatrixA& a, const MatrixB& b,
+            const value_type beta, MatrixC& c ) {
         BOOST_STATIC_ASSERT( (boost::is_same< typename traits::matrix_traits<
                 MatrixA >::value_type, typename traits::matrix_traits<
                 MatrixB >::value_type >::value) );
         BOOST_STATIC_ASSERT( (boost::is_same< typename traits::matrix_traits<
                 MatrixA >::value_type, typename traits::matrix_traits<
                 MatrixC >::value_type >::value) );
-        detail::gemm( transa, transb, traits::matrix_num_rows(c),
-                traits::matrix_num_columns(c), k, alpha,
-                traits::matrix_storage(a), traits::leading_dimension(a),
-                traits::matrix_storage(b), traits::leading_dimension(b), beta,
-                traits::matrix_storage(c), traits::leading_dimension(c) );
+        detail::gemm( transa, transb, , traits::matrix_num_columns(c),
+                (transa=='N'?traits::matrix_num_columns(a),
+                traits::matrix_num_rows(a)), alpha, traits::matrix_storage(a),
+                traits::leading_dimension(a), traits::matrix_storage(b),
+                traits::leading_dimension(b), beta, traits::matrix_storage(c),
+                traits::leading_dimension(c) );
     }
 };
 
@@ -98,14 +99,14 @@
 template< typename MatrixA, typename MatrixB, typename MatrixC >
 inline typename gemm_impl< typename traits::matrix_traits<
         MatrixA >::value_type >::return_type
-gemm( const char transa, const char transb, const integer_t k,
+gemm( const char transa, const char transb,
         const typename traits::matrix_traits< MatrixA >::value_type alpha,
         const MatrixA& a, const MatrixB& b,
         const typename traits::matrix_traits< MatrixA >::value_type beta,
         MatrixC& c ) {
     typedef typename traits::matrix_traits< MatrixA >::value_type value_type;
-    gemm_impl< value_type >::invoke( transa, transb, k, alpha, a, b,
-            beta, c );
+    gemm_impl< value_type >::invoke( transa, transb, alpha, a, b, beta,
+            c );
 }
 
 }}}}} // namespace boost::numeric::bindings::blas::level3
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	2009-09-17 07:18:43 EDT (Thu, 17 Sep 2009)
@@ -31,8 +31,8 @@
 namespace detail {
     inline void hemm( const char side, const char uplo, const integer_t m,
             const integer_t n, const traits::complex_f alpha,
-            traits::complex_f const* a, const integer_t lda,
-            traits::complex_f const* b, const integer_t ldb,
+            const traits::complex_f* a, const integer_t lda,
+            const traits::complex_f* b, const integer_t ldb,
             const traits::complex_f beta, traits::complex_f* c,
             const integer_t ldc ) {
         BLAS_CHEMM( &side, &uplo, &m, &n, traits::complex_ptr(&alpha),
@@ -41,8 +41,8 @@
     }
     inline void hemm( const char side, const char uplo, const integer_t m,
             const integer_t n, const traits::complex_d alpha,
-            traits::complex_d const* a, const integer_t lda,
-            traits::complex_d const* b, const integer_t ldb,
+            const traits::complex_d* a, const integer_t lda,
+            const traits::complex_d* b, const integer_t ldb,
             const traits::complex_d beta, traits::complex_d* c,
             const integer_t ldc ) {
         BLAS_ZHEMM( &side, &uplo, &m, &n, traits::complex_ptr(&alpha),
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	2009-09-17 07:18:43 EDT (Thu, 17 Sep 2009)
@@ -31,8 +31,8 @@
 namespace detail {
     inline void her2k( const char uplo, const char trans, const integer_t n,
             const integer_t k, const traits::complex_f alpha,
-            traits::complex_f const* a, const integer_t lda,
-            traits::complex_f const* b, const integer_t ldb, const float beta,
+            const traits::complex_f* a, const integer_t lda,
+            const traits::complex_f* b, const integer_t ldb, const float beta,
             traits::complex_f* c, const integer_t ldc ) {
         BLAS_CHER2K( &uplo, &trans, &n, &k, traits::complex_ptr(&alpha),
                 traits::complex_ptr(a), &lda, traits::complex_ptr(b), &ldb,
@@ -40,8 +40,8 @@
     }
     inline void her2k( const char uplo, const char trans, const integer_t n,
             const integer_t k, const traits::complex_d alpha,
-            traits::complex_d const* a, const integer_t lda,
-            traits::complex_d const* b, const integer_t ldb,
+            const traits::complex_d* a, const integer_t lda,
+            const traits::complex_d* b, const integer_t ldb,
             const double beta, traits::complex_d* c, const integer_t ldc ) {
         BLAS_ZHER2K( &uplo, &trans, &n, &k, traits::complex_ptr(&alpha),
                 traits::complex_ptr(a), &lda, traits::complex_ptr(b), &ldb,
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	2009-09-17 07:18:43 EDT (Thu, 17 Sep 2009)
@@ -30,14 +30,14 @@
 // overloaded functions to call blas
 namespace detail {
     inline void herk( const char uplo, const char trans, const integer_t n,
-            const integer_t k, const float alpha, traits::complex_f const* a,
+            const integer_t k, const float alpha, const traits::complex_f* a,
             const integer_t lda, const float beta, traits::complex_f* c,
             const integer_t ldc ) {
         BLAS_CHERK( &uplo, &trans, &n, &k, &alpha, traits::complex_ptr(a),
                 &lda, &beta, traits::complex_ptr(c), &ldc );
     }
     inline void herk( const char uplo, const char trans, const integer_t n,
-            const integer_t k, const double alpha, traits::complex_d const* a,
+            const integer_t k, const double alpha, const traits::complex_d* a,
             const integer_t lda, const double beta, traits::complex_d* c,
             const integer_t ldc ) {
         BLAS_ZHERK( &uplo, &trans, &n, &k, &alpha, traits::complex_ptr(a),
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	2009-09-17 07:18:43 EDT (Thu, 17 Sep 2009)
@@ -30,23 +30,23 @@
 // overloaded functions to call blas
 namespace detail {
     inline void symm( const char side, const char uplo, const integer_t m,
-            const integer_t n, const float alpha, float const* a,
-            const integer_t lda, float const* b, const integer_t ldb,
+            const integer_t n, const float alpha, const float* a,
+            const integer_t lda, const float* b, const integer_t ldb,
             const float beta, float* c, const integer_t ldc ) {
         BLAS_SSYMM( &side, &uplo, &m, &n, &alpha, a, &lda, b, &ldb, &beta, c,
                 &ldc );
     }
     inline void symm( const char side, const char uplo, const integer_t m,
-            const integer_t n, const double alpha, double const* a,
-            const integer_t lda, double const* b, const integer_t ldb,
+            const integer_t n, const double alpha, const double* a,
+            const integer_t lda, const double* b, const integer_t ldb,
             const double beta, double* c, const integer_t ldc ) {
         BLAS_DSYMM( &side, &uplo, &m, &n, &alpha, a, &lda, b, &ldb, &beta, c,
                 &ldc );
     }
     inline void symm( const char side, const char uplo, const integer_t m,
             const integer_t n, const traits::complex_f alpha,
-            traits::complex_f const* a, const integer_t lda,
-            traits::complex_f const* b, const integer_t ldb,
+            const traits::complex_f* a, const integer_t lda,
+            const traits::complex_f* b, const integer_t ldb,
             const traits::complex_f beta, traits::complex_f* c,
             const integer_t ldc ) {
         BLAS_CSYMM( &side, &uplo, &m, &n, traits::complex_ptr(&alpha),
@@ -55,8 +55,8 @@
     }
     inline void symm( const char side, const char uplo, const integer_t m,
             const integer_t n, const traits::complex_d alpha,
-            traits::complex_d const* a, const integer_t lda,
-            traits::complex_d const* b, const integer_t ldb,
+            const traits::complex_d* a, const integer_t lda,
+            const traits::complex_d* b, const integer_t ldb,
             const traits::complex_d beta, traits::complex_d* c,
             const integer_t ldc ) {
         BLAS_ZSYMM( &side, &uplo, &m, &n, traits::complex_ptr(&alpha),
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	2009-09-17 07:18:43 EDT (Thu, 17 Sep 2009)
@@ -30,23 +30,23 @@
 // overloaded functions to call blas
 namespace detail {
     inline void syr2k( const char uplo, const char trans, const integer_t n,
-            const integer_t k, const float alpha, float const* a,
-            const integer_t lda, float const* b, const integer_t ldb,
+            const integer_t k, const float alpha, const float* a,
+            const integer_t lda, const float* b, const integer_t ldb,
             const float beta, float* c, const integer_t ldc ) {
         BLAS_SSYR2K( &uplo, &trans, &n, &k, &alpha, a, &lda, b, &ldb, &beta,
                 c, &ldc );
     }
     inline void syr2k( const char uplo, const char trans, const integer_t n,
-            const integer_t k, const double alpha, double const* a,
-            const integer_t lda, double const* b, const integer_t ldb,
+            const integer_t k, const double alpha, const double* a,
+            const integer_t lda, const double* b, const integer_t ldb,
             const double beta, double* c, const integer_t ldc ) {
         BLAS_DSYR2K( &uplo, &trans, &n, &k, &alpha, a, &lda, b, &ldb, &beta,
                 c, &ldc );
     }
     inline void syr2k( const char uplo, const char trans, const integer_t n,
             const integer_t k, const traits::complex_f alpha,
-            traits::complex_f const* a, const integer_t lda,
-            traits::complex_f const* b, const integer_t ldb,
+            const traits::complex_f* a, const integer_t lda,
+            const traits::complex_f* b, const integer_t ldb,
             const traits::complex_f beta, traits::complex_f* c,
             const integer_t ldc ) {
         BLAS_CSYR2K( &uplo, &trans, &n, &k, traits::complex_ptr(&alpha),
@@ -55,8 +55,8 @@
     }
     inline void syr2k( const char uplo, const char trans, const integer_t n,
             const integer_t k, const traits::complex_d alpha,
-            traits::complex_d const* a, const integer_t lda,
-            traits::complex_d const* b, const integer_t ldb,
+            const traits::complex_d* a, const integer_t lda,
+            const traits::complex_d* b, const integer_t ldb,
             const traits::complex_d beta, traits::complex_d* c,
             const integer_t ldc ) {
         BLAS_ZSYR2K( &uplo, &trans, &n, &k, traits::complex_ptr(&alpha),
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	2009-09-17 07:18:43 EDT (Thu, 17 Sep 2009)
@@ -30,20 +30,20 @@
 // overloaded functions to call blas
 namespace detail {
     inline void syrk( const char uplo, const char trans, const integer_t n,
-            const integer_t k, const float alpha, float const* a,
+            const integer_t k, const float alpha, const float* a,
             const integer_t lda, const float beta, float* c,
             const integer_t ldc ) {
         BLAS_SSYRK( &uplo, &trans, &n, &k, &alpha, a, &lda, &beta, c, &ldc );
     }
     inline void syrk( const char uplo, const char trans, const integer_t n,
-            const integer_t k, const double alpha, double const* a,
+            const integer_t k, const double alpha, const double* a,
             const integer_t lda, const double beta, double* c,
             const integer_t ldc ) {
         BLAS_DSYRK( &uplo, &trans, &n, &k, &alpha, a, &lda, &beta, c, &ldc );
     }
     inline void syrk( const char uplo, const char trans, const integer_t n,
             const integer_t k, const traits::complex_f alpha,
-            traits::complex_f const* a, const integer_t lda,
+            const traits::complex_f* a, const integer_t lda,
             const traits::complex_f beta, traits::complex_f* c,
             const integer_t ldc ) {
         BLAS_CSYRK( &uplo, &trans, &n, &k, traits::complex_ptr(&alpha),
@@ -52,7 +52,7 @@
     }
     inline void syrk( const char uplo, const char trans, const integer_t n,
             const integer_t k, const traits::complex_d alpha,
-            traits::complex_d const* a, const integer_t lda,
+            const traits::complex_d* a, const integer_t lda,
             const traits::complex_d beta, traits::complex_d* c,
             const integer_t ldc ) {
         BLAS_ZSYRK( &uplo, &trans, &n, &k, traits::complex_ptr(&alpha),
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	2009-09-17 07:18:43 EDT (Thu, 17 Sep 2009)
@@ -31,21 +31,21 @@
 namespace detail {
     inline void trmm( const char side, const char uplo, const char transa,
             const char diag, const integer_t m, const integer_t n,
-            const float alpha, float const* a, const integer_t lda, float* b,
+            const float alpha, const float* a, const integer_t lda, float* b,
             const integer_t ldb ) {
         BLAS_STRMM( &side, &uplo, &transa, &diag, &m, &n, &alpha, a, &lda, b,
                 &ldb );
     }
     inline void trmm( const char side, const char uplo, const char transa,
             const char diag, const integer_t m, const integer_t n,
-            const double alpha, double const* a, const integer_t lda,
+            const double alpha, const double* a, const integer_t lda,
             double* b, const integer_t ldb ) {
         BLAS_DTRMM( &side, &uplo, &transa, &diag, &m, &n, &alpha, a, &lda, b,
                 &ldb );
     }
     inline void trmm( const char side, const char uplo, const char transa,
             const char diag, const integer_t m, const integer_t n,
-            const traits::complex_f alpha, traits::complex_f const* a,
+            const traits::complex_f alpha, const traits::complex_f* a,
             const integer_t lda, traits::complex_f* b, const integer_t ldb ) {
         BLAS_CTRMM( &side, &uplo, &transa, &diag, &m, &n,
                 traits::complex_ptr(&alpha), traits::complex_ptr(a), &lda,
@@ -53,7 +53,7 @@
     }
     inline void trmm( const char side, const char uplo, const char transa,
             const char diag, const integer_t m, const integer_t n,
-            const traits::complex_d alpha, traits::complex_d const* a,
+            const traits::complex_d alpha, const traits::complex_d* a,
             const integer_t lda, traits::complex_d* b, const integer_t ldb ) {
         BLAS_ZTRMM( &side, &uplo, &transa, &diag, &m, &n,
                 traits::complex_ptr(&alpha), traits::complex_ptr(a), &lda,
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	2009-09-17 07:18:43 EDT (Thu, 17 Sep 2009)
@@ -31,21 +31,21 @@
 namespace detail {
     inline void trsm( const char side, const char uplo, const char transa,
             const char diag, const integer_t m, const integer_t n,
-            const float alpha, float const* a, const integer_t lda, float* b,
+            const float alpha, const float* a, const integer_t lda, float* b,
             const integer_t ldb ) {
         BLAS_STRSM( &side, &uplo, &transa, &diag, &m, &n, &alpha, a, &lda, b,
                 &ldb );
     }
     inline void trsm( const char side, const char uplo, const char transa,
             const char diag, const integer_t m, const integer_t n,
-            const double alpha, double const* a, const integer_t lda,
+            const double alpha, const double* a, const integer_t lda,
             double* b, const integer_t ldb ) {
         BLAS_DTRSM( &side, &uplo, &transa, &diag, &m, &n, &alpha, a, &lda, b,
                 &ldb );
     }
     inline void trsm( const char side, const char uplo, const char transa,
             const char diag, const integer_t m, const integer_t n,
-            const traits::complex_f alpha, traits::complex_f const* a,
+            const traits::complex_f alpha, const traits::complex_f* a,
             const integer_t lda, traits::complex_f* b, const integer_t ldb ) {
         BLAS_CTRSM( &side, &uplo, &transa, &diag, &m, &n,
                 traits::complex_ptr(&alpha), traits::complex_ptr(a), &lda,
@@ -53,7 +53,7 @@
     }
     inline void trsm( const char side, const char uplo, const char transa,
             const char diag, const integer_t m, const integer_t n,
-            const traits::complex_d alpha, traits::complex_d const* a,
+            const traits::complex_d alpha, const traits::complex_d* a,
             const integer_t lda, traits::complex_d* b, const integer_t ldb ) {
         BLAS_ZTRSM( &side, &uplo, &transa, &diag, &m, &n,
                 traits::complex_ptr(&alpha), traits::complex_ptr(a), &lda,