Lines Matching refs:m
21 PetscInt m; /* the number of mesh points in x-direction */ member
60 data.m = 9; in main()
64 data.dx = 1.0 / (data.m + 1.0); in main()
66 mn = (data.m) * (data.n); in main()
198 PetscInt m, row, col; in Initial() local
205 m = data->m; in Initial()
217 row = 1 + (mybase + i) - ((mybase + i) / m) * m; in Initial()
218 col = (mybase + i) / m + 1; in Initial()
276 PetscInt m, n, mn; in RHSJacobian() local
280 m = data->m; in RHSJacobian()
282 mn = m * n; in RHSJacobian()
300 idx[2] = m; in RHSJacobian()
303 row = m - 1; in RHSJacobian()
307 idx[0] = m - 2; in RHSJacobian()
308 idx[1] = m - 1; in RHSJacobian()
309 idx[2] = m - 1 + m; in RHSJacobian()
312 for (i = 1; i < m - 1; i++) { in RHSJacobian()
321 idx[3] = i + m; in RHSJacobian()
326 row = j * m; in RHSJacobian()
331 idx[0] = j * m; in RHSJacobian()
332 idx[1] = j * m; in RHSJacobian()
333 idx[2] = j * m - m; in RHSJacobian()
334 idx[3] = j * m + m; in RHSJacobian()
337 row = j * m + m - 1; in RHSJacobian()
342 idx[0] = j * m + m - 1; in RHSJacobian()
343 idx[1] = j * m + m - 1 - 1; in RHSJacobian()
344 idx[2] = j * m + m - 1 - m; in RHSJacobian()
345 idx[3] = j * m + m - 1 + m; in RHSJacobian()
348 for (i = 1; i < m - 1; i++) { in RHSJacobian()
349 row = j * m + i; in RHSJacobian()
355 idx[0] = j * m + i; in RHSJacobian()
356 idx[1] = j * m + i - 1; in RHSJacobian()
357 idx[2] = j * m + i + 1; in RHSJacobian()
358 idx[3] = j * m + i - m; in RHSJacobian()
359 idx[4] = j * m + i + m; in RHSJacobian()
364 row = mn - m; in RHSJacobian()
368 idx[0] = mn - m; in RHSJacobian()
369 idx[1] = mn - m + 1; in RHSJacobian()
370 idx[2] = mn - m - m; in RHSJacobian()
379 idx[2] = mn - 1 - m; in RHSJacobian()
382 for (i = 1; i < m - 1; i++) { in RHSJacobian()
383 row = mn - m + i; in RHSJacobian()
388 idx[0] = mn - m + i - 1; in RHSJacobian()
389 idx[1] = mn - m + i; in RHSJacobian()
390 idx[2] = mn - m + i + 1; in RHSJacobian()
391 idx[3] = mn - m + i - m; in RHSJacobian()
404 PetscInt m, n, mn; in RHSFunction() local
417 m = data->m; in RHSFunction()
419 mn = m * n; in RHSFunction()
456 outptr[0] = xc * inptr[0] + xr * inptr[1] + yr * inptr[m]; in RHSFunction()
457 outptr[m - 1] = 2.0 * xl * inptr[m - 2] + xc * inptr[m - 1] + yr * inptr[m - 1 + m]; in RHSFunction()
458 …for (i = 1; i < m - 1; i++) outptr[i] = xc * inptr[i] + xl * inptr[i - 1] + xr * inptr[i + 1] + yr… in RHSFunction()
461 …outptr[j * m] = xc * inptr[j * m] + xr * inptr[j * m + 1] + yl * inptr[j * m - m] + yr * i… in RHSFunction()
462 …tr[j * m + m - 1] = xc * inptr[j * m + m - 1] + 2.0 * xl * inptr[j * m + m - 1 - 1] + yl * inptr[j… in RHSFunction()
463 …m - 1; i++) outptr[j * m + i] = xc * inptr[j * m + i] + xl * inptr[j * m + i - 1] + xr * inptr[j *… in RHSFunction()
466 outptr[mn - m] = xc * inptr[mn - m] + xr * inptr[mn - m + 1] + 2.0 * yl * inptr[mn - m - m]; in RHSFunction()
467 outptr[mn - 1] = 2.0 * xl * inptr[mn - 2] + xc * inptr[mn - 1] + 2.0 * yl * inptr[mn - 1 - m]; in RHSFunction()
468 …i < m - 1; i++) outptr[mn - m + i] = xc * inptr[mn - m + i] + xl * inptr[mn - m + i - 1] + xr * in… in RHSFunction()