<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Sun, Dec 28, 2014 at 11:02 AM, Umut Tabak <span dir="ltr"><<a href="mailto:u.tabak@tudelft.nl" target="_blank">u.tabak@tudelft.nl</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
  
    
  
  <div text="#000000" bgcolor="#FFFFFF"><span class="">
    <div>On 12/28/2014 05:54 PM, Alp Kalpalp
      wrote:<br>
    </div>
    <blockquote type="cite">
      
      <div dir="ltr">
        <div>
          <div>
            <div>
              <div>
                <div>Hi, <br>
                  <br>
                </div>
                Thank you Mark. <br>
                <br>
              </div>
              Let me clarify my questions;<br>
              <br>
            </div>
            1-)How to implement or activate a Reorthogonalization
            procedure for KSPCG..<br>
          </div>
          <div>As you know, search directions can be found more rapidly
            (with less numer of iterations) by using previous successive
            directions<br>
          </div>
        </div>
      </div>
    </blockquote></span>
    Without answering the PETSc related questions, interesting
    discussion, <br>
    <br>
    indeed, but at the cost of purging the previous directions(which
    means explicit orthogonalizations with respect to these vectors
    also), so I am not sure if you can gain something with this, cost
    wise...</div></blockquote><div><br></div><div>This has been proposed many times, but it has never been shown to work. I have tried every variant I could</div><div>find and it did not work. You can try LGMRES, which is the closest one to working in my opinion. There is</div><div>definitely no theoretical relation between Krylov directions from subsequent solves unless the operator is</div><div>identical.</div><div><br></div><div>   Matt</div><div><br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div text="#000000" bgcolor="#FFFFFF"><span class=""><blockquote type="cite">
      <div dir="ltr">
        <div>
          2-) How to implement or activate a projection space over CG. A
          sample projection can be;<br>
          P = I - G*((G'*G)\G'). <br>
        </div>
        I need to insert project,scale,precondition,re-scae,re-project
        steps during each KSPCG iteration. How can I utilize this?<br>
        <div><br>
        </div>
      </div>
    </blockquote></span>
    Just a side note, I had previous experience on this that these kinds
    of practice increase the cost more...<br>
    BR,<br>
    Umut<br>
  </div>

</blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature">What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>-- Norbert Wiener</div>
</div></div>