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