[Yuima-commits] r727 - in pkg/yuima: . src

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

Author: iacus
Date: 2020-03-05 14:36:59 +0100 (Thu, 05 Mar 2020)
New Revision: 727

Modified:
   pkg/yuima/DESCRIPTION
   pkg/yuima/NEWS
   pkg/yuima/src/euler.c
Log:
fixed potential PROTECT problems in euler.c

Modified: pkg/yuima/DESCRIPTION
===================================================================
--- pkg/yuima/DESCRIPTION	2020-03-05 11:24:54 UTC (rev 726)
+++ pkg/yuima/DESCRIPTION	2020-03-05 13:36:59 UTC (rev 727)
@@ -1,7 +1,7 @@
 Package: yuima
 Type: Package
 Title: The YUIMA Project Package for SDEs
-Version: 1.9.9
+Version: 1.10.0
 Depends: R(>= 2.10.0), methods, zoo, stats4, utils, expm, cubature,
         mvtnorm
 Imports: Rcpp (>= 0.12.1), boot (>= 1.3-2), glassoFast

Modified: pkg/yuima/NEWS
===================================================================
--- pkg/yuima/NEWS	2020-03-05 11:24:54 UTC (rev 726)
+++ pkg/yuima/NEWS	2020-03-05 13:36:59 UTC (rev 727)
@@ -64,4 +64,5 @@
 2019/7/19: added cce.factor.R, cce.factor.Rd
            modified llag.test.Rd, bns.test.Rd, hyavar.R, mpv.Rd, mllag.Rd, llag.Rd, cce.Rd, spectralcov.R, JBtest.Rd, snr.Rd
 2020/2/6: fixed length(class(matrix))>2 issue. Code for simPoi fastened.
+2020/3/5: should have fixed qmleLevy if() condition and rhck PROTECT problems
 

Modified: pkg/yuima/src/euler.c
===================================================================
--- pkg/yuima/src/euler.c	2020-03-05 11:24:54 UTC (rev 726)
+++ pkg/yuima/src/euler.c	2020-03-05 13:36:59 UTC (rev 727)
@@ -70,8 +70,14 @@
         /* PROTECT(b0 = eval(drift, rho));
         PROTECT(sigma0 = eval(diffusion, rho));
         AS_NUMERIC is added by YK (Dec 4, 2018) */
-        PROTECT(b0 = AS_NUMERIC(eval(drift, rho)));
+        /* PROTECT(b0 = AS_NUMERIC(eval(drift, rho)));
         PROTECT(sigma0 = AS_NUMERIC(eval(diffusion, rho)));
+        */
+        PROTECT(b0 = allocVector(REALSXP, d));
+        b0 = AS_NUMERIC(eval(drift, rho));
+        PROTECT(sigma0 = allocVector(REALSXP, d*r));
+        sigma0 = AS_NUMERIC(eval(diffusion, rho));
+        
         b = REAL(b0);
         sigma = REAL(sigma0);