GET_LOC_NODE_NUM.f90 File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	get_loc_node_num (nnz_loc, cs_loc, nn_loc, id)
	Computes number of local nodes.

Function/Subroutine Documentation

get_loc_node_num()

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

Computes number of local nodes.

Author

Ilario Mazzieri

Date

September, 2013

Version

1.0

Parameters

[in]	nnz_loc	length of cs_loc
[in]	cs_loc	local connectivity vector
[in]	id	number of mpi process
[out]	nn_loc	number of local unknowns

Definition at line 29 of file GET_LOC_NODE_NUM.f90.

 
      use qsort_c_module
      
      implicit none
 
      integer*4 :: nnz_loc
      integer*4 :: ie, ne, i, j, ind_nd, id
 
      integer*4, intent(out) :: nn_loc
 
      integer*4, dimension(0:nnz_loc), intent(in) :: cs_loc
      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
      !------------------------------------------------------------
      nn_loc = 0
      do ie = 0, cs_loc_copy(0)
         cs_loc_copy(ie) = 0
      enddo         
      
      call qsortc(cs_loc_copy)
      
      
      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
                  nn_loc = nn_loc + 1
               endif                
         endif           
      enddo
            
      !-------------------------------------------------------------
  
      !-------------------------------------------------------------- 
!      do i = cs_loc_copy(0) + 1, nnz_loc
!         ind_nd = cs_loc_copy(i)
!      
!         do j = i + 1, nnz_loc
!            if (ind_nd .eq. cs_loc_copy(j)) cs_loc_copy(j) = 0          
!         enddo
!      enddo
!      
!      do ie = 0, cs_loc_copy(0)
!         cs_loc_copy(ie) = 0
!      enddo         
!      
!      
!      nn_loc = 0
!      do i = 0, nnz_loc - 1
!         if(cs_loc_copy(i) .ne. 0) nn_loc = nn_loc + 1
!      enddo   
      !-------------------------------------------------------------
      
            

References qsort_c_module::qsortc().

Here is the call graph for this function: