Re: [R] Trouble returning 2D array into R from Fortran

2012-10-23 Thread paulfjbrowne
Thank you for your help. In any case, I eventually realized that I had reversed the Fortran/R column major order in the Fortran, which is why elements were returned incorrectly. As you point out, I did not run into memory access issues because of the self-similar array sizes. -- View this messa

Re: [R] Trouble returning 2D array into R from Fortran

2012-10-23 Thread Berend Hasselman
On 23-10-2012, at 12:33, paulfjbrowne wrote: > Yes, I had stripped out the Fortran from the geta subroutine downwards > because is wasn't very important in returning a test output matrix. > But then any results are not reproducible. So there is no purpose in trying to reproduce your results. S

Re: [R] Trouble returning 2D array into R from Fortran

2012-10-23 Thread paulfjbrowne
Yes, I had stripped out the Fortran from the geta subroutine downwards because is wasn't very important in returning a test output matrix. I think that I should make explicit what problem I believe I am seeing in what I have coded, and what form of the output matrix I would prefer. When I run th

Re: [R] Trouble returning 2D array into R from Fortran

2012-10-22 Thread Berend Hasselman
Your attachments did not make it through to the official R-help mailing list. Nabble R is not R-help. I downloaded the files from the links to Nabble. See further down for my comments. On 22-10-2012, at 14:24, paulfjbrowne wrote: > I haven't determined the root cause of the numerical accuracy bu

Re: [R] Trouble returning 2D array into R from Fortran

2012-10-22 Thread paulfjbrowne
I haven't determined the root cause of the numerical accuracy bug yet (I haven't determined a way to compare the values used when the code is run from R as opposed to from a Fortran test), but I did have another issue with retrieving the matrix from the Fortran code that I hoped you might be able t

Re: [R] Trouble returning 2D array into R from Fortran

2012-10-20 Thread Berend Hasselman
Look at my comments in between your post. On 20-10-2012, at 19:18, paulfjbrowne wrote: > I will look into using inline, but since the Fortran code is several thousand > lines long & is comprised of multiple subroutines, compiling it into a > shared object & dynamically loading it into R is prob

Re: [R] Trouble returning 2D array into R from Fortran

2012-10-20 Thread paulfjbrowne
I will look into using inline, but since the Fortran code is several thousand lines long & is comprised of multiple subroutines, compiling it into a shared object & dynamically loading it into R is probably the easier solution. I have also noticed a strange numerical problem when calling the routi

Re: [R] Trouble returning 2D array into R from Fortran

2012-10-20 Thread Berend Hasselman
See inline. On 20-10-2012, at 17:18, paulfjbrowne wrote: > Hello, > > I have been trying to use a collection of Fortran subroutines to return a 2D > array of calculated values to my R code, calling a Fortran wrapper > subroutine from R. I've done this successfully before with C & C++ code. >

[R] Trouble returning 2D array into R from Fortran

2012-10-20 Thread paulfjbrowne
Hello, I have been trying to use a collection of Fortran subroutines to return a 2D array of calculated values to my R code, calling a Fortran wrapper subroutine from R. I've done this successfully before with C & C++ code. The Fortran wrapper subroutine which is to be called by R takes a set of