Research Article

Using Coarrays to Parallelize Legacy Fortran Applications: Strategy and Case Study

Listing 2

Legacy particle loop.
l = 0 ! Global particle number
do k = 1, nb ! Loop over the bands
    do m = 1, k ! Loop over the particles in band
    ! First octant
   l = l + 1
    ! Do some computations
    ! Second octant
   l = l + 1
    ! Do some computations
    ! Third octant
   l = l + 1
    ! Do some computations
    ! Fourth octant
   l = l + 1
    ! Do some computations
    end do
end do