| Line 219... |
Line 219... |
| 219 |
|
219 |
|
| 220 |
## axis() when in subnormal range (x is subnormal if 0 < x < .Machine$double.xmin):
|
220 |
## axis() when in subnormal range (x is subnormal if 0 < x < .Machine$double.xmin):
|
| 221 |
(.min.exp.subnormal <- with(.Machine, double.min.exp + double.ulp.digits)) # -1074
|
221 |
(.min.exp.subnormal <- with(.Machine, double.min.exp + double.ulp.digits)) # -1074
|
| 222 |
(ry <- range(y74 <- 2^(.min.exp.subnormal + 0:50))) # all subnormal
|
222 |
(ry <- range(y74 <- 2^(.min.exp.subnormal + 0:50))) # all subnormal
|
| 223 |
plotNchk <- function(y) {
|
223 |
plotNchk <- function(y) {
|
| 224 |
plot(y, log = "y")
|
224 |
plot(y, log = "y", xlab="", ylab="", xaxt="n") # (as small pdf as possible)
|
| 225 |
ry <- range(y)
|
225 |
ry <- range(y)
|
| 226 |
xp <- par("yaxp")
|
226 |
xp <- par("yaxp")
|
| 227 |
## and indeed yaxp do cover the range(y) :
|
227 |
## and indeed yaxp do cover the range(y) :
|
| 228 |
stopifnot(xp[1] <= ry[1], ry[2] <= xp[2], xp[3] == 1)
|
228 |
stopifnot(xp[1] <= ry[1], ry[2] <= xp[2], xp[3] == 1)
|
| 229 |
invisible(xp)
|
229 |
invisible(xp)
|
| 230 |
}
|
230 |
}
|
| - |
|
231 |
if(interactive()) # not regularly, where pdf is stored
|
| 231 |
plotNchk(y74) # gives 3 warnings; 1. from pretty(): "very small range"
|
232 |
plotNchk(y74) # gives 3 warnings; 1. from pretty(): "very small range"
|
| 232 |
plotNchk(y74[1:8]) # 3 warnings *and* no error anymore
|
233 |
plotNchk(y74[1:8]) # 3 warnings *and* no error anymore
|
| 233 |
plotNchk(y74[1:2]) # (ditto)
|
234 |
plotNchk(y74[1:2]) # (ditto)
|