diff --git a/bfps/cpp/field.cpp b/bfps/cpp/field.cpp
index 6535c07aca723b161bdc6c1b265b48883bae2c8a..2558552dc5c35e70769eda3631323bca689f0fae 100644
--- a/bfps/cpp/field.cpp
+++ b/bfps/cpp/field.cpp
@@ -1102,14 +1102,8 @@ void field<rnumber, be, fc>::symmetrize()
             ptrdiff_t cindex1 = this->get_cindex(0, 0, this->clayout->sizes[1] - iz);
             for (int cc = 0; cc < int(ncomp(fc)); cc++)
             {
-                rnumber rtmp = ((*(cdata + cc + ncomp(fc)*cindex0))[0] +
-                                (*(cdata + cc + ncomp(fc)*cindex1))[0])/2;
-                rnumber ctmp = ((*(cdata + cc + ncomp(fc)*cindex0))[1] -
-                                (*(cdata + cc + ncomp(fc)*cindex1))[1])/2;
-                (*(cdata + cc + ncomp(fc)*cindex0))[0] =  rtmp;
-                (*(cdata + cc + ncomp(fc)*cindex1))[0] =  rtmp;
-                (*(cdata + cc + ncomp(fc)*cindex0))[1] =  ctmp;
-                (*(cdata + cc + ncomp(fc)*cindex1))[1] = -ctmp;
+                (*(cdata + cc + ncomp(fc)*cindex1))[0] =  (*(cdata + cc + ncomp(fc)*cindex0))[0];
+                (*(cdata + cc + ncomp(fc)*cindex1))[1] = -(*(cdata + cc + ncomp(fc)*cindex0))[1];
             }
         }
     }
diff --git a/bfps/test/test_bfps_NSVEparticles.py b/bfps/test/test_bfps_NSVEparticles.py
index 33212e7670728cfeb1f180d2dc51d37653724e86..7ceab8e540038878212b20858af83829dce33d96 100644
--- a/bfps/test/test_bfps_NSVEparticles.py
+++ b/bfps/test/test_bfps_NSVEparticles.py
@@ -42,13 +42,16 @@ def main():
     for iteration in [0, 32, 64]:
         field0 = f0['vorticity/complex/{0}'.format(iteration)].value
         field1 = f1['vorticity/complex/{0}'.format(iteration)].value
-        assert(np.max(np.abs(field0 - field1)) < 1e-5)
+        field_error = np.max(np.abs(field0 - field1))
         x0 = f0['tracers0/state/{0}'.format(iteration)].value
         x1 = f1['tracers0/state/{0}'.format(iteration)].value
-        assert(np.max(np.abs(x0 - x1)) < 1e-5)
+        traj_error = np.max(np.abs(x0 - x1))
         y0 = f0['tracers0/rhs/{0}'.format(iteration)].value
         y1 = f1['tracers0/rhs/{0}'.format(iteration)].value
-        assert(np.max(np.abs(y0 - y1)) < 1e-5)
+        rhs_error = np.max(np.abs(y0 - y1))
+        assert(field_error < 1e-5)
+        assert(traj_error < 1e-5)
+        assert(rhs_error < 1e-5)
     print('SUCCESS! Basic test passed.')
     return None