
<html><body><div style="font-family: times new roman, new york, times, serif; font-size: 12pt; color: #000000"><div>If I use MatISSetLocalMat with a preallocated MatCreateSeqAIJ local matrix, do I need to use MatISSetPreallocation for the global matrix ?</div><div><br></div><div>Here is the pseudo-code:<br></div><div>MatCreateIS(PETSC_COMM_WORLD, ..., &globalMat)<br></div><div>MatISSetPreallocation(globalMatrix, ...) // Is this necessary ?<br></div><div>MatCreateSeqAIJ(PETSC_COMM_SELF, ..., &localMatrix) // Prealloc done on the fly<br></div><div>MatSetValues(localMatrix, ...)</div><div>MatISSetLocalMat(globalMatrix, localMatrix)</div><div><br></div><div>Is it necessary to call MatISSetPreallocation for globalMatrix ? (prealloc should have been done locally for each local matrix, no ?)<br></div><div><br></div><div>Franck<br></div></div></body></html>