No longer sorts indices before selecting vols - naughty..

fslselectvols.cc

@@ -55,10 +55,10 @@ int get_vols(Matrix& id_vols,const string& vols){
     id_vols = read_ascii_matrix(vols);
     if(id_vols.Ncols()==1){id_vols=id_vols.t();}
 
-    comps.resize(id_vols.Ncols());
-    for(int i=0;i<(int)comps.size();i++){
-      comps[i]=id_vols(1,i+1);
-    }
+    //comps.resize(id_vols.Ncols());
+    //for(int i=0;i<(int)comps.size();i++){
+    //  comps[i]=id_vols(1,i+1);
+    //}
   }
   else{
     if(vols.length()>0 && parse_filterstring(comps,vols)){
@@ -67,14 +67,14 @@ int get_vols(Matrix& id_vols,const string& vols){
   }
 
   //sort and remove duplicates 
-  sort (comps.begin(), comps.end());
-  vector<int>::iterator it = unique (comps.begin(), comps.end());
-  comps.resize( it - comps.begin() );
+  //sort (comps.begin(), comps.end());
+  //vector<int>::iterator it = unique (comps.begin(), comps.end());
+  //comps.resize( it - comps.begin() );
     
-  id_vols.ReSize(1,comps.size());
-  for(int i=0;i<(int)comps.size();i++){
-    id_vols(1,i+1)=comps[i];
-  }
+  //id_vols.ReSize(1,comps.size());
+  //for(int i=0;i<(int)comps.size();i++){
+  //  id_vols(1,i+1)=comps[i];
+  //}
   
   
   return 0;