63   PetscReal         ftime, hx, dt, xmax, xmin;  in Brusselator()  local
124   hx = 1.0 / (PetscReal)(mx / 2 - 1);  in Brusselator()
125   dt = 0.4 * PetscSqr(hx) / user.alpha; /* Diffusive stability limit */  in Brusselator()
164   PetscReal     hx;  in FormIFunction()  local
170   hx = 1.0 / (PetscReal)(info.mx - 1);  in FormIFunction()
190       f[i].u = hx * (x[i].u - user->uleft);  in FormIFunction()
191       f[i].v = hx * (x[i].v - user->vleft);  in FormIFunction()
193       f[i].u = hx * (x[i].u - user->uright);  in FormIFunction()
194       f[i].v = hx * (x[i].v - user->vright);  in FormIFunction()
196       f[i].u = hx * xdot[i].u - user->alpha * (x[i - 1].u - 2. * x[i].u + x[i + 1].u) / hx;  in FormIFunction()
197       f[i].v = hx * xdot[i].v - user->alpha * (x[i - 1].v - 2. * x[i].v + x[i + 1].v) / hx;  in FormIFunction()
215   PetscReal     hx;  in FormRHSFunction()  local
221   hx = 1.0 / (PetscReal)(info.mx - 1);  in FormRHSFunction()
230     f[i].u = hx * (user->A + u * u * v - (user->B + 1) * u);  in FormRHSFunction()
231     f[i].v = hx * (user->B * u - u * u * v);  in FormRHSFunction()
249   PetscReal     hx;  in FormIJacobian()  local
256   hx = 1.0 / (PetscReal)(info.mx - 1);  in FormIJacobian()
267         {hx, 0 },  in FormIJacobian()
268         {0,  hx}  in FormIJacobian()
273 …const PetscScalar dxxL = -user->alpha / hx, dxx0 = 2. * user->alpha / hx, dxxR = -user->alpha / hx;  in FormIJacobian()
275         {{dxxL, 0}, {a * hx + dxx0, 0}, {dxxR, 0}},  in FormIJacobian()
276         {{0, dxxL}, {0, a * hx + dxx0}, {0, dxxR}}  in FormIJacobian()
302   PetscReal     hx;  in FormInitialSolution()  local
307   hx = 1.0 / (PetscReal)(info.mx - 1);  in FormInitialSolution()
314     PetscReal xi = i * hx;  in FormInitialSolution()