Discussion:
[petsc-users] petsc-users Digest, Vol 83, Issue 34
Åsmund Ervik
2015-11-10 12:51:30 UTC
Permalink
Subject: [petsc-users] Memory usage with DMDACreate3d and
DMDAGetCorners
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
Loading...