static char help[]= "Tests ISLocalToGlobalMappingCreateIS() for bs > 1.\n\n";

#include <petscis.h>
#include <petscviewer.h>

int main(int argc,char **argv)
{
  PetscInt               bs = 2,n = 3,ix[3] = {1,7,9},iy[2] = {0,2},mp[2];
  IS                     isx;
  ISLocalToGlobalMapping ltog;

  PetscCall(PetscInitialize(&argc,&argv,(char*)0,help));

  PetscCall(ISCreateBlock(PETSC_COMM_SELF,bs,n,ix,PETSC_COPY_VALUES,&isx));
  PetscCall(ISLocalToGlobalMappingCreateIS(isx,&ltog));

  PetscCall(PetscIntView(2,iy,PETSC_VIEWER_STDOUT_WORLD));
  PetscCall(ISLocalToGlobalMappingApply(ltog,2,iy,mp));
  PetscCall(PetscIntView(2,mp,PETSC_VIEWER_STDOUT_WORLD));

  PetscCall(PetscIntView(2,iy,PETSC_VIEWER_STDOUT_WORLD));
  PetscCall(ISLocalToGlobalMappingApplyBlock(ltog,2,iy,mp));
  PetscCall(PetscIntView(2,mp,PETSC_VIEWER_STDOUT_WORLD));

  PetscCall(ISLocalToGlobalMappingDestroy(&ltog));
  PetscCall(ISDestroy(&isx));

  PetscCall(PetscFinalize());
  return 0;
}

/*TEST

   test:

TEST*/