[petsc-dev] What should MatDuplicate_SeqAIJMKL() do?

Richard Tran Mills rtmills at anl.gov
Wed Jul 19 10:55:35 CDT 2017

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

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.

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20170719/e6c565c2/attachment.html>

More information about the petsc-dev mailing list