anchor.diff:              deallocate(m_i)
               deallocate(m_np_i)
   
  +
  +         !---------------------------------------------------------------------------
  +         ! Compute the force using only own particles and communicating a subset to
  +         ! each other to spare memory
  +         !---------------------------------------------------------------------------
            else !FIXME: Does not work with different masses
               s = N / memory_factor
               if (s*memory_factor /= N) then
  
positive.code:             allocate(m_np_i(N))
  
              !--------------------------------
              ! Scatter positions across all processes
              ! and compute interactions
              !--------------------------------
  
              do i = 0, nprocs - 1
  
negative.code:                   m_np_i = m
                 end if
  
                 !--------------------------------
                 ! Scatter positions and masses
                 !--------------------------------
                 if (communicate_right) then
                    i_to_translate(1) = i
                    call mpi_group_translate_ranks(wgroup, 1, i_to_translate, mpi_group_to_right(i), i_translated, err)
  