<div dir="ltr"><div>Thank you for the replay.</div><div><br></div><div>Now I have set the precondition on the command line with "-fieldsplit_1_lsc_pc_type" , but is it also possible to set it from within the
      code instead?<br></div><div>What is the equivalent code to obtain the effect of
      "-fieldsplit_1_lsc_pc_type mat" ?<font color="#888888"><br></font></div><div><font color="#888888"><br></font></div><div><font color="#888888">Elena<br></font></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">Il giorno mer 24 feb 2021 alle ore 06:04 Jed Brown <<a href="mailto:jed@jedbrown.org">jed@jedbrown.org</a>> ha scritto:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">If you've already attached a MatShell, you could presumably use -fieldsplit_1_lsc_pc_type mat to just call its MatMult.<br>
<br>
The sense I've gotten when I wrote PCLSC and was experimenting with these methods is that the main selling point of LSC (for most discretizations) is that it's more algebraic than the cheaper PCD methods.<br>
<br>
Elena Travaglia <<a href="mailto:elena.travaglia@edu.unito.it" target="_blank">elena.travaglia@edu.unito.it</a>> writes:<br>
<br>
> 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>
> 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>
> attached to the MATOP_SOLVE of the shell matrix?<br>
><br>
> Thanks in advance,<br>
><br>
> Elena and Matteo<br>
><br>
> -- <br>
> ------------------------<br>
><br>
><br>
><br>
> Indirizzo istituzionale di posta elettronica <br>
> degli studenti e dei laureati dell'Università degli Studi di TorinoOfficial <br>
> University of Turin email address for students and graduates <br>
</blockquote></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>