dolfin team mailing list archive

[Anders Logg <logg@xxxxxxxxx>]

On Mon, Jul 07, 2008 at 09:38:15AM +0100, Garth N. Wells wrote:
> I've had a look at the code for applying Dirichlet boundary conditions 
> by eliminating constrained dofs, and I think that we need to do it 
> differently. A constrained dof should be indicated by a negative index 
> in the DofMap, and then eliminated during the main assembly loop. We 
> then need to do something for DiscreteFunction so that it can pick up 
> the boundary values, which should be pretty easy.

PETSc handles this (skipping insertion for negative indices).

http://www-unix.mcs.anl.gov/petsc/petsc-as/snapshots/petsc-current/docs/manualpages/Mat/MatSetValues.html

Does Epetra do the same?

-- 
Anders

Attachment: signature.asc
Description: Digital signature