<div dir="ltr"><div class="gmail_extra">Hong,</div><div class="gmail_extra"><br>&gt; Did you follow same procedural call as petsc examples, e.g., petsc-dev/src/ksp/ksp/examples/tutorials/ex2.c:</div><div class="gmail_extra">

&gt;   ierr = MatCreate(PETSC_COMM_WORLD,&amp;A);CHKERRQ(ierr);</div><div class="gmail_extra">&gt;  ierr = MatSetSizes(A,PETSC_DECIDE,PETSC_DECIDE,m*n,m*n);CHKERRQ(ierr);</div><div class="gmail_extra">&gt;  ierr = MatSetFromOptions(A);CHKERRQ(ierr);</div>

<div class="gmail_extra">&gt;  ierr = MatMPIAIJSetPreallocation(A,5,PETSC_NULL,5,PETSC_NULL);CHKERRQ(ierr);</div><div class="gmail_extra">&gt;  ierr = MatSeqAIJSetPreallocation(A,5,PETSC_NULL);CHKERRQ(ierr);</div><div class="gmail_extra">

<br></div><div class="gmail_extra">Well what I do is *similar*</div><div class="gmail_extra"><br></div>    MatCreate(comm, &amp;A);<br>MatSetSizes(A, localRowSize, localColumnSize, globalRowSize, globalColumnSize);<br>    MatSetType(A, MATMPIAIJ);<br>

    MatMPIAIJSetPreallocation(A, 0, d_nnz, 0, o_nnz);<br>    MatSetFromOptions(A);<br>MatGetOwnershipRange(A, &amp;rStart, &amp;rEnd);<div><br></div><div><br></div><div>This (even without MatSetType(A, MATMPIAIJ);) works with 3.2-p6 but not dev. The only difference I can see is 1) the order of MatSetFromOptions and 2) I do not call MatSeqAIJSetPreallocation which I think I do not need anyway. Is there something I&#39;m doing wrong?</div>

<div><br></div><div>Mohammd</div></div>