[petsc-users] Automatic partition of Cartesian grids

[TAY wee-beng]

Hi,

I am currently using PETSc for my Fortran CFD code and I am manually 
partitioning my Cartesian grids. So in 2D, it will be something like 
u(1:size_x,jstart:jend), where the y component is partitioned into 2,4,8 
etc parts depending on the no. of processors.

However, it seems that there are better ways to do it to get a more 
balanced load. Do I use DMDA in PETSc to do it? I'm now using staggered 
Cartesian grids for my u,v,p. Is there an example to construct a 
Laplace/Poisson equation using DMDA?

It is mentioned in the manual that DMMG infrastructure will be replaced 
in the next release and we should not use it. Is this related to DMDA?

Yours sincerely,

TAY wee-beng