[petsc-users] MATNEST with shell matrices

[Jed Brown]

On Thu, Feb 21, 2013 at 5:48 PM, Amneet Bhalla <mail2amneet at gmail.com>wrote:

> Jed, so to use VecGetSubVector() with a regular Vec; it would need to set
> Vec->map->bs = no_of_components
> that I want to solve for? I understand that a wrapper code needs to be
> written for gelling SAMRAI data field and PETSc Vec; I am referring to such
> vector wrapper here.
>


No, you should never touch vec->map->bs directly in any circumstance and
it's totally wrong semantically here. You just need to define index sets
for the pieces.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20130221/4057dad4/attachment.html>