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()
32   hx     = 1.0/(PetscReal)(Mx-1);  in FormInitGuess()
66         if (i == 0 || j == 0 || k == 0 || i == Mx-1 || j == My-1 || k == Mz-1) {  in FormInitGuess()
70           x[k][j][i] = temp1*sqrt(PetscMin((PetscReal)(PetscMin(i,Mx-i-1))*hx,tempj));  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()
93   hx     = 1.0/(PetscReal)(Mx-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()
174   hx     = 1.0/(PetscReal)(Mx-1);  in FormJacobian()
226         if (i == 0 || j == 0 || k == 0|| i == Mx-1 || j == My-1 || k == Mz-1) {  in FormJacobian()