On Jun 12, 2012, at 4:25 PM, Jed Brown wrote: > On Tue, Jun 12, 2012 at 3:18 PM, Mark F. Adams <mark.adams at 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.
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); -------------- next part -------------- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20120612/b51bdfc1/attachment.html>
