# Tagged Questions

An eigenvector of an operator is a vector such that the action of the operator is the same as multiplication by a constant, called the eigenvalue. The eigensystem of an operator is the set of all such eigenvectors and their associated eigenvalues.

94 views

### finite difference frequency domain Eigenvalue matrix to get eigenmodes

I am using the FDFD method to calculate eigenmodes for an empty wave guide. It's a 1x1meter structure with a PEC boundary. (Here I have 6x6 points to make it simple). Sounds simple, but I can't get it ...
53 views

17 views

### Computing modes with fourier beam propagation

For my final optics projects I have spent the last few weeks implementing the beam proportionate method with Fourier split steps. This now works really well. Now I am trying to compute the modes of ...
32 views

### Learning Lanczos eigenvalue iteration - what is wrong?

I try to familiarize myself with iterative eigenvalue solvers such as Lanczos. So I tried rewrite it to python directly according to wiki. But it doesn't seem to work. The problem: it approximates ...
57 views

### Using DSYEV in LAPACK

I am calculating eigenvalues of a symmetric matrix using DSYEV in LAPACK.The concerned matrix is in a lower triangular form in a .dat file. I was searching for a simple explanation of the parameters ...
58 views

### Fastest way to diagonalize square matrices with c++

I am currently writing software in c++ which solves the eigenvalue problem of sparse hermitian matrices. The size of the matrix depends on the user input, but as an estimation it will roughly be ...
281 views

### Smallest eigenvalue without inverse

Suppose $A\in\mathbb{R}^{n\times n}$ is a symmetric, positive definite matrix. $A$ is big enough that it's expensive to solve $Ax=b$ directly. Is there an iterative algorithm for finding the ...
99 views

### Numerical computation of Perron-Frobenius eigenvector

I would like to efficiently and (to the extent possible) reliably find the Perron-Frobenius eigenvector of non-negative matrices. These are not stochastic matrices, they are typically dense, and their ...
41 views

### Memory requirement to find eigenvalues and -vectors of large sparse matrix

How can I estimate how much memory will be needed to find eigenvalues and eigenvectors of a given large sparse matrix? I have a real symmetric matrix with roughly $5 \times 10^4$ rows and columns, ...
36 views

### Rotate complex eigenvectors

I'm computing the eigenvectors of a real non-symmetric Matrix. I know that complex eigenvectors would come in conjugate. However, does anybody know a algorithm to rotate such complex eigenvectors ...
103 views

### Eigenvalue-style optimization with quadratic constraints

Suppose $A\in\mathbb{R}^{n\times n}$ is symmetric and positive definite and that we have several symmetric matrices $B_i\in\mathbb{R}^{n\times n}$ that are low-rank and indefinite. I need an ...
115 views

### Does this partial eigen-expansion have a name?

This question is a follow-up to this one. Let $A\in \mathbb{R}^{n\times n}$ be large, sparse, symmetric and positive definite. Suppose for I already know $m<n$ eigenpairs of $A$, corresponding to ...
101 views

### Julia: ordschur command

I just started leaning Julia, and my command of it is still very preliminary. I am trying to reorder a Schur decomposition. In Matlab, I could use the ordqz command and just had to specify the ...
78 views

146 views

### Appropriate iterative linear solver for an eigenvalue problem

I'm trying to solve a generalized eigenvalue problem $$Ax = \lambda Bx, \quad A = A^\top > 0,\; B = B^\top > 0$$ with $\lambda \approx \sigma$ using Rayleigh Quotient Iteration (RQI) (RQI is ...
58 views

### Is there an efficient $O(n^2)$ way to get the eigen decomposition given a LDL factorization?

Let's say I have a LDL factorization of a matrix A. Is there an efficient $O(n^2)$ way to get the eigen decomposition of A given it's LDL factorization? Is there a more efficient way, in case L and ...
285 views

### Verification in Eigenvalue problems

Let us start with a problem of the form $$(\mathcal{L} + k^2) u=0$$ with a set of given boundary conditions (Dirichlet, Neumann, Robin, Periodic, Bloch-Periodic). This corresponds with finding the ...
136 views

### How to correctly normalize modulus and phase of an eigenvector?

I am solving a linear stability problem using finite element discretization. Then, I have a generalised eigenvalue problem: $$\lambda M x = J x.$$ I obtain complex eigenvalue and eigenvectors from ...
18 views

61 views

### large eigenvalues with LAPACK

I have question about LAPACK. I calculate eigenvalues of a $16\times16$ Hermitian complex matrix with small entries by ZHEEV ...
252 views

### Armadillo eig_sym() for extracting eigenvalues. Is it parallel at all? [closed]

After wasting 3 days with scalapack, I gave up and moved to Armadillo, considering it uses lapack underneath its beatiful and easy interface. I would like to calculate the eigen values and eigen ...
123 views

109 views

### Finding Interior eigenvalues using Davidson algorithm

Is it possible to find interior eigenvalues closer to some lambda using Davidson method. I was searching online but found that most people use Jacobi-Davidson method for that. Thanks
80 views

### Computing eigendecomposition of a Hermitian matrix that is almost unitary

I have a dense Hermitian matrix that is approximately unitary, so it has eigenvalues that are $\sim \pm1$. I would like to compute all the eigenvectors corresponding to the $+1$ eigenvalue (not ...
278 views

### How to find the smallest positive eigenvalue of a large general system if they are all in +/- pairs of real eigenvalues

I used MKL Lapack SBGVX because it can solve for "selected" eigenvalues/modes (both positive and negative), thinking it would be efficient, but it is extremely slow when compared with Bathe's subspace ...
Let $\Sigma$ be a covariance matrix (e.g. symmetric positive definite). For arbitrary vectors $\epsilon$, I need to compute $\chi^2 \equiv \epsilon^\top\Sigma^{-1}\epsilon$, which I do using a ...