On Nov 7, 2013, at 4:31 PM, Jed Brown <[email protected]> wrote: > Barry Smith <[email protected]> writes: >> There is no easier way. Just live with the memory usage. > > DMDA with stencil width of 0 could use contiguous index sets (ISStride) > so that the scatters take up no memory, right? Why should stencil width > 0 use any problem-sized memory?
Maybe it doesn’t. Just give it a try.
