Re: [Rd] Philosophy behind converting Fortran to C for use in R

2017-06-07 Thread Martyn Byng
Hi, Just a quick comment on (1). The C-Fortran interface has been standardized since Fortran 2003. However, it does require the Fortran interface that is being called from C to have been written with C operability in mind as specific C interoperable types etc. must be used. Trying to call a

Re: [Rd] ALLOCATE in a FORTRAN subroutine

2016-03-07 Thread Martyn Byng
blems at compile time than the C ones). Martyn -Original Message- From: MAURICE Jean - externe [mailto:jean-externe.maur...@edf.fr] Sent: 07 March 2016 15:46 To: Martyn Byng Cc: r-devel@r-project.org Subject: RE: ALLOCATE in a FORTRAN subroutine Hi Martyn, Many thanks for your answer.

Re: [Rd] ALLOCATE in a FORTRAN subroutine

2016-03-04 Thread Martyn Byng
Hi, Until you get a more definitive answer, I will make an attempt to give some advice. When using an assumed sized array (i.e. REAL*8 array1(*)) you still need to allocate the memory prior to calling the Fortran subroutine, so you would still need to know its maximum length. Arrays created i

Re: [Rd] Error in simulation. NAN

2013-08-27 Thread Martyn Byng
Hi, You are going to have to track down exactly where the NaN is being produced. These usually occur due to performing invalid numerical operations, like trying to take the square-root of a negative number, trying to calculate 0/0 etc. Try adding a series of if / then blocks prior to performing a

Re: [Rd] a^b when a is large and b < 1 (64bit R on windows 7)

2011-11-22 Thread Martyn Byng
have specified I only checked on Windows. So: On Windows 7 > 64bit, the R-32bit works fine, the R-64bit gives the behaviour Martyn > reported. > > Cheers > Joris > > On Mon, Nov 21, 2011 at 3:58 PM, Brian G. Peterson wrote: >> On Mon, 2011-11-21 at 14:41 +, Martyn Byng w

[Rd] a^b when a is large and b < 1 (64bit R on windows 7)

2011-11-21 Thread Martyn Byng
Hi, I'm getting some strange behaviour when trying to use the power operator (a^b) when a is large and b is less than one: big <- .Machine$double.xmax big big^0.5 sqrt(big) > big <- 1.797693134862315708384e+308 > big^0.5 [1] Inf > sqrt(big) [1] 1.340781e+154 I'm guessing that this behaviour is

Re: [Rd] Randomness not due to seed

2011-07-20 Thread Martyn Byng
Hi, Even using the same math libraries you can get different results, depending on what sorts of instructions those libraries use, see the following (none R related) blog article: http://blog.nag.com/2011/02/wandering-precision.html. Martyn -Original Message- From: r-devel-boun...@r-proj

Re: [Rd] Calling Rscript from Makevars

2011-05-23 Thread Martyn Byng
-project.org] Sent: 20 May 2011 20:17 To: Martyn Byng Cc: r-devel@r-project.org Subject: Re: [Rd] Calling Rscript from Makevars On May 20, 2011, at 12:04 PM, Martyn Byng wrote: > Hi, > > I am trying to package some code to use with R and wanted to call > Rscript from within the Makevar

[Rd] Calling Rscript from Makevars

2011-05-20 Thread Martyn Byng
Hi, I am trying to package some code to use with R and wanted to call Rscript from within the Makevars file (I am trying to automate the setting of the location of a third party library depending on what is available / the system the package is being installed on). If I just have a simple Makevar