[petsc-users] Question about Nested matrices

[Alejandro D Otero]

Hi all,
I am trying to nest to square matrices to form a rectangular one with dims
2n x n (with n the size of the square matrices)
I want the two matrices to have interlaced rows so I passed the following 2
list of indices to MatCreateNest:

[[ 0  2  4  6  8 10 12 14 16 18 20 22 24 26 28 30 32 34 36 38 40 42 44 46
48]
[ 1  3  5  7  9 11 13 15 17 19 21 23 25 27 29 31 33 35 37 39 41 43 45 47
49]]

And we get the following error message:
Global sizes (25,25) of nested submatrix (0,0) do not agree with space
defined by index sets (50,25)

The size of the 2 square matrices is 25 x 25.

What am I doing wrong?

Thanks in advance,
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20150714/e15fbb28/attachment.html>