Lines Matching refs:k
24 PetscInt i,j,k,Mx,My,Mz,xs,ys,zs,xm,ym,zm; in FormInitGuess() local
61 for (k=zs; k<zs+zm; k++) { 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()
68 x[k][j][i] = 0.0; 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
129 for (k=zs; k<zs+zm; k++) { in FormFunction()
132 if (i == 0 || j == 0 || k == 0 || i == Mx-1 || j == My-1 || k == Mz-1) { in FormFunction()
134 f[k][j][i] = x[k][j][i] - 0.0; in FormFunction()
137 u = x[k][j][i]; in FormFunction()
138 u_east = x[k][j][i+1]; in FormFunction()
139 u_west = x[k][j][i-1]; in FormFunction()
140 u_north = x[k][j+1][i]; in FormFunction()
141 u_south = x[k][j-1][i]; in FormFunction()
142 u_up = x[k+1][j][i]; in FormFunction()
143 u_down = x[k-1][j][i]; in FormFunction()
147 f[k][j][i] = u_xx + u_yy + u_zz - sc*PetscExpScalar(u); in FormFunction()
165 PetscInt i,j,k,Mx,My,Mz,xs,ys,zs,xm,ym,zm; in FormJacobian() local
221 for (k=zs; k<zs+zm; k++) { in FormJacobian()
224 row.k = k; row.j = j; row.i = i; in FormJacobian()
226 if (i == 0 || j == 0 || k == 0|| i == Mx-1 || j == My-1 || k == Mz-1) { in FormJacobian()
231 v[0] = -hxhydhz; col[0].k=k-1;col[0].j=j; col[0].i = i; in FormJacobian()
232 v[1] = -hxhzdhy; col[1].k=k; col[1].j=j-1;col[1].i = i; in FormJacobian()
233 v[2] = -hyhzdhx; col[2].k=k; col[2].j=j; col[2].i = i-1; in FormJacobian()
234 …v[3] = 2.0*(hyhzdhx+hxhzdhy+hxhydhz)-sc*PetscExpScalar(x[k][j][i]);col[3].k=row.k;col[3].j=row.j;c… in FormJacobian()
235 v[4] = -hyhzdhx; col[4].k=k; col[4].j=j; col[4].i = i+1; in FormJacobian()
236 v[5] = -hxhzdhy; col[5].k=k; col[5].j=j+1;col[5].i = i; in FormJacobian()
237 v[6] = -hxhydhz; col[6].k=k+1;col[6].j=j; col[6].i = i; in FormJacobian()