Revision 9e5b9bb8a0564f8bf629df4ca5e96d3078f50363 authored by jmc on 14 July 2008, 00:00:00 UTC, committed by Gabor Csardi on 14 July 2008, 00:00:00 UTC
1 parent a3b36b1
Raw File
triDiag.R
triDiag <- function(diagonal, upper, lower,
                    nrow = length(diagonal), ncol = nrow) {
    value <- diag(diagonal, nrow, ncol)
    R <- row(value)
    C <- col(value)
    value[C == R + 1] <- upper
    value[C == R - 1] <- lower
    value
}
back to top