Do you use the same "mass matrix" in both cases, meaning do you call
eig(A, M) where A is the stiffness matrix and M is the mass matrix in
matlab?
You will need to compute the spectrum of the eigenvalue problem
A x = \lambda M x
Best,
Matthias
On Thu, May 30, 2024, at 10:55 CDT, Najwa Alsh
On 5/30/24 09:55, Najwa Alshehri wrote:
The problem arises when comparing the inf-sup constants obtained using Matlab
and dealii (solving on the same mesh ):
1. Matlab: Using the "eig" solver, the inf-sup constant decays to zero as the
mesh is refined.
2. deal.II: Using the "ArpackSolve
Dear all,
Thank you for your always support and help.
I am encountering an issue while solving an eigenvalue problem related to
the computation of the discrete inf-sup constant for a saddle point
problem. Specifically, I am solving the following system:
Bt A^-1 B eigenvector = eigenvalue
