<div dir="ltr"><div dir="ltr">
<div>Hi Barry, <br></div><div><br></div><div>Thanks for creating the new function. I'm somewhat confused as to how I'd use it. Given an MPIAIJ matrix, is one supposed to extract the local SeqAIJ matrix and set the preallocation on each mpi-rank independently followed by MatSetValues (on the MPIAIJ matrix) to fill the rows? Or, does one create SeqAIJ matrices on each rank and then combine them into a parallel MPIAIJ matrix say by using MatCreateMPIMatConcatenateSeqMat?<br><br></div><div>I tried the second approach but leaving the "number of local columns" for 
MatCreateMPIMatConcatenateSeqMat as PETSC_DECIDE causes a crash (when running with 1 mpi rank). Is this the correct approach to take and if yes what does "number of local columns" mean when combining the seqaij matrices ?<br></div><div><br></div><div><div>Thank You,<br></div><div><div dir="ltr"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div style="font-size:12.8px">Sajid Ali | PhD Candidate<br></div><div style="font-size:12.8px">Applied Physics<br></div><div style="font-size:12.8px">Northwestern University</div><div style="font-size:12.8px"><a href="http://s-sajid-ali.github.io" target="_blank">s-sajid-ali.github.io</a></div></div></div></div></div></div></div></div></div>

</div><br></div>