[Boost-commit] svn:boost r48263 - sandbox/boost/numeric/bindings/lapack

From: karl.meerbergen_at_[hidden]
Date: 2008-08-20 10:27:03

Author: karlmeerbergen
Date: 2008-08-20 10:27:03 EDT (Wed, 20 Aug 2008)
New Revision: 48263
URL: http://svn.boost.org/trac/boost/changeset/48263

Log:
bug in manipulation of array w

Text files modified: 
   sandbox/boost/numeric/bindings/lapack/geev.hpp |     6 +++---                                  
   1 files changed, 3 insertions(+), 3 deletions(-)

Modified: sandbox/boost/numeric/bindings/lapack/geev.hpp
==============================================================================
--- sandbox/boost/numeric/bindings/lapack/geev.hpp	(original)
+++ sandbox/boost/numeric/bindings/lapack/geev.hpp	2008-08-20 10:27:03 EDT (Wed, 20 Aug 2008)
@@ -171,7 +171,7 @@
                               work.storage(), &lwork);
 
         for (int i = 0; i < n; i++)
-	  w[i] = std::complex<value_type>(wr[i], wi[i]);
+	  traits::vector_storage(w)[i] = std::complex<value_type>(wr[i], wi[i]);
         return result;
       }
 
@@ -226,14 +226,14 @@
         
         for (int i = 0; i < n; i++)
         {
-	  w[i] = std::complex<value_type>(wr[i], wi[i]);
+          traits::vector_storage(w)[i] = std::complex<value_type>(wr[i], wi[i]);
           if (wi[i] != 0)
           {
             assert(i+1 < n);
             assert(wr[i+1] == wr[i]);
             assert(wi[i+1] == -wi[i]);
 
-	    w[i+1] = std::complex<value_type>(wr[i+1], wi[i+1]);
+            traits::vector_storage(w)[i+1] = std::complex<value_type>(wr[i+1], wi[i+1]);
             for (int j = 0; j < n; j++)
             {
               if (vl)