Re: [petsc-users] DMDACoor2d and DMDACoor3d in Fortran

[Danyang Su]

Hi Jed,

Thanks, the problem is solved.
        call DMGetCoordinateDM(dmda_flow%da,cda,ierr)
        call DMGetCoordinates(dmda_flow%da,gc,ierr)

        call DMDAGetLocalBoundingBox(cda,lmin,lmax,ierr)             
    !this will return node index value in x,y,z dim
        call DMDAGetLocalBoundingBox(dmda_flow%da,lmin,lmax,ierr)   
!this will return coordinate value

Thanks and regards,

Danyang


On 27/04/2014 7:00 AM, Jed Brown wrote:
Danyang Su <danyang...@gmail.com> writes:

Hi Barry,

Another question is about DMDAGetLocalBoundingBox in fortran.

PetscErrorCode  DMDAGetLocalBoundingBox(DM da,PetscReal lmin[],PetscReal
lmax[])

The first value of lmin and lmax (lmin(0)) are always zero,
lmin and lmax are normal arrays of length 3.  In Fortran, you don't
access lmin(0).

and lmin(1), lmin(2), and lmin(3) are for x, y, and z dimension,
respectively.  And the returned value is index (local node index in x,
y, z dim), not the coordinate. Correct?
It is the coordinate value, assuming you have set DMSetCoordinates.