211   PetscInt  i, j, mx, xs, ys, xm, ym;  in FormInitialSolution()  local
241   for (j = ys; j < ys + ym; j++) {  in FormInitialSolution()
243       x[j][i].u     = 0.0;  in FormInitialSolution()
244       x[j][i].v     = 0.0;  in FormInitialSolution()
245       x[j][i].omega = 0.0;  in FormInitialSolution()
246       x[j][i].temp  = (grashof > 0) * i * dx;  in FormInitialSolution()
260   PetscInt    xints, xinte, yints, yinte, i, j;  in FormIFunctionLocal()  local
291     j     = 0;  in FormIFunctionLocal()
295       f[j][i].u     = x[j][i].u;  in FormIFunctionLocal()
296       f[j][i].v     = x[j][i].v;  in FormIFunctionLocal()
297       f[j][i].omega = x[j][i].omega + (x[j + 1][i].u - x[j][i].u) * dhy;  in FormIFunctionLocal()
298       f[j][i].temp  = x[j][i].temp - x[j + 1][i].temp;  in FormIFunctionLocal()
304     j     = info->my - 1;  in FormIFunctionLocal()
308       f[j][i].u     = x[j][i].u - lid;  in FormIFunctionLocal()
309       f[j][i].v     = x[j][i].v;  in FormIFunctionLocal()
310       f[j][i].omega = x[j][i].omega + (x[j][i].u - x[j - 1][i].u) * dhy;  in FormIFunctionLocal()
311       f[j][i].temp  = x[j][i].temp - x[j - 1][i].temp;  in FormIFunctionLocal()
320     for (j = info->ys; j < info->ys + info->ym; j++) {  in FormIFunctionLocal()
321       f[j][i].u     = x[j][i].u;  in FormIFunctionLocal()
322       f[j][i].v     = x[j][i].v;  in FormIFunctionLocal()
323       f[j][i].omega = x[j][i].omega - (x[j][i + 1].v - x[j][i].v) * dhx;  in FormIFunctionLocal()
324       f[j][i].temp  = x[j][i].temp;  in FormIFunctionLocal()
333     for (j = info->ys; j < info->ys + info->ym; j++) {  in FormIFunctionLocal()
334       f[j][i].u     = x[j][i].u;  in FormIFunctionLocal()
335       f[j][i].v     = x[j][i].v;  in FormIFunctionLocal()
336       f[j][i].omega = x[j][i].omega - (x[j][i].v - x[j][i - 1].v) * dhx;  in FormIFunctionLocal()
337       f[j][i].temp  = x[j][i].temp - (PetscReal)(grashof > 0);  in FormIFunctionLocal()
342   for (j = yints; j < yinte; j++) {  in FormIFunctionLocal()
347       vx  = x[j][i].u;  in FormIFunctionLocal()
351       vy  = x[j][i].v;  in FormIFunctionLocal()
357       u         = x[j][i].u;  in FormIFunctionLocal()
358       udot      = user->parabolic ? xdot[j][i].u : 0.;  in FormIFunctionLocal()
359       uxx       = (2.0 * u - x[j][i - 1].u - x[j][i + 1].u) * hydhx;  in FormIFunctionLocal()
360       uyy       = (2.0 * u - x[j - 1][i].u - x[j + 1][i].u) * hxdhy;  in FormIFunctionLocal()
361       f[j][i].u = udot + uxx + uyy - .5 * (x[j + 1][i].omega - x[j - 1][i].omega) * hx;  in FormIFunctionLocal()
364       u         = x[j][i].v;  in FormIFunctionLocal()
365       udot      = user->parabolic ? xdot[j][i].v : 0.;  in FormIFunctionLocal()
366       uxx       = (2.0 * u - x[j][i - 1].v - x[j][i + 1].v) * hydhx;  in FormIFunctionLocal()
367       uyy       = (2.0 * u - x[j - 1][i].v - x[j + 1][i].v) * hxdhy;  in FormIFunctionLocal()
368       f[j][i].v = udot + uxx + uyy + .5 * (x[j][i + 1].omega - x[j][i - 1].omega) * hy;  in FormIFunctionLocal()
371       u   = x[j][i].omega;  in FormIFunctionLocal()
372       uxx = (2.0 * u - x[j][i - 1].omega - x[j][i + 1].omega) * hydhx;  in FormIFunctionLocal()
373       uyy = (2.0 * u - x[j - 1][i].omega - x[j + 1][i].omega) * hxdhy;  in FormIFunctionLocal()
374 …j][i].omega = (xdot[j][i].omega + uxx + uyy + (vxp * (u - x[j][i - 1].omega) + vxm * (x[j][i + 1].…  in FormIFunctionLocal()
377       u            = x[j][i].temp;  in FormIFunctionLocal()
378       uxx          = (2.0 * u - x[j][i - 1].temp - x[j][i + 1].temp) * hydhx;  in FormIFunctionLocal()
379       uyy          = (2.0 * u - x[j - 1][i].temp - x[j + 1][i].temp) * hxdhy;  in FormIFunctionLocal()
380 …j][i].temp = (xdot[j][i].temp + uxx + uyy + prandtl * ((vxp * (u - x[j][i - 1].temp) + vxm * (x[j]…  in FormIFunctionLocal()