simest (version 0.4)

solve.pentadiag:

Description

A function to solve pentadiagonal system of linear equations.

Usage

# S3 method for 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

Value

A vector containing the solution.

Details

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

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)