# HG changeset patch
# User Balazs Dezso <deba@inf.elte.hu>
# Date 1206029735 -3600
# Node ID 81563e019fa4cd5ce484d3eecc7e30da2e23fa30
# Parent  70f3967ca6eb8c62d5194a2ad078dd36d9a36374
Seeding random sequence

diff -r 70f3967ca6eb -r 81563e019fa4 lemon/random.h
--- a/lemon/random.h	Thu Mar 20 12:12:24 2008 +0000
+++ b/lemon/random.h	Thu Mar 20 17:15:35 2008 +0100
@@ -575,6 +575,26 @@
       return *this;
     }
 
+    /// \brief Seeding random sequence
+    ///
+    /// Seeding the random sequence. The current number type will be
+    /// converted to the architecture word type.
+    template <typename Number>
+    void seed(Number seed) { 
+      _random_bits::Initializer<Number, Word>::init(core, seed);
+    }
+
+    /// \brief Seeding random sequence
+    ///
+    /// Seeding the random sequence. The given range should contain
+    /// any number type and the numbers will be converted to the
+    /// architecture word type.
+    template <typename Iterator>
+    void seed(Iterator begin, Iterator end) { 
+      typedef typename std::iterator_traits<Iterator>::value_type Number;
+      _random_bits::Initializer<Number, Word>::init(core, begin, end);
+    }
+
     /// \brief Returns a random real number from the range [0, 1)
     ///
     /// It returns a random real number from the range [0, 1). The
diff -r 70f3967ca6eb -r 81563e019fa4 test/random_test.cc
--- a/test/random_test.cc	Thu Mar 20 12:12:24 2008 +0000
+++ b/test/random_test.cc	Thu Mar 20 17:15:35 2008 +0100
@@ -24,6 +24,8 @@
 ///\todo To be extended
 ///
 
+int seed_array[] = {1, 2};
+
 int main()
 {
   double a=lemon::rnd();
@@ -34,4 +36,10 @@
   //Does gamma work with integer k?
   a=lemon::rnd.gamma(4.0,0);
   a=lemon::rnd.poisson(.5);
+
+  lemon::rnd.seed(100);
+  lemon::rnd.seed(seed_array, seed_array + 
+		  (sizeof(seed_array) / sizeof(seed_array[0])));
+
+  return 0;
 }