Lines Matching refs:j
187 PetscInt i, j, mx, xs, ys, xm, ym; in FormInitialGuess() local
217 for (j = ys; j < ys + ym; j++) { in FormInitialGuess()
219 x[j][i].u = 0.0; in FormInitialGuess()
220 x[j][i].v = 0.0; in FormInitialGuess()
221 x[j][i].omega = 0.0; in FormInitialGuess()
222 x[j][i].temp = (grashof > 0) * i * dx; in FormInitialGuess()
236 PetscInt xints, xinte, yints, yinte, i, j; in FormFunctionLocal() local
267 j = 0; in FormFunctionLocal()
271 f[j][i].u = x[j][i].u; in FormFunctionLocal()
272 f[j][i].v = x[j][i].v; in FormFunctionLocal()
273 f[j][i].omega = x[j][i].omega + (x[j + 1][i].u - x[j][i].u) * dhy; in FormFunctionLocal()
274 f[j][i].temp = x[j][i].temp - x[j + 1][i].temp; in FormFunctionLocal()
280 j = info->my - 1; in FormFunctionLocal()
284 f[j][i].u = x[j][i].u - lid; in FormFunctionLocal()
285 f[j][i].v = x[j][i].v; in FormFunctionLocal()
286 f[j][i].omega = x[j][i].omega + (x[j][i].u - x[j - 1][i].u) * dhy; in FormFunctionLocal()
287 f[j][i].temp = x[j][i].temp - x[j - 1][i].temp; in FormFunctionLocal()
296 for (j = info->ys; j < info->ys + info->ym; j++) { in FormFunctionLocal()
297 f[j][i].u = x[j][i].u; in FormFunctionLocal()
298 f[j][i].v = x[j][i].v; in FormFunctionLocal()
299 f[j][i].omega = x[j][i].omega - (x[j][i + 1].v - x[j][i].v) * dhx; in FormFunctionLocal()
300 f[j][i].temp = x[j][i].temp; in FormFunctionLocal()
309 for (j = info->ys; j < info->ys + info->ym; j++) { in FormFunctionLocal()
310 f[j][i].u = x[j][i].u; in FormFunctionLocal()
311 f[j][i].v = x[j][i].v; in FormFunctionLocal()
312 f[j][i].omega = x[j][i].omega - (x[j][i].v - x[j][i - 1].v) * dhx; in FormFunctionLocal()
313 f[j][i].temp = x[j][i].temp - (PetscReal)(grashof > 0); in FormFunctionLocal()
318 for (j = yints; j < yinte; j++) { in FormFunctionLocal()
323 vx = x[j][i].u; in FormFunctionLocal()
327 vy = x[j][i].v; in FormFunctionLocal()
333 u = x[j][i].u; in FormFunctionLocal()
334 uxx = (2.0 * u - x[j][i - 1].u - x[j][i + 1].u) * hydhx; in FormFunctionLocal()
335 uyy = (2.0 * u - x[j - 1][i].u - x[j + 1][i].u) * hxdhy; in FormFunctionLocal()
336 f[j][i].u = uxx + uyy - .5 * (x[j + 1][i].omega - x[j - 1][i].omega) * hx; in FormFunctionLocal()
339 u = x[j][i].v; in FormFunctionLocal()
340 uxx = (2.0 * u - x[j][i - 1].v - x[j][i + 1].v) * hydhx; in FormFunctionLocal()
341 uyy = (2.0 * u - x[j - 1][i].v - x[j + 1][i].v) * hxdhy; in FormFunctionLocal()
342 f[j][i].v = uxx + uyy + .5 * (x[j][i + 1].omega - x[j][i - 1].omega) * hy; in FormFunctionLocal()
345 u = x[j][i].omega; in FormFunctionLocal()
346 uxx = (2.0 * u - x[j][i - 1].omega - x[j][i + 1].omega) * hydhx; in FormFunctionLocal()
347 uyy = (2.0 * u - x[j - 1][i].omega - x[j + 1][i].omega) * hxdhy; in FormFunctionLocal()
348 …j][i].omega = uxx + uyy + (vxp * (u - x[j][i - 1].omega) + vxm * (x[j][i + 1].omega - u)) * hy + (… in FormFunctionLocal()
351 u = x[j][i].temp; in FormFunctionLocal()
352 uxx = (2.0 * u - x[j][i - 1].temp - x[j][i + 1].temp) * hydhx; in FormFunctionLocal()
353 uyy = (2.0 * u - x[j - 1][i].temp - x[j + 1][i].temp) * hxdhy; in FormFunctionLocal()
354 …f[j][i].temp = uxx + uyy + prandtl * ((vxp * (u - x[j][i - 1].temp) + vxm * (x[j][i + 1].temp - u)… in FormFunctionLocal()
374 PetscInt xints, xinte, yints, yinte, i, j, k, l; in NonlinearGS() local
431 j = 0; in NonlinearGS()
435 bjiu = b[j][i].u; in NonlinearGS()
436 bjiv = b[j][i].v; in NonlinearGS()
441 x[j][i].u = 0.0 + bjiu; in NonlinearGS()
442 x[j][i].v = 0.0 + bjiv; in NonlinearGS()
448 j = info.my - 1; in NonlinearGS()
452 bjiu = b[j][i].u; in NonlinearGS()
453 bjiv = b[j][i].v; in NonlinearGS()
458 x[j][i].u = lid + bjiu; in NonlinearGS()
459 x[j][i].v = bjiv; in NonlinearGS()
467 for (j = info.ys; j < info.ys + info.ym; j++) { in NonlinearGS()
469 bjiu = b[j][i].u; in NonlinearGS()
470 bjiv = b[j][i].v; in NonlinearGS()
475 x[j][i].u = 0.0 + bjiu; in NonlinearGS()
476 x[j][i].v = 0.0 + bjiv; in NonlinearGS()
484 for (j = info.ys; j < info.ys + info.ym; j++) { in NonlinearGS()
486 bjiu = b[j][i].u; in NonlinearGS()
487 bjiv = b[j][i].v; in NonlinearGS()
492 x[j][i].u = 0.0 + bjiu; in NonlinearGS()
493 x[j][i].v = 0.0 + bjiv; in NonlinearGS()
498 for (j = info.ys; j < info.ys + info.ym; j++) { in NonlinearGS()
509 bjiu = b[j][i].u; in NonlinearGS()
510 bjiv = b[j][i].v; in NonlinearGS()
511 bjiomega = b[j][i].omega; in NonlinearGS()
512 bjitemp = b[j][i].temp; in NonlinearGS()
520 if (i != 0 && i != info.mx - 1 && j != 0 && j != info.my - 1) { in NonlinearGS()
522 u = x[j][i].u; in NonlinearGS()
523 uxx = (2.0 * u - x[j][i - 1].u - x[j][i + 1].u) * hydhx; in NonlinearGS()
524 uyy = (2.0 * u - x[j - 1][i].u - x[j + 1][i].u) * hxdhy; in NonlinearGS()
525 fu = uxx + uyy - .5 * (x[j + 1][i].omega - x[j - 1][i].omega) * hx - bjiu; in NonlinearGS()
528 u = x[j][i].v; in NonlinearGS()
529 uxx = (2.0 * u - x[j][i - 1].v - x[j][i + 1].v) * hydhx; in NonlinearGS()
530 uyy = (2.0 * u - x[j - 1][i].v - x[j + 1][i].v) * hxdhy; in NonlinearGS()
531 fv = uxx + uyy + .5 * (x[j][i + 1].omega - x[j][i - 1].omega) * hy - bjiv; in NonlinearGS()
536 vx = x[j][i].u; in NonlinearGS()
540 vy = x[j][i].v; in NonlinearGS()
545 u = x[j][i].omega; in NonlinearGS()
546 uxx = (2.0 * u - x[j][i - 1].omega - x[j][i + 1].omega) * hydhx; in NonlinearGS()
547 uyy = (2.0 * u - x[j - 1][i].omega - x[j + 1][i].omega) * hxdhy; in NonlinearGS()
548 …j][i - 1].omega) + vxm * (x[j][i + 1].omega - u)) * hy + (vyp * (u - x[j - 1][i].omega) + vym * (x… in NonlinearGS()
551 if (PetscRealPart(vx) > 0.0) dfodu = (u - x[j][i - 1].omega) * hy; in NonlinearGS()
552 else dfodu = (x[j][i + 1].omega - u) * hy; in NonlinearGS()
554 if (PetscRealPart(vy) > 0.0) dfodv = (u - x[j - 1][i].omega) * hx; in NonlinearGS()
555 else dfodv = (x[j + 1][i].omega - u) * hx; in NonlinearGS()
558 u = x[j][i].temp; in NonlinearGS()
559 uxx = (2.0 * u - x[j][i - 1].temp - x[j][i + 1].temp) * hydhx; in NonlinearGS()
560 uyy = (2.0 * u - x[j - 1][i].temp - x[j + 1][i].temp) * hxdhy; in NonlinearGS()
561 … + prandtl * ((vxp * (u - x[j][i - 1].temp) + vxm * (x[j][i + 1].temp - u)) * hy + (vyp * (u - x[j… in NonlinearGS()
563 if (PetscRealPart(vx) > 0.0) dftdu = prandtl * (u - x[j][i - 1].temp) * hy; in NonlinearGS()
564 else dftdu = prandtl * (x[j][i + 1].temp - u) * hy; in NonlinearGS()
566 if (PetscRealPart(vy) > 0.0) dftdv = prandtl * (u - x[j - 1][i].temp) * hx; in NonlinearGS()
567 else dftdv = prandtl * (x[j + 1][i].temp - u) * hx; in NonlinearGS()
581 x[j][i].u = x[j][i].u - yu; in NonlinearGS()
582 x[j][i].v = x[j][i].v - yv; in NonlinearGS()
583 x[j][i].temp = x[j][i].temp - yt; in NonlinearGS()
584 x[j][i].omega = x[j][i].omega - yo; in NonlinearGS()
587 fomega = x[j][i].omega - (x[j][i + 1].v - x[j][i].v) * dhx - bjiomega; in NonlinearGS()
588 ftemp = x[j][i].temp - bjitemp; in NonlinearGS()
591 x[j][i].omega = x[j][i].omega - fomega; in NonlinearGS()
592 x[j][i].temp = x[j][i].temp - ftemp; in NonlinearGS()
595 fomega = x[j][i].omega - (x[j][i].v - x[j][i - 1].v) * dhx - bjiomega; in NonlinearGS()
596 ftemp = x[j][i].temp - (PetscReal)(grashof > 0) - bjitemp; in NonlinearGS()
599 x[j][i].omega = x[j][i].omega - fomega; in NonlinearGS()
600 x[j][i].temp = x[j][i].temp - ftemp; in NonlinearGS()
602 if (j == 0) { in NonlinearGS()
603 fomega = x[j][i].omega + (x[j + 1][i].u - x[j][i].u) * dhy - bjiomega; in NonlinearGS()
604 ftemp = x[j][i].temp - x[j + 1][i].temp - bjitemp; in NonlinearGS()
607 x[j][i].omega = x[j][i].omega - fomega; in NonlinearGS()
608 x[j][i].temp = x[j][i].temp - ftemp; in NonlinearGS()
610 if (j == info.my - 1) { in NonlinearGS()
611 fomega = x[j][i].omega + (x[j][i].u - x[j - 1][i].u) * dhy - bjiomega; in NonlinearGS()
612 ftemp = x[j][i].temp - x[j - 1][i].temp - bjitemp; in NonlinearGS()
615 x[j][i].omega = x[j][i].omega - fomega; in NonlinearGS()
616 x[j][i].temp = x[j][i].temp - ftemp; in NonlinearGS()
623 …pxnorm = PetscRealPart(x[j][i].u * x[j][i].u + x[j][i].v * x[j][i].v + x[j][i].omega * x[j][i].ome… in NonlinearGS()