The following example just "hangs", i.e. the final contour() does not terminate. x <- pi*(0.5+((-13):12))/36 w <- 0.377*cos(2*x)*sqrt(abs(cos(2*x))) v <- c(rep(-1.47,13),rep(1.722,13))*abs(sin(x))^1.5 cr <- outer(w,w)+outer(v,v) s <- sqrt(diag(cr)) cr <- cr/outer(s,s) contour(cr[1:26,26:1],levels=seq(-1,1,0.1),axes=T) --please do not edit the information below-- Version: platform = i686-pc-linux-gnu arch = i686 os = linux-gnu system = i686, linux-gnu status = major = 1 minor = 2.2 year = 2001 month = 02 day = 26 language = R Search Path: .GlobalEnv, package:ctest, package:SfS, Autoloads, package:base -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- 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 _._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._
kuensch@stat.math.ethz.ch wrote:> > The following example just "hangs", i.e. > the final contour() does not terminate. > > x <- pi*(0.5+((-13):12))/36 > w <- 0.377*cos(2*x)*sqrt(abs(cos(2*x))) > v <- c(rep(-1.47,13),rep(1.722,13))*abs(sin(x))^1.5 > cr <- outer(w,w)+outer(v,v) > s <- sqrt(diag(cr)) > cr <- cr/outer(s,s) > contour(cr[1:26,26:1],levels=seq(-1,1,0.1),axes=T) > > --please do not edit the information below-- > > Version: > platform = i686-pc-linux-gnu > arch = i686 > os = linux-gnu > system = i686, linux-gnu > status > major = 1 > minor = 2.2 > year = 2001 > month = 02 > day = 26 > language = R... and R crashes on WinNT 4.0. Uwe Ligges -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- 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 _._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._
On Fri, 6 Apr 2001 kuensch@stat.math.ethz.ch wrote:> The following example just "hangs", i.e. > the final contour() does not terminate. > > x <- pi*(0.5+((-13):12))/36 > w <- 0.377*cos(2*x)*sqrt(abs(cos(2*x))) > v <- c(rep(-1.47,13),rep(1.722,13))*abs(sin(x))^1.5 > cr <- outer(w,w)+outer(v,v) > s <- sqrt(diag(cr)) > cr <- cr/outer(s,s) > contour(cr[1:26,26:1],levels=seq(-1,1,0.1),axes=T)It's following a list that ends up being circular (loop of about 10 entries) around plot3d.c:697. I don't understand the code that sets up the segment database, but presumably there's a bug there. The loop happens in the same place under gctorture, so I don't think it's memory management. -thomas> --please do not edit the information below-- > > Version: > platform = i686-pc-linux-gnu > arch = i686 > os = linux-gnu > system = i686, linux-gnu > status > major = 1 > minor = 2.2 > year = 2001 > month = 02 > day = 26 > language = R > > Search Path: > .GlobalEnv, package:ctest, package:SfS, Autoloads, package:base > > -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- > 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 > _._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._ >Thomas Lumley Asst. Professor, Biostatistics tlumley@u.washington.edu University of Washington, Seattle ^^^^^^^^^^^^^^^^^^^^^^^^ NOTE NEW EMAIL ADDRESS -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- 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 _._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._