Re: [Rd] In function isum in summary.c, k should be R_xlen_t

2020-12-31 Thread Jan Gorecki
Don't know this piece well but I am guessing that you haven't found an example because the iterator is going up to the length of a vector anymore but only to the number of batches, which is unlikely to be more than 2^31. On Tue, Dec 22, 2020 at 12:30 PM Suharto Anggono Suharto Anggono via R-devel

[Rd] In function isum in summary.c, k should be R_xlen_t

2020-12-31 Thread Suharto Anggono Suharto Anggono via R-devel
In summary.c, in function 'isum', the loop is 'ITERATE_BY_REGION' that contains 'for' loop for (int k = 0; k < nbatch; k++) It is since SVN revision 73445, in released R since version 3.5.0. Previously, the loop is for (R_xlen_t i = 0; i < n; i++) Inside 'ITERATE_BY_REGION', the type of the ind

Re: [Rd] From .Fortran to .Call?

2020-12-31 Thread Bill Dunlap
As the proverbial naive R (ab)user I’m left wondering: o if I updated my quantreg_init.c file in accordance with Bill’s suggestion could I then simply change my .Fortran calls to .Call? No. .Call(C_func, arg1, arg2) expects C_func's arguments to all be SEXP* (pointers to

Re: [Rd] From .Fortran to .Call?

2020-12-31 Thread Bill Dunlap
As for speed, I suspect much of the speedup when converting from .Fortran to .Call will be due to checking the validity of the arguments (lengths, types, ranges, etc.) in C code instead of in R code. There can also be less copying of the arguments and the returned objects will tend to be smaller.

Re: [Rd] From .Fortran to .Call?

2020-12-31 Thread Erin Hodgess
Will there be any interest in using Coarrays Fortran, please? Thanks, Erin Erin Hodgess, PhD mailto: erinm.hodg...@gmail.com On Wed, Dec 23, 2020 at 10:29 AM Martin Maechler wrote: > > Balasubramanian Narasimhan > > on Wed, 23 Dec 2020 08:34:40 -0800 writes: > > > I think it s

Re: [Rd] From .Fortran to .Call?

2020-12-31 Thread Erin Hodgess
Is .Fortran going to be deprecated, please? I have gotten amazing speed up with geostatistics processes using HPC type tools. Thanks On Sat, Dec 26, 2020 at 9:48 AM Koenker, Roger W wrote: > I’ve recoded a version of one of my quantile regression fitting functions > to use .C64 from dotCall64

Re: [Rd] From .Fortran to .Call?

2020-12-31 Thread Avraham Adler
I’ve tried recoding some of Delaporte to use the .Fortran interface and I don’t know what I’m doing wrong but it either doesn’t work or crashes my R instance completely. Avi On Sat, Dec 26, 2020 at 11:48 AM Koenker, Roger W wrote: > I’ve recoded a version of one of my quantile regression fittin

Re: [Rd] From .Fortran to .Call?

2020-12-31 Thread Koenker, Roger W
Erin, I think rumors of the deprecation of .Fortran are greatly exaggerated, but I’d welcome some confirmation of this from someone in R core. There is quite a lot of .Fortran usage in packages, and perhaps even in base R... > On Dec 26, 2020, at 7:57 PM, Erin Hodgess wrote: > > Is .Fortran