[Pomp-commits] r649 - in pkg: inst src

[noreply at r-forge.r-project.org]

Author: kingaa
Date: 2012-04-10 15:51:56 +0200 (Tue, 10 Apr 2012)
New Revision: 649

Modified:
   pkg/inst/NEWS
   pkg/src/trajectory.c
Log:
- update NEWS file
- minor tweak to trajectory.c


Modified: pkg/inst/NEWS
===================================================================
--- pkg/inst/NEWS	2012-04-10 12:45:52 UTC (rev 648)
+++ pkg/inst/NEWS	2012-04-10 13:51:56 UTC (rev 649)
@@ -1,4 +1,8 @@
 NEWS
+0.41-3
+     o	In 'trajectory', all vectorfield and map evaluation is now done in C for speed.
+     	For continuous-time dynamical systems, this gives an approximately 4-fold speedup.
+
 0.41-2
      o	Fix bug in 'bbs' example.
 

Modified: pkg/src/trajectory.c
===================================================================
--- pkg/src/trajectory.c	2012-04-10 12:45:52 UTC (rev 648)
+++ pkg/src/trajectory.c	2012-04-10 13:51:56 UTC (rev 649)
@@ -76,11 +76,8 @@
   }
 
   if (nzeros>0) {
-    for (j = 0; j < nreps; j++) {
-      for (i = 0; i < nzeros; i++) {
-	xp[zidx[i]+nvars*j] = 0.0;
-      }
-    }
+    for (j = 0; j < nreps; j++)
+      for (i = 0; i < nzeros; i++) xp[zidx[i]+nvars*j] = 0.0;
   }
 
   for (k = 0; k < ntimes; k++, tp++) {