functionwrapper_xBase
Compute the eigenvector bases according to Kautsky algorithm
Inputs
| Type | Name | Default | Description |
|---|---|---|---|
| Real[:,size(A, 1)] | A | Real square system matrix | |
| Real[size(A, 1),:] | B | Real input matrix | |
| Real[size(A, 1)] | gamma_real | Eigenvalue vector, real part | |
| Real[size(A, 1)] | gamma_imag | Eigenvalue vector, imaginary part | |
| Integer | ncp | Number of complex pairs of eigenvalues |
Outputs
| Type | Name | Default | Description |
|---|---|---|---|
| Real[size(A, 1),size(B, 2)] | U0 | U0 | |
| Real[size(B, 2),size(B, 2)] | Z | Z | |
| Real[size(A, 1),size(B, 2)*(size(A, 1) - ncp)] | S_real | Eigenvector bases, real part | |
| Real[size(A, 1),size(B, 2)*(size(A, 1) - ncp)] | S_imag | Eigenvector bases, imaginary part | |
| Integer | rankB | Rank of matrix B |