Would there be any objections to changing the dof indices used in DofMap from global to local indices (with a local-to-global map if required)? The motivation is to provide good and transparent block matrix and block preconditioning support. Assembling PETSc block matrices without any modification of the assembly code requires a switch to local indices (this is because local-to-global maps are cheap, whereas global-to-local is expensive). I'm not aware of any effect that this would have on user code.
Garth _______________________________________________ fenics mailing list [email protected] http://fenicsproject.org/mailman/listinfo/fenics
