anchor.diff:                    call minmax(value, min, max)
                 end if
   
  -              call centroid(pos, min, max, dim)
  +              call centroid(pos, min, dim)
   
   
             !! Contract from old position if minima value is still minima value
  
positive.code:       real*8, dimension(:,:), intent(inout) :: pos
        integer, intent(in)                   :: max, dim
        integer                               :: i
        real*8                                :: xsum, ysum
  
        do i = 1, dim - 1
            if (i.NE.max) then
  
negative.code:       implicit none
        integer                :: dim
        real*8, dimension(:,:) :: pos
        integer                :: i
  
        !!call init_random_seed()
  
  