anchor.diff:    allocate(pos1(3, N))
     allocate(pos2(3, N))
   
  -  call tester% init()
  +  call test% init()
   
     call random_seed(size = seed_size)
     allocate(seed(seed_size))
  
positive.code:   implicit none
  
    type(cell_system_t) :: solvent_cells
    type(tester_t) :: tester
  
    integer, parameter :: N = 1000000
    double precision, target, allocatable :: pos1(:, :), pos2(:, :)
  
negative.code:      end do
       c1 = floor( (pos(:, i) - solvent_cells% origin) / solvent_cells% a )
       c2 = compact_h_to_p(h-1, solvent_cells% M)
       call test% assert_equal(c1, c2)
    end do
  
    call test% print()
  
    deallocate(pos1)
    deallocate(pos2)
  