<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 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>