Skip to content
Snippets Groups Projects
Commit cab0e288 authored by Cristian Lalescu's avatar Cristian Lalescu
Browse files

use own symmetrize, which is fine

parent a6247072
Branches
Tags
1 merge request!23WIP: Feature/use cmake
Pipeline #
...@@ -1023,6 +1023,7 @@ void field<rnumber, be, fc>::symmetrize() ...@@ -1023,6 +1023,7 @@ void field<rnumber, be, fc>::symmetrize()
{ {
TIMEZONE("field::symmetrize"); TIMEZONE("field::symmetrize");
assert(!this->real_space_representation); assert(!this->real_space_representation);
// for debugging, just use FFTW
//this->ift(); //this->ift();
//this->dft(); //this->dft();
//this->normalize(); //this->normalize();
...@@ -1151,18 +1152,6 @@ void field<rnumber, be, fc>::symmetrize() ...@@ -1151,18 +1152,6 @@ void field<rnumber, be, fc>::symmetrize()
(*(cdata + cc + ncomp(fc)*cindex))[1] = 0.0; (*(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, template <typename rnumber,
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment