anchor.diff:                    end if
                  end if
               end do
  -
  +            deallocate(a_comm)
  +            deallocate(a_right)
  +            deallocate(r_i)
  +            deallocate(r_np_i)
  +            deallocate(r_right)
            case default
               stop "Unknown value of flag_compute_mpi"
         end select
   
  -      deallocate(a_comm)
  -      deallocate(a_right)
  -      deallocate(a_reduced)
  -      deallocate(r_i)
  -      deallocate(r_np_i)
  -      deallocate(r_right)
  -
      end subroutine compute_force_wrap
   
      subroutine compute_energy (m, r, v, Ec, Ep, E)
  
positive.code:             !--------------------------------
              call mpi_allreduce(a, a_reduced, npoints*3, MPI_REAL_XP, MPI_SUM, MPI_COMM_WORLD, err)
              a = a_reduced
           case(1)
              !--------------------------------
              ! Scatter positions across all processes
              ! and compute interactions
  
negative.code:       integer         :: istart, iend, istart2, iend2
  !      integer         :: jstart, jend
  
        select case(flag_compute_mpi)
           case(0)
              allocate(a_reduced(3, npoints))
              !--------------------------------
              ! Get the domain for integration from the number of nodes
              !--------------------------------
  