MatStencil only makes sense if you are using a distributed grid (DA), where it corresponds to physical field locations. You probably just want MatSetValuesBlocked ( http://www.mcs.anl.gov/petsc/petsc-as/snapshots/petsc-current/docs/manualpages/Mat/MatSetValuesBlocked.html )
Warm Regards, Aron On Mon, Dec 27, 2010 at 10:09 AM, Raeth, Peter <PRaeth at hpti.com> wrote: > Am a new PETSc user trying to make use of the matrix level to create > and operate on matrices whose memory exceeds that available on any one node. > To populate a distributed dense matrix with results of other matrix > calculations we are trying to use MatSetValuesBlockedStencil. Two of the > inputs to that function require structures of type MatStencil. After > searching the archives, tutorials, examples, and Google, I can not find > anything that explains what the values of MatStencil are meant to do > relative to where in the target matrix to place the block of values. > Would someone please point me in the right direction? > > Thanks, > > Peter. > > -------------- next part -------------- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20101227/12ac82bf/attachment.htm>
