Lines Matching refs:j
159 PetscInt i, j, Mx, My, xs, ys, xm, ym; in FormIFunction() local
192 for (j = ys; j < ys + ym; j++) { in FormIFunction()
195 if (i == 0 || j == 0 || i == Mx - 1 || j == My - 1) { in FormIFunction()
197 f[j][i] = uarray[j][i]; /* F = U */ in FormIFunction()
199 if (i == 0 && j == 0) { /* SW corner */ in FormIFunction()
200 f[j][i] = uarray[j][i] - uarray[j + 1][i + 1]; in FormIFunction()
201 } else if (i == Mx - 1 && j == 0) { /* SE corner */ in FormIFunction()
202 f[j][i] = uarray[j][i] - uarray[j + 1][i - 1]; in FormIFunction()
203 } else if (i == 0 && j == My - 1) { /* NW corner */ in FormIFunction()
204 f[j][i] = uarray[j][i] - uarray[j - 1][i + 1]; in FormIFunction()
205 } else if (i == Mx - 1 && j == My - 1) { /* NE corner */ in FormIFunction()
206 f[j][i] = uarray[j][i] - uarray[j - 1][i - 1]; in FormIFunction()
208 f[j][i] = uarray[j][i] - uarray[j][i + 1]; in FormIFunction()
210 f[j][i] = uarray[j][i] - uarray[j][i - 1]; in FormIFunction()
211 } else if (j == 0) { /* Bottom */ in FormIFunction()
212 f[j][i] = uarray[j][i] - uarray[j + 1][i]; in FormIFunction()
213 } else if (j == My - 1) { /* Top */ in FormIFunction()
214 f[j][i] = uarray[j][i] - uarray[j - 1][i]; in FormIFunction()
218 u = uarray[j][i]; in FormIFunction()
220 uxx = (-2.0 * u + uarray[j][i - 1] + uarray[j][i + 1]); in FormIFunction()
221 uyy = (-2.0 * u + uarray[j - 1][i] + uarray[j + 1][i]); in FormIFunction()
224 …uxx = 2.0 * uxx / 3.0 + (0.5 * (uarray[j - 1][i - 1] + uarray[j - 1][i + 1] + uarray[j + 1][i - 1]… in FormIFunction()
225 …uyy = 2.0 * uyy / 3.0 + (0.5 * (uarray[j - 1][i - 1] + uarray[j - 1][i + 1] + uarray[j + 1][i - 1]… in FormIFunction()
227 f[j][i] = udot[j][i] - (uxx * sx + uyy * sy); in FormIFunction()
248 PetscInt i, j, Mx, My, xs, ys, xm, ym, nc; in FormIJacobian() local
263 for (j = ys; j < ys + ym; j++) { in FormIJacobian()
266 row.j = j; in FormIJacobian()
268 if (user->boundary == 0 && (i == 0 || i == Mx - 1 || j == 0 || j == My - 1)) { in FormIJacobian()
269 col[nc].j = j; in FormIJacobian()
274 col[nc].j = j; in FormIJacobian()
277 col[nc].j = j; in FormIJacobian()
281 col[nc].j = j; in FormIJacobian()
284 col[nc].j = j; in FormIJacobian()
287 } else if (user->boundary > 0 && j == 0) { /* Bottom Neumann */ in FormIJacobian()
288 col[nc].j = j; in FormIJacobian()
291 col[nc].j = j + 1; in FormIJacobian()
294 } else if (user->boundary > 0 && j == My - 1) { /* Top Neumann */ in FormIJacobian()
295 col[nc].j = j; in FormIJacobian()
298 col[nc].j = j - 1; in FormIJacobian()
302 col[nc].j = j - 1; in FormIJacobian()
305 col[nc].j = j; in FormIJacobian()
308 col[nc].j = j; in FormIJacobian()
311 col[nc].j = j; in FormIJacobian()
314 col[nc].j = j + 1; in FormIJacobian()
341 PetscInt i, j, xs, ys, xm, ym, Mx, My; in FormInitialSolution() local
358 for (j = ys; j < ys + ym; j++) { in FormInitialSolution()
359 y = j * hy; in FormInitialSolution()
363 if (r < .125) u[j][i] = PetscExpReal(c * r * r * r); in FormInitialSolution()
364 else u[j][i] = 0.0; in FormInitialSolution()