anchor.diff:    type(cell_system_t) :: solvent_cells
     type(tester_t) :: tester
   
  -  integer, parameter :: N = 2000
  -  double precision, target :: pos1(3, N), pos2(3, N)
  +  integer, parameter :: N = 1000000
  +  double precision, target, allocatable :: pos1(:, :), pos2(:, :)
     double precision, pointer :: pos(:,:), pos_old(:,:), pos_pointer(:,:)
   
     integer :: i, L(3), seed_size, clock
     integer, allocatable :: seed(:)
     integer :: h, c1(3), c2(3)
  +
  +  allocate(pos1(3, N))
  +  allocate(pos2(3, N))
   
     call tester% init()
   
  
positive.code: 
    call solvent_cells%count_particles(pos)
  
    call sort
    call solvent_cells%count_particles(pos)
  
    call tester% assert_equal(solvent_cells%cell_start(1), 1)
  
negative.code: 
    call tester% print()
  
    deallocate(pos1)
    deallocate(pos2)
  
  end program try_all
  