Matthew Knepley <[email protected]> writes: > On Tue, Oct 17, 2017 at 5:01 PM, Jed Brown <[email protected]> wrote: > >> Tobin Isaac <[email protected]> writes: >> >> > Which is something so rote that creating DTFD to convert the pointwise >> > residuals/jacobians into stencils would go along way towards unifying >> > DMDA and DMPlex in a single useful interface. >> >> On a related topic, I have reason to add Fornberg's algorithm for >> computing finite difference weights to DMDT. (It's more stable than the >> usual Vandermonde approach.) I'd imagine some users would appreciate >> that to take some of the mystery/effort out of the design of higher >> order FD methods. >> > > That would be good. > > Is there a reason not to do all DMDA computation with p4est? We could > prescribe the refinement level and have the forest make the > decomposition.
GPL? Are you setting up indexing like DMDA so the same code will work for boundary conditions?
