Karl's geev fix.

inducer · Aug 21, 2008 · b90458e · b90458e
1 parent e72b237
commit b90458e
Showing 1 changed file with 5 additions and 3 deletions.
diff --git a/boost/numeric/bindings/lapack/geev.hpp b/boost/numeric/bindings/lapack/geev.hpp
@@ -171,7 +171,7 @@ namespace boost { namespace numeric { namespace bindings {
 			      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,16 @@ namespace boost { namespace numeric { namespace bindings {
 
 	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)