Lines Matching refs:yy
462 …PetscReal x, y, hx, pause, sx, len, max, xx[4], yy[4], xx_netforce, yy_netforce, yup, ydown, y2… in MyMonitor() local
512 …if (ctx->degenerate) yy[0] = PetscRealPart(.25 * (1. - u[i] * u[i]) * ctx->kappa * (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 …if (PetscRealPart(u[i]) < -1.0 + 2.0 * tol) yy[2] = .5 * theta * (2.0 * tol * PetscLogReal(tol) + … in MyMonitor()
526 …else if (PetscRealPart(u[i]) > 1.0 - 2.0 * tol) yy[2] = .5 * theta * (PetscRealPart(1.0 + u[i]) * … in MyMonitor()
527 …else yy[2] = .5 * theta * (PetscRealPart(1.0 + u[i]) * PetscLogReal(PetscRealPart(1.0 + u[i]) / 2.… in MyMonitor()
530 yy[1] = .5 * theta_c * PetscRealPart(1.0 - u[i] * u[i]); in MyMonitor()
531 …if (PetscRealPart(u[i]) < -1.0 + 2.0 * tol) yy[2] = .5 * theta * (2.0 * tol * PetscLogReal(tol) + … in MyMonitor()
532 …else if (PetscRealPart(u[i]) > 1.0 - 2.0 * tol) yy[2] = .5 * theta * (PetscRealPart(1.0 + u[i]) * … in MyMonitor()
533 …else yy[2] = .5 * theta * (PetscRealPart(1.0 + u[i]) * PetscLogReal(PetscRealPart(1.0 + u[i]) / 2.… in MyMonitor()
539 PetscCall(PetscDrawLGAddPoint(lg, xx, yy)); in MyMonitor()
569 yy[0] = PetscRealPart(-ctx->kappa * (l + r - 2.0 * c) * sx); in MyMonitor()
570 yy_netforce = yy[0]; in MyMonitor()
571 max = PetscMax(max, PetscAbs(yy[0])); in MyMonitor()
575 …yy[1] = PetscRealPart(6. * .25 * u[i] * (u[i + 1] - u[i - 1]) * (u[i + 1] - u[i - 1]) * sx + (3. *… in MyMonitor()
578 yy[1] = -PetscRealPart(u[i - 1] + u[i + 1] - 2.0 * u[i]) * sx; in MyMonitor()
581 …yy[1] = PetscRealPart(6. * .25 * u[i] * (u[i + 1] - u[i - 1]) * (u[i + 1] - u[i - 1]) * sx + (3. *… in MyMonitor()
583 …if (PetscRealPart(u[i]) < -1.0 + 2.0 * tol) yy[2] = (.25 * theta / (tol - tol * tol)) * PetscRealP… in MyMonitor()
584 …else if (PetscRealPart(u[i]) > 1.0 - 2.0 * tol) yy[2] = (.25 * theta / (tol - tol * tol)) * PetscR… in MyMonitor()
585 …else yy[2] = PetscRealPart(2.0 * theta * u[i] / ((1.0 - u[i] * u[i]) * (1.0 - u[i] * u[i])) * .25 … in MyMonitor()
589 …if (PetscRealPart(u[i]) < -1.0 + 2.0 * tol) yy[2] = PetscRealPart(-1.0 * a * .25 * (u[i + 1] - u[i… in MyMonitor()
590 …else if (PetscRealPart(u[i]) > 1.0 - 2.0 * tol) yy[2] = PetscRealPart(1.0 * a * .25 * (u[i + 1] - … in MyMonitor()
591 …else yy[2] = PetscRealPart(2.0 * theta * u[i] / ((1.0 - u[i] * u[i]) * (1.0 - u[i] * u[i])) * .25 … 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 …else if (PetscRealPart(u[i]) > 1.0 - 2.0 * tol) yy[2] = (.25 * theta / (tol - tol * tol)) * PetscR… in MyMonitor()
599 …else yy[2] = PetscRealPart(2.0 * theta * u[i] / ((1.0 - u[i] * u[i]) * (1.0 - u[i] * u[i])) * .25 … in MyMonitor()
603 …if (PetscRealPart(u[i]) < -1.0 + 2.0 * tol) yy[2] = PetscRealPart(-1.0 * a * .25 * (u[i + 1] - u[i… in MyMonitor()
604 …else if (PetscRealPart(u[i]) > 1.0 - 2.0 * tol) yy[2] = PetscRealPart(1.0 * a * .25 * (u[i + 1] - … in MyMonitor()
605 …else yy[2] = PetscRealPart(2.0 * theta * u[i] / ((1.0 - u[i] * u[i]) * (1.0 - u[i] * u[i])) * .25 … in MyMonitor()
612 max = PetscMax(max, PetscAbs(yy[1])); in MyMonitor()
613 yy[2] = yy[0] + yy[1]; in MyMonitor()
614 yy_netforce = yy[2]; in MyMonitor()
616 max = PetscMax(max, PetscAbs(yy[1] + yy[2])); in MyMonitor()
617 yy[3] = yy[0] + yy[1] + yy[2]; in MyMonitor()
618 yy_netforce = yy[3]; in MyMonitor()
624 PetscCall(PetscDrawLGAddPoint(lg, xx, yy)); in MyMonitor()
645 yy[0] = PetscRealPart(u[i]); in MyMonitor()
646 PetscCall(PetscDrawLGAddPoint(lg, xx, yy)); in MyMonitor()