Hi,
I need to verify if a matrix is positive (semi)-definite, ideally without computing the eigenvalues. Is there a standard way?
Thanks in advance!
Hi,
I need to verify if a matrix is positive (semi)-definite, ideally without computing the eigenvalues. Is there a standard way?
Thanks in advance!
I don’t know of any without computing the lowest eigenvalue.
Sorry, can you elaborate on how one could do that. I read the documentation and it is not very clear to me.
I forgot to ask whether your matrix is symmetric or not, but here is an example for Stokes equation FreeFem-sources/examples/hpddm/stokes-2d-SLEPc.edp at 2b2052ae9ec86ae7f0048607b20084234e03858f · FreeFem/FreeFem-sources · GitHub (determining whether the FE pair is inf-sup stable using the lowest eigenvalue).