MAKE_REN_LOC_NODE.f90 File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	make_ren_loc_node (nnz_loc, cs_loc, nn_loc, id, loc_n_nu
	Renumerates the grid nodes.

Function/Subroutine Documentation

make_ren_loc_node()

subroutine make_ren_loc_node	(	integer*4	nnz_loc,
		integer*4, dimension(0:nnz_loc), intent(in)	cs_loc,
		integer*4, intent(in)	nn_loc,
		integer*4	id,
			loc_n_nu
	)

Renumerates the grid nodes.

Author

Ilario Mazzieri

Date

September, 2013

Version

1.0

Parameters

[in]	nnz_loc	dimnsion for cs_loc
[in]	cs_loc	local connectivity vector
[in]	nn_loc	number of local unknowns
[in]	id	number of mpi process
[in,out]	loc_n_num	local node numbering (renumeration from local to global, loc_n_num(i) = global id for the i-th node)

Definition at line 31 of file MAKE_REN_LOC_NODE.f90.

 
      
      use qsort_c_module
      
      implicit none
 
      integer*4 :: nnz_loc
      integer*4 :: ie, ne, i, j, ind_nd, id, ic
 
      integer*4, intent(in) :: nn_loc
 
      integer*4, dimension(0:nnz_loc), intent(in) :: cs_loc
      integer*4, dimension(nn_loc), intent(inout) :: loc_n_num
      integer*4, dimension(0:nnz_loc) :: cs_loc_copy
 
 
            
      cs_loc_copy = cs_loc
      ne = cs_loc_copy(0) - 1
           
      
      do ie = 1,ne
         cs_loc_copy(cs_loc_copy(ie-1)) = 0       
      enddo
 
      !QUICKSORT FOR REORDING
      !------------------------------------------------------------
      do ie = 0, cs_loc_copy(0)
         cs_loc_copy(ie) = 0
      enddo         
      
      call qsortc(cs_loc_copy)
      
 
      ic = 0      
      do i = 1, nnz_loc
         if (cs_loc_copy(i) .ne. 0) then
               if(cs_loc_copy(i) .ne. cs_loc_copy(i-1)) then
                  ic = ic + 1
                  loc_n_num(ic) = cs_loc_copy(i)               
               endif                
         endif           
      enddo
            
 
            

References qsort_c_module::qsortc().

Here is the call graph for this function: