[petsc-users] Dirichlet boundary conditions for MATMPISBAIJ

[Aurelia Cuba Ramos]

Thank you Matt, I changed the format to BAIJ as you recommended. It
seems to work fine now.
 

On 04.12.2014 22:16, Matthew Knepley wrote:
> On Thu, Dec 4, 2014 at 7:14 AM, Aurelia Cuba Ramos
> <aurelia.cubaramos at epfl.ch <mailto: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/20141205/a7ed14ec/attachment.html>