Changeset 1468 for applications/pmsm/simulator_zdenek/ekf_example/mpf_double.cpp

Timestamp:

12/16/12 23:05:25 (11 years ago)

Author:

smidl

Message:

uprava generatoru

Files:

• applications/pmsm/simulator_zdenek/ekf_example/mpf_double.cpp (modified) (5 diffs)

applications/pmsm/simulator_zdenek/ekf_example/mpf_double.cpp

@@ -4 +4 @@
 #include "fastexp.h"
 #include "mpf_double.h"
-#include "RArandn.h"
 
 static floatx om[N];
@@ -21 +20 @@
 static floatx rinv;
 
-
+#if 1
+#include "RArandn.h"
 floatx qrandu() {
         if ( randu_i==randu_last )
@@ -39 +39 @@
         return *randn_i;
 }
+
+void rng_init()
+{
+randu_i=&RArandu[0];
+randn_i=&RArandn[0];
+}
+
+#else
+
+#include "rnor.h"
+floatx qrandu() {
+        return UNI;
+}
+
+floatx qrandn() {
+        return PolarRNG();
+}
+
+void rng_init()
+{
+zigset(86947731 );
+}
+
+
+#endif
+
 
 void resample() {
@@ -213 +239 @@
 
 void mpf_init(floatx qom0, floatx qth0, floatx r0) {
+        
+        rng_init();
+        
         int i;
         r=r0;
@@ -252 +281 @@
         for (i=0;i<N;i++) th1[i]=th[i];
 }
+void mpf_om(floatx om1[N]){
+        int i;
+        for (i=0;i<N;i++) om1[i]=om[i];
+}
 
 floatx kalman_om( const floatx isa, const floatx isb , const floatx usa, const floatx usb, const floatx th ) {