See do_summary() in summary.c, rsum() for doubles. R uses long double
type as accumulator on systems where available.
Best,
Tomas
On 2/14/19 2:08 PM, Rampal Etienne wrote:
Hello,
I am trying to write FORTRAN code to do the same as some R code I
have. I get (small) differences when using the sum function in R. I
know there are numerical routines to improve precision, but I have not
been able to figure out what algorithm R is using. Does anyone know
this? Or where can I find the code for the sum function?
Regards,
Rampal Etienne
______________________________________________
R-devel@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-devel
______________________________________________
R-devel@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-devel
