Dear all,
I try to create a compatible sparse MPI matrix A with dmplex global vector x, so I can do matrix-vector multiplication y = A*x. I think I can first get the local and global sizes of x on comm, say n and N, also sizes of y, m, M, then create A by using MatCreate(comm, &A), set the sizes using MatSetSizes(A, m, n, M, N), set the type using MatSetType(A, MATMPIAIJ). Is this process correct? Another question is: Do the entries not filled automatically compressed out? Thanks! Junming
