<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Fri, Nov 21, 2014 at 3:12 PM, Luc Berger-Vergiat <span dir="ltr"><<a href="mailto:lb2653@columbia.edu" target="_blank">lb2653@columbia.edu</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>Yes good call, I pulled master and got
      the patch.<br>
      I don't have errors anymore.<br></div></div></blockquote><div><br></div><div>+1 Lawrence</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"><div>
      Do this mean that when a dm is active, the dm is assumed to create
      hold the information about the vectors (residual, solution and
      internal work vectors) for the problem, whereas when the dm is not
      active or is not used, the ksp is assumed to hold the information
      about vectors format?<br></div></div></blockquote><div><br></div><div>Yes, exactly.</div><div><br></div><div>We wanted a way to have a DM associated with a problem, and yet have the Vec/Mat passed in</div><div>define it. We were trying to force the DM organization down the throat of users.</div><div><br></div><div>  Thanks,</div><div><br></div><div>    Matt</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"><div>
      <pre cols="72">Best,
Luc</pre><div><div class="h5">
      On 11/21/2014 12:08 PM, Lawrence Mitchell wrote:<br>
    </div></div></div><div><div class="h5">
    <blockquote type="cite">
      
      <div><br>
      </div>
      <div>On 21 Nov 2014, at 16:55, Luc Berger-Vergiat <<a href="mailto:lb2653@columbia.edu" target="_blank">lb2653@columbia.edu</a>>
        wrote:<br>
        <br>
      </div>
      <blockquote type="cite">Sure, so with gamg (which was my first
        choice too), I get the following:<br>
        <br>
        <font face="Courier New, Courier, monospace">[0]PETSC ERROR:
          --------------------- Error Message
          --------------------------------------------------------------<br>
          [0]PETSC ERROR:   <br>
          [0]PETSC ERROR: Must call DMShellSetGlobalVector() or
          DMShellSetCreateGlobalVector()<br>
          [0]PETSC ERROR: See <a href="http://www.mcs.anl.gov/petsc/documentation/faq.html" target="_blank">http://www.mcs.anl.gov/petsc/documentation/faq.html</a>
          for trouble shooting.<br>
          [0]PETSC ERROR: Petsc Release Version 3.5.2, Sep, 08, 2014 <br>
          [0]PETSC ERROR:
          /home/luc/research/feap_repo/ShearBands/parfeap/feap on a
          arch-opt named euler by</font></blockquote>
      <br>
      <div>I wonder if this is the same problem described here <a href="https://bitbucket.org/petsc/petsc/pull-request/221" target="_blank">https://bitbucket.org/petsc/petsc/pull-request/221</a></div>
      <div><br>
      </div>
      <div>Where the ksp contains an inactive dmshell that doesn't know
        how to create vectors </div>
      <div><br>
      </div>
      <div>Lawrence</div>
    </blockquote>
    <br>
  </div></div></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>