From 872b993bbadd9c22a581be819c8a03e6fb13b95c Mon Sep 17 00:00:00 2001
From: Cristian C Lalescu <Cristian.Lalescu@mpcdf.mpg.de>
Date: Mon, 13 Sep 2021 20:47:43 +0200
Subject: [PATCH] fixes symmetrize simd issue
Vectorization was wrongly requested for a few symmetrize loops, where
different arrays are accessed non-sequentially.
For safety, I removed all of them.
---
cpp/field.cpp | 4 ----
1 file changed, 4 deletions(-)
diff --git a/cpp/field.cpp b/cpp/field.cpp
index f6a65bde..49e976e3 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
--
GitLab