# Questions tagged [eigenvalues]

Questions on symbolically or numerically determining the eigenvalues of matrices (Eigenvalues, Eigensystem) or differential equations (DEigenvalues, DEigensystem, NDEigenvalues, NDEigensystem) in Mathematica. Also includes determining the eigenvalues of differential equations with DSolve or NDSolve.

111 questions with no upvoted or accepted answers
Filter by
Sorted by
Tagged with
1k views

### What algorithm is Mathematica using to find the smallest eigenvalue so quickly?

My question is what kind of black magic is Mathematica doing to obtain the correct answer so quickly compared to other programming languages? Details: I've written a Mathematica notebook to find the ...
559 views

### Numeric Solution Hydrogen Atom

I'm trying to solve the Schrödinger equation for the hydrogen atom without made the variable separation of the polar and radial coordinate. It is my test code to extrapolate to another system with its ...
341 views

### DEigenvalues and NDEigenvalues return different values

In the following example, DEigenvalues and NDEigenvalues return different results despite having identical arguments. Does anyone know why? (I use Mathematica 11.3) ...
268 views

### How to compute eigenvalues of linear function (not matrix)?

How to compute eigenvalues of a known linear function? In Julia, there is a package https://jutho.github.io/LinearMaps.jl/dev/ to compute the matrix representation of given function, then we can ...
132 views

### Differing behavior of Eigenvalues and Eigensystem

