cube_add__segment Subroutine

   pure subroutine cube_add__segment(self, segment)
      implicit none
      class(cube_t), intent(inout) :: self
      type(segment_t), intent(in) :: segment

      integer :: cp_min, cp_max, sps_size, i, j
      type(segment_t), allocatable :: tmp(:)
      type(segment_t) :: what_to_add

      if (segment == SEG_EPSILON) then
         call self%flag_epsilon()
         return
      end if

      cp_min = segment%min
      cp_max = segment%max

      if (cp_max < ASCII_SIZE_BIT .and. .not. self%is_switched_to_bmp) then
         call self%ascii%add(cp_min, cp_max)
         if (self%single_flag) self%single_flag = self%num() == 1
         return
      else
         call self%switch_bmp()
         call self%bmp%add(cp_min, cp_max)
      end if

      if (cp_max > BMP_SIZE_BIT) then

         what_to_add = segment_t(max(cp_min, BMP_SIZE_BIT), cp_max)

         if (allocated(self%sps)) then
            sps_size = size(self%sps, dim=1) + 1
            allocate(tmp(sps_size))
            j = 0
            do i = 1, size(self%sps)
               j = j + 1
               if (self%sps(i)%min < what_to_add%min) then
                  tmp(j) = self%sps(i)
               else
                  tmp(j) = what_to_add
               end if
            end do
            self%sps = tmp(1:sps_size) ! implicit reallocation

         else
            self%sps = [segment]
         end if

      end if

      if (self%single_flag) self%single_flag = self%num() == 1

   end subroutine cube_add__segment