Changeset 579 for library/bdm/itpp_ext.cpp

Timestamp:

08/22/09 14:50:03 (16 years ago)

Author:

smidl

Message:

Implementation of mixtools randun function

Files:

• library/bdm/itpp_ext.cpp (modified) (1 diff)

library/bdm/itpp_ext.cpp

@@ -376 +376 @@
         }
 }
-}
+
+class RandunStorage{
+        const double A;
+        const double M;
+        static double seed;
+        static int counter;
+        public:
+                RandunStorage(): A(16807), M(2147483647) {};
+                void set_seed(double seed0){seed=seed0;}
+                double get(){seed=mod(A*seed,M); counter++; return seed/M;}
+};
+static RandunStorage randun_global_storage;
+double RandunStorage::seed=1111111;
+int RandunStorage::counter=0;
+double randun(){return randun_global_storage.get();};
+vec randun(int n){vec res(n); for(int i=0;i<n;i++){res(i)=randun();}; return res;};
+mat randun(int n, int m){mat res(n,m); for(int i=0;i<n*m;i++){res(i)=randun();}; return res;};
+
+}