Problem creating a non-square MPIAIJ Matrix

Hi there, I am found a very mysterious problem of creating a MxN matrix when MN. Please take a look at the attached short test code I wrote to demonstrate the problem. Basically, I want to create a 8x6 matrix. A={1 2 3 4 0 0; 0 0 0 0 0 0; 0 0 0 0 0 0; 0 0 0 0 0 0; 0 0 0 0 0 0; 0

Problem creating a non-square MPIAIJ Matrix

On Fri, 16 Feb 2007, Shi Jin wrote: I actually used MatGetLocalSize(A,m,n) in the code. They give me m=4,n=3, as expected. I can also specify m=4,n=3 in MatCreateMPIAIJ() which is exactly identical to the previous code. If I specify anything else, I get error saying that they don't agree

Problem creating a non-square MPIAIJ Matrix

We are stoing the diagonal block and offdiagonal block separately. However both blocks are on the same processor. i.e each processor stores m*N values - in 2 submatrices m*n, m*(N-n). To understand this better - check manpage for MatCreateMPIAIJ(). Thanks. But this is completely different

Problem creating a non-square MPIAIJ Matrix

On Fri, 16 Feb 2007, Shi Jin wrote: We are stoing the diagonal block and offdiagonal block separately. However both blocks are on the same processor. i.e each processor stores m*N values - in 2 submatrices m*n, m*(N-n). To understand this better - check manpage for MatCreateMPIAIJ().