<div dir="ltr"><div>I think what confused me was the fact that using <span class="gmail-pl-c1">DMCreateMatrix</span>(da,&A) created a 12x12 matrix with 144 elements but summing up nel*nen from each rank gives only 2*2+3*2+3*2+3*2=20 elements. So this means that DMDAGetElements returns the elements for the vector created on the mesh and this happens to be used for indexing the matrix (created using the same DM object) via the local indices. Luckily since this code just creates a tridiagonal matrix, this works here (since if I wanted to create a dense matrix I'd want to have access to all the indices which would happen if the addressable indices at each rank (for the submatrix stored locally) add up to the size of the total matrix). Is my understanding correct? <br><br></div>Thanks for the help! <br><div><div><br>-- <br><div dir="ltr" class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div style="font-size:12.8px">Sajid Ali<br></div><div style="font-size:12.8px">Applied Physics<br></div><div style="font-size:12.8px">Northwestern University</div></div></div></div></div></div>