Lines Matching refs:i
200 PetscInt i, Mx, xs, xm; in FormFunction() local
238 for (i = xs; i < xs + xm; i++) { in FormFunction()
240 c = (1. - x[i] * x[i]) * (x[i - 1] + x[i + 1] - 2.0 * x[i]) * sx; in FormFunction()
241 r = (1. - x[i + 1] * x[i + 1]) * (x[i] + x[i + 2] - 2.0 * x[i + 1]) * sx; in FormFunction()
242 l = (1. - x[i - 1] * x[i - 1]) * (x[i - 2] + x[i] - 2.0 * x[i - 1]) * sx; in FormFunction()
244 c = (x[i - 1] + x[i + 1] - 2.0 * x[i]) * sx; in FormFunction()
245 r = (x[i] + x[i + 2] - 2.0 * x[i + 1]) * sx; in FormFunction()
246 l = (x[i - 2] + x[i] - 2.0 * x[i - 1]) * sx; in FormFunction()
248 f[i] = -ctx->kappa * (l + r - 2.0 * c) * sx; in FormFunction()
252 …f[i] += 6. * .25 * x[i] * (x[i + 1] - x[i - 1]) * (x[i + 1] - x[i - 1]) * sx + (3. * x[i] * x[i] -… in FormFunction()
255 f[i] += -(x[i - 1] + x[i + 1] - 2.0 * x[i]) * sx; in FormFunction()
258 …f[i] += 6. * .25 * x[i] * (x[i + 1] - x[i - 1]) * (x[i + 1] - x[i - 1]) * sx + (3. * x[i] * x[i] -… in FormFunction()
260 …if (PetscRealPart(x[i]) < -1.0 + 2.0 * tol) f[i] += (.25 * theta / (tol - tol * tol)) * (x[i - 1] … in FormFunction()
261 …else if (PetscRealPart(x[i]) > 1.0 - 2.0 * tol) f[i] += (.25 * theta / (tol - tol * tol)) * (x[i -… in FormFunction()
262 …i] += 2.0 * theta * x[i] / ((1.0 - x[i] * x[i]) * (1.0 - x[i] * x[i])) * .25 * (x[i + 1] - x[i - 1… in FormFunction()
266 …i]) < -1.0 + 2.0 * tol) f[i] += -1.0 * a * .25 * (x[i + 1] - x[i - 1]) * (x[i + 1] - x[i - 1]) * s… in FormFunction()
267 …x[i]) > 1.0 - 2.0 * tol) f[i] += 1.0 * a * .25 * (x[i + 1] - x[i - 1]) * (x[i + 1] - x[i - 1]) * s… in FormFunction()
268 …i] += 2.0 * theta * x[i] / ((1.0 - x[i] * x[i]) * (1.0 - x[i] * x[i])) * .25 * (x[i + 1] - x[i - 1… in FormFunction()
272 f[i] += -theta_c * (x[i - 1] + x[i + 1] - 2.0 * x[i]) * sx; in FormFunction()
274 …if (PetscRealPart(x[i]) < -1.0 + 2.0 * tol) f[i] += (.25 * theta / (tol - tol * tol)) * (x[i - 1] … in FormFunction()
275 …else if (PetscRealPart(x[i]) > 1.0 - 2.0 * tol) f[i] += (.25 * theta / (tol - tol * tol)) * (x[i -… in FormFunction()
276 …i] += 2.0 * theta * x[i] / ((1.0 - x[i] * x[i]) * (1.0 - x[i] * x[i])) * .25 * (x[i + 1] - x[i - 1… in FormFunction()
280 …i]) < -1.0 + 2.0 * tol) f[i] += -1.0 * a * .25 * (x[i + 1] - x[i - 1]) * (x[i + 1] - x[i - 1]) * s… in FormFunction()
281 …x[i]) > 1.0 - 2.0 * tol) f[i] += 1.0 * a * .25 * (x[i + 1] - x[i - 1]) * (x[i + 1] - x[i - 1]) * s… in FormFunction()
282 …i] += 2.0 * theta * x[i] / ((1.0 - x[i] * x[i]) * (1.0 - x[i] * x[i])) * .25 * (x[i + 1] - x[i - 1… in FormFunction()
306 PetscInt i, Mx, xs, xm; in FormJacobian() local
343 for (i = xs; i < xs + xm; i++) { in FormJacobian()
344 row.i = i; in FormJacobian()
351 cols[0].i = i - 2; in FormJacobian()
353 cols[1].i = i - 1; in FormJacobian()
355 cols[2].i = i; in FormJacobian()
357 cols[3].i = i + 1; in FormJacobian()
359 cols[4].i = i + 2; in FormJacobian()
396 PetscInt i, xs, xm, Mx, N, scale; in FormInitialSolution() local
432 for (i = xs; i < xs + xm; i++) { in FormInitialSolution()
433 x = i * hx; in FormInitialSolution()
435 if (r < .125) u[i] = 1.0; in FormInitialSolution()
436 else u[i] = -.5; in FormInitialSolution()
441 u[i] = f[scale * i]; in FormInitialSolution()
461 PetscInt Mx, i, xs, xm, cnt; in MyMonitor() local
510 for (i = xs; i < xs + xm; i++) { in MyMonitor()
512 …nerate) yy[0] = PetscRealPart(.25 * (1. - u[i] * u[i]) * ctx->kappa * (u[i - 1] - u[i + 1]) * (u[i… in MyMonitor()
513 … else yy[0] = PetscRealPart(.25 * ctx->kappa * (u[i - 1] - u[i + 1]) * (u[i - 1] - u[i + 1]) * sx); in MyMonitor()
518 yy[1] = .25 * PetscRealPart((1. - u[i] * u[i]) * (1. - u[i] * u[i])); in MyMonitor()
521 yy[1] = .5 * PetscRealPart(1. - u[i] * u[i]); in MyMonitor()
524 yy[1] = .25 * PetscRealPart((1. - u[i] * u[i]) * (1. - u[i] * u[i])); in MyMonitor()
525 …art(u[i]) < -1.0 + 2.0 * tol) yy[2] = .5 * theta * (2.0 * tol * PetscLogReal(tol) + PetscRealPart(… in MyMonitor()
526 …e if (PetscRealPart(u[i]) > 1.0 - 2.0 * tol) yy[2] = .5 * theta * (PetscRealPart(1.0 + u[i]) * Pet… in MyMonitor()
527 …RealPart(1.0 + u[i]) * PetscLogReal(PetscRealPart(1.0 + u[i]) / 2.0) + PetscRealPart(1.0 - u[i]) *… in MyMonitor()
530 yy[1] = .5 * theta_c * PetscRealPart(1.0 - u[i] * u[i]); in MyMonitor()
531 …art(u[i]) < -1.0 + 2.0 * tol) yy[2] = .5 * theta * (2.0 * tol * PetscLogReal(tol) + PetscRealPart(… in MyMonitor()
532 …e if (PetscRealPart(u[i]) > 1.0 - 2.0 * tol) yy[2] = .5 * theta * (PetscRealPart(1.0 + u[i]) * Pet… in MyMonitor()
533 …RealPart(1.0 + u[i]) * PetscLogReal(PetscRealPart(1.0 + u[i]) / 2.0) + PetscRealPart(1.0 - u[i]) *… in MyMonitor()
557 for (i = xs; i < xs + xm; i++) { in MyMonitor()
561 c = (1. - u[i] * u[i]) * (u[i - 1] + u[i + 1] - 2.0 * u[i]) * sx; in MyMonitor()
562 r = (1. - u[i + 1] * u[i + 1]) * (u[i] + u[i + 2] - 2.0 * u[i + 1]) * sx; in MyMonitor()
563 l = (1. - u[i - 1] * u[i - 1]) * (u[i - 2] + u[i] - 2.0 * u[i - 1]) * sx; in MyMonitor()
565 c = (u[i - 1] + u[i + 1] - 2.0 * u[i]) * sx; in MyMonitor()
566 r = (u[i] + u[i + 2] - 2.0 * u[i + 1]) * sx; in MyMonitor()
567 l = (u[i - 2] + u[i] - 2.0 * u[i - 1]) * sx; in MyMonitor()
575 …Part(6. * .25 * u[i] * (u[i + 1] - u[i - 1]) * (u[i + 1] - u[i - 1]) * sx + (3. * u[i] * u[i] - 1.… in MyMonitor()
578 yy[1] = -PetscRealPart(u[i - 1] + u[i + 1] - 2.0 * u[i]) * sx; in MyMonitor()
581 …Part(6. * .25 * u[i] * (u[i + 1] - u[i - 1]) * (u[i + 1] - u[i - 1]) * sx + (3. * u[i] * u[i] - 1.… in MyMonitor()
583 …if (PetscRealPart(u[i]) < -1.0 + 2.0 * tol) yy[2] = (.25 * theta / (tol - tol * tol)) * PetscRealP… in MyMonitor()
584 …if (PetscRealPart(u[i]) > 1.0 - 2.0 * tol) yy[2] = (.25 * theta / (tol - tol * tol)) * PetscRealPa… in MyMonitor()
585 …i] / ((1.0 - u[i] * u[i]) * (1.0 - u[i] * u[i])) * .25 * (u[i + 1] - u[i - 1]) * (u[i + 1] - u[i -… in MyMonitor()
589 …i]) < -1.0 + 2.0 * tol) yy[2] = PetscRealPart(-1.0 * a * .25 * (u[i + 1] - u[i - 1]) * (u[i + 1] -… in MyMonitor()
590 …i]) > 1.0 - 2.0 * tol) yy[2] = PetscRealPart(1.0 * a * .25 * (u[i + 1] - u[i - 1]) * (u[i + 1] - u… in MyMonitor()
591 …i] / ((1.0 - u[i] * u[i]) * (1.0 - u[i] * u[i])) * .25 * (u[i + 1] - u[i - 1]) * (u[i + 1] - u[i -… in MyMonitor()
595 yy[1] = theta_c * PetscRealPart(-(u[i - 1] + u[i + 1] - 2.0 * u[i])) * sx; in MyMonitor()
597 …if (PetscRealPart(u[i]) < -1.0 + 2.0 * tol) yy[2] = (.25 * theta / (tol - tol * tol)) * PetscRealP… in MyMonitor()
598 …if (PetscRealPart(u[i]) > 1.0 - 2.0 * tol) yy[2] = (.25 * theta / (tol - tol * tol)) * PetscRealPa… in MyMonitor()
599 …i] / ((1.0 - u[i] * u[i]) * (1.0 - u[i] * u[i])) * .25 * (u[i + 1] - u[i - 1]) * (u[i + 1] - u[i -… in MyMonitor()
603 …i]) < -1.0 + 2.0 * tol) yy[2] = PetscRealPart(-1.0 * a * .25 * (u[i + 1] - u[i - 1]) * (u[i + 1] -… in MyMonitor()
604 …i]) > 1.0 - 2.0 * tol) yy[2] = PetscRealPart(1.0 * a * .25 * (u[i + 1] - u[i - 1]) * (u[i + 1] - u… in MyMonitor()
605 …i] / ((1.0 - u[i] * u[i]) * (1.0 - u[i] * u[i])) * .25 * (u[i + 1] - u[i - 1]) * (u[i + 1] - u[i -… in MyMonitor()
643 for (i = xs; i < xs + xm; i++) { in MyMonitor()
645 yy[0] = PetscRealPart(u[i]); in MyMonitor()
659 for (i = xs; i < xs + xm; i += cnt) { in MyMonitor()
660 y = yup = ydown = PetscRealPart(u[i]); in MyMonitor()
661 c = (u[i - 1] + u[i + 1] - 2.0 * u[i]) * sx; in MyMonitor()
662 r = (u[i] + u[i + 2] - 2.0 * u[i + 1]) * sx; in MyMonitor()
663 l = (u[i - 2] + u[i] - 2.0 * u[i - 1]) * sx; in MyMonitor()
672 …Part(6. * .25 * u[i] * (u[i + 1] - u[i - 1]) * (u[i + 1] - u[i - 1]) * sx + (3. * u[i] * u[i] - 1.… in MyMonitor()
675 len = -.5 * PetscRealPart(u[i - 1] + u[i + 1] - 2.0 * u[i]) * sx / max; in MyMonitor()
678 …Part(6. * .25 * u[i] * (u[i + 1] - u[i - 1]) * (u[i + 1] - u[i - 1]) * sx + (3. * u[i] * u[i] - 1.… in MyMonitor()
683 … (PetscRealPart(u[i]) < -1.0 + 2.0 * tol) len2 = .5 * (.25 * theta / (tol - tol * tol)) * PetscRea… in MyMonitor()
684 … (PetscRealPart(u[i]) > 1.0 - 2.0 * tol) len2 = .5 * (.25 * theta / (tol - tol * tol)) * PetscReal… in MyMonitor()
685 …i] / ((1.0 - u[i] * u[i]) * (1.0 - u[i] * u[i])) * .25 * (u[i + 1] - u[i - 1]) * (u[i + 1] - u[i -… in MyMonitor()
689 …i]) < -1.0 + 2.0 * tol) len2 = PetscRealPart(.5 * (-1.0 * a * .25 * (u[i + 1] - u[i - 1]) * (u[i +… in MyMonitor()
690 …i]) > 1.0 - 2.0 * tol) len2 = PetscRealPart(.5 * (a * .25 * (u[i + 1] - u[i - 1]) * (u[i + 1] - u[… in MyMonitor()
691 …i] / ((1.0 - u[i] * u[i]) * (1.0 - u[i] * u[i])) * .25 * (u[i + 1] - u[i - 1]) * (u[i + 1] - u[i -… in MyMonitor()
697 len = -.5 * theta_c * PetscRealPart(-(u[i - 1] + u[i + 1] - 2.0 * u[i]) * sx / max); in MyMonitor()
702 … (PetscRealPart(u[i]) < -1.0 + 2.0 * tol) len2 = .5 * (.25 * theta / (tol - tol * tol)) * PetscRea… in MyMonitor()
703 … (PetscRealPart(u[i]) > 1.0 - 2.0 * tol) len2 = .5 * (.25 * theta / (tol - tol * tol)) * PetscReal… in MyMonitor()
704 …i] / ((1.0 - u[i] * u[i]) * (1.0 - u[i] * u[i])) * .25 * (u[i + 1] - u[i - 1]) * (u[i + 1] - u[i -… in MyMonitor()
708 …i]) < -1.0 + 2.0 * tol) len2 = .5 * PetscRealPart(-1.0 * a * .25 * (u[i + 1] - u[i - 1]) * (u[i + … in MyMonitor()
709 …i]) > 1.0 - 2.0 * tol) len2 = .5 * PetscRealPart(a * .25 * (u[i + 1] - u[i - 1]) * (u[i + 1] - u[i… in MyMonitor()
710 …i] / ((1.0 - u[i] * u[i]) * (1.0 - u[i] * u[i])) * .25 * (u[i + 1] - u[i - 1]) * (u[i + 1] - u[i -… in MyMonitor()