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

From: Alireza Mahani <>
Date: Mon, 18 Jul 2011 21:07:04 -0700 (PDT)


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 error accumulation occurs. For example, in differential equations where boundary values are integrated to arrive at interior values, etc. On the other hand, in my personal line of work (Hierarchical Bayesian models for quantitative marketing), we have so much inherent uncertainty and noise at so many levels in the problem (and no significant error accumulation sources) that single vs double precision issue is often inconsequential for us. So I think it really depends on the field as well as the nature of the problem.


View this message in context:
Sent from the R devel mailing list archive at

______________________________________________ mailing list
Received on Tue 19 Jul 2011 - 04:08:58 GMT

This quarter's messages: by month, or sorted: [ by date ] [ by thread ] [ by subject ] [ by author ]

All messages

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.2.0, at Tue 19 Jul 2011 - 08:30:10 GMT.

Mailing list information is available at Please read the posting guide before posting to the list.

list of date sections of archive