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().