<p dir="ltr">Hi Jose and Hong,</p>
<p dir="ltr">Thanks so much for the detailed response. I am glad to know the reason behind it--hopefully we eventually figure out why the solvers have this behavior! Hong, I really appreciate you working on a patch to throw an error in this case. It definitely bit me and some people using my code... Hopefully it doesn't happen to anyone else! :)</p>
<p dir="ltr">Best,<br>
Greg</p>
<br><div class="gmail_quote"><div dir="ltr">On Mon, Sep 25, 2017, 7:44 AM Hong <<a href="mailto:hzhang@mcs.anl.gov">hzhang@mcs.anl.gov</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Greg and Jose,<div>I'll check this case, at least have petsc provide error message.</div><div>Please give me some time because I'm in the middle of several tasks.</div><div>I'll let you know once I add this support.</div><div><br></div><div></div></div><div dir="ltr"><div>Hong</div></div><div dir="ltr"><div><br><div class="gmail_extra"><br><div class="gmail_quote">On Mon, Sep 25, 2017 at 3:46 AM, Jose E. Roman <span dir="ltr"><<a href="mailto:jroman@dsic.upv.es" target="_blank">jroman@dsic.upv.es</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Greg,<br>
<br>
The linear solver table probably needs to be updated. I have tried several Cholesky solvers. With mkl_pardiso I get an explicit error message that it does not support Cholesky with complex scalars. The rest (PETSc, MUMPS, CHOLMOD) give a wrong answer (without error message). The problem is not related to your matrix, nor to shift-and-invert in SLEPc. I tried with ex1.c under PETSC_DIR/src/ksp/ksp/examples/tutorials. The example works in complex scalars, but the matrix is real. As soon as you add complex entries Cholesky fails, for instance adding this:<br>
  ierr = MatSetValue(A,0,1,1.0+PETSC_i,INSERT_VALUES);CHKERRQ(ierr);<br>
  ierr = MatSetValue(A,1,0,1.0-PETSC_i,INSERT_VALUES);CHKERRQ(ierr);<br>
<br>
I don't know if it is a bug or that the algorithm cannot support complex Hermitian matrices. Maybe Hong can confirm any of these. In the latter case, I agree that all packages should give an error message, as mkl_pardiso does.<br>
<br>
As a side comment, I would suggest using LU instead of Cholesky. Cholesky performs less flops but it does not mean it will be faster - I have seen many cases where it is slower than LU, maybe because in shift-and-invert computations the matrix is often indefinite, so an indefinite factorization is computed rather than Cholesky. Some SLEPc eigensolvers (e.g. LOBPCG) require that the preconditioner is symmetric (Hermitian), but the default solver (Krylov-Schur) is quite robust when you use LU instead of Cholesky in Hermitian problems. And you can always solve the problem as non-Hermitian (the difference in accuracy should not be too noticeable).<br>
<span class="m_-4754538946754754778HOEnZb"><font color="#888888"><br>
Jose<br>
</font></span><div class="m_-4754538946754754778HOEnZb"><div class="m_-4754538946754754778h5"><br>
<br>
> El 25 sept 2017, a las 7:17, Greg Meyer <<a href="mailto:gregory.meyer@gmail.com" target="_blank">gregory.meyer@gmail.com</a>> escribió:<br>
><br>
> Hi all,<br>
><br>
> Hong--looking at your link, there may be no special algorithm for Hermitian matrices in MUMPS, but that doesn't mean it can't solve them like it would any matrix. Furthermore it appears that Cholesky of complex matrices is supported from this link: <a href="https://www.mcs.anl.gov/petsc/documentation/linearsolvertable.html" rel="noreferrer" target="_blank">https://www.mcs.anl.gov/petsc/documentation/linearsolvertable.html</a><br>
><br>
> So do you or anyone have any idea why I get incorrect eigenvalues?<br>
><br>
> Thanks,<br>
> Greg<br>
><br>
> On Thu, Sep 21, 2017 at 5:51 PM Greg Meyer <<a href="mailto:gregory.meyer@gmail.com" target="_blank">gregory.meyer@gmail.com</a>> wrote:<br>
> Ok, thanks. It seems that PETSc clearly should throw an error in this case instead of just giving incorrect answers? I am surprised that it does not throw an error...<br>
><br>
> On Thu, Sep 21, 2017 at 5:24 PM Hong <<a href="mailto:hzhang@mcs.anl.gov" target="_blank">hzhang@mcs.anl.gov</a>> wrote:<br>
> Greg :<br>
> Yes, they are Hermitian.<br>
><br>
> PETSc does not support  Cholesky factorization for Hermitian.<br>
> It seems mumps does not support Hermitian either<br>
> <a href="https://lists.mcs.anl.gov/mailman/htdig/petsc-users/2015-November/027541.html" rel="noreferrer" target="_blank">https://lists.mcs.anl.gov/mailman/htdig/petsc-users/2015-November/027541.html</a><br>
><br>
> Hong<br>
><br>
><br>
> On Thu, Sep 21, 2017 at 3:43 PM Hong <<a href="mailto:hzhang@mcs.anl.gov" target="_blank">hzhang@mcs.anl.gov</a>> wrote:<br>
> Greg:<br>
><br>
> OK, the difference is whether LU or Cholesky factorization is used. But I would hope that neither one should give incorrect eigenvalues, and when I run with the latter it does!<br>
> Are your matrices symmetric/Hermitian?<br>
> Hong<br>
><br>
> On Thu, Sep 21, 2017 at 2:05 PM Hong <<a href="mailto:hzhang@mcs.anl.gov" target="_blank">hzhang@mcs.anl.gov</a>> wrote:<br>
> Gregory :<br>
> Use '-eps_view' for both runs to check the algorithms being used.<br>
> Hong<br>
><br>
> Hi all,<br>
><br>
> I'm using shift-invert with EPS to solve for eigenvalues. I find that if I do only<br>
><br>
> ...<br>
>   ierr = EPSGetST(eps,&st);CHKERRQ(ierr);<br>
>   ierr = STSetType(st,STSINVERT);CHKERRQ(ierr);<br>
> ...<br>
><br>
> in my code I get correct eigenvalues. But if I do<br>
><br>
> ...<br>
>   ierr = EPSGetST(eps,&st);CHKERRQ(ierr);<br>
>   ierr = STSetType(st,STSINVERT);CHKERRQ(ierr);<br>
>   ierr = STGetKSP(st,&ksp);CHKERRQ(ierr);<br>
>   ierr = KSPGetPC(ksp,&pc);CHKERRQ(ierr);<br>
>   ierr = KSPSetType(ksp,KSPPREONLY);CHKERRQ(ierr);<br>
>   ierr = PCSetType(pc,PCCHOLESKY);CHKERRQ(ierr);<br>
> ...<br>
><br>
> the eigenvalues found by EPS are completely wrong! Somehow I thought I was supposed to do the latter, from the examples etc, but I guess that was not correct? I attach the full piece of test code and a test matrix.<br>
><br>
> Best,<br>
> Greg<br>
><br>
<br>
</div></div></blockquote></div><br></div></div></div></blockquote></div>