Lines Matching refs:k
24 PetscInt i,j,k,Mx,My,Mz,xs,ys,zs,xm,ym,zm;
60 for (k=zs; k<zs+zm; k++) {
61 tempk = (PetscReal)(PetscMin(k,Mz-k-1))*hz;
65 if (i == 0 || j == 0 || k == 0 || i == Mx-1 || j == My-1 || k == Mz-1) {
67 x[k][j][i] = 0.0;
69 x[k][j][i] = temp1*sqrt(PetscMin((PetscReal)(PetscMin(i,Mx-i-1))*hx,tempj));
84 PetscInt i,j,k,Mx,My,Mz,xs,ys,zs,xm,ym,zm;
125 for (k=zs; k<zs+zm; k++) {
128 if (i == 0 || j == 0 || k == 0 || i == Mx-1 || j == My-1 || k == Mz-1) {
130 f[k][j][i] = x[k][j][i] - 0.0;
133 u = x[k][j][i];
134 u_east = x[k][j][i+1];
135 u_west = x[k][j][i-1];
136 u_north = x[k][j+1][i];
137 u_south = x[k][j-1][i];
138 u_up = x[k+1][j][i];
139 u_down = x[k-1][j][i];
143 f[k][j][i] = u_xx + u_yy + u_zz - sc*PetscExpScalar(u);
161 PetscInt i,j,k,Mx,My,Mz,xs,ys,zs,xm,ym,zm;
214 for (k=zs; k<zs+zm; k++) {
217 row.k = k; row.j = j; row.i = i;
219 if (i == 0 || j == 0 || k == 0|| i == Mx-1 || j == My-1 || k == Mz-1) {
224 v[0] = -hxhydhz; col[0].k=k-1;col[0].j=j; col[0].i = i;
225 v[1] = -hxhzdhy; col[1].k=k; col[1].j=j-1;col[1].i = i;
226 v[2] = -hyhzdhx; col[2].k=k; col[2].j=j; col[2].i = i-1;
227 v[3] = 2.0*(hyhzdhx+hxhzdhy+hxhydhz)-sc*PetscExpScalar(x[k][j][i]);col[3].k=row.k;col[3].j=row.j;col[3].i = row.i;
228 v[4] = -hyhzdhx; col[4].k=k; col[4].j=j; col[4].i = i+1;
229 v[5] = -hxhzdhy; col[5].k=k; col[5].j=j+1;col[5].i = i;
230 v[6] = -hxhydhz; col[6].k=k+1;col[6].j=j; col[6].i = i;