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
signature.asc
Description: Digital signature
_______________________________________________ DOLFIN-dev mailing list [email protected] http://www.fenics.org/mailman/listinfo/dolfin-dev
