Lines Matching refs:sx
201 PetscReal hx, sx; in FormFunction() local
213 sx = 1.0 / (hx * hx); 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 …1] - x[i - 1]) * (x[i + 1] - x[i - 1]) * sx + (3. * x[i] * x[i] - 1.) * (x[i - 1] + x[i + 1] - 2.0… in FormFunction()
255 f[i] += -(x[i - 1] + x[i + 1] - 2.0 * x[i]) * sx; in FormFunction()
258 …1] - x[i - 1]) * (x[i + 1] - x[i - 1]) * sx + (3. * x[i] * x[i] - 1.) * (x[i - 1] + x[i + 1] - 2.0… in FormFunction()
260 …0 + 2.0 * tol) f[i] += (.25 * theta / (tol - tol * tol)) * (x[i - 1] + x[i + 1] - 2.0 * x[i]) * sx; in FormFunction()
261 …0 - 2.0 * tol) f[i] += (.25 * theta / (tol - tol * tol)) * (x[i - 1] + x[i + 1] - 2.0 * x[i]) * sx; in FormFunction()
262 …- x[i - 1]) * (x[i + 1] - x[i - 1]) * sx + (theta / (1.0 - x[i] * x[i])) * (x[i - 1] + x[i + 1] - … in FormFunction()
266 … 1] - x[i - 1]) * (x[i + 1] - x[i - 1]) * sx + (-1.0 * a * x[i] + b) * (x[i - 1] + x[i + 1] - 2.0 … in FormFunction()
267 …[i + 1] - x[i - 1]) * (x[i + 1] - x[i - 1]) * sx + (a * x[i] + b) * (x[i - 1] + x[i + 1] - 2.0 * x… in FormFunction()
268 …- x[i - 1]) * (x[i + 1] - x[i - 1]) * sx + (theta / (1.0 - x[i] * x[i])) * (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 …0 + 2.0 * tol) f[i] += (.25 * theta / (tol - tol * tol)) * (x[i - 1] + x[i + 1] - 2.0 * x[i]) * sx; in FormFunction()
275 …0 - 2.0 * tol) f[i] += (.25 * theta / (tol - tol * tol)) * (x[i - 1] + x[i + 1] - 2.0 * x[i]) * sx; in FormFunction()
276 …- x[i - 1]) * (x[i + 1] - x[i - 1]) * sx + (theta / (1.0 - x[i] * x[i])) * (x[i - 1] + x[i + 1] - … in FormFunction()
280 … 1] - x[i - 1]) * (x[i + 1] - x[i - 1]) * sx + (-1.0 * a * x[i] + b) * (x[i - 1] + x[i + 1] - 2.0 … in FormFunction()
281 …[i + 1] - x[i - 1]) * (x[i + 1] - x[i - 1]) * sx + (a * x[i] + b) * (x[i - 1] + x[i + 1] - 2.0 * x… in FormFunction()
282 …- x[i - 1]) * (x[i + 1] - x[i - 1]) * sx + (theta / (1.0 - x[i] * x[i])) * (x[i - 1] + x[i + 1] - … in FormFunction()
308 PetscReal hx, sx; in FormJacobian() local
319 sx = 1.0 / (hx * hx); in FormJacobian()
352 vals[0] = -ctx->kappa * sx * sx; in FormJacobian()
354 vals[1] = 4.0 * ctx->kappa * sx * sx; in FormJacobian()
356 vals[2] = -6.0 * ctx->kappa * sx * sx; in FormJacobian()
358 vals[3] = 4.0 * ctx->kappa * sx * sx; in FormJacobian()
360 vals[4] = -ctx->kappa * sx * sx; in FormJacobian()
462 …PetscReal x, y, hx, pause, sx, len, max, xx[4], yy[4], xx_netforce, yy_netforce, yup, ydown, y2… in MyMonitor() local
483 sx = 1.0 / (hx * hx); in MyMonitor()
512 …alPart(.25 * (1. - u[i] * u[i]) * ctx->kappa * (u[i - 1] - u[i + 1]) * (u[i - 1] - u[i + 1]) * sx); 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()
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()
569 yy[0] = PetscRealPart(-ctx->kappa * (l + r - 2.0 * c) * sx); in MyMonitor()
575 …1] - u[i - 1]) * (u[i + 1] - u[i - 1]) * sx + (3. * u[i] * u[i] - 1.) * (u[i - 1] + u[i + 1] - 2.0… in MyMonitor()
578 yy[1] = -PetscRealPart(u[i - 1] + u[i + 1] - 2.0 * u[i]) * sx; in MyMonitor()
581 …1] - u[i - 1]) * (u[i + 1] - u[i - 1]) * sx + (3. * u[i] * u[i] - 1.) * (u[i - 1] + u[i + 1] - 2.0… in MyMonitor()
583 …) yy[2] = (.25 * theta / (tol - tol * tol)) * PetscRealPart(u[i - 1] + u[i + 1] - 2.0 * u[i]) * sx; in MyMonitor()
584 …) yy[2] = (.25 * theta / (tol - tol * tol)) * PetscRealPart(u[i - 1] + u[i + 1] - 2.0 * u[i]) * sx; in MyMonitor()
585 …- u[i - 1]) * (u[i + 1] - u[i - 1]) * sx + (theta / (1.0 - u[i] * u[i])) * (u[i - 1] + u[i + 1] - … in MyMonitor()
589 … 1] - u[i - 1]) * (u[i + 1] - u[i - 1]) * sx + (-1.0 * a * u[i] + b) * (u[i - 1] + u[i + 1] - 2.0 … in MyMonitor()
590 …[i + 1] - u[i - 1]) * (u[i + 1] - u[i - 1]) * sx + (a * u[i] + b) * (u[i - 1] + u[i + 1] - 2.0 * u… in MyMonitor()
591 …- u[i - 1]) * (u[i + 1] - u[i - 1]) * sx + (theta / (1.0 - u[i] * u[i])) * (u[i - 1] + u[i + 1] - … in MyMonitor()
595 yy[1] = theta_c * PetscRealPart(-(u[i - 1] + u[i + 1] - 2.0 * u[i])) * sx; in MyMonitor()
597 …) yy[2] = (.25 * theta / (tol - tol * tol)) * PetscRealPart(u[i - 1] + u[i + 1] - 2.0 * u[i]) * sx; in MyMonitor()
598 …) yy[2] = (.25 * theta / (tol - tol * tol)) * PetscRealPart(u[i - 1] + u[i + 1] - 2.0 * u[i]) * sx; in MyMonitor()
599 …- u[i - 1]) * (u[i + 1] - u[i - 1]) * sx + (theta / (1.0 - u[i] * u[i])) * (u[i - 1] + u[i + 1] - … in MyMonitor()
603 … 1] - u[i - 1]) * (u[i + 1] - u[i - 1]) * sx + (-1.0 * a * u[i] + b) * (u[i - 1] + u[i + 1] - 2.0 … in MyMonitor()
604 …[i + 1] - u[i - 1]) * (u[i + 1] - u[i - 1]) * sx + (a * u[i] + b) * (u[i - 1] + u[i + 1] - 2.0 * u… in MyMonitor()
605 …- u[i - 1]) * (u[i + 1] - u[i - 1]) * sx + (theta / (1.0 - u[i] * u[i])) * (u[i - 1] + u[i + 1] - … 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()
664 len = -.5 * PetscRealPart(ctx->kappa * (l + r - 2.0 * c) * sx) / max; in MyMonitor()
672 …1] - u[i - 1]) * (u[i + 1] - u[i - 1]) * sx + (3. * u[i] * u[i] - 1.) * (u[i - 1] + u[i + 1] - 2.0… in MyMonitor()
675 len = -.5 * PetscRealPart(u[i - 1] + u[i + 1] - 2.0 * u[i]) * sx / max; in MyMonitor()
678 …1] - u[i - 1]) * (u[i + 1] - u[i - 1]) * sx + (3. * u[i] * u[i] - 1.) * (u[i - 1] + u[i + 1] - 2.0… in MyMonitor()
683 …5 * (.25 * theta / (tol - tol * tol)) * PetscRealPart(u[i - 1] + u[i + 1] - 2.0 * u[i]) * sx / max; in MyMonitor()
684 …5 * (.25 * theta / (tol - tol * tol)) * PetscRealPart(u[i - 1] + u[i + 1] - 2.0 * u[i]) * sx / max; in MyMonitor()
685 …- u[i - 1]) * (u[i + 1] - u[i - 1]) * sx + (theta / (1.0 - u[i] * u[i])) * (u[i - 1] + u[i + 1] - … in MyMonitor()
689 … 1] - u[i - 1]) * (u[i + 1] - u[i - 1]) * sx + (-1.0 * a * u[i] + b) * (u[i - 1] + u[i + 1] - 2.0 … in MyMonitor()
690 …[i + 1] - u[i - 1]) * (u[i + 1] - u[i - 1]) * sx + (a * u[i] + b) * (u[i - 1] + u[i + 1] - 2.0 * u… in MyMonitor()
691 …- u[i - 1]) * (u[i + 1] - u[i - 1]) * sx + (theta / (1.0 - u[i] * u[i])) * (u[i - 1] + u[i + 1] - … in MyMonitor()
697 len = -.5 * theta_c * PetscRealPart(-(u[i - 1] + u[i + 1] - 2.0 * u[i]) * sx / max); in MyMonitor()
702 …5 * (.25 * theta / (tol - tol * tol)) * PetscRealPart(u[i - 1] + u[i + 1] - 2.0 * u[i]) * sx / max; in MyMonitor()
703 …5 * (.25 * theta / (tol - tol * tol)) * PetscRealPart(u[i - 1] + u[i + 1] - 2.0 * u[i]) * sx / max; in MyMonitor()
704 …- u[i - 1]) * (u[i + 1] - u[i - 1]) * sx + (theta / (1.0 - u[i] * u[i])) * (u[i - 1] + u[i + 1] - … in MyMonitor()
708 … 1] - u[i - 1]) * (u[i + 1] - u[i - 1]) * sx + (-1.0 * a * u[i] + b) * (u[i - 1] + u[i + 1] - 2.0 … in MyMonitor()
709 …[i + 1] - u[i - 1]) * (u[i + 1] - u[i - 1]) * sx + (a * u[i] + b) * (u[i - 1] + u[i + 1] - 2.0 * u… in MyMonitor()
710 …- u[i - 1]) * (u[i + 1] - u[i - 1]) * sx + (theta / (1.0 - u[i] * u[i])) * (u[i - 1] + u[i + 1] - … in MyMonitor()