petsc and x-fem

2008-08-13 Thread Farshid Mossaiby
Hi, I always liked idea of a *simple* parallel FEM code in Petsc, which shows most of tricks used in various programs, like parallel assembly of elements and imposing Dirichlet BCs in a proper way. I think this can eliminate many of How should I do that questions from newbies like me. I think

petsc and x-fem

2008-08-13 Thread Barry Smith
On Aug 13, 2008, at 1:50 AM, Farshid Mossaiby wrote: Hi, I always liked idea of a *simple* parallel FEM code in Petsc, which shows most of tricks used in various programs, like parallel assembly of elements and imposing Dirichlet BCs in a proper way. I think this can eliminate many

PCApply_Shell ghosts...

2008-08-13 Thread Eric Chamberland
Hi, we are using PETSc in our code and we have a problem with, I think, the ghosted values that we expect. We developed our own pre-conditioner, and try to use it in a parallel environment. With other precondtioners (PETSc built-in), everything works fine. But with our home-made one, here

PCApply_Shell ghosts...

2008-08-13 Thread Matthew Knepley
PETSc knows nothing about your domain, so it can't know what you might want ghosted. I think the thing to do is use a VecScatter to map the input Vec to a ghosted Vec (called a local vector in DA language). Matt On Wed, Aug 13, 2008 at 3:41 PM, Eric Chamberland Eric.Chamberland at