added SIMD directive to code to enable vectorization