[petsc-users] ILU preconditioner

[Danesh Daroui]

Yes that's exactly what I want to do. I can also use PETSc to calculate
the preconditioner, but can I use it to calculate real part's ILU while
I am using complex solver?



On Fri, 2011-06-03 at 16:12 +0200, Jed Brown wrote:
> On Fri, Jun 3, 2011 at 16:07, Barry Smith <bsmith at mcs.anl.gov> wrote:
>         Unfortunately PETSc is restricted to working with only real or
>         complex objects at one time and one cannot mix some real and
>         some complex.
> 
> If I understand correctly, Danesh intends to use an external package
> (MKL) for the ILU on the real part. That will work fine with PCShell.
> His setup function will pull out the real part of the matrix, put it
> into an MKL data structure, and do the factorization. Then the apply
> function will take the input vector apart, solve with the MKL
> factorization, and put the result back together in the output vector.
> A bit awkward perhaps, but it should work fine.