unconstrained/tutorials/burgers_spectral.c

100   PetscInt     i, xs, xm, ind, j, lenglob;
161   PetscCall(DMDAGetCorners(appctx.da, &xs, NULL, NULL, &xm, NULL, NULL));
167   xs = xs / (appctx.param.N - 1);
174   for (i = xs; i < xs + xm; i++) {
298   PetscInt           i, xs, xn;
303   PetscCall(DMDAGetCorners(appctx->da, &xs, NULL, NULL, &xn, NULL, NULL));
304   for (i = xs; i < xs + xn; i++) s[i] = 2.0 * appctx->param.mu * PETSC_PI * PetscSinScalar(PETSC_PI * xg[i]) / (2.0 + PetscCosScalar(PETSC_PI * xg[i])) + 0.25 * PetscExpReal(-4.0 * PetscPowReal(xg[i] - 2.0, 2.0));
326   PetscInt           i, xs, xn;
331   PetscCall(DMDAGetCorners(appctx->da, &xs, NULL, NULL, &xn, NULL, NULL));
332   for (i = xs; i < xs + xn; i++) s[i] = 2.0 * appctx->param.mu * PETSC_PI * PetscSinScalar(PETSC_PI * xg[i]) / (2.0 + PetscCosScalar(PETSC_PI * xg[i]));
351   PetscInt           i, xs, xn;
356   PetscCall(DMDAGetCorners(appctx->da, &xs, NULL, NULL, &xn, NULL, NULL));
357   for (i = xs; i < xs + xn; i++) {
430   PetscInt    i, xs, xn, l, j;
448   PetscCall(DMDAGetCorners(appctx->da, &xs, NULL, NULL, &xn, NULL, NULL));
450   xs = xs / (appctx->param.N - 1);
457   for (j = xs; j < xs + xn; j++) {
458     for (l = 0; l < appctx->param.N; l++) rowsDM[l] = 1 + (j - xs) * (appctx->param.N - 1) + l;
491   PetscInt    xs, xn, l, j;
501   PetscCall(DMDAGetCorners(appctx->da, &xs, NULL, NULL, &xn, NULL, NULL));
503   xs = xs / (appctx->param.N - 1);
507   for (j = xs; j < xs + xn; j++) {
508     for (l = 0; l < appctx->param.N; l++) rowsDM[l] = 1 + (j - xs) * (appctx->param.N - 1) + l;