program main
#include <petsc/finclude/petscdmplex.h>
      use petscdmplex
      use petscsys
      implicit none
!
!
      DM dm
      PetscInt, target, dimension(4) :: EC
      PetscInt, pointer :: pEC(:)
      PetscInt, pointer :: pES(:)
      PetscInt c, firstCell, numCells
      PetscInt v, numVertices, numPoints
      PetscInt i0,i4
      PetscErrorCode ierr

      i0 = 0
      i4 = 4

      PetscCallA(PetscInitialize(ierr))

      PetscCallA(DMPlexCreate(PETSC_COMM_WORLD, dm, ierr))
      firstCell = 0
      numCells = 2
      numVertices = 6
      numPoints = numCells+numVertices
      PetscCallA(DMPlexSetChart(dm, i0, numPoints, ierr))
      do c=firstCell,numCells-1
         PetscCallA(DMPlexSetConeSize(dm, c, i4, ierr))
      end do
      PetscCallA(DMSetUp(dm, ierr))

      EC(1) = 2
      EC(2) = 3
      EC(3) = 4
      EC(4) = 5
      pEC => EC
      c = 0
      write(*,1000) 'cell',c,pEC
 1000 format (a,i4,50i4)
      PetscCallA(DMPlexSetCone(dm, c , pEC, ierr))
      PetscCallA(DMPlexGetCone(dm, c , pEC, ierr))
      write(*,1000) 'cell',c,pEC
      EC(1) = 4
      EC(2) = 5
      EC(3) = 6
      EC(4) = 7
      pEC => EC
      c = 1
      write(*,1000) 'cell',c,pEC
      PetscCallA(DMPlexSetCone(dm, c , pEC, ierr))
      PetscCallA(DMPlexGetCone(dm, c , pEC, ierr))
      write(*,1000) 'cell',c,pEC
      PetscCallA(DMPlexRestoreCone(dm, c , pEC, ierr))

      PetscCallA(DMPlexSymmetrize(dm, ierr))
      PetscCallA(DMPlexStratify(dm, ierr))

      v = 4
      PetscCallA(DMPlexGetSupport(dm, v , pES, ierr))
      write(*,1000) 'vertex',v,pES
      PetscCallA(DMPlexRestoreSupport(dm, v , pES, ierr))

      PetscCallA(DMDestroy(dm,ierr))
      PetscCallA(PetscFinalize(ierr))
      end

! /*TEST
!
! test:
!   suffix: 0
!
! TEST*/