anchor.diff:                    call mpi_bcast(m_np_i,   N, MPI_REAL_XP, i_translated(1), mpi_comm_to_left(i), err)
                  end if
   
  -               ! compute own interactions
  +               !--------------------------------
  +               ! Compute interactions
  +               !--------------------------------
                  a_comm_i = 0._xp
                  a_comm_np_i = 0._xp
                  if (rank == i) then
  
positive.code: 
                 end if
  
                 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)
  
negative.code: 
        if (flag_memory) then
  
           !---------------------------------------------------------------------------
           ! Compute the force using only own particles and communicating each other
           !---------------------------------------------------------------------------
           if (flag_diag) then
              !--------------------------------
              ! Allocate and initialize
  