FindConnectF.F

Go to the documentation of this file.

c FindConnect: Interacting with iMesh
c 
c This program shows how to get more information about a mesh, by
c getting connectivity two different ways (as connectivity and as
c adjacent 0-dimensional entities).
 
c Usage: FindConnect
 
      program findconnect
#include "iMesh_f.h"
 
c declarations
      imesh_instance mesh
      ibase_entityhandle  ents
      ibase_handle_t rpverts, rpallverts, ipoffsets
      pointer(rpents, ents(0:*))
      pointer(rpverts, verts(0:*))
      pointer(rpallverts, allverts(0:*))
      pointer(ipoffsets, ioffsets(0,*))
      integer ierr, ents_alloc, ents_size
      integer iverts_alloc, iverts_size
      integer allverts_alloc, allverts_size
      integer offsets_alloc, offsets_size
 
c create the Mesh instance
      call imesh_newmesh("MOAB", mesh, ierr)
 
c load the mesh
      call imesh_load(%VAL(mesh), %VAL(0), "125hex.vtk", "", ierr)
 
c get all 3d elements
      ents_alloc = 0
      call imesh_getentities(%VAL(mesh), %VAL(0), %VAL(ibase_region), 
     1     %VAL(imesh_all_topologies), rpents, ents_alloc, ents_size,
     1     ierr)
 
      ivert_uses = 0
 
c iterate through them; 
      do i = 0, ents_size-1
c get connectivity
         iverts_alloc = 0
         call imesh_getentadj(%VAL(mesh), %VAL(ents(i)), 
     1        %VAL(ibase_vertex), rpverts, iverts_alloc, iverts_size,
     1        ierr)
c sum number of vertex uses
 
         vert_uses = vert_uses + iverts_size
 
         call free(rpverts)
      end do
 
c now get adjacencies in one big block
      allverts_alloc = 0
      offsets_alloc = 0
      call imesh_getentarradj(%VAL(mesh), %VAL(rpents), 
     1     %VAL(ents_size), %VAL(ibase_vertex), rpallverts, 
     1     allverts_alloc, allverts_size, ipoffsets, offsets_alloc, 
     1     offsets_size, ierr)
 
      call free(rpallverts);
      call free(ipoffsets);
      call free(rpents);
 
c compare results of two calling methods
      if (allverts_size .ne. vert_uses) then
         write(*,'("Sizes didn''t agree!")')
      else 
         write(*,'("Sizes did agree!")')
      endif
 
      call imesh_dtor(%VAL(mesh), ierr)
 
      end