24   PetscInt       i,j,k,Mx,My,Mz,xs,ys,zs,xm,ym,zm;  in FormInitGuess()  local
29 …PetscCall(DMDAGetInfo(da,PETSC_IGNORE,&Mx,&My,&Mz,PETSC_IGNORE,PETSC_IGNORE,PETSC_IGNORE,PETSC_IGN…  in FormInitGuess()
34   hz     = 1.0/(PetscReal)(Mz-1);  in FormInitGuess()
62     tempk = (PetscReal)(PetscMin(k,Mz-k-1))*hz;  in FormInitGuess()
66         if (i == 0 || j == 0 || k == 0 || i == Mx-1 || j == My-1 || k == Mz-1) {  in FormInitGuess()
85   PetscInt       i,j,k,Mx,My,Mz,xs,ys,zs,xm,ym,zm;  in FormFunction()  local
91 …PetscCall(DMDAGetInfo(da,PETSC_IGNORE,&Mx,&My,&Mz,PETSC_IGNORE,PETSC_IGNORE,PETSC_IGNORE,PETSC_IGN…  in FormFunction()
95   hz     = 1.0/(PetscReal)(Mz-1);  in FormFunction()
132         if (i == 0 || j == 0 || k == 0 || i == Mx-1 || j == My-1 || k == Mz-1) {  in FormFunction()
165   PetscInt       i,j,k,Mx,My,Mz,xs,ys,zs,xm,ym,zm;  in FormJacobian()  local
172 …PetscCall(DMDAGetInfo(da,PETSC_IGNORE,&Mx,&My,&Mz,PETSC_IGNORE,PETSC_IGNORE,PETSC_IGNORE,PETSC_IGN…  in FormJacobian()
176   hz     = 1.0/(PetscReal)(Mz-1);  in FormJacobian()
226         if (i == 0 || j == 0 || k == 0|| i == Mx-1 || j == My-1 || k == Mz-1) {  in FormJacobian()