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.
