From cab0e2886eb96552f8e8dd08a9b11f662a9e00a4 Mon Sep 17 00:00:00 2001
From: Cristian C Lalescu <Cristian.Lalescu@ds.mpg.de>
Date: Tue, 29 May 2018 15:02:16 +0200
Subject: [PATCH] use own symmetrize, which is fine
---
bfps/cpp/field.cpp | 13 +------------
1 file changed, 1 insertion(+), 12 deletions(-)
diff --git a/bfps/cpp/field.cpp b/bfps/cpp/field.cpp
index e3c35a51..6535c07a 100644
--- a/bfps/cpp/field.cpp
+++ b/bfps/cpp/field.cpp
@@ -1023,6 +1023,7 @@ void field<rnumber, be, fc>::symmetrize()
{
TIMEZONE("field::symmetrize");
assert(!this->real_space_representation);
+ // for debugging, just use FFTW
//this->ift();
//this->dft();
//this->normalize();
@@ -1151,18 +1152,6 @@ void field<rnumber, be, fc>::symmetrize()
(*(cdata + cc + ncomp(fc)*cindex))[1] = 0.0;
}
}
- ///* put asymmetric data to 0 */
- //if (this->clayout->myrank == this->clayout->rank[0][this->clayout->sizes[0]/2])
- //{
- // ptrdiff_t tindex = ncomp(fc)*(this->clayout->sizes[0]/2 - this->clayout->starts[0])*this->clayout->sizes[1]*this->clayout->sizes[2];
- // for (ii = 0; ii < ptrdiff_t(this->clayout->sizes[1]); ii++)
- // {
- // std::fill_n((rnumber*)(data + tindex), ncomp(fc)*2*this->clayout->sizes[2], 0.0);
- // tindex += ncomp(fc)*this->clayout->sizes[2];
- // }
- //}
- ////tindex = ncomp(fc)*();
- ////std::fill_n((rnumber*)(data + tindex), ncomp(fc)*2, 0.0);
}
template <typename rnumber,
--
GitLab