Re: [Rd] [R] aligned memory allocation in C

2008-08-14 Thread christophe dutang
Thanks for your remarks. According to http://www.math.sci.hiroshima-u.ac.jp/~m-mat/MT/SFMT/howto-compile.html, I need 16 bit aligned memory when using fill_array64. So I suppose I need 8 bit aligned memory. I will test what you advise me and will come back to R-devel list after. Thanks again

Re: [Rd] extending the derivs table/fools rushing in

2008-08-14 Thread Gabor Grothendieck
While you are at it could you add { to the table so that this works: # this is ok f - function(x) x*x D(body(f), x) x + x # but not g which is same as f # except it has { ... } surrounding its body g - function(x) { x*x } D(body(g), x) Error in D(body(g), x) : Function '`{`' is not in the

Re: [Rd] extending the derivs table/fools rushing in

2008-08-14 Thread Ben Bolker
http://www.zoology.ufl.edu/bolker/deriv_patch2.txt has this change added as well. However, I'm not as confident that this is the right thing to do? Should curly brackets even be appearing in mathematical expressions? Ben Gabor Grothendieck wrote: While you are at it could you add { to

Re: [Rd] extending the derivs table/fools rushing in

2008-08-14 Thread Prof Brian Ripley
The derivative of plogis is surely dlogis. (And yes, there is a good reason why we have such a function: take a look at its C code.) That means we would need an entry for dlogis too, I guess. I am not convinced that there is a real need for these (and where does this stop?) What would be

Re: [Rd] extending the derivs table/fools rushing in

2008-08-14 Thread Ben Bolker
Prof Brian Ripley wrote: The derivative of plogis is surely dlogis. (And yes, there is a good reason why we have such a function: take a look at its C code.) Doh. That means we would need an entry for dlogis too, I guess. I am not convinced that there is a real need for these (and where

Re: [Rd] extending the derivs table/fools rushing in

2008-08-14 Thread Gabor Grothendieck
You could see if Ryacas package can do what you want: http://ryacas.googlecode.com 2008/8/14 Ben Bolker [EMAIL PROTECTED]: Prof Brian Ripley wrote: The derivative of plogis is surely dlogis. (And yes, there is a good reason why we have such a function: take a look at its C code.) Doh.

[Rd] Typo in proc.time.Rd

2008-08-14 Thread Stephen Weigand
There's a small typo in proc.time.Rd: \sQuote{system lime} should be \sQuote{system time} Thanks, Stephen -- Rochester, Minn. USA __ R-devel@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-devel

Re: [Rd] [R] RNG Cycle and Duplication (PR#12538)

2008-08-14 Thread murdoch
Shengqiao Li wrote: Hello all, I am generating large samples of random numbers. The RNG help page says: All the supplied uniform generators return 32-bit integer values that are converted to doubles, so they take at most 2^32 distinct values and long runs will return duplicated values.

Re: [Rd] [R] RNG Cycle and Duplication (PR#12540)

2008-08-14 Thread shli
This message is in MIME format. The first part should be readable text, while the remaining parts are likely unreadable without MIME-aware tools. ---559023410-851401618-1218751024=:15885 Content-Type: TEXT/PLAIN; charset=ISO-8859-1; format=flowed Content-Transfer-Encoding: QUOTED-PRINTABLE

Re: [Rd] [R] RNG Cycle and Duplication (PR#12540)

2008-08-14 Thread Prof Brian Ripley
Remember Wichmann-Hill is a composite generator. Its composition does take more than 2^32 distinct values. You still haven't identifed a problem here. The note is to warn that runif() does repeat within a cycle, because people wrote code assuming otherwise. It would be a poor use of