Lines Matching refs:j
111 PetscInt i, j, xs, ys, xm, ym; in FormInitialGuess() local
125 for (j = ys; j < ys + ym; j++) { in FormInitialGuess()
126 for (i = xs; i < xs + xm; i++) x[j][i] = tleft; in FormInitialGuess()
135 PetscInt i, j, mx, my, xs, ys, xm, ym; in FormFunction() local
164 for (j = ys; j < ys + ym; j++) { in FormFunction()
166 t0 = x[j][i]; in FormFunction()
168 if (i > 0 && i < mx - 1 && j > 0 && j < my - 1) { in FormFunction()
171 tw = x[j][i - 1]; in FormFunction()
176 te = x[j][i + 1]; in FormFunction()
181 ts = x[j - 1][i]; in FormFunction()
186 tn = x[j + 1][i]; in FormFunction()
198 te = x[j][i + 1]; in FormFunction()
203 if (j > 0) { in FormFunction()
204 ts = x[j - 1][i]; in FormFunction()
210 if (j < my - 1) { in FormFunction()
211 tn = x[j + 1][i]; in FormFunction()
219 tw = x[j][i - 1]; in FormFunction()
229 if (j > 0) { in FormFunction()
230 ts = x[j - 1][i]; in FormFunction()
236 if (j < my - 1) { in FormFunction()
237 tn = x[j + 1][i]; in FormFunction()
243 } else if (j == 0) { in FormFunction()
245 tw = x[j][i - 1]; in FormFunction()
250 te = x[j][i + 1]; in FormFunction()
257 tn = x[j + 1][i]; in FormFunction()
262 } else if (j == my - 1) { in FormFunction()
264 tw = x[j][i - 1]; in FormFunction()
269 te = x[j][i + 1]; in FormFunction()
274 ts = x[j - 1][i]; in FormFunction()
282 f[j][i] = -hydhx * (fe - fw) - hxdhy * (fn - fs); in FormFunction()
295 PetscInt i, j, mx, my, xs, ys, xm, ym; in FormJacobian() local
325 for (j = ys; j < ys + ym; j++) { in FormJacobian()
327 t0 = x[j][i]; in FormJacobian()
329 if (i > 0 && i < mx - 1 && j > 0 && j < my - 1) { in FormJacobian()
332 tw = x[j][i - 1]; in FormJacobian()
339 te = x[j][i + 1]; in FormJacobian()
346 ts = x[j - 1][i]; in FormJacobian()
353 tn = x[j + 1][i]; in FormJacobian()
361 col[0].j = j - 1; in FormJacobian()
364 col[1].j = j; in FormJacobian()
367 col[2].j = row.j = j; in FormJacobian()
370 col[3].j = j; in FormJacobian()
373 col[4].j = j + 1; in FormJacobian()
386 te = x[j][i + 1]; in FormJacobian()
394 if (j == 0) { in FormJacobian()
395 tn = x[j + 1][i]; in FormJacobian()
403 col[0].j = row.j = j; in FormJacobian()
406 col[1].j = j; in FormJacobian()
409 col[2].j = j + 1; in FormJacobian()
414 } else if (j < my - 1) { in FormJacobian()
415 ts = x[j - 1][i]; in FormJacobian()
422 tn = x[j + 1][i]; in FormJacobian()
430 col[0].j = j - 1; in FormJacobian()
433 col[1].j = row.j = j; in FormJacobian()
436 col[2].j = j; in FormJacobian()
439 col[3].j = j + 1; in FormJacobian()
444 ts = x[j - 1][i]; in FormJacobian()
452 col[0].j = j - 1; in FormJacobian()
455 col[1].j = row.j = j; in FormJacobian()
458 col[2].j = j; in FormJacobian()
465 tw = x[j][i - 1]; in FormJacobian()
480 if (j == 0) { in FormJacobian()
481 tn = x[j + 1][i]; in FormJacobian()
489 col[0].j = j; in FormJacobian()
492 col[1].j = row.j = j; in FormJacobian()
495 col[2].j = j + 1; in FormJacobian()
500 } else if (j < my - 1) { in FormJacobian()
501 ts = x[j - 1][i]; in FormJacobian()
508 tn = x[j + 1][i]; in FormJacobian()
516 col[0].j = j - 1; in FormJacobian()
519 col[1].j = j; in FormJacobian()
522 col[2].j = row.j = j; in FormJacobian()
525 col[3].j = j + 1; in FormJacobian()
530 ts = x[j - 1][i]; in FormJacobian()
538 col[0].j = j - 1; in FormJacobian()
541 col[1].j = j; in FormJacobian()
544 col[2].j = row.j = j; in FormJacobian()
550 } else if (j == 0) { in FormJacobian()
551 tw = x[j][i - 1]; in FormJacobian()
558 te = x[j][i + 1]; in FormJacobian()
565 tn = x[j + 1][i]; in FormJacobian()
573 col[0].j = j; in FormJacobian()
576 col[1].j = row.j = j; in FormJacobian()
579 col[2].j = j; in FormJacobian()
582 col[3].j = j + 1; in FormJacobian()
587 } else if (j == my - 1) { in FormJacobian()
588 tw = x[j][i - 1]; in FormJacobian()
595 te = x[j][i + 1]; in FormJacobian()
602 ts = x[j - 1][i]; in FormJacobian()
610 col[0].j = j - 1; in FormJacobian()
613 col[1].j = j; in FormJacobian()
616 col[2].j = row.j = j; in FormJacobian()
619 col[3].j = j; in FormJacobian()