wip move to addons

CMakeLists.txt

@@ -214,7 +214,8 @@ set(cpp_for_lib
     ${PROJECT_SOURCE_DIR}/cpp/particles/particles_inner_computer.cpp
     ${PROJECT_SOURCE_DIR}/cpp/full_code/ornstein_uhlenbeck_process.cpp
     ${PROJECT_SOURCE_DIR}/cpp/full_code/ornstein_uhlenbeck_test.cpp
-    ${PROJECT_SOURCE_DIR}/cpp/full_code/ou_vorticity_equation.cpp)
+    ${PROJECT_SOURCE_DIR}/cpp/full_code/ou_vorticity_equation.cpp
+    ${PROJECT_SOURCE_DIR}/cpp/full_code/NSVE_ou_forcing.cpp)
 
 set(hpp_for_lib
     ${PROJECT_SOURCE_DIR}/cpp/full_code/code_base.hpp
@@ -292,6 +293,7 @@ set(hpp_for_lib
     ${PROJECT_SOURCE_DIR}/cpp/full_code/ornstein_uhlenbeck_process.hpp
     ${PROJECT_SOURCE_DIR}/cpp/full_code/ornstein_uhlenbeck_test.hpp
     ${PROJECT_SOURCE_DIR}/cpp/full_code/ou_vorticity_equation.hpp
+    ${PROJECT_SOURCE_DIR}/cpp/full_code/NSVE_ou_forcing.hpp
     )
 #file(GLOB_RECURSE hpp_for_lib ${PROJECT_SOURCE_DIR}/*.hpp)
 LIST(APPEND source_files ${hpp_for_lib} ${cpp_for_lib})

TurTLE/DNS.py

@@ -655,6 +655,10 @@ class DNS(_code):
                 'NSVEp_extra_sampling',
                 help = 'plain Navier-Stokes vorticity formulation, with basic fluid tracers, that sample velocity gradient, as well as pressure and its derivatives.')
 
+        parser_NSVE_ou = subparsers.add_parser(
+                'NSVE_ou_forcing',
+                help = 'plain Navier-Stokes vorticity formulation, with ornstein-uhlenbeck forcing')
+
         for parser in ['NSVEparticles_no_output', 'NSVEp2', 'NSVEp2p', 'NSVEp_extra']:
             eval('self.simulation_parser_arguments({0})'.format('parser_' + parser))
             eval('self.job_parser_arguments({0})'.format('parser_' + parser))
@@ -1057,4 +1061,3 @@ class DNS(_code):
                 no_submit = opt.no_submit,
                 no_debug = opt.no_debug)
         return None
-

TurTLE/ou_quick_plot.py

@@ -56,9 +56,9 @@ class plot_field:
         
 
 
-# plot_spectrum()
-# plot_histogram()
+plot_spectrum()
+plot_histogram()
 p = plot_field()
-p.plot('ou_field.h5')
-p.plot('add_ou_test_field.h5')
+#p.plot('ou_field.h5')
+#p.plot('add_ou_test_field.h5')
 plt.show()

cpp/vorticity_equation.hpp

@@ -87,7 +87,7 @@ class vorticity_equation
         ~vorticity_equation(void);
 
         /* solver essential methods */
-        void omega_nonlin(int src);
+        virtual void omega_nonlin(int src);
         void step(double dt);
         void impose_zero_modes(void);
 
@@ -166,4 +166,3 @@ class vorticity_equation
 };
 
 #endif//VORTICITY_EQUATION
-