update replace method

cpp/full_code/ornstein_uhlenbeck_process.cpp

@@ -177,15 +177,16 @@ void ornstein_uhlenbeck_process<rnumber,be>::add_to_field_replace(
             ptrdiff_t zindex,
             double k2){
 
-              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);
+              if (k2 <= this->ou_kmax_squ && k2 >= this->ou_kmin_squ){
 
-                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);
+                rnumber tmp;
+
+                for(int cc=0; cc < 3; cc++){
+                  for(int imag=0; imag < 2; imag++){
+                    tmp = this->ou_field->cval(cindex,cc,imag);
+                    src->cval(cindex,cc,imag) = tmp;
+                  }
+                }
               }
 
       }