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

Subject: [Boost-commit] svn:boost r58928 - sandbox/numeric_bindings/libs/numeric/bindings/lapack/test
From: rutger_at_[hidden]
Date: 2010-01-12 06:09:52

Author: rutger
Date: 2010-01-12 06:09:52 EST (Tue, 12 Jan 2010)
New Revision: 58928
URL: http://svn.boost.org/trac/boost/changeset/58928

Log:
ublas_getrf_getrs is now passing regression

Text files modified: 
   sandbox/numeric_bindings/libs/numeric/bindings/lapack/test/ublas_getrf_getrs.cc |     9 +++++----                               
   1 files changed, 5 insertions(+), 4 deletions(-)

Modified: sandbox/numeric_bindings/libs/numeric/bindings/lapack/test/ublas_getrf_getrs.cc
==============================================================================
--- sandbox/numeric_bindings/libs/numeric/bindings/lapack/test/ublas_getrf_getrs.cc	(original)
+++ sandbox/numeric_bindings/libs/numeric/bindings/lapack/test/ublas_getrf_getrs.cc	2010-01-12 06:09:52 EST (Tue, 12 Jan 2010)
@@ -9,12 +9,13 @@
 #include <boost/numeric/bindings/lapack/computational/getri.hpp>
 #include <boost/numeric/bindings/lapack/computational/getrs.hpp>
 #include <boost/numeric/bindings/ublas/matrix.hpp>
+#include <boost/numeric/bindings/ublas/matrix_proxy.hpp>
 #include <boost/numeric/bindings/std/vector.hpp>
-#include <boost/numeric/ublas/matrix_proxy.hpp>
+#include <boost/numeric/bindings/trans.hpp>
 #include "utils.h"
 
 namespace ublas = boost::numeric::ublas;
-namespace traits = boost::numeric::bindings::traits;
+namespace bindings = boost::numeric::bindings;
 namespace lapack = boost::numeric::bindings::lapack;
 
 using std::size_t; 
@@ -82,7 +83,7 @@
 
   lapack::getrf (a, ipiv);      // factor a
   m_t ia (a);
-  lapack::getrs ('N', a, ipiv, b);   // solve from factorization 
+  lapack::getrs ( a, ipiv, b);   // solve from factorization 
   print_m (b, "X"); 
   cout << endl; 
   lapack::getri (ia, ipiv);     // invert a
@@ -104,7 +105,7 @@
   
   lapack::getrf (a2, ipiv); // factor a
   m_t ia2 (a2);
-  lapack::getrs ('T', a2, ipiv, b2); // solve 
+  lapack::getrs ( bindings::trans( a2 ), ipiv, b2); // solve 
   print_m (b2, "X"); 
   cout << endl; 
   lapack::getri (ia2, ipiv); // invert a2