Re: [Rd] convolve: request for "usual" behaviour + some improvements + some fixes

2007-02-02 Thread Herve Pages
Last but not least: convolve2 can be made 100 times or 1000 times faster than convolve by choosing a power of 2 for the length of the fft-buffer (a length of 2^n is the best case for the fft, the worst case being when the length is a prime number): > x <- 1:13 > y <- 1:1 > system.time(cc <- co

[Rd] convolve: request for "usual" behaviour + some improvements + some fixes

2007-02-02 Thread Herve Pages
Hi again, There are many problems with current 'convolve' function. The author of the man page seems to be aware that 'convolve' does _not_ the usual thing: Note that the usual definition of convolution of two sequences 'x' and 'y' is given by 'convolve(x, rev(y), type = "o")'. and indeed, i

[Rd] Inaccuracy in ?convolve

2007-02-02 Thread hpages
Hi, Man page for 'convolve' says: conj: logical; if 'TRUE', take the complex _conjugate_ before back-transforming (default, and used for usual convolution). The complex conjugate of 'x', of 'y', of both? In fact it seems that it takes the complex conjugate of 'y' only which is OK

Re: [Rd] SEXP i/o, .Call(), and garbage collection.

2007-02-02 Thread Jeffrey Horner
K. B. Udaya wrote: > Apologies for any obtuseness in the following. We have been working > on Version 2.0 of the randomSurvivalForest CRAN package and we're > encountering a perplexing 'memory not mapped' segfault that we believe > is "influenced" by GC. [...] > We are wondering if there is someth

Re: [Rd] Help with OS X (BSD) ps command

2007-02-02 Thread Prof Brian Ripley
Most systems do not have pgrep: it is not POSIX. From Linux: STANDARDS pkill and pgrep were introduced in Sun's Solaris 7. This implementation is fully compatible. The man pages of Darwin and many other systems are online, and linked from developer.r-project.org. You could do

[Rd] Help with OS X (BSD) ps command

2007-02-02 Thread Jeffrey J. Hallman
My fame package has a function that checks to see if a FAME SERVER process is already running. On Linux, I can do this in one of two ways: pid <- Sys.getpid() user <- Sys.info()["user"] cmd <- paste("pgrep -fU", user, "-P", pid, "'FAME SERVER'") fameRunning <- as.logical(length(system(cm