unique64 Function

private function unique64(vec) result(vec_unique)

Return only the unique values from vec.

Arguments

Type IntentOptional Attributes Name
integer(kind=ip), intent(in), dimension(:) :: vec

Return Value integer(kind=ip), dimension(:), allocatable


Calls

proc~~unique64~~CallsGraph proc~unique64 unique64 interface~sort sort proc~unique64->interface~sort proc~sort_ascending sort_ascending interface~sort->proc~sort_ascending proc~sort_ascending_64 sort_ascending_64 interface~sort->proc~sort_ascending_64 interface~swap swap proc~sort_ascending->interface~swap proc~swap64 swap64 proc~sort_ascending_64->proc~swap64 interface~swap->proc~swap64 proc~swap32 swap32 interface~swap->proc~swap32 proc~swap_str swap_str interface~swap->proc~swap_str

Called by

proc~~unique64~~CalledByGraph proc~unique64 unique64 interface~unique unique interface~unique->proc~unique64

Source Code

    function unique64(vec) result(vec_unique)

        integer(ip),dimension(:),intent(in) :: vec
        integer(ip),dimension(:),allocatable :: vec_unique

        integer(ip) :: i,num
        logical,dimension(size(vec)) :: mask

        mask = .false.

        do i=1,size(vec)

            !count the number of occurrences of this element:
            num = count( vec(i)==vec )

            if (num==1) then
                !there is only one, flag it:
                mask(i) = .true.
            else
                !flag this value only if it hasn't already been flagged:
                if (.not. any(vec(i)==vec .and. mask) ) mask(i) = .true.
            end if

        end do

        !return only flagged elements:
        allocate( vec_unique(count(mask)) )
        vec_unique = pack( vec, mask )

        ! also sort it:
        call sort(vec_unique)

    end function unique64