From tlumley@u.washington.edu  Fri Apr  6 16:49:34 2001
From: tlumley@u.washington.edu (Thomas Lumley)
Date: Fri, 6 Apr 2001 08:49:34 -0700 (PDT)
Subject: [Rd] bug in contour() ? (PR#897)
In-Reply-To: <200104061421.QAA26081@pubhealth.ku.dk>
Message-ID: <Pine.A41.4.33.0104060843400.89396-100000@homer01.u.washington.edu>

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


[Rd] bug in contour() ? (PR#897)
Thomas Lumley