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>
