<div dir="ltr"><div><div><div><div><div><div>Fellow PETSc developers,<br><br></div>I'm trying to figure out what my MatDuplicate() implementation for AIJMKL should do when the sparse inspector-executor model is being used. I think that the behavior should be<br><br></div>* If MAT_COPY_VALUES has been specified, then mkl_sparse_optimize() should be executed for the new matrix inside MatDuplicate_SeqAIJMKL().<br></div>* If MAT_COPY_VALUES has not been specified, then don't do the mkl_sparse_optimize(), as the user is going to have to specify the nonzero entries at some point, and then do MatAssemblyBegin/End (which calls mkl_sparse_optimize).<br><br></div>I'm a little unclear on all of the situations in which MatDuplicate gets used, however, so let me know if I'm misunderstanding something.<br><br></div>Thanks,<br></div>Richard<br></div>