<div class="gmail_quote">On Fri, Jun 3, 2011 at 16:07, Barry Smith <span dir="ltr"><<a href="mailto:bsmith@mcs.anl.gov">bsmith@mcs.anl.gov</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
Unfortunately PETSc is restricted to working with only real or complex objects at one time and one cannot mix some real and some complex.</blockquote></div><br><div>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.</div>