Many thanks for your advice. Another little question: Are there alternate methods for renumbering the matrix (instead of MatGetSubMatrix())? Thanks again for your time. Pierre.