Dear PETSc users,<br /><br />I implemented PCCOMPOSITE preconditioner of type PC_COMPOSITE_ADDITIVE. It has two components:<br /><br /> - pc1: PCASM          1-level additive Schwarz preconditioning<br /> - pc2: PCSHELL      2-level additive Schwarz preconditioning<br /><br />so basically the whole preconditioner is a 2-level additive Schwarz preconditioner, where the second level is hand-made in a MATSHELL matrix structure. For PCASM, I use defalut options. I run tests with PCSHELL activated and deactivated.<br /><br />I use the default stopping criterion with this setting:<br /><br />ierr=KSPSetTolerances(ksp,1e-7,PETSC_DEFAULT,PETSC_DEFAULT,PETSC_DEFAULT);CHKERRQ(ierr);<br /><br />In the end I get these convergence results with 256 MPI tasks (1054729 total dofs) using KSPCG<br /><br /> PSASM only:<br /> ...<br />134 KSP preconditioned resid norm 4.680038240771e-06 true resid norm 5.417972183036e-07 ||r(i)||/||b|| 2.402197263846e-06<br />135 KSP preconditioned resid norm 3.681603361690e-06 true resid norm 4.091445069661e-07 ||r(i)||/||b|| 1.814047363013e-06<br />136 KSP preconditioned resid norm 2.938575711894e-06 true resid norm 3.280865690135e-07 ||r(i)||/||b|| 1.454656155040e-06<br /> - convergence detected<br /><br /> PCASM + PCSHELL<br />  ...<br /> 22 KSP preconditioned resid norm 6.610753500258e-05 true resid norm 4.725684013113e-06 ||r(i)||/||b|| 2.095253504927e-05<br /> 23 KSP preconditioned resid norm 1.173818907551e-04 true resid norm 2.984650385207e-06 ||r(i)||/||b|| 1.323321483881e-05<br /> 24 KSP preconditioned resid norm 4.409443889428e-05 true resid norm 1.959757043931e-06 ||r(i)||/||b|| 8.689086709369e-06<br /> - convergence detected<br /><br />The preconditioner, in both (1- and 2-level) versions, works fine, but the relative residual norm stopping criterion is not achieved. Observing the residual norms, I don't think the next iteration would achieve the requested 10^-7 tolerance.<br /><br />Can anyone explain me why?<br /><br />Thank you in advance,<br />Federico Marini
<HR><BR>
<TABLE cellpadding="10">
<TR><TD>
<A HREF="http://www.unimi.it/13084.htm?utm_source=firmaMail&utm_medium=email&utm_content=linkFirmaEmail&utm_campaign=5xmille">
<img alt="5xmilleUniMi" src="http://www.unimi.it/img/banner_190x85_5mille.jpg">
</A>
</TD>
<TD>
<I>
<FONT color="#FC7F01" face="Verdana,Arial,Helvetica,sans-serif" size="2">
Il 5 x mille</FONT> <FONT color="#535354" face="Verdana,Arial,Helvetica,sans-serif" size="2">alla nostra Università è un investimento sui giovani,<BR>
sui loro migliori progetti.
<P>
Sostiene la libera ricerca.<BR>
Alimenta le loro speranze nel futuro.
</P>
<P>
Investi il tuo <FONT color="#FC7F01" face="Verdana,Arial,Helvetica,sans-serif" size="2">
5 x mille</FONT> <FONT color="#535354" face="Verdana,Arial,Helvetica,sans-serif" size="2">sui giovani.
</FONT>
</P>
</I>
</TD>
</TR>
</TABLE>