solve.pentadiag function

Pentadiagonal Linear Solver.

A function to solve pentadiagonal system of linear equations. data

## S3 method for class 'pentadiag'
solve(a, b, ...)

Arguments

• a: a numeric square matrix with pentadiagonal rows. The function does NOT check for pentadiagonal matrix.
• b: a numeric vector of the same length as nrows(a). This argument cannot be a matrix.
• ...: any additional arguments

Details

This function is written mainly for use in this package. It may not be the most efficient code.

Returns

A vector containing the solution.

Author(s)

Arun Kumar Kuchibhotla, arunku@wharton.upenn.edu

Examples

A <- matrix(c(2,1,1,0,0,
                          1,2,1,1,0,
                          1,1,2,1,1,
                          0,1,1,2,1,
                          0,0,1,1,2),nrow = 5)
b <- rnorm(5)
tmp <- solve.pentadiag(A, b)