“overallocating” is exactly what we can live on at the moment, as long as it is easier to work with coloring on dmstag.
So it sounds like if we can provide a preallocated matrix with a proper stencil through DMCreateMatrix, then it should work with dmstag and coloring already. Most APIs are already there. Qi On Dec 13, 2021, at 12:13 PM, Matthew Knepley <knep...@gmail.com<mailto:knep...@gmail.com>> wrote: Yes, and would not handle higher order stencils.I think the overallocating is livable for the first imeplementation.