[Boost-commit] svn:boost r58960 - sandbox/numeric_bindings/libs/numeric/bindings/lapack/test

Subject: [Boost-commit] svn:boost r58960 - sandbox/numeric_bindings/libs/numeric/bindings/lapack/test
From: rutger_at_[hidden]
Date: 2010-01-13 05:04:22

Author: rutger
Date: 2010-01-13 05:04:21 EST (Wed, 13 Jan 2010)
New Revision: 58960
URL: http://svn.boost.org/trac/boost/changeset/58960

Log:
ublas_geqrf regression compiles now, but doesn't pass yet

Text files modified: 
   sandbox/numeric_bindings/libs/numeric/bindings/lapack/test/ublas_geqrf.cpp |    26 ++++++++++++--------------              
   1 files changed, 12 insertions(+), 14 deletions(-)

Modified: sandbox/numeric_bindings/libs/numeric/bindings/lapack/test/ublas_geqrf.cpp
==============================================================================
--- sandbox/numeric_bindings/libs/numeric/bindings/lapack/test/ublas_geqrf.cpp	(original)
+++ sandbox/numeric_bindings/libs/numeric/bindings/lapack/test/ublas_geqrf.cpp	2010-01-13 05:04:21 EST (Wed, 13 Jan 2010)
@@ -32,31 +32,29 @@
 
 struct apply_real {
   template< typename MatrixA, typename VectorTAU, typename Workspace >
-  static inline integer_t orgqr( const integer_t m, const integer_t n,
-        const integer_t k, MatrixA& a, const VectorTAU& tau, Workspace work ) {
-    return lapack::orgqr( m, n, k, a, tau, work );
+  static inline integer_t orgqr(  MatrixA& a, const VectorTAU& tau, Workspace work ) {
+    return lapack::orgqr( a, tau, work );
   }
   template< typename MatrixA, typename VectorTAU, typename MatrixC,
         typename Workspace >
-  static inline integer_t ormqr( const char side, const char trans,
-        const integer_t k, const MatrixA& a, const VectorTAU& tau, MatrixC& c,
+  static inline integer_t ormqr( const char side, 
+        const MatrixA& a, const VectorTAU& tau, MatrixC& c,
         Workspace work ) {
-    return lapack::ormqr( side, trans, k, a, tau, c, work );
+    return lapack::ormqr( side, a, tau, c, work );
   }
 };
 
 struct apply_complex {
   template< typename MatrixA, typename VectorTAU, typename Workspace >
-  static inline integer_t orgqr( const integer_t m, const integer_t n,
-        const integer_t k, MatrixA& a, const VectorTAU& tau, Workspace work ) {
-    return lapack::ungqr( m, n, k, a, tau, work );
+  static inline integer_t orgqr( MatrixA& a, const VectorTAU& tau, Workspace work ) {
+    return lapack::ungqr( a, tau, work );
   }
   template< typename MatrixA, typename VectorTAU, typename MatrixC,
         typename Workspace >
-  static inline integer_t ormqr( const char side, const char trans,
-        const integer_t k, const MatrixA& a, const VectorTAU& tau, MatrixC& c,
+  static inline integer_t ormqr( const char side, 
+        const MatrixA& a, const VectorTAU& tau, MatrixC& c,
         Workspace work ) {
-    return lapack::unmqr( side, trans, k, a, tau, c, work );
+    return lapack::unmqr( side, a, tau, c, work );
   }
 };
 
@@ -122,14 +120,14 @@
    lapack::geqrf( a, tau, workspace ) ;
 
    // Apply the orthogonal transformations to a2
-   apply_t::ormqr( 'L', transpose<T>::value, bindings::size (tau), a, tau, a2, workspace );
+   apply_t::ormqr( 'L', a, tau, a2, workspace );
 
    // The upper triangular parts of a and a2 must be equal.
    if (norm_frobenius( upper_part( a - a2 ) )
             > std::numeric_limits<real_type>::epsilon() * 10.0 * norm_frobenius( upper_part( a ) ) ) return 255 ;
 
    // Generate orthogonal matrix
-   apply_t::orgqr( bindings::size_row (a), bindings::size_column (a), bindings::size (tau), a, tau, workspace );
+   apply_t::orgqr( a, tau, workspace );
 
    // The result of lapack::ormqr and the equivalent matrix product must be equal.
    if (norm_frobenius( a2 - prod(herm(a), a3) )