<div class="gmail_quote">On Fri, Aug 17, 2012 at 7:11 AM, Alexander Grayver <span dir="ltr"><<a href="mailto:agrayver@gfz-potsdam.de" target="_blank">agrayver@gfz-potsdam.de</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">

  
    
  
  <div bgcolor="#FFFFFF" text="#000000"><div class="im">
    <div>On 17.08.2012 14:52, Jed Brown wrote:<br>
    </div>
    <blockquote type="cite">
      <div class="gmail_quote"><span dir="ltr"></span><br>
        <div>See this <b>ADDITIVE</b>? You are preconditioning with
          I+P_{SOR}. You probably want multiplicative.</div>
      </div>
    </blockquote></div>
    Jed,<br>
    <br>
    It didn't work.<br>
    <br>
    -pc_type composite -pc_composite_pcs shell,sor -pc_composite_type
    multiplicative <br>
    -----><br>
    Linear solve did not converge due to DIVERGED_ITS iterations 250<br>
    <br>
    KSP Object:() 1 MPI processes<br>
      type: tfqmr<br>
      maximum iterations=250, initial guess is zero<br>
      tolerances:  relative=1e-07, absolute=1e-50, divergence=10000<br>
      left preconditioning<br>
      using PRECONDITIONED norm type for convergence test<div class="im"><br>
    PC Object:() 1 MPI processes<br>
      type: composite<br></div>
      Composite PC type - MULTIPLICATIVE<div class="im"><br>
      PCs on composite preconditioner follow<br></div>
    ....<br>
    <br>
    Is this something unexpected?<br></div></blockquote><div><br></div><div>Sorry, I misinterpreted what you want. PC_COMPOSITE_MULTIPLICATIVE re-evaluates the residual after applying the first stage of the preconditioner. Do you _really_ want to just apply one preconditioner after the next? You can use PC_COMPOSITE_SPECIAL, but read this man page and understand why your preconditioner has that form. <a href="http://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/PC/PCCompositeType.html">http://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/PC/PCCompositeType.html</a></div>
<div><br></div><div><br></div><div>Barry, it looks like despite the interface, "alpha" has never been used by this implementation.</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div bgcolor="#FFFFFF" text="#000000">
    I can try to write minimal example to make sure that this is not
    related to my app itself.<span class="HOEnZb"><font color="#888888"><br>
    <br>
    <pre cols="72">-- 
Regards,
Alexander</pre>
  </font></span></div>

</blockquote></div><br>