[petsc-users] High-dimensional DMDA

[Jed Brown]

Matthew Knepley <knepley at gmail.com> writes:

> On Tue, Oct 17, 2017 at 5:01 PM, Jed Brown <jed at jedbrown.org> wrote:
>
>> Tobin Isaac <tisaac at cc.gatech.edu> 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?