The R Project SVN R

Rev

Rev 26196 | Rev 28529 | Go to most recent revision | Show entire file | Ignore whitespace | Details | Blame | Last modification | View Log | RSS feed

Rev 26196 Rev 27960
Line 123... Line 123...
123
n <-  120; all(abs(rErr( gamma(1:n), cumprod(c(1,1:(n-1))))) < 1000*Meps)
123
n <-  120; all(abs(rErr( gamma(1:n), cumprod(c(1,1:(n-1))))) < 1000*Meps)
124
n <- 10000;all(abs(rErr(lgamma(1:n),cumsum(log(c(1,1:(n-1)))))) < 100*Meps)
124
n <- 10000;all(abs(rErr(lgamma(1:n),cumsum(log(c(1,1:(n-1)))))) < 100*Meps)
125
 
125
 
126
all(is.nan(gamma(0:-4))) # + warn.
126
all(is.nan(gamma(0:-4))) # + warn.
127
 
127
 
-
 
128
## choose() {and lchoose}:
-
 
129
n51 <- c(196793068630200, 229591913401900, 247959266474052)
-
 
130
abs(c(n51, rev(n51))- choose(51, 23:28)) <= 2
-
 
131
all(choose(0:4,2) == c(0,0,1,3,6))
-
 
132
## 3 to 8 units off and two NaN's in 1.8.1
128
 
133
 
129
## fft():
134
## fft():
130
ok <- TRUE
135
ok <- TRUE
131
##test EXTENSIVELY:	for(N in 1:100) {
136
##test EXTENSIVELY:	for(N in 1:100) {
132
    cat(".")
137
    cat(".")