diff --git a/cpp/full_code/ornstein_uhlenbeck_process.cpp b/cpp/full_code/ornstein_uhlenbeck_process.cpp
index 42cdfdcb868b6969271bc2dfc377f13ff02bc29f..101093e940b2f06e6acafd3fde75458303c524f6 100644
--- a/cpp/full_code/ornstein_uhlenbeck_process.cpp
+++ b/cpp/full_code/ornstein_uhlenbeck_process.cpp
@@ -174,9 +174,6 @@ void ornstein_uhlenbeck_process<rnumber,be>::add_to_field_replace(
                 src->cval(cindex,1,1) = this->ou_field->cval(cindex,1,1);
                 src->cval(cindex,2,1) = this->ou_field->cval(cindex,2,1);
               }
-              // else {
-              //   assert(src->cval(cindex,0,0) == this->ou_field->cval(cindex,0,0));
-              // }
 
       }
 
@@ -197,7 +194,13 @@ void ornstein_uhlenbeck_process<rnumber,be>::add_to_field_sharp(
 
               if (k2 <= this->ou_kmax_squ && k2 >= this->ou_kmin_squ)
               {
+                src->cval(cindex,0,0) += this->ou_field->cval(cindex,0,0);
+                src->cval(cindex,1,0) += this->ou_field->cval(cindex,1,0);
+                src->cval(cindex,2,0) += this->ou_field->cval(cindex,2,0);
 
+                src->cval(cindex,0,1) += this->ou_field->cval(cindex,0,1);
+                src->cval(cindex,1,1) += this->ou_field->cval(cindex,1,1);
+                src->cval(cindex,2,1) += this->ou_field->cval(cindex,2,1);
               }
 
       }