#include "petscdmda.h" PetscErrorCode DMDAGetPreallocationCenterDimension(DM dm, PetscInt *preallocCenterDim)
FEM: Two points p and q are adjacent if q \in closure(star(p)), preallocCenterDim = dim
FVM: Two points p and q are adjacent if q \in star(cone(p)), preallocCenterDim = dim-1
FVM++: Two points p and q are adjacent if q \in star(closure(p)), preallocCenterDim = 0