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 > <[email protected] <mailto:[email protected]>> 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
