mmhh, i think i am missing something. Doesn't PCApply() apply the preconditioner to a vector? So how would that work (easily) with a matrix? matt On Thursday 09 October 2008, Matthew Knepley wrote: > CApply().