On Jun 12, 2012, at 4:56 PM, Jed Brown wrote:

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.

There is a vector load, at least ex10 loads a vec and its seems to work (so I assume Blaise put one in there).

Where is this input file? I'll grab it sort out the issue.

You can get an archive containing the KSP and the data files at https://filestogeaux.lsu.edu/public/download.php?FILE=bourdin/52338Q0at2Y These are just standard "binaryoutput" files I think, Mark