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,<og);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(<og);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