<br><br><div class="gmail_quote">Thomas :<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">In my FETI-DP code, each rank creates a SEQAIJ matrix that represents the discretization of the interior domain. Just for debugging, I would like to join these sequential matrices to one global MPIAIJ matrix. This matrix has no off diagonal nnzs and should be stored corresponding to the ranks unknowns, thus, first all rows of the first rank and so on. What's the most efficient way to do this? Is it possible to create this parallel matrix just as a view of the sequential ones, so without copying the data? Thanks for any advise.</blockquote>
<div><br></div><div><a href="http://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/Mat/MatCreateMPIAIJConcatenateSeqAIJ.html">http://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/Mat/MatCreateMPIAIJConcatenateSeqAIJ.html</a></div>
<div><br></div><div>Note: entries in seqaij matrices are copied into a mpiaij matrix without</div><div>inter-processor communication. Use petsc-3.3 for this function.</div><div><br></div><div>Hong</div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<span class="HOEnZb"><font color="#888888"><br>
<br>
Thomas<br>
</font></span></blockquote></div><br>