[Rd] Look, Watson! La.svd & ATLAS

[Douglas Bates]

"Liaw, Andy" <andy_liaw@merck.com> writes:

> Dear R-devel,
> 
> I had attempted to compile r-devel (dated Oct. 31, 2001) on WinNT with link
> to ATLAS, with mostly success.  However, when I tried the following, I got a
> visit from Dr. Watson:
> 
> 
> R : Copyright 2001, The R Development Core Team
> Version 1.4.0 Under development (unstable) (2001-10-31)
> 
> R is free software and comes with ABSOLUTELY NO WARRANTY.
> You are welcome to redistribute it under certain conditions.
> Type `license()' or `licence()' for distribution details.
> 
> R is a collaborative project with many contributors.
> Type `contributors()' for more information.
> 
> Type `demo()' for some demos, `help()' for on-line help, or
> `help.start()' for a HTML browser interface to help.
> Type `q()' to quit R.
> 
> > set.seed(1)
> > x <- matrix(runif(1e4), 1e2, 1e2)
> > invisible(La.svd(x))                         # This ran fine.
> > x <- matrix(runif(1e5), 1e3, 1e2)
> > invisible(La.svd(x))
> 
> (Dr. Watson says there's an access violation.)

I was able to run those examples without ill effects on an R-1.4.0
(development) (2001-10-30) with Atlas 3.2.1 on a Linux system (Debian
testing)

> set.seed(1)
> x <- matrix(runif(1e4), 1e2, 1e2)
> invisible(La.svd(x))
> x <- matrix(runif(1e5), 1e3, 1e2)
> invisible(La.svd(x))
> x <- matrix(runif(1e4), 1e2, 1e2)
> system.time(invisible(La.svd(x)))
[1] 0.04 0.00 0.04 0.00 0.00
> x <- matrix(runif(1e5), 1e3, 1e2)
> system.time(invisible(La.svd(x)))
[1] 1.15 0.01 2.39 0.00 0.00
> system.time(La.svd(x, nu = 0, nv = 0))
[1] 0.33 0.00 0.63 0.00 0.00
-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
r-devel mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html
Send "info", "help", or "[un]subscribe"
(in the "body", not the subject !)  To: r-devel-request@stat.math.ethz.ch
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._