Undo changes of 3560f54e3eb3 and change to result a little bit to avoid allocations. Trying different signatres reveals some strange interactions with julias dispatch: allocate: getrange(gridsize::Integer, closuresize::Integer, region::Type{<:Region}) no allocate: getrange(gridsize::Integer, closuresize::Integer, region::R) where R allocate: getrange(gridsize::Integer, closuresize::Integer, region::R) where R <: Type{<:Region} no allocate: getrange(gridsize::Integer, closuresize::Integer, region::DataType)

# D2 order 4

boundary_stencils
     2    -5       4      -1     0     0
     1    -2       1       0     0     0
 -4/43 59/43 -110/43   59/43 -4/43     0
 -1/49     0   59/49 -118/49 64/49 -4/49

inner_stencil
-1/12     4/3  -5/2   4/3 -1/12

e
1

d1
-11/6 3 -3/2 1/3