Prof. Brian Ripley wrote: >Why not use exp(y*log(x)) if it is adequate for your purposes? It is >faster under Windows.
With the CRAN version of R 1.8.0 under Win Xp: > system.time(exp(3.45*log(runif(1000000)))) [1] 1.30 0.06 1.49 NA NA > system.time(runif(1000000)^3.45) [1] 7.14 0.03 7.20 NA NA With R compiled using MingW 2.0.1 on the same machine: > system.time(exp(3.45*log(runif(1000000)))) [1] 1.31 0.02 1.35 NA NA > system.time(runif(1000000)^3.45) [1] 1.04 0.00 1.04 NA NA OK, this is fine for me. I define: "%^%" <- function(a, b) return(exp(b*log(a))) which I use as a substitute for ^ to make sure performance does not drop too much under Windows with the current CRAN version of R in my application. Thanks, Philippe Grosjean ______________________________________________ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-devel
