[petsc-users] Slepc: Nonlinear eigenvalue problem

[Jose E. Roman]

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: http://slepc.upv.es/documentation/current/src/nep/examples/tutorials/ex21.c.html

Jose


> El 22 oct 2018, a las 21:40, Manav Bhatia <bhatiamanav at gmail.com> escribió:
> 
> 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 <bhatiamanav at gmail.com> wrote:
>> 
>> 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 include support for complex eigenpairs in a library complied with real scalars? 
>> 
>> Regards,
>> Manav
>> 
>> 
>