OK Barry, I'll have a look - thanks. Should I use DMShellSetCreateMatrix to set the routine I write for matrix allocation? -------------- next part -------------- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20171031/3d92e276/attachment.html>