Re: [Rd] Manipulating single-precision (float) arrays in .Call functions

2011-07-18 Thread Prof Brian Ripley
On Mon, 18 Jul 2011, Alireza Mahani wrote: Simon, Thank you for elaborating on the limitations of R in handling float types. I think I'm pretty much there with you. As for the insufficiency of single-precision math (and hence limitations of GPU), my personal take so far has been that double-pr

Re: [Rd] [R] Problem compiling in extra/xdr

2011-07-18 Thread Prof Brian Ripley
On Mon, 18 Jul 2011, Allin Cottrell wrote: On Mon, 18 Jul 2011, Prof Brian Ripley wrote: 1) R-help is the wrong list: see the posting guide. I've moved this to R-devel. 2) A glibc system should not be compiling in that directory. glibc 2.14 is rather recent and NEWS does say * The RPC i

Re: [Rd] Manipulating single-precision (float) arrays in .Call functions

2011-07-18 Thread Alireza Mahani
Simon, Thank you for elaborating on the limitations of R in handling float types. I think I'm pretty much there with you. As for the insufficiency of single-precision math (and hence limitations of GPU), my personal take so far has been that double-precision becomes crucial when some sort of erro

Re: [Rd] Manipulating single-precision (float) arrays in .Call functions

2011-07-18 Thread Simon Urbanek
On Jul 18, 2011, at 6:15 PM, Alireza Mahani wrote: > Duncan, > > Thank you for your reply. This is a rather unfortunate limitation, because > for large data sizes there is a significant difference between the > performance of '.C' and '.Call'. I think you may have missed the main point - R do

Re: [Rd] Manipulating single-precision (float) arrays in .Call functions

2011-07-18 Thread Alireza Mahani
Duncan, Thank you for your reply. This is a rather unfortunate limitation, because for large data sizes there is a significant difference between the performance of '.C' and '.Call'. I will have to do some tests to see what sort of penalty I incur for copying from double to float inside my C++ cod

Re: [Rd] [R] Problem compiling in extra/xdr

2011-07-18 Thread Prof Brian Ripley
1) R-help is the wrong list: see the posting guide. I've moved this to R-devel. 2) A glibc system should not be compiling in that directory. glibc 2.14 is rather recent and NEWS does say * The RPC implementation in libc is obsoleted. Old programs keep working but new programs cannot be l

[Rd] filter(x, method="recursive") when there are NA's in x messes up

2011-07-18 Thread William Dunlap
If x contains an NA then filter(x, f, method="recursive") puts NA's and 0.0's in the output in an odd pattern. It puts NA's in the length(f) positions ending at the position of the input NA and it puts a run of length(f) 0.0's starting at the position immediately after the input NA. The values aft

Re: [Rd] Fwd: Understanding R's "Environment" concept

2011-07-18 Thread Greg Snow
Here is an attempt at the general concept without getting technical. How many people in the world answer to the name/title "Dad"? Yet based on context we can usually tell who someone is talking about when they use "Dad". It is the same in programming, I may write a function which includes a v

Re: [Rd] Fwd: Understanding R's "Environment" concept

2011-07-18 Thread Davor Cubranic
On 2011-07-18, at 11:52 AM, Nipesh Bajaj wrote: > Hi all, I am trying to understand the R's "environment" concept > however the underlying help files look quite technical to me. Can > experts here provide me some more intuitive ideas behind this concept > like, why it is there, what exactly it is

[Rd] Fwd: Understanding R's "Environment" concept

2011-07-18 Thread Nipesh Bajaj
*Initially, I posted this topic in R-help however, folks there suggested me to post this in R-devel forum. Here is my problem* Hi all, I am trying to understand the R's "environment" concept however the underlying help files look quite technical to me. Can experts here provide me

Re: [Rd] Manipulating single-precision (float) arrays in .Call functions

2011-07-18 Thread Duncan Murdoch
On 18/07/2011 11:52 AM, Alireza Mahani wrote: I am writing a wrapper function in C++ that calls a GPU kernel. My array type for the GPU kernel is float, so I would like my wrapper function to receive float arrays from R. I understand that I can use 'as.single' in R to copy a double-precision vect

Re: [Rd] Confusing inheritance problem

2011-07-18 Thread Terry Therneau
I've packaged the test library up as a tar file at ftp.mayo.edu directory therneau, file ktest.tar login username: mayoftp password: KPlLiFoz This will disappear in 3 days (Mayo is very fussy about outside access). In response to Uwe's comments 1. "2.13.0"

[Rd] Manipulating single-precision (float) arrays in .Call functions

2011-07-18 Thread Alireza Mahani
I am writing a wrapper function in C++ that calls a GPU kernel. My array type for the GPU kernel is float, so I would like my wrapper function to receive float arrays from R. I understand that I can use 'as.single' in R to copy a double-precision vector from R in single-precision format while using

Re: [Rd] image adds lines

2011-07-18 Thread Uwe Ligges
Yes, thanks. I think I know where this is coming from and will discuss with Duncan Murdoch who was working on this topic recently. Best, Uwe On 13.07.2011 20:06, Robert J. Hijmans wrote: There seems to be a bug in "image" in R 13.1 (on windows 32& 64 bits) and on R-devel that is not pre