diff --git a/bfps/cpp/field.cpp b/bfps/cpp/field.cpp
index e3c35a517f092b06f63b4c91754d313f7665a18e..6535c07aca723b161bdc6c1b265b48883bae2c8a 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,