|
Vectorize Eigen Solver (VES)
Calling VES Library
VES Library consist of One Call: 1. VES_LIB (Solve Eigen Value Problem) CALL VSS_LIB(NEQ,KMAX,MMAX,
DIAG_K, COEFS_K, KPT_K, KIND_K,
DIAG_M, COEFS_M, KPT_M, KIND_M,
IPT, RPT, EIG, XEIG)
NEQ [I] : number of equations
KMAX [I] : number of non-zero coefficients in the Stiffness Matrix
MMAX [I] : number of non-zero coefficients in the Mass Matrix
0 for Lumped Mass Problem.
DIAG_K [I] : diagonal terms in the Stiffness Matrix [size NEQ]
COEFS_K [I] : (off-diagonal) coefficients in row format in the Stiffness Matrix [size KMAX ]
KPT_K [I] : number of (off- diagonal) coefficients in each row in the Stiffness Matrix [size NEQ]
KIND_K [I] : column location of each coefficient in the Stiffness Matrix [size KMAX ]
DIAG_M [I] : diagonal terms [size NEQ] in the Mass Matrix [size NEQ ]
COEFS_M [I] : (off-diagonal) coefficients in row format in the Mass Matrix [size MMAX , OR 0 for Lumped Mass Problem]
KPT_M [I] : number of (off-diagonal) coefficients in each row in the Mass Matrix [size NEQ , OR 0 for Lumped Mass Problem]
KIND_M [I] : column location of each coefficient in the Mass Matrix [size NCOEFF, OR 0 for Lumped Mass Problem]
IPT [I] : Integer Parameter Array [size 20 ]
IPT(1) = Number of Equations
IPT(2) = Number of Eigen Value required
RPT [I] : Real Parameter Array [size 20 ]
EIG [O] : Vector Contains Eigen Values [size (NEQ,2)]
XEIG [O] : Eigen Value Vector [size (NEQ,Number of Eigenvalue Required]
[I] Denotes Input
[O] Denotes Output
Example for the data layout for the stiffness or mass matrix:
1 2 3 4 5 6
1 | 100 1 2 5 | | X1 | | 201 |
2 | 200 6 7 9 | | X2 | | 202 |
3 | 300 10 11 12 | | X3 | | 203 |
A = 4 | 400 13 14 | | X4 | = | 204 |
5 | 500 15 | | X5 | | 205 |
6 | 600 | | X6 | | 206 |
NUMBER OF EQUATIONS = 6
NUMBER OF COEFFICIENTS = 12
NEQ = 6
NCOEFF = 12
DIAG = { 100, 200, 300, 400, 500, 600}
KPT = { 3, 3, 3, 2, 1, 0}
KIND = { 2, 3, 6, 3, 4, 6, 4, 5, 6, 5, 6, 6}
AMAT = { 1, 2, 5, 6, 7, 9, 10,11,12, 13,14, 15}
RHS = { 201, 202, 203, 204, 205, 206}
|