Feng Shi <[email protected]> writes: > Hi Folks, > > I get hard times in using the routine MatSetValuesStencil for > dof>1. There numerous examples for dof=1 but dof>1, neither in the > Petsc manual. The routine on-line manual > (http://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/Mat/MatSetValuesStencil.html) > says, idxm/idxn are grid coordinates (and component number when dof > 1) for > matrix columns being entered. What is the general form of these indices?
If you have dof>1, please use MatSetValuesBlockedStencil. See src/snes/examples/tutorials/ex48.c for an example usage for a finite-element method. > Am I right, for example, for idxn, if I write idxn[ix, dof].i? That's not valid syntax, so I don't know what you're asking.
pgp0wtwbry5Vc.pgp
Description: PGP signature
