[petsc-users] Error to TaoSolve() call

[Jed Brown]

Khai Hong Pham <khaipham at utexas.edu> writes:

> Jed,
>
> Thank you for your advise. I'm solving multi-field problem which were
> solved by static-splitting schema. I tried to reused my code, that's why I
> tried to use the nested matrix and vector. I'm trying to use the contiguous
> vector and matrix storage now. I appreciate if you could give me any advise
> on managing the matrix and vector (based on the skeleton of the existed
> code for separate problems) to get good performance.

I recommend using MatGetLocalSubMatrix() similar to
src/snes/examples/tutorials/ex28.c.  That way your old assembly should
work (provided you use local indices instead of global, but local
indices are better for this purpose anyway).
-------------- next part --------------
A non-text attachment was scrubbed...
Name: not available
Type: application/pgp-signature
Size: 818 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20140812/c49ddbac/attachment.pgp>