I have further checked this function. <br>In MatMatMult(Mat A,Mat B,MatReuse scall,PetscReal fill,Mat *C) <br>I am wondering why&nbsp;the&nbsp;type&nbsp;of&nbsp;C&nbsp;is&nbsp;MATAIJ&nbsp;when&nbsp;the&nbsp;types&nbsp;of&nbsp;A&nbsp;and&nbsp;B&nbsp;are&nbsp;MATAIJ.&nbsp;although&nbsp;A&nbsp;and&nbsp;B&nbsp;is&nbsp;MATAIJ,&nbsp;C&nbsp;should&nbsp;be&nbsp;dense.&nbsp;If&nbsp;C&nbsp;uses&nbsp;MATAIJ&nbsp;type,&nbsp;it&nbsp;should&nbsp;take&nbsp;more&nbsp;memory,&nbsp;is&nbsp;it&nbsp;right?<br>
<br>thanks&nbsp;a&nbsp;lot.<br><br>Regards,<br>Yujie&nbsp;<br>