MAKE_SPX_NODES.f90 File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	make_spx_nodes (nb_node, spx_con_nnz, spx_con, node_wgt, nnz
	Makes a pointer for spectral connectivity vector.

Function/Subroutine Documentation

make_spx_nodes()

subroutine make_spx_nodes	(	integer*4	nb_node,
		integer*4	spx_con_nnz,
		integer*4, dimension(0:spx_con_nnz)	spx_con,
		integer*4, dimension(nb_node)	node_wgt,
		integer*4	nnz
	)

Makes a pointer for spectral connectivity vector.

Author

Ilario Mazzieri

Date

September, 2013

Version

1.0

Parameters

[in]	nb_node	GLL nodes
[in]	spx_con	spectral connectivity vector
[in]	node_wgt	node_wgt(i) contains the multiplicity of the GLL node i (e.g. if the GLL node i is shared by 4 elements node_wgt(i) = 4
[in]	nnz	number of GLL nodes + number of GLL nodes including repetitions + 1
[out]	node_pointer	as explained in MAKE_GRID_NODES.f90

Definition at line 30 of file MAKE_SPX_NODES.f90.

      
  
      
      implicit none
      
      integer*4 :: nb_node,spx_con_nnz,nnz
      integer*4 :: i,j,ie,nb_elem,nn3
 
      integer*4, dimension(:), allocatable :: i4count
      integer*4, dimension(0:spx_con_nnz) :: spx_con
      integer*4, dimension(nb_node) :: node_wgt
      integer*4, dimension(0:nnz) :: node_pointer
      
      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
      
      nb_elem = spx_con(0) - 1
      do ie = 1, nb_elem
      
         nn3 = spx_con(ie) - spx_con(ie -1) -1
         do i = 1, nn3
      
            j = spx_con(spx_con(ie -1) +i)       
            node_pointer(i4count(j)) = ie
            i4count(j) = i4count(j) + 1
         enddo
      enddo
      
      deallocate(i4count)
      
      return