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
