populate_coloring_result Subroutine

    subroutine populate_coloring_result(self)
        implicit none
        class(type_coloring), intent(inout) :: self
        integer(int32) :: i, j, counts, num_nodes

        num_nodes = size(self%color)
        if (num_nodes == 0) then
            self%num_colors = 0
            return
        end if

        self%num_colors = maxval(self%color)
        if (self%num_colors == 0) return

        if (allocated(self%colored)) deallocate (self%colored)
        allocate (self%colored(self%num_colors))

        do i = 1, self%num_colors
            counts = count(self%color == i)
            self%colored(i)%num_elements = counts
            if (counts > 0) then
                allocate (self%colored(i)%elements(counts))
                self%colored(i)%elements = pack([(j, j=1, num_nodes)], self%color == i)
            end if
        end do
    end subroutine populate_coloring_result