Changeset 63

Timestamp:

Apr 19, 2004, 1:21:28 PM (17 years ago)

Author:

Peter

Message:

minmax function added

Location:

trunk/src

Files:

• vector.cc (modified) (2 diffs)
• vector.h (modified) (2 diffs)

trunk/src/vector.cc

@@ -5 +5 @@
 #include <string>
 #include <vector>
+#include <utility>
 
 #include "vector.h"
@@ -125 +126 @@
   }
 
+//Peter, add default checking all elements using GSL
+  std::pair<u_int,u_int> vector::minmax_index(const std::vector<u_int>& subset ) const
+  {
+    u_int min_index = subset[0];
+    u_int max_index = subset[0];
+    for (unsigned int i=1; i<subset.size(); i++) {
+      if (gsl_vector_get( v_, subset[i]) < gsl_vector_get( v_, min_index))
+        min_index = subset[i];
+      else if (gsl_vector_get( v_, subset[i]) > gsl_vector_get( v_, max_index))
+        max_index = subset[i];
+    }
+    return std::pair<u_int,u_int>(min_index, max_index);
+  }
+
 
 

trunk/src/vector.h

@@ -5 +5 @@
 
 #include <iostream>
+#include <vector>
+#include <utility>
 
 #include <gsl/gsl_vector.h>
@@ -86 +88 @@
     /// @return A pointer to the internal GSL vector,
     ///
-    inline gsl_vector* gsl_vector_pointer(void) { return v_; };
+    inline gsl_vector* gsl_vector_pointer(void) { return v_; }
+
+    ///This function returns the indices of the minimum and maximum values in
+    ///the sub-vector (defined by subset), storing them in imin and imax. When
+    ///there are several equal minimum or maximum elements then the lowest
+    ///indices are returned. The returned index is the index from the complete vector (not the sub-vector)  @return Index corresponding to the smallest
+    ///and largest value.
+    ///
+    std::pair<u_int,u_int> vector::minmax_index(const std::vector<u_int>& subset ) const;
 
     ///