Hi,
suppose you have two versions of the same algorithm: one in pure R, the
other one in C/C++ called via .Call().
Assuming there is no bug in the implementations (i.e. they both do the
same thing), is there any well known reason why the C/C++ implementation
could return numerical results non
Thank you Duncan for your reply.
Currently I am using 'double' for the computations.
What type should I use for extended real in my intermediate computations?
The result will still be 'double' anyway right?
On 10/09/2010 13:00, Duncan Murdoch wrote:
On 10/09/2010 6:46 AM, Renaud Gaujoux
On Fri, Sep 10, 2010 at 11:46 AM, Renaud Gaujoux
ren...@mancala.cbio.uct.ac.za wrote:
Hi,
suppose you have two versions of the same algorithm: one in pure R, the
other one in C/C++ called via .Call().
Assuming there is no bug in the implementations (i.e. they both do the same
thing), is
Ok.
I quickly tried using LDOUBLE wherever I could, but it did not changed
the results. I might try harder...
I agree with you Barry, and I will re-double re-check my code.
Thank you both for your help.
Bests,
Renaud
On 10/09/2010 13:24, Barry Rowlingson wrote:
On Fri, Sep 10, 2010 at 11:46
of the same problems.
Paul
-Original Message-
From: r-devel-boun...@r-project.org [mailto:r-devel-boun...@r-
project.org] On Behalf Of Renaud Gaujoux
Sent: September 10, 2010 6:47 AM
To: r-devel@r-project.org
Subject: [Rd] Non identical numerical results from R code vs C/C++
code?
Hi,
suppose
[mailto:r-devel-boun...@r-
project.org] On Behalf Of Renaud Gaujoux
Sent: September 10, 2010 6:47 AM
To: r-devel@r-project.org
Subject: [Rd] Non identical numerical results from R code vs C/C++
code?
Hi,
suppose you have two versions of the same algorithm: one in pure R, the
other one in C
On Fri, 10 Sep 2010, Duncan Murdoch wrote:
On 10/09/2010 7:07 AM, Renaud Gaujoux wrote:
Thank you Duncan for your reply.
Currently I am using 'double' for the computations.
What type should I use for extended real in my intermediate computations?
I think it depends on compiler details. On