<div dir="ltr">







<p class="">I am getting a new-non-zero error in PCComputeExplicitOperator at the bottom of precon.c.  it looks like the matrix is set to have not memory:</p><p class=""><br></p><p class="">PetscErrorCode  <span class="">PCComputeExplicitOperator</span>(PC pc,Mat *mat)</p>

<p class=""><br></p><p class=""> ....</p><p class=""><br></p><p class="">  if (size == <span class="">1</span>) {    ierr = MatSetType(*mat,MATSEQDENSE);CHKERRQ(ierr);</p>

<p class="">    ierr = MatSeqDenseSetPreallocation(*mat,<span class="">NULL</span>);CHKERRQ(ierr);</p>
<p class="">  } else {</p>
<p class="">    ierr = MatSetType(*mat,MATMPIAIJ);CHKERRQ(ierr);</p>
<p class="">    ierr = MatMPIAIJSetPreallocation(*mat,<span class="">0</span>,<span class="">NULL</span>,<span class="">0</span>,<span class="">NULL</span>);CHKERRQ(ierr);</p>
<p class="">  }</p>
<p class=""> .....</p>
<p class="">    ierr = MatSetValues(*mat,m,rows,<span class="">1</span>,&i,array,INSERT_VALUES);CHKERRQ(ierr);</p>
<p class=""><br></p><p class="">I get the error on this last line (#1835 precon.c).  Matt mentioned that this code needed work.  I can work on this but am not sure about the best way to proceed.  I am doing a full Schur complement and don't mind simply going through and solving for each column of A_01 in A_00^-1 * A_01.  I guess the pc->pmat is a matrix shell with (A_11 - A_10 * A_00^-1 * A_01).</p>
<p class="">should this just be a dense matrix in parallel as well as not?</p><p class="">Mark</p><p class=""><br></p></div>