Lines Matching refs:j
286 PetscInt i, j, Mx, My, xs, ys, xm, ym; in FormInitialGuess() local
317 for (j = ys; j < ys + ym; j++) { in FormInitialGuess()
318 temp = (PetscReal)(PetscMin(j, My - j - 1)) * hy; in FormInitialGuess()
320 if (i == 0 || j == 0 || i == Mx - 1 || j == My - 1) { in FormInitialGuess()
322 x[j][i] = 0.0; in FormInitialGuess()
326 … x[j][i] = temp1 * PetscSqrtReal(PetscMin((PetscReal)(PetscMin(i, Mx - i - 1)) * hx, temp)); in FormInitialGuess()
330 … const PetscReal xx = 2 * (PetscReal)i / (Mx - 1) - 1, yy = 2 * (PetscReal)j / (My - 1) - 1; in FormInitialGuess()
331 x[j][i] = (1 - xx * xx) * (1 - yy * yy) * xx * yy; in FormInitialGuess()
334 x[j][i] = 0.; in FormInitialGuess()
336 … const PetscReal xx = 2 * (PetscReal)i / (Mx - 1) - 1, yy = 2 * (PetscReal)j / (My - 1) - 1; in FormInitialGuess()
337 x[j][i] = (1 - xx * xx) * (1 - yy * yy) * xx * yy; in FormInitialGuess()
340 … x[j][i] = temp1 * PetscSqrtReal(PetscMin((PetscReal)(PetscMin(i, Mx - i - 1)) * hx, temp)); in FormInitialGuess()
342 … x[j][i] = 0.5 * PetscSqrtReal(PetscMin((PetscReal)(PetscMin(i, Mx - i - 1)) * hx, temp)); in FormInitialGuess()
367 PetscInt i, j, Mx, My, xs, ys, xm, ym; in FormRHS() local
378 for (j = ys; j < ys + ym; j++) { in FormRHS()
380 if (i == 0 || j == 0 || i == Mx - 1 || j == My - 1) { in FormRHS()
381 b[j][i] = 0.0; in FormRHS()
383 b[j][i] = hx * hy * user->source; in FormRHS()
412 PetscInt i, j; in FormFunctionLocal() local
424 for (j = info->ys; j < info->ys + info->ym; j++) { in FormFunctionLocal()
426 PetscReal xx = i * hx, yy = j * hy; in FormFunctionLocal()
427 if (i == 0 || j == 0 || i == info->mx - 1 || j == info->my - 1) { in FormFunctionLocal()
428 f[j][i] = x[j][i]; in FormFunctionLocal()
430 …j][i], ux_E = dhx * (x[j][i + 1] - x[j][i]), uy_E = 0.25 * dhy * (x[j + 1][i] + x[j + 1][i + 1] - … in FormFunctionLocal()
437 f[j][i] = uxx + uyy - sc * eu; in FormFunctionLocal()
453 PetscInt i, j; in FormFunctionPicardLocal() local
462 for (j = info->ys; j < info->ys + info->ym; j++) { in FormFunctionPicardLocal()
464 if (!(i == 0 || j == 0 || i == info->mx - 1 || j == info->my - 1)) { in FormFunctionPicardLocal()
465 const PetscScalar u = x[j][i]; in FormFunctionPicardLocal()
466 f[j][i] = sc * PetscExpScalar(u); in FormFunctionPicardLocal()
468 … f[j][i] = 0.0; /* this is zero because the A(x) x term forces the x to be zero on the boundary */ in FormFunctionPicardLocal()
481 PetscInt i, j; in FormJacobianLocal() local
505 for (j = info->ys; j < info->ys + info->ym; j++) { in FormJacobianLocal()
507 PetscReal xx = i * hx, yy = j * hy; in FormJacobianLocal()
508 row.j = j; in FormJacobianLocal()
511 if (i == 0 || j == 0 || i == info->mx - 1 || j == info->my - 1) { in FormJacobianLocal()
516 …j][i + 1] - x[j][i]), uy_E = 0.25 * dhy * (x[j + 1][i] + x[j + 1][i + 1] - x[j - 1][i] - x[j - 1][… in FormJacobianLocal()
522 col[0].j = j - 1; in FormJacobianLocal()
525 col[1].j = j; in FormJacobianLocal()
528 col[2].j = row.j; in FormJacobianLocal()
531 col[3].j = j; in FormJacobianLocal()
534 col[4].j = j + 1; in FormJacobianLocal()
541 col[0].j = j - 1; in FormJacobianLocal()
544 col[1].j = j; in FormJacobianLocal()
547 col[2].j = row.j; in FormJacobianLocal()
550 col[3].j = j; in FormJacobianLocal()
553 col[4].j = j + 1; in FormJacobianLocal()
559 col[0].j = j - 1; in FormJacobianLocal()
561 col[1].j = j; in FormJacobianLocal()
563 col[2].j = j; in FormJacobianLocal()
565 col[3].j = j; in FormJacobianLocal()
567 col[4].j = j + 1; in FormJacobianLocal()
582 col[0].j = j - 1; in FormJacobianLocal()
584 col[1].j = j - 1; in FormJacobianLocal()
586 col[2].j = j - 1; in FormJacobianLocal()
588 col[3].j = j; in FormJacobianLocal()
590 col[4].j = j; in FormJacobianLocal()
592 col[5].j = j; in FormJacobianLocal()
594 col[6].j = j + 1; in FormJacobianLocal()
596 col[7].j = j + 1; in FormJacobianLocal()
598 col[8].j = j + 1; in FormJacobianLocal()
728 PetscInt i, j, k, xs, ys, xm, ym, its, tot_its, sweeps, l, m; in NonlinearGS() local
770 for (j = ys; j < ys + ym; j++) { in NonlinearGS()
771 for (i = xs + (m + j) % 2; i < xs + xm; i += 2) { in NonlinearGS()
772 PetscReal xx = i * hx, yy = j * hy; in NonlinearGS()
773 if (B) bij = b[j][i]; in NonlinearGS()
776 if (i == 0 || j == 0 || i == info.mx - 1 || j == info.my - 1) { in NonlinearGS()
778 x[j][i] = 0.0 + bij; in NonlinearGS()
780 … const PetscScalar u_E = x[j][i + 1], u_W = x[j][i - 1], u_N = x[j + 1][i], u_S = x[j - 1][i]; in NonlinearGS()
781 …j + 1][i] + x[j + 1][i + 1] - x[j - 1][i] - x[j - 1][i + 1]), uy_W = 0.25 * dhy * (x[j + 1][i - 1]… in NonlinearGS()
782 u = x[j][i]; in NonlinearGS()
797 x[j][i] = u; in NonlinearGS()