diff --git a/cpp/full_code/ornstein_uhlenbeck_test.cpp b/cpp/full_code/ornstein_uhlenbeck_test.cpp
index 64bb1ecb178c45de23b5d1c57de690934938df17..b9e485a3504cb2f9404f71fc9e508ceca30cfc53 100644
--- a/cpp/full_code/ornstein_uhlenbeck_test.cpp
+++ b/cpp/full_code/ornstein_uhlenbeck_test.cpp
@@ -77,6 +77,7 @@ int ornstein_uhlenbeck_test<rnumber>::do_work(void)
       this->ou->ou_field->get_cdata(),
       out_file,"ou_spectra",0);
 
+    this->ou->ou_field->symmetrize();
     this->ou->ou_field->ift();
     std::vector<double> me;
     me.resize(3);
@@ -84,6 +85,8 @@ int ornstein_uhlenbeck_test<rnumber>::do_work(void)
     // this->ou->ou_field->template compute_rspace_stats<rnumber, FFTW, THREE>(out_file, "ou", 0, me);
     this->ou->ou_field->compute_rspace_stats(out_file, "ou", 0, me);
 
+    this->ou->ou_field->io("ou_field.h5", "ou_field", 0, false);
+
     return EXIT_SUCCESS;
 }