Fellow PETSc developers, 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
* If MAT_COPY_VALUES has been specified, then mkl_sparse_optimize() should be executed for the new matrix inside MatDuplicate_SeqAIJMKL(). * 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). 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. Thanks, Richard
