| aab60f1b | 27-Oct-2017 |
Richard Tran Mills <rmills@rmills.org> |
MatSeqAIJMKL_create_from_mkl_handle() now simply destroys the old matrix if MAT_REUSE_MATRIX is specified.
Previously, we attempted to save some expense by directly copying exported CSR arrays from
MatSeqAIJMKL_create_from_mkl_handle() now simply destroys the old matrix if MAT_REUSE_MATRIX is specified.
Previously, we attempted to save some expense by directly copying exported CSR arrays from MKL, but this causes problems because the MKL and PETSc representations differ in small ways (e.g., more explicit nonzeros in the PETSc representation due to preallocation.)
show more ...
|