Lines Matching refs:j
112 PetscInt i, j, k, xs, ys, xm, ym, zs, zm; in FormInitialGuess() local
124 for (j = ys; j < ys + ym; j++) { in FormInitialGuess()
125 for (i = xs; i < xs + xm; i++) x[k][j][i] = user->tleft; in FormInitialGuess()
135 PetscInt i, j, k, mx, my, mz, xs, ys, zs, xm, ym, zm; in FormFunction() local
167 for (j = ys; j < ys + ym; j++) { in FormFunction()
169 t0 = x[k][j][i]; in FormFunction()
171 if (i > 0 && i < mx - 1 && j > 0 && j < my - 1 && k > 0 && k < mz - 1) { in FormFunction()
174 tw = x[k][j][i - 1]; in FormFunction()
179 te = x[k][j][i + 1]; in FormFunction()
184 ts = x[k][j - 1][i]; in FormFunction()
189 tn = x[k][j + 1][i]; in FormFunction()
194 td = x[k - 1][j][i]; in FormFunction()
199 tu = x[k + 1][j][i]; in FormFunction()
211 te = x[k][j][i + 1]; in FormFunction()
216 if (j > 0) { in FormFunction()
217 ts = x[k][j - 1][i]; in FormFunction()
225 if (j < my - 1) { in FormFunction()
226 tn = x[k][j + 1][i]; in FormFunction()
235 td = x[k - 1][j][i]; in FormFunction()
244 tu = x[k + 1][j][i]; in FormFunction()
254 tw = x[k][j][i - 1]; in FormFunction()
264 if (j > 0) { in FormFunction()
265 ts = x[k][j - 1][i]; in FormFunction()
273 if (j < my - 1) { in FormFunction()
274 tn = x[k][j + 1][i]; in FormFunction()
283 td = x[k - 1][j][i]; in FormFunction()
292 tu = x[k + 1][j][i]; in FormFunction()
300 } else if (j == 0) { in FormFunction()
302 tw = x[k][j][i - 1]; in FormFunction()
307 te = x[k][j][i + 1]; in FormFunction()
314 tn = x[k][j + 1][i]; in FormFunction()
320 td = x[k - 1][j][i]; in FormFunction()
329 tu = x[k + 1][j][i]; in FormFunction()
337 } else if (j == my - 1) { in FormFunction()
339 tw = x[k][j][i - 1]; in FormFunction()
344 te = x[k][j][i + 1]; in FormFunction()
349 ts = x[k][j - 1][i]; in FormFunction()
357 td = x[k - 1][j][i]; in FormFunction()
366 tu = x[k + 1][j][i]; in FormFunction()
376 tw = x[k][j][i - 1]; in FormFunction()
381 te = x[k][j][i + 1]; in FormFunction()
386 ts = x[k][j - 1][i]; in FormFunction()
391 tn = x[k][j + 1][i]; in FormFunction()
398 tu = x[k + 1][j][i]; in FormFunction()
405 tw = x[k][j][i - 1]; in FormFunction()
410 te = x[k][j][i + 1]; in FormFunction()
415 ts = x[k][j - 1][i]; in FormFunction()
420 tn = x[k][j + 1][i]; in FormFunction()
425 td = x[k - 1][j][i]; in FormFunction()
433 f[k][j][i] = -hyhzdhx * (fe - fw) - hzhxdhy * (fn - fs) - hxhydhz * (fu - fd); in FormFunction()
447 PetscInt i, j, k, mx, my, mz, xs, ys, zs, xm, ym, zm; in FormJacobian() local
481 for (j = ys; j < ys + ym; j++) { in FormJacobian()
483 t0 = x[k][j][i]; in FormJacobian()
485 row.j = j; in FormJacobian()
487 if (i > 0 && i < mx - 1 && j > 0 && j < my - 1 && k > 0 && k < mz - 1) { in FormJacobian()
490 tw = x[k][j][i - 1]; in FormJacobian()
497 te = x[k][j][i + 1]; in FormJacobian()
504 ts = x[k][j - 1][i]; in FormJacobian()
511 tn = x[k][j + 1][i]; in FormJacobian()
518 td = x[k - 1][j][i]; in FormJacobian()
525 tu = x[k + 1][j][i]; in FormJacobian()
533 c[0].j = j; in FormJacobian()
537 c[1].j = j - 1; in FormJacobian()
541 c[2].j = j; in FormJacobian()
545 c[3].j = j; in FormJacobian()
549 c[4].j = j; in FormJacobian()
553 c[5].j = j + 1; in FormJacobian()
557 c[6].j = j; in FormJacobian()
571 te = x[k][j][i + 1]; in FormJacobian()
579 if (j == 0) { in FormJacobian()
580 tn = x[k][j + 1][i]; in FormJacobian()
589 tu = x[k + 1][j][i]; in FormJacobian()
597 c[0].j = j; in FormJacobian()
601 c[1].j = j; in FormJacobian()
605 c[2].j = j + 1; in FormJacobian()
609 c[3].j = j; in FormJacobian()
616 tu = x[k + 1][j][i]; in FormJacobian()
623 td = x[k - 1][j][i]; in FormJacobian()
631 c[0].j = j; in FormJacobian()
635 c[1].j = j; in FormJacobian()
639 c[2].j = j; in FormJacobian()
643 c[3].j = j + 1; in FormJacobian()
647 c[4].j = j; in FormJacobian()
654 td = x[k - 1][j][i]; in FormJacobian()
662 c[0].j = j; in FormJacobian()
666 c[1].j = j; in FormJacobian()
670 c[2].j = j; in FormJacobian()
674 c[3].j = j + 1; in FormJacobian()
681 } else if (j == my - 1) { in FormJacobian()
682 ts = x[k][j - 1][i]; in FormJacobian()
691 tu = x[k + 1][j][i]; in FormJacobian()
699 c[0].j = j - 1; in FormJacobian()
703 c[1].j = j; in FormJacobian()
707 c[2].j = j; in FormJacobian()
711 c[3].j = j; in FormJacobian()
718 tu = x[k + 1][j][i]; in FormJacobian()
725 td = x[k - 1][j][i]; in FormJacobian()
733 c[0].j = j; in FormJacobian()
737 c[1].j = j - 1; in FormJacobian()
741 c[2].j = j; in FormJacobian()
745 c[3].j = j; in FormJacobian()
749 c[4].j = j; in FormJacobian()
756 td = x[k - 1][j][i]; in FormJacobian()
764 c[0].j = j; in FormJacobian()
768 c[1].j = j - 1; in FormJacobian()
772 c[2].j = j; in FormJacobian()
776 c[3].j = j; in FormJacobian()
783 ts = x[k][j - 1][i]; in FormJacobian()
790 tn = x[k][j + 1][i]; in FormJacobian()
799 tu = x[k + 1][j][i]; in FormJacobian()
807 c[0].j = j - 1; in FormJacobian()
811 c[1].j = j; in FormJacobian()
815 c[2].j = j; in FormJacobian()
819 c[3].j = j + 1; in FormJacobian()
823 c[4].j = j; in FormJacobian()
830 td = x[k - 1][j][i]; in FormJacobian()
838 c[0].j = j; in FormJacobian()
842 c[1].j = j - 1; in FormJacobian()
846 c[2].j = j; in FormJacobian()
850 c[3].j = j; in FormJacobian()
854 c[4].j = j + 1; in FormJacobian()
860 td = x[k - 1][j][i]; in FormJacobian()
867 tu = x[k + 1][j][i]; in FormJacobian()
875 c[0].j = j; in FormJacobian()
879 c[1].j = j - 1; in FormJacobian()
883 c[2].j = j; in FormJacobian()
887 c[3].j = j; in FormJacobian()
891 c[4].j = j + 1; in FormJacobian()
895 c[5].j = j; in FormJacobian()
904 tw = x[k][j][i - 1]; in FormJacobian()
919 if (j == 0) { in FormJacobian()
920 tn = x[k][j + 1][i]; in FormJacobian()
929 tu = x[k + 1][j][i]; in FormJacobian()
937 c[0].j = j; in FormJacobian()
941 c[1].j = j; in FormJacobian()
945 c[2].j = j + 1; in FormJacobian()
949 c[3].j = j; in FormJacobian()
956 tu = x[k + 1][j][i]; in FormJacobian()
963 td = x[k - 1][j][i]; in FormJacobian()
971 c[0].j = j; in FormJacobian()
975 c[1].j = j; in FormJacobian()
979 c[2].j = j; in FormJacobian()
983 c[3].j = j + 1; in FormJacobian()
987 c[4].j = j; in FormJacobian()
994 td = x[k - 1][j][i]; in FormJacobian()
1002 c[0].j = j; in FormJacobian()
1006 c[1].j = j; in FormJacobian()
1010 c[2].j = j; in FormJacobian()
1014 c[3].j = j + 1; in FormJacobian()
1021 } else if (j == my - 1) { in FormJacobian()
1022 ts = x[k][j - 1][i]; in FormJacobian()
1031 tu = x[k + 1][j][i]; in FormJacobian()
1039 c[0].j = j - 1; in FormJacobian()
1043 c[1].j = j; in FormJacobian()
1047 c[2].j = j; in FormJacobian()
1051 c[3].j = j; in FormJacobian()
1058 tu = x[k + 1][j][i]; in FormJacobian()
1065 td = x[k - 1][j][i]; in FormJacobian()
1073 c[0].j = j; in FormJacobian()
1077 c[1].j = j - 1; in FormJacobian()
1081 c[2].j = j; in FormJacobian()
1085 c[3].j = j; in FormJacobian()
1089 c[4].j = j; in FormJacobian()
1096 td = x[k - 1][j][i]; in FormJacobian()
1104 c[0].j = j; in FormJacobian()
1108 c[1].j = j - 1; in FormJacobian()
1112 c[2].j = j; in FormJacobian()
1116 c[3].j = j; in FormJacobian()
1123 ts = x[k][j - 1][i]; in FormJacobian()
1130 tn = x[k][j + 1][i]; in FormJacobian()
1139 tu = x[k + 1][j][i]; in FormJacobian()
1147 c[0].j = j - 1; in FormJacobian()
1151 c[1].j = j; in FormJacobian()
1155 c[2].j = j; in FormJacobian()
1159 c[3].j = j + 1; in FormJacobian()
1163 c[4].j = j; in FormJacobian()
1170 td = x[k - 1][j][i]; in FormJacobian()
1178 c[0].j = j; in FormJacobian()
1182 c[1].j = j - 1; in FormJacobian()
1186 c[2].j = j; in FormJacobian()
1190 c[3].j = j; in FormJacobian()
1194 c[4].j = j + 1; in FormJacobian()
1201 td = x[k - 1][j][i]; in FormJacobian()
1208 tu = x[k + 1][j][i]; in FormJacobian()
1216 c[0].j = j; in FormJacobian()
1220 c[1].j = j - 1; in FormJacobian()
1224 c[2].j = j; in FormJacobian()
1228 c[3].j = j; in FormJacobian()
1232 c[4].j = j + 1; in FormJacobian()
1236 c[5].j = j; in FormJacobian()
1243 } else if (j == 0) { in FormJacobian()
1244 tw = x[k][j][i - 1]; in FormJacobian()
1251 te = x[k][j][i + 1]; in FormJacobian()
1258 tn = x[k][j + 1][i]; in FormJacobian()
1267 tu = x[k + 1][j][i]; in FormJacobian()
1275 c[0].j = j; in FormJacobian()
1279 c[1].j = j; in FormJacobian()
1283 c[2].j = j; in FormJacobian()
1287 c[3].j = j + 1; in FormJacobian()
1291 c[4].j = j; in FormJacobian()
1298 td = x[k - 1][j][i]; in FormJacobian()
1306 c[0].j = j; in FormJacobian()
1310 c[1].j = j; in FormJacobian()
1314 c[2].j = j; in FormJacobian()
1318 c[3].j = j; in FormJacobian()
1322 c[4].j = j + 1; in FormJacobian()
1329 tu = x[k + 1][j][i]; in FormJacobian()
1336 td = x[k - 1][j][i]; in FormJacobian()
1344 c[0].j = j; in FormJacobian()
1348 c[1].j = j; in FormJacobian()
1352 c[2].j = j; in FormJacobian()
1356 c[3].j = j; in FormJacobian()
1360 c[4].j = j + 1; in FormJacobian()
1364 c[5].j = j; in FormJacobian()
1370 } else if (j == my - 1) { in FormJacobian()
1371 tw = x[k][j][i - 1]; in FormJacobian()
1378 te = x[k][j][i + 1]; in FormJacobian()
1385 ts = x[k][j - 1][i]; in FormJacobian()
1394 tu = x[k + 1][j][i]; in FormJacobian()
1402 c[0].j = j - 1; in FormJacobian()
1406 c[1].j = j; in FormJacobian()
1410 c[2].j = j; in FormJacobian()
1414 c[3].j = j; in FormJacobian()
1418 c[4].j = j; in FormJacobian()
1425 td = x[k - 1][j][i]; in FormJacobian()
1433 c[0].j = j; in FormJacobian()
1437 c[1].j = j - 1; in FormJacobian()
1441 c[2].j = j; in FormJacobian()
1445 c[3].j = j; in FormJacobian()
1449 c[4].j = j; in FormJacobian()
1456 tu = x[k + 1][j][i]; in FormJacobian()
1463 td = x[k - 1][j][i]; in FormJacobian()
1471 c[0].j = j; in FormJacobian()
1475 c[1].j = j - 1; in FormJacobian()
1479 c[2].j = j; in FormJacobian()
1483 c[3].j = j; in FormJacobian()
1487 c[4].j = j; in FormJacobian()
1491 c[5].j = j; in FormJacobian()
1500 tw = x[k][j][i - 1]; in FormJacobian()
1507 te = x[k][j][i + 1]; in FormJacobian()
1514 ts = x[k][j - 1][i]; in FormJacobian()
1521 tn = x[k][j + 1][i]; in FormJacobian()
1528 tu = x[k + 1][j][i]; in FormJacobian()
1536 c[0].j = j - 1; in FormJacobian()
1540 c[1].j = j; in FormJacobian()
1544 c[2].j = j; in FormJacobian()
1548 c[3].j = j; in FormJacobian()
1552 c[4].j = j + 1; in FormJacobian()
1556 c[5].j = j; in FormJacobian()
1564 tw = x[k][j][i - 1]; in FormJacobian()
1571 te = x[k][j][i + 1]; in FormJacobian()
1578 ts = x[k][j - 1][i]; in FormJacobian()
1585 tn = x[k][j + 1][i]; in FormJacobian()
1592 td = x[k - 1][j][i]; in FormJacobian()
1600 c[0].j = j; in FormJacobian()
1604 c[1].j = j - 1; in FormJacobian()
1608 c[2].j = j; in FormJacobian()
1612 c[3].j = j; in FormJacobian()
1616 c[4].j = j; in FormJacobian()
1620 c[5].j = j + 1; in FormJacobian()