On Thu, Jul 13, 2023 at 12:43 PM Cea, Stephen M <stephen.m....@intel.com> wrote:
> Hi, > > > > I had a questions about PETSc and pardiso integration. Can > you use any of the built into PETSc matrix reordering routines with pardiso > or just the built in options in pardiso? > Currently, you can only use the built-in reorderings. However, it would be simple to extend. We already support using a user-supplied permutation -mat_mkl_pardiso_5 but there is no interface for setting the "perm" array in the Pardiso struct. We could add an interface to take the IS from MatGetOrdering() and stick it into this array. Note that Pardiso only supports symmetric reorderings. Thanks, Matt > Thanks! > > steve > -- What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead. -- Norbert Wiener https://www.cse.buffalo.edu/~knepley/ <http://www.cse.buffalo.edu/~knepley/>