MAKE_WGT_SPX_NODES.f90 File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	make_wgt_spx_nodes (nb_node, spx_con_nnz, spx_con, node_wgt
	Computes multeplicity for spectral nodes.

Function/Subroutine Documentation

make_wgt_spx_nodes()

subroutine make_wgt_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
	)

Computes multeplicity for spectral nodes.

Author

Ilario Mazzieri

Date

September, 2013

Version

1.0

Parameters

[in]	nb_node	GLL nodes
[in]	spx_con	spectral connectivity vector
[out]	node_wgt	node_wgt(i) contains the multiplicity of the i-th GLL node (e.g. if the node i is shared by 4 elements node_wgt(i) = 4
[out]	nb_nz	number of GLL nodes + number of GLL nodes including repetitions + 1

Definition at line 29 of file MAKE_WGT_SPX_NODES.f90.

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