<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Wed, May 21, 2014 at 2:09 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>So I just pulled an updated version of
      petsc-dev today (I switched from the <b>next</b> branch to the <b>master</b>
      branch due to some compilation error existing with the last commit
      on <b>next</b>).<br>
      I still have the same error and I believe this is the whole error
      message I have.<br>
      I mean I am running multiple time steps for my simulation so I
      have the same message at each time step, but I don't think that it
      is important to report these duplicates, is it?<br></div></div></blockquote><div><br></div><div>Sorry I am just getting back to this. I asked for more of the error message, because I assumed that some other</div><div>
PETSc function was calling DMGetGlobalVector(), not you, since you say it goes away with a different PC. If</div><div>you are calling it and its generating this error, it is because you are calling it before setting a Vec in the DMShell.</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>
      On 05/20/2014 09:14 PM, Matthew Knepley wrote:<br>
    </div>
    <blockquote type="cite">
      <div dir="ltr"><br>
        <div class="gmail_extra">
          <div class="gmail_quote">On Tue, May 20, 2014 at 4:33 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"> <font color="#ff0000"><font color="#000000">Hi all,<br>
                    I am running an FEM simulation that uses Petsc as a
                    linear solver.<br>
                    I am setting up ISs and pass them to a DMShell in
                    order to use the FieldSplit capabilities of Petsc.<br>
                    <br>
                    When I pass the following options to Petsc:<br>
                  </font></font>
                <blockquote>" -ksp_type gmres -pc_type fieldsplit
                  -pc_fieldsplit_type schur
                  -pc_fieldsplit_schur_factorization_type full
                  -pc_fieldsplit_schur_precondition selfp
                  -pc_fieldsplit_0_fields 1,2 -pc_fieldsplit_1_fields 0
                  -fieldsplit_0_ksp_type preonly -fieldsplit_0_pc_type
                  ilu -fieldsplit_Field_0_ksp_type gmres
                  -fieldsplit_Field_0_pc_type mg -malloc_log mlog
                  -log_summary time.log"<br>
                </blockquote>
                I get an error message:<br>
                <font color="#ff0000"><br>
                  [0]PETSC ERROR: --------------------- Error Message
                  --------------------------------------------------------------</font><br>
                [0]PETSC ERROR:   <br>
                [0]PETSC ERROR: Must call DMShellSetGlobalVector() or
                DMShellSetCreateGlobalVector()<br>
                [0]PETSC ERROR: See <a href="http://http://www.mcs.anl.gov/petsc/documentation/faq.html" target="_blank">http://http://www.mcs.anl.gov/petsc/documentation/faq.html</a>
                for trouble shooting.<br>
                [0]PETSC ERROR: Petsc Development GIT revision:
                v3.4.4-5071-g1163a46  GIT Date: 2014-03-26 22:20:51
                -0500<br>
                [0]PETSC ERROR:
                /home/luc/research/feap_repo/ShearBands/parfeap-dev/feap
                on a arch-linux2-c-opt named euler by luc Tue May 20
                11:31:11 2014<br>
                [0]PETSC ERROR: Configure options --download-cmake
                --download-hypre --download-metis --download-mpich
                --download-parmetis --with-debugging=no
                --with-share-libraries=no<br>
                [0]PETSC ERROR: #1 DMCreateGlobalVector_Shell() line 245
                in /home/luc/research/petsc/src/dm/impls/shell/dmshell.c<br>
                [0]PETSC ERROR: #2 DMCreateGlobalVector() line 669 in
                /home/luc/research/petsc/src/dm/interface/dm.c<br>
                [0]PETSC ERROR: #3 DMGetGlobalVector() line 154 in
                /home/luc/research/petsc/src/dm/interface/dmget.c<br>
              </div>
            </blockquote>
            <div><br>
            </div>
            <div>Always always always give the entire error message.</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"> I am not really
                sure why this happens but it only happens when
                -fieldsplit_Field_0_pc_type mg, with other
                preconditioners, I have no problems. I attached the
                ksp_view in case that's any use.<span class="HOEnZb"><font color="#888888"><span><font color="#888888"><br>
                    <pre cols="72">-- 
Best,
Luc</pre>
                  </font></span></font></span></div><span class="HOEnZb"><font color="#888888">
            </font></span></blockquote><span class="HOEnZb"><font color="#888888">
          </font></span></div><span class="HOEnZb"><font color="#888888">
          <br>
          <br clear="all">
          <div><br>
          </div>
          -- <br>
          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
        </font></span></div>
      </div>
    </blockquote>
    <br>
  </div>

</blockquote></div><br><br clear="all"><div><br></div>-- <br>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>