<div dir="ltr"><span class="gmail-im">Dear PETSc users,<br>
<br>
we would like to compare our preconditioner for the Schur complement<br>
of a Stokes system, with the LSC preconditioner already implemented in<br></span>
PETSc. Following the example in the PETSc manual, we've tried<br>
   -fieldsplit_1_pc_type lsc -fieldsplit_1_lsc_pc_type ml<br>
but this is not working (properly) on our problem.<br>
<br>
On the other hand we think we have a good preconditioner for A10*A01,<br>
so  we'd like to try<br>
   -fieldsplit_1_pc_type lsc -fieldsplit_1_lsc_pc_type shell<br>
but we cannot figure out how to attach our apply() routine to<br>
the pc object of fieldsplit_1_lsc.<br>
<br>
Can this be done in the current interface?<br>
Or perhaps, should we call KSPGetOperators on the fieldsplit_1 solver<br>
and attach to its Sp operator a "LSC_Lp" of type MATSHELL with our routine<br><div>
attached to the MATOP_SOLVE of the shell matrix?</div><div><br></div><div>Thanks in advance,</div><div> <br></div><div>Elena and Matteo<br></div></div>

<br>
<div style="font-size:1.3em">------------------------</div><span style="font-size:small"><div><img src="http://www.unito.it/sites/default/files/logounito.png"></div><div><span style="font-size:small"><br></span></div>Indirizzo istituzionale di posta elettronica degli studenti e dei laureati dell'Università degli Studi di Torino</span><div><font size="2">Official University of Turin email address for students and graduates </font></div>