On Tue, Oct 23, 2018 at 10:53 AM Manav Bhatia wrote:
> Really interesting!
>
> So this is a limitation of the algorithm and not the implementation.
>
> The challenge is that the eigenvalue solution in my workflow is a small
> component of a large computation done with real numbers in an
Really interesting!
So this is a limitation of the algorithm and not the implementation.
The challenge is that the eigenvalue solution in my workflow is a small
component of a large computation done with real numbers in an optimization
problem. I could do the whole thing with complex
> El 23 oct 2018, a las 16:10, Manav Bhatia escribió:
>
> Thanks for the clarification.
>
> Does this also apply to the standard non-hermitian eigenvalue problem? Do I
> need to compile with complex numbers if I want to capture the complex
> eigenvalues? Or does it work with real number
Thanks for the clarification.
Does this also apply to the standard non-hermitian eigenvalue problem? Do I
need to compile with complex numbers if I want to capture the complex
eigenvalues? Or does it work with real number support?
Thanks
Manav
Sent from my iPhone
> On Oct 23, 2018, at 3:43
If eigenvalues are complex then NLEIGS also needs to work in complex arithmetic
because it needs a region of the complex plane containing the wanted
eigenvalues. It seems that complex arithmetic is the only change in your
problem.
Jose
> El 22 oct 2018, a las 22:01, Manav Bhatia escribió:
>
It depends on the solver. For instance, NEPRII builds the matrix T(lambda) and
then uses it for matrix-vector multiplications and also for linear solves. So
the required operations depend on which preconditioner you use for the linear
solves. This example can use Jacobi preconditioner:
As a followup to this, if I am using a shell matrix for eigensolution (linear
or nonlinear eigenproblems), what operations should be defined for the shell
matrix?
-Manav
> On Oct 22, 2018, at 2:05 PM, Manav Bhatia wrote:
>
> Hi,
>
> I am exploring the nonlinear eigenvalue problem solver
> El 22 oct 2018, a las 21:05, Manav Bhatia escribió:
>
> Hi,
>
> I am exploring the nonlinear eigenvalue problem solver in Slepc.
>
> From the notes in "Sec 6.4: Retrieving the Solution”, it appears that if I
> expect to find complex eigenpairs then I must compile the library (and
Hi,
I am exploring the nonlinear eigenvalue problem solver in Slepc.
From the notes in "Sec 6.4: Retrieving the Solution”, it appears that if I
expect to find complex eigenpairs then I must compile the library (and Petsc)
with complex scalars. Is that correct?
Is there a way to