<html>
  <head>
    <meta content="text/html; charset=utf-8" http-equiv="Content-Type">
  </head>
  <body bgcolor="#FFFFFF" text="#000000">
    Thank you Matt, I changed the format to BAIJ as you recommended. It
    seems to work fine now. <br>
     <br>
    <br>
    <div class="moz-cite-prefix">On 04.12.2014 22:16, Matthew Knepley
      wrote:<br>
    </div>
    <blockquote
cite="mid:CAMYG4GmpYqgn9pBsyaXnLSXNo90mcbgtWuV2eU++zDo8TBFn3w@mail.gmail.com"
      type="cite">
      <meta http-equiv="Content-Type" content="text/html; charset=utf-8">
      <div dir="ltr">
        <div class="gmail_extra">
          <div class="gmail_quote">On Thu, Dec 4, 2014 at 7:14 AM,
            Aurelia Cuba Ramos <span dir="ltr"><<a
                moz-do-not-send="true"
                href="mailto:aurelia.cubaramos@epfl.ch" target="_blank">aurelia.cubaramos@epfl.ch</a>></span>
            wrote:<br>
            <blockquote class="gmail_quote" style="margin:0 0 0
              .8ex;border-left:1px #ccc solid;padding-left:1ex">Hi all,<br>
              <br>
              I recently started using PETSc and I am working with the
              MATMPISBAIJ<br>
              matrix format. I am trying to apply Dirichlet boundary
              conditions by<br>
              calling MatZeroRowsColumns() but this function seems to be
              not available<br>
              for symmetric block matrices in parallel. I saw on the FAQ
              website that<br>
              another possible way is to use MatZeroRows() and -ksp_type
              preonly<br>
              -pc_type redistribute. But when I try calling
              MatZeroRows() I get again<br>
              a PetscErrorCode 56, so it doesn't seem to be supported
              neither. I am<br>
              currently using PETSc 3.4. Is MatZeroRowsColumns()
              avaibale for<br>
              MATMPISBAIJ in PETSc 3.5 and if not, what is the fastest
              way to apply<br>
              Dirichlet BC's for this matrix format?<br>
            </blockquote>
            <div><br>
            </div>
            <div>1) SBAIJ is just an optimization to save memory. I
              would get everything running</div>
            <div>    with BAIJ, and then if you run out of memory switch
              to SBAIJ</div>
            <div><br>
            </div>
            <div>2) We have only coded that up for the sequential
              version. I will add it to our list.</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">
              Many thanks,<br>
              <br>
              Aurelia<br>
            </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>
    </blockquote>
    <br>
  </body>
</html>