Actual source code: ex8.c


  2: static char help[] = "Tests ISLocalToGlobalMappingSetBlockSize.\n\n";

  4: #include <petscis.h>
  5: #include <petscviewer.h>

  7: int main(int argc, char **argv)
  8: {
  9:   PetscInt               indices[]  = {0, 1, 2, 3, -1, -1, -1, -1, 4, 5, 6, 7};
 10:   PetscInt               indices2[] = {0, 1, 2, 3, 4, 5, -1, -1, -1, -1, -1, -1, 6, 7, 8, 9, 10, 11};
 11:   ISLocalToGlobalMapping map;

 13:   PetscFunctionBeginUser;
 14:   PetscCall(PetscInitialize(&argc, &argv, (char *)0, help));
 15:   PetscCall(ISLocalToGlobalMappingCreate(PETSC_COMM_WORLD, 1, 12, indices, PETSC_COPY_VALUES, &map));
 16:   PetscCall(ISLocalToGlobalMappingView(map, NULL));
 17:   PetscCall(ISLocalToGlobalMappingSetBlockSize(map, 2));
 18:   PetscCall(ISLocalToGlobalMappingView(map, NULL));
 19:   PetscCall(ISLocalToGlobalMappingSetBlockSize(map, 4));
 20:   PetscCall(ISLocalToGlobalMappingView(map, NULL));
 21:   PetscCall(ISLocalToGlobalMappingSetBlockSize(map, 2));
 22:   PetscCall(ISLocalToGlobalMappingView(map, NULL));
 23:   PetscCall(ISLocalToGlobalMappingSetBlockSize(map, 1));
 24:   PetscCall(ISLocalToGlobalMappingView(map, NULL));
 25:   PetscCall(ISLocalToGlobalMappingDestroy(&map));
 26:   PetscCall(ISLocalToGlobalMappingCreate(PETSC_COMM_WORLD, 1, 18, indices2, PETSC_COPY_VALUES, &map));
 27:   PetscCall(ISLocalToGlobalMappingView(map, NULL));
 28:   PetscCall(ISLocalToGlobalMappingSetBlockSize(map, 3));
 29:   PetscCall(ISLocalToGlobalMappingView(map, NULL));
 30:   PetscCall(ISLocalToGlobalMappingSetBlockSize(map, 6));
 31:   PetscCall(ISLocalToGlobalMappingView(map, NULL));
 32:   PetscCall(ISLocalToGlobalMappingSetBlockSize(map, 3));
 33:   PetscCall(ISLocalToGlobalMappingView(map, NULL));
 34:   PetscCall(ISLocalToGlobalMappingSetBlockSize(map, 1));
 35:   PetscCall(ISLocalToGlobalMappingView(map, NULL));
 36:   PetscCall(ISLocalToGlobalMappingDestroy(&map));
 37:   PetscCall(ISLocalToGlobalMappingCreate(PETSC_COMM_WORLD, 5, 2, indices2, PETSC_COPY_VALUES, &map));
 38:   PetscCall(ISLocalToGlobalMappingView(map, NULL));
 39:   PetscCall(ISLocalToGlobalMappingSetBlockSize(map, 2));
 40:   PetscCall(ISLocalToGlobalMappingView(map, NULL));
 41:   PetscCall(ISLocalToGlobalMappingDestroy(&map));
 42:   PetscCall(PetscFinalize());
 43:   return 0;
 44: }

 46: /*TEST

 48:    test:

 50: TEST*/