Actual source code: ex1.c
2: static char help[]= "Tests ISLocalToGlobalMappingCreateIS() for bs > 1.\n\n";
4: #include <petscis.h>
5: #include <petscviewer.h>
7: int main(int argc,char **argv)
8: {
9: PetscInt bs = 2,n = 3,ix[3] = {1,7,9},iy[2] = {0,2},mp[2];
10: IS isx;
11: ISLocalToGlobalMapping ltog;
13: PetscInitialize(&argc,&argv,(char*)0,help);
15: ISCreateBlock(PETSC_COMM_SELF,bs,n,ix,PETSC_COPY_VALUES,&isx);
16: ISLocalToGlobalMappingCreateIS(isx,<og);
18: PetscIntView(2,iy,PETSC_VIEWER_STDOUT_WORLD);
19: ISLocalToGlobalMappingApply(ltog,2,iy,mp);
20: PetscIntView(2,mp,PETSC_VIEWER_STDOUT_WORLD);
22: PetscIntView(2,iy,PETSC_VIEWER_STDOUT_WORLD);
23: ISLocalToGlobalMappingApplyBlock(ltog,2,iy,mp);
24: PetscIntView(2,mp,PETSC_VIEWER_STDOUT_WORLD);
26: ISLocalToGlobalMappingDestroy(<og);
27: ISDestroy(&isx);
29: PetscFinalize();
30: return 0;
31: }
33: /*TEST
35: test:
37: TEST*/