On Tue, Jun 12, 2012 at 3:18 PM, Mark F. Adams <mark.adams@columbia.edu> wrote:

This is the calling sequence in ex10:

ierr = VecCreate(PETSC_COMM_WORLD,&b);CHKERRQ(ierr);
ierr = VecSetFromOptions(b);CHKERRQ(ierr);
ierr = VecLoad(b,fd);CHKERRQ(ierr);

I set the block size for matrices before SetFromOptions.

How is that? You don't even know the Vec type.