137   PetscScalar   hx, hy, hxdhy, hydhx;  in FormFunction()  local
151   hydhx  = hy / hx;  in FormFunction()
282       f[j][i] = -hydhx * (fe - fw) - hxdhy * (fn - fs);  in FormFunction()
296   PetscScalar one = 1.0, hx, hy, hxdhy, hydhx, t0, tn, ts, te, tw;  in FormJacobian()  local
311   hydhx  = hy / hx;  in FormJacobian()
363         v[1]     = -hydhx * (dw - gw);  in FormJacobian()
366         v[2]     = hxdhy * (ds + dn + gs - gn) + hydhx * (dw + de + gw - ge);  in FormJacobian()
369         v[3]             = -hydhx * (de + ge);  in FormJacobian()
402           v[0]     = hxdhy * (dn - gn) + hydhx * (dw + de + gw - ge);  in FormJacobian()
405           v[1]             = -hydhx * (de + ge);  in FormJacobian()
432           v[1]     = hxdhy * (ds + dn + gs - gn) + hydhx * (dw + de + gw - ge);  in FormJacobian()
435           v[2]             = -hydhx * (de + ge);  in FormJacobian()
454           v[1]     = hxdhy * (ds + gs) + hydhx * (dw + de + gw - ge);  in FormJacobian()
457           v[2]             = -hydhx * (de + ge);  in FormJacobian()
488           v[0]     = -hydhx * (dw - gw);  in FormJacobian()
491           v[1]     = hxdhy * (dn - gn) + hydhx * (dw + de + gw - ge);  in FormJacobian()
518           v[1]     = -hydhx * (dw - gw);  in FormJacobian()
521           v[2]     = hxdhy * (ds + dn + gs - gn) + hydhx * (dw + de + gw - ge);  in FormJacobian()
540           v[1]     = -hydhx * (dw - gw);  in FormJacobian()
543           v[2]     = hxdhy * (ds + gs) + hydhx * (dw + de + gw - ge);  in FormJacobian()
572         v[0]     = -hydhx * (dw - gw);  in FormJacobian()
575         v[1]     = hxdhy * (dn - gn) + hydhx * (dw + de + gw - ge);  in FormJacobian()
578         v[2]             = -hydhx * (de + ge);  in FormJacobian()
612         v[1]     = -hydhx * (dw - gw);  in FormJacobian()
615         v[2]     = hxdhy * (ds + gs) + hydhx * (dw + de + gw - ge);  in FormJacobian()
618         v[3]             = -hydhx * (de + ge);  in FormJacobian()