1 
2 static char help[]= "Tests ISLocalToGlobalMappingCreateIS() for bs > 1.\n\n";
3 
4 #include <petscis.h>
5 #include <petscviewer.h>
6 
main(int argc,char ** argv)7 int main(int argc,char **argv)
8 {
9   PetscErrorCode         ierr;
10   PetscInt               bs = 2,n = 3,ix[3] = {1,7,9},iy[2] = {0,2},mp[2];
11   IS                     isx;
12   ISLocalToGlobalMapping ltog;
13 
14   ierr = PetscInitialize(&argc,&argv,(char*)0,help);if (ierr) return ierr;
15 
16   ierr = ISCreateBlock(PETSC_COMM_SELF,bs,n,ix,PETSC_COPY_VALUES,&isx);CHKERRQ(ierr);
17   ierr = ISLocalToGlobalMappingCreateIS(isx,&ltog);CHKERRQ(ierr);
18 
19   ierr = PetscIntView(2,iy,PETSC_VIEWER_STDOUT_WORLD);CHKERRQ(ierr);
20   ierr = ISLocalToGlobalMappingApply(ltog,2,iy,mp);CHKERRQ(ierr);
21   ierr = PetscIntView(2,mp,PETSC_VIEWER_STDOUT_WORLD);CHKERRQ(ierr);
22 
23   ierr = PetscIntView(2,iy,PETSC_VIEWER_STDOUT_WORLD);CHKERRQ(ierr);
24   ierr = ISLocalToGlobalMappingApplyBlock(ltog,2,iy,mp);CHKERRQ(ierr);
25   ierr = PetscIntView(2,mp,PETSC_VIEWER_STDOUT_WORLD);CHKERRQ(ierr);
26 
27   ierr = ISLocalToGlobalMappingDestroy(&ltog);CHKERRQ(ierr);
28   ierr = ISDestroy(&isx);CHKERRQ(ierr);
29 
30   ierr = PetscFinalize();
31   return ierr;
32 }
33 
34 
35 
36 /*TEST
37 
38    test:
39 
40 TEST*/
41