45   PetscReal         ftime, hx, dt;  in main()  local
106   hx = 1.0 / (PetscReal)(mx / 2 - 1);  in main()
107   dt = 0.4 * PetscSqr(hx) / user.alpha; /* Diffusive stability limit */  in main()
142   PetscReal     hx;  in FormIFunction()  local
148   hx = 1.0 / (PetscReal)(info.mx - 1);  in FormIFunction()
168       f[i].u = hx * (x[i].u - user->uleft);  in FormIFunction()
169       f[i].v = hx * (x[i].v - user->vleft);  in FormIFunction()
171       f[i].u = hx * (x[i].u - user->uright);  in FormIFunction()
172       f[i].v = hx * (x[i].v - user->vright);  in FormIFunction()
174       f[i].u = hx * xdot[i].u - user->alpha * (x[i - 1].u - 2. * x[i].u + x[i + 1].u) / hx;  in FormIFunction()
175       f[i].v = hx * xdot[i].v - user->alpha * (x[i - 1].v - 2. * x[i].v + x[i + 1].v) / hx;  in FormIFunction()
193   PetscReal     hx;  in FormRHSFunction()  local
199   hx = 1.0 / (PetscReal)(info.mx - 1);  in FormRHSFunction()
208     f[i].u = hx * (user->A + u * u * v - (user->B + 1) * u);  in FormRHSFunction()
209     f[i].v = hx * (user->B * u - u * u * v);  in FormRHSFunction()
227   PetscReal     hx;  in FormIJacobian()  local
234   hx = 1.0 / (PetscReal)(info.mx - 1);  in FormIJacobian()
245         {hx, 0 },  in FormIJacobian()
246         {0,  hx}  in FormIJacobian()
251 …const PetscScalar dxxL = -user->alpha / hx, dxx0 = 2. * user->alpha / hx, dxxR = -user->alpha / hx;  in FormIJacobian()
253         {{dxxL, 0}, {a * hx + dxx0, 0}, {dxxR, 0}},  in FormIJacobian()
254         {{0, dxxL}, {0, a * hx + dxx0}, {0, dxxR}}  in FormIJacobian()
280   PetscReal     hx;  in FormInitialSolution()  local
285   hx = 1.0 / (PetscReal)(info.mx - 1);  in FormInitialSolution()
292     PetscReal xi = i * hx;  in FormInitialSolution()