With the update to v12.0, I seem to be getting different behavior of eigenvalues returned by Eigenvalues and Eigensystem (oddly, ...
167 views

### Memory usage for smallest eigenvalues

I have a bunch of hermitian matrices which are huge (of order 2^17 x 2^17) but extremely sparse so that, when I build the matrices, the usage of RAM is low (say of order 1 GB or similar). The ...
157 views

### Solving a matrix pencil (quadratic eigenvalue) problem with Mathematica

According to Wikipedia The matrix pencil of degree $\ell$ is the matrix-valued function defined on the complex numbers $L(k) = \sum_{i=0}^{\ell} k^{i} A_{i}$. Here $A_{\ell}$ are non-zero $n\times n$...
110 views

### Does current Mathematica capability allow solving PDEs (e.g. reaction diffusion equation) over surface of a mesh or surfaces

I came across this nice paper (https://arxiv.org/pdf/1605.01583.pdf) where the authors simulated a reaction diffusion system over the surface of a gecko, primarily to understand how various patterns ...
80 views

171 views

### Finding matrix in Krylov subspace (Lanczos method)

The Lanczos method for finding the smallest eigenvalue of a hermiteian matrix $H$ is based on the construction of a vector subspace (Krylov space) where one can build a matrix $H_{Krylov}$ which is ...
50 views

### Difference in computation speed of Eigensystem not to expectation for ParallelTable and Table

I was trying to compute the time it takes for Eigensystem to evaluate while being inside a ParallelTable, as it is well-known LAPACK subroutine has an inbuilt Parallelization to it. And the difference ...
502 views

### Orthogonal matrix decomposition of symmetric matrix?

If matrix mat is symmetric, we should be able to decompose it into eigenvalue matrix matJ and orthogonal matrix ...
226 views

### What real symmetric matrices of this type can Mathematica find symbolic eigenvalues for?

I'm working on a problem calculating symbolic eigenvalues of matrices that always have a very simple form: they are real and symmetric and usually sparse. They have two distinct symbolic parameters. (...
78 views

### First few smallest eigenvalues of a large dense symmetric matrix

I construct a large (say 2000x2000) matrix M whose entries are real random variables drawn from a certain distribution. Most of these values will be nonzero, so <...
108 views

### (Possibly bugs? ) Wrong results provided by DEigensystem

I was trying DEigensystem with the following code: ...
330 views

I refurbished my post in order to be more understandable. After computing simulations of Chladni patterns with Mathematica (see my previous topics), I finally went to practice. I realized my own ...
146 views

### Derivative of eigenvalues

I work with 4x4 Hermitian matrices (r). I want to calculate a derivative of a function f[t,r] (ff[t_,r_]=1/2*D[f[t,r],t]), where the function f depends on the absolute value of the eigenvalues of r. ...
84 views

### Interface points of NDEigensystem

When solving an eigenvalue problem with "NDEigensystem", e.g. a 1D Eigenvalue problem with the interval composed of different materials, which should be solved by the pure numerical method such as FEM,...
125 views

100 views

71 views

### What kind of performance should I expect out of Eigensystem using FEAST?

I'm numerically solving a time-independent Schrödinger equation using Eigensystem's FEAST method. It takes a lot longer than I ...
249 views

### Block diagonalizing a complex anti-symmetric matrix

I am going to evaluate the block diagonal form of few skew-matrices. When matrix elements are real I can simply follow the approach suggested in this thread which I have implemented that as ...
252 views

### Boundary value problem, multiple dimensional shooting, coupled eigenvalue problem

Following the one dimensional boundary value problem here, I would like to understand the easiest way to solve a BVP for a coupled system. In the 1D case, BVP can be converted to an initial value ...
134 views

### DEigenvalues with Robin B.C. sign problem

To find eigenvalues for $y''=\lambda y$ with robin boundary conditions on one end, and Dirichlet on the other end, I am getting correct value when robin B.C. on the right side, but when I flip things, ...
81 views

### Eigenvector matrix of a real positive matrix to be from $SO(n)$. How?

Suppose I have a positive real matrix of dimension $n$ and as such there exists a diagonalizing (rotation) matrix that belongs to $SO(n)$. How do I force Eigensystem...
1 vote
56 views

### Eigenvectors of a matrix (Solving and Plotting)

Given a nxn matrix h[k] ...
1 vote
33 views

### Wrong eigenvalues for 2D QHO using DEigensystem[]

I try to solve 2D Quantum Harmonic Oscillator using DEigensystem[] in Mathematica 13.0. Here is my code: ...
1 vote
82 views

### Complex eigenvalue of a symmetric matrix

I am working on a eigenvalue problem. I am using Eigensystem However, I am facing issues when I change the input parameters: I get complex eigenvalues. I couldn't understand the reason although my ...
1 vote
53 views

### Finding the (symbolic) eigenvalues of a matrix with some assumption

I have the following 3$\times$3 symmetric matrix with symbolic entries. Meff = \begin{array}{ccc} 2 \beta \text{c12}^2 \text{c13}^2 \text{$\eta$11} \text{m1}+\alpha \text{c13}^2 \text{D31} \text{...
1 vote
71 views

### Why does the minimum eigenvalue sharply decrease when the number of basis functions increases from 37 to 38 and more?

I have the following Hamiltonian: H = -1/2 * Laplacian -1/r - 2 * r/5 * (Exp[-r * 1.6] + Exp[-r * 3.1]) I'm trying to find the minimum eigenvalue of this Hamiltonian using the matrix method. As basis ...
1 vote
50 views

### Calculating eigenvalues and solving inequalities with parameters

I would like to calculate b in terms of c and a, in order to satisfy the following conditions: ...
1 vote
28 views

### Why does Sparsearray diagonalization by Krylov technique become slower for very large sizes?

Here is an example sparse matrix, taken from example Hamiltonian, which was a question I asked before in this forum. ...
1 vote
203 views

### Intel MKL ERROR: Parameter 2 was incorrect on entry to ZGEHD2

I have some code which is diagonalising very large random matrices, and storing certain statistical properties of their eigenvalues and eigenvectors. The code runs fine on my computer. When running it ...
1 vote
42 views

### Eigenvalues broken for nonsymmetric matrices

It looks that eigensystem calculation can easily go broken with real nonsymmetric/general complex matrices. Here are two examples with the possibly complex spectra plotted in the complex plane. I'm ...
1 vote
Background I am using NDEigensystem to solve the following eigenvalue problem:  \left( \begin{matrix} m&-i\partial_x \\ -i\partial_x & -m\end{matrix}\right) \left( \begin{matrix} u_u(x) \\ ...