The missed vectorization seems related to the loop index being unknown: if I replace 'do i = 1, n' with 'do i = 1, nl' with nl set to some given value (say 100) the loop vectorize with -fgraphite-identity. Dominique