[petsc-users] Eigenvalues always converge to zero when using slepc4py-complex

Yelyzaveta Velizhanina velizhaninae at gmail.com
Thu Oct 7 04:07:06 CDT 2021


Hello Lisandro,

I am testing EPS for a random non-Hermitian complex matrix in shift-invert mode and comparing with SciPy’s wrapper of Lapack – seems to be working perfectly. Thanks.

Regards,
Yelyzaveta.

From: Lisandro Dalcin <dalcinl at gmail.com>
Date: Wednesday, 6 October 2021 at 20:52
To: Yelyzaveta Velizhanina <velizhaninae at gmail.com>
Cc: Jose E. Roman <jroman at dsic.upv.es>, petsc-users at mcs.anl.gov <petsc-users at mcs.anl.gov>
Subject: Re: [petsc-users] Eigenvalues always converge to zero when using slepc4py-complex
The new builds are up. Jose added instructions in the SLEPc FAQ. Could you try again?

Regards,

On Mon, 4 Oct 2021 at 16:13 Yelyzaveta Velizhanina <velizhaninae at gmail.com<mailto:velizhaninae at gmail.com>> wrote:
I see. Thanks, much appreciated.

Best regards,
Yelyzaveta Velizhanina.

From: Jose E. Roman <jroman at dsic.upv.es<mailto:jroman at dsic.upv.es>>
Date: Monday, 4 October 2021 at 15:10
To: Yelyzaveta Velizhanina <velizhaninae at gmail.com<mailto:velizhaninae at gmail.com>>
Cc: petsc-users at mcs.anl.gov<mailto:petsc-users at mcs.anl.gov> <petsc-users at mcs.anl.gov<mailto:petsc-users at mcs.anl.gov>>
Subject: Re: [petsc-users] Eigenvalues always converge to zero when using slepc4py-complex
Conda supports complex scalars for petsc4py. However, this is not implemented in slepc4py. Lisandro is trying to get this fixed, so if no issues arise this will be available in a couple of days, with slepc4py-3.16.0.

Jose


> El 3 oct 2021, a las 22:45, Yelyzaveta Velizhanina <velizhaninae at gmail.com<mailto:velizhaninae at gmail.com>> escribió:
>
> Dear all,
>
> I am having a problem to get EPS run properly with PETSc and SLEPc build with scalar_value=complex. I am using petsc4py and slepc4py. Installed everything, including PETSc and SLEPc, with conda. While real scalar value build works well, when using the complex one, all the eigenvalues always converge to 0 for any matrix and any solver. I’ve tried running examples given in this repo https://github.com/myousefi2016/slepc4py as well - same outcome, only zero eigenvalues. I am running MacOSX BigSur.
>
> Will appreciate any help,
>
> Best regards,
> Yelyzaveta Velizhanina.
--
Lisandro Dalcin
============
Senior Research Scientist
Extreme Computing Research Center (ECRC)
King Abdullah University of Science and Technology (KAUST)
http://ecrc.kaust.edu.sa/
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20211007/03c61260/attachment.html>


More information about the petsc-users mailing list