MAKE_GRID_NODES.f90 File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	make_grid_nodes (nb_node, nb_elem, con_matrix, node_wgt, nnz
	Makes pointer for connectivity of mesh nodes.

Function/Subroutine Documentation

make_grid_nodes()

subroutine make_grid_nodes	(	integer*4	nb_node,
		integer*4	nb_elem,
		integer*4, dimension(nb_elem,9)	con_matrix,
		integer*4, dimension(nb_node)	node_wgt,
		integer*4	nnz
	)

Makes pointer for connectivity of mesh nodes.

Author

Ilario Mazzieri

Date

September, 2013

Version

1.0

Parameters

[in]	nb_node	grid nodes
[in]	nb_elem	number of elements (hexes)
[in]	con_matrix	con_matrix(i,1) -> element con_matrix(i,2),...,con_matrix(i,9) -> grid nodes of element con_matrix(i,1)
[in]	node_wgt	node_wgt(i) contains the multiplicity of the grid node i (e.g. if the node node i is shared by 4 elements node_wgt(i) = 4
[in]	nnz	number of grid nodes + number of grid nodes including repetitions + 1
[out]	node_pointer	pointer for mesh nodes connectivity

Definition at line 72 of file MAKE_GRID_NODES.f90.

      
      
      implicit none
      
      integer*4 :: nb_node,nb_elem,nnz
      integer*4 :: i,j,ie
      
      integer*4, dimension(:), allocatable :: i4count
      integer*4, dimension(nb_node) :: node_wgt
      integer*4, dimension(0:nnz) :: node_pointer
      
      integer*4, dimension(nb_elem,9) :: con_matrix
      
 
      allocate(i4count(nb_node))
      
      node_pointer(0) = nb_node +1
      do j = 1,nb_node
         node_pointer(j) = node_pointer(j-1) + node_wgt(j)
      enddo
      
      do j = 1,nb_node
         i4count(j) = node_pointer(j-1)
      enddo
      
      
      do ie = 1,nb_elem
         do i = 1,8
            j = con_matrix(ie,i +1)
            node_pointer(i4count(j)) = ie      
            i4count(j) = i4count(j) +1
         enddo
      enddo
      
      deallocate(i4count)
      
      return