zhaohf_2000@yahoo.com science forum beginner
Joined: 06 Jul 2006
Posts: 5

Posted: Thu Jul 06, 2006 6:02 pm Post subject:
how to efficiently solve a relatively large scale linear equation system?



Hello, guys,
I am doing a project relating to solve relatively large scalce linear
equation system, here is my equations:
(A + T' *B * T) * X = (c + T' *d)
Unknowns: X (N*1 vector, N is between 100 and 1000)
Knowns:
AN*N symmetric narrow band matrix with band width= 5 (actually it
comes from assembling stiffness matrix of scale 6*6)
BN*N symmetric narrow band matrix with band width= 2
TN*N lower triangular matrix with first 3 elements in diagnol are
zeros. (this means this
is a singular matrix)
T'transpose of T
cN*1 vector
dN*1 vector
I need to solve this equation system in C++ efficiently and the speed
for solving the equation is CRITICAL for me.
Anybody has any idea?
Thanks in advance!
Helen 
