[petsc-users] Dirichlet boundary conditions for MATMPISBAIJ

[Matthew Knepley]

On Thu, Dec 4, 2014 at 7:14 AM, Aurelia Cuba Ramos <
aurelia.cubaramos at epfl.ch> wrote:

> Hi all,
>
> I recently started using PETSc and I am working with the MATMPISBAIJ
> matrix format. I am trying to apply Dirichlet boundary conditions by
> calling MatZeroRowsColumns() but this function seems to be not available
> for symmetric block matrices in parallel. I saw on the FAQ website that
> another possible way is to use MatZeroRows() and -ksp_type preonly
> -pc_type redistribute. But when I try calling MatZeroRows() I get again
> a PetscErrorCode 56, so it doesn't seem to be supported neither. I am
> currently using PETSc 3.4. Is MatZeroRowsColumns() avaibale for
> MATMPISBAIJ in PETSc 3.5 and if not, what is the fastest way to apply
> Dirichlet BC's for this matrix format?
>

1) SBAIJ is just an optimization to save memory. I would get everything
running
    with BAIJ, and then if you run out of memory switch to SBAIJ

2) We have only coded that up for the sequential version. I will add it to
our list.

  Thanks,

    Matt


> Many thanks,
>
> Aurelia
>



-- 
What most experimenters take for granted before they begin their
experiments is infinitely more interesting than any results to which their
experiments lead.
-- Norbert Wiener
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20141204/cf8298ca/attachment.html>