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