How well is PETSC suited for solving elliptic PDEs on four-dimensional domains? Looking at the list of commands, I see e.g. DMDACreate{1d,2d,3d} with no 4d version. On the other hand, I see you have VecGetArray{1d,2d,3d,4d}.
Is there e.g. a convenient way of using the PETSCs data structures for solving the Poisson equation on a 4d "cube"? Best regards Torquil Sørensen