755   PetscInt     ii, i0;  in SetupDS()  local
760   for (ii = ctx->species_offset[grid], i0 = 0; ii < ctx->species_offset[grid + 1]; ii++, i0++) {  in SetupDS()
766     PetscCall(DMSetField(ctx->plex[grid], i0, NULL, (PetscObject)ctx->fe[ii]));  in SetupDS()
770 …for (PetscInt ii = ctx->species_offset[grid], i0 = 0; ii < ctx->species_offset[grid + 1]; ii++, i0…  in SetupDS()  local
773     PetscCall(PetscSectionSetComponentName(section, i0, 0, buf));  in SetupDS()
859 …for (PetscInt ii = ctx->species_offset[grid], i0 = 0; ii < ctx->species_offset[grid + 1]; ii++, i0…  in DMPlexLandauAddMaxwellians()  local
860     mctxs[i0]      = &data[i0];  in DMPlexLandauAddMaxwellians()
861     data[i0].v_0   = ctx->v_0;                             // v_0 same for all grids  in DMPlexLandauAddMaxwellians()
862     data[i0].kT_m  = ctx->k * temps[ii] / ctx->masses[ii]; /* kT/m */  in DMPlexLandauAddMaxwellians()
863     data[i0].n     = ns[ii];  in DMPlexLandauAddMaxwellians()
864     initu[i0]      = maxwellian;  in DMPlexLandauAddMaxwellians()
865     data[i0].shift = 0;  in DMPlexLandauAddMaxwellians()
2226 …for (PetscInt sp = ctx->species_offset[grid], i0 = 0; sp < ctx->species_offset[grid + 1]; sp++, i0…  in DMPlexLandauAccess()  local
2228       PetscInt vf[1] = {i0};  in DMPlexLandauAccess()
2239         PetscCall(func(vdm, vec, i0, grid, b_id, user_ctx));  in DMPlexLandauAccess()
2397   PetscInt    cStart, cEnd, dim, ii, i0, nDMs;  in DMPlexLandauPrintNorms()  local
2417     for (ii = ctx->species_offset[grid], i0 = 0; ii < ctx->species_offset[grid + 1]; ii++, i0++) {  in DMPlexLandauPrintNorms()
2418       PetscScalar user[2] = {(PetscScalar)i0, ctx->charges[ii]};  in DMPlexLandauPrintNorms()