swh:1:snp:a4c99a50dc49f82b591f268001b320f8c3ca0041
Tip revision: ce60f670aac0a708d3a9af5f0cf46a752d46ba8a authored by John M Chambers on 28 October 2020, 07:59:48 UTC
version 1.0-6.1
version 1.0-6.1
Tip revision: ce60f67
triDiag2.R
triDiag2 <- function(diagonal, upper, lower,
nrow = length(diagonal), ncol = nrow) {
value <- diag(diagonal, nrow = nrow, ncol = ncol)
n <- min(nrow, ncol-1)
if(n>0) {
rseq <- 1:n
value[cbind(rseq, rseq+1)] <- upper
}
n <- min(nrow-1, ncol)
if(n > 0) {
rseq <- 1:n
value[cbind(rseq+1, rseq)] <- lower
}
value
}