anchor.diff:              iend2   = iend*2
               select case (flag_compute_force)
                  case(0)
  -                  call compute_force(m, r, istart2, iend2, r, istart2, iend2, a)          ! Compute a(t+dt) with sequential version
  +                  call compute_force(m, r, istart2, iend2, r, istart2, iend2, a)         ! Compute a(t+dt) with sequential version
                  case(1)
  -                  call compute_force_diag(m, r, istart, iend, npoints, a) ! Compute a(t+dt) with fast OpenMP version
  +                  call compute_force_diag(m, r, istart, iend, npoints, a)                ! Compute a(t+dt) with fast OpenMP version
                  case(2)
  -                  call compute_force_omp(m, r, istart2, iend2, r, istart2, iend2, a)      ! Compute a(t+dt) with naive OpenMP version
  +                  call compute_force_omp(m, r, istart2, iend2, r, istart2, iend2, a)     ! Compute a(t+dt) with naive OpenMP version
                  case(3)
  -                  call compute_force_omp_diag(m, r, istart, iend, npoints, a) ! Compute a(t+dt) with fast OpenMP version
  +                  call compute_force_omp_diag(m, r, istart, iend, npoints, a)            ! Compute a(t+dt) with fast OpenMP version
                  case default
                     stop "Unknown value of flag_compute_force"
               end select
  
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:                   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
  
     end subroutine compute_force_wrap
  
     subroutine compute_energy (m, r, v, Ec, Ep, E)
  