PowerfulMaxEigenpair (version 0.1.1)

thomas.tri.sol: Solve the linear equation (-Q-zI)w=v.

Description

Construct the solution of linear equation (-Q-zI)w=v.

Usage

thomas.tri.sol(Q, v, z)

Arguments

Q

The given tridiagonal matrix.

v

The column vector on the right hand of equation.

z

The Rayleigh shift.

Value

A solution sequence \(w\) to the equation (-Q-zI)w=v.

Examples

# NOT RUN {
nn = 8
a = c(1:(nn - 1))^2
b = c(1:(nn - 1))^2
C = c(b[1], a[1:(nn - 2)] + b[2:(nn - 1)], a[nn - 1] + nn^2)
Q = tridiag(b, a, -C)
zstart = 6
thomas.tri.sol(Q, z=zstart, v=rep(1,dim(Q)[1]))
# }