expand_ga_ellipsoid Subroutine

private subroutine expand_ga_ellipsoid(me, pos, length)

Arguments

Type IntentOptional Attributes Name
type(grw_ellipsoidarray_type) :: me
integer, intent(in), optional :: pos
integer, intent(in), optional :: length

optional length to expand the array


Called by

proc~~expand_ga_ellipsoid~~CalledByGraph proc~expand_ga_ellipsoid expand_ga_ellipsoid interface~expand~27 expand interface~expand~27->proc~expand_ga_ellipsoid proc~placeat_ga_ellipsoid placeat_ga_ellipsoid proc~placeat_ga_ellipsoid->interface~expand~27 proc~placeat_ga_ellipsoid_vec placeat_ga_ellipsoid_vec proc~placeat_ga_ellipsoid_vec->interface~expand~27 proc~append_ga_ellipsoid append_ga_ellipsoid proc~append_ga_ellipsoid->interface~expand~27 proc~append_ga_ellipsoid_vec append_ga_ellipsoid_vec proc~append_ga_ellipsoid_vec->interface~expand~27 interface~placeat~17 placeat interface~placeat~17->proc~placeat_ga_ellipsoid interface~placeat~17->proc~placeat_ga_ellipsoid_vec interface~append~28 append interface~append~28->proc~append_ga_ellipsoid interface~append~28->proc~append_ga_ellipsoid_vec

Contents

Source Code


Source Code

  subroutine expand_ga_ellipsoid(me, pos, length)
    type(grw_ellipsoidarray_type) :: me !< array to resize
    integer, intent(in), optional :: pos !< optional predefined position
    !> optional length to expand the array
    integer, intent(in), optional :: length

    type(tem_ellipsoid_type), allocatable :: swpval(:)
    integer :: explen, ii

    explen = 0
    ! increase the container by the requested length of double it
    if( present(length) ) then
      explen = max( length, minlength )
    else
      ! set the global minimum length, if doubling would be smaller than that
      explen = max(me%containersize, minlength)
    end if

    ! if a position is given, increase the container to at least the size to
    ! fit the position.
    if( present(pos) ) explen = max(explen, pos-me%containersize)

    ! if the current size plus explen exceeds the max container size,
    ! reduce the size to the max container size.
    if( (huge(me%containersize) - explen) <= me%containersize) then
      ! set max container size
      me%containersize = huge(me%containersize)
    else
      ! set the new container size
      me%containersize = me%containersize + explen
    end if

    if ( me%nvals > 0 ) then
      allocate(swpval(me%containersize))
      do ii = 1, me%nvals
        swpval(ii) = me%val(ii)
      end do
      call move_alloc( swpval, me%val )
    else ! me%nvals == 0
      if ( allocated(me%val) ) deallocate( me%val )
      allocate( me%val(me%containersize) )
    end if

  end subroutine expand_ga_ellipsoid