Lines Matching refs:j
162 PetscInt i, j, mx, xs, ys, xm, ym; in FormInitialGuess() local
192 for (j = ys; j < ys + ym; j++) { in FormInitialGuess()
194 x[j][i].u = 0.0; in FormInitialGuess()
195 x[j][i].v = 0.0; in FormInitialGuess()
196 x[j][i].omega = 0.0; in FormInitialGuess()
197 x[j][i].temp = (grashof > 0) * i * dx; in FormInitialGuess()
211 PetscInt xints, xinte, yints, yinte, i, j; in FormFunctionLocal() local
248 j = 0; in FormFunctionLocal()
252 f[j][i].u = x[j][i].u; in FormFunctionLocal()
253 f[j][i].v = x[j][i].v; in FormFunctionLocal()
254 f[j][i].omega = x[j][i].omega + (x[j + 1][i].u - x[j][i].u) * dhy; in FormFunctionLocal()
255 f[j][i].temp = x[j][i].temp - x[j + 1][i].temp; in FormFunctionLocal()
261 j = info->my - 1; in FormFunctionLocal()
265 f[j][i].u = x[j][i].u - lid; in FormFunctionLocal()
266 f[j][i].v = x[j][i].v; in FormFunctionLocal()
267 f[j][i].omega = x[j][i].omega + (x[j][i].u - x[j - 1][i].u) * dhy; in FormFunctionLocal()
268 f[j][i].temp = x[j][i].temp - x[j - 1][i].temp; in FormFunctionLocal()
277 for (j = info->ys; j < info->ys + info->ym; j++) { in FormFunctionLocal()
278 f[j][i].u = x[j][i].u; in FormFunctionLocal()
279 f[j][i].v = x[j][i].v; in FormFunctionLocal()
280 f[j][i].omega = x[j][i].omega - (x[j][i + 1].v - x[j][i].v) * dhx; in FormFunctionLocal()
281 f[j][i].temp = x[j][i].temp; in FormFunctionLocal()
290 for (j = info->ys; j < info->ys + info->ym; j++) { in FormFunctionLocal()
291 f[j][i].u = x[j][i].u; in FormFunctionLocal()
292 f[j][i].v = x[j][i].v; in FormFunctionLocal()
293 f[j][i].omega = x[j][i].omega - (x[j][i].v - x[j][i - 1].v) * dhx; in FormFunctionLocal()
294 f[j][i].temp = x[j][i].temp - (PetscReal)(grashof > 0); in FormFunctionLocal()
299 for (j = yints; j < yinte; j++) { in FormFunctionLocal()
304 vx = x[j][i].u; in FormFunctionLocal()
308 vy = x[j][i].v; in FormFunctionLocal()
314 u = x[j][i].u; in FormFunctionLocal()
315 uxx = (2.0 * u - x[j][i - 1].u - x[j][i + 1].u) * hydhx; in FormFunctionLocal()
316 uyy = (2.0 * u - x[j - 1][i].u - x[j + 1][i].u) * hxdhy; in FormFunctionLocal()
317 f[j][i].u = uxx + uyy - .5 * (x[j + 1][i].omega - x[j - 1][i].omega) * hx; in FormFunctionLocal()
320 u = x[j][i].v; in FormFunctionLocal()
321 uxx = (2.0 * u - x[j][i - 1].v - x[j][i + 1].v) * hydhx; in FormFunctionLocal()
322 uyy = (2.0 * u - x[j - 1][i].v - x[j + 1][i].v) * hxdhy; in FormFunctionLocal()
323 f[j][i].v = uxx + uyy + .5 * (x[j][i + 1].omega - x[j][i - 1].omega) * hy; in FormFunctionLocal()
326 u = x[j][i].omega; in FormFunctionLocal()
327 uxx = (2.0 * u - x[j][i - 1].omega - x[j][i + 1].omega) * hydhx; in FormFunctionLocal()
328 uyy = (2.0 * u - x[j - 1][i].omega - x[j + 1][i].omega) * hxdhy; in FormFunctionLocal()
329 …j][i].omega = uxx + uyy + (vxp * (u - x[j][i - 1].omega) + vxm * (x[j][i + 1].omega - u)) * hy + (… in FormFunctionLocal()
332 u = x[j][i].temp; in FormFunctionLocal()
333 uxx = (2.0 * u - x[j][i - 1].temp - x[j][i + 1].temp) * hydhx; in FormFunctionLocal()
334 uyy = (2.0 * u - x[j - 1][i].temp - x[j + 1][i].temp) * hxdhy; in FormFunctionLocal()
335 …f[j][i].temp = uxx + uyy + prandtl * ((vxp * (u - x[j][i - 1].temp) + vxm * (x[j][i + 1].temp - u)… in FormFunctionLocal()