added saving of vel moments and histogramm to ou_test

TurTLE/TEST.py

@@ -380,7 +380,13 @@ class TEST(_code):
             if (self.dns_type == 'ornstein_uhlenbeck_test'):
                 with h5py.File(os.path.join(self.work_dir, self.simname + '_ou_output.h5'), 'a') as ofile:
                     ofile.require_group('spectra')
+                    ofile.require_group('moments')
+                    ofile.require_group('histograms')
                     ofile['spectra'].create_dataset('ou_spectra', (1,300,3,3),dtype=np.float)
+                    ofile['moments'].create_dataset('ou', (1,10,3),dtype=np.float)
+                    ofile['histograms'].create_dataset('ou', (1,128,3),dtype=np.float)
+
+
 
         self.run(
                 nb_processes = opt.nb_processes,

cpp/full_code/ornstein_uhlenbeck_test.cpp

@@ -4,7 +4,6 @@
 #include "scope_timer.hpp"
 
 
-
 template <typename rnumber>
 int ornstein_uhlenbeck_test<rnumber>::initialize(void)
 {
@@ -67,19 +66,24 @@ int ornstein_uhlenbeck_test<rnumber>::do_work(void)
     {
         this->ou->step(1e-3);
     }
-    hid_t gg;
+    hid_t out_file;
 
     if (this->myrank == 0) {
-      gg = H5Fopen("test_ou_output.h5", H5F_ACC_RDWR, H5P_DEFAULT);
-      std::cerr << "peep" << std::endl; //DELETE
+      out_file = H5Fopen("test_ou_output.h5", H5F_ACC_RDWR, H5P_DEFAULT);
     }
 
-    // if (this->myrank == 0) gt = H5Fopen((this->simname + std::string(".h5")).c_str(), H5F_ACC_RDWR, H5P_DEFAULT);
-    // this->ou->kk->store(gt);
     this->ou->kk->template cospectrum<rnumber,THREE>(
       this->ou->ou_field->get_cdata(),
       this->ou->ou_field->get_cdata(),
-      gg,"ou_spectra",0);
+      out_file,"ou_spectra",0);
+
+    this->ou->ou_field->ift();
+    std::vector<double> me;
+    me.resize(3);
+    std::fill(me.begin(),me.end(),30);
+    // 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);
+
     return EXIT_SUCCESS;
 }