2052   Vec        Xsub[LANDAU_MAX_GRIDS];  in DMPlexLandauCreateVelocitySpace()  local
2072     PetscCall(DMCreateGlobalVector(ctx->plex[grid], &Xsub[grid]));  in DMPlexLandauCreateVelocitySpace()
2073     PetscCall(PetscObjectSetName((PetscObject)Xsub[grid], "u_orig"));  in DMPlexLandauCreateVelocitySpace()
2075     PetscCall(LandauSetInitialCondition(ctx->plex[grid], Xsub[grid], grid, 0, 1, ctx));  in DMPlexLandauCreateVelocitySpace()
2079       PetscCall(adapt(grid, ctx, &Xsub[grid])); // forest goes in, plex comes out  in DMPlexLandauCreateVelocitySpace()
2086       PetscCall(LandauSetInitialCondition(ctx->plex[grid], Xsub[grid], grid, 0, 1, ctx));  in DMPlexLandauCreateVelocitySpace()
2090       PetscCall(VecSetOptionsPrefix(Xsub[grid], prefix));  in DMPlexLandauCreateVelocitySpace()
2091       PetscCall(VecViewFromOptions(Xsub[grid], NULL, "-dm_landau_amr_vec_view"));  in DMPlexLandauCreateVelocitySpace()
2112     PetscCall(VecGetLocalSize(Xsub[grid], &n));  in DMPlexLandauCreateVelocitySpace()
2127     PetscCall(VecGetLocalSize(Xsub[grid], &n));  in DMPlexLandauCreateVelocitySpace()
2131 … PetscCall(LandauSetInitialCondition(ctx->plex[grid], Xsub[grid], grid, b_id, ctx->batch_sz, ctx));  in DMPlexLandauCreateVelocitySpace()
2132       PetscCall(VecGetArrayRead(Xsub[grid], &values)); // Drop whole grid in Plex ordering  in DMPlexLandauCreateVelocitySpace()
2134       PetscCall(VecRestoreArrayRead(Xsub[grid], &values));  in DMPlexLandauCreateVelocitySpace()
2138   for (PetscInt grid = 0; grid < ctx->num_grids; grid++) PetscCall(VecDestroy(&Xsub[grid]));  in DMPlexLandauCreateVelocitySpace()