Khai Hong Pham <[email protected]> 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).
pgphJ0JcTygmg.pgp
Description: PGP signature
