On Tue, Jun 12, 2012 at 3:54 PM, Mark F. Adams <[email protected]>wrote:

How is what? The ex10 code or "I set the block size for matrices before SetFromOptions".

This how I create matrices with blocks sizes:

ierr = MatCreate( wcomm, &Prol ); CHKERRQ(ierr); ierr = MatSetSizes(Prol,nloc*bs,nLocalSelected*col_bs,PETSC_DETERMINE,PETSC_DETERMINE); CHKERRQ(ierr); ierr = MatSetBlockSizes( Prol, bs, col_bs ); CHKERRQ(ierr);

Sorry, I thought you were talking about vectors because the error is coming from a vector load. Where is this input file? I'll grab it sort out the issue.