<html><head></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><br><div><div>On Jun 12, 2012, at 4:25 PM, Jed Brown wrote:</div><br class="Apple-interchange-newline"><blockquote type="cite"><div class="gmail_quote">On Tue, Jun 12, 2012 at 3:18 PM, Mark F. Adams <span dir="ltr"><<a href="mailto:mark.adams@columbia.edu" target="_blank">mark.adams@columbia.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0.8ex; border-left-width: 1px; border-left-color: rgb(204, 204, 204); border-left-style: solid; padding-left: 1ex; position: static; z-index: auto; ">
<div>This is the calling sequence in ex10:</div><div><br></div><div><div> ierr = VecCreate(PETSC_COMM_WORLD,&b);CHKERRQ(ierr);</div><div> ierr = VecSetFromOptions(b);CHKERRQ(ierr);</div><div> ierr = VecLoad(b,fd);CHKERRQ(ierr);</div>
</div><div><br></div><div>I set the block size for matrices before SetFromOptions. </div></blockquote><div><br></div><div>How is that? You don't even know the Vec type.</div></div>
</blockquote></div><br><div>How is what? The ex10 code or "I set the block size for matrices before SetFromOptions".</div><div><br></div><div>This how I create matrices with blocks sizes:</div><div><br></div><div><div> ierr = MatCreate( wcomm, &Prol ); CHKERRQ(ierr);</div><div> ierr = MatSetSizes(Prol,nloc*bs,nLocalSelected*col_bs,PETSC_DETERMINE,PETSC_DETERMINE); </div><div> CHKERRQ(ierr);</div><div> ierr = MatSetBlockSizes( Prol, bs, col_bs ); CHKERRQ(ierr);</div></div><div><br></div><div><br></div></body></html>