dropping columns/rows from matrix?

Peter Schröder ps at cs.caltech.edu
Fri Nov 23 20:13:20 CST 2007

Barry Smith wrote:
>    We use MatGetSubMatrix() to do this. Rather than dropping/adding 
> rows and columns just grab the
> part you want and use it then destroy it and grab another part you 
> want. We've done this with active set
> methods and the time to do the MatGetSubMatrix() has been surprisingly 
> small percentage of the time.
Aaaah. Ok. Even if the submatrix is everyone but one column/row? Very well.

The one other way I had thought might be to fix one of the variables 
(some entry of the x vector is forced to zero) and ignore the 
corresponding part of the residual. But this may be too hacky an 
approach given the overall structure of Petsc.

I'll use MatGetSubMatrix() for now.


More information about the petsc-users mailing list