<div dir="ltr">Sorry, I was getting "composite" from DMCompositeGetAccessArray, I am talking about MatNest.  I can get the matrices with <span style="color:rgb(0,0,0)">MatNestSetSubMats....</span><div class="gmail_extra"><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex"><span class=""><br>
>   ! copy in linear operator<br>
>   call MatCopy(a_ts%FJacobean,Jpre,ierr);CHKERRQ(ierr)<br>
>   ! shift 1 & 2 by 'shift'<br>
>   call MatShift(mat00,shift,ierr);CHKERRQ(ierr)  ????<br>
>   call MatShift(mat11,shift,ierr);CHKERRQ(ierr)  ????<br>
<br>
</span>What does this mean?  Why not include the shift while assembling?<br></blockquote><div><br></div><div>Perhaps I am using this incorrectly but the shift is added in "FormIJacobian" and only 2 of the 3 fields have u_t.  The rest of the matrix is linear and not time dependant so I create it once in the setup and copy it in in FormIJacobian.</div><div><br></div><div>Thanks,</div><div>Mark</div></div><br></div></div>