77   PetscScalar u, v, omega, temp;  member
221       x[j][i].omega = 0.0;  in FormInitialGuess()
273       f[j][i].omega = x[j][i].omega + (x[j + 1][i].u - x[j][i].u) * dhy;  in FormFunctionLocal()
286       f[j][i].omega = x[j][i].omega + (x[j][i].u - x[j - 1][i].u) * dhy;  in FormFunctionLocal()
299       f[j][i].omega = x[j][i].omega - (x[j][i + 1].v - x[j][i].v) * dhx;  in FormFunctionLocal()
312       f[j][i].omega = x[j][i].omega - (x[j][i].v - x[j][i - 1].v) * dhx;  in FormFunctionLocal()
336       f[j][i].u = uxx + uyy - .5 * (x[j + 1][i].omega - x[j - 1][i].omega) * hx;  in FormFunctionLocal()
342       f[j][i].v = uxx + uyy + .5 * (x[j][i + 1].omega - x[j][i - 1].omega) * hy;  in FormFunctionLocal()
345       u             = x[j][i].omega;  in FormFunctionLocal()
346       uxx           = (2.0 * u - x[j][i - 1].omega - x[j][i + 1].omega) * hydhx;  in FormFunctionLocal()
347       uyy           = (2.0 * u - x[j - 1][i].omega - x[j + 1][i].omega) * hxdhy;  in FormFunctionLocal()
348 …].omega = uxx + uyy + (vxp * (u - x[j][i - 1].omega) + vxm * (x[j][i + 1].omega - u)) * hy + (vyp …  in FormFunctionLocal()
511             bjiomega = b[j][i].omega;  in NonlinearGS()
525             fu    = uxx + uyy - .5 * (x[j + 1][i].omega - x[j - 1][i].omega) * hx - bjiu;  in NonlinearGS()
531             fv    = uxx + uyy + .5 * (x[j][i + 1].omega - x[j][i - 1].omega) * hy - bjiv;  in NonlinearGS()
545             u      = x[j][i].omega;  in NonlinearGS()
546             uxx    = (2.0 * u - x[j][i - 1].omega - x[j][i + 1].omega) * hydhx;  in NonlinearGS()
547             uyy    = (2.0 * u - x[j - 1][i].omega - x[j + 1][i].omega) * hxdhy;  in NonlinearGS()
548 …(vxp * (u - x[j][i - 1].omega) + vxm * (x[j][i + 1].omega - u)) * hy + (vyp * (u - x[j - 1][i].ome…  in NonlinearGS()
551             if (PetscRealPart(vx) > 0.0) dfodu = (u - x[j][i - 1].omega) * hy;  in NonlinearGS()
552             else dfodu = (x[j][i + 1].omega - u) * hy;  in NonlinearGS()
554             if (PetscRealPart(vy) > 0.0) dfodv = (u - x[j - 1][i].omega) * hx;  in NonlinearGS()
555             else dfodv = (x[j + 1][i].omega - u) * hx;  in NonlinearGS()
584             x[j][i].omega = x[j][i].omega - yo;  in NonlinearGS()
587             fomega        = x[j][i].omega - (x[j][i + 1].v - x[j][i].v) * dhx - bjiomega;  in NonlinearGS()
591             x[j][i].omega = x[j][i].omega - fomega;  in NonlinearGS()
595             fomega        = x[j][i].omega - (x[j][i].v - x[j][i - 1].v) * dhx - bjiomega;  in NonlinearGS()
599             x[j][i].omega = x[j][i].omega - fomega;  in NonlinearGS()
603             fomega        = x[j][i].omega + (x[j + 1][i].u - x[j][i].u) * dhy - bjiomega;  in NonlinearGS()
607             x[j][i].omega = x[j][i].omega - fomega;  in NonlinearGS()
611             fomega        = x[j][i].omega + (x[j][i].u - x[j - 1][i].u) * dhy - bjiomega;  in NonlinearGS()
615             x[j][i].omega = x[j][i].omega - fomega;  in NonlinearGS()
623 …cRealPart(x[j][i].u * x[j][i].u + x[j][i].v * x[j][i].v + x[j][i].omega * x[j][i].omega + x[j][i].…  in NonlinearGS()