From bates@stat.wisc.edu  Fri Nov  2 14:27:01 2001
From: bates@stat.wisc.edu (Douglas Bates)
Date: 02 Nov 2001 08:27:01 -0600
Subject: [Rd] Look, Watson!  La.svd & ATLAS
In-Reply-To: <51F9C42DA15CD311BD220008C707D81903DC8350@usrymx10.merck.com>
References: <51F9C42DA15CD311BD220008C707D81903DC8350@usrymx10.merck.com>
Message-ID: <6rwv198czu.fsf@franz.stat.wisc.edu>

"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
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._


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