> From: TAY wee-beng <[email protected]> > To: PETSc users list <[email protected]> > Subject: [petsc-users] Memory usage with DMDACreate3d and > DMDAGetCorners > Message-ID: <[email protected]> > Content-Type: text/plain; charset=utf-8; format=flowed > > Hi, > > I need a subroutine in Fortran to partition a subset of my grid in the 3 > x,y,z directions for MPI. I thought of using DMDACreate3d and > DMDAGetCorners to get the starting and width of the partitioned grid.
You might want to look at dm/examples/ex13f90.F90 (and ex13f90aux.F90) which does a detailed walkthrough of the DMDA routines for a very simple case in 3D. http://www.mcs.anl.gov/petsc/petsc-current/src/dm/examples/tutorials/ex13f90.F90.html > > Because I need to partition at every time step and the subset grid > changes dimension and index at every time step, so I will also need to > use DMDestroy after each time step > > Will that use alot of memory? Will the grid actually be created? So I > wonder if this DMDACreate3d and DMDestroy calls will take a lot of time. > > > -- > Thank you. > > Yours sincerely, > > TAY wee-beng
signature.asc
Description: OpenPGP digital signature
