Lines Matching refs:k
112 PetscInt i, j, k, xs, ys, xm, ym, zs, zm; in FormInitialGuess() local
123 for (k = zs; k < zs + zm; k++) { 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
166 for (k = zs; k < zs + zm; k++) { 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()
217 ts = x[k][j - 1][i]; in FormFunction()
226 tn = x[k][j + 1][i]; in FormFunction()
234 if (k > 0) { in FormFunction()
235 td = x[k - 1][j][i]; in FormFunction()
243 if (k < mz - 1) { in FormFunction()
244 tu = x[k + 1][j][i]; in FormFunction()
254 tw = x[k][j][i - 1]; in FormFunction()
265 ts = x[k][j - 1][i]; in FormFunction()
274 tn = x[k][j + 1][i]; in FormFunction()
282 if (k > 0) { in FormFunction()
283 td = x[k - 1][j][i]; in FormFunction()
291 if (k < mz - 1) { in FormFunction()
292 tu = x[k + 1][j][i]; 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()
319 if (k > 0) { in FormFunction()
320 td = x[k - 1][j][i]; in FormFunction()
328 if (k < mz - 1) { in FormFunction()
329 tu = x[k + 1][j][i]; 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()
356 if (k > 0) { in FormFunction()
357 td = x[k - 1][j][i]; in FormFunction()
365 if (k < mz - 1) { in FormFunction()
366 tu = x[k + 1][j][i]; in FormFunction()
374 } else if (k == 0) { 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()
403 } else if (k == mz - 1) { 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
480 for (k = zs; k < zs + zm; k++) { in FormJacobian()
483 t0 = x[k][j][i]; in FormJacobian()
484 row.k = k; 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()
532 c[0].k = k - 1; in FormJacobian()
536 c[1].k = k; in FormJacobian()
540 c[2].k = k; in FormJacobian()
544 c[3].k = k; in FormJacobian()
548 c[4].k = k; in FormJacobian()
552 c[5].k = k; in FormJacobian()
556 c[6].k = k + 1; in FormJacobian()
571 te = x[k][j][i + 1]; in FormJacobian()
580 tn = x[k][j + 1][i]; in FormJacobian()
588 if (k == 0) { in FormJacobian()
589 tu = x[k + 1][j][i]; in FormJacobian()
596 c[0].k = k; in FormJacobian()
600 c[1].k = k; in FormJacobian()
604 c[2].k = k; in FormJacobian()
608 c[3].k = k + 1; in FormJacobian()
615 } else if (k < mz - 1) { in FormJacobian()
616 tu = x[k + 1][j][i]; in FormJacobian()
623 td = x[k - 1][j][i]; in FormJacobian()
630 c[0].k = k - 1; in FormJacobian()
634 c[1].k = k; in FormJacobian()
638 c[2].k = k; in FormJacobian()
642 c[3].k = k; in FormJacobian()
646 c[4].k = k + 1; in FormJacobian()
654 td = x[k - 1][j][i]; in FormJacobian()
661 c[0].k = k - 1; in FormJacobian()
665 c[1].k = k; in FormJacobian()
669 c[2].k = k; in FormJacobian()
673 c[3].k = k; in FormJacobian()
682 ts = x[k][j - 1][i]; in FormJacobian()
690 if (k == 0) { in FormJacobian()
691 tu = x[k + 1][j][i]; in FormJacobian()
698 c[0].k = k; in FormJacobian()
702 c[1].k = k; in FormJacobian()
706 c[2].k = k; in FormJacobian()
710 c[3].k = k + 1; in FormJacobian()
717 } else if (k < mz - 1) { in FormJacobian()
718 tu = x[k + 1][j][i]; in FormJacobian()
725 td = x[k - 1][j][i]; in FormJacobian()
732 c[0].k = k - 1; in FormJacobian()
736 c[1].k = k; in FormJacobian()
740 c[2].k = k; in FormJacobian()
744 c[3].k = k; in FormJacobian()
748 c[4].k = k + 1; in FormJacobian()
756 td = x[k - 1][j][i]; in FormJacobian()
763 c[0].k = k - 1; in FormJacobian()
767 c[1].k = k; in FormJacobian()
771 c[2].k = k; in FormJacobian()
775 c[3].k = k; in FormJacobian()
783 ts = x[k][j - 1][i]; in FormJacobian()
790 tn = x[k][j + 1][i]; in FormJacobian()
798 if (k == 0) { in FormJacobian()
799 tu = x[k + 1][j][i]; in FormJacobian()
806 c[0].k = k; in FormJacobian()
810 c[1].k = k; in FormJacobian()
814 c[2].k = k; in FormJacobian()
818 c[3].k = k; in FormJacobian()
822 c[4].k = k + 1; in FormJacobian()
828 } else if (k == mz - 1) { /* left-hand up interior edge */ in FormJacobian()
830 td = x[k - 1][j][i]; in FormJacobian()
837 c[0].k = k - 1; in FormJacobian()
841 c[1].k = k; in FormJacobian()
845 c[2].k = k; in FormJacobian()
849 c[3].k = k; in FormJacobian()
853 c[4].k = k; in FormJacobian()
860 td = x[k - 1][j][i]; in FormJacobian()
867 tu = x[k + 1][j][i]; in FormJacobian()
874 c[0].k = k - 1; in FormJacobian()
878 c[1].k = k; in FormJacobian()
882 c[2].k = k; in FormJacobian()
886 c[3].k = k; in FormJacobian()
890 c[4].k = k; in FormJacobian()
894 c[5].k = k + 1; in FormJacobian()
904 tw = x[k][j][i - 1]; in FormJacobian()
920 tn = x[k][j + 1][i]; in FormJacobian()
928 if (k == 0) { in FormJacobian()
929 tu = x[k + 1][j][i]; in FormJacobian()
936 c[0].k = k; in FormJacobian()
940 c[1].k = k; in FormJacobian()
944 c[2].k = k; in FormJacobian()
948 c[3].k = k + 1; in FormJacobian()
955 } else if (k < mz - 1) { in FormJacobian()
956 tu = x[k + 1][j][i]; in FormJacobian()
963 td = x[k - 1][j][i]; in FormJacobian()
970 c[0].k = k - 1; in FormJacobian()
974 c[1].k = k; in FormJacobian()
978 c[2].k = k; in FormJacobian()
982 c[3].k = k; in FormJacobian()
986 c[4].k = k + 1; in FormJacobian()
994 td = x[k - 1][j][i]; in FormJacobian()
1001 c[0].k = k - 1; in FormJacobian()
1005 c[1].k = k; in FormJacobian()
1009 c[2].k = k; in FormJacobian()
1013 c[3].k = k; in FormJacobian()
1022 ts = x[k][j - 1][i]; in FormJacobian()
1030 if (k == 0) { in FormJacobian()
1031 tu = x[k + 1][j][i]; in FormJacobian()
1038 c[0].k = k; in FormJacobian()
1042 c[1].k = k; in FormJacobian()
1046 c[2].k = k; in FormJacobian()
1050 c[3].k = k + 1; in FormJacobian()
1057 } else if (k < mz - 1) { in FormJacobian()
1058 tu = x[k + 1][j][i]; in FormJacobian()
1065 td = x[k - 1][j][i]; in FormJacobian()
1072 c[0].k = k - 1; in FormJacobian()
1076 c[1].k = k; in FormJacobian()
1080 c[2].k = k; in FormJacobian()
1084 c[3].k = k; in FormJacobian()
1088 c[4].k = k + 1; in FormJacobian()
1096 td = x[k - 1][j][i]; in FormJacobian()
1103 c[0].k = k - 1; in FormJacobian()
1107 c[1].k = k; in FormJacobian()
1111 c[2].k = k; in FormJacobian()
1115 c[3].k = k; in FormJacobian()
1123 ts = x[k][j - 1][i]; in FormJacobian()
1130 tn = x[k][j + 1][i]; in FormJacobian()
1138 if (k == 0) { in FormJacobian()
1139 tu = x[k + 1][j][i]; in FormJacobian()
1146 c[0].k = k; in FormJacobian()
1150 c[1].k = k; in FormJacobian()
1154 c[2].k = k; in FormJacobian()
1158 c[3].k = k; in FormJacobian()
1162 c[4].k = k + 1; in FormJacobian()
1168 } else if (k == mz - 1) { /* right-hand up interior edge */ in FormJacobian()
1170 td = x[k - 1][j][i]; in FormJacobian()
1177 c[0].k = k - 1; in FormJacobian()
1181 c[1].k = k; in FormJacobian()
1185 c[2].k = k; in FormJacobian()
1189 c[3].k = k; in FormJacobian()
1193 c[4].k = k; in FormJacobian()
1201 td = x[k - 1][j][i]; in FormJacobian()
1208 tu = x[k + 1][j][i]; in FormJacobian()
1215 c[0].k = k - 1; in FormJacobian()
1219 c[1].k = k; in FormJacobian()
1223 c[2].k = k; in FormJacobian()
1227 c[3].k = k; in FormJacobian()
1231 c[4].k = k; in FormJacobian()
1235 c[5].k = k + 1; 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()
1266 if (k == 0) { in FormJacobian()
1267 tu = x[k + 1][j][i]; in FormJacobian()
1274 c[0].k = k; in FormJacobian()
1278 c[1].k = k; in FormJacobian()
1282 c[2].k = k; in FormJacobian()
1286 c[3].k = k; in FormJacobian()
1290 c[4].k = k + 1; in FormJacobian()
1296 } else if (k == mz - 1) { /* bottom up interior edge */ in FormJacobian()
1298 td = x[k - 1][j][i]; in FormJacobian()
1305 c[0].k = k - 1; in FormJacobian()
1309 c[1].k = k; in FormJacobian()
1313 c[2].k = k; in FormJacobian()
1317 c[3].k = k; in FormJacobian()
1321 c[4].k = k; in FormJacobian()
1329 tu = x[k + 1][j][i]; in FormJacobian()
1336 td = x[k - 1][j][i]; in FormJacobian()
1343 c[0].k = k - 1; in FormJacobian()
1347 c[1].k = k; in FormJacobian()
1351 c[2].k = k; in FormJacobian()
1355 c[3].k = k; in FormJacobian()
1359 c[4].k = k; in FormJacobian()
1363 c[5].k = k + 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()
1393 if (k == 0) { in FormJacobian()
1394 tu = x[k + 1][j][i]; in FormJacobian()
1401 c[0].k = k; in FormJacobian()
1405 c[1].k = k; in FormJacobian()
1409 c[2].k = k; in FormJacobian()
1413 c[3].k = k; in FormJacobian()
1417 c[4].k = k + 1; in FormJacobian()
1423 } else if (k == mz - 1) { /* top up interior edge */ in FormJacobian()
1425 td = x[k - 1][j][i]; in FormJacobian()
1432 c[0].k = k - 1; in FormJacobian()
1436 c[1].k = k; in FormJacobian()
1440 c[2].k = k; in FormJacobian()
1444 c[3].k = k; in FormJacobian()
1448 c[4].k = k; in FormJacobian()
1456 tu = x[k + 1][j][i]; in FormJacobian()
1463 td = x[k - 1][j][i]; in FormJacobian()
1470 c[0].k = k - 1; in FormJacobian()
1474 c[1].k = k; in FormJacobian()
1478 c[2].k = k; in FormJacobian()
1482 c[3].k = k; in FormJacobian()
1486 c[4].k = k; in FormJacobian()
1490 c[5].k = k + 1; in FormJacobian()
1497 } else if (k == 0) { 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()
1535 c[0].k = k; in FormJacobian()
1539 c[1].k = k; in FormJacobian()
1543 c[2].k = k; in FormJacobian()
1547 c[3].k = k; in FormJacobian()
1551 c[4].k = k; in FormJacobian()
1555 c[5].k = k + 1; in FormJacobian()
1561 } else if (k == mz - 1) { 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()
1599 c[0].k = k - 1; in FormJacobian()
1603 c[1].k = k; in FormJacobian()
1607 c[2].k = k; in FormJacobian()
1611 c[3].k = k; in FormJacobian()
1615 c[4].k = k; in FormJacobian()
1619 c[5].k = k; in FormJacobian()