DirectAccessNoHolesF90.F90 File Reference

Go to the source code of this file.

Functions/Subroutines
program	directaccessnoholesf90
subroutine	create_mesh_no_holes (imesh, nquads, ierr)

Function/Subroutine Documentation

create_mesh_no_holes()

subroutine create_mesh_no_holes	(		imesh,
		integer	nquads,
		integer	ierr
	)

Definition at line 196 of file DirectAccessNoHolesF90.F90.

 use, intrinsic :: iso_c_binding
 implicit none
#include "iMesh_f.h"
 
 imesh_instance imesh
 integer nquads, ierr
 
 real(C_DOUBLE), pointer :: coords(:,:)
 TYPE(C_PTR) tmp_ents_ptr, stat_ptr
 ibase_entityhandle, pointer :: connect(:), tmp_ents(:)
 integer nverts, tmp_size, stat_size, i
 
 ! first make the mesh, a 1d array of quads with left hand side x = elem_num;
 ! vertices are numbered in layers
 ! create verts, num is 2(nquads+1) because they're in a 1d row
 nverts = 2*(nquads+1)
 allocate(coords(0:2, 0:nverts-1))
 do i = 0, nquads-1
 coords(0,2*i) = i
 coords(0,2*i+1) = i ! x values are all i
 coords(1,2*i) = 0.0
 coords(1,2*i+1) = 1.0 ! y coords
 coords(2,2*i) = 0.0
 coords(2,2*i+1) = 0.0 ! z values, all zero (2d mesh)
 end do
 ! last two vertices
 coords(0, nverts-2) = nquads
 coords(0, nverts-1) = nquads
 coords(1, nverts-2) = 0.0
 coords(1, nverts-1) = 1.0 ! y coords
 coords(2, nverts-2) = 0.0
 coords(2, nverts-1) = 0.0 ! z values, all zero (2d mesh)
 tmp_size = 0
 !!!!!! VSM: This is the culprit call that screws up IBM xlf compiler
 call imesh_createvtxarr(%val(imesh), %val(nverts), %val(ibase_interleaved), &
 coords, %val(3*nverts), tmp_ents_ptr, tmp_size, tmp_size, ierr)
 check(ierr)
 call c_f_pointer(tmp_ents_ptr, tmp_ents, [nverts])
 deallocate(coords)
 allocate(connect(0:4*nquads-1))
 do i = 0, nquads-1
 connect(4*i+0) = tmp_ents(1+2*i)
 connect(4*i+1) = tmp_ents(1+2*i+2)
 connect(4*i+2) = tmp_ents(1+2*i+3)
 connect(4*i+3) = tmp_ents(1+2*i+1)
 end do
 stat_size = 0
 stat_ptr = c_null_ptr
 ! re-use tmp_ents here;
 ! we know it'll always be larger than nquads so it's ok
 call imesh_createentarr(%val(imesh), %val(imesh_quadrilateral), connect, %val(4*nquads), &
 tmp_ents_ptr, tmp_size, tmp_size, stat_ptr, stat_size, stat_size, ierr)
 check(ierr)
 
 ierr = ibase_success
 
 call imesh_freememory(%val(imesh), tmp_ents_ptr)
 call imesh_freememory(%val(imesh), stat_ptr)
 deallocate(connect)
 
 return

References check().

Referenced by directaccessnoholesf90().

directaccessnoholesf90()

program directaccessnoholesf90

Definition at line 45 of file DirectAccessNoHolesF90.F90.

References check(), and create_mesh_no_holes().