imeshp_test.F90 File Reference

#include "moab/MOABConfig.h"
#include "mpif.h"
#include "iMeshP_f.h"

Include dependency graph for imeshp_test.F90:

Go to the source code of this file.

Macros
#define	ERROR(rval)   if (0 .ne. rval) call exit(1)

Functions/Subroutines
program	imeshp_test
subroutine	create_mesh (imesh, imeshp, part, comm, numv, nume, vgids, nvpe, tp, posn, iconn, vertsPtr, entsPtr, ierr)

Macro Definition Documentation

ERROR

#define ERROR

(

rval

)

if (0 .ne. rval) call exit(1)

Function/Subroutine Documentation

create_mesh()

subroutine create_mesh	(		imesh,
			imeshp,
			part,
		integer	comm,
		integer	numv,
		integer	nume,
		integer, dimension(0:*)	vgids,
		integer	nvpe,
		integer	tp,
		real*8, dimension(0:*)	posn,
		integer, dimension(0:*)	iconn,
		type(c_ptr)	vertsPtr,
		type(c_ptr)	entsPtr,
		integer	ierr
	)

Definition at line 176 of file imeshp_test.F90.

 !
 ! create a mesh with numv vertices and nume elements, with elements of type tp
 ! vertices have positions in posn (3 coordinates each, interleaved xyzxyz...), indexed from 0
 ! elements have nvpe vertices per entity, with connectivity indices stored in iconn, referencing
 ! vertices using 0-based indices; vertex and entity handles are output in arrays passed in
 !
 ! if imesh/imeshp are 0, imesh/imeshp are initialized in this subroutine
 !
 
 use iso_c_binding
 implicit none
 
#ifdef MOAB_HAVE_MPI
# include "iMeshP_f.h"
# include "mpif.h"
#else
# include "iMesh_f.h"
#endif
 
 ! subroutine arguments
 imesh_instance imesh
 TYPE(C_PTR) :: vertsPtr, entsPtr
 integer numv, nume, nvpe, vgids(0:*), iconn(0:*), ierr, tp
 real*8 posn(0:*)
 
 imeshp_partitionhandle imeshp
 integer comm
 
 
 ! local variables
 integer comm_sz, comm_rank, numa, numo, iv, ie
 TYPE(C_PTR) :: statsPtr
 integer, allocatable, target :: stats(:)
 ibase_taghandle tagh
 integer i
 ibase_entityhandle, pointer :: verts(:), ents(:)
 ibase_entityhandle, allocatable :: conn(:)
 ibase_entitysethandle root_set
 ibase_entitysethandle file_set
 
 ibase_handle_t mpi_comm_c
 TYPE(C_PTR) :: partsPtr
 imeshp_parthandle, pointer :: parts(:)
 imeshp_parthandle part
 integer partsa, partso
 character (len=10) filename
 
 ! create the Mesh instance
 if (imesh .eq. 0) then
 call imesh_newmesh("MOAB", imesh, ierr)
 end if
 
 
 if (imeshp .eq. 0) then
 call imeshp_getcommunicator(%VAL(imesh), mpi_comm_world, mpi_comm_c, ierr)
 error(ierr)
 call imeshp_createpartitionall(%VAL(imesh), %VAL(mpi_comm_c), imeshp, ierr)
 error(ierr)
 call imeshp_createpart(%VAL(imesh), %VAL(imeshp), part, ierr)
 error(ierr)
 else
 partsa = 0
 call imeshp_getlocalparts(%VAL(imesh), %VAL(imeshp), partsptr, partsa, partso, ierr)
 error(ierr)
 call c_f_pointer(partsptr, parts, [partso])
 part = parts(1)
 end if
 call mpi_comm_rank(comm, comm_rank, ierr)
 error(ierr)
 call mpi_comm_size(comm, comm_sz, ierr)
 error(ierr)
 
 ! create the vertices, all in one call
 numa = 0
 call imesh_createvtxarr(%VAL(imesh), %VAL(numv), %VAL(ibase_interleaved), posn, %VAL(3*numv), &
 vertsptr, numa, numo, ierr)
 error(ierr)
 
 ! fill in the connectivity array, based on indexing from iconn
 allocate (conn(0:nvpe*nume-1))
 call c_f_pointer(vertsptr, verts, [numv])
 do i = 0, nvpe*nume-1
 conn(i) = verts(1+iconn(i))
 end do
 ! create the elements
 numa = 0
 allocate(stats(0:nume-1))
 statsptr = c_loc(stats(0))
 call imesh_createentarr(%VAL(imesh), %VAL(tp), conn, %VAL(nvpe*nume), &
 entsptr, numa, numo, statsptr, numa, numo, ierr)
 deallocate(stats)
 deallocate(conn)
 
 ! take care of parallel stuff
 
 ! add entities to part, using iMesh
 call c_f_pointer(entsptr, ents, [numo])
 call imesh_addentarrtoset(%VAL(imesh), ents, %VAL(numo), %VAL(part), ierr)
 error(ierr)
 ! set global ids on vertices, needed for sharing between procs
 call imesh_gettaghandle(%VAL(imesh), "GLOBAL_ID", tagh, ierr, %VAL(9))
 if (ibase_success .ne. ierr) then
 ! didn't get handle, need to create the tag
 call imesh_createtag(%VAL(imesh), "GLOBAL_ID", %VAL(ibase_integer), tagh, ierr)
 error(ierr)
 end if
 call imesh_setintarrdata(%VAL(imesh), verts, %VAL(numv), %VAL(tagh), vgids, %VAL(numv), ierr)
 error(ierr)
 ! now resolve shared verts and exchange ghost cells
 call imeshp_syncmeshall(%VAL(imesh), %VAL(imeshp), ierr)
 error(ierr)
 
 call imeshp_createghostentsall(%VAL(imesh), %VAL(imeshp), %VAL(2), %VAL(1), %VAL(1), %VAL(0), ierr)
 error(ierr)
 
 call imesh_freememory(%VAL(imesh), vertsptr);
 call imesh_freememory(%VAL(imesh), entsptr);
 
 return

References error().

imeshp_test()

program imeshp_test

Definition at line 11 of file imeshp_test.F90.

References create_mesh(), and error().