diff --git a/cpp/field.cpp b/cpp/field.cpp
index f6a65bde9e65a853832d31765849e281eb42454d..49e976e3745e6e2ffd9422934e796734c7e64475 100644
--- a/cpp/field.cpp
+++ b/cpp/field.cpp
@@ -2180,7 +2180,6 @@ void field<rnumber, be, fc>::symmetrize()
             if (this->clayout->myrank == rankp)
             {
                 ptrdiff_t iyy = iy - this->clayout->starts[0];
-                #pragma omp simd
                 for (ptrdiff_t iz = zstart; iz < zend; iz++)
                 {
                     ptrdiff_t cindex = this->get_cindex(0, iyy, iz);
@@ -2194,7 +2193,6 @@ void field<rnumber, be, fc>::symmetrize()
             if (this->clayout->myrank == rankm)
             {
                 ptrdiff_t iyy = (this->clayout->sizes[0] - iy) - this->clayout->starts[0];
-                #pragma omp simd
                 for (ptrdiff_t iz = zstart; iz < zend; iz++)
                 {
                     ptrdiff_t cindex = this->get_cindex(0, iyy, iz);
@@ -2249,7 +2247,6 @@ void field<rnumber, be, fc>::symmetrize()
             if (this->clayout->myrank == rankp)
             {
                 ptrdiff_t iyy = iy - this->clayout->starts[0];
-                #pragma omp simd
                 for (ptrdiff_t iz = zstart; iz < zend; iz++)
                 {
                     // modulo operation makes sense only for slab decomposition
@@ -2266,7 +2263,6 @@ void field<rnumber, be, fc>::symmetrize()
             if (this->clayout->myrank == rankm)
             {
                 ptrdiff_t iyy = (this->clayout->sizes[0] - iy) - this->clayout->starts[0];
-                #pragma omp simd
                 for (ptrdiff_t iz = zstart; iz < zend; iz++)
                 {
                     // modulo operation makes sense only for slab decomposition