Changeset 877 for trunk/yat/utility/vector.cc

Timestamp:

Sep 20, 2007, 9:01:57 PM (16 years ago)

Author:

Markus Ringnér

Message:

This revision breaks the vector_test. I know this is against the rules but it is intented for getting help from Jari to identify linking problem. A new function sort_index is added to vector, and tested in vector_test. Everything compiles but vector_test fails to link and I cannot see why

File:

• trunk/yat/utility/vector.cc (modified) (1 diff)

trunk/yat/utility/vector.cc

@@ -513 +513 @@
 
 
+  void sort_index(std::vector<size_t>& sort_index, const vector& invec)
+  {
+    assert(invec.gsl_vector_p());
+    gsl_permutation* p = gsl_permutation_alloc(invec.size());
+    int status=gsl_sort_vector_index(p,invec.gsl_vector_p());
+    if (status) {
+      gsl_permutation_free(p);
+      throw utility::GSL_error(std::string("sort_index(vector&,const vector&)",status));      
+    }
+    for(size_t i=0;i<p->size;i++) {
+      sort_index.resize(0);
+      sort_index.push_back(gsl_permutation_get(p,i));
+    }
+    gsl_permutation_free(p);
+  }
+
+
+
   double sum(const vector& v)
   {