From 4b534fe30641c07d64b000e566b0dac3bd94f857 Mon Sep 17 00:00:00 2001 From: Oliver Hahn Date: Sun, 23 Aug 2020 22:07:59 +0200 Subject: [PATCH] hard coded phase flip in panphasia --- src/plugins/random_ngenic.cc | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/src/plugins/random_ngenic.cc b/src/plugins/random_ngenic.cc index f8ecd38..6ced6c3 100644 --- a/src/plugins/random_ngenic.cc +++ b/src/plugins/random_ngenic.cc @@ -117,7 +117,7 @@ public: if (i == nres_ / 2 || j == nres_ / 2 || k == nres_ / 2) continue; if (i == 0 && j == 0 && k == 0) continue; - ampl = -std::sqrt(-std::log(ampl)); + ampl = std::sqrt(-std::log(ampl)); ccomplex_t zrand(ampl*std::cos(phase),ampl*std::sin(phase)); if (k > 0) { @@ -143,7 +143,7 @@ public: } } else if (i < nres_ / 2) { if(i_in_range) g.kelem(ip,j,k) = zrand; - if (ii_in_range) g.kelem(iip,jj,k) = std::conj(zrand); + if(ii_in_range) g.kelem(iip,jj,k) = std::conj(zrand); } } }