ut_sortOnProcs.F90

Go to the documentation of this file.

!! NOTICE
!!  Copyright 2022 UChicago Argonne, LLC and contributors
!!
!!  Licensed under the Apache License, Version 2.0 (the "License");
!!  you may not use this file except in compliance with the License.
!!
!!  Unless required by applicable law or agreed to in writing, software
!!  distributed under the License is distributed on an "AS IS" BASIS,
!!  WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
!!  See the License for the specific language governing permissions and
!!  limitations under the License.
!!
 
 
subroutine ut_sortOnProcs(count, props, attrib, numProcs,&
     storage, workspace, &
     perProc, ifNonZeroProc,nonZeroProcsCount)
 
#include "constants.h"
#include "Simulation.h"
 
  implicit none
 
  integer, intent(IN) :: count, props, attrib, numProcs
  real,dimension(props,count),intent(INOUT) ::storage,workspace
  integer,dimension(numProcs),intent(OUT) :: perProc, ifNonZeroProc
  integer, intent(OUT) :: nonZeroProcsCount
 
  integer :: i,j,k,n
 
  integer, dimension(numProcs) :: pntr
 
  perProc=0
  nonZeroProcsCount=0
  ifNonZeroProc=0
 
  if(count> 0) then
     
     k=1
     n=0
     do i = 1,count
        j=int(storage(attrib,i))
        if(j/=NONEXISTENT) then
           n=n+1
           workspace(1:props,n)=storage(1:props,i)
           j=j+1
           perProc(j)=perProc(j)+1
           if(ifNonZeroProc(j)==0) then
              ifNonZeroProc(j)=1
              nonZeroProcsCount=nonZeroProcsCount+1
           end if
        end if
     end do
     
     pntr(1)=1
     do i = 2,numProcs
        pntr(i)=pntr(i-1)+perProc(i-1)
     end do
     
     do i = 1,n
        j=int(workspace(attrib,i))+1
        k = pntr(j)
        storage(1:props,k)=workspace(1:props,i)
        pntr(j)=pntr(j)+1
     end do
  end if
end subroutine ut_sortOnProcs