From: Michael S. <m-s...@us...> - 2005-03-25 13:07:45
|
On 25.03.05, Benjamin S. Kirk wrote: > A fairly arbitrary sparsity pattern can be computed if you use the > DofMap::dof_coupling array. By default, all DOFs are considered coupled > to each other, but it sounds like for your case (in 2D) you probably > want something like this: > > u v P > |0 0 1| > |0 0 1| > |1 1 0| I needed something like this for testing purposes some time ago. Because I was not interested in a very efficient program I let Petsc compute the sparsity pattern for me which is sub-optimal. I think, a framework for this is of broader interest. Best greetings, Michael -- "A mathematician is a device for turning coffee into theorems" Paul Erdös. |