<div dir="ltr"><div><div><div><div><div><div><div><div>Thanks for the answers,<br><br></div>Please forgive me, I forgot to say that my stiffness matrix is not changing during time steps. I could not remember directly but just after a google search..I just hit this <br><br><a href="http://web.stanford.edu/group/frg/publications/recent/FETI-stoch.pdf">http://web.stanford.edu/group/frg/publications/recent/FETI-stoch.pdf</a><br><br></div><div>please look around eq37<br></div><div><br></div>My problem is not related to this random paper I found. But, I think I can find several others that shows the enhancing power of orthogonalization with successive directions when the system's behaviour is not changing rapidly. In my current sample case a gradually increasing force is applied to a linear system.<br><br></div>Since I use FETIDP, <span class="im">preconditioned conjugate projected gradient (PCPG)</span> is crucial in order to select any generalized inverse for the system.<br><br></div>So, any suggestions on how to complete these tasks?<br><br></div>For example anyway of obtaining search direction from KSPCG?<br><br></div>or <br><br></div>how to implement a projection space?<br><br></div>Is it posible or too difficuly to code a variant of a KSPCG that meets my requirements?<br></div><div class="gmail_extra"><br><div class="gmail_quote">On Sun, Dec 28, 2014 at 7:08 PM, Matthew Knepley <span dir="ltr"><<a href="mailto:knepley@gmail.com" target="_blank">knepley@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div class="gmail_extra"><div class="gmail_quote"><span class="">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>
    <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></span><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><span class=""><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><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></span></div><span class="HOEnZb"><font color="#888888"><br><br clear="all"><div><br></div>-- <br><div>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>
</font></span></div></div>
</blockquote></div><br></div>