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,&ltog);

 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(&ltog);
 27:   ISDestroy(&isx);

 29:   PetscFinalize();
 30:   return 0;
 31: }

 33: /*TEST

 35:    test:

 37: TEST*/