<div dir="ltr"><div class="gmail_quote"><div dir="ltr">On Tue, Oct 23, 2018 at 10:53 AM Manav Bhatia <<a href="mailto:bhatiamanav@gmail.com">bhatiamanav@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Really interesting! <br>
<br>
So this is a limitation of the algorithm and not the implementation. <br>
<br>
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 numbers, but would be wasting a lot of cpu cycles doing that. <br>
<br>
Do you know if there is anyway to link to both real and complex versions of the library or switch between them at runtime? <br></blockquote><div><br></div><div>Unfortunately, that is a huge amount of work. One group has done it, but we do not support it.</div><div><br></div><div>  Thanks,</div><div><br></div><div>    Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
-Manav<br>
<br>
Sent from my iPhone<br>
<br>
> On Oct 23, 2018, at 9:42 AM, Jose E. Roman <<a href="mailto:jroman@dsic.upv.es" target="_blank">jroman@dsic.upv.es</a>> wrote:<br>
> <br>
> <br>
> <br>
>> El 23 oct 2018, a las 16:10, Manav Bhatia <<a href="mailto:bhatiamanav@gmail.com" target="_blank">bhatiamanav@gmail.com</a>> escribió:<br>
>> <br>
>> Thanks for the clarification. <br>
>> <br>
>> 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? <br>
> <br>
> No, linear eigenproblems (EPS) can be solved with real scalars for complex eigenvalues, but nonlinear eigenproblems (NEP) cannot.<br>
> <br>
> Jose<br>
> <br>
>> <br>
>> Thanks<br>
>> Manav<br>
>> <br>
>> Sent from my iPhone<br>
>> <br>
>>> On Oct 23, 2018, at 3:43 AM, Jose E. Roman <<a href="mailto:jroman@dsic.upv.es" target="_blank">jroman@dsic.upv.es</a>> wrote:<br>
>>> <br>
>>> 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.<br>
>>> <br>
>>> Jose<br>
>>> <br>
>>> <br>
>>>> El 22 oct 2018, a las 22:01, Manav Bhatia <<a href="mailto:bhatiamanav@gmail.com" target="_blank">bhatiamanav@gmail.com</a>> escribió:<br>
>>>> <br>
>>>> Thanks, Jose. <br>
>>>> <br>
>>>> How difficult would it be to add the support for the general case (if at all possible)? <br>
>>>> <br>
>>>> My eigenvalue problem is of the form shown in the attachment. Beta is the eigenvalue and X_s^\Delta is the eigenvector. While some of the matrices are known, others are defined only as matrix vector products. <br>
>>>> <br>
>>>> I am interested in eigenvalues with the largest real part. I expect to find complex eigenvalues, although for a small subset of cases these will be real. <br>
>>>> <br>
>>>> What is your recommendation for attacking this problem with the nonlinear eigenvalue support in Slepc? <br>
>>>> <br>
>>>> Would appreciate your guidance. <br>
>>>> <br>
>>>> Regards,<br>
>>>> Manav<br>
>>>> <br>
>>>> <PastedGraphic-1.pdf><br>
>>>> <br>
>>>> <br>
>>>>> On Oct 22, 2018, at 2:40 PM, Jose E. Roman <<a href="mailto:jroman@dsic.upv.es" target="_blank">jroman@dsic.upv.es</a>> wrote:<br>
>>>>> <br>
>>>>> <br>
>>>>> <br>
>>>>>> El 22 oct 2018, a las 21:05, Manav Bhatia <<a href="mailto:bhatiamanav@gmail.com" target="_blank">bhatiamanav@gmail.com</a>> escribió:<br>
>>>>>> <br>
>>>>>> Hi, <br>
>>>>>> <br>
>>>>>> I am exploring the nonlinear eigenvalue problem solver in Slepc. <br>
>>>>>> <br>
>>>>>> 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? <br>
>>>>>> <br>
>>>>>> Is there a way to include support for complex eigenpairs in a library complied with real scalars? <br>
>>>>>> <br>
>>>>>> Regards,<br>
>>>>>> Manav<br>
>>>>>> <br>
>>>>>> <br>
>>>>> <br>
>>>>> Currently, the only combination that supports complex eigenpairs with real scalars is the split form for the nonlinear function with the NLEIGS solver.<br>
>>>>> <br>
>>>>> Jose<br>
>>>> <br>
>>> <br>
> <br>
</blockquote></div><br clear="all"><div><br></div>-- <br><div dir="ltr" class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div>What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>-- Norbert Wiener</div><div><br></div><div><a href="http://www.cse.buffalo.edu/~knepley/" target="_blank">https://www.cse.buffalo.edu/~knepley/</a><br></div></div></div></div></div></div></div></div>