strip ou from field

cpp/full_code/ornstein_uhlenbeck_process.cpp

@@ -175,6 +175,34 @@ void ornstein_uhlenbeck_process<rnumber, be>::let_converge(void)
   }
 }
 
+template <class rnumber, field_backend be>
+void ornstein_uhlenbeck_process<rnumber, be>::strip_from_field(
+  field<rnumber, be, THREE> *src)
+{
+  assert(src->real_space_representation==false);
+
+  this->kk->CLOOP_K2(
+          [&](ptrdiff_t cindex,
+            ptrdiff_t xindex,
+            ptrdiff_t yindex,
+            ptrdiff_t zindex,
+            double k2){
+
+              if (k2 <= this->ou_kmax_squ && k2 >= this->ou_kmin_squ){
+
+                for(int cc=0; cc < 3; cc++){
+                  for(int imag=0; imag < 2; imag++){
+                    src->cval(cindex,cc,imag) = 0;
+                  }
+                }
+              }
+
+      }
+
+  );
+
+}
+
 template <class rnumber, field_backend be>
 void ornstein_uhlenbeck_process<rnumber,be>::add_to_field_replace(
   field<rnumber, be, THREE> *src, std::string uv)

cpp/full_code/ornstein_uhlenbeck_process.hpp

@@ -62,6 +62,9 @@ class ornstein_uhlenbeck_process{
         void add_to_field_replace(
           field<rnumber,be,THREE> *src, std::string uv);
 
+        void strip_from_field(
+          field<rnumber, be, THREE> *src);
+
         void calc_ou_vorticity(void);