search for: r_pow

In R, 0 ^ NaN yields Inf. I would have expected NaN or perhaps 0. Is this behaviour intended? >sessionInfo() R version 2.8.0 (2008-10-20) i386-pc-mingw32 locale: LC_COLLATE=English_United States.1252;LC_CTYPE=English_United States.1252;LC_MONETARY=English_United States.1252;LC_NUMERIC=C;LC_TIME=English_United States.1252 attached base packages: [1] stats graphics grDevices utils

at lines 651 & 653 (integer_binary function): if (code == DIVOP || code == POWOP) ans = allocVector(REALSXP, n); else ans = allocVector(INTSXP, n); There are calls to warningcall() later in the function, which can trigger garbbage collection. Looks like the typical scenario where it seemed pretty safe to not PROTECT in the original version of the function but

..., but this introduces a whole lot of problems > that were not there before. In particular R crashes when using the > graphics device, which it did not without this fix. I found that the problem goes away when we use isnanl() instead of isnan(). This is very similar to the problem we have for R_pow() where we need to use powl() instead of pow() for recent versions mingw-w64 (this still needs to be patched in r-devel). So the full solution is: #if (defined(_WIN32) || defined(_WIN64)) && defined(__GNUC__) # define R_sqrt(x) (isnanl(x) ? x : sqrt(x)) #else # define R_sqrt sqrt #endif...

On 14/09/2015 9:36 AM, luke-tierney at uiowa.edu wrote: > I believe the issue is that on Windows the sqrt function when called > with a NaN does not return the same NaN, as it does on other platforms. > We have > > #if (defined(_WIN32) || defined(_WIN64)) && defined(__GNUC__) && \ > __GNUC__ <= 4 > # define R_sqrt(x) (ISNAN(x) ? x : sqrt(x)) > #else

I've continued to work on speeding up R, and now have a collection of fourteen patches, some of which speed up particular functions, and some of which reduce general interpretive overhead. The total speed improvement from these patches is substantial. It varies a lot from one R program to the next, of course, and probably from one machine to the next, but speedups of 25% can be expected in

...e/v4.2.1/lib -lreadline -ltermcap -lm -lUfor -lfor -lFutil -lots -lm_c32 -lmld -lexc ld: Warning: Unresolved: R_registerRoutines R_NaReal Rf_error Rf_warning R_alloc R_chk_calloc R_chk_free rwarn_ interv_ d1mach_ dsvdc_ dqrsl_ dqrdc_ ddot_ idamax_ dnrm2_ dblepr_ intpr_ R_finite Rprintf REprintf R_pow Rf_fsign dpbsl_ dpbfa_ mkdir -p -- ../../../../library/modreg/libs gmake[5]: Leaving directory `/soft/R/R-1.5.1/src/library/modreg/src' gmake[4]: Leaving directory `/soft/R/R-1.5.1/src/library/modreg/src' gmake[3]: Leaving directory `/soft/R/R-1.5.1/src/library/modreg' gmake[3]: Enterin...